Direct high-throughput amplification and sequencing of immunoglobulin genes from single human B cells

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.

Clonal evolution and TCR specificity of the human TFH cell response to Plasmodium falciparum CSP

T follicular helper (T_FH) cells play a crucial role in the development of long-lived, high-quality B cell responses after infection and vaccination. However, little is known about how antigen-specific T_FH cells clonally evolve in response to complex pathogens and what guides the targeting of different epitopes. Here, we assessed the cell phenotype, clonal dynamics, and T cell receptor (TCR) specificity of human circulating T_FH (cT_FH) cells during successive malaria immunizations with radiation-attenuated Plasmodium falciparum (Pf) sporozoites. Repeated parasite exposures induced a dynamic, polyclonal cT_FH response with high frequency of cells specific to a small number of epitopes in Pf circumsporozoite protein (PfCSP), the primary sporozoite surface protein and well-defined vaccine target. Human leukocyte antigen (HLA) restrictions and differences in TCR generation probability were associated with differences in the epitope targeting frequency and indicated the potential of amino acids 311 to 333 in the Th2R/T* region as a T cell supertope. But most of vaccine-induced anti-amino acid 311 to 333 TCRs, including convergent TCRs with high sequence similarity, failed to tolerate natural polymorphisms in their target peptide sequence, thus demonstrating that the T_FH cell response was limited to the vaccine strain. These data suggest that the high parasite diversity in endemic areas will limit boosting of the vaccine-induced T_FH cell response by natural infections. Our findings may guide the further design of PfCSP-based malaria vaccines able to induce potent T helper cell responses for broad, long-lasting antibody responses.

An efficient single-cell based method for linking human T cell phenotype to T cell receptor sequence and specificity

Single-cell antigen-receptor gene amplification and sequencing platforms have been used to characterize T cell receptor (TCR) repertoires but typically fail to generate paired full-length gene products for direct expression cloning and do not enable linking this data to cell phenotype information. To overcome these limitations, we established a high-throughput platform for the quantitative and qualitative analysis of human TCR repertoires that provides insights into the clonal and functional composition of human CD4⁺ and CD8⁺ αβ T cells at the molecular and cellular level. The strategy is a powerful tool to qualitatively assess differences between antigen receptors of phenotypically defined αβ T cell subsets, e.g. in immune responses to cancer, vaccination, or infection, and in autoimmune diseases.

ALDH4A1 is an atherosclerosis auto-antigen targeted by protective antibodies

Cardiovascular disease (CVD) is the leading cause of mortality in the world, with most CVD-related deaths resulting from myocardial infarction or stroke. The main underlying cause of thrombosis and cardiovascular events is atherosclerosis, an inflammatory disease that can remain asymptomatic for long periods. There is an urgent need for therapeutic and diagnostic options in this area. Atherosclerotic plaques contain autoantibodies^1,2, and there is a connection between atherosclerosis and autoimmunity³. However, the immunogenic trigger and the effects of the autoantibody response during atherosclerosis are not well understood^3-5. Here we performed high-throughput single-cell analysis of the atherosclerosis-associated antibody repertoire. Antibody gene sequencing of more than 1,700 B cells from atherogenic Ldlr^-/- and control mice identified 56 antibodies expressed by in-vivo-expanded clones of B lymphocytes in the context of atherosclerosis. One-third of the expanded antibodies were reactive against atherosclerotic plaques, indicating that various antigens in the lesion can trigger antibody responses. Deep proteomics analysis identified ALDH4A1, a mitochondrial dehydrogenase involved in proline metabolism, as a target antigen of one of these autoantibodies, A12. ALDH4A1 distribution is altered during atherosclerosis, and circulating ALDH4A1 is increased in mice and humans with atherosclerosis, supporting the potential use of ALDH4A1 as a disease biomarker. Infusion of A12 antibodies into Ldlr^-/- mice delayed plaque formation and reduced circulating free cholesterol and LDL, suggesting that anti-ALDH4A1 antibodies can protect against atherosclerosis progression and might have therapeutic potential in CVD.

IGLV3-21*01 is an inherited risk factor for CLL through the acquisition of a single-point mutation enabling autonomous BCR signaling

CLL is characterized by autonomous B cell receptor (BCR) signaling. CLL subsets are empirically defined by sequence similarities of the BCR heavy chain. However, in the unfavorable subset 2, an acquired mutation (termed R110) in the light chain stimulates autonomous BCR signaling. This study demonstrates that the oncogenic R110 mutation dictates the unfavorable prognosis and is not restricted to the conventional subset 2. Interestingly, carriers of a particular light-chain allele (IGLV3-21*01) are predisposed to develop CLL because this allele enables autonomous BCR signaling by R110 as a single-point mutation. Monoclonal antibodies permit convenient screening for R110-expressing CLL, showing that it is the largest immunologically defined CLL subset and an example of functional rather than empirical CLL subclassification.

