The inference of antigen selection on Ig genes

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.

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.

Signatures of selection in the human antibody repertoire: Selective sweeps, competing subclones, and neutral drift

Antibodies are created and refined by somatic evolution in B cell populations, which endows the human immune system with the ability to recognize and eliminate diverse pathogens. However, the evolutionary processes that sculpt antibody repertoires remain poorly understood. Here, using an unbiased repertoire-scale approach, we show that the population genetic signatures of evolution are evident in human B cell lineages and reveal how antibodies evolve somatically. We measured the dynamics and genetic diversity of B cell responses in five adults longitudinally before and after influenza vaccination using high-throughput antibody repertoire sequencing. We identified vaccine-responsive B cell lineages that carry signatures of selective sweeps driven by positive selection, and discovered that they often display evidence for selective sweeps favoring multiple subclones. We also found persistent B cell lineages that exhibit stable population dynamics and carry signatures of neutral drift. By exploiting the relationship between B cell fitness and antibody binding affinity, we demonstrate the potential for using phylogenetic approaches to identify antibodies with high binding affinity. This quantitative characterization reveals that antibody repertoires are shaped by an unexpectedly broad spectrum of evolutionary processes and shows how signatures of evolutionary history can be harnessed for antibody discovery and engineering.

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.

Antibody Heavy Chain Variable Domains of Different Germline Gene Origins Diversify through Different Paths

B cells produce antibodies, key effector molecules in health and disease. They mature their properties, including their affinity for antigen, through hypermutation events; processes that involve, e.g., base substitution, codon insertion and deletion, often in association with an isotype switch. Investigations of antibody evolution define modes whereby particular antibody responses are able to form, and such studies provide insight important for instance for development of efficient vaccines. Antibody evolution is also used in vitro for the design of antibodies with improved properties. To better understand the basic concepts of antibody evolution, we analyzed the mutational paths, both in terms of amino acid substitution and insertions and deletions, taken by antibodies of the IgG isotype. The analysis focused on the evolution of the heavy chain variable domain of sets of antibodies, each with an origin in 1 of 11 different germline genes representing six human heavy chain germline gene subgroups. Investigated genes were isolated from cells of human bone marrow, a major site of antibody production, and characterized by next-generation sequencing and an in-house bioinformatics pipeline. Apart from substitutions within the complementarity determining regions, multiple framework residues including those in protein cores were targets of extensive diversification. Diversity, both in terms of substitutions, and insertions and deletions, in antibodies is focused to different positions in the sequence in a germline gene-unique manner. Altogether, our findings create a framework for understanding patterns of evolution of antibodies from defined germline genes.

Rare, high-affinity anti-pathogen antibodies from human repertoires, discovered using microfluidics and molecular genomics

Affinity-matured, functional anti-pathogen antibodies are present at low frequencies in natural human repertoires. These antibodies are often excellent candidates for therapeutic monoclonal antibodies. However, mining natural human antibody repertoires is a challenge. In this study, we demonstrate a new method that uses microfluidics, yeast display, and deep sequencing to identify 247 natively paired anti-pathogen single-chain variable fragments (scFvs), which were initially as rare as 1 in 100,000 in the human repertoires. Influenza A vaccination increased the frequency of influenza A antigen-binding scFv within the peripheral B cell repertoire from <0.1% in non-vaccinated donors to 0.3-0.4% in vaccinated donors, whereas pneumococcus vaccination did not increase the frequency of antigen-binding scFv. However, the pneumococcus scFv binders from the vaccinated library had higher heavy and light chain Replacement/Silent mutation (R/S) ratios, a measure of affinity maturation, than the pneumococcus binders from the corresponding non-vaccinated library. Thus, pneumococcus vaccination may increase the frequency of affinity-matured antibodies in human repertoires. We synthesized 10 anti-influenza A and nine anti-pneumococcus full-length antibodies that were highly abundant among antigen-binding scFv. All 10 anti-influenza A antibodies bound the appropriate antigen at K_D<10 nM and neutralized virus in cellular assays. All nine anti-pneumococcus full-length antibodies bound at least one polysaccharide serotype, and 71% of the anti-pneumococcus antibodies that we tested were functional in cell killing assays. Our approach has future application in a variety of fields, including the development of therapeutic antibodies for emerging viral diseases, autoimmune disorders, and cancer.

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.

Characterization of T and B cell repertoire diversity in patients with RAG deficiency

