Individual variation in the germline Ig gene repertoire inferred from variable region gene rearrangements

Humanized mouse G6 anti-idiotypic monoclonal antibody has therapeutic potential against IGHV1-69 germline gene-based B-CLL

In 10-20% of the cases of chronic lymphocytic leukemia of B-cell phenotype (B-CLL), the IGHV1-69 germline is utilized as VH gene of the B cell receptor (BCR). Mouse G6 (MuG6) is an anti-idiotypic monoclonal antibody discovered in a screen against rheumatoid factors (RFs) that binds with high affinity to an idiotope expressed on the 51p1 alleles of IGHV1-69 germline gene encoded antibodies (G6-id(+)). The finding that unmutated IGHV1-69 encoded BCRs are frequently expressed on B-CLL cells provides an opportunity for anti-idiotype monoclonal antibody immunotherapy. In this study, we first showed that MuG6 can deplete B cells encoding IGHV1-69 BCRs using a novel humanized GTL mouse model. Next, we humanized MuG6 and demonstrated that the humanized antibodies (HuG6s), especially HuG6.3, displayed ∼2-fold higher binding affinity for G6-id(+) antibody compared to the parental MuG6. Additional studies showed that HuG6.3 was able to kill G6-id(+) BCR expressing cells and patient B-CLL cells through antibody-dependent cell-mediated cytotoxicity (ADCC) and complement-dependent cytotoxicity (CDC). Finally, both MuG6 and HuG6.3 mediate in vivo depletion of B-CLL cells in NSG mice. These data suggest that HuG6.3 may provide a new precision medicine to selectively kill IGHV1-69-encoding G6-id(+) B-CLL cells.

Multi-loci diagnosis of acute lymphoblastic leukaemia with high-throughput sequencing and bioinformatics analysis

High-throughput sequencing (HTS) is considered a technical revolution that has improved our knowledge of lymphoid and autoimmune diseases, changing our approach to leukaemia both at diagnosis and during follow-up. As part of an immunoglobulin/T cell receptor-based minimal residual disease (MRD) assessment of acute lymphoblastic leukaemia patients, we assessed the performance and feasibility of the replacement of the first steps of the approach based on DNA isolation and Sanger sequencing, using a HTS protocol combined with bioinformatics analysis and visualization using the Vidjil software. We prospectively analysed the diagnostic and relapse samples of 34 paediatric patients, thus identifying 125 leukaemic clones with recombinations on multiple loci (TRG, TRD, IGH and IGK), including Dd2/Dd3 and Intron/KDE rearrangements. Sequencing failures were halved (14% vs. 34%, P = 0.0007), enabling more patients to be monitored. Furthermore, more markers per patient could be monitored, reducing the probability of false negative MRD results. The whole analysis, from sample receipt to clinical validation, was shorter than our current diagnostic protocol, with equal resources. V(D)J recombination was successfully assigned by the software, even for unusual recombinations. This study emphasizes the progress that HTS with adapted bioinformatics tools can bring to the diagnosis of leukaemia patients.

Isolation of potent neutralizing antibodies from a survivor of the 2014 Ebola virus outbreak

Antibodies targeting the Ebola virus surface glycoprotein (EBOV GP) are implicated in protection against lethal disease, but the characteristics of the human antibody response to EBOV GP remain poorly understood. We isolated and characterized 349 GP-specific monoclonal antibodies (mAbs) from the peripheral B cells of a convalescent donor who survived the 2014 EBOV Zaire outbreak. Remarkably, 77% of the mAbs neutralize live EBOV, and several mAbs exhibit unprecedented potency. Structures of selected mAbs in complex with GP reveal a site of vulnerability located in the GP stalk region proximal to the viral membrane. Neutralizing antibodies targeting this site show potent therapeutic efficacy against lethal EBOV challenge in mice. The results provide a framework for the design of new EBOV vaccine candidates and immunotherapies.

High-Throughput DNA Sequencing Analysis of Antibody Repertoires

New high-throughput DNA sequencing (HTS) technologies developed in the past decade have begun to be applied to the study of the complex gene rearrangements that encode human antibodies. This article first reviews the genetic features of Ig loci and the HTS technologies that have been applied to human repertoire studies, then discusses key choices for experimental design and data analysis in these experiments and the insights gained in immunological and infectious disease studies with the use of these approaches.

Efficient construct of a large and functional scFv yeast display library derived from the ascites B cells of ovarian cancer patients by three-fragment transformation-associated recombination

Over the past decade, yeast display technology has emerged as a powerful tool for the isolation of high-affinity immunoglobulin fragments with potential utility as clinical diagnostic and therapeutic reagents. Despite significant refinement of the various methodologies underpinning library construction and selections, certain aspects remain challenging and process limiting. We have sought to significantly improve the robustness of the single-chain Fv (scFv) library construction step by overcoming the technical inefficiencies frequently encountered during the PCR-mediated assembly of scFvs from the discrete heavy and light V-domain repertoires. Using a novel primer set designed to provide maximum amplification coverage of the known germ-line V-domain repertoire, we have exploited the potential of the in vivo homologous gap-repair apparatus of Saccharomyces cerevisiae to assemble intact scFvs directly from co-transformed PBMC-derived VH, VL, and linearized vector component fragments. We have successfully applied this three-fragment assembly strategy to construct a large (>10(9)) scFv yeast display library from the ascites immune repertoire of ovarian cancer patients and validated the approach by applying FACS-based sorting to readily isolate scFvs that recognize various tumor marker antigens (TMAs). It is expected that this simplified construction method may find general utility, both for de novo scFv library construction and for subsequent combinatorial affinity maturation manipulations that require more than two fragments.

Multiple Myeloma Minimal Residual Disease

Assessment of minimal residual disease (MRD) is becoming standard diagnostic care for potentially curable neoplasms such as some acute leukemias as well as chronic myeloid and lymphocytic leukemia. Although multiple myeloma (MM) remains as an incurable disease, around half of the patients achieve complete remission (CR), and recent data suggests increasing rates of curability with "total-therapy-like" programs. This landscape is likely to be improved with the advent of new antibodies and small molecules. Therefore, conventional serological and morphological techniques have become suboptimal for sensitive evaluation of highly effective treatment strategies. Although, existing data suggests that MRD could be used as a biomarker to evaluate treatment efficacy, help on therapeutic decisions, and act as surrogate for overall survival, the role of MRD in MM is still a matter of extensive debate. Here, we review the different levels of remission used to define depth of response in MM and their clinical significance, as well as the prognostic value and unique characteristics of MRD detection using immunophenotypic, molecular, and imaging techniques. Key facts The higher efficacy of new treatment strategies for MM demand the incorporation of highly sensitive techniques to monitor treatment efficacy MRD could be used as a more potent surrogate biomarker for survival than standard CR We need to understand the pros and cons of the different MRD techniques The time has come to incorporate highly sensitive, cost-effective, readily available, and standardized MRD techniques into clinical trials to assess its role in therapeutic decisions.

DJ Pairing during VDJ Recombination Shows Positional Biases That Vary among Individuals with Differing IGHD Locus Immunogenotypes

