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

Computational detection of antigen-specific B cell receptors following immunization

B cell receptors (BCRs) play a crucial role in recognizing and fighting foreign antigens. High-throughput sequencing enables in-depth sampling of the BCRs repertoire after immunization. However, only a minor fraction of BCRs actively participate in any given infection. To what extent can we accurately identify antigen-specific sequences directly from BCRs repertoires? We present a computational method grounded on sequence similarity, aimed at identifying statistically significant responsive BCRs. This method leverages well-known characteristics of affinity maturation and expected diversity. We validate its effectiveness using longitudinally sampled human immune repertoire data following influenza vaccination and SARS-CoV-2 infections. We show that different lineages converge to the same responding Complementarity Determining Region 3, demonstrating convergent selection within an individual. The outcomes of this method hold promise for application in vaccine development, personalized medicine, and antibody-derived therapeutics.

Biomarkers of immunotherapy for non-small cell lung cancer

Immunotherapy is revolutionizing the treatment of non-small cell lung cancer by targeting immune checkpoint proteins, including programmed death-1, programmed death ligand 1 and cytotoxic T-lymphocyte-associated antigen 4. Several immune checkpoint inhibitors, including programmed death ligand 1 inhibitors, programmed death-1 inhibitors and cytotoxic T-lymphocyte-associated antigen 4 inhibitors, were approved for the treatment of patients with advanced non-small cell lung cancer. Programmed death ligand 1 expression is currently the only predictive biomarker for immune checkpoint inhibitors to guide the treatment strategy in these patients. However, programmed death ligand 1 expression is not a perfect biomarker for predicting the efficacy of immunotherapy. Therefore, various biomarkers such as tumour mutation burden, tumour microenvironment, gut microbiome and T-cell receptor repertoire have been proposed to predict the efficacy of immunotherapy more accurately. Additionally, combining different biomarkers may provide a more accurate prediction of response to immunotherapy. This article reports the review of the latest evidence of the predictive marker of immunotherapy in patients with advanced non-small cell lung cancer.

An open protocol for modeling T Cell Clonotype repertoires using TCRβ CDR3 sequences

T cell receptor repertoires can be profiled using next generation sequencing (NGS) to measure and monitor adaptive dynamical changes in response to disease and other perturbations. Genomic DNA-based bulk sequencing is cost-effective but necessitates multiplex target amplification using multiple primer pairs with highly variable amplification efficiencies. Here, we utilize an equimolar primer mixture and propose a single statistical normalization step that efficiently corrects for amplification bias post sequencing. Using samples analyzed by both our open protocol and a commercial solution, we show high concordance between bulk clonality metrics. This approach is an inexpensive and open-source alternative to commercial solutions.

An open protocol for modeling T Cell Clonotype repertoires using TCRβ CDR3 sequences

T cell receptor repertoires can be profiled using next generation sequencing (NGS) to measure and monitor adaptive dynamical changes in response to disease and other perturbations. Genomic DNA-based bulk sequencing is cost-effective but necessitates multiplex target amplification using multiple primer pairs with highly variable amplification efficiencies. Here, we utilize an equimolar primer mixture and propose a single statistical normalization step that efficiently corrects for amplification bias post sequencing. Using samples analyzed by both our open protocol and a commercial solution, we show high concordance between bulk clonality metrics. This approach is an inexpensive and open-source alternative to commercial solutions.

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.

Circulating T-cell receptor diversity as predictive biomarker for PARP inhibitors maintenance therapy in high grade serous ovarian cancer

OBJECTIVE

T-cell receptor (TCR) repertoire diversity is getting increasing attention as a predictive biomarker in cancer patients. However, the characteristics of the TCR together with its predictive significance for high grade serous ovarian cancer (HGSOC) patients receiving poly (ADP-ribose) polymerase inhibitor (PARPi) maintenance therapy remain unknown.

METHODS

Twenty-seven patients with HGSOC were analyzed including 22 patients receiving PARPi maintenance therapy and 5 untreated patients as control. Peripheral blood samples were collected for TCR sequencing at baseline as well as one month and three months after the exposure to PARPi. To determine whether TCR diversity was related to PARPi efficacy, we compared the TCR repertoire between patients who had received PARPi and those who had not.

RESULTS

