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

The application of next generation sequencing to the understanding of antibody repertoires

In the decade since the human genome sequence was declared complete, the development of next generation sequencing (NGS) or "deep" sequencing to deliver cost-effective genomic sequencing has influenced advances beyond its primary application and changed the research landscape in many other areas. This review will survey recent applications of NGS which have broadened the understanding of natural antibody repertoires (the "antibodyome") and how these evolve in response to viral infection. We will also report examples where deep sequencing of binding populations, derived from both natural and synthetic repertoires, have been used to benefit antibody engineering. This knowledge will ultimately lead to the design of more effective biological drugs and vaccines.

The shape of the lymphocyte receptor repertoire: lessons from the B cell receptor

Both the B cell receptor (BCR) and the T cell receptor (TCR) repertoires are generated through essentially identical processes of V(D)J recombination, exonuclease trimming of germline genes, and the random addition of non-template encoded nucleotides. The naïve TCR repertoire is constrained by thymic selection, and TCR repertoire studies have therefore focused strongly on the diversity of MHC-binding complementarity determining region (CDR) CDR3. The process of somatic point mutations has given B cell studies a major focus on variable (IGHV, IGLV, and IGKV) genes. This in turn has influenced how both the naïve and memory BCR repertoires have been studied. Diversity (D) genes are also more easily identified in BCR VDJ rearrangements than in TCR VDJ rearrangements, and this has allowed the processes and elements that contribute to the incredible diversity of the immunoglobulin heavy chain CDR3 to be analyzed in detail. This diversity can be contrasted with that of the light chain where a small number of polypeptide sequences dominate the repertoire. Biases in the use of different germline genes, in gene processing, and in the addition of non-template encoded nucleotides appear to be intrinsic to the recombination process, imparting "shape" to the repertoire of rearranged genes as a result of differences spanning many orders of magnitude in the probabilities that different BCRs will be generated. This may function to increase the precursor frequency of naïve B cells with important specificities, and the likely emergence of such B cell lineages upon antigen exposure is discussed with reference to public and private T cell clonotypes.

Automated analysis of immunoglobulin genes from high-throughput sequencing: life without a template

Background

Immunoglobulin (that is, antibody) and T cell receptor genes are created through somatic gene rearrangement from gene segment libraries. Immunoglobulin genes are further diversified by somatic hypermutation and selection during the immune response. Studying the repertoires of these genes yields valuable insights into immune system function in infections, aging, autoimmune diseases and cancers. The introduction of high throughput sequencing has generated unprecedented amounts of repertoire and mutation data from immunoglobulin genes. However, common analysis programs are not appropriate for pre-processing and analyzing these data due to the lack of a template or reference for the whole gene.

Results

We present here the automated analysis pipeline we created for this purpose, which integrates various software packages of our own development and others’, and demonstrate its performance.

Conclusions

Our analysis pipeline presented here is highly modular, and makes it possible to analyze the data resulting from high-throughput sequencing of immunoglobulin genes, in spite of the lack of a template gene. An executable version of the Automation program (and its source code) is freely available for downloading from our website: http://immsilico2.lnx.biu.ac.il/Software.html.

Network properties derived from deep sequencing of human B-cell receptor repertoires delineate B-cell populations

The adaptive immune response selectively expands B- and T-cell clones following antigen recognition by B- and T-cell receptors (BCR and TCR), respectively. Next-generation sequencing is a powerful tool for dissecting the BCR and TCR populations at high resolution, but robust computational analyses are required to interpret such sequencing. Here, we develop a novel computational approach for BCR repertoire analysis using established next-generation sequencing methods coupled with network construction and population analysis. BCR sequences organize into networks based on sequence diversity, with differences in network connectivity clearly distinguishing between diverse repertoires of healthy individuals and clonally expanded repertoires from individuals with chronic lymphocytic leukemia (CLL) and other clonal blood disorders. Network population measures defined by the Gini Index and cluster sizes quantify the BCR clonality status and are robust to sampling and sequencing depths. BCR network analysis therefore allows the direct and quantifiable comparison of BCR repertoires between samples and intra-individual population changes between temporal or spatially separated samples and over the course of therapy.

Conserved variation: identifying patterns of stability and variability in BCR and TCR V genes with different diversity and richness metrics

The immune system can detect most invading pathogens. The potential for detection of pathogens is dependent on the somatic diversity of the immune repertoires. While it is known that this somatic diversity is carefully generated, it is unclear how the diversity is distributed in the different genes encoding receptors of immune cells. Utilizing different metrics for richness and diversity at the level of small sequence fragments, we present here an analysis of the entire known human germline repertoire as represented by the sequences from the ImMunoGeneTics database of immune receptors. We have developed a fragment sequence quantification analysis to track variation of repertoires with different degrees of precision. Somatic diversity has previously been functionally characterized mostly by division of the V gene sequences into the more conserved and invariant framework (FR) of the receptor and more varied complementarity determining regions (CDR), that interact with the antigen. We find that CDR and FR can be explicitly identified with our sequence fragment diversity quantification technique. In terms of diversity, CDR and FR are especially distinct in B cell V genes. T cell V genes show less of the CDR/FR periodicity but are more diverse overall. Our analysis further shows that there are other areas of diversity outside the CDR and FR that are found widely dispersed in T cell receptor V genes and more tightly focused in FR1 and FR3 in the B cell receptor V genes. The diversity we observe is not dependent on allelic differences nor is this diversity segregated by individual V gene families. We would thus expect that each individual exhibit a diversity equivalent to that of the entire potential repertoire.

Complete haplotype sequence of the human immunoglobulin heavy-chain variable, diversity, and joining genes and characterization of allelic and copy-number variation

