High-throughput sequencing of the zebrafish antibody repertoire

Challenges and promise for the development of human immune monitoring

The immune system is critical for protection and health maintenance and is likely required for a long lifespan. Yet, despite its importance for health, the ability to assess its quality of function has been poor, nor is much known on its variation between individuals. Hence direct assessment of immune health has largely been missing from medicine, and metrics of immune health are not well defined, especially in non-extreme states. This is chiefly due to the high complexity of the immune system. Recently emerging technologies now enable broad surveying of many immune system components at high resolution, setting forth a transformation of immunology and, through it, medicine. Such technologies enable, for the first time, high-resolution monitoring of an individual's immune system. The resulting information can be used for diagnostic and prognostic purposes, as well as to provide a quantitative, global view of the immune system, i.e. "systems immunology." This is especially relevant in the context of aging, in which the immune system exhibits profound alterations in state and function.

Chromatin conformation governs T-cell receptor Jβ gene segment usage

T cells play fundamental roles in adaptive immunity, relying on a diverse repertoire of T-cell receptor (TCR) α and β chains. Diversity of the TCR β chain is generated in part by a random yet intrinsically biased combinatorial rearrangement of variable (Vβ), diversity (Dβ), and joining (Jβ) gene segments. The mechanisms that determine biases in gene segment use remain unclear. Here we show, using a high-throughput TCR sequencing approach, that a physical model of chromatin conformation at the DJβ genomic locus explains more than 80% of the biases in Jβ use that we measured in murine T cells. This model also predicts correctly how differences in intersegment genomic distances between humans and mice translate into differences in Jβ bias between TCR repertoires of these two species. As a consequence of these structural and other biases, TCR sequences are produced with different a priori frequencies, thus affecting their probability of becoming public TCRs that are shared among individuals. Surprisingly, we find that many more TCR sequences are shared among all five mice we studied than among only subgroups of three or four mice. We derive a necessary mathematical condition explaining this finding, which indicates that the TCR repertoire contains a core set of receptor sequences that are highly abundant among individuals, if their a priori probability of being produced by the recombination process is higher than a defined threshold. Our results provide evidence for an expanded role of chromatin conformation in VDJ rearrangement, from control of gene accessibility to precise determination of gene segment use.

Statistical inference of the generation probability of T-cell receptors from sequence repertoires

Stochastic rearrangement of germline V-, D-, and J-genes to create variable coding sequence for certain cell surface receptors is at the origin of immune system diversity. This process, known as "VDJ recombination", is implemented via a series of stochastic molecular events involving gene choices and random nucleotide insertions between, and deletions from, genes. We use large sequence repertoires of the variable CDR3 region of human CD4+ T-cell receptor beta chains to infer the statistical properties of these basic biochemical events. Because any given CDR3 sequence can be produced in multiple ways, the probability distribution of hidden recombination events cannot be inferred directly from the observed sequences; we therefore develop a maximum likelihood inference method to achieve this end. To separate the properties of the molecular rearrangement mechanism from the effects of selection, we focus on nonproductive CDR3 sequences in T-cell DNA. We infer the joint distribution of the various generative events that occur when a new T-cell receptor gene is created. We find a rich picture of correlation (and absence thereof), providing insight into the molecular mechanisms involved. The generative event statistics are consistent between individuals, suggesting a universal biochemical process. Our probabilistic model predicts the generation probability of any specific CDR3 sequence by the primitive recombination process, allowing us to quantify the potential diversity of the T-cell repertoire and to understand why some sequences are shared between individuals. We argue that the use of formal statistical inference methods, of the kind presented in this paper, will be essential for quantitative understanding of the generation and evolution of diversity in the adaptive immune system.

Application of circular consensus sequencing and network analysis to characterize the bovine IgG repertoire

BACKGROUND

Vertebrate immune systems generate diverse repertoires of antibodies capable of mediating response to a variety of antigens. Next generation sequencing methods provide unique approaches to a number of immuno-based research areas including antibody discovery and engineering, disease surveillance, and host immune response to vaccines. In particular, single-molecule circular consensus sequencing permits the sequencing of antibody repertoires at previously unattainable depths of coverage and accuracy. We approached the bovine immunoglobulin G (IgG) repertoire with the objective of characterizing diversity of expressed IgG transcripts. Here we present single-molecule real-time sequencing data of expressed IgG heavy-chain repertoires of four individual cattle. We describe the diversity observed within antigen binding regions and visualize this diversity using a network-based approach.

RESULTS

We generated 49,945 high quality cDNA sequences, each spanning the entire IgG variable region from four Bos taurus calves. From these sequences we identified 49,521 antigen binding regions using the automated Paratome web server. Approximately 9% of all unique complementarity determining 2 (CDR2) sequences were of variable lengths. A bimodal distribution of unique CDR3 sequence lengths was observed, with common lengths of 5-6 and 21-25 amino acids. The average number of cysteine residues in CDR3s increased with CDR3 length and we observed that cysteine residues were centrally located in CDR3s. We identified 19 extremely long CDR3 sequences (up to 62 amino acids in length) within IgG transcripts. Network analyses revealed distinct patterns among the expressed IgG antigen binding repertoires of the examined individuals.

