High-throughput sequencing of the zebrafish antibody repertoire

Tracking global changes induced in the CD4 T-cell receptor repertoire by immunization with a complex antigen using short stretches of CDR3 protein sequence

Motivation: The clonal theory of adaptive immunity proposes that immunological responses are encoded by increases in the frequency of lymphocytes carrying antigen-specific receptors. In this study, we measure the frequency of different T-cell receptors (TcR) in CD4 + T cell populations of mice immunized with a complex antigen, killed Mycobacterium tuberculosis, using high throughput parallel sequencing of the TcRβ chain. Our initial hypothesis that immunization would induce repertoire convergence proved to be incorrect, and therefore an alternative approach was developed that allows accurate stratification of TcR repertoires and provides novel insights into the nature of CD4 + T-cell receptor recognition.

Results: To track the changes induced by immunization within this heterogeneous repertoire, the sequence data were classified by counting the frequency of different clusters of short (3 or 4) continuous stretches of amino acids within the antigen binding complementarity determining region 3 (CDR3) repertoire of different mice. Both unsupervised (hierarchical clustering) and supervised (support vector machine) analyses of these different distributions of sequence clusters differentiated between immunized and unimmunized mice with 100% efficiency. The CD4 + TcR repertoires of mice 5 and 14 days postimmunization were clearly different from that of unimmunized mice but were not distinguishable from each other. However, the repertoires of mice 60 days postimmunization were distinct both from naive mice and the day 5/14 animals. Our results reinforce the remarkable diversity of the TcR repertoire, resulting in many diverse private TcRs contributing to the T-cell response even in genetically identical mice responding to the same antigen. However, specific motifs defined by short stretches of amino acids within the CDR3 region may determine TcR specificity and define a new approach to TcR sequence classification.

Availability and implementation: The analysis was implemented in R and Python, and source code can be found in Supplementary Data.

Contact:b.chain@ucl.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Systematic characterization and comparative analysis of the rabbit immunoglobulin repertoire

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

Quantifying selection in immune receptor repertoires

The efficient recognition of pathogens by the adaptive immune system relies on the diversity of receptors displayed at the surface of immune cells. T-cell receptor diversity results from an initial random DNA editing process, called VDJ recombination, followed by functional selection of cells according to the interaction of their surface receptors with self and foreign antigenic peptides. Using high-throughput sequence data from the β-chain of human T-cell receptors, we infer factors that quantify the overall effect of selection on the elements of receptor sequence composition: the V and J gene choice and the length and amino acid composition of the variable region. We find a significant correlation between biases induced by VDJ recombination and our inferred selection factors together with a reduction of diversity during selection. Both effects suggest that natural selection acting on the recombination process has anticipated the selection pressures experienced during somatic evolution. The inferred selection factors differ little between donors or between naive and memory repertoires. The number of sequences shared between donors is well-predicted by our model, indicating a stochastic origin of such public sequences. Our approach is based on a probabilistic maximum likelihood method, which is necessary to disentangle the effects of selection from biases inherent in the recombination process.

Antibody repertoire deep sequencing reveals antigen-independent selection in maturing B cells

Antibody repertoires are known to be shaped by selection for antigen binding. Unexpectedly, we now show that selection also acts on a non-antigen-binding antibody region: the heavy-chain variable (VH)-encoded "elbow" between variable and constant domains. By sequencing 2.8 million recombined heavy-chain genes from immature and mature B-cell subsets in mice, we demonstrate a striking gradient in VH gene use as pre-B cells mature into follicular and then into marginal zone B cells. Cells whose antibodies use VH genes that encode a more flexible elbow are more likely to mature. This effect is distinct from, and exceeds in magnitude, previously described maturation-associated changes in heavy-chain complementarity determining region 3, a key antigen-binding region, which arise from junctional diversity rather than differential VH gene use. Thus, deep sequencing reveals a previously unidentified mode of B-cell selection.

Immunoglobulin gene repertoire diversification and selection in the stomach - from gastritis to gastric lymphomas

Chronic gastritis is characterized by gastric mucosal inflammation due to autoimmune responses or infection, frequently with Helicobacter pylori. Gastritis with H. pylori background can cause gastric mucosa-associated lymphoid tissue lymphoma (MALT-L), which sometimes further transforms into diffuse large B-cell lymphoma (DLBCL). However, gastric DLBCL can also be initiated de novo. The mechanisms underlying transformation into DLBCL are not completely understood. We analyzed immunoglobulin repertoires and clonal trees to investigate whether and how immunoglobulin gene repertoires, clonal diversification, and selection in gastritis, gastric MALT-L, and DLBCL differ from each other and from normal responses. The two gastritis types (positive or negative for H. pylori) had similarly diverse repertoires. MALT-L dominant clones (defined as the largest clones in each sample) presented higher diversification and longer mutational histories compared with all other conditions. DLBCL dominant clones displayed lower clonal diversification, suggesting the transforming events are triggered by similar responses in different patients. These results are surprising, as we expected to find similarities between the dominant clones of gastritis and MALT-L and between those of MALT-L and DLBCL.

Characterization of immunoglobulin by mass spectrometry with applications for the clinical laboratory

Studies monitoring immunoglobulin (Ig) antigen specificity have brought to light key Ig biomarkers for immunity, autoimmunity, cancer detection, and immune system function evaluation. A fundamentally new approach to the detection of Igs based on the primary structure of the Ig is beginning to emerge in the literature. This approach has only become feasible in light of advances in proteomics and rapid improvements in mass spectrometry (MS). Driven primarily by the development of Ig pharmaceuticals, Ig MS-based proteomic methods are revealing structural features which were previously unavailable with other characterization techniques. The task of adapting these techniques to clinical chemistry is in its infancy, but these methods have the potential to dramatically alter testing for Ig biomarkers. The purpose of this article is to review the advances that have been made in proteomic characterization of Igs by MS and the early attempts to apply these methods to clinical samples.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Studying the antibody repertoire after vaccination: practical applications

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

