Statistical inference of a convergent antibody repertoire response to influenza vaccine

Extraction of the CDRH3 sequence of the mouse antibody repertoire selected upon influenza virus infection by subtraction of the background antibody repertoire

Historically, antibody reactivity to pathogens and vaccine antigens has been evaluated using serological measurements of antigen-specific antibodies. However, it is difficult to evaluate all antibodies that contribute to various functions in a single assay, such as the measurement of the neutralizing antibody titer. Bulk antibody repertoire analysis using next-generation sequencing is a comprehensive method for analyzing the overall antibody response; however, it is unreliable for estimating antigen-specific antibodies due to individual variation. To address this issue, we propose a method to subtract the background signal from the repertoire of data of interest. In this study, we analyzed changes in antibody diversity and inferred the heavy-chain complementarity-determining region 3 (CDRH3) sequences of antibody clones that were selected upon influenza virus infection in a mouse model using bulk repertoire analysis. A decrease in the diversity of the antibody repertoire was observed upon viral infection, along with an increase in neutralizing antibody titers. Using kernel density estimation of sequences in a high-dimensional sequence space with background signal subtraction, we identified several clusters of CDRH3 sequences induced upon influenza virus infection. Most of these repertoires were detected more frequently in infected mice than in uninfected control mice, suggesting that infection-specific antibody sequences can be extracted using this method. Such an accurate extraction of antigen- or infection-specific repertoire information will be a useful tool for vaccine evaluation in the future.

IMPORTANCE

As specific interactions between antigens and cell-surface antibodies trigger the proliferation of B-cell clones, the frequency of each antibody sequence in the samples reflects the size of each clonal population. Nevertheless, it is extremely difficult to extract antigen-specific antibody sequences from the comprehensive bulk antibody sequences obtained from blood samples due to repertoire bias influenced by exposure to dietary antigens and other infectious agents. This issue can be addressed by subtracting the background noise from the post-immunization or post-infection repertoire data. In the present study, we propose a method to quantify repertoire data from comprehensive repertoire data. This method allowed subtraction of the background repertoire, resulting in more accurate extraction of expanded antibody repertoires upon influenza virus infection. This accurate extraction of antigen- or infection-specific repertoire information is a useful tool for vaccine evaluation.

Effect of Ocrelizumab on B- and T-Cell Receptor Repertoire Diversity in Patients With Relapsing Multiple Sclerosis From the Randomized Phase III OPERA Trial

BACKGROUND AND OBJECTIVES

The B cell-depleting anti-CD20 antibody ocrelizumab (OCR) effectively reduces MS disease activity and slows disability progression. Given the role of B cells as antigen-presenting cells, the primary goal of this study was to evaluate the effect of OCR on the T-cell receptor repertoire diversity.

METHODS

To examine whether OCR substantially alters the molecular diversity of the T-cell receptor repertoire, deep immune repertoire sequencing (RepSeq) of CD4+ and CD8+ T-cell receptor β-chain variable regions was performed on longitudinal blood samples. The IgM and IgG heavy chain variable region repertoire was also analyzed to characterize the residual B-cell repertoire under OCR treatment.

RESULTS

Peripheral blood samples for RepSeq were obtained from 8 patients with relapsing MS enrolled in the OPERA I trial over a period of up to 39 months. Four patients each were treated with OCR or interferon β1-a during the double-blind period of OPERA I. All patients received OCR during the open-label extension. The diversity of the CD4+/CD8+ T-cell repertoires remained unaffected in OCR-treated patients. The expected OCR-associated B-cell depletion was mirrored by reduced B-cell receptor diversity in peripheral blood and a shift in immunoglobulin gene usage. Despite deep B-cell depletion, longitudinal persistence of clonally related B-cells was observed.

DISCUSSION

Our data illustrate that the diversity of CD4+/CD8+ T-cell receptor repertoires remained unaltered in OCR-treated patients with relapsing MS. Persistence of a highly diverse T-cell repertoire suggests that aspects of adaptive immunity remain intact despite extended anti-CD20 therapy.

TRIAL REGISTRATION INFORMATION

This is a substudy (BE29353) of the OPERA I (WA21092; NCT01247324) trial. Date of registration, November 23, 2010; first patient enrollment, August 31, 2011.

Reference-based comparison of adaptive immune receptor repertoires