Human IgH diversity is influenced by biases in the pairing of IGHD and IGHJ genes, but these biases have not been described in detail. We used high-throughput sequencing of VDJ rearrangements to explore DJ pairing biases in 29 individuals. It was possible to infer three contrasting IGHD-IGHJ haplotypes in nine of these individuals, and two of these haplotypes include deletion polymorphisms involving multiple contiguous IGHD genes. Therefore, we were able to explore how the underlying genetic makeup of the H chain locus influences the formation of particular DJ pairs. Analysis of nonproductive rearrangements demonstrates that 3' IGHD genes tend to pair preferentially with 5' IGHJ genes, whereas 5' IGHD genes pair preferentially with 3' IGHJ genes; the relationship between IGHD gene pairing frequencies and IGHD gene position is a near linear one for each IGHJ gene. However, striking differences are seen in individuals who carry deletion polymorphisms in the D locus. The absence of different blocks of IGHD genes leads to increases in the utilization frequencies of just a handful of genes, and these genes have no clear positional relationships to the deleted genes. This suggests that pairing frequencies may be influenced by additional complex positional relationships that perhaps arise from chromatin structure. In contrast to IGHD gene usage, IGHJ gene usage is unaffected by the IGHD gene-deletion polymorphisms. Such an outcome would be expected if the recombinase complex associates with an IGHJ gene before associating with an IGHD gene partner.

Immune Repertoire Profiling Reveals that Clonally Expanded B and T Cells Infiltrating Diseased Human Kidneys Can Also Be Tracked in Blood

Recent advances in high-throughput sequencing allow for the competitive analysis of the human B and T cell immune repertoire. In this study we compared Immunoglobulin and T cell receptor repertoires of lymphocytes found in kidney and blood samples of 10 patients with various renal diseases based on next-generation sequencing data. We used Biomed-2 primer panels and ImmunExplorer software to sequence, analyze and compare complementarity determining regions and V-(D)-J elements. While generally an individual's renal receptor repertoire is different from the repertoire present in blood, 94% (30/32) of the lymphocytes with clonal expansion in kidney can also be traced in blood however, not all of these clonotypes are equally abundant. Summarizing the data of all analyzed patients, 68% of highly expanded T cell clonotypes and 30% of the highly expanded B cell clonotypes that have infiltrated the kidney can be found amongst the five most abundant clonotypes in blood. In addition, complementarity determining region 3 sequences of the immunoglobulin heavy chains are on average more diverse than T cell receptor beta chains. Immune repertoire analysis of tissue infiltrating B and T cells adds new approaches to the assessment of adaptive immune response in kidney diseases. Our data suggest that expanded clonotypes in the tissues might be traceable in blood samples in the course of treatment or the natural history of the disease.

Practical guidelines for B-cell receptor repertoire sequencing analysis

High-throughput sequencing of B-cell immunoglobulin repertoires is increasingly being applied to gain insights into the adaptive immune response in healthy individuals and in those with a wide range of diseases. Recent applications include the study of autoimmunity, infection, allergy, cancer and aging. As sequencing technologies continue to improve, these repertoire sequencing experiments are producing ever larger datasets, with tens- to hundreds-of-millions of sequences. These data require specialized bioinformatics pipelines to be analyzed effectively. Numerous methods and tools have been developed to handle different steps of the analysis, and integrated software suites have recently been made available. However, the field has yet to converge on a standard pipeline for data processing and analysis. Common file formats for data sharing are also lacking. Here we provide a set of practical guidelines for B-cell receptor repertoire sequencing analysis, starting from raw sequencing reads and proceeding through pre-processing, determination of population structure, and analysis of repertoire properties. These include methods for unique molecular identifiers and sequencing error correction, V(D)J assignment and detection of novel alleles, clonal assignment, lineage tree construction, somatic hypermutation modeling, selection analysis, and analysis of stereotyped or convergent responses. The guidelines presented here highlight the major steps involved in the analysis of B-cell repertoire sequencing data, along with recommendations on how to avoid common pitfalls.

Persistence and evolution of allergen-specific IgE repertoires during subcutaneous specific immunotherapy

BACKGROUND

Specific immunotherapy (SIT) is the only treatment with proved long-term curative potential in patients with allergic disease. Allergen-specific IgE is the causative agent of allergic disease, and antibodies contribute to SIT, but the effects of SIT on aeroallergen-specific B-cell repertoires are not well understood.

OBJECTIVE

We sought to characterize the IgE sequences expressed by allergen-specific B cells and track the fate of these B-cell clones during SIT.

METHODS

We used high-throughput antibody gene sequencing and identification of allergen-specific IgE with combinatorial antibody fragment library technology to analyze immunoglobulin repertoires of blood and the nasal mucosa from aeroallergen-sensitized subjects before and during the first year of subcutaneous SIT.

RESULTS

Of 52 distinct allergen-specific IgE heavy chains from 8 allergic donors, 37 were also detected by using high-throughput antibody gene sequencing of blood samples, nasal mucosal samples, or both. The allergen-specific clones had increased persistence, higher likelihood of belonging to clones expressing other switched isotypes, and possibly larger clone size than the rest of the IgE repertoire. Clone members in nasal tissue showed close mutational relationships.

CONCLUSION

In the future, combining functional binding studies, deep antibody repertoire sequencing, and information on clinical outcomes in larger studies might aid assessment of SIT mechanisms and efficacy.

Survival of human lymphoma cells requires B-cell receptor engagement by self-antigens

The activated B-cell-like (ABC) subtype of diffuse large B-cell lymphoma (DLBCL) relies on chronic active B-cell receptor (BCR) signaling. BCR pathway inhibitors induce remissions in a subset of ABC DLBCL patients. BCR microclusters on the surface of ABC cells resemble those generated following antigen engagement of normal B cells. We speculated that binding of lymphoma BCRs to self-antigens initiates and maintains chronic active BCR signaling in ABC DLBCL. To assess whether antigenic engagement of the BCR is required for the ongoing survival of ABC cells, we developed isogenic ABC cells that differed solely with respect to the IgH V region of their BCRs. In competitive assays with wild-type cells, substitution of a heterologous V region impaired the survival of three ABC lines. The viability of one VH4-34(+) ABC line and the ability of its BCR to bind to its own cell surface depended on V region residues that mediate the intrinsic autoreactivity of VH4-34 to self-glycoproteins. The BCR of another ABC line reacted with self-antigens in apoptotic debris, and the survival of a third ABC line was sustained by reactivity of its BCR to an idiotypic epitope in its own V region. Hence, a diverse set of self-antigens is responsible for maintaining the malignant survival of ABC DLBCL cells. IgH V regions used by the BCRs of ABC DLBCL biopsy samples varied in their ability to sustain survival of these ABC lines, suggesting a screening procedure to identify patients who might benefit from BCR pathway inhibition.

LymAnalyzer: a tool for comprehensive analysis of next generation sequencing data of T cell receptors and immunoglobulins

The adaptive immune system includes populations of B and T cells capable of binding foreign epitopes via antigen specific receptors, called immunoglobulin (IG) for B cells and the T cell receptor (TCR) for T cells. In order to provide protection from a wide range of pathogens, these cells display highly diverse repertoires of IGs and TCRs. This is achieved through combinatorial rearrangement of multiple gene segments in addition, for B cells, to somatic hypermutation. Deep sequencing technologies have revolutionized analysis of the diversity of these repertoires; however, accurate TCR/IG diversity profiling requires specialist bioinformatics tools. Here we present LymAnalzyer, a software package that significantly improves the completeness and accuracy of TCR/IG profiling from deep sequence data and includes procedures to identify novel alleles of gene segments. On real and simulated data sets LymAnalyzer produces highly accurate and complete results. Although, to date we have applied it to TCR/IG data from human and mouse, it can be applied to data from any species for which an appropriate database of reference genes is available. Implemented in Java, it includes both a command line version and a graphical user interface and is freely available at https://sourceforge.net/projects/lymanalyzer/.