Comprehensive evaluation and optimization of amplicon library preparation methods for high-throughput antibody sequencing

High-throughput sequencing (HTS) of antibody repertoire libraries has become a powerful tool in the field of systems immunology. However, numerous sources of bias in HTS workflows may affect the obtained antibody repertoire data. A crucial step in antibody library preparation is the addition of short platform-specific nucleotide adapter sequences. As of yet, the impact of the method of adapter addition on experimental library preparation and the resulting antibody repertoire HTS datasets has not been thoroughly investigated. Therefore, we compared three standard library preparation methods by performing Illumina HTS on antibody variable heavy genes from murine antibody-secreting cells. Clonal overlap and rank statistics demonstrated that the investigated methods produced equivalent HTS datasets. PCR-based methods were experimentally superior to ligation with respect to speed, efficiency, and practicality. Finally, using a two-step PCR based method we established a protocol for antibody repertoire library generation, beginning from inputs as low as 1 ng of total RNA. In summary, this study represents a major advance towards a standardized experimental framework for antibody HTS, thus opening up the potential for systems-based, cross-experiment meta-analyses of antibody repertoires.

Peptides from the variable region of specific antibodies are shared among lung cancer patients

Late diagnosis of lung cancer is still the main reason for high mortality rates in lung cancer. Lung cancer is a heterogeneous disease which induces an immune response to different tumor antigens. Several methods for searching autoantibodies have been described that are based on known purified antigen panels. The aim of our study is to find evidence that parts of the antigen-binding-domain of antibodies are shared among lung cancer patients. This was investigated by a novel approach based on sequencing antigen-binding-fragments (Fab) of immunoglobulins using proteomic techniques without the need of previously known antigen panels. From serum of 93 participants of the NELSON trial IgG was isolated and subsequently digested into Fab and Fc. Fab was purified from the digested mixture by SDS-PAGE. The Fab containing gel-bands were excised, tryptic digested and measured on a nano-LC-Orbitrap-Mass-spectrometry system. Multivariate analysis of the mass spectrometry data by linear canonical discriminant analysis combined with stepwise logistic regression resulted in a 12-antibody-peptide model which was able to distinguish lung cancer patients from controls in a high risk population with a sensitivity of 84% and specificity of 90%. With our Fab-purification combined Orbitrap-mass-spectrometry approach, we found peptides from the variable-parts of antibodies which are shared among lung cancer patients.

Proteomic identification of monoclonal antibodies from serum

Characterizing the in vivo dynamics of the polyclonal antibody repertoire in serum, such as that which might arise in response to stimulation with an antigen, is difficult due to the presence of many highly similar immunoglobulin proteins, each specified by distinct B lymphocytes. These challenges have precluded the use of conventional mass spectrometry for antibody identification based on peptide mass spectral matches to a genomic reference database. Recently, progress has been made using bottom-up analysis of serum antibodies by nanoflow liquid chromatography/high-resolution tandem mass spectrometry combined with a sample-specific antibody sequence database generated by high-throughput sequencing of individual B cell immunoglobulin variable domains (V genes). Here, we describe how intrinsic features of antibody primary structure, most notably the interspersed segments of variable and conserved amino acid sequences, generate recurring patterns in the corresponding peptide mass spectra of V gene peptides, greatly complicating the assignment of correct sequences to mass spectral data. We show that the standard method of decoy-based error modeling fails to account for the error introduced by these highly similar sequences, leading to a significant underestimation of the false discovery rate. Because of these effects, antibody-derived peptide mass spectra require increased stringency in their interpretation. The use of filters based on the mean precursor ion mass accuracy of peptide-spectrum matches is shown to be particularly effective in distinguishing between "true" and "false" identifications. These findings highlight important caveats associated with the use of standard database search and error-modeling methods with nonstandard data sets and custom sequence databases.

T cells translate individual, quantal activation into collective, analog cytokine responses via time-integrated feedbacks

Variability within isogenic T cell populations yields heterogeneous ‘local’ signaling responses to shared antigenic stimuli, but responding clones may communicate ‘global’ antigen load through paracrine messengers, such as cytokines. Such coordination of individual cell responses within multicellular populations is critical for accurate collective reactions to shared environmental cues. However, cytokine production may saturate as a function of antigen input, or be dominated by the precursor frequency of antigen-specific T cells. Surprisingly, we found that T cells scale their collective output of IL-2 to total antigen input over a large dynamic range, independently of population size. Through experimental quantitation and computational modeling, we demonstrate that this scaling is enforced by an inhibitory cross-talk between antigen and IL-2 signaling, and a nonlinear acceleration of IL-2 secretion per cell. Our study reveals how time-integration of these regulatory loops within individual cell signaling generates scaled collective responses and can be leveraged for immune monitoring.

DOI:http://dx.doi.org/10.7554/eLife.01944.001

The cells of the immune system face the challenge of removing viruses and other pathogens without endangering healthy tissues. Cells called T cells plays a variety of roles in the immune response: some T cells directly destroy infected cells, some recruit other cells called phagocytes to the site of infection, and some release small proteins called cytokines. These cytokines help cells to communicate with other cells and, therefore, to tailor the overall immune responses to deal with a particular pathogen.

It is known that mammals are capable of adjusting the T cell response to match the overall severity of an infection. However, it is not clear how individual T cells coordinate their seemingly binary response—they are either activated when they recognize a pathogen, or they are not activated—into a response at the collective cell level that can be varied continuously over a wide range of values.

Here, Tkach et al. show that T cell populations match their production of the cytokine interleukin 2 (IL-2) to the abundance of antigens—molecules released by the pathogen—over an unexpectedly large range of concentrations. Through a combination of experimental and computational analyses, Tkach et al. identified two novel IL-2 feedback loops that help to generate the correct quantity of cytokine, irrespective of the total number of T cells. Furthermore, this model can be used to estimate antigen quantities within diseased tissues. The work of Tkach et al. illustrates the potential of feedback integration in cell signalling and gene regulation as a mechanism to allow cellular populations to respond to environmental stimuli in a graded, collective fashion.