The immunoglobulin heavy-chain locus (IGH) encodes variable (IGHV), diversity (IGHD), joining (IGHJ), and constant (IGHC) genes and is responsible for antibody heavy-chain biosynthesis, which is vital to the adaptive immune response. Programmed V-(D)-J somatic rearrangement and the complex duplicated nature of the locus have impeded attempts to reconcile its genomic organization based on traditional B-lymphocyte derived genetic material. As a result, sequence descriptions of germline variation within IGHV are lacking, haplotype inference using traditional linkage disequilibrium methods has been difficult, and the human genome reference assembly is missing several expressed IGHV genes. By using a hydatidiform mole BAC clone resource, we present the most complete haplotype of IGHV, IGHD, and IGHJ gene regions derived from a single chromosome, representing an alternate assembly of ∼1 Mbp of high-quality finished sequence. From this we add 101 kbp of previously uncharacterized sequence, including functional IGHV genes, and characterize four large germline copy-number variants (CNVs). In addition to this germline reference, we identify and characterize eight CNV-containing haplotypes from a panel of nine diploid genomes of diverse ethnic origin, discovering previously unmapped IGHV genes and an additional 121 kbp of insertion sequence. We genotype four of these CNVs by using PCR in 425 individuals from nine human populations. We find that all four are highly polymorphic and show considerable evidence of stratification (Fst = 0.3-0.5), with the greatest differences observed between African and Asian populations. These CNVs exhibit weak linkage disequilibrium with SNPs from two commercial arrays in most of the populations tested.

JEC. Secondary mechanisms of diversification in the human antibody repertoire

V(D)J recombination and somatic hypermutation (SHM) are the primary mechanisms for diversification of the human antibody repertoire. These mechanisms allow for rapid humoral immune responses to a wide range of pathogenic challenges. V(D)J recombination efficiently generate a virtually limitless diversity through random recombination of variable (V), diversity (D), and joining (J) genes with diverse non-templated junctions between the selected gene segments. Following antigen stimulation, affinity maturation by SHM produces antibodies with refined specificity mediated by mutations typically focused in complementarity determining regions (CDRs), which form the bulk of the antigen recognition site. While V(D)J recombination and SHM are responsible for much of the diversity of the antibody repertoire, there are several secondary mechanisms that, while less frequent, make substantial contributions to antibody diversity including V(DD)J recombination (or D-D fusion), SHM-associated insertions and deletions, and affinity maturation and antigen contact by non-CDR regions of the antibody. In addition to enhanced diversity, these mechanisms allow the production of antibodies that are critical to response to a variety of viral and bacterial pathogens but that would be difficult to generate using only the primary mechanisms of diversification.

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.

Automated cleaning and pre-processing of immunoglobulin gene sequences from high-throughput sequencing

High-throughput sequencing (HTS) yields tens of thousands to millions of sequences that require a large amount of pre-processing work to clean various artifacts. Such cleaning cannot be performed manually. Existing programs are not suitable for immunoglobulin (Ig) genes, which are variable and often highly mutated. This paper describes Ig High-Throughput Sequencing Cleaner (Ig-HTS-Cleaner), a program containing a simple cleaning procedure that successfully deals with pre-processing of Ig sequences derived from HTS, and Ig Insertion—Deletion Identifier (Ig-Indel-Identifier), a program for identifying legitimate and artifact insertions and/or deletions (indels). Our programs were designed for analyzing Ig gene sequences obtained by 454 sequencing, but they are applicable to all types of sequences and sequencing platforms. Ig-HTS-Cleaner and Ig-Indel-Identifier have been implemented in Java and saved as executable JAR files, supported on Linux and MS Windows. No special requirements are needed in order to run the programs, except for correctly constructing the input files as explained in the text. The programs' performance has been tested and validated on real and simulated data sets.

Characterization of human IgG repertoires in an acute HIV-1 infection

All known broadly neutralizing antibodies (bnAbs) are highly somatically mutated and therefore significantly differ from their germline predecessors. Thus although the mature bnAbs bind to conserved epitopes of the HIV-1 envelope glycoprotein (Env) with high affinity their germline predecessors do not or weakly bind Envs failing to initiate an effective immune response. The identification of less somatically mutated bnAbs and/or antibody maturation intermediates that are clonally related to bnAbs may be useful to circumvent the major problem of initiating immune responses leading to elicitation of bnAbs. Here, we describe the identification of IgG antibodies from an acutely HIV-1-infected patient using a combination of phage display and high-throughput sequencing. We found two antibodies with only a single point mutation in the V region of their heavy chain variable domains compared to their putative germline predecessors which bound with high affinity to several Envs. They targeted the Env gp41 and did not neutralize HIV-1. Using high-throughput sequencing, we identified several highly abundant CDR3s, germline-like as well as somatically mutated V genes in the VH/VL repertoires of the patient which may provide antibody intermediates corresponding to known bnAbs as templates for design of novel HIV-1 vaccine immunogens.

Onset of immune senescence defined by unbiased pyrosequencing of human immunoglobulin mRNA repertoires

The immune system protects us from foreign substances or pathogens by generating specific antibodies. The variety of immunoglobulin (Ig) paratopes for antigen recognition is a result of the V(D)J rearrangement mechanism, while a fast and efficient immune response is mediated by specific immunoglobulin isotypes obtained through class switch recombination (CSR). To get a better understanding on how antibody-based immune protection works and how it changes with age, the interdependency between these two parameters need to be addressed. Here, we have performed an in depth analysis of antibody repertoires of 14 healthy donors representing different gender and age groups. For this task, we developed a unique pyrosequencing approach, which is able to monitor the expression levels of all immunoglobulin V(D)J recombinations of all isotypes including subtypes in an unbiased and quantitative manner. Our results show that donors have individual immunoglobulin repertoires and cannot be clustered according to V(D)J recombination patterns, neither by age nor gender. However, after incorporating isotype-specific analysis and considering CSR information into hierarchical clustering the situation changes. For the first time the donors cluster according to age and separate into young adults and elderly donors (>50). As a direct consequence, this clustering defines the onset of immune senescence at the age of fifty and beyond. The observed age-dependent reduction of CSR ability proposes a feasible explanation why reduced efficacy of vaccination is seen in the elderly and implies that novel vaccine strategies for the elderly should include the "Golden Agers".

B cell exchange across the blood-brain barrier in multiple sclerosis