B and T cell receptor (immune) repertoires can represent an individual's immune history. While current repertoire analysis methods aim to discriminate between health and disease states, they are typically based on only a limited number of parameters. Here, we introduce immuneREF: a quantitative multidimensional measure of adaptive immune repertoire (and transcriptome) similarity that allows interpretation of immune repertoire variation by relying on both repertoire features and cross-referencing of simulated and experimental datasets. To quantify immune repertoire similarity landscapes across health and disease, we applied immuneREF to >2,400 datasets from individuals with varying immune states (healthy, [autoimmune] disease, and infection). We discovered, in contrast to the current paradigm, that blood-derived immune repertoires of healthy and diseased individuals are highly similar for certain immune states, suggesting that repertoire changes to immune perturbations are less pronounced than previously thought. In conclusion, immuneREF enables the population-wide study of adaptive immune response similarity across immune states.

B-Cell Receptor Repertoire Sequencing: Deeper Digging into the Mechanisms and Clinical Aspects of Immune-mediated Diseases

B cells play an essential role in adaptive immunity and are intimately correlated with pleiotropic immune-mediated diseases. Each B cell occupies a unique B cell receptor (BCR), and all BCRs throughout our body form “BCR repertoire.” With the development of sequencing technology and coupled bioinformatics, accumulating evidence indicates that BCR repertoire largely varies under physiological and pathological conditions. Therefore, comprehensive grasp of BCR repertoire will provide new insights into the pathogenesis of immune-mediated diseases and help exploit efficient diagnostic and treatment strategies. In this review, we start with an overview of BCR repertoire and related sequencing technologies and summarize their current applications in immune-mediated diseases. We also underscore the challenges of this emerging field and propose promising future directions in advancing BCR repertoire exploration.

A Pan-Cancer Analysis of Tumor-Infiltrating B Cell Repertoires

Tumor-infiltrating B cells can play an important role in anti-tumor responses but their presence is not well understood. In this study, we extracted the B cell receptor repertoires from 9522 tumor and adjacent non-tumor samples across 28 tumor types in the Cancer Genome Atlas project and performed diversity and network analysis. We identified differences in diversity and network statistics across tumor types and subtypes and observed a trend towards increased clonality in primary tumors compared to adjacent non-tumor tissues. We also found significant associations between the repertoire features and mutation load, tumor stage, and age. Our V-gene usage analysis identified similar V-gene usage patterns in colorectal and endometrial cancers. Lastly, we evaluated the prognostic value of the repertoire features and identified significant associations with survival in seven tumor types. This study warrants further research into better understanding the role of tumor-infiltrating B cells across a wide range of tumor types.

The Future of Blood Testing Is the Immunome

It is increasingly clear that an extraordinarily diverse range of clinically important conditions—including infections, vaccinations, autoimmune diseases, transplants, transfusion reactions, aging, and cancers—leave telltale signatures in the millions of V(D)J-rearranged antibody and T cell receptor [TR per the Human Genome Organization (HUGO) nomenclature but more commonly known as TCR] genes collectively expressed by a person’s B cells (antibodies) and T cells. We refer to these as the immunome. Because of its diversity and complexity, the immunome provides singular opportunities for advancing personalized medicine by serving as the substrate for a highly multiplexed, near-universal blood test. Here we discuss some of these opportunities, the current state of immunome-based diagnostics, and highlight some of the challenges involved. We conclude with a call to clinicians, researchers, and others to join efforts with the Adaptive Immune Receptor Repertoire Community (AIRR-C) to realize the diagnostic potential of the immunome.

Convergent antibody evolution and clonotype expansion following influenza virus vaccination

Recent advances in high-throughput single cell sequencing have opened up new avenues into the investigation of B cell receptor (BCR) repertoires. In this study, PBMCs were collected from 17 human participants vaccinated with the split-inactivated influenza virus vaccine during the 2016-2017 influenza season. A combination of Immune Repertoire Capture (IRCTM) technology and IgG sequencing was performed on ~7,800 plasmablast (PB) cells and preferential IgG heavy-light chain pairings were investigated. In some participants, a single expanded clonotype accounted for ~22% of their PB BCR repertoire. Approximately 60% (10/17) of participants experienced convergent evolution, possessing public PBs that were elicited independently in multiple participants. Binding profiles of one private and three public PBs confirmed they were all subtype-specific, cross-reactive hemagglutinin (HA) head-directed antibodies. Collectively, this high-resolution antibody repertoire analysis demonstrated the impact evolution can have on BCRs in response to influenza virus vaccination, which can guide future universal influenza prophylactic approaches.

Data Integration for Immunology

Over the last several years, next-generation sequencing and its recent push toward single-cell resolution have transformed the landscape of immunology research by revealing novel complexities about all components of the immune system. With the vast amounts of diverse data currently being generated, and with the methods of analyzing and combining diverse data improving as well, integrative systems approaches are becoming more powerful. Previous integrative approaches have combined multiple data types and revealed ways that the immune system, both as a whole and as individual parts, is affected by genetics, the microbiome, and other factors. In this review, we explore the data types that are available for studying immunology with an integrative systems approach, as well as the current strategies and challenges for conducting such analyses.