DOI:http://dx.doi.org/10.7554/eLife.01944.002

High-resolution antibody dynamics of vaccine-induced immune responses

The adaptive immune system confers protection by generating a diverse repertoire of antibody receptors that are rapidly expanded and contracted in response to specific targets. Next-generation DNA sequencing now provides the opportunity to survey this complex and vast repertoire. In the present work, we describe a set of tools for the analysis of antibody repertoires and their application to elucidating the dynamics of the response to viral vaccination in human volunteers. By analyzing data from 38 separate blood samples across 2 y, we found that the use of the germ-line library of V and J segments is conserved between individuals over time. Surprisingly, there appeared to be no correlation between the use level of a particular VJ combination and degree of expansion. We found the antibody RNA repertoire in each volunteer to be highly dynamic, with each individual displaying qualitatively different response dynamics. By using combinatorial phage display, we screened selected VH genes paired with their corresponding VL library for affinity against the vaccine antigens. Altogether, this work presents an additional set of tools for profiling the human antibody repertoire and demonstrates characterization of the fast repertoire dynamics through time in multiple individuals responding to an immune challenge.

The functional repertoire of rabbit antibodies and antibody discovery via next-generation sequencing

To gain insight into the functional antibody repertoire of rabbits, the VH and VL repertoires of bone marrow (BM) and spleen (SP) of a naïve New Zealand White rabbit (NZW; Oryctolagus cuniculus) and that of lymphocytes collected from a NZW rabbit immunized (IM) with a 16-mer peptide were deep-sequenced. Two closely related genes, IGHV1S40 (VH1a3) and IGHV1S45 (VH4), were found to dominate (~90%) the VH repertoire of BM and SP, whereas, IGHV1S69 (VH1a1) contributed significantly (~40%) to IM. BM and SP antibodies recombined predominantly with IGHJ4. A significant proportion (~30%) of IM sequences recombined with IGHJ2. The VK repertoire was encoded by nine IGKV genes recombined with one IGKJ gene, IGKJ1. No significant bias in the VK repertoire of the BM, SP and IM samples was observed. The complementarity-determining region (CDR)-H3 and -L3 length distributions were similar in the three samples following a Gaussian curve with average length of 12.2 ± 2.4 and 11.1 ± 1.1 amino acids, respectively. The amino acid composition of the predominant CDR-H3 and -L3 loop lengths was similar to that of humans and mice, rich in Tyr, Gly, Ser and, in some specific positions, Asp. The average number of mutations along the IGHV/KV genes was similar in BM, SP and IM; close to 12 and 15 mutations for VH and VL, respectively. A monoclonal antibody specific for the peptide used as immunogen was obtained from the IM rabbit. The CDR-H3 sequence was found in 1,559 of 61,728 (2.5%) sequences, at position 10, in the rank order of the CDR-H3 frequencies. The CDR-L3 was found in 24 of 11,215 (0.2%) sequences, ranking 102. No match was found in the BM and SP samples, indicating positive selection for the hybridoma sequence. Altogether, these findings lay foundations for engineering of rabbit V regions to enhance their potential as therapeutics, i.e., design of strategies for selection of specific rabbit V regions from NGS data mining, humanization and design of libraries for affinity maturation campaigns.

Profiling the IgOme: meeting the challenge

The entire repertoire of antibodies in our serum, the IgOme, is a historical record of our past experiences and a reflection of our immune status at any given moment. Understanding the dynamics of the IgOme and how the diversity and specificities of serum antibodies change in response to disease and maintenance of homeostasis can directly impact the ability to design and develop novel vaccines, diagnostics and therapeutics. Here we review both direct and indirect methodologies that are being developed to map the complexity and specificities of the antibodies in polyclonal serum - the IgOme.

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

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

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

Lineage structure of the human antibody repertoire in response to influenza vaccination

The human antibody repertoire is one of the most important defenses against infectious disease, and the development of vaccines has enabled the conferral of targeted protection to specific pathogens. However, there are many challenges to measuring and analyzing the immunoglobulin sequence repertoire, including that each B cell's genome encodes a distinct antibody sequence, that the antibody repertoire changes over time, and the high similarity between antibody sequences. We have addressed these challenges by using high-throughput long read sequencing to perform immunogenomic characterization of expressed human antibody repertoires in the context of influenza vaccination. Informatic analysis of 5 million antibody heavy chain sequences from healthy individuals allowed us to perform global characterizations of isotype distributions, determine the lineage structure of the repertoire, and measure age- and antigen-related mutational activity. Our analysis of the clonal structure and mutational distribution of individuals' repertoires shows that elderly subjects have a decreased number of lineages but an increased prevaccination mutation load in their repertoire and that some of these subjects have an oligoclonal character to their repertoire in which the diversity of the lineages is greatly reduced relative to younger subjects. We have thus shown that global analysis of the immune system's clonal structure provides direct insight into the effects of vaccination and provides a detailed molecular portrait of age-related effects.

The immune system as a biomonitor: explorations in innate and adaptive immunity

The human immune system has a highly complex, multi-layered structure which has evolved to detect and respond to changes in the internal microenvironment of the body. Recognition occurs at the molecular or submolecular scale, via classical reversible receptor-ligand interactions, and can lead to a response with great sensitivity and speed. Remarkably, recognition is coupled to memory, such that responses are modulated by events which occurred years or even decades before. Although the immune system in general responds differently and more vigorously to stimuli entering the body from the outside (e.g. infections), this is an emergent property of the system: many of the recognition molecules themselves have no inherent bias towards external stimuli (non-self) but also bind targets found within the body (self). It is quite clear that the immune response registers pathophysiological changes in general. Cancer, wounding and chronic tissue injury are some obvious examples. Against this background, the immune system 'state' tracks the internal processes of the body, and is likely to encode information regarding both current and past disease processes. Moreover, the distributed nature of most immune responses (e.g. typically involving lymphoid tissue, non-lymphoid tissue, bone marrow, blood, extracellular interstitial spaces, etc.) means that many of the changes associated with immune responses are manifested systemically, and specifically can be detected in blood. This provides a very convenient route to sampling immune cells. We consider two different and complementary ways of querying the human immune 'state' using high-dimensional genomic screening methodologies, and discuss the potentials of these approaches and some of the technological and computational challenges to be overcome.