For patients receiving PARPi treatment or not, we evaluated changes in clone abundance during PARPi maintenance and the similarity of the TCR repertoire before and after the treatment. The results revealed that patients receiving PARPi had TCR repertoires that were more stable than those of untreated cases. We next correlated TCR diversity with the efficacy of PARPi in the treatment group. The rising trend of TCR diversity after three months with PARPi treatment was associated with a longer PFS (21.7 vs 7.4 months, hazard ratio = 0.19, p < 0.001) and a better response to PARPi (91.7% vs 25.0%, p = 0.004). Furthermore, we discovered that the primary characteristic with predictive value for the effectiveness of PARPi is the considerable reduction of the high-frequency T cell clones.

CONCLUSION

We suggested that the circulating TCR diversity could be a potential predictive biomarker for PARPi maintenance therapy in HGSOC.

Immune environment and antigen specificity of the T cell receptor repertoire of malignant ascites in ovarian cancer

We evaluated the association of disease outcome with T cell immune-related characteristics and T cell receptor (TCR) repertoire in malignant ascites from patients with high-grade epithelial ovarian cancer. Ascitic fluid samples were collected from 47 high-grade epithelial ovarian cancer patients and analyzed using flow cytometry and TCR sequencing to characterize the complementarity determining region 3 TCR β-chain. TCR functions were analyzed using the McPAS-TCR and VDJ databases. TCR clustering was implemented using Grouping of Lymphocyte Interactions by Paratope Hotspots software. Patients with poor prognosis had ascites characterized by an increased ratio of CD8+ T cells to regulatory T cells, which correlated with an increased productive frequency of the top 100 clones and decreased productive entropy. TCRs enriched in patients with an excellent or good prognosis were more likely to recognize cancer antigens and contained more TCR reads predicted to recognize epithelial ovarian cancer antigens. In addition, a TCR motif that is predicted to bind the TP53 neoantigen was identified, and this motif was enriched in patients with an excellent or good prognosis. Ascitic fluid in high-grade epithelial ovarian cancer patients with an excellent or good prognosis is enriched with TCRs that may recognize ovarian cancer-specific neoantigens, including mutated TP53 and TEAD1. These results suggest that an effective antigen-specific immune response in ascites is vital for a good outcome in high-grade epithelial ovarian cancer.

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.

High-throughput retrieval of physical DNA for NGS-identifiable clones in phage display library

In antibody discovery, in-depth analysis of an antibody library and high-throughput retrieval of clones in the library are crucial to identifying and exploiting rare clones with different properties. However, existing methods have technical limitations, such as low process throughput from the laborious cloning process and waste of the phenotypic screening capacity from unnecessary repetitive tests on the dominant clones. To overcome the limitations, we developed a new high-throughput platform for the identification and retrieval of clones in the library, TrueRepertoire™. This new platform provides highly accurate sequences of the clones with linkage information between heavy and light chains of the antibody fragment. Additionally, the physical DNA of clones can be retrieved in high throughput based on the sequence information. We validated the high accuracy of the sequences and demonstrated that there is no platform-specific bias. Moreover, the applicability of TrueRepertoire™ was demonstrated by a phage-displayed single-chain variable fragment library targeting human hepatocyte growth factor protein.

A comparison of unamplified and massively multiplexed PCR amplification for murine antibody repertoire sequencing

Sequencing antibody repertoires has steadily become cheaper and easier. Sequencing methods usually rely on some form of amplification, often a massively multiplexed PCR prior to sequencing. To eliminate potential biases and create a data set that could be used for other studies, our laboratory compared unamplified sequencing results from the splenic heavy-chain repertoire in the mouse to those processed through two commercial applications. We also compared the use of mRNA vs total RNA, reverse transcriptase, and primer usage for cDNA synthesis and submission. The use of mRNA for cDNA synthesis resulted in higher read counts but reverse transcriptase and primer usage had no statistical effects on read count. Although most of the amplified data sets contained more antibody reads than the unamplified data set, we detected more unique variable (V)-gene segments in the unamplified data set. Although unique CDR3 detection was much lower in the unamplified data set, RNASeq detected 98% of the high-frequency CDR3s. We have shown that unamplified profiling of the antibody repertoire is possible, detects more V-gene segments, and detects high-frequency clones in the repertoire.

High-throughput sequencing of murine immunoglobulin heavy chain repertoires using single side unique molecular identifiers on an Ion Torrent PGM