CONCLUSIONS

We utilized circular consensus sequencing technology to provide baseline data of the expressed bovine IgG repertoire that can be used for future studies important to livestock research. Somatic mutation resulting in base insertions and deletions in CDR2 further diversifies the bovine antibody repertoire. In contrast to previous studies, our data indicate that unusually long CDR3 sequences are not unique to IgM antibodies in cattle. Centrally located cysteine residues in bovine CDR3s provide further evidence that disulfide bond formation is likely of structural importance. We hypothesize that network or cluster-based analyses of expressed antibody repertoires from controlled challenge experiments will help identify novel natural antigen binding solutions to specific pathogens of interest.

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.

The porcine antibody repertoire: variations on the textbook theme

The genes encoding the heavy and light chains of swine antibodies are organized in the same manner as in other eutherian mammals. There are ∼30 VH genes, two functional DH genes and one functional JH gene, 14-60 Vκ genes, 5 Jκ segments, 12-13 functional Vλ genes, and two functional Jλ genes. The heavy chain constant regions encode the same repertoire of isotypes common to other eutherian mammals. The piglet models offers advantage over rodent models since the fetal repertoire develops without maternal influences and the precocial nature of their multiple offspring allows the experimenter to control the influences of environmental and maternal factors on repertoire development postnatally. B cell lymphogenesis in swine begins in the fetal yolk sac at 20 days of gestation (DG), moves to the fetal liver at 30 DG and eventually to the bone marrow which dominates until birth (114 DG) and to at least 5 weeks postpartum. There is no evidence that the ileal Peyers patches are a site of B cell lymphogenesis or are required for B cell maintenance. Unlike rodents and humans, light chain rearrangement begins first in the lambda locus; kappa rearrangements are not seen until late gestation. Dissimilar to lab rodents and more in the direction of the rabbit, swine utilize a small number of VH genes to form >90% of their pre-immune repertoire. Diversification in response to environmental antigen does not alter this pattern and is achieved by somatic hypermutation (SHM) of the same small number of VH genes. The situation for light chains is less well studied, but certain Vκ and Jκ and Vλ and Jλ are dominant in transcripts and in contrast to rearranged heavy chains, there is little junctional diversity, less SHM, and mutations are not concentrated in CDR regions. The transcribed and secreted pre-immune antibodies of the fetus include mainly IgM, IgA, and IgG3; this last isotype may provide a type of first responder mucosal immunity. Development of functional adaptive immunity is dependent on bacterial MAMPs or MAMPs provided by viral infections, indicating the importance of innate immunity for development of adaptive immunity. The structural analysis of Ig genes of this species indicate that especially the VH and Cγ gene are the result of tandem gene duplication in the context of genomic gene conversion. Since only a few of these duplicated VH genes substantially contribute to the antibody repertoire, polygeny may be a vestige from a time before somatic processes became prominently evolved to generate the antibody repertoire. In swine we believe such duplications within the genome have very limited functional significance and their occurrence is therefore overrated.

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.

Evolution of the immune system in the lower vertebrates

The evolutionary emergence of vertebrates was accompanied by the invention of adaptive immunity. This is characterized by extraordinarily diverse repertoires of somatically assembled antigen receptors and the facility of antigen-specific memory, leading to more rapid and efficient secondary immune responses. Adaptive immunity emerged twice during early vertebrate evolution, once in the lineage leading to jawless fishes (such as lamprey and hagfish) and, independently, in the lineage leading to jawed vertebrates (comprising the overwhelming majority of extant vertebrates, from cartilaginous fishes to mammals). Recent findings on the immune systems of jawless and jawed fishes (here referred to as lower vertebrates) impact on the identification of general principles governing the structure and function of adaptive immunity and its coevolution with innate defenses. The discovery of conserved features of adaptive immunity will guide attempts to generate synthetic immunological functionalities and thus provide new avenues for intervening with faulty immune functions in humans.

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.

Molecular tools for companion diagnostics

The heterogeneous nature of cancer results in highly variable therapeutic responses even among patients with identical stages and grades of a malignancy. The move towards personalised medicine in cancer therapy has therefore been motivated by a need to customise therapy according to molecular features of individual tumours. Companion diagnostics serves to support early drug development, it can provide surrogate markers in clinical trials, and also guide selection of individual therapies and monitoring of responses in routine clinical care. The era of companion diagnostics can be said to have begun with the introduction of the HercepTest - a first-of-a-kind diagnostic tool developed by DakoCytomation in 1998 to select patients for therapy with the anticancer drug Herceptin (trastuzumab). Herceptin and the paired test proved that companion diagnostics can help guide patient-tailored therapies. We will discuss herein technologies to analyse companion diagnostics markers at the level of DNA, RNA or protein, focusing on a series of methods developed in our laboratory that can facilitate drug development and help stratify patients for therapy.

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.