Convergent antibody signatures in human dengue

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

Error analysis of deep sequencing of phage libraries: peptides censored in sequencing

Next-generation sequencing techniques empower selection of ligands from phage-display libraries because they can detect low abundant clones and quantify changes in the copy numbers of clones without excessive selection rounds. Identification of errors in deep sequencing data is the most critical step in this process because these techniques have error rates >1%. Mechanisms that yield errors in Illumina and other techniques have been proposed, but no reports to date describe error analysis in phage libraries. Our paper focuses on error analysis of 7-mer peptide libraries sequenced by Illumina method. Low theoretical complexity of this phage library, as compared to complexity of long genetic reads and genomes, allowed us to describe this library using convenient linear vector and operator framework. We describe a phage library as N × 1 frequency vector n = ||ni||, where ni is the copy number of the ith sequence and N is the theoretical diversity, that is, the total number of all possible sequences. Any manipulation to the library is an operator acting on n. Selection, amplification, or sequencing could be described as a product of a N × N matrix and a stochastic sampling operator (Sa). The latter is a random diagonal matrix that describes sampling of a library. In this paper, we focus on the properties of Sa and use them to define the sequencing operator (Seq). Sequencing without any bias and errors is Seq = Sa IN, where IN is a N × N unity matrix. Any bias in sequencing changes IN to a nonunity matrix. We identified a diagonal censorship matrix (CEN), which describes elimination or statistically significant downsampling, of specific reads during the sequencing process.

Deep sequencing and Circos analyses of antibody libraries reveal antigen-driven selection of Ig VH genes during HIV-1 infection

The vast diversity of antibody repertoires is largely attributed to heavy chain (V(H)) recombination of variable (V), diversity (D) and joining (J) gene segments. We used 454 sequencing information of the variable domains of the antibody heavy chain repertoires from neonates, normal adults and an HIV-1-infected individual, to analyze, with Circos software, the VDJ pairing patterns at birth, adulthood and a time-dependent response to HIV-1 infection. Our comparative analyses of the Ig VDJ repertoires from these libraries indicated that, from birth to adulthood, VDJ recombination patterns remain the same with some slight changes, whereas some V(H) families are selected and preferentially expressed after long-term infection with HIV-1. We also demonstrated that the immune system responds to HIV-1 chronic infection by selectively expanding certain HV families in an attempt to combat infection. Our findings may have implications for understanding immune responses in pathology as well as for development of new therapeutics and vaccines.

The past, present, and future of immune repertoire biology - the rise of next-generation repertoire analysis

T and B cell repertoires are collections of lymphocytes, each characterized by its antigen-specific receptor. We review here classical technologies and analysis strategies developed to assess immunoglobulin (IG) and T cell receptor (TR) repertoire diversity, and describe recent advances in the field. First, we describe the broad range of available methodological tools developed in the past decades, each of which answering different questions and showing complementarity for progressive identification of the level of repertoire alterations: global overview of the diversity by flow cytometry, IG repertoire descriptions at the protein level for the identification of IG reactivities, IG/TR CDR3 spectratyping strategies, and related molecular quantification or dynamics of T/B cell differentiation. Additionally, we introduce the recent technological advances in molecular biology tools allowing deeper analysis of IG/TR diversity by next-generation sequencing (NGS), offering systematic and comprehensive sequencing of IG/TR transcripts in a short amount of time. NGS provides several angles of analysis such as clonotype frequency, CDR3 diversity, CDR3 sequence analysis, V allele identification with a quantitative dimension, therefore requiring high-throughput analysis tools development. In this line, we discuss the recent efforts made for nomenclature standardization and ontology development. We then present the variety of available statistical analysis and modeling approaches developed with regards to the various levels of diversity analysis, and reveal the increasing sophistication of those modeling approaches. To conclude, we provide some examples of recent mathematical modeling strategies and perspectives that illustrate the active rise of a "next-generation" of repertoire analysis.

IgM Repertoire Biodiversity is Reduced in HIV-1 Infection and Systemic Lupus Erythematosus

Background: HIV-1 infection or systemic lupus erythematosus (SLE) disrupt B cell homeostasis, reduce memory B cells, and impair function of IgG and IgM antibodies.

Objective: To determine how disturbances in B cell populations producing polyclonal antibodies relate to the IgM repertoire, the IgM transcriptome in health and disease was explored at the complementarity determining region 3 (CDRH3) sequence level.

Methods: 454-deep pyrosequencing in combination with a novel analysis pipeline was applied to define populations of IGHM CDRH3 sequences based on absence or presence of somatic hypermutations (SHM) in peripheral blood B cells.

Results: HIV or SLE subjects have reduced biodiversity within their IGHM transcriptome compared to healthy subjects, mainly due to a significant decrease in the number of unique combinations of alleles, although recombination machinery was intact. While major differences between sequences without or with SHM occurred among all groups, IGHD and IGHJ allele use, CDRH3 length distribution, or generation of SHM were similar among study cohorts. Antiretroviral therapy failed to normalize IGHM biodiversity in HIV-infected individuals. All subjects had a low frequency of allelic combinations within the IGHM repertoire similar to known broadly neutralizing HIV-1 antibodies.

Conclusion: Polyclonal expansion would decrease overall IgM biodiversity independent of other mechanisms for development of the B cell repertoire. Applying deep sequencing as a strategy to follow development of the IgM repertoire in health and disease provides a novel molecular assessment of multiple points along the B cell differentiation pathway that is highly sensitive for detecting perturbations within the repertoire at the population level.