Recombination-activating genes 1 and 2 (RAG1 and RAG2) play a critical role in T and B cell development by initiating the recombination process that controls the expression of T cell receptor (TCR) and immunoglobulin genes. Mutations in the RAG1 and RAG2 genes in humans cause a broad spectrum of phenotypes, including severe combined immunodeficiency (SCID) with lack of T and B cells, Omenn syndrome, leaky SCID, and combined immunodeficiency with granulomas or autoimmunity (CID-G/AI). Using next-generation sequencing, we analyzed the TCR and B cell receptor (BCR) repertoire in 12 patients with RAG mutations presenting with Omenn syndrome (n = 5), leaky SCID (n = 3), or CID-G/AI (n = 4). Restriction of repertoire diversity skewed usage of variable (V), diversity (D), and joining (J) segment genes, and abnormalities of CDR3 length distribution were progressively more prominent in patients with a more severe phenotype. Skewed usage of V, D, and J segment genes was present also within unique sequences, indicating a primary restriction of repertoire. Patients with Omenn syndrome had a high proportion of class-switched immunoglobulin heavy chain transcripts and increased somatic hypermutation rate, suggesting in vivo activation of these B cells. These data provide a framework to better understand the phenotypic heterogeneity of RAG deficiency.

Non-infectious environmental antigens as a trigger for the initiation of an autoimmune skin disease

Pemphigus represents a group of organ specific autoimmune blistering disorders of the skin mediated by pathogenic autoantibodies with well-defined antigenic targets. While most of these diseases are sporadic, endemic forms of disease do exist. The endemic form of pemphigus foliaceus (also known as fogo selvagem, FS) exhibits epidemiological features that suggest exposure to hematophagous insect bites are a possible precipitating factor of this autoimmune disease, and provides a unique opportunity to study how environmental factors contribute to autoimmune disease development. FS patients and healthy individuals from endemic regions show an autoreactive IgM response that starts in early childhood and becomes restricted to IgG4 autoantibodies in FS patients. In searching for triggering environmental antigens, we have found that IgG4 and IgE autoantibodies from FS patients cross-react with a salivary antigen from sand flies. The presence of these cross-reactive antibodies and antibody genetic analysis confirming that these antibodies evolve from the same naïve B cells provides compelling evidence that this non-infectious environmental antigen could be the initial target of the autoantibody response in FS. Consequently, FS serves as an ideal model to study the impact of environmental antigens in the development of autoimmune disease.

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.

Inferring processes underlying B-cell repertoire diversity

We quantify the VDJ recombination and somatic hypermutation processes in human B cells using probabilistic inference methods on high-throughput DNA sequence repertoires of human B-cell receptor heavy chains. Our analysis captures the statistical properties of the naive repertoire, first after its initial generation via VDJ recombination and then after selection for functionality. We also infer statistical properties of the somatic hypermutation machinery (exclusive of subsequent effects of selection). Our main results are the following: the B-cell repertoire is substantially more diverse than T-cell repertoires, owing to longer junctional insertions; sequences that pass initial selection are distinguished by having a higher probability of being generated in a VDJ recombination event; somatic hypermutations have a non-uniform distribution along the V gene that is well explained by an independent site model for the sequence context around the hypermutation site.

Use of the VH6-1 gene segment to code for anti-interleukin-18 autoantibodies in multiple sclerosis

We investigated whether levels and repertoires of anti-interleukin-18 (IL-18) autoantibodies (auto-Abs) differ in multiple sclerosis (MS) patients and healthy donors (HDs). IL-18 concentration in MS patients' sera was higher than in HD, but the level of anti-IL-18 auto-Abs was lower in MS patients. Correlation patterns of IL-18/anti-IL-18 auto-Abs system differed in HD and MS patients, so we have compared segment composition of the anti-IL-18 single-chain variable fragments (scFvs) selected from MS and naïve phage display libraries. Considerable differences between anti-IL-18 auto-Abs of these libraries were found. MS panel contained auto-Abs displaying both signs of "fetal" and somatically hypermutated repertoires. Naïve panel mainly contained the naïve antibodies. These variations from the norm are possible results of abnormal regulation of the repertoire perhaps determined by remodeling of the molecular mechanisms involved in the V(D)J recombination and/or abnormal selection by antigen in MS pathogenesis.

Clonality Analysis of Immunoglobulin Gene Rearrangement by Next-Generation Sequencing in Endemic Burkitt Lymphoma Suggests Antigen Drive Activation of BCR as Opposed to Sporadic Burkitt Lymphoma