Antibody repertoires in humanized NOD-scid-IL2Rγ(null) mice and human B cells reveals human-like diversification and tolerance checkpoints in the mouse

Immunodeficient mice reconstituted with human hematopoietic stem cells enable the in vivo study of human hematopoiesis. In particular, NOD-scid-IL2Rγ^null engrafted mice have been shown to have reasonable levels of T and B cell repopulation and can mount T-cell dependent responses; however, antigen-specific B-cell responses in this model are generally poor. We explored whether developmental defects in the immunoglobulin gene repertoire might be partly responsible for the low level of antibody responses in this model. Roche 454 sequencing was used to obtain over 685,000 reads from cDNA encoding immunoglobulin heavy (IGH) and light (IGK and IGL) genes isolated from immature, naïve, or total splenic B cells in engrafted NOD-scid-IL2Rγ^null mice, and compared with over 940,000 reads from peripheral B cells of two healthy volunteers. We find that while naïve B-cell repertoires in humanized mice are chiefly indistinguishable from those in human blood B cells, and display highly correlated patterns of immunoglobulin gene segment use, the complementarity-determining region H3 (CDR-H3) repertoires are nevertheless extremely diverse and are specific for each individual. Despite this diversity, preferential D_H-J_H pairings repeatedly occur within the CDR-H3 interval that are strikingly similar across all repertoires examined, implying a genetic constraint imposed on repertoire generation. Moreover, CDR-H3 length, charged amino-acid content, and hydropathy are indistinguishable between humans and humanized mice, with no evidence of global autoimmune signatures. Importantly, however, a statistically greater usage of the inherently autoreactive IGHV4-34 and IGKV4-1 genes was observed in the newly formed immature B cells relative to naïve B or total splenic B cells in the humanized mice, a finding consistent with the deletion of autoreactive B cells in humans. Overall, our results provide evidence that key features of the primary repertoire are shaped by genetic factors intrinsic to human B cells and are principally unaltered by differences between mouse and human stromal microenvironments.

Rep-Seq: uncovering the immunological repertoire through next-generation sequencing

Recent scientific discoveries fuelled by the application of next-generation DNA and RNA sequencing technologies highlight the striking impact of these platforms in characterizing multiple aspects in genomics research. This technology has been used in the study of the B-cell and T-cell receptor repertoire. The novelty of immunosequencing comes from the recent rapid development of techniques and the exponential reduction in cost of sequencing. Here, we describe some of the technologies, which we collectively refer to as Rep-Seq (repertoire sequencing), to portray achievements in the field and to present the essential and inseparable role of next-generation sequencing to the understanding of entities in immune response. The large Rep-Seq data sets that should be available in the near future call for new computational algorithms to segue the transition from 'classic' molecular-based analysis to system-wide analysis. The combination of new algorithms with high-throughput data will form the basis for possible new clinical implications in personalized medicine and deeper understanding of immune behaviour and immune response.

A novel DNAseq program for enhanced analysis of Illumina GAII data: a case study on antibody complementarity-determining regions

High-throughput DNA sequencing technologies are increasingly becoming powerful systems for the comprehensive analysis of variations in whole genomes or various DNA libraries. As they are capable of producing massive collections of short sequences with varying lengths, a major challenge is how to turn these reads into biologically meaningful information. The first stage is to assemble the short reads into longer sequences through an in silico process. However, currently available software/programs allow only the assembly of abundant sequences, which apparently results in the loss of highly variable (or rare) sequences or creates artefact assemblies. In this paper, we describe a novel program (DNAseq) that is capable of assembling highly variable sequences and displaying them directly for phylogenetic analysis. In addition, this program is Microsoft Windows-based and runs by a normal PC with 700MB RAM for a general use. We have applied it to analyse a human naive single-chain antibody (scFv) library, comprehensively revealing the diversity of antibody variable complementarity-determining regions (CDRs) and their families. Although only a scFv library was exemplified here, we envisage that this program could be applicable to other genome libraries.

Structural elucidation of the mechanistic basis of degeneracy in the primary humoral response

The mechanistic basis for efficient combating of the infinite range of foreign Ags by the limited repertoire of naive Abs expressed on primary B cell surfaces during their first encounter was addressed through elegantly designed crystallographic analyses. Resolution of the discrepancy arising from the limited number of possible germline Ab receptors on primary B cells for recognizing the unlimited pool of possible Ags has been attempted by invoking the degenerate recognition potential of the germline Abs. Structural analyses of germline mAb BBE6.12H3 in an Ag-free state, as well as bound to four different peptide Ags, established the correlation of its degenerate specificity with conformational versatility of the paratope. Six distinct paratope topologies observed for a single germline mAb provided a quantitative description of the primary Ag recognition repertoire at the tertiary structural level. Each of the four different peptide Ags was bound specifically to a distinct conformation of the paratope, which was also different from that of the Ag-free states of the same germline mAb. A minimal conserved motif in the pristine Ag-combining site essential for multispecificity and Ag binding-mediated change in the elbow angle of Fab was also discernible. It is proposed that the generation of a primary Ab repertoire involves large, yet finite, germline Ab clones, each capable of adopting discrete conformations, which in turn exhibit diverse binding modes.