Origin and evolution of adaptive immunity

The evolutionary emergence of vertebrates was accompanied by major morphological and functional innovations, including the development of an adaptive immune system. Vertebrate adaptive immunity is based on the clonal expression of somatically diversifying antigen receptors on lymphocytes. This is a common feature of both the jawless and jawed vertebrates , although these two groups of extant vertebrates employ structurally different types of antigen receptors and principal mechanisms for their somatic diversification . These observations suggest that the common vertebrate ancestor must have already possessed a complex immune system, including B- and T-like lymphocyte lineages and primary lymphoid organs, such as the thymus, but possibly lacked the facilities for somatic diversification of antigen receptors. Interestingly, memory formation, previously considered to be a defining feature of adaptive immunity, also occurs in the context of innate immune responses and can even be observed in unicellular organisms, attesting to the convergent evolutionary history of distinct aspects of adaptive immunity.

IgG variable region and VH CDR3 diversity in unimmunized mice analyzed by massively parallel sequencing

Most antigen-specific mouse antibodies have been derived by hybridoma technology, predominantly through use of the Balb/c strain. Much of the Balb/c germline repertoire of variable genes (V regions) is known. However, there is little information about the background expressed repertoire of IgG antibodies in mice, which reflects the baseline against which antigen-specific antibodies are generated through immunization. To assess this baseline repertoire, RNA was isolated from splenic B-cells enriched for expression of IgG from three mice. The RNA was individually amplified with three distinct PCR primer sets for comprehensive recovery of the heavy and light chain variable regions. Each PCR product was independently subjected to deep sequencing using 454 pyro-sequencing technology and analysed for redundancy, open reading frame, germline representation, and CDR3 sequence of the heavy chain variable region (VH CDR3) within and across the primer sets and mice. A highly skewed abundance of heavy and light chain variable gene usage was observed for all three primers in all three mice. While showing considerable overlap, there were differences among these profiles indicative of primer bias and animal-to-animal variation. VH CDR3 sequences were likewise highly skewed indicating that the heavy chain genes profiles substantially reflected individual antibodies. This observation was confirmed through analysis of randomly selected complete heavy chain variable sequences. However, there was very little redundancy in VH CDR3 sequences across the different mice. We conclude that the background IgG repertoire in young, unimmunized mice is highly skewed within individual mice and is diverse among them, a pattern similar to that observed in highly immunized mice.

Single-cell based high-throughput sequencing of full-length immunoglobulin heavy and light chain genes

Single-cell PCR and sequencing of full-length Ig heavy (Igh) and Igk and Igl light chain genes is a powerful tool to measure the diversity of antibody repertoires and allows the functional assessment of B-cell responses through direct Ig gene cloning and the generation of recombinant mAbs. However, the current methodology is not high-throughput compatible. Here we developed a two-dimensional bar-coded primer matrix to combine Igh and Igk/Igl chain gene single-cell PCR with next-generation sequencing for the parallel analysis of the antibody repertoire of over 46 000 individual B cells. Our approach provides full-length Igh and corresponding Igk/Igl chain gene-sequence information and permits the accurate correction of sequencing errors by consensus building. The use of indexed cell sorting for the isolation of single B cells enables the integration of flow cytometry and Ig gene sequence information. The strategy is fully compatible with established protocols for direct antibody gene cloning and expression and therefore advances over previously described high-throughput approaches to assess antibody repertoires at the single-cell level.

Genomic and systems approaches to translational biomarker discovery in immunological diseases

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

Expression of immunoglobulin receptors with distinctive features indicating antigen selection by marginal zone B cells from human spleen

Marginal zone (MZ) B cells, identified as surface (s)IgM(high)sIgD(low)CD23(low/-)CD21(+)CD38(-) B cells, were purified from human spleens, and the features of their V(D)J gene rearrangements were investigated and compared with those of germinal center (GC), follicular mantle (FM) and switched memory (SM) B cells. Most MZ B cells were CD27(+) and exhibited somatic hypermutations (SHM), although to a lower extent than SM B cells. Moreover, among MZ B-cell rearrangements, recurrent sequences were observed, some of which displayed intraclonal diversification. The same diversifying sequences were detected in very low numbers in GC and FM B cells and only when a highly sensitive, gene-specific polymerase chain reaction was used. This result indicates that MZ B cells could expand and diversify in situ and also suggested the presence of a number of activation-induced cytidine deaminase (AID)-expressing B cells in the MZ. The notion of antigen-driven expansion/selection in situ is further supported by the VH CDR3 features of MZ B cells with highly conserved amino acids at specific positions and by the finding of shared ("stereotyped") sequences in two different spleens. Collectively, the data are consistent with the notion that MZ B cells are a special subset selected by in situ antigenic stimuli.

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.

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.

Cerebrospinal-fluid-derived immunoglobulin G of different multiple sclerosis patients shares mutated sequences in complementarity determining regions

B lymphocytes play a pivotal role in multiple sclerosis pathology, possibly via both antibody-dependent and -independent pathways. Intrathecal immunoglobulin G in multiple sclerosis is produced by clonally expanded B-cell populations. Recent studies indicate that the complementarity determining regions of immunoglobulins specific for certain antigens are frequently shared between different individuals. In this study, our main objective was to identify specific proteomic profiles of mutated complementarity determining regions of immunoglobulin G present in multiple sclerosis patients but absent in healthy controls. To achieve this objective, we purified immunoglobulin G from the cerebrospinal fluid of 29 multiple sclerosis patients and 30 healthy controls and separated the corresponding heavy and light chains via SDS-PAGE. Subsequently, bands were excised, trypsinized, and measured with high-resolution mass spectrometry. We sequenced 841 heavy and 771 light chain variable region peptides. We observed 24 heavy and 26 light chain complementarity determining regions that were solely present in a number of multiple sclerosis patients. Using stringent criteria for the identification of common peptides, we found five complementarity determining regions shared in three or more patients and not in controls. Interestingly, one complementarity determining region with a single mutation was found in six patients. Additionally, one other patient carrying a similar complementarity determining region with another mutation was observed. In addition, we found a skew in the κ-to-λ ratio and in the usage of certain variable heavy regions that was previously observed at the transcriptome level. At the protein level, cerebrospinal fluid immunoglobulin G shares common characteristics in the antigen binding region among different multiple sclerosis patients. The indication of a shared fingerprint may indicate common antigens for B-cell activation.