Objectives: Recent studies using next-generation sequencing (NGS) analysis disclosed the importance of the intrinsic activation of the B-cell receptor (BCR) pathway in the pathogenesis of sporadic Burkitt lymphoma (sBL) due to mutations of TCF3/ID3 genes. Since no definitive data are available on the genetic landscape of endemic Burkitt (eBL), we first assessed the mutation frequency of TCF3/ID3 in eBL compared with sBL and subsequently the somatic hypermutation status of the BCR to answer whether an extrinsic activation of BCR signaling could also be demonstrated in Burkitt lymphoma.

Methods: We assessed the mutations of TCF3/ID3 by RNAseq and the BCR status by NGS analysis of the immunoglobulin genes (IGs).

Results: We detected mutations of TCF3/ID3 in about 30% of the eBL cases. This rate is significantly lower than that detected in sBL (64%). The NGS analysis of IGs revealed intraclonal diversity, suggesting an active targeted somatic hypermutation process in eBL compared with sBL.

Conclusions: These findings support the view that the antigenic pressure plays a key role in the pathogenetic pathways of eBL, which may be partially distinct from those driving sBL development.

Aging affects B-cell antigen receptor repertoire diversity in primary and secondary lymphoid tissues

The elderly immune system is characterized by reduced responses to infections and vaccines, and an increase in the incidence of autoimmune diseases and cancer. Age-related deficits in the immune system may be caused by peripheral homeostatic pressures that limit bone marrow B-cell production or migration to the peripheral lymphoid tissues. Studies of peripheral blood B-cell receptor spectratypes have shown that those of the elderly are characterized by reduced diversity, which is correlated with poor health status. In the present study, we performed for the first time high-throughput sequencing of immunoglobulin genes from archived biopsy samples of primary and secondary lymphoid tissues in old (74 ± 7 years old, range 61-89) versus young (24 ± 5 years old, range 18-45) individuals, analyzed repertoire diversities and compared these to results in peripheral blood. We found reduced repertoire diversity in peripheral blood and lymph node repertoires from old people, while in the old spleen samples the diversity was larger than in the young. There were no differences in somatic hypermutation characteristics between age groups. These results support the hypothesis that age-related immune frailty stems from altered B-cell homeostasis leading to narrower memory B-cell repertoires, rather than changes in somatic hypermutation mechanisms.

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.

Vimentin is a dominant target of in situ humoral immunity in human lupus tubulointerstitial nephritis

OBJECTIVE

In lupus nephritis (LN), severe tubulointerstitial inflammation (TII) predicts progression to renal failure. Severe TII is associated with tertiary lymphoid neogenesis and in situ antigen-driven clonal B cell selection. The autoantigen(s) driving in situ B cell selection in TII are not known. This study was undertaken to identify the dominant driving autoantigen(s).

METHODS

Single CD38+ or Ki-67+ B cells were laser captured from 7 biopsy specimens that were diagnostic for LN. Eighteen clonally expanded immunoglobulin heavy- and light-chain variable region pairs were cloned and expressed as monoclonal antibodies. Seven more antibodies were cloned from flow-sorted CD38+ cells from an eighth biopsy specimen. Antigen characterization was performed using a combination of confocal microscopy, enzyme-linked immunosorbent assay, screening protoarrays, immunoprecipitation, and mass spectrometry. Serum IgG titers to the dominant antigen in 48 LN and 35 non-nephritic lupus samples were determined using purified antigen-coated arrays. Autoantigen expression on normal and LN kidney was localized by immunohistochemistry and immunofluorescence.

RESULTS

Eleven of 25 antibodies reacted with cytoplasmic structures, 4 reacted with nuclei, and none reacted with double-stranded DNA. Vimentin was the only autoantigen identified by both mass spectrometry and protoarray. Ten of the 11 anticytoplasmic TII antibodies directly bound vimentin. Vimentin was highly expressed by tubulointerstitial inflammatory cells, and the TII antibodies tested preferentially bound inflamed tubulointerstitium. Finally, high titers of serum antivimentin antibodies were associated with severe TII (P = 0.0001).

CONCLUSION

Vimentin, an antigenic feature of inflammation, is a dominant autoantigen targeted in situ in LN TII. This adaptive autoimmune response likely feeds forward to worsen TII and renal damage.

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.

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.