The evolution within us

The B-cell immune response is a remarkable evolutionary system found in jawed vertebrates. B-cell receptors, the membrane-bound form of antibodies, are capable of evolving high affinity to almost any foreign protein. High germline diversity and rapid evolution upon encounter with antigen explain the general adaptability of B-cell populations, but the dynamics of repertoires are less well understood. These dynamics are scientifically and clinically important. After highlighting the remarkable characteristics of naive and experienced B-cell repertoires, especially biased usage of genes encoding the B-cell receptors, we contrast methods of sequence analysis and their attempts to explain patterns of B-cell evolution. These phylogenetic approaches are currently unlinked to explicit models of B-cell competition, which analyse repertoire evolution at the level of phenotype, the affinities and specificities to particular antigenic sites. The models, in turn, suggest how chance, infection history and other factors contribute to different patterns of immunodominance and protection between people. Challenges in rational vaccine design, specifically vaccines to induce broadly neutralizing antibodies to HIV, underscore critical gaps in our understanding of B cells' evolutionary and ecological dynamics.

HIV-1 VACCINES. Diversion of HIV-1 vaccine-induced immunity by gp41-microbiota cross-reactive antibodies

An HIV-1 DNA prime vaccine, with a recombinant adenovirus type 5 (rAd5) boost, failed to protect from HIV-1 acquisition. We studied the nature of the vaccine-induced antibody (Ab) response to HIV-1 envelope (Env). HIV-1-reactive plasma Ab titers were higher to Env gp41 than to gp120, and repertoire analysis demonstrated that 93% of HIV-1-reactive Abs from memory B cells responded to Env gp41. Vaccine-induced gp41-reactive monoclonal antibodies were non-neutralizing and frequently polyreactive with host and environmental antigens, including intestinal microbiota (IM). Next-generation sequencing of an immunoglobulin heavy chain variable region repertoire before vaccination revealed an Env-IM cross-reactive Ab that was clonally related to a subsequent vaccine-induced gp41-reactive Ab. Thus, HIV-1 Env DNA-rAd5 vaccine induced a dominant IM-polyreactive, non-neutralizing gp41-reactive Ab repertoire response that was associated with no vaccine efficacy.

Recovery of a human natural antibody against the noncollagenous-1 domain of type IV collagen using humanized models

BACKGROUND

Anti-glomerular basement membrane nephritis and Goodpasture syndrome result from autoantibody (Ab)-mediated destruction of kidney and lung. Ab target the noncollagenous 1 (NC1) domain of alpha3(IV) collagen, but little is known about Ab origins or structure. This ignorance is due in part to the inability to recover monoclonal Ab by transformation of patients' blood cells. The aim of this study was to assess the suitability of two humanized models for this purpose.

METHODS

NOD-scid-gamma immunodeficient mice were engrafted either with human CD34+ hematopoietic stem cells (HSC) (Hu-HSC mice) and immunized with alpha3(IV)NC1 collagen containing the Goodpasture epitopes or with nephritis patients' peripheral blood leukocytes (PBL) (Hu-PBL mice). After in vivo immune cell development and/or expansion, recovered human B cells were Epstein Barr virus (EBV)-transformed, screened for antigen (Ag) binding, electrofused with a mouse-human heterohybridoma, subcloned, and human Ab RNA sequenced by PCR after reverse transcription to cDNA. Flow cytometry was used to assess human B cell markers and differentiation in Hu-PBL mice.

RESULTS

Sequence analysis of a human Ab derived from an immunized Hu-HSC mouse and reactive with alpha3(IV)NC1 collagen reveals that it is encoded by unmutated heavy and light chain genes. The heavy chain complementarity determining region 3, a major determinant of Ag binding, contains uncommon motifs, including an N-region somatically-introduced highly hydrophobic tetrapeptide and dual cysteines encoded by a uniquely human IGHD2-2 Ab gene segment that lacks a murine counterpart. Comparison of human and mouse autoantibodies suggests that structurally similar murine Ab may arise by convergent selection. In contrast to the Hu-HSC model, transformed human B cells are rarely recovered from Hu-PBL mice, in which human B cells terminally differentiate and lose expression of EBV receptor CD21, thus precluding their transformation and recovery.

CONCLUSIONS

Hu-HSC mice reveal that potentially pathogenic B cells bearing unmutated Ig receptors reactive with the NC1 domain on alpha3(IV) collagen can be generated in, and not purged from, the human preimmune repertoire. Uniquely human gene elements are recruited to generate the antigen binding site in at least a subset of these autoantibodies, indicating that humanized models may provide insights inaccessible using conventional mouse models.

HTJoinSolver: Human immunoglobulin VDJ partitioning using approximate dynamic programming constrained by conserved motifs

Background

Partitioning the human immunoglobulin variable region into variable (V), diversity (D), and joining (J) segments is a common sequence analysis step. We introduce a novel approximate dynamic programming method that uses conserved immunoglobulin gene motifs to improve performance of aligning V-segments of rearranged immunoglobulin (Ig) genes. Our new algorithm enhances the former JOINSOLVER algorithm by processing sequences with insertions and/or deletions (indels) and improves the efficiency for large datasets provided by high throughput sequencing.

Results

In our simulations, which include rearrangements with indels, the V-matching success rate improved from 61% for partial alignments of sequences with indels in the original algorithm to over 99% in the approximate algorithm. An improvement in the alignment of human VDJ rearrangements over the initial JOINSOLVER algorithm was also seen when compared to the Stanford.S22 human Ig dataset with an online VDJ partitioning software evaluation tool.

Conclusions

HTJoinSolver can rapidly identify V- and J-segments with indels to high accuracy for mutated sequences when the mutation probability is around 30% and 20% respectively. The D-segment is much harder to fit even at 20% mutation probability. For all segments, the probability of correctly matching V, D, and J increases with our alignment score.

Soma-to-germline feedback is implied by the extreme polymorphism at IGHV relative to MHC: The manifest polymorphism of the MHC appears greatly exceeded at Immunoglobulin loci, suggesting antigen-selected somatic V mutants penetrate Weismann's Barrier

Soma-to-germline feedback is forbidden under the neo-Darwinian paradigm. Nevertheless, there is a growing realization it occurs frequently in immunoglobulin (Ig) variable (V) region genes. This is a surprising development. It arises from a most unlikely source in light of the exposure of co-author EJS to the haplotype data of RL Dawkins and others on the polymorphism of the Major Histocompatibility Complex, which is generally assumed to be the most polymorphic region in the genome (spanning ∼4 Mb). The comparison between the magnitude of MHC polymorphism with estimates for the human heavy chain immunoglobulin V locus (spanning ∼1 Mb), suggests IGHV could be many orders of magnitude more polymorphic than the MHC. This conclusion needs airing in the literature as it implies generational churn and soma-to-germline gene feedback. Pedigree-based experimental strategies to resolve the IGHV issue are outlined.

Sequencing the functional antibody repertoire--diagnostic and therapeutic discovery

The development of high-throughput DNA sequencing technologies has enabled large-scale characterization of functional antibody repertoires, a new method of understanding protective and pathogenic immune responses. Important parameters to consider when sequencing antibody repertoires include the methodology, the B-cell population and clinical characteristics of the individuals analysed, and the bioinformatic analysis. Although focused sequencing of immunoglobulin heavy chains or complement determining regions can be utilized to monitor particular immune responses and B-cell malignancies, high-fidelity analysis of the full-length paired heavy and light chains expressed by individual B cells is critical for characterizing functional antibody repertoires. Bioinformatic identification of clonal antibody families and recombinant expression of representative members produces recombinant antibodies that can be used to identify the antigen targets of functional immune responses and to investigate the mechanisms of their protective or pathogenic functions. Integrated analysis of coexpressed functional genes provides the potential to further pinpoint the most important antibodies and clonal families generated during an immune response. Sequencing antibody repertoires is transforming our understanding of immune responses to autoimmunity, vaccination, infection and cancer. We anticipate that antibody repertoire sequencing will provide next-generation biomarkers, diagnostic tools and therapeutic antibodies for a spectrum of diseases, including rheumatic diseases.