In multiple sclerosis (MS) pathogenic B cells likely act on both sides of the blood-brain barrier (BBB). However, it is unclear whether antigen-experienced B cells are shared between the CNS and the peripheral blood (PB) compartments. We applied deep repertoire sequencing of IgG heavy chain variable region genes (IgG-VH) in paired cerebrospinal fluid and PB samples from patients with MS and other neurological diseases to identify related B cells that are common to both compartments. For the first time to our knowledge, we found that a restricted pool of clonally related B cells participated in robust bidirectional exchange across the BBB. Some clusters of related IgG-VH appeared to have undergone active diversification primarily in the CNS, while others have undergone active diversification in the periphery or in both compartments in parallel. B cells are strong candidates for autoimmune effector cells in MS, and these findings suggest that CNS-directed autoimmunity may be triggered and supported on both sides of the BBB. These data also provide a powerful approach to identify and monitor B cells in the PB that correspond to clonally amplified populations in the CNS in MS and other inflammatory states.

HLA DNA typing: past, present, and future

The human leukocyte antigen (HLA) class I and class II loci are the most polymorphic genes in the human genome, with a highly clustered and patchwork pattern of sequence motifs. In the three decades since the first HLA gene was isolated by molecular cloning (a cDNA clone of HLA-B7), thousands of alleles have been identified and the names and sequences of all known alleles have been curated in the IMGT/HLA database. This extensive allelic diversity made and continues to make high-resolution HLA DNA typing very challenging. The first attempt at HLA DNA typing involved restriction fragment length polymorphism (RFLP) analysis, but this approach had many limitations. The development of PCR in 1985 allowed for the amplification of the polymorphic exons of the HLA class I and class II genes and the analysis of the polymorphic sequence motifs with sequence-specific oligonucleotide (SSO) hybridization probes. The immobilization of these probes on membranes and later on beads, along with primer sets for sequence-specific priming (SSP), gave rise to the current set of HLA typing reagents. Sanger sequencing has provided high-resolution typing but, in many cases, genotyping 'ambiguity' remains an issue. In the past few years, the introduction of next-generation sequencing, with the critical properties of massively parallel and clonal sequencing, has significantly reduced HLA genotyping ambiguity. Here, our lab's efforts to develop high-resolution and high-throughput HLA DNA typing using the 454 Sequencing System are reviewed, and the potential future developments and applications of HLA DNA typing are discussed.

Comprehensive assessment of potential multiple myeloma immunoglobulin heavy chain V-D-J intraclonal variation using massively parallel pyrosequencing

Multiple myeloma (MM) is characterized by the accumulation of malignant plasma cells (PCs) in the bone marrow (BM). MM is viewed as a clonal disorder due to lack of verified intraclonal sequence diversity in the immunoglobulin heavy chain variable region gene (IGHV). However, this conclusion is based on analysis of a very limited number of IGHV subclones and the methodology employed did not permit simultaneous analysis of the IGHV repertoire of non-malignant PCs in the same samples. Here we generated genomic DNA and cDNA libraries from purified MM BMPCs and performed massively parallel pyrosequencing to determine the frequency of cells expressing identical IGHV sequences. This method provided an unprecedented opportunity to interrogate the presence of clonally related MM cells and evaluate the IGHV repertoire of non-MM PCs. Within the MM sample, 37 IGHV genes were expressed, with 98.9% of all immunoglobulin sequences using the same IGHV gene as the MM clone and 83.0% exhibiting exact nucleotide sequence identity in the IGHV and heavy chain complementarity determining region 3 (HCDR3). Of interest, we observed in both genomic DNA and cDNA libraries 48 sets of identical sequences with single point mutations in the MM clonal IGHV or HCDR3 regions. These nucleotide changes were suggestive of putative subclones and therefore were subjected to detailed analysis to interpret: 1) their legitimacy as true subclones; and 2) their significance in the context of MM. Finally, we report for the first time the IGHV repertoire of normal human BMPCs and our data demonstrate the extent of IGHV repertoire diversity as well as the frequency of clonally-related normal BMPCs. This study demonstrates the power and potential weaknesses of in-depth sequencing as a tool to thoroughly investigate the phylogeny of malignant PCs in MM and the IGHV repertoire of normal BMPCs.

Human antibodies that neutralize HIV-1: identification, structures, and B cell ontogenies

Antibodies that neutralize diverse strains of HIV-1 develop in ∼20% of HIV-1-infected individuals, and isolation and structural characterization of these antibodies are revealing how they recognize the envelope glycoprotein spike. Broadly reactive neutralizing antibodies utilize just a few sites of spike vulnerability and converge on select modes of recognition. These antibodies have unusual features: uncommonly long complementarity-determining loops, extensive somatic mutation, or both. Recent advances in next-generation sequencing of antibody-gene transcripts are providing genetic records of the development of neutralizing antibodies. These records inform an understanding of the naive B cell repertoire, of somatic mutation, and of the resulting antibody features that are critical to effective HIV-1 neutralization; based on these, we propose an ontogeny and structure-based system of antibody classification. The human immune system is capable of developing antibodies that broadly neutralize HIV-1--and an increasingly detailed view is accumulating for how effective immunity against HIV-1 can be generated.

Somatic Populations of PGT135-137 HIV-1-Neutralizing Antibodies Identified by 454 Pyrosequencing and Bioinformatics

Select HIV-1-infected individuals develop sera capable of neutralizing diverse viral strains. The molecular basis of this neutralization is currently being deciphered by the isolation of HIV-1-neutralizing antibodies. In one infected donor, three neutralizing antibodies, PGT135–137, were identified by assessment of neutralization from individually sorted B cells and found to recognize an epitope containing an N-linked glycan at residue 332 on HIV-1 gp120. Here we use next-generation sequencing and bioinformatics methods to interrogate the B cell record of this donor to gain a more complete understanding of the humoral immune response. PGT135–137-gene family specific primers were used to amplify heavy-chain and light-chain variable-domain sequences. Pyrosequencing produced 141,298 heavy-chain sequences of IGHV4-39 origin and 87,229 light-chain sequences of IGKV3-15 origin. A number of heavy and light-chain sequences of ∼90% identity to PGT137, several to PGT136, and none of high identity to PGT135 were identified. After expansion of these sequences to include close phylogenetic relatives, a total of 202 heavy-chain sequences and 72 light-chain sequences were identified. These sequences were clustered into populations of 95% identity comprising 15 for heavy chain and 10 for light chain, and a select sequence from each population was synthesized and reconstituted with a PGT137-partner chain. Reconstituted antibodies showed varied neutralization phenotypes for HIV-1 clade A and D isolates. Sequence diversity of the antibody population represented by these tested sequences was notably higher than observed with a 454 pyrosequencing-control analysis on 10 antibodies of defined sequence, suggesting that this diversity results primarily from somatic maturation. Our results thus provide an example of how pathogens like HIV-1 are opposed by a varied humoral immune response, derived from intrinsic mechanisms of antibody development, and embodied by somatic populations of diverse antibodies.

