454 antibody sequencing - error characterization and correction

Bioinformatic Analysis of Natively Paired VH:VL Antibody Repertoires for Antibody Discovery

Next-generation DNA sequencing (NGS) of human antibody repertoires has been extensively implemented to discover novel antibody drugs, to analyze B-cell developmental features, and to investigate antibody responses to infectious diseases and vaccination. Because the antibody repertoire encoded by human B cells is highly diverse, NGS analyses of antibody genes have provided a new window into understanding antibody responses for basic immunology, biopharmaceutical drug discovery, and immunotherapy. However, many antibody discovery protocols analyze the heavy and light chains separately due to the short-read nature of most NGS technologies, whereas paired heavy and light chain data are required for complete antibody characterization. Here, we describe a computational workflow to process millions of paired antibody heavy and light chain DNA sequence reads using the Illumina MiSeq 2x300 NGS platform. In this workflow, we describe raw NGS read processing and initial quality filtering, the annotation and assembly of antibody clonotypes relating to paired heavy and light chain antibody lineages, and the generation of complete heavy+light consensus sequences for the downstream cloning and expression of human antibody proteins.

Adaptive immune receptor repertoires, an overview of this exciting field

The adaptive immune response in jawed vertebrates relies on the huge diversity and specificity of the B cell and T cell antigen receptors, the immunoglobulins (IG) or antibodies and the T cell receptors (TR), respectively. The high level of diversity has represented a barrier to a comprehensive analysis of the adaptive immune response before the emergence of high-throughput sequencing (HTS) technologies. The size and complexity of HTS data requires the generation of novel computational and analytical approaches, which are transforming how the adaptive immune responses are deciphered to understand the clonal dynamics and properties of antigen-specific B and T cells in response to different kind of antigens. This exciting and rapidly evolving field is not only impacting human and clinical immunology but also comparative immunology. We are now closer to understanding the evolution of adaptive immune response in jawed vertebrates. This review provides an overview about classical and current strategies developed to assess the IG/TR diversity and their applications in basic and clinical immunology.

Combined Influence of B-Cell Receptor Rearrangement and Somatic Hypermutation on B-Cell Class-Switch Fate in Health and in Chronic Lymphocytic Leukemia

A diverse B-cell receptor (BCR) repertoire is required to bind a wide range of antigens. BCRs are generated through genetic recombination and can be diversified through somatic hypermutation (SHM) or class-switch recombination (CSR). Patterns of repertoire diversity can vary substantially between different health conditions. We use isotype-resolved BCR sequencing to compare B-cell evolution and class-switch fate in healthy individuals and in patients with chronic lymphocytic leukemia (CLL). We show that the patterns of SHM and CSR in B-cells from healthy individuals are distinct from CLL. We identify distinct properties of clonal expansion that lead to the generation of antibodies of different classes in healthy, malignant, and non-malignant CLL BCR repertoires. We further demonstrate that BCR diversity is affected by relationships between antibody variable and constant regions leading to isotype-specific signatures of variable gene usage. This study provides powerful insights into the mechanisms underlying the evolution of the adaptive immune responses in health and their aberration during disease.

History, applications, and challenges of immune repertoire research

The diversity of T and B cells in terms of their receptor sequences is huge in the vertebrate's immune system and provides broad protection against the vast diversity of pathogens. Immune repertoire is defined as the sum of T cell receptors and B cell receptors (also named immunoglobulin) that makes the organism's adaptive immune system. Before the emergence of high-throughput sequencing, the studies on immune repertoire were limited by the underdeveloped methodologies, since it was impossible to capture the whole picture by the low-throughput tools. The massive paralleled sequencing technology suits perfectly the researches on immune repertoire. In this article, we review the history of immune repertoire studies, in terms of technologies and research applications. Particularly, we discuss several aspects of challenges in this field and highlight the efforts to develop potential solutions, in the era of high-throughput sequencing of the immune repertoire.

Deciphering antigen-responding antibody repertoires by using next-generation sequencing and confirming them through antibody-gene synthesis