The prognosis of chronic lymphocytic leukemia (CLL) depends on different markers, including cytogenetic aberrations, oncogenic mutations, and mutational status of the immunoglobulin (Ig) heavy-chain variable (IGHV) gene. The number of IGHV mutations distinguishes mutated (M) CLL with a markedly superior prognosis from unmutated (UM) CLL cases. In addition, B cell antigen receptor (BCR) stereotypes as defined by IGHV usage and complementarity-determining regions (CDRs) classify ∼30% of CLL cases into prognostically important subsets. Subset 2 expresses a BCR with the combination of IGHV3-21–derived heavy chains (HCs) with IGLV3-21–derived light chains (LCs), and is associated with an unfavorable prognosis. Importantly, the subset 2 LC carries a single-point mutation, termed R110, at the junction between the variable and constant LC regions. By analyzing 4 independent clinical cohorts through BCR sequencing and by immunophenotyping with antibodies specifically recognizing wild-type IGLV3-21 and R110-mutated IGLV3-21 (IGLV3-21^R110), we show that IGLV3-21^R110–expressing CLL represents a distinct subset with poor prognosis independent of IGHV mutations. Compared with other alleles, only IGLV3-21*01 facilitates effective homotypic BCR–BCR interaction that results in autonomous, oncogenic BCR signaling after acquiring R110 as a single-point mutation. Presumably, this mutation acts as a standalone driver that transforms IGLV3-21*01–expressing B cells to develop CLL. Thus, we propose to expand the conventional definition of CLL subset 2 to subset 2L by including all IGLV3-21^R110–expressing CLL cases regardless of IGHV mutational status. Moreover, the generation of monoclonal antibodies recognizing IGLV3-21 or mutated IGLV3-21^R110 facilitates the recognition of B cells carrying this mutation in CLL patients or healthy donors.

openPrimeR for multiplex amplification of highly diverse templates

To study the diversity of immune receptors and pathogens, multiplex PCR has become a central approach in research and diagnostics. However, insufficient primer design against highly diverse templates often prevents amplification and therefore limits the correct understanding of biological processes. Here, we present openPrimeR, an R-based tool for evaluating and designing multiplex PCR primers. openPrimeR provides a functional and intuitive interface and uses either a greedy algorithm or an integer linear program to compute the minimal set of primers that performs full target coverage. As proof of concept, we used openPrimeR to find optimal primer sets for the amplification of highly mutated immunoglobulins. Comprehensive analyses on specifically generated immunoglobulin variable gene segment libraries resulted in the composition of highly effective primer sets (oPR-IGHV, oPR-IGKV and oPR-IGLV) that demonstrated to be particularly suitable for the isolation of novel human antibodies.

High-throughput single-cell sequencing of paired TCRα and TCRβ genes for the direct expression-cloning and functional analysis of murine T-cell receptors

Precise clonal and functional assessments of the T cell receptor (TCR) repertoire diversity require paired TCRα and TCRβ gene sequence information at monoclonal level. However, available single‐cell strategies are typically limited in throughput and often do not provide full‐length DNA templates for direct gene cloning. Here, we describe a high‐throughput strategy for the unbiased amplification and automated sequence analysis of paired TCRα and TCRβ genes from primary mouse T cells. The platform links cell phenotype and TCR gene sequence information at single‐cell level. Furthermore, it enables direct functional analyses through the efficient cloning of both genes and the generation of stable TCR expressing cell lines. This highly efficient workflow is a powerful tool to determine the diversity and quality of the murine T‐cell repertoire in various settings, for example in vaccine development, infectious diseases, autoimmunity, or cancer.

Expression Cloning of Antibodies from Single Human B Cells

The majority of lymphomas originate from B cells at the germinal center stage. Preferential selection of B-cell clones by a limited set of antigens has been suggested to drive lymphoma development. While recent studies in chronic lymphocytic leukemia have shown that self-reactive B-cell receptors (BCR) can generate cell-autonomous signaling and proliferation, our knowledge about the role of BCRs for the development or survival of other lymphomas remains limited. Here, we describe a strategy to characterize the antibody reactivity of human B cells. The approach allows the unbiased characterization of the human antibody repertoire at single-cell level through the generation of recombinant monoclonal antibodies from single primary human B cells of defined origin. This protocol offers a detailed description of the method starting from the flow-cytometric isolation of single human B cells to the reverse transcription-polymerase chain reaction (RT-PCR)-based amplification of the expressed immunoglobulin (Ig) transcripts (IGH, IGK, and IGL) and their subsequent cloning into expression vectors for the in vitro production of recombinant monoclonal antibodies. The strategy may be used to obtain information on the clonal evolution of B-cell lymphomas by single-cell sequencing of Ig transcripts and on the antibody reactivity of human lymphoma B cells.