Age, microbiota, and T cells shape diverse individual IgA repertoires in the intestine

Intestinal immunoglobulin A (IgA) ensures host defense and symbiosis with our commensal microbiota. Yet previous studies hint at a surprisingly low diversity of intestinal IgA, and it is unknown to what extent the diverse Ig arsenal generated by somatic recombination and diversification is actually used. In this study, we analyze more than one million mouse IgA sequences to describe the shaping of the intestinal IgA repertoire, its determinants, and stability over time. We show that expanded and infrequent clones combine to form highly diverse polyclonal IgA repertoires with very little overlap between individual mice. Selective homing allows expanded clones to evenly seed the small but not large intestine. Repertoire diversity increases during aging in a dual process. On the one hand, microbiota-, T cell-, and transcription factor RORγt-dependent but Peyer's patch-independent somatic mutations drive the diversification of expanded clones, and on the other hand, new clones are introduced into the repertoire of aged mice. An individual's IgA repertoire is stable and recalled after plasma cell depletion, which is indicative of functional memory. These data provide a conceptual framework to understand the dynamic changes in the IgA repertoires to match environmental and intrinsic stimuli.

Determinants of public T cell responses

Historically, sharing T cell receptors (TCRs) between individuals has been speculated to be impossible, considering the dramatic discrepancy between the potential enormity of the TCR repertoire and the limited number of T cells generated in each individual. However, public T cell response, in which multiple individuals share identical TCRs in responding to a same antigenic epitope, has been extensively observed in a variety of immune responses across many species. Public T cell responses enable individuals within a population to generate similar antigen-specific TCRs against certain ubiquitous pathogens, leading to favorable biological outcomes. However, the relatively concentrated feature of TCR repertoire may limit T cell response in a population to some other pathogens. It could be a great benefit for human health if public T cell responses can be manipulated. Therefore, the mechanistic insight of public TCR generation is important to know. Recently, high-throughput DNA sequencing has revolutionized the study of immune receptor repertoires, which allows a much better understanding of the factors that determine the overlap of TCR repertoire among individuals. Here, we summarize the current knowledge on public T-cell response and discuss future challenges in this field.

Mass spectrometry analyses of κ and λ fractions result in increased number of complementarity-determining region identifications

Sera from lung cancer patients contain antibodies against tumor-associated antigens. Specific amino acid sequences of the complementarity-determining regions (CDRs) in the antigen-binding fragment (Fab) of these antibodies have potential as lung cancer biomarkers. Detection and identification of CDRs by mass spectrometry can significantly be improved by reduction of the complexity of the immunoglobulin molecule. Our aim was to molecular dissect IgG into κ and λ fragments to reduce the complexity and thereby identify substantially more CDRs than by just total Fab isolation. We purified Fab, Fab-κ, Fab-λ, κ and λ light chains from serum from 10 stage I lung adenocarcinoma patients and 10 matched controls from the current and former smokers. After purification, the immunoglobulin fragments were enzymatically digested and measured by high-resolution mass spectrometry. Finally, we compared the number of CDRs identified in these immunoglobulin fragments with that in the Fab fragments. Twice as many CDRs were identified when Fab-κ, Fab-λ, κ and λ (3330) were combined than in the Fab fraction (1663) alone. The number of CDRs and κ:λ ratio was statistically similar in both cases and controls. Molecular dissection of IgG identifies significantly more CDRs, which increases the likelihood of finding lung cancer-related CDR sequences.

Coupling mammalian cell surface display with somatic hypermutation for the discovery and maturation of human antibodies

A novel approach has been developed for the isolation and maturation of human antibodies that replicates key features of the adaptive immune system by coupling in vitro somatic hypermutation (SHM) with mammalian cell display. SHM is dependent on the action of the B cell specific enzyme, activation-induced cytidine deaminase (AID), and can be replicated in non-B cells through expression of recombinant AID. A library of human antibodies, based on germline V-gene segments with recombined human regions was used to isolate low-affinity antibodies to human β nerve growth factor (hβNGF). These antibodies, initially naïve to SHM, were subjected to AID-directed SHM in vitro and selected using the same mammalian cell display system, as illustrated by the maturation of one of the antibodies to low pM K(D). This approach overcomes many of the previous limitations of mammalian cell display, enabling direct selection and maturation of antibodies as full-length, glycosylated IgGs.

On-chip screening of experimental conditions for the synthesis of noble-metal nanostructures with different morphologies