Vast diversity and high specificity of antigen recognition by antibodies are hallmarks of the acquired immune system. Although the molecular mechanisms that yield the extremely large antibody repertoires are precisely understood, comprehensive description of the global antibody repertoire generated in individual bodies has been hindered by the lack of powerful measures. To obtain holistic understanding of the antibody-repertoire space, we used next-generation sequencing (NGS) to analyze the deep profiles of naive and antigen-responding repertoires of the IgM, IgG1, and IgG2c classes formed in individual mice. The overall landscapes of naive IgM repertoires were almost the same for each mouse, whereas those of IgG1 and IgG2c differed considerably among naive individuals. Next, we immunized mice with a model antigen, nitrophenol (NP)-hapten linked to chicken γ-globulin (CGG) carrier, and compared the antigen-responding repertoires in individual mice. To extract the complete antigen response, we developed an intelligible method for detecting common components of antigen-responding repertoires. The major responding antibodies were IGHV1-72/IGHD1-1/IGHJ2 for NP-hapten and IGHV9-3/IGHD3-1/IGHJ2 for CGG-carrier protein. The antigen-binding specificities of the identified antibodies were confirmed through ELISA after antibody-gene synthesis and expression of the corresponding NGS reads. Thus, we deciphered antigen-responding antibody repertoires by inclusively analyzing the antibody-repertoire space generated in individual bodies by using NGS, which avoided inadvertent omission of key antibody repertoires.

The analysis of clonal expansions in normal and autoimmune B cell repertoires

Clones are the fundamental building blocks of immune repertoires. The number of different clones relates to the diversity of the repertoire, whereas their size and sequence diversity are linked to selective pressures. Selective pressures act both between clones and within different sequence variants of a clone. Understanding how clonal selection shapes the immune repertoire is one of the most basic questions in all of immunology. But how are individual clones defined? Here we discuss different approaches for defining clones, starting with how antibodies are diversified during different stages of B cell development. Next, we discuss how clones are defined using different experimental methods. We focus on high-throughput sequencing datasets, and the computational challenges and opportunities that these data have for mining the antibody repertoire landscape. We discuss methods that visualize sequence variants within the same clone and allow us to consider collections of shared mutations to determine which sequences share a common ancestry. Finally, we comment on features of frequently encountered expanded B cell clones that may be of particular interest in the setting of autoimmunity and other chronic conditions.

Immunoglobulin class-switched B cells form an active immune axis between CNS and periphery in multiple sclerosis

In multiple sclerosis (MS), lymphocyte--in particular B cell--transit between the central nervous system (CNS) and periphery may contribute to the maintenance of active disease. Clonally related B cells exist in the cerebrospinal fluid (CSF) and peripheral blood (PB) of MS patients; however, it remains unclear which subpopulations of the highly diverse peripheral B cell compartment share antigen specificity with intrathecal B cell repertoires and whether their antigen stimulation occurs on both sides of the blood-brain barrier. To address these questions, we combined flow cytometric sorting of PB B cell subsets with deep immune repertoire sequencing of CSF and PB B cells. Immunoglobulin (IgM and IgG) heavy chain variable (VH) region repertoires of five PB B cell subsets from MS patients were compared with their CSF Ig-VH transcriptomes. In six of eight patients, we identified peripheral CD27(+)IgD(-) memory B cells, CD27(hi)CD38(hi) plasma cells/plasmablasts, or CD27(-)IgD(-) B cells that had an immune connection to the CNS compartment. Pinpointing Ig class-switched B cells as key component of the immune axis thought to contribute to ongoing MS disease activity strengthens the rationale of current B cell-targeting therapeutic strategies and may lead to more targeted approaches.

Discovery of novel candidate therapeutics and diagnostics based on engineered human antibody domains

The smallest independently folded antibody fragments, the domains, are emerging as promising scaffolds for candidate therapeutics and diagnostics that bind specifically targets of interest. The discovery of such binders is based on several technologies including structure-based design and generation of libraries of mutants displayed on phage or yeast, next-generation sequencing for diversity analysis, panning and screening of the libraries, affinity maturation of selected binders, and their expression, purification, and characterization for specific binding, function, and aggregation propensity. In this review, we describe these technologies as applied for the generation of engineered antibody domains (eAds), especially those derived from the human immunoglobulin heavy chain variable region (VH) and the second domain of IgG1 heavy chain constant region (CH2) as potential candidate therapeutics and diagnostics, and discuss examples of eAds against HIV-1 and cancer-related proteins.

Analysis of plant microbe interactions in the era of next generation sequencing technologies