Violation of an evolutionarily conserved immunoglobulin diversity gene sequence preference promotes production of dsDNA-specific IgG antibodies

Variability in the developing antibody repertoire is focused on the third complementarity determining region of the H chain (CDR-H3), which lies at the center of the antigen binding site where it often plays a decisive role in antigen binding. The power of VDJ recombination and N nucleotide addition has led to the common conception that the sequence of CDR-H3 is unrestricted in its variability and random in its composition. Under this view, the immune response is solely controlled by somatic positive and negative clonal selection mechanisms that act on individual B cells to promote production of protective antibodies and prevent the production of self-reactive antibodies. This concept of a repertoire of random antigen binding sites is inconsistent with the observation that diversity (D_H) gene segment sequence content by reading frame (RF) is evolutionarily conserved, creating biases in the prevalence and distribution of individual amino acids in CDR-H3. For example, arginine, which is often found in the CDR-H3 of dsDNA binding autoantibodies, is under-represented in the commonly used D_H RFs rearranged by deletion, but is a frequent component of rarely used inverted RF1 (iRF1), which is rearranged by inversion. To determine the effect of altering this germline bias in D_H gene segment sequence on autoantibody production, we generated mice that by genetic manipulation are forced to utilize an iRF1 sequence encoding two arginines. Over a one year period we collected serial serum samples from these unimmunized, specific pathogen-free mice and found that more than one-fifth of them contained elevated levels of dsDNA-binding IgG, but not IgM; whereas mice with a wild type D_H sequence did not. Thus, germline bias against the use of arginine enriched D_H sequence helps to reduce the likelihood of producing self-reactive antibodies.

Automated analysis of high-throughput B-cell sequencing data reveals a high frequency of novel immunoglobulin V gene segment alleles

Individual variation in germline and expressed B-cell immunoglobulin (Ig) repertoires has been associated with aging, disease susceptibility, and differential response to infection and vaccination. Repertoire properties can now be studied at large-scale through next-generation sequencing of rearranged Ig genes. Accurate analysis of these repertoire-sequencing (Rep-Seq) data requires identifying the germline variable (V), diversity (D), and joining (J) gene segments used by each Ig sequence. Current V(D)J assignment methods work by aligning sequences to a database of known germline V(D)J segment alleles. However, existing databases are likely to be incomplete and novel polymorphisms are hard to differentiate from the frequent occurrence of somatic hypermutations in Ig sequences. Here we develop a Tool for Ig Genotype Elucidation via Rep-Seq (TIgGER). TIgGER analyzes mutation patterns in Rep-Seq data to identify novel V segment alleles, and also constructs a personalized germline database containing the specific set of alleles carried by a subject. This information is then used to improve the initial V segment assignments from existing tools, like IMGT/HighV-QUEST. The application of TIgGER to Rep-Seq data from seven subjects identified 11 novel V segment alleles, including at least one in every subject examined. These novel alleles constituted 13% of the total number of unique alleles in these subjects, and impacted 3% of V(D)J segment assignments. These results reinforce the highly polymorphic nature of human Ig V genes, and suggest that many novel alleles remain to be discovered. The integration of TIgGER into Rep-Seq processing pipelines will increase the accuracy of V segment assignments, thus improving B-cell repertoire analyses.

Sequencing of the human IG light chain loci from a hydatidiform mole BAC library reveals locus-specific signatures of genetic diversity

Germline variation at immunoglobulin (IG) loci is critical for pathogen-mediated immunity, but establishing complete haplotype sequences in these regions has been problematic because of complex sequence architecture and diploid source DNA. We sequenced BAC clones from the effectively haploid human hydatidiform mole cell line, CHM1htert, across the light chain IG loci, kappa (IGK) and lambda (IGL), creating single haplotype representations of these regions. The IGL haplotype generated here is 1.25 Mb of contiguous sequence, including four novel IGLV alleles, one novel IGLC allele, and an 11.9-kb insertion. The CH17 IGK haplotype consists of two 644 kb proximal and 466 kb distal contigs separated by a large gap of unknown size; these assemblies added 49 kb of unique sequence extending into this gap. Our analysis also resulted in the characterization of seven novel IGKV alleles and a 16.7-kb region exhibiting signatures of interlocus sequence exchange between distal and proximal IGKV gene clusters. Genetic diversity in IGK/IGL was compared with that of the IG heavy chain (IGH) locus within the same haploid genome, revealing threefold (IGK) and sixfold (IGL) higher diversity in the IGH locus, potentially associated with increased levels of segmental duplication and the telomeric location of IGH.

HIV antibodies. Antigen modification regulates competition of broad and narrow neutralizing HIV antibodies

Some HIV-infected individuals develop broadly neutralizing antibodies (bNAbs), whereas most develop antibodies that neutralize only a narrow range of viruses (nNAbs). bNAbs, but not nNAbs, protect animals from experimental infection and are likely a key component of an effective vaccine. nNAbs and bNAbs target the same regions of the viral envelope glycoprotein (Env), but for reasons that remain unclear only nNAbs are elicited by Env immunization. We show that in contrast to germline-reverted (gl) bNAbs, glnNAbs recognized diverse recombinant Envs. Moreover, owing to binding affinity differences, nNAb B cell progenitors had an advantage in becoming activated and internalizing Env compared with bNAb B cell progenitors. We then identified an Env modification strategy that minimized the activation of nNAb B cells targeting epitopes that overlap those of bNAbs.

B-cell repertoire responses to varicella-zoster vaccination in human identical twins

Adaptive immune responses in humans rely on somatic genetic rearrangements of Ig and T-cell receptor loci to generate diverse antigen receptors. It is unclear to what extent an individual's genetic background affects the characteristics of the antibody repertoire used in responding to vaccination or infection. We studied the B-cell repertoires and clonal expansions in response to attenuated varicella-zoster vaccination in four pairs of adult identical twins and found that the global antibody repertoires of twin pair members showed high similarity in antibody heavy chain V, D, and J gene segment use, and in the length and features of the complementarity-determining region 3, a major determinant of antigen binding. These twin similarities were most pronounced in the IgM-expressing B-cell pools, but were seen to a lesser extent in IgG-expressing B cells. In addition, the degree of antibody somatic mutation accumulated in the B-cell repertoire was highly correlated within twin pair members. Twin pair members had greater numbers of shared convergent antibody sequences, including mutated sequences, suggesting similarity among memory B-cell clonal lineages. Despite these similarities in the memory repertoire, the B-cell clones used in acute responses to ZOSTAVAX vaccination were largely unique to each individual. Taken together, these results suggest that the overall B-cell repertoire is significantly shaped by the underlying germ-line genome, but that stochastic or individual-specific effects dominate the selection of clones in response to an acute antigenic stimulus.

Characterizing immune repertoires by high throughput sequencing: strategies and applications