The applications of nanostructures critically depend on their morphologies. Although significant progress has been made in the chemical synthesis of nanostructures with a variety of different morphologies, it is still highly desired to develop an approach that allows one to quickly identify the best set of parameters for nanostructure syntheses. Herein, an on-chip approach to the rapid screening of experimental conditions pivotal to the production of nanostructures with different morphologies is reported. The key component of this approach is an array of reactors containing solutions with a one- or two-dimensional gradient in reagent concentration, pH value, or reaction temperature. In the proof-of-concept experiments, the parameters needed for the production of Au and Pd nanostructures with various morphologies are quickly identified. In principle, this approach can be extended to other systems for rapid screening and optimization of experimental conditions involved in the syntheses of different types of nanostructures.

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.

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.

A multi-scale model for correlation in B cell VDJ usage of zebrafish

The zebrafish (Danio rerio) is one of the model animals used for the study of immunology because the dynamics in the adaptive immune system of zebrafish are similar to that in higher animals. In this work, we built a multi-scale model to simulate the dynamics of B cells in the primary and secondary immune responses of zebrafish. We use this model to explain the reported correlation between VDJ usage of B cell repertoires in individual zebrafish. We use a delay ordinary differential equation (ODE) system to model the immune responses in the 6-month lifespan of a zebrafish. This mean field theory gives the number of high-affinity B cells as a function of time during an infection. The sequences of those B cells are then taken from a distribution calculated by a 'microscopic' random energy model. This generalized NK model shows that mature B cells specific to one antigen largely possess a single VDJ recombination. The model allows first-principle calculation of the probability, p, that two zebrafish responding to the same antigen will select the same VDJ recombination. This probability p increases with the B cell population size and the B cell selection intensity. The probability p decreases with the B cell hypermutation rate. The multi-scale model predicts correlations in the immune system of the zebrafish that are highly similar to that from experiment.

High-resolution description of antibody heavy-chain repertoires in humans

Antibodies' protective, pathological, and therapeutic properties result from their considerable diversity. This diversity is almost limitless in potential, but actual diversity is still poorly understood. Here we use deep sequencing to characterize the diversity of the heavy-chain CDR3 region, the most important contributor to antibody binding specificity, and the constituent V, D, and J segments that comprise it. We find that, during the stepwise D-J and then V-DJ recombination events, the choice of D and J segments exert some bias on each other; however, we find the choice of the V segment is essentially independent of both. V, D, and J segments are utilized with different frequencies, resulting in a highly skewed representation of VDJ combinations in the repertoire. Nevertheless, the pattern of segment usage was almost identical between two different individuals. The pattern of V, D, and J segment usage and recombination was insufficient to explain overlap that was observed between the two individuals' CDR3 repertoires. Finally, we find that while there are a near-infinite number of heavy-chain CDR3s in principle, there are about 3–9 million in the blood of an adult human being.

Pumping-induced perturbation of flow in microfluidic channels and its implications for on-chip cell culture

We study the rate of response to changes in the rate of flow and the perturbations in flow in polydimethylsiloxane (PDMS) microfluidic chips that are subjected to several common flow-control systems. We find that the flow rate of liquid delivered from a syringe pump equipped with a glass syringe responds faster to the changes in the conditions of flow than the same liquid delivered from a plastic syringe; and the rate of flow delivered from compressed air responds faster than that from a glass syringe. We discover that the rate of flow that is driven by a syringe pump and regulated by an integrated pneumatic valve responds even faster, but this flow-control method is characterized by large perturbations. We also examine the possible effects of these large perturbations on NIH 3T3 cells in microfluidic channels and find that they could cause the detachment of NIH 3T3 cells in the microchannels.

Detecting selection in immunoglobulin sequences

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

Systems analysis of adaptive immunity by utilization of high-throughput technologies

A new generation of high-throughput technologies for quantitative and clonal analysis of adaptive immune responses have been developed. Functional analysis of lymphocyte populations has been accomplished via microfluidic assay systems. Additionally, lymphocyte receptor repertoires have been characterized on proteomic and genomic levels with multiplexed protein microarrays and high-throughput DNA sequencing. These tools are providing an unprecedented level of information depth on the distribution of adaptive immune cell (B and T cell) functionalities and repertoires, which develop upon activation following vaccination, pathogenic infection, or in disease states. These various high-throughput technologies have unlocked the potential to transform immunology into an information-rich science that will enable rapid expansion of the field of experimental systems immunology.

Deep mutational scanning: assessing protein function on a massive scale

Analysis of protein mutants is an effective means to understand their function. Protein display is an approach that allows large numbers of mutants of a protein to be selected based on their activity, but only a handful with maximal activity have been traditionally identified for subsequent functional analysis. However, the recent application of high-throughput sequencing (HTS) to protein display and selection has enabled simultaneous assessment of the function of hundreds of thousands of mutants that span the activity range from high to low. Such deep mutational scanning approaches are rapid and inexpensive with the potential for broad utility. In this review, we discuss the emergence of deep mutational scanning, the challenges associated with its use and some of its exciting applications.

Determinism and stochasticity during maturation of the zebrafish antibody repertoire