Profiling immunoglobulin repertoires across multiple human tissues using RNA sequencing

Profiling immunoglobulin (Ig) receptor repertoires with specialized assays can be cost-ineffective and time-consuming. Here we report ImReP, a computational method for rapid and accurate profiling of the Ig repertoire, including the complementary-determining region 3 (CDR3), using regular RNA sequencing data such as those from 8,555 samples across 53 tissues types from 544 individuals in the Genotype-Tissue Expression (GTEx v6) project. Using ImReP and GTEx v6 data, we generate a collection of 3.6 million Ig sequences, termed the atlas of immunoglobulin repertoires (TAIR), across a broad range of tissue types that often do not have reported Ig repertoires information. Moreover, the flow of Ig clonotypes and inter-tissue repertoire similarities across immune-related tissues are also evaluated. In summary, TAIR is one of the largest collections of CDR3 sequences and tissue types, and should serve as an important resource for studying immunological diseases.

Information on immune receptor repertoire provides important insights on disease progression and therapy development, but can be expensive and time-consuming to obtain. Here the authors report ImReP, a computational method that can extract detailed immune repertoire information from existing tissue-specific RNA sequencing data.

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.

Inferring B cell specificity for vaccines using a Bayesian mixture model

Background

Vaccines have greatly reduced the burden of infectious disease, ranking in their impact on global health second only after clean water. Most vaccines confer protection by the production of antibodies with binding affinity for the antigen, which is the main effector function of B cells. This results in short term changes in the B cell receptor (BCR) repertoire when an immune response is launched, and long term changes when immunity is conferred. Analysis of antibodies in serum is usually used to evaluate vaccine response, however this is limited and therefore the investigation of the BCR repertoire provides far more detail for the analysis of vaccine response.

Results

Here, we introduce a novel Bayesian model to describe the observed distribution of BCR sequences and the pattern of sharing across time and between individuals, with the goal to identify vaccine-specific BCRs. We use data from two studies to assess the model and estimate that we can identify vaccine-specific BCRs with 69% sensitivity.

Conclusion

Our results demonstrate that statistical modelling can capture patterns associated with vaccine response and identify vaccine specific B cells in a range of different data sets. Additionally, the B cells we identify as vaccine specific show greater levels of sequence similarity than expected, suggesting that there are additional signals of vaccine response, not currently considered, which could improve the identification of vaccine specific B cells.

Characterizing pre-transplant and post-transplant kidney rejection risk by B cell immune repertoire sequencing

Studying immune repertoire in the context of organ transplant provides important information on how adaptive immunity may contribute and modulate graft rejection. Here we characterize the peripheral blood immune repertoire of individuals before and after kidney transplant using B cell receptor sequencing in a longitudinal clinical study. Individuals who develop rejection after transplantation have a more diverse immune repertoire before transplant, suggesting a predisposition for post-transplant rejection risk. Additionally, over 2 years of follow-up, patients who develop rejection demonstrate a specific set of expanded clones that persist after the rejection. While there is an overall reduction of peripheral B cell diversity, likely due to increased general immunosuppression exposure in this cohort, the detection of specific IGHV gene usage across all rejecting patients supports that a common pool of immunogenic antigens may drive post-transplant rejection. Our findings may have clinical implications for the prediction and clinical management of kidney transplant rejection.

Current and future influenza vaccines

Although antiviral drugs and vaccines have reduced the economic and healthcare burdens of influenza, influenza epidemics continue to take a toll. Over the past decade, research on influenza viruses has revealed a potential path to improvement. The clues have come from accumulated discoveries from basic and clinical studies. Now, virus surveillance allows researchers to monitor influenza virus epidemic trends and to accumulate virus sequences in public databases, which leads to better selection of candidate viruses for vaccines and early detection of drug-resistant viruses. Here we provide an overview of current vaccine options and describe efforts directed toward the development of next-generation vaccines. Finally, we propose a plan for the development of an optimal influenza vaccine.

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

Engineering Improved Antiphosphotyrosine Antibodies Based on an Immunoconvergent Binding Motif