Candidate antibody-based therapeutics against HIV-1

Antibody-based therapeutics have been successfully used for the treatment of various diseases and as research tools. Several well characterized, broadly neutralizing monoclonal antibodies (bnmAbs) targeting HIV-1 envelope glycoproteins or related host cell surface proteins show sterilizing protection of animals, but they are not effective when used for therapy of an established infection in humans. Recently, a number of novel bnmAbs, engineered antibody domains (eAds), and multifunctional fusion proteins have been reported which exhibit exceptionally potent and broad neutralizing activity against a wide range of HIV-1 isolates from diverse genetic subtypes. eAds could be more effective in vivo than conventional full-size antibodies generated by the human immune system. Because of their small size (12∼15 kD), they can better access sterically restricted epitopes and penetrate densely packed tissue where HIV-1 replicates than the larger full-size antibodies. HIV-1 possesses a number of mechanisms to escape neutralization by full-size antibodies but could be less likely to develop resistance to eAds. Here, we review the in vitro and in vivo antiviral efficacies of existing HIV-1 bnmAbs, summarize the development of eAds and multispecific fusion proteins as novel types of HIV-1 inhibitors, and discuss possible strategies to generate more potent antibody-based candidate therapeutics against HIV-1, including some that could be used to eradicate the virus.

Next generation deep sequencing and vaccine design: today and tomorrow

Next generation sequencing (NGS) technologies have redefined the modus operandi in both human and microbial genetics research, allowing the unprecedented generation of very large sequencing datasets on a short time scale and at affordable costs. Vaccine development research is rapidly taking full advantage of the advent of NGS. This review provides a concise summary of the current applications of NGS in relation to research seeking to develop vaccines for human infectious diseases, incorporating studies of both the pathogen and the host. We focus on rapidly mutating viral pathogens, which are major targets in current vaccine research. NGS is unraveling the complex dynamics of viral evolution and host responses against these viruses, thus contributing substantially to the likelihood of successful vaccine development.

Massive evolution of the immunoglobulin heavy chain locus in children with B precursor acute lymphoblastic leukemia

The ability to distinguish clonal B-cell populations based on the sequence of their rearranged immunoglobulin heavy chain (IgH) locus is an important tool for diagnosing B-cell neoplasms and monitoring treatment response. Leukemic precursor B cells may continue to undergo recombination of the IgH gene after malignant transformation; however, the magnitude of evolution at the IgH locus is currently unknown. We used next-generation sequencing to characterize the repertoire of IgH sequences in diagnostic samples of 51 children with B precursor acute lymphoblastic leukemia (B-ALL). We identified clonal IgH rearrangements in 43 of 51 (84%) cases and found that the number of evolved IgH sequences per patient ranged dramatically from 0 to 4024. We demonstrate that the evolved IgH sequences are not the result of amplification artifacts and are unique to leukemic precursor B cells. In addition, the evolution often follows an allelic exclusion pattern, where only 1 of 2 rearranged IgH loci exhibit ongoing recombination. Thus, precursor B-cell leukemias maintain evolution at the IgH locus at levels that were previously underappreciated. This finding sheds light on the mechanisms associated with leukemic clonal evolution and may fundamentally change approaches for monitoring minimal residual disease burden.

PCR amplification and high throughput sequencing of immunoglobulin heavy chain genes from formalin-fixed paraffin-embedded human biopsies

The use of high throughput sequencing (HTS) technologies in biomedicine is expanding in a variety of fields in recent years. The 454 system is an HTS platform that is ideally suited to characterize B cell receptor (BCR) repertoires by sequencing of immunoglobulin (Ig) genes, as it is able to sequence stretches of several hundred nucleotides. Most studies that used this platform for antibody repertoire analyses have started from fresh or frozen tissues or peripheral blood samples, and rely on starting with optimal quality DNA. In this paper we demonstrate that BCR repertoire analysis can be done using DNA from formalin-fixed paraffin-embedded (FFPE) human tissue samples. The heterogeneity of BCR repertoires we obtained confirms the plausibility of HTS of DNA from FFPE specimens. The establishment of experimental protocols and computational tools that enable sequence data analysis from the low quality DNA of FFPE tissues is important for enabling research, as it would enable the use of the rich source of preserved samples in clinical biobanks and biopsy archives.

Models and methods for analysis of lymphocyte repertoire generation, development, selection and evolution

T and B cell receptor repertoires are diversified by variable region gene rearrangement and selected based on functionality and lack of self-reactivity. Repertoires can also be defined based on phenotype and function rather than receptor specificity - such as the diversity of T helper cell subsets. Natural killer (NK) cell repertoires, in which each cell expresses a randomly chosen subset of its inhibitory receptor genes, and is educated based on self-MHC recognition by yet unknown mechanisms, are also phenotypic repertoires. Studying the generation, development and selection of lymphocyte repertoires, and their functions during immune responses, is essential for understanding the function of the immune system in healthy individuals and in immune deficient, autoimmune or cancer patients. The study of lymphocyte repertoires will enable clinical immunologists to develop better therapeutic monoclonal antibodies, vaccines, transplantation donor-recipient matching protocols, and other immune intervention strategies. The recent development of high-throughput methods for repertoire data collection - from multicolor flow cytometry through single-cell imaging to deep sequencing - presents us now, for the first time, with the ability to analyze and compare large samples of lymphocyte repertoires in health, aging and disease. The exponential growth of these datasets, however, challenges the theoretical immunology community to develop methods for data organization and analysis. Furthermore, the need to test hypotheses regarding immune function, and generate predictions regarding the outcomes of medical interventions, necessitates the development of complex mathematical and computational models, covering processes on multiple scales, from the genetic and molecular to the cellular and system scales.