It is thought that the adaptive immune system of immature organisms follows a more deterministic program of antibody creation than is found in adults. We used high-throughput sequencing to characterize the diversifying antibody repertoire in zebrafish over five developmental time points. We found that the immune system begins in a highly stereotyped state with preferential use of a small number of V (variable) D (diverse) J (joining) gene segment combinations, but that this stereotypy decreases dramatically as the zebrafish mature, with many of the top VDJ combinations observed in 2-wk-old zebrafish virtually disappearing by 1 mo. However, we discovered that, in the primary repertoire, there are strong correlations in VDJ use that increase with zebrafish maturity, suggesting that VDJ recombination involves a level of deterministic programming that is unexpected. This stereotypy is masked by the complex diversification processes of antibody maturation; the variation and lack of correlation in full repertoires between individuals appears to be derived from randomness in clonal expansion during the affinity maturation process. These data provide a window into the mechanisms of VDJ recombination and diversity creation and allow us to better understand how the adaptive immune system achieves diversity.

The whole-organism heavy chain B cell repertoire from Zebrafish self-organizes into distinct network features

BACKGROUND

The adaptive immune system is based on selected populations of molecularly distinct individual B and T cell clones. However, it has not been possible to characterize these clones in a comprehensive and informatics manner to date; attempts have been limited by the number of cells in the adaptive immune system and an inability to quantify them. Recently, using the Zebrafish (ZF) Danio rerio as a model organism and parallel sequencing as the quantifying technology, Weinstein et al. overcame this major hurdle and quantified the entire heavy chain B-cell repertoire in ZF. Here, we present a novel network analysis of the data from the Weinstein group, providing new insights into the network structure of the B-cell repertoire.

RESULTS

Using a collection of computational methods, the IgM sequences from 14 fish were analyzed. This analysis demonstrated that the B-cell repertoire of the ZF is structured along similar lines to those previously detected in limited parts of the human B-cell immune system. The analysis confirms the validity of the global data and the evolutionary placement of the ZF based on known sequence motifs. Recombination events in the repertoire were quantified, and demonstrated a lack of shared recombined V, J groups across fish. Nevertheless, it was demonstrated that a similar network architecture is shared among fish. However, the network analysis identified two distinct populations within the group; these findings are compatible with the occurrence of an immune response in a subset of the fish. The emerging connectivity network was demonstrated and quantified, and mutation drifts within the groups were characterized. Dissection of sequence data revealed common network features of the B-cell repertoire as well as individual differences.

CONCLUSION

The ZF B-cell repertoire reveals an underlying order that is compatible with self-organization representing every portion of the sequence-based network. This pattern varies in individual specimens, perhaps as a response to an immune challenge. However, a sequence-non-specific network that maintains a common architecture of sequence diversity was detected.The common feature among different individuals can be captured by the network architecture and characteristics, rather than specific clones. We believe that further study of the dynamics of this network could provide insight into modes of operation of the immune system.

Sequencing antibody repertoires: the next generation

Genomic studies have been revolutionized by the use of next generation sequencing (NGS) that delivers huge amounts of sequence information in a short span of time. The number of applications for NGS is rapidly expanding and significantly transforming many areas of life sciences. The field of antibody research and discovery is no exception. Several recent studies have harnessed the power of NGS for analyzing natural or synthetic immunoglobulin repertoires with unprecedented resolution and exploring alternative paths for antibody discovery. Thus, appreciating, and then exploiting, these advances is essential for staying at the edge of antibody innovation.

Factors important in evolutionary shaping of immunoglobulin gene loci

Background

The extraordinary diversity characterizing the antibody repertoire is generated by both evolution and lymphocyte development. Much of this diversity is due to the existence of immunoglobulin (Ig) variable region gene segment libraries, which were diversified during evolution and, in higher vertebrates, are used in generating the combinatorial diversity of antibody genes. The aim of the present study was to address the following questions: What evolutionary parameters affect the size and structure of gene libraries? Are the number of genes in libraries of contemporary species, and the corresponding gene locus structure, a random result of evolutionary history, or have these properties been optimized with respect to individual or population fitness? If a larger number of genes or different genome structures do not increase the fitness, then the current structure is probably optimized.

Results

We used a simulation of variable region gene library evolution. We measured the effect of different parameters on gene library size and diversity, and the corresponding fitness. We found compensating relationships between parameters, which optimized Ig library size and diversity.

Conclusions

We conclude that contemporary species' Ig libraries have been optimized by evolution in terms of Ig sequence lengths, the number and diversity of Ig genes, and antibody-antigen affinities.

An antibody-based biomarker discovery method by mass spectrometry sequencing of complementarity determining regions