With the advent of high-throughput sequencing (HTS), profiling immunoglobulin (IG) repertoires has become an essential part of immunological research. Advances in sequencing technology enable the IonTorrent Personal Genome Machine (PGM) to cover the full-length of IG mRNA transcripts. Nucleotide insertions and deletions (indels) are the dominant errors of the PGM sequencing platform and can critically influence IG repertoire assessments. Here, we present a PGM-tailored IG repertoire sequencing approach combining error correction through unique molecular identifier (UID) barcoding and indel detection through ImMunoGeneTics (IMGT), the most commonly used sequence alignment database for IG sequences. Using artificially falsified sequences for benchmarking, we found that IMGT's underlying algorithms efficiently detect 98% of the introduced indels. Undetected indels are either located at the end of the sequences or produce masked frameshifts with an insertion and deletion in close proximity. The complementary determining regions 3 (CDR3s) are returned correct for up to 3 insertions or 3 deletions through conservative culling. We further show, that our PGM-tailored unique molecular identifiers result in highly accurate HTS data if combined with the presented processing strategy. In this regard, considering sequences with at least two copies from datasets with UID families of minimum 3 reads result in correct sequences with over 99% confidence. Finally, we show that the protocol can readily be used to generate homogenous datasets for bulk sequencing of murine bone marrow samples. Taken together, this approach will help to establish benchtop-scale sequencing of IG heavy chain transcripts in the field of IG repertoire research.

IRProfiler - a software toolbox for high throughput immune receptor profiling

Background

The study of the huge diversity of immune receptors, often referred to as immune repertoire profiling, is a prerequisite for diagnosis, prognostication and monitoring of hematological disorders. In the era of high-throughput sequencing (HTS), the abundance of immunogenetic data has revealed unprecedented opportunities for the thorough profiling of T-cell receptors (TR) and B-cell receptors (BcR). However, the volume of the data to be analyzed mandates for efficient and ease-to-use immune repertoire profiling software applications.

Results

This work introduces Immune Repertoire Profiler (IRProfiler), a novel software pipeline that delivers a number of core receptor repertoire quantification and comparison functionalities on high-throughput TR and BcR sequencing data. Adopting 5 alternative clonotype definitions, IRProfiler implements a series of algorithms for 1) data filtering, 2) calculation of clonotype diversity and expression, 3) calculation of gene usage for the V and J subgroups, 4) detection of shared and exclusive clonotypes among multiple repertoires, and 5) comparison of gene usage for V and J subgroups among multiple repertoires. IRProfiler has been implemented as a toolbox of the Galaxy bioinformatics platform, comprising 6 tools. Theoretical and experimental evaluation has shown that the tools of IRProfiler are able to scale well with respect to the size of input dataset(s). IRProfiler has been utilized by a number of recently published studies concerning hematological disorders.

Conclusion

IRProfiler is made freely available via 3 distribution channels, including the Galaxy Tool Shed. Despite being a new entry in a crowded ecosystem of immune repertoire profiling software, IRProfiler founds its added value on its support for alternative clonotype definitions in conjunction with a combination of properties stemming from its user-centric design, namely ease-of-use, ease-of-access, exploitability of the output data, and analysis flexibility.

Computational Strategies for Dissecting the High-Dimensional Complexity of Adaptive Immune Repertoires

The adaptive immune system recognizes antigens via an immense array of antigen-binding antibodies and T-cell receptors, the immune repertoire. The interrogation of immune repertoires is of high relevance for understanding the adaptive immune response in disease and infection (e.g., autoimmunity, cancer, HIV). Adaptive immune receptor repertoire sequencing (AIRR-seq) has driven the quantitative and molecular-level profiling of immune repertoires, thereby revealing the high-dimensional complexity of the immune receptor sequence landscape. Several methods for the computational and statistical analysis of large-scale AIRR-seq data have been developed to resolve immune repertoire complexity and to understand the dynamics of adaptive immunity. Here, we review the current research on (i) diversity, (ii) clustering and network, (iii) phylogenetic, and (iv) machine learning methods applied to dissect, quantify, and compare the architecture, evolution, and specificity of immune repertoires. We summarize outstanding questions in computational immunology and propose future directions for systems immunology toward coupling AIRR-seq with the computational discovery of immunotherapeutics, vaccines, and immunodiagnostics.

Functionally Convergent B Cell Receptor Sequences in Transgenic Rats Expressing a Human B Cell Repertoire in Response to Tetanus Toxoid and Measles Antigens