In vivo VL-targeted microbial superantigen induced global shifts in the B cell repertoire

To subvert host defenses, some microbial pathogens produce proteins that interact with conserved motifs in V regions of B cell Ag receptor shared by large sets of lymphocytes, which define the properties of a superantigen. Because the clonal composition of the lymphocyte pool is a major determinant of immune responsiveness, this study was undertaken to examine the in vivo effect on the host immune system of exposure to a B cell superantigen, protein L (PpL), a product of the common commensal bacterial species, Finegoldia magna, which is one of the most common pathogenic species among Gram-positive anaerobic cocci. Libraries of Vκ L chain transcripts were generated from the spleens of control and PpL-exposed mice, and the expressed Vκ rearrangements were characterized by high-throughput sequencing. A total of 120,855 sequencing reads could be assigned to a germline Vκ gene, with all 20 known Vκ subgroups represented. In control mice, we found a recurrent and consistent hierarchy of Vκ gene usage, as well as patterns of preferential Vκ-Jκ pairing. PpL exposure induced significant targeted global shifts in repertoire with reduction of Vκ that contain the superantigen binding motif in all exposed mice. We found significant targeted reductions in the expression of clonotypes encoded by 14 specific Vκ genes with the predicted PpL binding motif. These rigorous surveys document the capacity of a microbial protein to modulate the composition of the expressed lymphocyte repertoire, which also has broad potential implications for host-microbiome and host-pathogen relationships.

New tools for classification and monitoring of autoimmune diseases

Rheumatologists see patients with a range of autoimmune diseases. Phenotyping these diseases for diagnosis, prognosis and selection of therapies is an ever increasing problem. Advances in multiplexed assay technology at the gene, protein, and cellular level have enabled the identification of 'actionable biomarkers'; that is, biological metrics that can inform clinical practice. Not only will such biomarkers yield insight into the development, remission, and exacerbation of a disease, they will undoubtedly improve diagnostic sensitivity and accuracy of classification, and ultimately guide treatment. This Review provides an introduction to these powerful technologies that could promote the identification of actionable biomarkers, including mass cytometry, protein arrays, and immunoglobulin and T-cell receptor high-throughput sequencing. In our opinion, these technologies should become part of routine clinical practice for the management of autoimmune diseases. The use of analytical tools to deconvolve the data obtained from use of these technologies is also presented here. These analyses are revealing a more comprehensive and interconnected view of the immune system than ever before and should have an important role in directing future treatment approaches for autoimmune diseases.

High-throughput antibody sequencing reveals genetic evidence of global regulation of the naïve and memory repertoires that extends across individuals

Vast diversity in the antibody repertoire is a key component of the adaptive immune response. This diversity is generated centrally through the assembly of variable, diversity and joining gene segments, and peripherally by somatic hypermutation and class-switch recombination. The peripheral diversification process is thought to only occur in response to antigenic stimulus, producing antigen-selected memory B cells. Surprisingly, analyses of the variable, diversity and joining gene segments have revealed that the naïve and memory subsets are composed of similar proportions of these elements. Lacking, however, is a more detailed study, analyzing the repertoires of naïve and memory subsets at the level of the complete V(D)J recombinant. This report presents a thorough examination of V(D)J recombinants in the human peripheral blood repertoire, revealing surprisingly large repertoire differences between circulating B-cell subsets and providing genetic evidence for global control of repertoire diversity in naïve and memory circulating B-cell subsets.

The immunoglobulin heavy chain locus: genetic variation, missing data, and implications for human disease

The immunoglobulin (IG) loci consist of repeated and highly homologous sets of genes of different types, variable (V), diversity (D) and junction (J), that rearrange in developing B cells to produce an individual's highly variable repertoire of expressed antibodies, designed to bind to a vast array of pathogens. This repeated structure makes these loci susceptible to a high frequency of insertion and deletion events through evolutionary time, and also makes them difficult to characterize at the genomic level or assay with high-throughput techniques. Given the central role of antibodies in the adaptive immune system, it is not surprising that early candidate gene approaches showed that germline polymorphisms in these regions correlated with susceptibility to both infectious and autoimmune diseases. However, more recent studies, particularly those using high-throughput genome-wide arrays, have failed to implicate these loci in disease. In this review of the IG heavy chain variable gene cluster (IGHV), we examine how poorly we understand the distribution of haplotype variation in this genomic region, and we argue that this lack of information may mask candidate loci in the IGHV gene cluster as causative factors for infectious and autoimmune diseases.

Expressed antibody repertoires in human cord blood cells: 454 sequencing and IMGT/HighV-QUEST analysis of germline gene usage, junctional diversity, and somatic mutations

Human cord blood cell-derived IgM antibodies are important for the neonate immune responses and construction of germline-based immunoglobulin libraries. Several previous studies of a relatively small number of sequences found that they exhibit restrictions in the usage of germline genes and in the diversity of the variable heavy chain complementarity determining region 3 compared to adults. To further characterize such restrictions on a larger scale and to compare the early B-cell diversity to adult IgM repertoires, we performed 454 sequencing and IMGT/HighV-QUEST analysis of cord blood IG libraries from two babies and determined germline gene usage, V-D-J rearrangement, VHCDR3 diversity, and somatic mutations to characterize human neonate repertoire. Most of the germline subgroups were identified with frequencies comparable to those present in the adult IgM repertoire except for the IGHV1-2 gene that was preferentially expressed in the cord blood cells. The gene usage diversity contributed to 1,430 unique IGH V-D-J rearrangement patterns while the exonuclease trimming and N region addition at the V-D-J junctions along with gene diversity created a wide range of VHCDR3 with different lengths and sequence variability. We observed a lower degree of somatic mutations in the CDR and framework regions of antibodies from cord blood cells compared to adults. These results provide insights into the characteristics of human cord blood antibody repertoires, which have gene usage diversity and VHCDR3 lengths similar to that of the adult IgM repertoire but differ significantly in some of the gene usages, V-D-J rearrangements, junctional diversity, and somatic mutations.

Systems biology approach for new target and biomarker identification