Autoantibodies are increasingly used as biomarkers in the detection of autoimmune disorders and cancer. Disease specific antibodies are generally detected by their binding to specific antigens. As an alternative approach, we propose to identify specific complementarity determining regions (CDR) of IgG that relate to an autoimmune disorder or cancer instead of the specific antigen(s). In this manuscript, we tested the technical feasibility to detect and identify CDRs of specific antibodies by mass spectrometry. We used a commercial pooled IgG preparation as well as purified serum IgG fractions that were spiked with different amounts of a fully human monoclonal antibody (adalimumab). These samples were enzymatically digested and analyzed by nanoLC Orbitrap mass spectrometry. In these samples, we were able to identify peptides derived from the CDRs of adalimumab. These peptides could be detected at an amount of 110 attomole, 5 orders of magnitude lower than the total IgG concentration in these samples. Using higher energy collision induced dissociation (HCD) fragmentation and subsequent de novo sequencing, we could successfully identify 50% of the detectable CDR peptides of adalimumab. In addition, we demonstrated that an affinity purification with anti-dinitrophenol (DNP) monoclonal antibody enhanced anti-DNP derived CDR detection in a serum IgG background. In conclusion, specific CDR peptides could be detected and sequenced at relatively low levels (attomole-femtomole range) which should allow the detection of clinically relevant CDR peptides in patient samples.

Sequencing and quantifying IgG fragments and antigen-binding regions by mass spectrometry

In cancer and autoimmune diseases, immunoglobulins with a specific molecular signature that could potentially be used as diagnostic or prognostic markers are released into body fluids. An immunomics approach based on this phenomenon relies on the ability to identify the specific amino acid sequences of the complementarity-determining regions (CDR) of these immunoglobulins, which in turn depends on the level of accuracy, resolution, and sensitivity that can be achieved by advanced mass spectrometry. Reproducible isolation and sequencing of antibody fragments (e.g., Fab) by high-resolution mass spectrometry (MS) from seven healthy donors revealed 43 217 MS signals: 225 could be associated with CDR1 peptides, 513 with CDR2 peptides, and 19 with CDR3 peptides. Seventeen percent of the 43 217 MS signals did not overlap between the seven donors. The Fab isolation method used is reproducible and fast, with a high yield. It provides only one Fab sample fraction for subsequent characterization by high-resolution MS. In 17% and 4% of these seven healthy donors, qualitative (presence/absence) and quantitative (intensity) differences in Fab fragments could be demonstrated, respectively. From these results, we conclude that the identification of a CDR signature as biomarker for autoimmune diseases and cancer without prior knowledge of the antigen is feasible.

Targets of somatic hypermutation within immunoglobulin light chain genes in zebrafish

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

Targeted resequencing of candidate genes using selector probes

Targeted genome enrichment is a powerful tool for making use of the massive throughput of novel DNA-sequencing instruments. We herein present a simple and scalable protocol for multiplex amplification of target regions based on the Selector technique. The updated version exhibits improved coverage and compatibility with next-generation-sequencing (NGS) library-construction procedures for shotgun sequencing with NGS platforms. To demonstrate the performance of the technique, all 501 exons from 28 genes frequently involved in cancer were enriched for and sequenced in specimens derived from cell lines and tumor biopsies. DNA from both fresh frozen and formalin-fixed paraffin-embedded biopsies were analyzed and 94% specificity and 98% coverage of the targeted region was achieved. Reproducibility between replicates was high (R² = 0, 98) and readily enabled detection of copy-number variations. The procedure can be carried out in <24 h and does not require any dedicated instrumentation.

Model for comparative analysis of antigen receptor repertoires

In modern molecular biology one of the standard ways of analyzing a vertebrate immune system is to sequence and compare the counts of specific antigen receptor clones (either immunoglobulins or T-cell receptors) derived from various tissues under different experimental or clinical conditions. The resulting statistical challenges are difficult and do not fit readily into the standard statistical framework of contingency tables primarily due to the serious under-sampling of the receptor populations. This under-sampling is caused, on one hand, by the extreme diversity of antigen receptor repertoires maintained by the immune system and, on the other, by the high cost and labor intensity of the receptor data collection process. In most of the recent immunological literature the differences across antigen receptor populations are examined via non-parametric statistical measures of the species overlap and diversity borrowed from ecological studies. While this approach is robust in a wide range of situations, it seems to provide little insight into the underlying clonal size distribution and the overall mechanism differentiating the receptor populations. As a possible alternative, the current paper presents a parametric method that adjusts for the data under-sampling as well as provides a unifying approach to a simultaneous comparison of multiple receptor groups by means of the modern statistical tools of unsupervised learning. The parametric model is based on a flexible multivariate Poisson-lognormal distribution and is seen to be a natural generalization of the univariate Poisson-lognormal models used in the ecological studies of biodiversity patterns. The procedure for evaluating a model's fit is described along with the public domain software developed to perform the necessary diagnostics. The model-driven analysis is seen to compare favorably vis a vis traditional methods when applied to the data from T-cell receptors in transgenic mice populations.

The origins of vertebrate adaptive immunity

Adaptive immunity is mediated through numerous genetic and cellular processes that generate favourable somatic variants of antigen-binding receptors under evolutionary selection pressure by pathogens and other factors. Advances in our understanding of immunity in mammals and other model organisms are revealing the underlying basis and complexity of this remarkable system. Although the evolution of adaptive immunity has been thought to occur by the acquisition of novel molecular capabilities, an increasing amount of information from new model systems suggest that co-option and redirection of pre-existing systems are the main source of innovation. We combine evidence from a wide range of organisms to obtain an integrated view of the origins and patterns of divergence in adaptive immunity.