Sequencing and Affinity Determination of Antigen-Specific B Lymphocytes from Peripheral Blood

Here we describe methods for screening human blood to isolate peripheral blood mononuclear cells (PBMCs) capable of binding fluorescently labeled antigen, as well as methods for the amplification and sequencing of B cell receptor (BCR) heavy and light chain genes. Detailed protocols are provided for transient mammalian expression in a hexahistidine-tagged Fab format, purification by immobilized metal affinity chromatography (IMAC), and affinity determination by BioLayer interferometry (BLI).

Using High-Throughput Sequencing to Characterize the Development of the Antibody Repertoire During Infections: A Case Study of HIV-1

High throughput sequencing (HTS) approaches have only recently been applied to describing the antibody/B-cell repertoire in fine detail, but these data sets have already become critical to the design of vaccines and therapeutics, and monitoring of cancer immunotherapy. As a case study, we describe the potential and present limitations of HTS studies of the Ab repertoire during infection with HIV-1. Most of the present studies restrict their analyses to lineages of specific bnAbs. We discuss future initiatives to expand this type of analysis to more complete repertoires and to improve comparing and sharing of these Ab repertoire data across studies and institutions.

Computational Evaluation of B-Cell Clone Sizes in Bulk Populations

B cell clones expand and contract during adaptive immune responses and can persist or grow uncontrollably in lymphoproliferative disorders. One way to monitor and track B cell clones is to perform large-scale sampling of bulk cell populations, amplifying, and sequencing antibody gene rearrangements by next-generation sequencing (NGS). Here, we describe a series of computational approaches for estimating B cell clone size in NGS immune repertoire profiling data of antibody heavy chain gene rearrangements. We define three different measures of B cell clone size-copy numbers, instances, and unique sequences-and show how these measures can be used to rank clones, analyze their diversity, and study their distribution within and between individuals. We provide a detailed, step-by-step procedure for performing these analyses using two different data sets of spleen samples from human organ donors. In the first data set, 19 independently generated biological replicates from a single individual are analyzed for B cell clone size, diversity and sampling sufficiency for clonal overlap analysis. In the second data set, B cell clones are compared in eight different organ donors. We comment upon frequently encountered pitfalls and offer practical advice with alternative approaches. Overall, we provide a series of pragmatic analytical approaches and show how different clone size measures can be used to study the clonal landscape in bulk B cell immune repertoire profiling data.

Cross-specificity of protective human antibodies against Klebsiella pneumoniae LPS O-antigen

Humoral immune responses to microbial polysaccharide surface antigens can prevent bacterial infection but are typically strain specific and fail to mediate broad protection against different serotypes. Here we describe a panel of affinity-matured monoclonal human antibodies from peripheral blood immunoglobulin M-positive (IgM⁺) and IgA⁺ memory B cells and clonally related intestinal plasmablasts, directed against the lipopolysaccharide (LPS) O-antigen of Klebsiella pneumoniae, an opportunistic pathogen and major cause of antibiotic-resistant nosocomial infections. The antibodies showed distinct patterns of in vivo cross-specificity and protection against different clinically relevant K. pneumoniae serotypes. However, cross-specificity was not limited to K. pneumoniae, as K. pneumoniae-specific antibodies recognized diverse intestinal microbes and neutralized not only K. pneumoniae LPS but also non-K. pneumoniae LPS. Our data suggest that the recognition of minimal glycan epitopes abundantly expressed on microbial surfaces might serve as an efficient humoral immunological mechanism to control invading pathogens and the large diversity of the human microbiota with a limited set of cross-specific antibodies.

Rare PfCSP C-terminal antibodies induced by live sporozoite vaccination are ineffective against malaria infection

Scally et al. show the molecular, structural, and functional characterization of human antibodies against the C-terminal domain of Plasmodium falciparum (Pf) circumsporozoite (CSP [C-PfCSP]) and reveal that its arrangement on the Pf sporozoite surface and epitope polymorphism contribute to poor C-PfCSP immunogenicity and ineffective humoral responses in volunteers protected against Pf malaria.