Biased immunoglobulin light chain use in the Chlamydophila psittaci negative ocular adnexal marginal zone lymphomas

Ocular adnexal mucosa associated lymphoid tissue lymphomas (OAMALTL) are the most common lymphomas of the eye. The potential roles for specific antigens in these lymphomas are still controversial. Previously we examined IGHV usage and mutations in Chlamydophila (C) psittaci-negative OAMALTL, demonstrating biased use of the IGHV4 family and IGHV4-34 gene and evidence for antigen selection. Herein, we examined the IGKV/IGLV gene usage and mutations in 34 C. psittaci-negative OAMALTL originating from the orbit (15), conjunctivae (14), and lacrimal gland (5). Clonal potentially functional IGKV/IGLV gene sequences were identified in 30 tumors (18 kappa and 12 lambda). An overrepresentation of the IGKV4 family (P < 0.01) was observed. The IGKV3-20*01 allele was used at a greater frequency than in normal peripheral blood B-lymphocytes (P = 0.02) and commonly paired with the IGHV4-34 allele. Twenty-seven of the 30 unique light chain sequences displayed mutations from germline and evidence for antigen selection. Overall our findings demonstrate that in C. psittaci-negative OAMALTL there is a biased usage of IGKV families and genes, which harbor somatic mutations. These findings and the specific paring between the IGKV3-20*01 and IGHV4-34 alleles suggest that specific antigens could play an important role in the pathogenesis of these lymphomas.

Comparison of two functional kappa light-chain transcripts amplified from a hybridoma

Three heavy-chain and three kappa (κ)-chain transcripts were amplified from hybridoma cells secreting a monoclonal antibody (mAb) against transferrin receptor. Sequence analysis via IMGT/V-QUEST yielded the functional/aberrant prediction. Two functional κ-chain transcripts, Vκ2 and Vκ3, and one functional VH1 were revealed. Comprehensive bioinformatics analyses including sequence alignment, phylogenetic tree, somatic hypermutation prediction, and three-dimensional-molecular structure modeling were used to predict the origin of the two κ-chain transcripts. The results of bioinformatics analysis suggest that Vκ3 is derived from the myeloma partner of the hybridoma; Vκ2 is derived from B-cell. Functional transcripts VH1 and Vκ2 and Vκ3 were then used to construct two chimeric antibodies chi-C2 (Vκ2-VH1) and chi-C3 (Vκ3-VH1), respectively. Antigen-binding experiments showed that only chi-C2 remained the same affinity as its parental mAb. Possible explanations for the coexistence of two functional κ-chain transcripts and the different affinity of the two chimeric antibodies are discussed.

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.

Immunoglobulin light chains in medaka (Oryzias latipes)

The gene segments encoding antibodies have been studied in many capacities and represent some of the best-characterized gene families in traditional animal disease models (mice and humans). To date, multiple immunoglobulin light chain (IgL) isotypes have been found in vertebrates and it is unclear as to which isotypes might be more primordial in nature. Sequence data emerging from an array of fish genome projects is a valuable resource for discerning complex multigene assemblages in this critical branch point of vertebrate phylogeny. Herein, we have analyzed the genomic organization of medaka (Oryzias latipes) IgL gene segments based on recently released genome data. The medaka IgL locus located on chromosome 11 contains at least three clusters of IgL gene segments comprised of multiple gene assemblages of the kappa light chain isotype. These data suggest that medaka IgL gene segments may undergo both intra- and inter-cluster rearrangements as a means to generate additional diversity. Alignments of expressed sequence tags to concordant gene segments which revealed each of the three IgL clusters are expressed. Collectively, these data provide a genomic framework for IgL genes in medaka and indicate that Ig diversity in this species is achieved from at least three distinct chromosomal regions.

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.

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.

Biased use of the IGHV4 family and evidence for antigen selection in Chlamydophila psittaci-negative ocular adnexal extranodal marginal zone lymphomas