By-passing in vitro screening--next generation sequencing technologies applied to antibody display and in silico candidate selection

In recent years, unprecedented DNA sequencing capacity provided by next generation sequencing (NGS) has revolutionized genomic research. Combining the Illumina sequencing platform and a scFv library designed to confine diversity to both CDR3, >1.9 × 10(7) sequences have been generated. This approach allowed for in depth analysis of the library's diversity, provided sequence information on virtually all scFv during selection for binding to two targets and a global view of these enrichment processes. Using the most frequent heavy chain CDR3 sequences, primers were designed to rescue scFv from the third selection round. Identification, based on sequence frequency, retrieved the most potent scFv and valuable candidates that were missed using classical in vitro screening. Thus, by combining NGS with display technologies, laborious and time consuming upfront screening can be by-passed or complemented and valuable insights into the selection process can be obtained to improve library design and understanding of antibody repertoires.

Regulation of the B cell receptor repertoire and self-reactivity by BAFF

The TNF-family cytokine BAFF (BLyS) promotes B lymphocyte survival and is overexpressed in individuals with systemic lupus erythematosus and Sjögren's Syndrome. BAFF can rescue anergic autoreactive B cells from death, but only when competition from nonautoreactive B cells is lacking. Yet, high BAFF levels promote autoantibody formation in individuals possessing diverse B cells. To better understand how excess BAFF promotes autoimmunity in a polyclonal immune system, Ig L chain usage was analyzed in 3H9 site-directed IgH chain transgenic mice, whose B cells recognize DNA and chromatin when they express certain endogenous L chains. BAFF levels were manipulated in 3H9 mice by introducing transgenes expressing either BAFF or its natural inhibitor ΔBAFF. B cells in BAFF/3H9 mice were elevated in number, used a broad L chain repertoire, including L chains generating high-affinity autoreactivity, and produced abundant autoantibodies. Comparison of spleen and lymph node B cells suggested that highly autoreactive B cells were expanded. By contrast, ΔBAFF/3H9 mice had reduced B cell numbers with a repertoire similar to that of 3H9 mice, but lacking usage of a subset of Vκ genes. The results show that limiting BAFF signaling only slightly selects against higher affinity autoreactive B cells, whereas its overexpression leads to broad tolerance escape and positive selection of autoreactive cells. The results have positive implications for the clinical use of BAFF-depleting therapy.

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.

Monoclonal antibodies isolated without screening by analyzing the variable-gene repertoire of plasma cells

Isolation of antigen-specific monoclonal antibodies (mAbs) and antibody fragments relies on high-throughput screening of immortalized B cells or recombinant antibody libraries. We bypassed the screening step by using high-throughput DNA sequencing and bioinformatic analysis to mine antibody variable region (V)-gene repertoires from bone marrow plasma cells (BMPC) of immunized mice. BMPCs, which cannot be immortalized, produce the vast majority of circulating antibodies. We found that the V-gene repertoire of BMPCs becomes highly polarized after immunization, with the most abundant sequences represented at frequencies between approximately 1% and >10% of the total repertoire. We paired the most abundant variable heavy (V(H)) and variable light (V(L)) genes based on their relative frequencies, reconstructed them using automated gene synthesis, and expressed recombinant antibodies in bacteria or mammalian cells. Antibodies generated in this manner from six mice, each immunized with one of three antigens were overwhelmingly antigen specific (21/27 or 78%). Those generated from a mouse with high serum titers had nanomolar binding affinities.

Rapid construction and characterization of synthetic antibody libraries without DNA amplification

We report on a simple method to rapidly generate very large libraries of genes encoding mutant proteins without the use of DNA amplification, and the application of this methodology in the construction of synthetic immunoglobulin variable heavy (V(H)) and light (V(kappa)) libraries. Four high quality, chemically synthesized polynucleotides (90-140 bases) were annealed and extended using T4 DNA polymerase. Following electroporation, >10(9) transformants could be synthesized within 1 day. Fusion to beta-lactamase and selection on ampicillin resulted in 3.7 x 10(8) V(H) and 6.9 x 10(8) V(kappa) clones highly enriched for full-length, in-frame genes. High-throughput 454 DNA sequencing of >250,000 V(H) and V(kappa) genes from the pre- and post-selection libraries revealed that, in addition to the expected reduction in reading-frame shifts and stop codons, selection for functional expression also resulted in a statistical decrease in the cysteine content. Apart from these differences, there was a good agreement between the expected and actual diversity, indicating that neither oligonucleotide synthesis nor biological constrains due to protein synthesis of V(H)/V(kappa)-beta-lactamase fusions introduce biases in the amino acid composition of the randomized regions. This methodology can be employed for the rapid construction of highly diverse libraries with the near elimination of PCR errors in invariant regions.