The pharmaceutical industry is spending increasingly large amounts of money on the discovery and development of novel medicines, but this investment is not adequately paying off in an increased rate of newly approved drugs by the FDA. The post-genomic era has provided a wealth of novel approaches for generating large, high-dimensional genetic and transcriptomic data sets from large cohorts of preclinical species as well as normal and diseased individuals. This systems biology approach to understanding disease-related biology is revolutionizing our understanding of the cellular pathways and gene networks underlying the onset of disease, and the mechanisms of pharmacological treatments that ameliorate disease phenotypes. In this article, we review a number of approaches being used by pharmaceutical and biotechnology companies, e.g., high-throughput DNA genotyping, sequencing, and genome-wide gene expression profiling, to enable drug discovery and development through the identification of new drug targets and biomarkers of disease progression, drug pharmacodynamics, and predictive markers for selecting the patients most likely to respond to therapy.

New frontiers of primary antibody deficiencies

Primary antibody deficiencies (PAD) form the largest group of inherited disorders of the immune system. They are characterized by a marked reduction or absence of serum immunoglobulins (Ig) due to disturbed B cell differentiation and by a poor response to vaccination. PAD can be divided into agammaglobulinemia, Ig class switch recombination deficiencies, and idiopathic hypogammaglobulinemia. Over the past 20 years, defects have been identified in 18 different genes, but in many PAD patients the underlying gene defects have not been found. Diagnosis of known PAD and discovery of new PAD is important for good patient care. In this review, we present the effects of genetic defects in the context of normal B cell differentiation, and we discuss how new technical developments can support understanding and discovering new genetic defects in PAD.

The inference of phased haplotypes for the immunoglobulin H chain V region gene loci by analysis of VDJ gene rearrangements

The existence of many highly similar genes in the lymphocyte receptor gene loci makes them difficult to investigate, and the determination of phased "haplotypes" has been particularly problematic. However, V(D)J gene rearrangements provide an opportunity to infer the association of Ig genes along the chromosomes. The chromosomal distribution of H chain genes in an Ig genotype can be inferred through analysis of VDJ rearrangements in individuals who are heterozygous at points within the IGH locus. We analyzed VDJ rearrangements from 44 individuals for whom sufficient unique rearrangements were available to allow comprehensive genotyping. Nine individuals were identified who were heterozygous at the IGHJ6 locus and for whom sufficient suitable VDJ rearrangements were available to allow comprehensive haplotyping. Each of the 18 resulting IGHV│IGHD│IGHJ haplotypes was unique. Apparent deletion polymorphisms were seen that involved as many as four contiguous, functional IGHV genes. Two deletion polymorphisms involving multiple contiguous IGHD genes were also inferred. Three previously unidentified gene duplications were detected, where two sequences recognized as allelic variants of a single gene were both inferred to be on a single chromosome. Phased genomic data brings clarity to the study of the contribution of each gene to the available repertoire of rearranged VDJ genes. Analysis of rearrangement frequencies suggests that particular genes may have substantially different yet predictable propensities for rearrangement within different haplotypes. Together with data highlighting the extent of haplotypic variation within the population, this suggests that there may be substantial variability in the available Ab repertoires of different individuals.

Naive antibody gene-segment frequencies are heritable and unaltered by chronic lymphocyte ablation

A diverse antibody repertoire is essential for an effective adaptive immune response to novel molecular surfaces. Although past studies have observed common patterns of V-segment use, as well as variation in V-segment use between individuals, the relative contributions to variance from genetics, disease, age, and environment have remained unclear. Using high-throughput sequence analysis of monozygotic twins, we show that variation in naive V(H) and D(H) segment use is strongly determined by an individual's germ-line genetic background. The inherited segment-use profiles are resilient to differential environmental exposure, disease processes, and chronic lymphocyte depletion therapy. Signatures of the inherited profiles were observed in class switched germ-line use of each individual. However, despite heritable segment use, the rearranged complementarity-determining region-H3 repertoires remained highly specific to the individual. As it has been previously demonstrated that certain V-segments exhibit biased representation in autoimmunity, lymphoma, and viral infection, we anticipate our findings may provide a unique mechanism for stratifying individual risk profiles in specific diseases.

Antibody repertoire development in fetal and neonatal piglets XXI. Usage of most VH genes remains constant during fetal and postnatal development

Usage of variable region gene segments during development of the antibody repertoire in mammals is unresolved in part because of the complexity of the locus in mice and humans and the difficulty of distinguishing intrinsic from extrinsic influences in these species. We present the first vertical studies on VH usage that spans the fetal and neonatal period using the piglet model. We tracked VH usage in DNA rearrangements and in VDJ transcripts throughout 75 days of gestation (DG) in outbred fetuses, thereafter in outbred germfree and colonized isolator piglets, isolator piglets infected with swine influenza and in conventionally reared nematode-infected adults. Seven VH genes account for >90% of the pre-immune repertoire which is the same among tissues and in both transcripts and DNA rearrangements. Statistical modeling supports the view that proportional usage of the major genes remains constant during fetal life and that postnatal usage ranking is similar to that during fetal life. Changes in usage ranking are developmental not antigen dependent. In this species exposure to environmental antigens results in diversification of the repertoire by somatic hypermutation of the same small number of VH genes that comprise the pre-immune repertoire, not by using other VH gene available in the germline. Therefore in swine a small number of VH genes shape the antibody repertoire throughout life questioning the need for extensive VH polygeny.

454 antibody sequencing - error characterization and correction

BACKGROUND

454 sequencing is currently the method of choice for sequencing of antibody repertoires and libraries containing large numbers (106 to 1012) of different molecules with similar frameworks and variable regions which poses significant challenges for identifying sequencing errors. Identification and correction of sequencing errors in such mixtures is especially important for the exploration of complex maturation pathways and identification of putative germline predecessors of highly somatically mutated antibodies. To quantify and correct errors incorporated in 454 antibody sequencing, we sequenced six antibodies at different known concentrations twice over and compared them with the corresponding known sequences as determined by standard Sanger sequencing.

RESULTS