Extranodal marginal zone lymphomas (EMZL) are the most common lymphomas in the ocular adnexa. The etiology and potential role for antigenic stimulation in these lymphomas are still controversial. We have examined IGHV gene usage and mutations in 67 Chlamydophila psittaci-negative ocular adnexal EMZL. Clonal IGHV gene sequences were identified in 43 tumors originating from the orbit (19), conjunctivae (18) and lacrimal gland (6). Forty four potentially functional clonal IGHV gene sequences were detected with overrepresentation of the IGHV4 family and IGHV4-34 gene. All but 3 sequences were mutated with the average percent homology to the germ line of 93.5±6.1. Multinomial model and Focused binomial test demonstrated evidence for positive and/or negative antigen selection in 59% of the potentially functional IGHV genes. Intraclonal variation was detected in 8 of 11 tumor specimens. Overall our findings demonstrate that C. psittaci-negative ocular adnexal EMZL exhibit biased usage of IGHV families and genes with evidence for intraclonal heterogeneity and antigen selection in multiple tumors, implicating B-cell receptor-mediated antigen stimulation in the pathogenesis of these lymphomas.

The repertoire of λ light chains causing predominant amyloid heart involvement and identification of a preferentially involved germline gene, IGLV1-44

Monoclonal Ig light chains (LC) can be responsible for pathologic conditions in humans, as in systemic amyloid light amyloidosis. Protean clinical manifestations characterize this disorder with the most varied combination of symptoms generated by different degrees of diverse organ involvement. Kidney and heart are most frequently interested, with major heart involvement as the most relevant prognostic factor. The identification of the underlying mechanism involved in organ targeting is of major relevance for the pathobiology of this disorder. To this aim, we characterized the repertoire of variable region germline genes of λ LC preferentially targeting the heart and compared it with the repertoire of LC that do not in a case-control study. We found that the repertoires were highly restricted, showing preferential use of the same few germline genes but with a different frequency pattern. A single gene, IGVL1-44, was found associated with a 5-fold increase in the odds of dominant heart involvement (after adjusting for confounders in a multivariable logistic model). These results support an involvement of LC genetics in the determination of organ targeting. Study of the characteristics of IGVL1-44-LC with, and of the minority without, heart involvement might lead to identification of LC/tissue interactions.

Antibody V(h) repertoire differences between resolving and chronically evolving hepatitis C virus infections

Despite the production of neutralizing antibodies to hepatitis C virus (HCV), many patients fail to clear the virus and instead develop chronic infection and long-term complications. To understand how HCV infection perturbs the antibody repertoire and to identify molecular features of antibody genes associated with either viral clearance or chronic infection, we sequenced the V(D)J region of naïve and memory B cells of 6 persons who spontaneously resolved an HCV infection (SR), 9 patients with a newly diagnosed chronically evolving infection (CE), and 7 healthy donors. In both naïve and memory B cells, the frequency of use of particular antibody gene subfamilies and segments varied among the three clinical groups, especially between SR and CE. Compared to CE, SR antibody genes used fewer VH, D and JH gene segments in naïve B cells and fewer VH segments in memory B cells. SR and CE groups significantly differed in the frequency of use of 7 gene segments in naïve B cell clones and 3 gene segments in memory clones. The nucleotide mutation rates were similar among groups, but the pattern of replacement and silent mutations in memory B cell clones indicated greater antigen selection in SR than CE. Greater clonal evolution of SR than CE memory B cells was revealed by analysis of phylogenetic trees and CDR3 lengths. Pauciclonality of the peripheral memory B cell population is a distinguishing feature of persons who spontaneously resolved an HCV infection. This finding, previously considered characteristic only of patients with HCV-associated lymphoproliferative disorders, suggests that the B cell clones potentially involved in clearance of the virus may also be those susceptible to abnormal expansion.

Immunologic aspects of monoclonal B-cell lymphocytosis

Monoclonal B-cell lymphocytosis (MBL) is a preclinical hematologic condition wherein small numbers of clonal B cells can be detected in the blood of otherwise healthy individuals. Most MBL have a surface immunophenotype nearly identical to that of chronic lymphocytic leukemia (CLL), though other phenotypes can also be identified. MBL has been shown to be a precursor state for CLL, but most MBL clones are quite small and apparently have minimal potential to progress of CLL or other B-cell lymphoproliferative disorder (B-LPD). The investigation of MBL as a precursor state for CLL will likely lead to important insights into mechanisms of disease pathogenesis. The review will cover clinical and translational aspects of MBL, with a particular emphasis on the prevalence of MBL; the relationship between MBL, CLL, and other B-LPDs; and the capacity of MBL to modulate the normal B- and T-cell compartments.

Splenic marginal zone lymphoma with VH1-02 gene rearrangement expresses poly- and self-reactive antibodies with similar reactivity