Antibodies against the central repeat of the Plasmodium falciparum (Pf) circumsporozoite protein (CSP) inhibit parasite activity and correlate with protection from malaria. However, the humoral response to the PfCSP C terminus (C-PfCSP) is less well characterized. Here, we describe B cell responses to C-PfCSP from European donors who underwent immunization with live Pf sporozoites (PfSPZ Challenge) under chloroquine prophylaxis (PfSPZ-CVac), and were protected against controlled human malaria infection. Out of 215 PfCSP-reactive monoclonal antibodies, only two unique antibodies were specific for C-PfCSP, highlighting the rare occurrence of C-PfCSP–reactive B cells in PfSPZ-CVac–induced protective immunity. These two antibodies showed poor sporozoite binding and weak inhibition of parasite traversal and development, and did not protect mice from infection with PfCSP transgenic Plasmodium berghei sporozoites. Structural analyses demonstrated that one antibody interacts with a polymorphic region overlapping two T cell epitopes, suggesting that variability in C-PfCSP may benefit parasite escape from humoral and cellular immunity. Our data identify important features underlying C-PfCSP shortcomings as a vaccine target.

Integrating high-throughput screening and sequencing for monoclonal antibody discovery and engineering

Monoclonal antibody discovery and engineering is a field that has traditionally been dominated by high-throughput screening platforms (e.g. hybridomas and surface display). In recent years the emergence of high-throughput sequencing has made it possible to obtain large-scale information on antibody repertoire diversity. Additionally, it has now become more routine to perform high-throughput sequencing on antibody repertoires to also directly discover antibodies. In this review, we provide an overview of the progress in this field to date and show how high-throughput screening and sequencing are converging to deliver powerful new workflows for monoclonal antibody discovery and engineering.

Next-generation sequencing in chronic lymphocytic leukemia: recent findings and new horizons

The rapid progress in next-generation sequencing technologies has significantly contributed to our knowledge of the genetic events associated with the development, progression and treatment resistance of chronic lymphocytic leukemia patients. Together with the discovery of new driver mutations, next-generation sequencing has revealed an immense degree of both intra- and inter-tumor heterogeneity and enabled us to describe marked clonal evolution. Advances in immunogenetics may be implemented to detect minimal residual disease more sensitively and to track clonal B cell populations, their dynamics and molecular characteristics. The interpretation of these aspects is indispensable to thoroughly examine the genetic background of chronic lymphocytic leukemia. We review and discuss the recent results provided by the different next-generation sequencing techniques used in studying the chronic lymphocytic leukemia genome, as well as future perspectives in the methodologies and applications.

Novel Approaches to Analyze Immunoglobulin Repertoires

Analysis of immunoglobulin (Ig) repertoires aims to comprehend Ig diversity with the goal of predicting humoral immune responses in the context of infection, vaccination, autoimmunity, and malignancies. The first next-generation sequencing (NGS) analyses of bulk B cell populations dramatically advanced sampling depth over previous low-throughput single-cell-based protocols, albeit at the expense of accuracy and loss of chain-pairing information. In recent years the field has substantially differentiated, with bulk analyses becoming more accurate while single-cell approaches have gained in throughput. Additionally, new platforms striving to combine high throughput and chain pairing have been developed as well as various computational tools for analysis. Here we review the developments of the past 4-5 years and discuss the open challenges.

Exploring Human Antimicrobial Antibody Responses on a Single B Cell Level

Analysis of monoclonal antibodies (MAbs) derived from single B cell cloning has been highly beneficial for antimicrobial immunotherapy, vaccine design, and advancing our understanding of pathogen-triggered effects on the human immunoglobulin repertoire. Sequencing of variable domains of single B cells, and characterization of binding and functional activities of MAbs derived from those sequences, provides in-depth insight not only into sites of susceptibility for antibody-mediated neutralization or opsonization of the pathogen but also into the dynamics of protective antibody evolution during infection. This information can be utilized to rapidly develop novel immunotherapies of completely human origin and provides a roadmap for structure-based vaccine design that aims to elicit similar protective antibody responses. Here, we summarize recent aspects of the single B cell cloning approach.

Advanced Methodologies in High-Throughput Sequencing of Immune Repertoires

In recent years, major efforts have been made to develop sophisticated experimental and bioinformatic workflows for sequencing adaptive immune repertoires. The immunological insight gained has been applied to fields as varied as lymphocyte biology, immunodiagnostics, vaccines, cancer immunotherapy, and antibody engineering. In this review, we provide a detailed overview of these advanced methodologies, focusing specifically on strategies to reduce sequencing errors and bias and to achieve high-throughput pairing of variable regions (e.g., heavy-light or alpha-beta chains). In addition, we highlight recent technologies for single-cell transcriptome sequencing that can be integrated with immune repertoires. Finally, we provide a perspective on advanced immune repertoire sequencing and its ability to impact basic immunology, biopharmaceutical drug discovery and development, and cancer immunotherapy.