We found that 454 antibody sequencing could lead to approximately 20% incorrect reads due to insertions that were mostly found at shorter homopolymer regions of 2-3 nucleotide length, and less so by insertions, deletions and other variants at random sites. Correction of errors might reduce this population of erroneous reads down to 5-10%. However, there are a certain number of errors accounting for 4-8% of the total reads that could not be corrected unless several repeated sequencing is performed, although this may not be possible for large diverse libraries and repertoires including complete sets of antibodies (antibodyomes).

CONCLUSIONS

The experimental test procedure carried out for assessing 454 antibody sequencing errors reveals high (up to 20%) incorrect reads; the errors can be reduced down to 5-10% but not less which suggests the use of caution to avoid false discovery of antibody variants and diversity.

Immunoglobulin A: a bridge between innate and adaptive immunity

PURPOSE OF REVIEW

The review summarizes the recent progress that has been made in understanding the function of immunoglobulin A (IgA) in promoting a healthy mutualism with the commensal microbiota and protecting against pathogens. Although IgA is by far the most abundant antibody produced by mammals, direct experimental evidence for its function is still lacking.

RECENT FINDINGS

IgA is the predominant antibody induced in response to intestinal colonization with commensal bacteria: even fish have been shown to have a mucosal immunoglobulin (IgT), which is produced in the mucosa and coats commensals in the intestinal lumen. Recent studies indicate that intestinal IgA can be highly specific to the inducing commensals. Priming of IgA also appears to be a long-lasting response dependent on the overall dose (integral) of the bacteria sampled rather than exhibiting prime-boost effects normally observed with systemic immunoglobulin responses. Not only is human IgA highly mutated, but a mouse model with deficient hypermutation but intact class-switch recombination also shows that this mutation process (presumably leading to better anticommensal affinities) is important for IgA protection at the mucosal surface. It has been shown that some IgA can be induced independently of T cells through stimulation by epithelial cell and plasmacytoid dendritic cell cytokines including BAFF and APRIL, although the relative roles of the T-dependent and T-independent IgA pathways in generating mucosal protection are still unclear.

SUMMARY

Protection at mucosal surfaces through the secretion of antibodies is a phylogenetically ancient function. Mammals can produce high and low-affinity IgA against their commensal microbes via T-cell-dependent and T-cell-independent pathways to contribute to host microbial mutualism. The process of improving IgA affinity to intestinal luminal contents through somatic hypermutation of immunoglobulin genes improves the level of protection at the mucosal surface and such mutations are abundant in human IgA sequences.

Ultra-sensitive detection of rare T cell clones

Advances in high-throughput sequencing have enabled technologies that probe the adaptive immune system with unprecedented depth. We have developed a multiplex PCR method to sequence tens of millions of T cell receptors (TCRs) from a single sample in a few days. A method is presented to test the precision, accuracy, and sensitivity of this assay. T cell clones, each with one fixed productive TCR rearrangement, are doped into complex blood cell samples. TCRs from a total of eleven samples are sequenced, with the doped T cell clones ranging from 10% of the total sample to 0.001% (one cell in 100,000). The assay is able to detect even the rarest clones. The precision of the assay is demonstrated across five orders of magnitude. The accuracy for each clone is within an overall factor of three across the 100,000 fold dynamic range. Additionally, the assay is shown to be highly repeatable.

Epitope-specific human influenza antibody repertoires diversify by B cell intraclonal sequence divergence and interclonal convergence

We generated from a single blood sample five independent human mAbs that recognized the Sa antigenic site on the head of influenza hemagglutinin and exhibited inhibitory activity against a broad panel of H1N1 strains. All five Abs used the V(H)3-7 and J(H)6 gene segments, but at least four independent clones were identified by junctional analysis. High-throughput sequence analysis of circulating B cells revealed that each of the independent clones were members of complex phylogenetic lineages that had diversified widely using a pattern of progressive diversification through somatic mutation. Unexpectedly, B cells encoding multiple diverging lineages of these clones, including many containing very few mutations in the Ab genes, persisted in the circulation. Conversely, we noted frequent instances of amino acid sequence convergence in the Ag combining sites exhibited by members of independent clones, suggesting a strong selection for optimal binding sites. We suggest that maintenance in circulation of a wide diversity of somatic variants of dominant clones may facilitate recognition of drift variant virus epitopes that occur in rapidly mutating virus Ags, such as influenza hemagglutinin. In fact, these Ab clones recognize an epitope that acquired three glycosylation sites mediating escape from previously isolated human Abs.

Divergent human populations show extensive shared IGK rearrangements in peripheral blood B cells

We have analysed the transcribed immunoglobulin kappa (IGK) repertoire of peripheral blood B cells from four individuals from two genetically distinct populations, Papua New Guinean and Australian, using high-throughput DNA sequencing. The depth of sequencing data for each individual averaged 5,548 high-quality IGK reads, and permitted genotyping of the inferred IGKV and IGKJ germline gene segments for each individual. All individuals were homozygous at each IGKJ locus and had highly similar inferred IGKV genotypes. Preferential gene usage was seen at both the IGKV and IGKJ loci, but only IGKV segment usage varied significantly between individuals. Despite the differences in IGKV gene utilisation, the rearranged IGK repertoires showed extensive identity at the amino acid level. Public rearrangements (those shared by two or more individuals) made up 60.2% of the total sequenced IGK rearrangements. The total diversity of IGK rearrangements of each individual was estimated to range from just 340 to 549 unique amino acid sequences. Thus, the repertoire of unique expressed IGK rearrangements is dramatically less than previous theoretical estimates of IGK diversity, and the majority of expressed IGK rearrangements are likely to be extensively shared in individual human beings.

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

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

Diversity of immunoglobulin E-encoding transcripts in sinus mucosa of subjects diagnosed with non-allergic fungal eosinophilic sinusitis

BACKGROUND

The role of allergy in the aetiopathogenesis of chronic rhinosinusitis (CRS) remains controversial. For example, in some cases with sinus fungal infections allergy can be demonstrated by standard tests. In other cases, such signs can be absent despite elevated levels of IgE-positive cells in sinus tissue and the presence of IgE and eosinophils in the sinus mucous.

OBJECTIVE

To define the nature of molecular diversity in antibodies of the IgE isotype at the site of local inflammation in subjects diagnosed with non-allergic fungal eosinophilic sinusitis (NAFES).