One-third of all splenic marginal zone lymphomas (SMZL) use the IgH VH1-02 gene. These cases are usually not associated with hepatitis C virus infection. Of interest, the rearranged VH1-02 genes display similar complementarity determining regions 3, a finding confirmed by our study. The latter suggests that these SMZL may produce antibodies with similar reactivity. We produced recombinant antibodies from 5 SMZL cases with VH1-02 gene rearrangement to study the binding reactivity of these antibodies. Surprisingly, the recombinant antibodies demonstrated poly- and self-reactivity as demonstrated by their reactivity with nuclear, cytoplasmic, as well as membranous antigens expressed by human cells and by reactivity with human serum. This polyreactivity was specific as demonstrated by ELISA. The antibodies did not react with proteins on the cell surface that are induced by apoptosis as shown for antibodies produced by chronic lymphatic leukemia with VH1-02 gene rearrangement. The results indicate that a common subset of SMZL arises from polyreactive B cells, a subset of marginal zone B cells that are important in the immunologic defense against infection.

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.

Human recombinant antibodies againstTrypanosoma cruziribosomal P2βprotein

Patients with chronic Chagas' Heart Disease (cChHD) develop an antibody response that is suspected to be involved in the cardiac pathogenesis. The response againstTrypanosoma cruziribosomal P proteins is of particular interest, as these antibodies can cross-react with host cardiac receptors causing electrophysiological alterations. To better understand the humoral anti-P response we constructed a single-chain variable fragment library derived from a cChHD patient. The variable heavy and light regions were amplified from bone-marrow RNA and subcloned into the vector pComb3X. The phage library was subsequently panned againstT. cruziribosomal P2βprotein (TcP2β). We obtained 3 different human recombinant antibodies that specifically reacted with TcP2βin ELISA and Western blots. Two of them reacted with the C-terminal region of TcP2β, peptide R13, as the recombinant autoanti-P antibodies from Systemic Lupus Erythematosus (SLE) patients. Interestingly, the third one was specific for TcP2βbut did not recognize R13, confirming the specific nature of the anti-P response in Chagas disease. Neither sequence nor VH usage similarities between Chagas and SLE anti-P autoantibodies were observed. Herein, the first human mAbs against TcP2βhave been obtained and characterized showing that the humoral anti-P response is directed against the parasite and does not include an autoimmune component.

Molecular characterization of the cervical and systemic B-cell repertoire: Unique, yet overlapping, immune compartments of an HIV-1 resistant individual

The cervical mucosa of women who are highly exposed to HIV-1, yet remain persistently seronegative (HEPS), presents a unique opportunity to study the dynamics of an immune compartment potentially capable of preventing HIV-1 infection. Herein, we provide a detailed characterization of the immunoglobulin repertoire of cervical and systemic B cells from one such HEPS individual from Nairobi, Kenya. Analysis was done on 512 VH sequences that were RT-PCR amplified from B cells in a paired sample from the cervix and peripheral blood. The VH3 and DH repertoire of class switched cervical B cells differs significantly from that of systemic B cells indicating that the cervical environment affects local B cell populations and hence VH gene expression. Six networks of clonally related, heavily mutated B cells were identified that spanned the systemic and cervical B cell compartments. Analysis of somatic mutations suggests this is likely the result of systemic, class switched B cells homing to the cervical mucosa. Multiple networks of somatically mutated V-gene sequences, unique to the cervical mucosa, were also identified. This supports the notion that site specific responses occur and have unique regulation of tolerance and recruitment into local memory or blast B cell compartments. We conclude that while the nature of the cervical environment shapes the local B cell repertoire, the infusion of post germinal center B cells to the human cervix is a common occurrence, and represents a means by which systemic immunization could provide the local antibodies necessary to prevent HIV-1 at the site of initial contact.

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.

Targets of somatic hypermutation within immunoglobulin light chain genes in zebrafish