As the key cellular effectors of adaptive immunity, T and B lymphocytes utilize specialized receptors to recognize, respond to, and neutralize a diverse array of extrinsic threats. These receptors (immunoglobulins in B lymphocytes, T cell receptors in T lymphocytes) are incredibly variable, the products of specialized genetic diversification mechanisms that generate complex lymphocyte repertoires with extensive collections of antigen specificities. Recent advances in high throughput sequencing (HTS) technologies have transformed our ability to examine antigen receptor repertoires at single nucleotide, and more recently, single cell, resolution. Here we review current approaches to examining antigen receptor repertoires by HTS, and discuss inherent biological and technical challenges. We further describe emerging applications of this powerful methodology for exploring the adaptive immune system.

Identification of antigen-specific B cell receptor sequences using public repertoire analysis

High-throughput sequencing allows detailed study of the BCR repertoire postimmunization, but it remains unclear to what extent the de novo identification of Ag-specific sequences from the total BCR repertoire is possible. A conjugate vaccine containing Haemophilus influenzae type b (Hib) and group C meningococcal polysaccharides, as well as tetanus toxoid (TT), was used to investigate the BCR repertoire of adult humans following immunization and to test the hypothesis that public or convergent repertoire analysis could identify Ag-specific sequences. A number of Ag-specific BCR sequences have been reported for Hib and TT, which made a vaccine containing these two Ags an ideal immunological stimulus. Analysis of identical CDR3 amino acid sequences that were shared by individuals in the postvaccine repertoire identified a number of known Hib-specific sequences but only one previously described TT sequence. The extension of this analysis to nonidentical, but highly similar, CDR3 amino acid sequences revealed a number of other TT-related sequences. The anti-Hib avidity index postvaccination strongly correlated with the relative frequency of Hib-specific sequences, indicating that the postvaccination public BCR repertoire may be related to more conventional measures of immunogenicity correlating with disease protection. Analysis of public BCR repertoire provided evidence of convergent BCR evolution in individuals exposed to the same Ags. If this finding is confirmed, the public repertoire could be used for rapid and direct identification of protective Ag-specific BCR sequences from peripheral blood.

Frequency Patterns of T-Cell Exposed Amino Acid Motifs in Immunoglobulin Heavy Chain Peptides Presented by MHCs

Immunoglobulins are highly diverse protein sequences that are processed and presented to T-cells by B-cells and other antigen presenting cells. We examined a large dataset of immunoglobulin heavy chain variable regions (IGHV) to assess the diversity of T-cell exposed motifs (TCEMs). TCEM comprise those amino acids in a MHC-bound peptide, which face outwards, surrounded by the MHC histotope, and which engage the T-cell receptor. Within IGHV there is a distinct pattern of predicted MHC class II binding and a very high frequency of re-use of the TCEMs. The re-use frequency indicates that only a limited number of different cognate T-cells are required to engage many different clonal B-cells. The amino acids in each outward-facing TCEM are intercalated with the amino acids of inward-facing MHC groove-exposed motifs (GEM). Different GEM may have differing, allele-specific, MHC binding affinities. The intercalation of TCEM and GEM in a peptide allows for a vast combinatorial repertoire of epitopes, each eliciting a different response. Outcome of T-cell receptor binding is determined by overall signal strength, which is a function of the number of responding T-cells and the duration of engagement. Hence, the frequency of TCEM re-use appears to be an important determinant of whether a T-cell response is stimulatory or suppressive. The frequency distribution of TCEMs implies that somatic hypermutation is followed by T-cell clonal expansion that develops along repeated pathways. The observations of TCEM and GEM derived from immunoglobulins suggest a relatively simple, yet powerful, mechanism to correlate T-cell polyspecificity, through re-use of TCEMs, with a very high degree of specificity achieved by combination with a diversity of GEMs. The frequency profile of TCEMs also points to an economical mechanism for maintaining T-cell memory, recall, and self-discrimination based on an endogenously generated profile of motifs.

Human IgE against the major allergen Bet v 1--defining an epitope with limited cross-reactivity between different PR-10 family proteins

Background

The interaction between IgE and allergen is a key event at the initiation of an allergic response, and its characteristics have substantial effects on the clinical manifestation. Despite this, the molecular details of the interaction between human IgE and the major birch allergen Bet v 1, one of the most potent tree allergens, still remain poorly investigated.

Objective

To isolate Bet v 1-specific human monoclonal IgE and characterize their interaction with the allergen.

Methods

Recombinant human IgE were isolated from a combinatorial antibody fragment library and their interaction with Bet v 1 assessed using various immunological assays. The structure of one such IgE in the single-chain fragment variable format was determined using X-ray crystallography.

Results

We present four novel Bet v 1-specific IgE, for one of which we solve the structure, all with their genetic origin in the IGHV5 germline gene, and demonstrate that they target two non-overlapping epitopes on the surface of Bet v 1, thereby fulfilling the basic criteria for FcεRI cross-linkage. We further define these epitopes and for one epitope pinpoint single amino acid residues important for the interaction with human IgE. This provides a potential explanation, at the molecular level, for the differences in recognition of isoforms of Bet v 1 and other allergens in the PR-10 protein family displayed by IgE targeting this epitope. Finally, we present the first high-resolution structure of a human allergen-specific IgE fragment in the single-chain fragment variable (scFv) format.

Conclusions and Clinical Relevance

We here display the usefulness of allergen-specific human monoclonal IgE as a tool in studies of the crucial molecular interaction taking place at the initiation of an allergic response. Such studies may aid us in development of better diagnostic tools and guide us in the development of new therapeutic compounds.

The antibody genetics of multiple sclerosis: comparing next-generation sequencing to sanger sequencing

We previously identified a distinct mutation pattern in the antibody genes of B cells isolated from cerebrospinal fluid (CSF) that can identify patients who have relapsing-remitting multiple sclerosis (RRMS) and patients with clinically isolated syndromes who will convert to RRMS. This antibody gene signature (AGS) was developed using Sanger sequencing of single B cells. While potentially helpful to patients, Sanger sequencing is not an assay that can be practically deployed in clinical settings. In order to provide AGS evaluations to patients as part of their diagnostic workup, we developed protocols to generate AGS scores using next-generation DNA sequencing (NGS) on CSF-derived cell pellets without the need to isolate single cells. This approach has the potential to increase the coverage of the B-cell population being analyzed, reduce the time needed to generate AGS scores, and may improve the overall performance of the AGS approach as a diagnostic test in the future. However, no investigations have focused on whether NGS-based repertoires will properly reflect antibody gene frequencies and somatic hypermutation patterns defined by Sanger sequencing. To address this issue, we isolated paired CSF samples from eight patients who either had MS or were at risk to develop MS. Here, we present data that antibody gene frequencies and somatic hypermutation patterns are similar in Sanger and NGS-based antibody repertoires from these paired CSF samples. In addition, AGS scores derived from the NGS database correctly identified the patients who initially had or subsequently converted to RRMS, with precision similar to that of the Sanger sequencing approach. Further investigation of the utility of the AGS in predicting conversion to MS using NGS-derived antibody repertoires in a larger cohort of patients is warranted.

Deep sequencing reveals clonal evolution patterns and mutation events associated with relapse in B-cell lymphomas

Background

Molecular mechanisms associated with frequent relapse of diffuse large B-cell lymphoma (DLBCL) are poorly defined. It is especially unclear how primary tumor clonal heterogeneity contributes to relapse. Here, we explore unique features of B-cell lymphomas - VDJ recombination and somatic hypermutation - to address this question.

Results