METHODS

The local occurrence and sequence characteristics of IgE-encoding transcripts in NAFES patients were investigated and compared with sequences found in subjects diagnosed with CRS featuring systemic allergy. These sequences have also been compared with other reported IgE-encoding transcriptomes. Results IGHV genes derived from major subgroups 1, 3, 4 and 5 and a diverse set of IGHD and IGHJ genes were shown to create the IgE repertoire in patients diagnosed with NAFES and CRS. The average lengths of the third hypervariable loop in these populations were 15.8 and 14.6 residues. The sequences showed evidence of extensive somatic hypermutation (mutation frequency: NAFES, 6.4 ± 3.2%; CRS, 7.0 ± 4.4%) with substitutions targeted to complementarity-determining regions. These sequence collections thus show extensive similarities to those found in other polyclonal Ig repertoires including those encoding IgE.

CONCLUSION AND CLINICAL RELEVANCE

We conclude that sinus IgE-encoding transcripts in subjects diagnosed with NAFES show evidence of conventional IgE responses and we suggest that allergens with characteristics of classical antigens should be investigated for a role in the local response occurring in NAFES. This investigation illustrates that assessment of local immunity might be an important diagnostic tool in conditions like NAFES with no systemic signs of allergy to identify or rule out an allergic component of the patient's disease.

Genomic screening by 454 pyrosequencing identifies a new human IGHV gene and sixteen other new IGHV allelic variants

Complete and accurate knowledge of the genes and allelic variants of the human immunoglobulin gene loci is critical for studies of B cell repertoire development and somatic point mutation, but evidence from studies of VDJ rearrangements suggests that our knowledge of the available immunoglobulin gene repertoire is far from complete. The reported repertoire has changed little over the last 15 years. This is, in part, a consequence of the inefficiencies involved in searching for new members of large, multigenic gene families by cloning and sequencing. The advent of high-throughput sequencing provides a new avenue by which the germline repertoire can be explored. In this report, we describe pyrosequencing studies of the heavy chain IGHV1, IGHV3 and IGHV4 gene subgroups in ten Papua New Guineans. Thousands of 454 reads aligned with complete identity to 51 previously reported functional IGHV genes and allelic variants. A new gene, IGHV3-NL1*01, was identified, which differs from the nearest previously reported gene by 15 nucleotides. Sixteen new IGHV alleles were also identified, 15 of which varied from previously reported functional IGHV genes by between one and four nucleotides, while one sequence appears to be a functional variant of the pseudogene IGHV3-25. BLAST searches suggest that at least six of these new genes are carried within the relatively well-studied populations of North America, Europe or Asia. This study substantially expands the known immunoglobulin gene repertoire and demonstrates that genetic variation of immunoglobulin genes can now be efficiently explored in different human populations using high-throughput pyrosequencing.

Benchmarking the performance of human antibody gene alignment utilities using a 454 sequence dataset

MOTIVATION

Immunoglobulin heavy chain genes are formed by recombination of genes randomly selected from sets of IGHV, IGHD and IGHJ genes. Utilities have been developed to identify genes that contribute to observed VDJ rearrangements, but in the absence of datasets of known rearrangements, the evaluation of these utilities is problematic. We have analyzed thousands of VDJ rearrangements from an individual (S22) whose IGHV, IGHD and IGHJ genotype can be inferred from the dataset. Knowledge of this genotype means that the Stanford_S22 dataset can serve to benchmark the performance of IGH alignment utilities.

RESULTS

We evaluated the performance of seven utilities. Failure to partition a sequence into genes present in the S22 genome was considered an error, and error rates for different utilities ranged from 7.1% to 13.7%.

AVAILABILITY

Supplementary data includes the S22 genotypes and alignments. The Stanford_S22 dataset and an evaluation tool is available at http://www.emi.unsw.edu.au/~ihmmune/IGHUtilityEval/.

Clustering-based identification of clonally-related immunoglobulin gene sequence sets

Background

Clonal expansion of B lymphocytes coupled with somatic mutation and antigen selection allow the mammalian humoral immune system to generate highly specific immunoglobulins (IG) or antibodies against invading bacteria, viruses and toxins. The availability of high-throughput DNA sequencing methods is providing new avenues for studying this clonal expansion and identifying the factors guiding the generation of antibodies. The identification of groups of rearranged immunoglobulin gene sequences descended from the same rearrangement (clonally-related sets) in very large sets of sequences is facilitated by the availability of immunoglobulin gene sequence alignment and partitioning software that can accurately predict component germline gene, but has required painstaking visual inspection and analysis of sequences.

Results

We have developed and implemented an algorithm for identifying sets of clonally-related sequences in large human immunoglobulin heavy chain gene variable region sequence sets. The program processes sequences that have been partitioned using iHMMune-align, and uses pairwise comparisons of CDR3 sequences and similarity in IGHV and IGHJ germline gene assignments to construct a distance matrix. Agglomerative hierarchical clustering is then used to identify likely groups of clonally-related sequences. The program is available for download from http://www.cse.unsw.edu.au/~ihmmune/ClonalRelate/ClonalRelate.zip.

Conclusions

The method was evaluated on several benchmark datasets and provided a more accurate and considerably faster identification of clonally-related immunoglobulin gene sequences than visual inspection by domain experts.

Overlap and effective size of the human CD8+ T cell receptor repertoire

Diversity in T lymphocyte antigen receptors is generated by somatic rearrangement of T cell receptor (TCR) genes and is concentrated within the third complementarity-determining region 3 (CDR3) of each chain of the TCR heterodimer. We sequenced the CDR3 regions from millions of rearranged TCR beta chain genes in naïve and memory CD8(+) T cells of seven adults. The CDR3 sequence repertoire realized in each individual is strongly biased toward specific V(beta)-J(beta) pair utilization, dominated by sequences containing few inserted nucleotides, and drawn from a defined subset comprising less than 0.1% of the estimated 5 x 10(11) possible sequences. Surprisingly, the overlap in the naïve CD8(+) CDR3 sequence repertoires of any two of the individuals is approximately 7000-fold larger than predicted and appears to be independent of the degree of human leukocyte antigen matching.