Genetic measurement of memory B-cell recall using antibody repertoire sequencing

Annual influenza vaccinations aim to protect against seasonal infections, and vaccine strain compositions are updated every year. This protection is based on antibodies that are produced by either newly activated or memory B cells recalled from previous encounters with influenza vaccination or infection. The extent to which the B-cell repertoire responds to vaccination and recalls antibodies has so far not been analyzed at a genetic level--which is to say, at the level of antibody sequences. Here, we developed a consensus read sequencing approach that incorporates unique barcode labels on each starting RNA molecule. These labels allow one to combine multiple sequencing reads covering the same RNA molecule to reduce the error rate to a desired level, and they also enable accurate quantification of RNA and isotype levels. We validated this approach and analyzed the differential response of the antibody repertoire to live-attenuated or trivalent-inactivated influenza vaccination. Additionally, we analyzed the antibody repertoire in response to repeated yearly vaccinations with trivalent-inactivated influenza vaccination. We found antibody sequences that were present in both years, providing a direct genetic measurement of B-cell recall.

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.

Deep mutational scanning of an antibody against epidermal growth factor receptor using mammalian cell display and massively parallel pyrosequencing

We developed a method for deep mutational scanning of antibody complementarity-determining regions (CDRs) that can determine in parallel the effect of every possible single amino acid CDR substitution on antigen binding. The method uses libraries of full length IgGs containing more than 1000 CDR point mutations displayed on mammalian cells, sorted by flow cytometry into subpopulations based on antigen affinity and analyzed by massively parallel pyrosequencing. Higher, lower and neutral affinity mutations are identified by their enrichment or depletion in the FACS subpopulations. We applied this method to a humanized version of the anti-epidermal growth factor receptor antibody cetuximab, generated a near comprehensive data set for 1060 point mutations that recapitulates previously determined structural and mutational data for these CDRs and identified 67 point mutations that increase affinity. The large-scale, comprehensive sequence-function data sets generated by this method should have broad utility for engineering properties such as antibody affinity and specificity and may advance theoretical understanding of antibody-antigen recognition.

Deep mRNA sequencing analysis to capture the transcriptome landscape of zebrafish embryos and larvae

Transcriptome analysis is a powerful tool to obtain large amount genome-scale gene expression profiles. Despite its extensive usage to diverse biological problems in the last decade, transcriptomic researches approaching the zebrafish embryonic development have been very limited. Several recent studies have made great progress in this direction, yet the large gap still exists, especially regarding to the transcriptome dynamics of embryonic stages from early gastrulation onwards. Here, we present a comprehensive analysis about the transcriptomes of 9 different stages covering 7 major periods (cleavage, blastula, gastrula, segmentation, pharyngula, hatching and early larval stage) in zebrafish development, by recruiting the RNA-sequencing technology. We detected the expression for at least 24,065 genes in at least one of the 9 stages. We identified 16,130 genes that were significantly differentially expressed between stages and were subsequently classified into six clusters. Each revealed gene cluster had distinct expression patterns and characteristic functional pathways, providing a framework for the understanding of the developmental transcriptome dynamics. Over 4000 genes were identified as preferentially expressed in one of the stages, which could be of high relevance to stage-specific developmental and molecular events. Among the 68 transcription factor families active during development, most had enhanced average expression levels and thus might be crucial for embryogenesis, whereas the inactivation of the other families was likely required by the activation of the zygotic genome. We discussed our RNA-seq data together with previous findings about the Wnt signaling pathway and some other genes with known functions, to show how our data could be used to advance our understanding about these developmental functional elements. Our study provides ample information for further study about the molecular and cellular mechanisms underlying vertebrate development.

Single-cell analysis and sorting using droplet-based microfluidics

We present a droplet-based microfluidics protocol for high-throughput analysis and sorting of single cells. Compartmentalization of single cells in droplets enables the analysis of proteins released from or secreted by cells, thereby overcoming one of the major limitations of traditional flow cytometry and fluorescence-activated cell sorting. As an example of this approach, we detail a binding assay for detecting antibodies secreted from single mouse hybridoma cells. Secreted antibodies are detected after only 15 min by co-compartmentalizing single mouse hybridoma cells, a fluorescent probe and single beads coated with anti-mouse IgG antibodies in 50-pl droplets. The beads capture the secreted antibodies and, when the captured antibodies bind to the probe, the fluorescence becomes localized on the beads, generating a clearly distinguishable fluorescence signal that enables droplet sorting at ∼200 Hz as well as cell enrichment. The microfluidic system described is easily adapted for screening other intracellular, cell-surface or secreted proteins and for quantifying catalytic or regulatory activities. In order to screen ∼1 million cells, the microfluidic operations require 2-6 h; the entire process, including preparation of microfluidic devices and mammalian cells, requires 5-7 d.

Rheumatoid arthritis synovial tissue harbours dominant B-cell and plasma-cell clones associated with autoreactivity

OBJECTIVE

To identify potential autoreactive B-cell and plasma-cell clones by quantitatively analysing the complete human B-cell receptor (BCR) repertoire in synovium and peripheral blood in early and established rheumatoid arthritis (RA).

METHODS

The BCR repertoire was screened in synovium and blood of six patients with early RA (ERA) (<6 months) and six with established RA (ESRA) (>20 months). In two patients, the repertoires in different joints were compared. Repertoires were analysed by next-generation sequencing from mRNA, generating >10 000 BCR heavy-chain sequence reads per sample. For each clone, the degree of expansion was calculated as the percentage of the total number of reads encoding the specific clonal sequence. Clones with a frequency ≥ 0.5% were considered dominant.