We performed high-throughput sequencing of rearranged VDJ junctions in 14 pairs of matched diagnosis-relapse tumors, among which 7 pairs were further characterized by exome sequencing. We identify two distinctive modes of clonal evolution of DLBCL relapse: an early-divergent mode in which clonally related diagnosis and relapse tumors diverged early and developed in parallel; and a late-divergent mode in which relapse tumors developed directly from diagnosis tumors with minor divergence. By examining mutation patterns in the context of phylogenetic information provided by VDJ junctions, we identified mutations in epigenetic modifiers such as KMT2D as potential early driving events in lymphomagenesis and immune escape alterations as relapse-associated events.

Conclusions

Altogether, our study for the first time provides important evidence that DLBCL relapse may result from multiple, distinct tumor evolutionary mechanisms, providing rationale for therapies for each mechanism. Moreover, this study highlights the urgent need to understand the driving roles of epigenetic modifier mutations in lymphomagenesis, and immune surveillance factor genetic lesions in relapse.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0432-0) contains supplementary material, which is available to authorized users.

Capturing needles in haystacks: a comparison of B-cell receptor sequencing methods

BACKGROUND

Deep-sequencing methods are rapidly developing in the field of B-cell receptor (BCR) and T-cell receptor (TCR) diversity. These promise to revolutionise our understanding of adaptive immune dynamics, identify novel antibodies, and allow monitoring of minimal residual disease. However, different methods for BCR and TCR enrichment and amplification have been proposed. Here we perform the first systematic comparison between different methods of enrichment, amplification and sequencing for generating BCR and TCR repertoires using large sample numbers.

RESULTS

Resampling from the same RNA or cDNA pool results in highly correlated and reproducible repertoires, but resampling low frequency clones leads to stochastic variance. Repertoires generated by different sequencing methods (454 Roche and Illumina MiSeq) and amplification methods (multiplex PCR, 5' Rapid amplification of cDNA ends (5'RACE), and RNA-capture) are highly correlated, and resulting IgHV gene frequencies between the different methods were not significantly different. Read length has an impact on captured repertoire structure, and ultimately full-length BCR sequences are most informative for repertoire analysis as diversity outside of the CDR is very useful for phylogenetic analysis. Additionally, we show RNA-based BCR repertoires are more informative than using DNA.

CONCLUSIONS

Repertoires generated by different sequencing and amplification methods are consistent, but we show that read lengths, depths and error profiles should be considered in experimental design, and multiple sampling approaches could be employed to minimise stochastic sampling variation. This detailed investigation of immune repertoire sequencing methods is essential for informing basic and clinical research.

Systematic characterization and comparative analysis of the rabbit immunoglobulin repertoire

Rabbits have been used extensively as a model system for the elucidation of the mechanism of immunoglobulin diversification and for the production of antibodies. We employed Next Generation Sequencing to analyze Ig germline V and J gene usage, CDR3 length and amino acid composition, and gene conversion frequencies within the functional (transcribed) IgG repertoire of the New Zealand white rabbit (Oryctolagus cuniculus). Several previously unannotated rabbit heavy chain variable (VH) and light chain variable (VL) germline elements were deduced bioinformatically using multidimensional scaling and k-means clustering methods. We estimated the gene conversion frequency in the rabbit at 23% of IgG sequences with a mean gene conversion tract length of 59±36 bp. Sequencing and gene conversion analysis of the chicken, human, and mouse repertoires revealed that gene conversion occurs much more extensively in the chicken (frequency 70%, tract length 79±57 bp), was observed to a small, yet statistically significant extent in humans, but was virtually absent in mice.

The evolution of clonality testing in the diagnosis and monitoring of hematological malignancies

Currently, distinguishing between benign and malignant lymphoid proliferations is based on a combination of clinical characteristics, cyto/histomorphology, immunophenotype and the identification of well-defined chromosomal aberrations. However, such diagnoses remain challenging in 10-15% of cases of lymphoproliferative disorders, and clonality assessments are often required to confirm diagnostic suspicions. In recent years, the development of new techniques for clonality detection has allowed researchers to better characterize, classify and monitor hematological neoplasms. In the past, clonality was primarily studied by performing Southern blotting analyses to characterize rearrangements in segments of the IG and TCR genes. Currently, the most commonly used method in the clinical molecular diagnostic laboratory is polymerase chain reaction (PCR), which is an extremely sensitive technique for detecting nucleic acids. This technique is rapid, accurate, specific, and sensitive, and it can be used to analyze small biopsies as well as formalin-fixed paraffin-embedded samples. These advantages make PCR-based approaches the current gold standard for IG/TCR clonality testing. Since the completion of the first human genome sequence, there has been a rapid development of technologies to facilitate high-throughput sequencing of DNA. These techniques have been applied to the deep characterization and classification of various diseases, patient stratification, and the monitoring of minimal residual disease. Furthermore, these novel approaches have the potential to significantly improve the sensitivity and cost of clonality assays and post-treatment monitoring of B- and T-cell malignancies. However, more studies will be required to demonstrate the utility, sensitivity, and benefits of these methods in order to warrant their adoption into clinical practice. In this review, recent developments in clonality testing are examined with an emphasis on highly sensitive systems for improving diagnostic workups and minimal residual disease assessments.

Tissue-specific expressed antibody variable gene repertoires

Recent developments in genetic technologies allow deep analysis of the sequence diversity of immune repertoires, but little work has been reported on the architecture of immune repertoires in mucosal tissues. Antibodies are the key to prevention of infections at the mucosal surface, but it is currently unclear whether the B cell repertoire at mucosal surfaces reflects the dominant antibodies found in the systemic compartment or whether mucosal tissues harbor unique repertoires. We examined the expressed antibody variable gene repertoires from 10 different human tissues using RNA samples derived from a large number of individuals. The results revealed that mucosal tissues such as stomach, intestine and lung possess unique antibody gene repertoires that differed substantially from those found in lymphoid tissues or peripheral blood. Mutation frequency analysis of mucosal tissue repertoires revealed that they were highly mutated, with little evidence for the presence of naïve B cells, in contrast to blood. Mucosal tissue repertoires possessed longer heavy chain complementarity determining region 3 loops than lymphoid tissue repertoires. We also noted a large increase in frequency of both insertions and deletions in the small intestine antibody repertoire. These data suggest that mucosal immune repertoires are distinct in many ways from the systemic compartment.

Fast multiclonal clusterization of V(D)J recombinations from high-throughput sequencing

BACKGROUND

V(D)J recombinations in lymphocytes are essential for immunological diversity. They are also useful markers of pathologies. In leukemia, they are used to quantify the minimal residual disease during patient follow-up. However, the full breadth of lymphocyte diversity is not fully understood.

RESULTS

We propose new algorithms that process high-throughput sequencing (HTS) data to extract unnamed V(D)J junctions and gather them into clones for quantification. This analysis is based on a seed heuristic and is fast and scalable because in the first phase, no alignment is performed with germline database sequences. The algorithms were applied to TR γ HTS data from a patient with acute lymphoblastic leukemia, and also on data simulating hypermutations. Our methods identified the main clone, as well as additional clones that were not identified with standard protocols.

CONCLUSIONS

The proposed algorithms provide new insight into the analysis of high-throughput sequencing data for leukemia, and also to the quantitative assessment of any immunological profile. The methods described here are implemented in a C++ open-source program called Vidjil.

Studying the antibody repertoire after vaccination: practical applications