The identification and tracking of antigen-specific immunoglobulin (Ig) sequences within total Ig repertoires is central to high-throughput sequencing (HTS) studies of infections or vaccinations. In this context, public Ig sequences shared by different individuals exposed to the same antigen could be valuable markers for tracing back infections, measuring vaccine immunogenicity, and perhaps ultimately allow the reconstruction of the immunological history of an individual. Here, we immunized groups of transgenic rats expressing human Ig against tetanus toxoid (TT), Modified Vaccinia virus Ankara (MVA), measles virus hemagglutinin and fusion proteins expressed on MVA, and the environmental carcinogen benzo[a]pyrene, coupled to TT. We showed that these antigens impose a selective pressure causing the Ig heavy chain (IgH) repertoires of the rats to converge toward the expression of antibodies with highly similar IgH CDR3 amino acid sequences. We present a computational approach, similar to differential gene expression analysis, that selects for clusters of CDR3s with 80% similarity, significantly overrepresented within the different groups of immunized rats. These IgH clusters represent antigen-induced IgH signatures exhibiting stereotypic amino acid patterns including previously described TT- and measles-specific IgH sequences. Our data suggest that with the presented methodology, transgenic Ig rats can be utilized as a model to identify antigen-induced, human IgH signatures to a variety of different antigens.

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.

Overview of methodologies for T-cell receptor repertoire analysis

Background

The T-cell receptor (TCR), located on the surface of T cells, is responsible for the recognition of the antigen-major histocompatibility complex, leading to the initiation of an inflammatory response. Analysing the TCR repertoire may help to gain a better understanding of the immune system features and of the aetiology and progression of diseases, in particular those with unknown antigenic triggers. The extreme diversity of the TCR repertoire represents a major analytical challenge; this has led to the development of specialized methods which aim to characterize the TCR repertoire in-depth. Currently, next generation sequencing based technologies are most widely employed for the high-throughput analysis of the immune cell repertoire.

Results

Here, we report on the latest methodological advancements in the field by describing and comparing the available tools; from the choice of the starting material and library preparation method, to the sequencing technologies and data analysis.

Finally, we provide a practical example and our own experience by reporting some exemplary results from a small internal benchmark study, where current approaches from the literature and the market are employed and compared.

Conclusions

Several valid methods for clonotype identification and TCR repertoire analysis exist, however, a gold standard method for the field has not yet been identified. Depending on the purpose of the scientific study, some approaches may be more suitable than others. Finally, due to possible method specific biases, scientists must be careful when comparing results obtained using different methods.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-017-0379-9) contains supplementary material, which is available to authorized users.

Insights into immune system development and function from mouse T-cell repertoires

The ability of the adaptive immune system to respond to arbitrary pathogens stems from the broad diversity of immune cell surface receptors. This diversity originates in a stochastic DNA editing process (VDJ recombination) that acts on the surface receptor gene each time a new immune cell is created from a stem cell. By analyzing T-cell receptor (TCR) sequence repertoires taken from the blood and thymus of mice of different ages, we quantify the changes in the VDJ recombination process that occur from embryo to young adult. We find a rapid increase with age in the number of random insertions and a dramatic increase in diversity. Because the blood accumulates thymic output over time, blood repertoires are mixtures of different statistical recombination processes, and we unravel the mixture statistics to obtain a picture of the time evolution of the early immune system. Sequence repertoire analysis also allows us to detect the statistical impact of selection on the output of the VDJ recombination process. The effects we find are nearly identical between thymus and blood, suggesting that our analysis mainly detects selection for proper folding of the TCR receptor protein. We further find that selection is weaker in laboratory mice than in humans and it does not affect the diversity of the repertoire.

RTCR: a pipeline for complete and accurate recovery of T cell repertoires from high throughput sequencing data

Motivation: High Throughput Sequencing (HTS) has enabled researchers to probe the human T cell receptor (TCR) repertoire, which consists of many rare sequences. Distinguishing between true but rare TCR sequences and variants generated by polymerase chain reaction (PCR) and sequencing errors remains a formidable challenge. The conventional approach to handle errors is to remove low quality reads, and/or rare TCR sequences. Such filtering discards a large number of true and often rare TCR sequences. However, accurate identification and quantification of rare TCR sequences is essential for repertoire diversity estimation.