RESULTS

Multiple dominant clones were found in inflamed synovium but hardly any in blood. Within an individual patient, the same dominant clones were detected in different joints. The majority of the synovial clones were class-switched; however, the fraction of clones that expressed IgM was higher in ESRA than ERA patients. Dominant synovial clones showed autoreactive features: in ERA in particular the clones were enriched for immunoglobulin heavy chain gene segment V4-34 (IGHV4-34) and showed longer CDR3 lengths. Dominant synovial clones that did not encode IGHV4-34 also had longer CDR3s than peripheral blood.

CONCLUSIONS

In RA, the synovium forms a niche where expanded--potentially autoreactive--B cells and plasma cells reside. The inflamed target tissue, especially in the earliest phase of disease, seems to be the most promising compartment for studying autoreactive cells.

Wrestling with the repertoire: the promise and perils of next generation sequencing for antigen receptors

Next generation sequencing technologies are revolutionizing the study of immune repertoires. These methods provide a previously unimaginable amount of sequence data, unfortunately accompanied by numerous challenges associated with error correction and interpretation that remain to be solved. For antigen receptors, these challenges will require dedicated solutions beyond those developed for genome sequencing, which may differ depending on the sequencing technology used and the purpose of the experiment. Many investigators are developing such methods, based on different sequencing platforms, but critical details of protocol and performance are proprietary. The field will move forward when these methods are shared and standardized, and when the accuracy, sensitivity and reproducibility of various sequencing, and analytic methods are evaluated using standardized samples in comparative experiments.

Mining the antibodyome for HIV-1-neutralizing antibodies with next-generation sequencing and phylogenetic pairing of heavy/light chains

Next-generation sequencing of antibody transcripts from HIV-1-infected individuals with broadly neutralizing antibodies could provide an efficient means for identifying somatic variants and characterizing their lineages. Here, we used 454 pyrosequencing and identity/divergence grid sampling to analyze heavy- and light-chain sequences from donor N152, the source of the broadly neutralizing antibody 10E8. We identified variants with up to 28% difference in amino acid sequence. Heavy- and light-chain phylogenetic trees of identified 10E8 variants displayed similar architectures, and 10E8 variants reconstituted from matched and unmatched phylogenetic branches displayed significantly lower autoreactivity when matched. To test the generality of phylogenetic pairing, we analyzed donor International AIDS Vaccine Initiative 84, the source of antibodies PGT141-145. Heavy- and light-chain phylogenetic trees of PGT141-145 somatic variants also displayed remarkably similar architectures; in this case, branch pairings could be anchored by known PGT141-145 antibodies. Altogether, our findings suggest that phylogenetic matching of heavy and light chains can provide a means to approximate natural pairings.

Deep sequencing of phage display libraries to support antibody discovery

The use of next generation sequencing (NGS) for the analysis of antibody sequences both in phage display libraries and during in vitro selection processes has become increasingly popular in the last few years. Here, our methods developed for DNA preparation, sequencing and data analysis are presented. A key parameter has also been to develop new software designed for high throughput antibody sequence analysis that is used in combination with publicly available tools. As an example of our methods, we provide data from the extensive analysis of five scFv libraries generated using different heavy chain CDR3 diversification strategies. The results not only confirm that the library designs were correct but also reveal differences in quality not easily identified by standard DNA sequencing approaches. The very large number of reads permits extensive sequence coverage after the selection process. Furthermore, as samples can be multiplexed, costs decrease and more information is gained per NGS run. Using examples of results obtained post phage display selections against two antigens, frequency and clustering analysis identified novel antibody fragments that were then shown to be specific for the target antigen. In summary, the methods described here demonstrate how NGS analysis enhances quality control of complex antibody libraries as well as facilitates the antibody discovery process.

Diversity-oriented approaches for interrogating T-cell receptor repertoire, ligand recognition, and function

Molecular diversity lies at the heart of adaptive immunity. T-cell receptors and peptide-major histocompatibility complex molecules utilize and rely upon an enormous degree of diversity at the levels of genetics, chemistry, and structure to engage one another and carry out their functions. This high level of diversity complicates the systematic study of important aspects of T-cell biology, but recent technical advances have allowed for the ability to study diversity in a comprehensive manner. In this review, we assess insights gained into T-cell receptor function and biology from our increasingly precise ability to assess the T-cell repertoire as a whole or to perturb individual receptors with engineered reagents. We conclude with a perspective on a new class of high-affinity, non-stimulatory peptide ligands we have recently discovered using diversity-oriented techniques that challenges notions for how we think about T-cell receptor signaling.

The astonishing diversity of Ig classes and B cell repertoires in teleost fish

With lymphoid tissue anatomy different than mammals, and diverse adaptations to all aquatic environments, fish constitute a fascinating group of vertebrate to study the biology of B cell repertoires in a comparative perspective. Fish B lymphocytes express immunoglobulin (Ig) on their surface and secrete antigen-specific antibodies in response to immune challenges. Three antibody classes have been identified in fish, namely IgM, IgD, and IgT, while IgG, IgA, and IgE are absent. IgM and IgD have been found in all fish species analyzed, and thus seem to be primordial antibody classes. IgM and IgD are normally co-expressed from the same mRNA through alternative splicing, as in mammals. Tetrameric IgM is the main antibody class found in serum. Some species of fish also have IgT, which seems to exist only in fish and is specialized in mucosal immunity. IgM/IgD and IgT are expressed by two different sub-populations of B cells. The tools available to investigate B cell responses at the cellular level in fish are limited, but the progress of fish genomics has started to unravel a rich diversity of IgH and immunoglobulin light chain locus organization, which might be related to the succession of genome remodelings that occurred during fish evolution. Moreover, the development of deep sequencing techniques has allowed the investigation of the global features of the expressed fish B cell repertoires in zebrafish and rainbow trout, in steady state or after infection. This review provides a description of the organization of fish Ig loci, with a particular emphasis on their heterogeneity between species, and presents recent data on the structure of the expressed Ig repertoire in healthy and infected fish.