Phosphotyrosine (pY) is one of the most highly studied posttranslational modifications that is responsible for tightly regulating many signaling pathways in eukaryotes. Pan-specific pY antibodies have emerged as powerful tools for understanding the role of these modifications. Nevertheless, structures have not been reported for pan-specific pY antibodies, greatly impeding the further development of tools for integrating this ubiquitous posttranslational modification using structure-guided designs. Here, we present the first crystal structures of two widely utilized pan-specific pY antibodies, PY20 and 4G10. The two antibodies, although developed independently from animal immunizations, have surprisingly similar modes of recognition of the phosphate group, implicating a generic binding structure among pan-specific pY antibodies. Sequence alignments revealed that many pY binding residues are predominant in the mouse V germline genes, which consequently led to the convergent antibodies. On the basis of the convergent structure, we designed a phage display library by lengthening the CDR-L3 loop with the aid of computational modeling. Panning with this library resulted in a series of 4G10 variants with 4 to 11-fold improvements in pY binding affinities. The crystal structure of one improved variant showed remarkable superposition to the computational model, where the lengthened CDR-L3 loop creates an additional hydrogen bond indirectly bound to the phosphate group via a water molecule. The engineered variants exhibited superior performance in Western blot and immunofluorescence.

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.

Immunoglobulin gene analysis as a tool for investigating human immune responses

The human immunoglobulin repertoire is a hugely diverse set of sequences that are formed by processes of gene rearrangement, heavy and light chain gene assortment, class switching and somatic hypermutation. Early B cell development produces diverse IgM and IgD B cell receptors on the B cell surface, resulting in a repertoire that can bind many foreign antigens but which has had self-reactive B cells removed. Later antigen-dependent development processes adjust the antigen affinity of the receptor by somatic hypermutation. The effector mechanism of the antibody is also adjusted, by switching the class of the antibody from IgM to one of seven other classes depending on the required function. There are many instances in human biology where positive and negative selection forces can act to shape the immunoglobulin repertoire and therefore repertoire analysis can provide useful information on infection control, vaccination efficacy, autoimmune diseases, and cancer. It can also be used to identify antigen-specific sequences that may be of use in therapeutics. The juxtaposition of lymphocyte development and numerical evaluation of immune repertoires has resulted in the growth of a new sub-speciality in immunology where immunologists and computer scientists/physicists collaborate to assess immune repertoires and develop models of immune action.

VDJServer: A Cloud-Based Analysis Portal and Data Commons for Immune Repertoire Sequences and Rearrangements

Background

Recent technological advances in immune repertoire sequencing have created tremendous potential for advancing our understanding of adaptive immune response dynamics in various states of health and disease. Immune repertoire sequencing produces large, highly complex data sets, however, which require specialized methods and software tools for their effective analysis and interpretation.

Results

VDJServer is a cloud-based analysis portal for immune repertoire sequence data that provide access to a suite of tools for a complete analysis workflow, including modules for preprocessing and quality control of sequence reads, V(D)J gene segment assignment, repertoire characterization, and repertoire comparison. VDJServer also provides sophisticated visualizations for exploratory analysis. It is accessible through a standard web browser via a graphical user interface designed for use by immunologists, clinicians, and bioinformatics researchers. VDJServer provides a data commons for public sharing of repertoire sequencing data, as well as private sharing of data between users. We describe the main functionality and architecture of VDJServer and demonstrate its capabilities with use cases from cancer immunology and autoimmunity.

Conclusion

VDJServer provides a complete analysis suite for human and mouse T-cell and B-cell receptor repertoire sequencing data. The combination of its user-friendly interface and high-performance computing allows large immune repertoire sequencing projects to be analyzed with no programming or software installation required. VDJServer is a web-accessible cloud platform that provides access through a graphical user interface to a data management infrastructure, a collection of analysis tools covering all steps in an analysis, and an infrastructure for sharing data along with workflows, results, and computational provenance. VDJServer is a free, publicly available, and open-source licensed resource.

ImmPort, toward repurposing of open access immunological assay data for translational and clinical research

Immunology researchers are beginning to explore the possibilities of reproducibility, reuse and secondary analyses of immunology data. Open-access datasets are being applied in the validation of the methods used in the original studies, leveraging studies for meta-analysis, or generating new hypotheses. To promote these goals, the ImmPort data repository was created for the broader research community to explore the wide spectrum of clinical and basic research data and associated findings. The ImmPort ecosystem consists of four components–Private Data, Shared Data, Data Analysis, and Resources—for data archiving, dissemination, analyses, and reuse. To date, more than 300 studies have been made freely available through the Shared Data portal (www.immport.org/immport-open), which allows research data to be repurposed to accelerate the translation of new insights into discoveries.

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.

ROP: dumpster diving in RNA-sequencing to find the source of 1 trillion reads across diverse adult human tissues