Nearly all licensed vaccines have been developed to confer protection against infectious diseases by stimulating the production of antibodies by B cells, but the nature of a successful antibody response has been difficult to capture. Recent advances in next-generation sequencing (NGS) technology have allowed high-resolution characterization of the antibody repertoire, and of the changes that occur following vaccination. These approaches have yielded important insights into the B cell response, and have raised the possibility of using specific antibody sequences as measures of vaccine immunogenicity. Here, we review recent findings based on antibody repertoire sequencing, and discuss potential applications of these new technologies and of the analyses of the increasing volume of antibody sequence data in the context of vaccine development.

Molecular signatures of hemagglutinin stem-directed heterosubtypic human neutralizing antibodies against influenza A viruses

Recent studies have shown high usage of the IGHV1-69 germline immunoglobulin gene for influenza hemagglutinin stem-directed broadly-neutralizing antibodies (HV1-69-sBnAbs). Here we show that a major structural solution for these HV1-69-sBnAbs is achieved through a critical triad comprising two CDR-H2 loop anchor residues (a hydrophobic residue at position 53 (Ile or Met) and Phe54), and CDR-H3-Tyr at positions 98±1; together with distinctive V-segment CDR amino acid substitutions that occur in positions sparse in AID/polymerase-η recognition motifs. A semi-synthetic IGHV1-69 phage-display library screen designed to investigate AID/polη restrictions resulted in the isolation of HV1-69-sBnAbs that featured a distinctive Ile52Ser mutation in the CDR-H2 loop, a universal CDR-H3 Tyr at position 98 or 99, and required as little as two additional substitutions for heterosubtypic neutralizing activity. The functional importance of the Ile52Ser mutation was confirmed by mutagenesis and by BCR studies. Structural modeling suggests that substitution of a small amino acid at position 52 (or 52a) facilitates the insertion of CDR-H2 Phe54 and CDR-H3-Tyr into adjacent pockets on the stem. These results support the concept that activation and expansion of a defined subset of IGHV1-69-encoded B cells to produce potent HV1-69-sBnAbs does not necessarily require a heavily diversified V-segment acquired through recycling/reentry into the germinal center; rather, the incorporation of distinctive amino acid substitutions by Phase 2 long-patch error-prone repair of AID-induced mutations or by random non-AID SHM events may be sufficient. We propose that these routes of B cell maturation should be further investigated and exploited as a pathway for HV1-69-sBnAb elicitation by vaccination.

The quest for universal influenza vaccine has gained wide interest with the discovery of human neutralizing antibodies that are able to variably cross neutralize and protect against different influenza strains, subtypes, groups and lineages. These antibodies, which bind to a highly conserved epitope in the hemagglutinin stem, are often encoded by rearranged IGHV1-69 germline genes that alone make contact with HA and prevent virus entry and emergence of escape mutants. Our study was undertaken to gain an understanding of what structural requirements enable a rearranged IGHV1-69 Ab to become a potent cross-neutralizing antibody. We found that in addition to a critical amino acid triad consisting of a pair of anchor residues in CDR-H2 and a properly positioned CDR-H3 Tyr, distinctive V-segment substitutions that arise in positions that are distinct from phase I AID somatic hypermutation (SHM) hotspot motifs are often required. As few as two V-segment SHM can fulfill this role which appears to facilitate the optimal binding of CDR-H2 Phe54 and CHR-H3-Tyr into adjacent hydrophobic pockets in the HA stem. These studies provide new information on the SHM requirements for IGHV1-69-encoded B cells to produce HV1-69-sBnAbs and suggest that there may exist alternative routes to their elicitation by vaccination.

Identification of human neutralizing antibodies against MERS-CoV and their role in virus adaptive evolution

The newly emerging Middle East Respiratory Syndrome coronavirus (MERS-CoV) causes a Severe Acute Respiratory Syndrome-like disease with ∼43% mortality. Given the recent detection of virus in dromedary camels, zoonotic transfer of MERS-CoV to humans is suspected. In addition, little is known about the role of human neutralizing Ab (nAb) pressure as a driving force in MERS-CoV adaptive evolution. Here, we used a well-characterized nonimmune human Ab-phage library and a panning strategy with proteoliposomes and cells to identify seven human nAbs against the receptor-binding domain (RBD) of the MERS-CoV Spike protein. These nAbs bind to three different epitopes in the RBD and human dipeptidyl peptidase 4 (hDPP4) interface with subnanomolar/nanomolar binding affinities and block the binding of MERS-CoV Spike protein with its hDPP4 receptor. Escape mutant assays identified five amino acid residues that are critical for neutralization escape. Despite the close proximity of the three epitopes on the RBD interface, escape from one epitope did not have a major impact on neutralization with Abs directed to a different epitope. Importantly, the majority of escape mutations had negative impacts on hDPP4 receptor binding and viral fitness. To our knowledge, these results provide the first report on human nAbs against MERS-CoV that may contribute to MERS-CoV clearance and evolution. Moreover, in the absence of a licensed vaccine or antiviral for MERS, this panel of nAbs offers the possibility of developing human mAb-based immunotherapy, especially for health-care workers.

Differences in the composition of the human antibody repertoire by B cell subsets in the blood

The vast initial diversity of the antibody repertoire is generated centrally by means of a complex series of V(D)J gene rearrangement events, variation in the site of gene segment joining, and TdT catalyzed N-region addition. Although the diversity is great, close inspection has revealed distinct and unique characteristics in the antibody repertoires expressed by different B cell developmental subsets. In order to illustrate our approach to repertoire analysis, we present an in-depth comparison of V(D)J gene usage, hydrophobicity, length, D_H reading frame, and amino acid usage between heavy chain repertoires expressed by immature, transitional, mature, memory IgD⁺, memory IgD⁻, and plasmacytes isolated from the blood of a single individual. Our results support the view that in both human and mouse, the H chain repertoires expressed by individual, developmental B cell subsets appear to differ in sequence content. Sequencing of unsorted B cells from the blood is thus likely to yield an incomplete or compressed view of what is actually happening in the immune response of the individual. Our findings support the view that studies designed to correlate repertoire expression with diseases of immune function will likely require deep sequencing of B cells sorted by subset.

IgE-associated IGHV genes from venom and peanut allergic individuals lack mutational evidence of antigen selection

Antigen selection of B cells within the germinal center reaction generally leads to the accumulation of replacement mutations in the complementarity-determining regions (CDRs) of immunoglobulin genes. Studies of mutations in IgE-associated VDJ gene sequences have cast doubt on the role of antigen selection in the evolution of the human IgE response, and it may be that selection for high affinity antibodies is a feature of some but not all allergic diseases. The severity of IgE-mediated anaphylaxis is such that it could result from higher affinity IgE antibodies. We therefore investigated IGHV mutations in IgE-associated sequences derived from ten individuals with a history of anaphylactic reactions to bee or wasp venom or peanut allergens. IgG sequences, which more certainly experience antigen selection, served as a control dataset. A total of 6025 unique IgE and 5396 unique IgG sequences were generated using high throughput 454 pyrosequencing. The proportion of replacement mutations seen in the CDRs of the IgG dataset was significantly higher than that of the IgE dataset, and the IgE sequences showed little evidence of antigen selection. To exclude the possibility that 454 errors had compromised analysis, rigorous filtering of the datasets led to datasets of 90 core IgE sequences and 411 IgG sequences. These sequences were present as both forward and reverse reads, and so were most unlikely to include sequencing errors. The filtered datasets confirmed that antigen selection plays a greater role in the evolution of IgG sequences than of IgE sequences derived from the study participants.

Qualifying high-throughput immune repertoire sequencing