In mammals, somatic hypermutation (SHM) of immunoglobulin (Ig) genes is critical for the generation of high-affinity antibodies and effective immune responses. Knowledge of sequence-specific biases in the targeting of somatic mutations can be useful for studies aimed at understanding antibody repertoires produced in response to infections, B-cell neoplasms, or autoimmune disease. To evaluate potential nucleotide targets of somatic mutation in zebrafish (Danio rerio), an enriched IgL cDNA library was constructed and > 250 randomly selected clones were sequenced and analysed. In total, 55 unique VJ-C sequences were identified encoding a total of 125 mutations. Mutations were most prevalent in V(L) with a bias towards single base transitions and increased mutation in the complementarity-determining regions (CDRs). Overall, mutations were overrepresented at WRCH/DGYW motifs suggestive of activation-induced cytidine deaminase (AID) targeting which is common in mice and humans. In contrast to mammalian models, N and P addition was not observed and mutations at AID hotspots were largely restricted to palindromic WRCH/DGYW motifs. Mutability indexes for di- and trinucleotide combinations confirmed C/G targets within WRCH/DGYW motifs to be statistically significant mutational hotspots and showed trinucleotides ATC and ATG to be mutation coldspots. Additive mutations in VJ-C sequences revealed patterns of clonal expansion consistent with affinity maturation responses seen in higher vertebrates. Taken together, the data reveal specific nucleotide targets of SHM in zebrafish and suggest that AID and affinity maturation contribute to antibody diversification in this emerging immunological model.

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.

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.

IgE transcripts in the circulation of allergic children reflect a classical antigen-driven B cell response and not a superantigen-like activation

Allergic asthma is the most frequent chronic disorder in childhood. Although IgE is a central effector molecule in allergic diseases, the nature of the IgE response is still under debate. The objective of our study was to clarify whether the IgE repertoire in the circulation of allergic children represents a classical Ag-driven and oligoclonal B cell response, a superantigen-like activation of a subset of B cells, or a polyclonal B-1 cell expansion. Using a highly sensitive RT-PCR method, we amplified, cloned, and sequenced IgE H chain transcripts from 13 children with allergic asthma. We gained 1366 functional IgE sequences, which currently represent the most extensive collection of human IgE transcripts. Compared to IgM transcripts from the same children, the somatic mutation rate was significantly enhanced in IgE transcripts (21 per thousand versus 72 per thousand; p < 0.001), which renders a polyclonal B-1 response unlikely. Moreover, IgE sequences displayed significantly enhanced Ag selection and hence were indicative of a classical Ag-driven immune response with affinity maturation (p < 0.001). In contrast to several recent studies, the usage pattern of variable gene segment of the H Ig chain in IgE transcripts followed the germline complexity, arguing against a superantigen-like interaction. We conclude that IgE transcripts in the circulation of children with allergic asthma reflect a classical adaptive B-2 cell response. This study provides reference data for a better characterization of the IgE response under immunomodulating therapies, such as anti-IgE therapy or allergen-specific immunotherapy.

Mathematical analysis of antigen selection in somatically mutated immunoglobulin genes associated with autoimmunity

Affinity maturation is a process by which low-affinity antibodies are transformed into highly specific antibodies in germinal centres. This process occurs by hypermutation of immunoglobulin heavy chain variable (IgH V) region genes followed by selection for high-affinity variants. It has been proposed that statistical tests can identify affinity maturation and antigen selection by analysing the frequency of replacement and silent mutations in the complementarity determining regions (CDRs) that contact antigen and the framework regions (FRs) that encode structural integrity. In this study three different methods that have been proposed for detecting selection: the binomial test, the multinomial test and the focused binomial test, have been assessed for their reliability and ability to detect selection in human IgH V genes. We observe first that no statistical test is able to identify selection in the CDR antigen-binding sites, second that tests can reliably detect selection in the FR and third that antibodies from nasal biopsies from patients with Wegener’s granulomatosis and pathogenic antibodies from systemic lupus erythematosus do not appear to be as stringently selected for structural integrity as other groups of functional sequences.

The heavy chain variable segment gene repertoire in chronic Chagas' heart disease

Patients chronically infected with Trypanosoma cruzi develop chronic Chagas' heart disease (cChHD). Their Ab response is suspected to be involved in the cardiac pathogenesis. Reactivity of serum Abs from these patients has been extensively studied but little is known about the diversity of the in vivo IgG repertoire. We analyzed 125 variable H chain (VH) genes and compared it to repertoires from healthy individuals, and patients with autoimmune processes and other infections. VH were from plasma cells isolated from heart tissue of three cChHD patients and from a Fab combinatorial library derived from bone marrow of another cChHD patient. The role of the parasite in shaping the Ab repertoire was assessed analyzing VH genes before and after panning against T. cruzi Ag. Among recovered VH genes, a significantly increased representation of VH4 was observed. Plasma cells at the site of cardiac infiltration showed an increased VH1 usage. CDR3 lengths were similar to the ones found in the healthy repertoire and significantly shorter than in other infections. VH derived from anti-T. cruzi Fab and plasma cells showed a higher proportion of hypermutated genes, 46.9% and 43.75%, respectively, vs 30.9% of the cChHD patient repertoire, pointing to the role of parasite Ags in the shaping of the humoral response in Chagas' disease. No histological evidence of germinal center-like structures was observed in heart tissue. In accordance, VH analysis of heart plasmocytes revealed no evidence of clonal B cell expansion, suggesting that they migrated into heart tissue from secondary lymphoid organs.