Results: We devised a pipeline, called Recover TCR (RTCR), that accurately recovers TCR sequences, including rare TCR sequences, from HTS data (including barcoded data) even at low coverage. RTCR employs a data-driven statistical model to rectify PCR and sequencing errors in an adaptive manner. Using simulations, we demonstrate that RTCR can easily adapt to the error profiles of different types of sequencers and exhibits consistently high recall and high precision even at low coverages where other pipelines perform poorly. Using published real data, we show that RTCR accurately resolves sequencing errors and outperforms all other pipelines.

Availability and Implementation: The RTCR pipeline is implemented in Python (v2.7) and C and is freely available at http://uubram.github.io/RTCR/along with documentation and examples of typical usage.

Contact:b.gerritsen@uu.nl

Dynamics of Individual T Cell Repertoires: From Cord Blood to Centenarians

The diversity, architecture, and dynamics of the TCR repertoire largely determine our ability to effectively withstand infections and malignancies with minimal mistargeting of immune responses. In this study, we have employed deep TCRβ repertoire sequencing with normalization based on unique molecular identifiers to explore the long-term dynamics of T cell immunity. We demonstrate remarkable stability of repertoire, where approximately half of all T cells in peripheral blood are represented by clones that persist and generally preserve their frequencies for 3 y. We further characterize the extremes of lifelong TCR repertoire evolution, analyzing samples ranging from umbilical cord blood to centenarian peripheral blood. We show that the fetal TCR repertoire, albeit structurally maintained within regulated borders due to the lower numbers of randomly added nucleotides, is not limited with respect to observed functional diversity. We reveal decreased efficiency of nonsense-mediated mRNA decay in umbilical cord blood, which may reflect specific regulatory mechanisms in development. Furthermore, we demonstrate that human TCR repertoires are functionally more similar at birth but diverge during life, and we track the lifelong behavior of CMV- and EBV-specific T cell clonotypes. Finally, we reveal gender differences in dynamics of TCR diversity constriction, which come to naught in the oldest age. Based on our data, we propose a more general explanation for the previous observations on the relationships between longevity and immunity.

Accurate and predictive antibody repertoire profiling by molecular amplification fingerprinting

High-throughput antibody repertoire sequencing (Ig-seq) provides quantitative molecular information on humoral immunity. However, Ig-seq is compromised by biases and errors introduced during library preparation and sequencing. By using synthetic antibody spike-in genes, we determined that primer bias from multiplex polymerase chain reaction (PCR) library preparation resulted in antibody frequencies with only 42 to 62% accuracy. Additionally, Ig-seq errors resulted in antibody diversity measurements being overestimated by up to 5000-fold. To rectify this, we developed molecular amplification fingerprinting (MAF), which uses unique molecular identifier (UID) tagging before and during multiplex PCR amplification, which enabled tagging of transcripts while accounting for PCR efficiency. Combined with a bioinformatic pipeline, MAF bias correction led to measurements of antibody frequencies with up to 99% accuracy. We also used MAF to correct PCR and sequencing errors, resulting in enhanced accuracy of full-length antibody diversity measurements, achieving 98 to 100% error correction. Using murine MAF-corrected data, we established a quantitative metric of recent clonal expansion-the intraclonal diversity index-which measures the number of unique transcripts associated with an antibody clone. We used this intraclonal diversity index along with antibody frequencies and somatic hypermutation to build a logistic regression model for prediction of the immunological status of clones. The model was able to predict clonal status with high confidence but only when using MAF error and bias corrected Ig-seq data. Improved accuracy by MAF provides the potential to greatly advance Ig-seq and its utility in immunology and biotechnology.

Next-Generation Sequencing of a Single Domain Antibody Repertoire Reveals Quality of Phage Display Selected Candidates

Next-Generation Sequencing and bioinformatics are powerful tools for analyzing the large number of DNA sequences present in an immune library. In this work, we constructed a cDNA library of single domain antibodies from a llama immunized with staphylococcal enterotoxin B. The resulting library was sequenced, resulting in approximately 8.5 million sequences with 5.4 million representing intact, useful sequences. The sequenced library was interrogated using sequences of known SEB-binding single domain antibodies from the library obtained through phage display panning methods in a previous study. New antibodies were identified, produced, and characterized, and were shown to have affinities and melting temperatures comparable to those obtained by traditional panning methods. This demonstrates the utility of using NGS as a complementary tool to phage-displayed biopanning as a means for rapidly obtaining additional antibodies from an immune library. It also shows that phage display, using a library of high diversity, is able to select high quality antibodies even when they are low in frequency.