High-throughput RNA-sequencing (RNA-seq) technologies provide an unprecedented opportunity to explore the individual transcriptome. Unmapped reads are a large and often overlooked output of standard RNA-seq analyses. Here, we present Read Origin Protocol (ROP), a tool for discovering the source of all reads originating from complex RNA molecules. We apply ROP to samples across 2630 individuals from 54 diverse human tissues. Our approach can account for 99.9% of 1 trillion reads of various read length. Additionally, we use ROP to investigate the functional mechanisms underlying connections between the immune system, microbiome, and disease. ROP is freely available at https://github.com/smangul1/rop/wiki.

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.

Quantitative Analysis of Repertoire-Scale Immunoglobulin Properties in Vaccine-Induced B-Cell Responses

Recent advances in the next-generation sequencing of B-cell receptors (BCRs) enable the characterization of humoral responses at a repertoire-wide scale and provide the capability for identifying unique features of immune repertoires in response to disease, vaccination, or infection. Immunosequencing now readily generates 10³–10⁵ sequences per sample; however, statistical analysis of these repertoires is challenging because of the high genetic diversity of BCRs and the elaborate clonal relationships among them. To date, most immunosequencing analyses have focused on reporting qualitative trends in immunoglobulin (Ig) properties, such as usage or somatic hypermutation (SHM) percentage of the Ig heavy chain variable (IGHV) gene segment family, and on reducing complex Ig property distributions to simple summary statistics. However, because Ig properties are typically not normally distributed, any approach that fails to assess the distribution as a whole may be inadequate in (1) properly assessing the statistical significance of repertoire differences, (2) identifying how two repertoires differ, and (3) determining appropriate confidence intervals for assessing the size of the differences and their potential biological relevance. To address these issues, we have developed a technique that uses Wilcox’ robust statistics toolbox to identify statistically significant vaccine-specific differences between Ig repertoire properties. The advantage of this technique is that it can determine not only whether but also where the distributions differ, even when the Ig repertoire properties are non-normally distributed. We used this technique to characterize murine germinal center (GC) B-cell repertoires in response to a complex Ebola virus-like particle (eVLP) vaccine candidate with known protective efficacy. The eVLP-mediated GC B-cell responses were highly diverse, consisting of thousands of clonotypes. Despite this staggering diversity, we identified statistically significant differences between non-immunized, vaccine only, and vaccine-plus-adjuvant groups in terms of Ig properties, including IGHV-family usage, SHM percentage, and characteristics of the BCR complementarity-determining region. Most notably, our analyses identified a robust eVLP-specific feature—enhanced IGHV8-family usage in B-cell repertoires. These findings demonstrate the utility of our technique in identifying statistically significant BCR repertoire differences following vaccination. More generally, our approach is potentially applicable to a wide range of studies in infection, vaccination, auto-immunity, and cancer.

The Individual and Population Genetics of Antibody Immunity

Antibodies (Abs) produced by immunoglobulin (IG) genes are the most diverse proteins expressed in humans. While part of this diversity is generated by recombination during B-cell development and mutations during affinity maturation, the germ-line IG loci are also diverse across human populations and ethnicities. Recently, proof-of-concept studies have demonstrated genotype–phenotype correlations between specific IG germ-line variants and the quality of Ab responses during vaccination and disease. However, the functional consequences of IG genetic variation in Ab function and immunological outcomes remain underexplored. In this opinion article, we outline interconnections between IG genomic diversity and Ab-expressed repertoires and structure. We further propose a strategy for integrating IG genotyping with functional Ab profiling data as a means to better predict and optimize humoral responses in genetically diverse human populations, with immediate implications for personalized medicine.

Genetic variation in human populations affects how individuals are able to mount functional antibody responses.

Different alleles can encode convergent binding motifs that result in successful Ab responses against specific infections and vaccinations.

Given the complexity of the IG loci and the diversity of the antibody repertoire, links between IG polymorphism and antibody repertoire variability have not been thoroughly explored.

We present a strategy to mine genotype–repertoire–disease associations.

A Public Database of Memory and Naive B-Cell Receptor Sequences

The vast diversity of B-cell receptors (BCR) and secreted antibodies enables the recognition of, and response to, a wide range of epitopes, but this diversity has also limited our understanding of humoral immunity. We present a public database of more than 37 million unique BCR sequences from three healthy adult donors that is many fold deeper than any existing resource, together with a set of online tools designed to facilitate the visualization and analysis of the annotated data. We estimate the clonal diversity of the naive and memory B-cell repertoires of healthy individuals, and provide a set of examples that illustrate the utility of the database, including several views of the basic properties of immunoglobulin heavy chain sequences, such as rearrangement length, subunit usage, and somatic hypermutation positions and dynamics.