Post-transplant lymphoproliferative disorders: role of viral infection, genetic lesions and antigen stimulation in the pathogenesis of the disease

Post-transplant lymphoproliferative disorders (PTLD) are a life-threatening complication of solid organ transplantation or, more rarely, hematopoietic stem cell transplantation. The majority of PTLD is of B-cell origin and associated with Epstein–Barr virus (EBV) infection. PTLD generally display involvement of extranodal sites, aggressive histology and aggressive clinical behavior. The molecular pathogenesis of PTLD involves infection by oncogenic viruses, namely EBV, as well as genetic or epigenetic alterations of several cellular genes. At variance with lymphoma arising in immunocompetent hosts, whose genome is relatively stable, a fraction of PTLD are characterized by microsatellite instability as a consequence of defects in the DNA mismatch repair mechanism. Apart from microsatellite instability, molecular alterations of cellular genes recognized in PTLD include alterations of cMYC, BCL6, TP53, DNA hypermethylation, and aberrant somatic hypermutation of protooncogenes. The occurrence of IGV mutations in the overwhelming majority of PTLD documents that malignant transformation targets germinal centre (GC) B-cells and their descendants both in EBV–positive and EBV–negative cases. Analysis of phenotypic markers of B-cell histogenesis, namely BCL6, MUM1 and CD138, allows further distinction of PTLD histogenetic categories. PTLD expressing the BCL6+/MUM1+/-/CD138− profile reflect B-cells actively experiencing the GC reaction, and comprise diffuse large B-cell lymphoma (DLBCL) centroblastic and Burkitt lymphoma. PTLD expressing the BCL6−/MUM1+/CD138− phenotype putatively derive from B-cells that have concluded the GC reaction, and comprise the majority of polymorphic PTLD and a fraction of DLBCL immunoblastic. A third group of PTLD is reminiscent of post-GC and preterminally differentiated B-cells that show the BCL6−/MUM1+/CD138+ phenotype, and are morphologically represented by either polymorphic PTLD or DLBCL immunoblastic.

Immunoglobulin heavy chain variable region gene usage and mutational status of the leukemic B cells in Iranian patients with chronic lymphocytic leukemia

The mutational status of the immunoglobulin variable region heavy chain genes (IGHV) is an important prognostic marker in chronic lymphocytic leukemia (CLL). The data accumulated in the literature has largely been derived from studies conducted on Caucasian Western populations. Little is known about Asian CLL patients. In this study the IGHV genes usage and somatic hypermutation status have been investigated in 87 Iranian CLL patients. Based on a cut-off of 98% nucleotide sequence homology, 64.4% and 35.6% of the patients expressed mutated and unmutated IGHV genes, respectively, with most non-progressive patients being in the mutated group (35/44 vs 19/40; P = 0.009). Progression-free survival (PFS) and time to first treatment (TTFT) were significantly higher in our mutated and non-progressive patients compared to unmutated and progressive subtypes, respectively. The most frequently used IGHV gene was IGHV3-7 (12.6%) followed by IGHV3-30 (11.4%), IGHV3-48 (9.2%), IGHV4-39 (6.9%), and IGHV1-8 (6.9%) genes, which taken together comprised nearly half of the IGHV genes expressed in the Iranian CLL patients. Of the IGHV genes, IGHV3-7 was significantly over-represented in non-progressive compared to progressive CLL patients (P = 0.036), whereas IGHV1-69 and IGHV1-2 were expressed at a higher frequency in unmutated compared to mutated CLL patients (P < 0.03). Comparison of IGHV gene usage in our patients with that of Western CLL patients revealed significant differences in expression of IGHV1-69, IGHV3-7, IGHV3-21, and IGHV4-34 genes. Analysis of the IGHV third complementary determining region (HCDR3) sequences revealed a high frequency use of certain HCDR3 motifs, such as YYYGMDV, in our samples. These findings imply contribution of antigen selection and regional (ethnic/geographic) parameters in the leukomogenesis of CLL.