Complexity of the human memory B-cell compartment is determined by the versatility of clonal diversification in germinal centers

Our knowledge about the clonal composition and intraclonal diversity of the human memory B-cell compartment and the relationship between memory B-cell subsets is still limited, although these are central issues for our understanding of adaptive immunity. We performed a deep sequencing analysis of rearranged immunoglobulin (Ig) heavy chain genes from biological replicates, covering more than 100,000 memory B lymphocytes from two healthy adults. We reveal a highly similar B-cell receptor repertoire among the four main human IgM(+) and IgG(+) memory B-cell subsets. Strikingly, in both donors, 45% of sequences could be assigned to expanded clones, demonstrating that the human memory B-cell compartment is characterized by many, often very large, B-cell clones. Twenty percent of the clones consisted of class switched and IgM(+)(IgD(+)) members, a feature that correlated significantly with clone size. Hence, we provide strong evidence that the vast majority of Ig mutated B cells--including IgM(+)IgD(+)CD27(+) B cells--are post-germinal center (GC) memory B cells. Clone members showed high intraclonal sequence diversity and high intraclonal versatility in Ig class and IgG subclass composition, with particular patterns of memory B-cell clone generation in GC reactions. In conclusion, GC produce amazingly large, complex, and diverse memory B-cell clones, equipping the human immune system with a versatile and highly diverse compartment of IgM(+)(IgD(+)) and class-switched memory B cells.

Systems integration of innate and adaptive immunity

The pathogens causing AIDS, malaria, and tuberculosis have proven too complex to be overcome by classical approaches to vaccination. The complexities of human immunology and pathogen-induced modulation of the immune system mandate new approaches to vaccine discovery and design. A new field, systems vaccinology, weds holistic analysis of innate and adaptive immunity within a quantitative framework to enable rational design of new vaccines that elicit tailored protective immune responses. A key step in the approach is to discover relationships between the earliest innate inflammatory responses to vaccination and the subsequent vaccine-induced adaptive immune responses and efficacy. Analysis of these responses in clinical studies is complicated by the inaccessibility of relevant tissue compartments (such as the lymph node), necessitating reliance upon peripheral blood responses as surrogates. Blood transcriptomes, although indirect to vaccine mechanisms, have proven very informative in systems vaccinology studies. The approach is most powerful when innate and adaptive immune responses are integrated with vaccine efficacy, which is possible for malaria with the advent of a robust human challenge model. This is more difficult for AIDS and tuberculosis, given that human challenge models are lacking and efficacy observed in clinical trials has been low or highly variable. This challenge can be met by appropriate clinical trial design for partially efficacious vaccines and by analysis of natural infection cohorts. Ultimately, systems vaccinology is an iterative approach in which mechanistic hypotheses-derived from analysis of clinical studies-are evaluated in model systems, and then used to guide the development of new vaccine strategies. In this review, we will illustrate the above facets of the systems vaccinology approach with case studies.

A bioinformatic framework for immune repertoire diversity profiling enables detection of immunological status

BACKGROUND

Lymphocyte receptor repertoires are continually shaped throughout the lifetime of an individual in response to environmental and pathogenic exposure. Thus, they may serve as a fingerprint of an individual's ongoing immunological status (e.g., healthy, infected, vaccinated), with far-reaching implications for immunodiagnostics applications. The advent of high-throughput immune repertoire sequencing now enables the interrogation of immune repertoire diversity in an unprecedented and quantitative manner. However, steadily increasing sequencing depth has revealed that immune repertoires vary greatly among individuals in their composition; correspondingly, it has been reported that there are few shared sequences indicative of immunological status ('public clones'). Disconcertingly, this means that the wealth of information gained from repertoire sequencing remains largely unused for determining the current status of immune responses, thereby hampering the implementation of immune-repertoire-based diagnostics.

METHODS

Here, we introduce a bioinformatics repertoire-profiling framework that possesses the advantage of capturing the diversity and distribution of entire immune repertoires, as opposed to singular public clones. The framework relies on Hill-based diversity profiles composed of a continuum of single diversity indices, which enable the quantification of the extent of immunological information contained in immune repertoires.

RESULTS