Natively glycosylated HIV-1 Env structure reveals new mode for antibody recognition of the CD4-binding site

HIV-1 vaccine design is informed by structural studies elucidating mechanisms by which broadly neutralizing antibodies (bNAbs) recognize and/or accommodate N-glycans on the trimeric envelope glycoprotein (Env). Variability in high-mannose and complex-type Env glycoforms leads to heterogeneity that usually precludes visualization of the native glycan shield. We present 3.5-Å- and 3.9-Å-resolution crystal structures of the HIV-1 Env trimer with fully processed and native glycosylation, revealing a glycan shield of high-mannose and complex-type N-glycans, which we used to define complete epitopes of two bNAbs. Env trimer was complexed with 10-1074 (against the V3-loop) and IOMA, a new CD4-binding site (CD4bs) antibody. Although IOMA derives from VH1-2*02, the germline gene of CD4bs-targeting VRC01-class bNAbs, its light chain lacks the short CDRL3 that defines VRC01-class bNAbs. Thus IOMA resembles 8ANC131-class/VH1-46-derived CD4bs bNAbs, which have normal-length CDRL3s. The existence of bNAbs that combine features of VRC01-class and 8ANC131-class antibodies has implications for immunization strategies targeting VRC01-like bNAbs.

sciReptor: analysis of single-cell level immunoglobulin repertoires

Background

The sequencing of immunoglobulin (Ig) transcripts from single B cells yields essential information about Ig heavy:light chain pairing, which is lost in conventional bulk sequencing experiments. The previously limited throughput of single-cell approaches has recently been overcome by the introduction of multiple next-generation sequencing (NGS)-based platforms. Furthermore, single-cell techniques allow the assignment of additional data types (e.g. cell surface marker expression), which are crucial for biological interpretation. However, the currently available computational tools are not designed to handle single-cell data and do not provide integral solutions for linking of sequence data to other biological data.

Results

Here we introduce sciReptor, a flexible toolkit for the processing and analysis of antigen receptor repertoire sequencing data at single-cell level. The software combines bioinformatics tools for immunoglobulin sequence annotation with a relational database, where raw data and analysis results are stored and linked. sciReptor supports attribution of additional data categories such as cell surface marker expression or immunological metadata. Furthermore, it comprises a quality control module as well as basic repertoire visualization tools.

Conclusion

sciReptor is a flexible framework for standardized sequence analysis of antigen receptor repertoires on single-cell level. The relational database allows easy data sharing and downstream analyses as well as immediate comparisons between different data sets.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-0920-1) contains supplementary material, which is available to authorized users.

Isolation of human monoclonal autoantibodies derived from pancreatic lymph node and peripheral blood B cells of islet autoantibody-positive patients

AIMS/HYPOTHESIS

Autoantibodies against pancreatic islets and infections by enteroviruses are associated with type 1 diabetes, but the specificity of immune responses within the type 1 diabetic pancreas is poorly characterised. We investigated whether pancreatic lymph nodes could provide a source of antigen-specific B cells for analysis of immune responses within the (pre)diabetic pancreas.

METHODS

Human IgG antibodies were cloned from single B lymphocytes sorted from pancreatic lymph node cells of three organ donors positive for islet autoantibodies, and from the peripheral blood of a patient with type 1 diabetes. Antibodies to insulinoma-associated antigen 2 (IA-2), GAD65, zinc transporter 8 (ZnT8) and Coxsackie B virus proteins were assayed by immunoprecipitation and by immunofluorescence on pancreatic sections.

RESULTS

Human IgG antibodies (863) were successfully cloned and produced from 4,092 single B cells from lymph nodes and peripheral blood. Reactivity to the protein tyrosine phosphatase domain of the IA-2 autoantigen was detected in two cloned antibodies: one derived from a pancreatic lymph node and one from peripheral blood. Epitopes for these two antibodies were similar to each other and to those for circulating antibodies in type 1 diabetes. The remaining 861 antibodies were negative for reactivity to IA-2, GAD65 or ZnT8 by both assays tested. Reactivity to a Coxsackie viral protein 2 was detected in one antibody derived from a peripheral blood B cell, but not from lymph nodes.

CONCLUSIONS/INTERPRETATION

We show evidence for the infrequent presence of autoantigen-specific IgG+ B lymphocytes in the pancreatic-draining lymph nodes of islet autoantibody-positive individuals.