Next generation sequencing (NGS) technologies have impressively accelerated research in biological science during the last years by enabling the production of large volumes of sequence data to a drastically lower price per base, compared to traditional sequencing methods. The recent and ongoing developments in the field allow addressing research questions in plant-microbe biology that were not conceivable just a few years ago. The present review provides an overview of NGS technologies and their usefulness for the analysis of microorganisms that live in association with plants. Possible limitations of the different sequencing systems, in particular sources of errors and bias, are critically discussed and methods are disclosed that help to overcome these shortcomings. A focus will be on the application of NGS methods in metagenomic studies, including the analysis of microbial communities by amplicon sequencing, which can be considered as a targeted metagenomic approach. Different applications of NGS technologies are exemplified by selected research articles that address the biology of the plant associated microbiota to demonstrate the worth of the new methods.

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.

Routine performance and errors of 454 HLA exon sequencing in diagnostics

Background

Next-generation sequencing (NGS) has changed genomics significantly. More and more applications strive for sequencing with different platforms. Now, in 2012, after a decade of development and evolution, NGS has been accepted for a variety of research fields. Determination of sequencing errors is essential in order to follow next-generation sequencing beyond research use only. This study describes the overall 454 system performance of using multiple GS Junior runs with an in-house established and validated diagnostic assay for human leukocyte antigen (HLA) exon sequencing. Based on this data, we extracted, evaluated and characterized errors and variants of 60 HLA loci per run with respect to their adjacencies.

Results

We determined an overall error rate of 0.18% in a total of 118,484,408 bases. 31.3% of all reads analyzed (n=349,503) contain one or more errors. The largest group are deletions that account for 50% of the errors. Incorrect bases are not distributed equally along sequences and tend to be more frequent at sequence ends. Certain sequence positions in the middle or at the beginning of the read accumulate errors. Typically, the corresponding quality score at the actual error position is lower than the adjacent scores.

Conclusions

Here we present the first error assessment in a human next-generation sequencing diagnostics assay in an amplicon sequencing approach. Improvements of sequence quality and error rate that have been made over the years are evident and it is shown that both have now reached a level where diagnostic applications become feasible. Our presented data are better than previously published error rates and we can confirm and quantify the often described relation of homopolymers and errors. Nevertheless, a certain depth of coverage is needed, in particular with challenging areas of the sequencing target. Furthermore, the usage of error correcting tools is not essential but might contribute towards the capacity and efficiency of a sequencing run.

Human monoclonal antibodies as candidate therapeutics against emerging viruses and HIV-1

More than 40 monoclonal antibodies (mAbs) have been approved for a number of disease indications with only one of these (Synagis) - for a viral disease, and not for therapy but for prevention. However, in the last decade novel potent mAbs have been discovered and characterized with potential as therapeutics against viruses of major importance for public health and biosecurity including Hendra virus (HeV), Nipah virus (NiV), severe acute respiratory syndrome coronavirus (SARS-CoV), Ebola virus (EBOV), West Nile virus (WNV), influenza virus (IFV) and human immunodeficiency virus type 1 (HIV-1). Here, we review such mAbs with an emphasis on antibodies of human origin, and highlight recent results as well as technologies and mechanisms related to their potential as therapeutics.

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.

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

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

Origin, diversity, and maturation of human antiviral antibodies analyzed by high-throughput sequencing

Our understanding of how antibodies are generated and function could help develop effective vaccines and antibody-based therapeutics against viruses such as HIV-1, SARS coronavirus (SARS CoV), and Hendra and Nipah viruses (henipaviruses). Although broadly neutralizing antibodies (bnAbs) against the HIV-1 were observed in patients, elicitation of such bnAbs remains a major challenge when compared to other viral targets. We previously hypothesized that HIV-1 could have evolved a strategy to evade the immune system due to absent or very weak binding of germline antibodies to the conserved epitopes that may not be sufficient to initiate and/or maintain an effective immune response. To further explore our hypothesis, we used the 454 sequence analysis of a large naïve library of human IgM antibodies which had been used for selecting antibodies against SARS CoV receptor-binding domain (RBD), and soluble G proteins (sG) of henipaviruses. We found that the human IgM repertoires from the 454 sequencing have diverse germline usages, recombination patterns, junction diversity, and a lower extent of somatic mutation. In this study, we identified antibody maturation intermediates that are related to bnAbs against the HIV-1 and other viruses as observed in normal individuals, and compared their genetic diversity and somatic mutation level along with available structural and functional data. Further computational analysis will provide framework for understanding the underlying genetic and molecular determinants related to maturation pathways of antiviral bnAbs that could be useful for applying novel approaches to the design of effective vaccine immunogens and antibody-based therapeutics.