High-throughput sequencing of the paired human immunoglobulin heavy and light chain repertoire

Each B-cell receptor consists of a pair of heavy and light chains. High-throughput sequencing can identify large numbers of heavy- and light-chain variable regions (V(H) and V(L)) in a given B-cell repertoire, but information about endogenous pairing of heavy and light chains is lost after bulk lysis of B-cell populations. Here we describe a way to retain this pairing information. In our approach, single B cells (>5 × 10(4) capacity per experiment) are deposited in a high-density microwell plate (125 pl/well) and lysed in situ. mRNA is then captured on magnetic beads, reverse transcribed and amplified by emulsion V(H):V(L) linkage PCR. The linked transcripts are analyzed by Illumina high-throughput sequencing. We validated the fidelity of V(H):V(L) pairs identified by this approach and used the method to sequence the repertoire of three human cell subsets-peripheral blood IgG(+) B cells, peripheral plasmablasts isolated after tetanus toxoid immunization and memory B cells isolated after seasonal influenza vaccination.

Fractal organization of the human T cell repertoire in health and after stem cell transplantation

T cell repertoire diversity is generated in part by recombination of variable (V), diversity (D), and joining (J) segments in the T cell receptor β (TCR) locus. T cell clonal frequency distribution determined by high-throughput sequencing of TCR β in 10 stem cell transplantation (SCT) donors revealed a fractal, self-similar frequency distribution of unique TCR bearing clones with respect to V, D, and J segment usage in the T cell repertoire of these individuals. Further, ranking of T cell clones by frequency of gene segment usage in the observed sequences revealed an ordered distribution of dominant clones conforming to a power law, with a fractal dimension of 1.6 and 1.8 in TCR β DJ and VDJ containing clones in healthy stem cell donors. This self-similar distribution was perturbed in the recipients after SCT, with patients demonstrating a lower level of complexity in their TCR repertoire at day 100 followed by a modest improvement by 1 year post-SCT. A large shift was observed in the frequency distribution of the dominant T cell clones compared to the donor, with fewer than one third of the VDJ-containing clones shared in the top 4 ranks. In conclusion, the normal T cell repertoire is highly ordered with a TCR gene segment usage that results in a fractal self-similar motif of pattern repetition across levels of organization. Fractal analysis of high-throughput TCR β sequencing data provides a comprehensive measure of immune reconstitution after SCT.

Teleost fish mount complex clonal IgM and IgT responses in spleen upon systemic viral infection

Upon infection, B-lymphocytes expressing antibodies specific for the intruding pathogen develop clonal responses triggered by pathogen recognition via the B-cell receptor. The constant region of antibodies produced by such responding clones dictates their functional properties. In teleost fish, the clonal structure of B-cell responses and the respective contribution of the three isotypes IgM, IgD and IgT remain unknown. The expression of IgM and IgT are mutually exclusive, leading to the existence of two B-cell subsets expressing either both IgM and IgD or only IgT. Here, we undertook a comprehensive analysis of the variable heavy chain (VH) domain repertoires of the IgM, IgD and IgT in spleen of homozygous isogenic rainbow trout (Onchorhynchus mykiss) before, and after challenge with a rhabdovirus, the Viral Hemorrhagic Septicemia Virus (VHSV), using CDR3-length spectratyping and pyrosequencing of immunoglobulin (Ig) transcripts. In healthy fish, we observed distinct repertoires for IgM, IgD and IgT, respectively, with a few amplified μ and τ junctions, suggesting the presence of IgM- and IgT-secreting cells in the spleen. In infected animals, we detected complex and highly diverse IgM responses involving all VH subgroups, and dominated by a few large public and private clones. A lower number of robust clonal responses involving only a few VH were detected for the mucosal IgT, indicating that both IgM(+) and IgT(+) spleen B cells responded to systemic infection but at different degrees. In contrast, the IgD response to the infection was faint. Although fish IgD and IgT present different structural features and evolutionary origin compared to mammalian IgD and IgA, respectively, their implication in the B-cell response evokes these mouse and human counterparts. Thus, it appears that the general properties of antibody responses were already in place in common ancestors of fish and mammals, and were globally conserved during evolution with possible functional convergences.

High-resolution deep sequencing reveals biodiversity, population structure, and persistence of HIV-1 quasispecies within host ecosystems

Background

Deep sequencing provides the basis for analysis of biodiversity of taxonomically similar organisms in an environment. While extensively applied to microbiome studies, population genetics studies of viruses are limited. To define the scope of HIV-1 population biodiversity within infected individuals, a suite of phylogenetic and population genetic algorithms was applied to HIV-1 envelope hypervariable domain 3 (Env V3) within peripheral blood mononuclear cells from a group of perinatally HIV-1 subtype B infected, therapy-naïve children.

Results

Biodiversity of HIV-1 Env V3 quasispecies ranged from about 70 to 270 unique sequence clusters across individuals. Viral population structure was organized into a limited number of clusters that included the dominant variants combined with multiple clusters of low frequency variants. Next generation viral quasispecies evolved from low frequency variants at earlier time points through multiple non-synonymous changes in lineages within the evolutionary landscape. Minor V3 variants detected as long as four years after infection co-localized in phylogenetic reconstructions with early transmitting viruses or with subsequent plasma virus circulating two years later.

Conclusions

Deep sequencing defines HIV-1 population complexity and structure, reveals the ebb and flow of dominant and rare viral variants in the host ecosystem, and identifies an evolutionary record of low-frequency cell-associated viral V3 variants that persist for years. Bioinformatics pipeline developed for HIV-1 can be applied for biodiversity studies of virome populations in human, animal, or plant ecosystems.

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".