Diversity of B and T cell receptors, achieved by gene recombination and somatic hypermutation, allows the immune system for recognition and targeted reaction against various threats. Next-generation sequencing for assessment of a cell's gene composition and variation makes deep analysis of one individual's immune spectrum feasible. An easy to apply but detailed analysis and visualization strategy is necessary to process all sequences generated. We performed sequencing utilizing the 454 system for CLL and control samples, utilized the IMGT database and applied the presented analysis tools. With the applied protocol, malignant clones are found and characterized, mutational status compared to germline identity is elaborated in detail showing that the CLL mutation status is not as monoclonal as generally thought. On the other hand, this strategy is not solely applicable to the 454 sequencing system but can easily be transferred to any other next-generation sequencing platform.

Heavy-light chain interrelations of MS-associated immunoglobulins probed by deep sequencing and rational variation

The mechanisms triggering most of autoimmune diseases are still obscure. Autoreactive B cells play a crucial role in the development of such pathologies and, in particular, production of autoantibodies of different specificities. The combination of deep-sequencing technology with functional studies of antibodies selected from highly representative immunoglobulin combinatorial libraries may provide unique information on specific features in the repertoires of autoreactive B cells. Here, we have analyzed cross-combinations of the variable regions of human immunoglobulins against the myelin basic protein (MBP) previously selected from a multiple sclerosis (MS)-related scFv phage-display library. On the other hand, we have performed deep sequencing of the sublibraries of scFvs against MBP, Epstein-Barr virus (EBV) latent membrane protein 1 (LMP1), and myelin oligodendrocyte glycoprotein (MOG). Bioinformatics analysis of sequencing data and surface plasmon resonance (SPR) studies have shown that it is the variable fragments of antibody heavy chains that mainly determine both the affinity of antibodies to the parent autoantigen and their cross-reactivity. It is suggested that LMP1-cross-reactive anti-myelin autoantibodies contain heavy chains encoded by certain germline gene segments, which may be a hallmark of the EBV-specific B cell subpopulation involved in MS triggering.

The promise and challenge of high-throughput sequencing of the antibody repertoire

Georgiou and colleagues discuss rapidly evolving methods for high-throughput sequencing of the antibody repertoire, and how the resulting data may be applied to answer basic and translational research questions.

Efforts to determine the antibody repertoire encoded by B cells in the blood or lymphoid organs using high-throughput DNA sequencing technologies have been advancing at an extremely rapid pace and are transforming our understanding of humoral immune responses. Information gained from high-throughput DNA sequencing of immunoglobulin genes (Ig-seq) can be applied to detect B-cell malignancies with high sensitivity, to discover antibodies specific for antigens of interest, to guide vaccine development and to understand autoimmunity. Rapid progress in the development of experimental protocols and informatics analysis tools is helping to reduce sequencing artifacts, to achieve more precise quantification of clonal diversity and to extract the most pertinent biological information. That said, broader application of Ig-seq, especially in clinical settings, will require the development of a standardized experimental design framework that will enable the sharing and meta-analysis of sequencing data generated by different laboratories.

Convergent antibody signatures in human dengue

Dengue is the most prevalent mosquito-borne viral disease in humans, and the lack of early prognostics, vaccines, and therapeutics contributes to immense disease burden. To identify patterns that could be used for sequence-based monitoring of the antibody response to dengue, we examined antibody heavy-chain gene rearrangements in longitudinal peripheral blood samples from 60 dengue patients. Comparing signatures between acute dengue, postrecovery, and healthy samples, we found increased expansion of B cell clones in acute dengue patients, with higher overall clonality in secondary infection. Additionally, we observed consistent antibody sequence features in acute dengue in the highly variable major antigen-binding determinant, complementarity-determining region 3 (CDR3), with specific CDR3 sequences highly enriched in acute samples compared to postrecovery, healthy, or non-dengue samples. Dengue thus provides a striking example of a human viral infection where convergent immune signatures can be identified in multiple individuals. Such signatures could facilitate surveillance of immunological memory in communities.

Next generation sequencing: new tools in immunology and hematology

One of the hallmarks of the adaptive immune system is the specificity of B and T cell receptors. Thanks to somatic recombination, a large repertoire of receptors can be generated within an individual that guarantee the recognition of a vast number of antigens. Monoclonal antibodies have limited applicability, given the high degree of diversity among these receptors, in BCR and TCR monitoring. Furthermore, with regard to cancer, better characterization of complex genomes and the ability to monitor tumor-specific cryptic mutations or translocations are needed to develop better tailored therapies. Novel technologies, by enhancing the ability of BCR and TCR monitoring, can help in the search for minimal residual disease during hematological malignancy diagnosis and follow-up, and can aid in improving bone marrow transplantation techniques. Recently, a novel technology known as next generation sequencing has been developed; this allows the recognition of unique sequences and provides depth of coverage, heterogeneity, and accuracy of sequencing. This provides a powerful tool that, along with microarray analysis for gene expression, may become integral in resolving the remaining key problems in hematology. This review describes the state of the art of this novel technology, its application in the immunological and hematological fields, and the possible benefits it will provide for the hematology and immunology community.

Effects of aging, cytomegalovirus infection, and EBV infection on human B cell repertoires

Elderly humans show decreased humoral immunity to pathogens and vaccines, yet the effects of aging on B cells are not fully known. Chronic viral infection by CMV is implicated as a driver of clonal T cell proliferations in some aging humans, but whether CMV or EBV infection contributes to alterations in the B cell repertoire with age is unclear. We have used high-throughput DNA sequencing of IGH gene rearrangements to study the BCR repertoires over two successive years in 27 individuals ranging in age from 20 to 89 y. Some features of the B cell repertoire remain stable with age, but elderly subjects show increased numbers of B cells with long CDR3 regions, a trend toward accumulation of more highly mutated IgM and IgG Ig genes, and persistent clonal B cell populations in the blood. Seropositivity for CMV or EBV infection alters B cell repertoires, regardless of the individual's age: EBV infection correlates with the presence of persistent clonal B cell expansions, whereas CMV infection correlates with the proportion of highly mutated Ab genes. These findings isolate effects of aging from those of chronic viral infection on B cell repertoires and provide a baseline for understanding human B cell responses to vaccination or infectious stimuli.

The human IgE repertoire

IgE is a key mediator in allergic diseases. However, in strong contrast to other antibody isotypes, many details of the composition of the human IgE repertoire are poorly defined. The low levels of human IgE in the circulation and the rarity of IgE-producing B cells are important reasons for this lack of knowledge. In this review, we summarize the current knowledge on these repertoires both in terms of their complexity and activity, i.e. knowledge which despite the difficulties encountered when studying the molecular details of human IgE has been acquired in recent years. We also take a look at likely future developments, for instance through improvements in sequencing technology and methodology that allow the isolation of additional allergen-specific human antibodies mimicking IgE, as this certainly will support our understanding of human IgE in the context of human disease in the years to come.

Genomic and systems approaches to translational biomarker discovery in immunological diseases

The high failure rate of new therapeutic mechanisms tested in clinical development has spurred an upsurge in research dedicated to discovering biomarker readouts that can improve decision-making. Increasingly, systems biology and genomic technologies, such as transcriptional profiling, are being leveraged to aid in the discovery of biomarker readouts. For inflammatory and immunological diseases, such as rheumatoid arthritis (RA) and asthma, progress has been made in developing biomarkers to monitor disease activity, prediction of response to therapy, and pharmacodynamic (PD) measurements. In this review, we discuss recent successes and challenges in these endeavors, highlighting the importance of human clinical studies of standard-of-care treatments in control subjects and patients with disease as the most direct path toward identifying useful translational biomarkers for clinical development.