We coupled diversity profiles with unsupervised (hierarchical clustering) and supervised (support vector machine and feature selection) machine learning approaches in order to correlate patients' immunological statuses with their B- and T-cell repertoire data. We could predict with high accuracy (greater than or equal to 80 %) a wide range of immunological statuses such as healthy, transplantation recipient, and lymphoid cancer, suggesting as a proof of principle that diversity profiling can recover a large amount of immunodiagnostic fingerprints from immune repertoire data. Our framework is highly scalable as it easily allowed for the analysis of 1000 simulated immune repertoires; this exceeds the size of published immune repertoire datasets by one to two orders of magnitude.

CONCLUSIONS

Our framework offers the possibility to advance immune-repertoire-based fingerprinting, which may in the future enable a systems immunogenomics approach for vaccine profiling and the accurate and early detection of disease and infection.

Characterization of T cell repertoire of blood, tumor, and ascites in ovarian cancer patients using next generation sequencing

Tumor-infiltrating lymphocytes (TILs) play an important role in regulating the host immune response and are one of key factors in defining tumor microenvironment. Some studies have indicated that T cell infiltration in malignant ascites is associated with clinical outcome, but few studies have performed detailed characterization of T cell diversity or clonality in malignant effusions. We have applied a next generation sequencing method to characterize T cell repertoire of a set of primary cancers, ascites, and blood from 12 ovarian cancer patients and also analyzed the T cell subtype populations in malignant fluids from 3 ovarian cancer patients. We observed enrichment of certain T cells in tumors and ascites, but most of the enriched T cell receptor (TCR) sequences in tumors and ascites were not common. Moreover, we analyzed TCR sequences of T cell subtypes (CD4⁺, CD8⁺, and regulatory T cells) isolated from malignant effusions and also found clonal expansion of certain T cell populations, but the TCR sequences were almost mutually exclusive among the three subgroups. Although functional studies of clonally expanded T cell populations are definitely required, our approach offers a detailed characterization of T cell immune microenvironment in tumors and ascites that might differently affect antitumor immune response.

Systems-level analysis of innate immunity

Systems-level analysis of biological processes strives to comprehensively and quantitatively evaluate the interactions between the relevant molecular components over time, thereby enabling development of models that can be employed to ultimately predict behavior. Rapid development in measurement technologies (omics), when combined with the accessible nature of the cellular constituents themselves, is allowing the field of innate immunity to take significant strides toward this lofty goal. In this review, we survey exciting results derived from systems biology analyses of the immune system, ranging from gene regulatory networks to influenza pathogenesis and systems vaccinology.

Toward a more accurate view of human B-cell repertoire by next-generation sequencing, unbiased repertoire capture and single-molecule barcoding

B-cell repertoire analysis using next-generation sequencing has become a valuable tool for interrogating the genetic record of humoral response to infection. However, key obstacles such as low throughput, short read length, high error rate, and undetermined bias of multiplex PCR method have hindered broader application of this technology. In this study, we report several technical advances in antibody repertoire sequencing. We first demonstrated the ability to sequence antibody variable domains using the Ion Torrent PGM platform. As a test case, we analyzed the PGT121 class of antibodies from IAVI donor 17, an HIV-1-infected individual. We then obtained "unbiased" antibody repertoires by sequencing the 5'-RACE PCR products of B-cell transcripts from IAVI donor 17 and two HIV-1-uninfected individuals. We also quantified the bias of previously published gene-specific primers by comparing the repertoires generated by 5'-RACE PCR and multiplex PCR. We further developed a single-molecule barcoding strategy to reduce PCR-based amplification noise. Lastly, we evaluated several new PGM technologies in the context of antibody sequencing. We expect that, based upon long-read and high-fidelity next-generation sequencing technologies, the unbiased analysis will provide a more accurate view of the overall antibody repertoire while the barcoding strategy will facilitate high-resolution analysis of individual antibody families.

Assessing T cell clonal size distribution: a non-parametric approach

Clonal structure of the human peripheral T-cell repertoire is shaped by a number of homeostatic mechanisms, including antigen presentation, cytokine and cell regulation. Its accurate tuning leads to a remarkable ability to combat pathogens in all their variety, while systemic failures may lead to severe consequences like autoimmune diseases. Here we develop and make use of a non-parametric statistical approach to assess T cell clonal size distributions from recent next generation sequencing data. For 41 healthy individuals and a patient with ankylosing spondylitis, who undergone treatment, we invariably find power law scaling over several decades and for the first time calculate quantitatively meaningful values of decay exponent. It has proved to be much the same among healthy donors, significantly different for an autoimmune patient before the therapy, and converging towards a typical value afterwards. We discuss implications of the findings for theoretical understanding and mathematical modeling of adaptive immunity.

The impact of "omic" and imaging technologies on assessing the host immune response to biodefence agents

Understanding the interactions between host and pathogen is important for the development and assessment of medical countermeasures to infectious agents, including potential biodefence pathogens such as Bacillus anthracis, Ebola virus, and Francisella tularensis. This review focuses on technological advances which allow this interaction to be studied in much greater detail. Namely, the use of “omic” technologies (next generation sequencing, DNA, and protein microarrays) for dissecting the underlying host response to infection at the molecular level; optical imaging techniques (flow cytometry and fluorescence microscopy) for assessing cellular responses to infection; and biophotonic imaging for visualising the infectious disease process. All of these technologies hold great promise for important breakthroughs in the rational development of vaccines and therapeutics for biodefence agents.

Next generation sequencing reveals skewing of the T and B cell receptor repertoires in patients with wiskott-Aldrich syndrome

The Wiskott–Aldrich syndrome (WAS) is due to mutations of the WAS gene encoding for the cytoskeletal WAS protein, leading to abnormal downstream signaling from the T cell and B cell antigen receptors (TCR and BCR). We hypothesized that the impaired signaling through the TCR and BCR in WAS would subsequently lead to aberrations in the immune repertoire of WAS patients. Using next generation sequencing (NGS), the T cell receptor β and B cell immunoglobulin heavy chain (IGH) repertoires of eight patients with WAS and six controls were sequenced. Clonal expansions were identified within memory CD4⁺ cells as well as in total, naïve and memory CD8⁺ cells from WAS patients. In the B cell compartment, WAS patient IGH repertoires were also clonally expanded and showed skewed usage of IGHV and IGHJ genes, and increased usage of IGHG constant genes, compared with controls. To our knowledge, this is the first study that demonstrates significant abnormalities of the immune repertoire in WAS patients using NGS.

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.

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.

Qualifying high-throughput immune repertoire sequencing

Diversity of B and T cell receptors, achieved by gene recombination and somatic hypermutation, allows the immune system for recognition and targeted reaction against various threats. Next-generation sequencing for assessment of a cell's gene composition and variation makes deep analysis of one individual's immune spectrum feasible. An easy to apply but detailed analysis and visualization strategy is necessary to process all sequences generated. We performed sequencing utilizing the 454 system for CLL and control samples, utilized the IMGT database and applied the presented analysis tools. With the applied protocol, malignant clones are found and characterized, mutational status compared to germline identity is elaborated in detail showing that the CLL mutation status is not as monoclonal as generally thought. On the other hand, this strategy is not solely applicable to the 454 sequencing system but can easily be transferred to any other next-generation sequencing platform.

Correlation of contact sensitizer potency with T cell frequency and TCR repertoire diversity

Allergic contact dermatitis is a T cell-mediated skin disease. Many hundreds of organic chemicals and some metal ions are contact sensitizers. They induce an innate inflammatory immune response in the skin that results in the priming of contact sensitizer-specific T cells by dendritic cells in the draining lymph nodes. The factors that determine the strength of this T cell response and thereby define the potency of a contact sensitizer are largely unknown. This chapter highlights different variables such as precursor frequency of antigen-specific T cells, possible bystander activation, and T cell receptor diversity or avidity of the TCR/peptide-MHC interactions, which might impact the quality and strength of T cell responses to contact sensitizers. In addition, different methods available to determine both the frequency of antigen-specific T cells and T cell receptor repertoires are discussed. Identification of the factors determining potency may allow for the development of suitable in vitro assays for potency assessment of contact sensitizers.

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.

Preparing unbiased T-cell receptor and antibody cDNA libraries for the deep next generation sequencing profiling

High-throughput sequencing has the power to reveal the nature of adaptive immunity as represented by the full complexity of T-cell receptor (TCR) and antibody (IG) repertoires, but is at present severely compromised by the quantitative bias, bottlenecks, and accumulated errors that inevitably occur in the course of library preparation and sequencing. Here we report an optimized protocol for the unbiased preparation of TCR and IG cDNA libraries for high-throughput sequencing, starting from thousands or millions of live cells in an investigated sample. Critical points to control are revealed, along with tips that allow researchers to minimize quantitative bias, accumulated errors, and cross-sample contamination at each stage, and to enhance the subsequent bioinformatic analysis. The protocol is simple, reliable, and can be performed in 1–2 days.

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.

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.