Quantitative assessment of the robustness of next-generation sequencing of antibody variable gene repertoires from immunized mice

Bulk T cell repertoire sequencing (TCR-Seq) is a powerful technology for understanding inflammation-mediated diseases

Multiple alterations in the T cell repertoire were identified in many chronic inflammatory diseases such as inflammatory bowel disease, multiple sclerosis, and rheumatoid arthritis, suggesting that T cells might, directly or indirectly, be implicated in these pathologies. This has sparked a deep interest in characterizing disease-associated T cell clonotypes as well as in identifying and quantifying their contribution to the pathophysiology of different autoimmune and inflammation-mediated diseases. Bulk T cell repertoire sequencing (TCR-Seq) has emerged as a powerful method to profile the T cell repertoire of a sample in a high throughput fashion. Given the increasing utilization of TCR-Seq, we aimed here to provide a comprehensive, up-to-date review of the method, its extensions, and its ability to investigate chronic and autoimmune diseases. Specifically, we started by introducing the immunological basis of TCR repertoire generation and features, followed by discussing different experimental approach to perform TCR-Seq, then we describe different methods and frameworks for analyzing the generated datasets. Subsequently, different experimental techniques for investigating the antigenicity of T cell clonotypes are described. Lastly, we discuss recent studies that utilized TCR-Seq to understand different inflammation-mediated diseases, discuss fallbacks of the technology and potential future directions to overcome current limitations.

Machine learning for precision diagnostics of autoimmunity

Early and accurate diagnosis is crucial to prevent disease development and define therapeutic strategies. Due to predominantly unspecific symptoms, diagnosis of autoimmune diseases (AID) is notoriously challenging. Clinical decision support systems (CDSS) are a promising method with the potential to enhance and expedite precise diagnostics by physicians. However, due to the difficulties of integrating and encoding multi-omics data with clinical values, as well as a lack of standardization, such systems are often limited to certain data types. Accordingly, even sophisticated data models fall short when making accurate disease diagnoses and presenting data analyses in a user-friendly form. Therefore, the integration of various data types is not only an opportunity but also a competitive advantage for research and industry. We have developed an integration pipeline to enable the use of machine learning for patient classification based on multi-omics data in combination with clinical values and laboratory results. The application of our framework resulted in up to 96% prediction accuracy of autoimmune diseases with machine learning models. Our results deliver insights into autoimmune disease research and have the potential to be adapted for applications across disease conditions.

Enhancing comparative T cell receptor repertoire analysis in small biological samples through pooling homologous cell samples from multiple mice

Accurate characterization and comparison of T cell receptor (TCR) repertoires from small biological samples present significant challenges. The main challenge is the low material input, which compromises the quality of bulk sequencing and hinders the recovery of sufficient TCR sequences for robust analyses. We aimed to address this limitation by implementing a strategic approach to pool homologous biological samples. Our findings demonstrate that such pooling indeed enhances the TCR repertoire coverage, particularly for cell subsets of constrained sizes, and enables accurate comparisons of TCR repertoires at different levels of complexity across T cell subsets with different sizes. This methodology holds promise for advancing our understanding of T cell repertoires in scenarios where sample size constraints are a prevailing concern.

Guidelines for reproducible analysis of adaptive immune receptor repertoire sequencing data

Enhancing the reproducibility and comprehension of adaptive immune receptor repertoire sequencing (AIRR-seq) data analysis is critical for scientific progress. This study presents guidelines for reproducible AIRR-seq data analysis, and a collection of ready-to-use pipelines with comprehensive documentation. To this end, ten common pipelines were implemented using ViaFoundry, a user-friendly interface for pipeline management and automation. This is accompanied by versioned containers, documentation and archiving capabilities. The automation of pre-processing analysis steps and the ability to modify pipeline parameters according to specific research needs are emphasized. AIRR-seq data analysis is highly sensitive to varying parameters and setups; using the guidelines presented here, the ability to reproduce previously published results is demonstrated. This work promotes transparency, reproducibility, and collaboration in AIRR-seq data analysis, serving as a model for handling and documenting bioinformatics pipelines in other research domains.

Application and challenges of TCR and BCR sequencing to investigate T- and B-cell clonality in elastase-induced experimental murine abdominal aortic aneurysm

Background

An abdominal aortic aneurysm (AAA) is a life-threatening cardiovascular disease. Although its pathogenesis is still poorly understood, recent evidence suggests that AAA displays autoimmune disease characteristics. Particularly, T cells responding to AAA-related antigens in the aortic wall may contribute to an initial immune response. Single-cell RNA (scRNA) T cell receptor (TCR) and B cell receptor (BCR) sequencing is a powerful tool for investigating clonality. However, difficulties such as limited numbers of isolated cells must be considered during implementation and data analysis, making biological interpretation challenging. Here, we perform a representative single-cell immune repertoire analysis in experimental murine AAA and show a reliable bioinformatic processing pipeline highlighting opportunities and limitations of this approach.

Methods

We performed scRNA TCR and BCR sequencing of isolated lymphocytes from the infrarenal aorta of male C57BL/6J mice 3, 7, 14, and 28 days after AAA induction via elastase perfusion of the aorta. Sham-operated mice at days 3 and 28 and non-operated mice served as controls.

Results

Comparison of complementarity-determining region (CDR3) length distribution of 179 B cells and 796 T cells revealed neither differences between AAA and control nor between the disease stages. We found no clonal expansion of B cells in AAA. For T cells, we identified several clones in 11 of 16 AAA samples and one of eight control samples. Immune receptor repertoire comparison indicated that only a few clones were shared between the individual AAA samples. The most frequently used V-genes in the TCR beta chain in AAA were TRBV3, TRBV19, and the splicing variant TRBV12-2 + TRBV13-2.

Conclusion

We found no clonal expansion of B cells but evidence for clonal expansion of T cells in elastase-induced AAA in mice. Our findings imply that a more precise characterization of TCR and BCR distribution requires a more extensive number of lymphocytes to prevent undersampling and potentially detect rare clones. Thus, further experiments are necessary to confirm our findings. In summary, this paper examines TCR and BCR sequencing results, identifies limitations and pitfalls, and offers guidance for future studies.

Apyrase-mediated amplification of secretory IgA promotes intestinal homeostasis

Secretory immunoglobulin A (SIgA) interaction with commensal bacteria conditions microbiota composition and function. However, mechanisms regulating reciprocal control of microbiota and SIgA are not defined. Bacteria-derived adenosine triphosphate (ATP) limits T follicular helper (Tfh) cells in the Peyer's patches (PPs) via P2X7 receptor (P2X7R) and thereby SIgA generation. Here we show that hydrolysis of extracellular ATP (eATP) by apyrase results in amplification of the SIgA repertoire. The enhanced breadth of SIgA in mice colonized with apyrase-releasing Escherichia coli influences topographical distribution of bacteria and expression of genes involved in metabolic versus immune functions in the intestinal epithelium. SIgA-mediated conditioning of bacteria and enterocyte function is reflected by differences in nutrient absorption in mice colonized with apyrase-expressing bacteria. Apyrase-induced SIgA improves intestinal homeostasis and attenuates barrier impairment and susceptibility to infection by enteric pathogens in antibiotic-induced dysbiosis. Therefore, amplification of SIgA by apyrase can be leveraged to restore intestinal fitness in dysbiotic conditions.

Characterization of B cell receptor H-CDR3 repertoire of spleen in PRV-infected mice

Pseudorabies virus (PRV), also known as suid Alphaherpesvirus 1 (SuHV-1), which is one of the most devastating infectious pathogen of swine industry worldwide. Vaccination is the safest and most effective PRV prevention and control strategy. B cell receptor (BCR) is membrane-bound immunoglobulin located on the surface of B cells capable of specifically binding foreign antigens, which is one of the most important molecules regulating the proliferation and function of B cells. Here, to assess the molecular diversity of BCR H-CDR3 repertoire after different PRV strains infection, we detected the IGHV, IGHD, IGHJ genes usage and CDR3 sequence changes of mice spleen with PRV vaccine strain (Bartha-K61), variant strain (XJ) and mock infection by high-throughput sequencing. We found that PRV-infected groups shared partial BCR sequences, which are most likely to be PRV-specific BCR candidates. However, there were still differences in the IGHV genes usage as well as the combined usage of IGHV and IGHJ genes between the Bartha-K61 strain and XJ strain infection groups. In addition, the CDR3 sequences exhibited large differences in the types and lengths in PRV infection groups. Our study contributes to a better understanding of the host adaptive immune response to PRV infection and provides a theoretical basis for further research on novel and efficient PRV vaccines in the future.

TCRpower: quantifying the detection power of T-cell receptor sequencing with a novel computational pipeline calibrated by spike-in sequences

T-cell receptor (TCR) sequencing has enabled the development of innovative diagnostic tests for cancers, autoimmune diseases and other applications. However, the rarity of many T-cell clonotypes presents a detection challenge, which may lead to misdiagnosis if diagnostically relevant TCRs remain undetected. To address this issue, we developed TCRpower, a novel computational pipeline for quantifying the statistical detection power of TCR sequencing methods. TCRpower calculates the probability of detecting a TCR sequence as a function of several key parameters: in-vivo TCR frequency, T-cell sample count, read sequencing depth and read cutoff. To calibrate TCRpower, we selected unique TCRs of 45 T-cell clones (TCCs) as spike-in TCRs. We sequenced the spike-in TCRs from TCCs, together with TCRs from peripheral blood, using a 5' RACE protocol. The 45 spike-in TCRs covered a wide range of sample frequencies, ranging from 5 per 100 to 1 per 1 million. The resulting spike-in TCR read counts and ground truth frequencies allowed us to calibrate TCRpower. In our TCR sequencing data, we observed a consistent linear relationship between sample and sequencing read frequencies. We were also able to reliably detect spike-in TCRs with frequencies as low as one per million. By implementing an optimized read cutoff, we eliminated most of the falsely detected sequences in our data (TCR α-chain 99.0% and TCR β-chain 92.4%), thereby improving diagnostic specificity. TCRpower is publicly available and can be used to optimize future TCR sequencing experiments, and thereby enable reliable detection of disease-relevant TCRs for diagnostic applications.

Adaptive Immune Receptor Repertoire (AIRR) Community Guide to Repertoire Analysis

Adaptive immune receptor repertoires (AIRRs) are rich with information that can be mined for insights into the workings of the immune system. Gene usage, CDR3 properties, clonal lineage structure, and sequence diversity are all capable of revealing the dynamic immune response to perturbation by disease, vaccination, or other interventions. Here we focus on a conceptual introduction to the many aspects of repertoire analysis and orient the reader toward the uses and advantages of each. Along the way, we note some of the many software tools that have been developed for these investigations and link the ideas discussed to chapters on methods provided elsewhere in this volume.

An Analysis of the Effects of Spaceflight and Vaccination on Antibody Repertoire Diversity

Ab repertoire diversity plays a critical role in the host's ability to fight pathogens. CDR3 is partially responsible for Ab-Ag binding and is a significant source of diversity in the repertoire. CDR3 diversity is generated during VDJ rearrangement because of gene segment selection, gene segment trimming and splicing, and the addition of nucleotides. We analyzed the Ab repertoire diversity across multiple experiments examining the effects of spaceflight on the Ab repertoire after vaccination. Five datasets from four experiments were analyzed using rank-abundance curves and Shannon indices as measures of diversity. We discovered a trend toward lower diversity as a result of spaceflight but did not find the same decrease in our physiological model of microgravity in either the spleen or bone marrow. However, the bone marrow repertoire showed a reduction in diversity after vaccination. We also detected differences in Shannon indices between experiments and tissues. We did not detect a pattern of CDR3 usage across the experiments. Overall, we were able to find differences in the Ab repertoire diversity across experimental groups and tissues.

Antibody Repertoire Analysis of Tumor-Infiltrating B Cells Reveals Distinct Signatures and Distributions Across Tissues

The role of B cells in the tumor microenvironment (TME) has largely been under investigated, and data regarding the antibody repertoire encoded by B cells in the TME and the adjacent lymphoid organs are scarce. Here, we utilized B cell receptor high-throughput sequencing (BCR-Seq) to profile the antibody repertoire signature of tumor-infiltrating lymphocyte B cells (TIL−Bs) in comparison to B cells from three anatomic compartments in a mouse model of triple-negative breast cancer. We found that TIL-Bs exhibit distinct antibody repertoire measures, including high clonal polarization and elevated somatic hypermutation rates, suggesting a local antigen-driven B-cell response. Importantly, TIL-Bs were highly mutated but non-class switched, suggesting that class-switch recombination may be inhibited in the TME. Tracing the distribution of TIL-B clones across various compartments indicated that they migrate to and from the TME. The data thus suggests that antibody repertoire signatures can serve as indicators for identifying tumor-reactive B cells.

Recapitulation of posttransfusion purpura by cross-strain platelet immunization in mice

Posttransfusion purpura (PTP) is an uncommon but life-threatening condition characterized by profound thrombocytopenia occurring ∼1 week after a blood transfusion. The hallmark of PTP is a potent immunoglobulin G antibody specific for a transfused platelet-specific alloantigen, usually located on glycoprotein IIb/IIIa (GPIIb/IIIa; αIIb/β3 integrin). It is widely thought that this alloantibody somehow causes the thrombocytopenia, despite absence from host platelets of the alloantigen for which it is specific. In studies described here, we found that cross-strain platelet immunization in mice commonly induces GPIIb/IIIa-specific alloantibodies combined with platelet-specific autoantibodies and varying degrees of thrombocytopenia, and we identified 1 strain combination (129S1Svlm/PWKPhJ) in which 95% of immunized mice made both types of antibody and developed severe thrombocytopenia. There was a strong inverse correlation between autoantibody strength and platelet decline (P < .0001) and plasma from mice that produced autoantibodies caused thrombocytopenia when transfused to syngeneic animals, arguing that autoantibodies were the cause of thrombocytopenia. The findings define a model in which a routine alloimmune response to platelets regularly transitions to an autoimmune reaction capable of causing severe thrombocytopenia and support the hypothesis that PTP is an autoimmune disorder.

A single donor is sufficient to produce a highly functional in vitro antibody library

Antibody complementarity determining region diversity has been considered to be the most important metric for the production of a functional antibody library. Generally, the greater the antibody library diversity, the greater the probability of selecting a diverse array of high affinity leads. According to this paradigm, the primary means of elevating library diversity has been by increasing the number of donors. In the present study we explored the possibility of creating an in vitro antibody library from a single healthy individual, showing that the number of lymphocytes, rather than the number of donors, is the key criterion in the production of a diverse and functional antibody library. We describe the construction of a high-quality phage display library comprising 5 × 10⁹ human antibodies by applying an efficient B cell extraction protocol from a single donor and a targeted V-gene amplification strategy favoring specific antibody families for their improved developability profiles. Each step of the library generation process was followed and validated by next generation sequencing to monitor the library quality and diversity. The functionality of the library was tested using several therapeutically relevant targets for which a vast number of different antibodies with desired biophysical properties were obtained.

PASA: Proteomic analysis of serum antibodies web server

MOTIVATION

A comprehensive characterization of the humoral response towards a specific antigen requires quantification of the B-cell receptor repertoire by next-generation sequencing (BCR-Seq), as well as the analysis of serum antibodies against this antigen, using proteomics. The proteomic analysis is challenging since it necessitates the mapping of antigen-specific peptides to individual B-cell clones.

RESULTS

The PASA web server provides a robust computational platform for the analysis and integration of data obtained from proteomics of serum antibodies. PASA maps peptides derived from antibodies raised against a specific antigen to corresponding antibody sequences. It then analyzes and integrates proteomics and BCR-Seq data, thus providing a comprehensive characterization of the humoral response. The PASA web server is freely available at https://pasa.tau.ac.il and open to all users without a login requirement.

FcRn Overexpression Expands Diversity of the Humoral Immune Response in bFcRn Transgenic Mice

The neonatal Fc receptor (FcRn) plays key roles in IgG and albumin homeostasis, maternal IgG transport, and antigen presentation of IgG-opsonized antigens. Previously, we reported that transgenic (Tg) mice that overexpress bovine FcRn (bFcRn) have augmented T-dependent humoral immune response with increased IgG protection, higher level of antigen-specific antibodies, greater number of antigen-specific B cells, and effective immune response even against weakly immunogenic epitopes. In this study we analyzed the diversity of the humoral immune response of bFcRn Tg mice, using a length distribution analysis (spectratyping) and next generation sequencing (NGS) of the immunoglobulin heavy chain variable regions. Our analysis showed that in response to immunization with ovalbumin or transfected cells that expressed a unique membrane protein, our Tg animals developed a more diverse plasma cell repertoire than controls, which manifested in greater numbers of different clones, and clusters with fewer highly expanded large clones, as identified by the variable region (CDR3) of the immunoglobulin heavy chain. The increased antibody diversity in Tg mice after immunization was observed at both IgM and IgG levels, indicating that the increased humoral immune diversity in Tg mice is due to a higher number of both activated, antigen-specific naïve and isotype switched B cells. We thus demonstrated that the BCR repertoire of the immunized bFcRn Tg animals is more diverse compared to wild type mice, which likely makes these Tg mice a better choice for monoclonal antibody production against challenging antigens, including the extracellular regions of cell membrane proteins.

Mucosal or systemic microbiota exposures shape the B cell repertoire

Colonization by the microbiota causes a marked stimulation of B cells and induction of immunoglobulin, but mammals colonized with many taxa have highly complex and individualized immunoglobulin repertoires^1,2. Here we use a simplified model of defined transient exposures to different microbial taxa in germ-free mice³ to deconstruct how the microbiota shapes the B cell pool and its functional responsiveness. We followed the development of the immunoglobulin repertoire in B cell populations, as well as single cells by deep sequencing. Microbial exposures at the intestinal mucosa generated oligoclonal responses that differed from those of germ-free mice, and from the diverse repertoire that was generated after intravenous systemic exposure to microbiota. The IgA repertoire-predominantly to cell-surface antigens-did not expand after dose escalation, whereas increased systemic exposure broadened the IgG repertoire to both microbial cytoplasmic and cell-surface antigens. These microbial exposures induced characteristic immunoglobulin heavy-chain repertoires in B cells, mainly at memory and plasma cell stages. Whereas sequential systemic exposure to different microbial taxa diversified the IgG repertoire and facilitated alternative specific responses, sequential mucosal exposure produced limited overlapping repertoires and the attrition of initial IgA binding specificities. This shows a contrast between a flexible response to systemic exposure with the need to avoid fatal sepsis, and a restricted response to mucosal exposure that reflects the generic nature of host-microbial mutualism in the mucosa.

Validation of Methods to Assess the Immunoglobulin Gene Repertoire in Tissues Obtained from Mice on the International Space Station

Spaceflight is known to affect immune cell populations. In particular, splenic B-cell numbers decrease during spaceflight and in ground-based physiological models. Although antibody isotype changes have been assessed during and after spaceflight, an extensive characterization of the impact of spaceflight on antibody composition has not been conducted in mice. Next Generation Sequencing and bioinformatic tools are now available to assess antibody repertoires. We can now identify immunoglobulin gene-segment usage, junctional regions, and modifications that contribute to specificity and diversity. Due to limitations on the International Space Station, alternate sample collection and storage methods must be employed. Our group compared Illumina MiSeq® sequencing data from multiple sample preparation methods in normal C57Bl/6J mice to validate that sample preparation and storage would not bias the outcome of antibody repertoire characterization. In this report, we also compared sequencing techniques and a bioinformatic workflow on the data output when we assessed the IgH and Igκ variable gene usage. Our bioinformatic workflow has been optimized for Illumina HiSeq® and MiSeq® datasets, and is designed specifically to reduce bias, capture the most information from Ig sequences, and produce a data set that provides other data mining options.

Uncoupling of invasive bacterial mucosal immunogenicity from pathogenicity

There is the notion that infection with a virulent intestinal pathogen induces generally stronger mucosal adaptive immunity than the exposure to an avirulent strain. Whether the associated mucosal inflammation is important or redundant for effective induction of immunity is, however, still unclear. Here we use a model of auxotrophic Salmonella infection in germ-free mice to show that live bacterial virulence factor-driven immunogenicity can be uncoupled from inflammatory pathogenicity. Although live auxotrophic Salmonella no longer causes inflammation, its mucosal virulence factors remain the main drivers of protective mucosal immunity; virulence factor-deficient, like killed, bacteria show reduced efficacy. Assessing the involvement of innate pathogen sensing mechanisms, we show MYD88/TRIF, Caspase-1/Caspase-11 inflammasome, and NOD1/NOD2 nodosome signaling to be individually redundant. In colonized animals we show that microbiota metabolite cross-feeding may recover intestinal luminal colonization but not pathogenicity. Consequent immunoglobulin A immunity and microbial niche competition synergistically protect against Salmonella wild-type infection.

Molecular T-Cell Repertoire Analysis as Source of Prognostic and Predictive Biomarkers for Checkpoint Blockade Immunotherapy

The T cells are key players of the response to checkpoint blockade immunotherapy (CBI) and monitoring the strength and specificity of antitumor T-cell reactivity remains a crucial but elusive component of precision immunotherapy. The entire assembly of T-cell receptor (TCR) sequences accounts for antigen specificity and strength of the T-cell immune response. The TCR repertoire hence represents a “footprint” of the conditions faced by T cells that dynamically evolves according to the challenges that arise for the immune system, such as tumor neo-antigenic load. Hence, TCR repertoire analysis is becoming increasingly important to comprehensively understand the nature of a successful antitumor T-cell response, and to improve the success and safety of current CBI.

Molecular Landscape of Anti-Drug Antibodies Reveals the Mechanism of the Immune Response Following Treatment With TNFα Antagonists

Drugs formulated from monoclonal antibodies (mAbs) are clinically effective in various diseases. Repeated administration of mAbs, however, elicits an immune response in the form of anti-drug-antibodies (ADA), thereby reducing the drug's efficacy. Notwithstanding their importance, the molecular landscape of ADA and the mechanisms involved in their formation are not fully understood. Using a newly developed quantitative bio-immunoassay, we found that ADA concentrations specific to TNFα antagonists can exceed extreme concentrations of 1 mg/ml with a wide range of neutralization capacity. Our data further suggest a preferential use of the λ light chain in a subset of neutralizing ADA. Moreover, we show that administration of TNFα antagonists result in a vaccine-like response whereby ADA formation is governed by the extrafollicular T cell-independent immune response. Our bio-immunoassay coupled with insights on the nature of the immune response can be leveraged to improve mAb immunogenicity assessment and facilitate improvement in therapeutic intervention strategies.

Resident memory CD8 T cells persist for years in human small intestine

Bartolomé-Casado et al. demonstrate that human gut contains large numbers of resident memory CD8 T cells that survive for years. Intestinal CD8 Trm cells have a clonally expanded immune repertoire that is stable over time and exhibit enhanced protective capabilities.

Resident memory CD8 T (Trm) cells have been shown to provide effective protective responses in the small intestine (SI) in mice. A better understanding of the generation and persistence of SI CD8 Trm cells in humans may have implications for intestinal immune-mediated diseases and vaccine development. Analyzing normal and transplanted human SI, we demonstrated that the majority of SI CD8 T cells were bona fide CD8 Trm cells that survived for >1 yr in the graft. Intraepithelial and lamina propria CD8 Trm cells showed a high clonal overlap and a repertoire dominated by expanded clones, conserved both spatially in the intestine and over time. Functionally, lamina propria CD8 Trm cells were potent cytokine producers, exhibiting a polyfunctional (IFN-γ⁺ IL-2⁺ TNF-α⁺) profile, and efficiently expressed cytotoxic mediators after stimulation. These results suggest that SI CD8 Trm cells could be relevant targets for future oral vaccines and therapeutic strategies for gut disorders.

Epithelial endoplasmic reticulum stress orchestrates a protective IgA response

Immunoglobulin A (IgA) is the major secretory immunoglobulin isotype found at mucosal surfaces, where it regulates microbial commensalism and excludes luminal factors from contacting intestinal epithelial cells (IECs). IgA is induced by both T cell-dependent and -independent (TI) pathways. However, little is known about TI regulation. We report that IEC endoplasmic reticulum (ER) stress induces a polyreactive IgA response, which is protective against enteric inflammation. IEC ER stress causes TI and microbiota-independent expansion and activation of peritoneal B1b cells, which culminates in increased lamina propria and luminal IgA. Increased numbers of IgA-producing plasma cells were observed in healthy humans with defective autophagy, who are known to exhibit IEC ER stress. Upon ER stress, IECs communicate signals to the peritoneum that induce a barrier-protective TI IgA response.

Comparative Analysis of B-Cell Receptor Repertoires Induced by Live Yellow Fever Vaccine in Young and Middle-Age Donors

Age-related changes can significantly alter the state of adaptive immune system and often lead to attenuated response to novel pathogens and vaccination. In present study we employed 5′RACE UMI-based full length and nearly error-free immunoglobulin profiling to compare plasma cell antibody repertoires in young (19–26 years) and middle-age (45–58 years) individuals vaccinated with a live yellow fever vaccine, modeling a newly encountered pathogen. Our analysis has revealed age-related differences in the responding antibody repertoire ranging from distinct IGH CDR3 repertoire properties to differences in somatic hypermutation intensity and efficiency and antibody lineage tree structure. Overall, our findings suggest that younger individuals respond with a more diverse antibody repertoire and employ a more efficient somatic hypermutation process than elder individuals in response to a newly encountered pathogen.

Observed Antibody Space: A Resource for Data Mining Next-Generation Sequencing of Antibody Repertoires

Abs are immune system proteins that recognize noxious molecules for elimination. Their sequence diversity and binding versatility have made Abs the primary class of biopharmaceuticals. Recently, it has become possible to query their immense natural diversity using next-generation sequencing of Ig gene repertoires (Ig-seq). However, Ig-seq outputs are currently fragmented across repositories and tend to be presented as raw nucleotide reads, which means nontrivial effort is required to reuse the data for analysis. To address this issue, we have collected Ig-seq outputs from 55 studies, covering more than half a billion Ab sequences across diverse immune states, organisms (primarily human and mouse), and individuals. We have sorted, cleaned, annotated, translated, and numbered these sequences and make the data available via our Observed Antibody Space (OAS) resource at http://antibodymap.org The data within OAS will be regularly updated with newly released Ig-seq datasets. We believe OAS will facilitate data mining of immune repertoires for improved understanding of the immune system and development of better biotherapeutics.

Comparison of methods for phylogenetic B-cell lineage inference using time-resolved antibody repertoire simulations (AbSim)

Motivation

The evolution of antibody repertoires represents a hallmark feature of adaptive B-cell immunity. Recent advancements in high-throughput sequencing have dramatically increased the resolution to which we can measure the molecular diversity of antibody repertoires, thereby offering for the first time the possibility to capture the antigen-driven evolution of B cells. However, there does not exist a repertoire simulation framework yet that enables the comparison of commonly utilized phylogenetic methods with regard to their accuracy in inferring antibody evolution.

Results

Here, we developed AbSim, a time-resolved antibody repertoire simulation framework, which we exploited for testing the accuracy of methods for the phylogenetic reconstruction of B-cell lineages and antibody molecular evolution. AbSim enables the (i) simulation of intermediate stages of antibody sequence evolution and (ii) the modeling of immunologically relevant parameters such as duration of repertoire evolution, and the method and frequency of mutations. First, we validated that our repertoire simulation framework recreates replicates topological similarities observed in experimental sequencing data. Second, we leveraged Absim to show that current methods fail to a certain extent to predict the true phylogenetic tree correctly. Finally, we formulated simulation-validated guidelines for antibody evolution, which in the future will enable the development of accurate phylogenetic methods.

Availability and implementation

https://cran.r-project.org/web/packages/AbSim/index.html.

Contact

sai.reddy@ethz.ch.

Supplementary information

Supplementary data are available at Bioinformatics online.

ASAP - A Webserver for Immunoglobulin-Sequencing Analysis Pipeline

Reproducible and robust data on antibody repertoires are invaluable for basic and applied immunology. Next-generation sequencing (NGS) of antibody variable regions has emerged as a powerful tool in systems immunology, providing quantitative molecular information on antibody polyclonal composition. However, major computational challenges exist when analyzing antibody sequences, from error handling to hypermutation profiles and clonal expansion analyses. In this work, we developed the ASAP (A webserver for Immunoglobulin-Seq Analysis Pipeline) webserver (https://asap.tau.ac.il). The input to ASAP is a paired-end sequence dataset from one or more replicates, with or without unique molecular identifiers. These datasets can be derived from NGS of human or murine antibody variable regions. ASAP first filters and annotates the sequence reads using public or user-provided germline sequence information. The ASAP webserver next performs various calculations, including somatic hypermutation level, CDR3 lengths, V(D)J family assignments, and V(D)J combination distribution. These analyses are repeated for each replicate. ASAP provides additional information by analyzing the commonalities and differences between the repeats (“joint” analysis). For example, ASAP examines the shared variable regions and their frequency in each replicate to determine which sequences are less likely to be a result of a sample preparation derived and/or sequencing errors. Moreover, ASAP clusters the data to clones and reports the identity and prevalence of top ranking clones (clonal expansion analysis). ASAP further provides the distribution of synonymous and non-synonymous mutations within the V genes somatic hypermutations. Finally, ASAP provides means to process the data for proteomic analysis of serum/secreted antibodies by generating a variable region database for liquid chromatography high resolution tandem mass spectrometry (LC-MS/MS) interpretation. ASAP is user-friendly, free, and open to all users, with no login requirement. ASAP is applicable for researchers interested in basic questions related to B cell development and differentiation, as well as applied researchers who are interested in vaccine development and monoclonal antibody engineering. By virtue of its user-friendliness, ASAP opens the antibody analysis field to non-expert users who seek to boost their research with immune repertoire analysis.

Synthetic Standards Combined With Error and Bias Correction Improve the Accuracy and Quantitative Resolution of Antibody Repertoire Sequencing in Human Naïve and Memory B Cells

High-throughput sequencing of immunoglobulin (Ig) repertoires (Ig-seq) is a powerful method for quantitatively interrogating B cell receptor sequence diversity. When applied to human repertoires, Ig-seq provides insight into fundamental immunological questions, and can be implemented in diagnostic and drug discovery projects. However, a major challenge in Ig-seq is ensuring accuracy, as library preparation protocols and sequencing platforms can introduce substantial errors and bias that compromise immunological interpretation. Here, we have established an approach for performing highly accurate human Ig-seq by combining synthetic standards with a comprehensive error and bias correction pipeline. First, we designed a set of 85 synthetic antibody heavy-chain standards (in vitro transcribed RNA) to assess correction workflow fidelity. Next, we adapted a library preparation protocol that incorporates unique molecular identifiers (UIDs) for error and bias correction which, when applied to the synthetic standards, resulted in highly accurate data. Finally, we performed Ig-seq on purified human circulating B cell subsets (naïve and memory), combined with a cellular replicate sampling strategy. This strategy enabled robust and reliable estimation of key repertoire features such as clonotype diversity, germline segment, and isotype subclass usage, and somatic hypermutation. We anticipate that our standards and error and bias correction pipeline will become a valuable tool for researchers to validate and improve accuracy in human Ig-seq studies, thus leading to potentially new insights and applications in human antibody repertoire profiling.

Microbial symbionts regulate the primary Ig repertoire

Symbiotic relationships help shape immune fitness. Chen et al. demonstrate that microbial symbionts influence host immunity by enriching frequencies of antibacterial specificities within the naive B cell receptor repertoire and that this may have consequences for mucosal and systemic immunity.

The ability of immunoglobulin (Ig) to recognize pathogens is critical for optimal immune fitness. Early events that shape preimmune Ig repertoires, expressed on IgM⁺ IgD⁺ B cells as B cell receptors (BCRs), are poorly defined. Here, we studied germ-free mice and conventionalized littermates to explore the hypothesis that symbiotic microbes help shape the preimmune Ig repertoire. Ig-binding assays showed that exposure to conventional microbial symbionts enriched frequencies of antibacterial IgM⁺ IgD⁺ B cells in intestine and spleen. This enrichment affected follicular B cells, involving a diverse set of Ig-variable region gene segments, and was T cell–independent. Functionally, enrichment of microbe reactivity primed basal levels of small intestinal T cell–independent, symbiont-reactive IgA and enhanced systemic IgG responses to bacterial immunization. These results demonstrate that microbial symbionts influence host immunity by enriching frequencies of antibacterial specificities within preimmune B cell repertoires and that this may have consequences for mucosal and systemic immunity.

Many Routes to an Antibody Heavy-Chain CDR3: Necessary, Yet Insufficient, for Specific Binding

Because of its great potential for diversity, the immunoglobulin heavy-chain complementarity-determining region 3 (HCDR3) is taken as an antibody molecule’s most important component in conferring binding activity and specificity. For this reason, HCDR3s have been used as unique identifiers to investigate adaptive immune responses in vivo and to characterize in vitro selection outputs where display systems were employed. Here, we show that many different HCDR3s can be identified within a target-specific antibody population after in vitro selection. For each identified HCDR3, a number of different antibodies bearing differences elsewhere can be found. In such selected populations, all antibodies with the same HCDR3 recognize the target, albeit at different affinities. In contrast, within unselected populations, the majority of antibodies with the same HCDR3 sequence do not bind the target. In one HCDR3 examined in depth, all target-specific antibodies were derived from the same VDJ rearrangement, while non-binding antibodies with the same HCDR3 were derived from many different V and D gene rearrangements. Careful examination of previously published in vivo datasets reveals that HCDR3s shared between, and within, different individuals can also originate from rearrangements of different V and D genes, with up to 26 different rearrangements yielding the same identical HCDR3 sequence. On the basis of these observations, we conclude that the same HCDR3 can be generated by many different rearrangements, but that specific target binding is an outcome of unique rearrangements and VL pairing: the HCDR3 is necessary, albeit insufficient, for specific antibody binding.

Characterization of the naive murine antibody repertoire using unamplified high-throughput sequencing

Antibody specificity and diversity are generated through the enzymatic splicing of genomic gene segments within each B cell. Antibodies are heterodimers of heavy- and light-chains encoded on separate loci. We studied the antibody repertoire from pooled, splenic tissue of unimmunized, adult female C57BL/6J mice, using high-throughput sequencing (HTS) without amplification of antibody transcripts. We recovered over 90,000 heavy-chain and over 135,000 light-chain immunoglobulin sequences. Individual V-, D-, and J-gene segment usage was uniform among the three mouse pools, particularly in highly abundant gene segments, with low frequency V-gene segments not being detected in all pools. Despite the similar usage of individual gene segments, the repertoire of individual B-cell CDR3 amino acid sequences in each mouse pool was highly varied, affirming the combinatorial diversity in the B-cell pool that has been previously demonstrated. There also was some skewing in the V-gene segments that were detected depending on chromosomal location. This study presents a unique, non-primer biased glimpse of the conventionally housed, unimmunized antibody repertoire of the C57BL6/J mouse.

VDJPipe: a pipelined tool for pre-processing immune repertoire sequencing data

Background

Pre-processing of high-throughput sequencing data for immune repertoire profiling is essential to insure high quality input for downstream analysis. VDJPipe is a flexible, high-performance tool that can perform multiple pre-processing tasks with just a single pass over the data files.

Results

Processing tasks provided by VDJPipe include base composition statistics calculation, read quality statistics calculation, quality filtering, homopolymer filtering, length and nucleotide filtering, paired-read merging, barcode demultiplexing, 5′ and 3′ PCR primer matching, and duplicate reads collapsing. VDJPipe utilizes a pipeline approach whereby multiple processing steps are performed in a sequential workflow, with the output of each step passed as input to the next step automatically. The workflow is flexible enough to handle the complex barcoding schemes used in many immunosequencing experiments. Because VDJPipe is designed for computational efficiency, we evaluated this by comparing execution times with those of pRESTO, a widely-used pre-processing tool for immune repertoire sequencing data. We found that VDJPipe requires <10% of the run time required by pRESTO.

Conclusions

VDJPipe is a high-performance tool that is optimized for pre-processing large immune repertoire sequencing data sets.

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.

Validation of Methods to Assess the Immunoglobulin Gene Repertoire in Tissues Obtained from Mice on the International Space Station

Spaceflight is known to affect immune cell populations. In particular, splenic B cell numbers decrease during spaceflight and in ground-based physiological models. Although antibody isotype changes have been assessed during and after space flight, an extensive characterization of the impact of spaceflight on antibody composition has not been conducted in mice. Next Generation Sequencing and bioinformatic tools are now available to assess antibody repertoires. We can now identify immunoglobulin gene- segment usage, junctional regions, and modifications that contribute to specificity and diversity. Due to limitations on the International Space Station, alternate sample collection and storage methods must be employed. Our group compared Illumina MiSeq sequencing data from multiple sample preparation methods in normal C57Bl/6J mice to validate that sample preparation and storage would not bias the outcome of antibody repertoire characterization. In this report, we also compared sequencing techniques and a bioinformatic workflow on the data output when we assessed the IgH and Igκ variable gene usage. This included assessments of our bioinformatic workflow on Illumina HiSeq and MiSeq datasets and is specifically designed to reduce bias, capture the most information from Ig sequences, and produce a data set that provides other data mining options. We validated our workflow by comparing our normal mouse MiSeq data to existing murine antibody repertoire studies validating it for future antibody repertoire studies.

Assessment of antibody library diversity through next generation sequencing and technical error compensation

Antibody libraries are important resources to derive antibodies to be used for a wide range of applications, from structural and functional studies to intracellular protein interference studies to developing new diagnostics and therapeutics. Whatever the goal, the key parameter for an antibody library is its complexity (also known as diversity), i.e. the number of distinct elements in the collection, which directly reflects the probability of finding in the library an antibody against a given antigen, of sufficiently high affinity. Quantitative evaluation of antibody library complexity and quality has been for a long time inadequately addressed, due to the high similarity and length of the sequences of the library. Complexity was usually inferred by the transformation efficiency and tested either by fingerprinting and/or sequencing of a few hundred random library elements. Inferring complexity from such a small sampling is, however, very rudimental and gives limited information about the real diversity, because complexity does not scale linearly with sample size. Next-generation sequencing (NGS) has opened new ways to tackle the antibody library complexity quality assessment. However, much remains to be done to fully exploit the potential of NGS for the quantitative analysis of antibody repertoires and to overcome current limitations. To obtain a more reliable antibody library complexity estimate here we show a new, PCR-free, NGS approach to sequence antibody libraries on Illumina platform, coupled to a new bioinformatic analysis and software (Diversity Estimator of Antibody Library, DEAL) that allows to reliably estimate the complexity, taking in consideration the sequencing error.

High-Throughput Immunogenetics for Clinical and Research Applications in Immunohematology: Potential and Challenges

Analysis and interpretation of Ig and TCR gene rearrangements in the conventional, low-throughput way have their limitations in terms of resolution, coverage, and biases. With the advent of high-throughput, next-generation sequencing (NGS) technologies, a deeper analysis of Ig and/or TCR (IG/TR) gene rearrangements is now within reach, which impacts on all main applications of IG/TR immunogenetic analysis. To bridge the generation gap from low- to high-throughput analysis, the EuroClonality-NGS Consortium has been formed, with the main objectives to develop, standardize, and validate the entire workflow of IG/TR NGS assays for 1) clonality assessment, 2) minimal residual disease detection, and 3) repertoire analysis. This concerns the preanalytical (sample preparation, target choice), analytical (amplification, NGS), and postanalytical (immunoinformatics) phases. Here we critically discuss pitfalls and challenges of IG/TR NGS methodology and its applications in hemato-oncology and immunology.

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.

Deep sequencing in library selection projects: what insight does it bring?

High throughput sequencing is poised to change all aspects of the way antibodies and other binders are discovered and engineered. Millions of available sequence reads provide an unprecedented sampling depth able to guide the design and construction of effective, high quality naïve libraries containing tens of billions of unique molecules. Furthermore, during selections, high throughput sequencing enables quantitative tracing of enriched clones and position-specific guidance to amino acid variation under positive selection during antibody engineering. Successful application of the technologies relies on specific PCR reagent design, correct sequencing platform selection, and effective use of computational tools and statistical measures to remove error, identify antibodies, estimate diversity, and extract signatures of selection from the clone down to individual structural positions. Here we review these considerations and discuss some of the remaining challenges to the widespread adoption of the technology.

Eye on the B-ALL: B-cell receptor repertoires reveal persistence of numerous B-lymphoblastic leukemia subclones from diagnosis to relapse

The strongest predictor of relapse in B-cell acute lymphoblastic leukemia (B-ALL) is the level of persistence of tumor cells after initial therapy. The high mutation rate of the B-cell receptor (BCR) locus allows high-resolution tracking of the architecture, evolution and clonal dynamics of B-ALL. Using longitudinal BCR repertoire sequencing, we find that the BCR undergoes an unexpectedly high level of clonal diversification in B-ALL cells through both somatic hypermutation and secondary rearrangements, which can be used for tracking the subclonal composition of the disease and detect minimal residual disease with unprecedented sensitivity. We go on to investigate clonal dynamics of B-ALL using BCR phylogenetic analyses of paired diagnosis-relapse samples and find that large numbers of small leukemic subclones present at diagnosis re-emerge at relapse alongside a dominant clone. Our findings suggest that in all informative relapsed patients, the survival of large numbers of clonogenic cells beyond initial chemotherapy is a surrogate for inherent partial chemoresistance or inadequate therapy, providing an increased opportunity for subsequent emergence of fully resistant clones. These results frame early cytoreduction as an important determinant of long-term outcome.

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

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

The mouse antibody heavy chain repertoire is germline-focused and highly variable between inbred strains

The human and mouse antibody repertoires are formed by identical processes, but like all small animals, mice only have sufficient lymphocytes to express a small part of the potential antibody repertoire. In this study, we determined how the heavy chain repertoires of two mouse strains are generated. Analysis of IgM- and IgG-associated VDJ rearrangements generated by high-throughput sequencing confirmed the presence of 99 functional immunoglobulin heavy chain variable (IGHV) genes in the C57BL/6 genome, and inferred the presence of 164 IGHV genes in the BALB/c genome. Remarkably, only five IGHV sequences were common to both strains. Compared with humans, little N nucleotide addition was seen in the junctions of mouse VDJ genes. Germline human IgG-associated IGHV genes are rare, but many murine IgG-associated IGHV genes were unmutated. Together these results suggest that the expressed mouse repertoire is more germline-focused than the human repertoire. The apparently divergent germline repertoires of the mouse strains are discussed with reference to reports that inbred mouse strains carry blocks of genes derived from each of the three subspecies of the house mouse. We hypothesize that the germline genes of BALB/c and C57BL/6 mice may originally have evolved to generate distinct germline-focused antibody repertoires in the different mouse subspecies.

IgRepertoireConstructor: a novel algorithm for antibody repertoire construction and immunoproteogenomics analysis

The analysis of concentrations of circulating antibodies in serum (antibody repertoire) is a fundamental, yet poorly studied, problem in immunoinformatics. The two current approaches to the analysis of antibody repertoires [next generation sequencing (NGS) and mass spectrometry (MS)] present difficult computational challenges since antibodies are not directly encoded in the germline but are extensively diversified by somatic recombination and hypermutations. Therefore, the protein database required for the interpretation of spectra from circulating antibodies is custom for each individual. Although such a database can be constructed via NGS, the reads generated by NGS are error-prone and even a single nucleotide error precludes identification of a peptide by the standard proteomics tools. Here, we present the IgRepertoireConstructor algorithm that performs error-correction of immunosequencing reads and uses mass spectra to validate the constructed antibody repertoires.

AVAILABILITY AND IMPLEMENTATION

IgRepertoireConstructor is open source and freely available as a C++ and Python program running on all Unix-compatible platforms. The source code is available from http://bioinf.spbau.ru/igtools.

CONTACT

ppevzner@ucsd.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

In-Depth Assessment of Within-Individual and Inter-Individual Variation in the B Cell Receptor Repertoire

High-throughput sequencing of the B cell receptor (BCR) repertoire can provide rapid characterization of the B cell response in a wide variety of applications in health, after vaccination and in infectious, inflammatory and immune-driven disease, and is starting to yield clinical applications. However, the interpretation of repertoire data is compromised by a lack of studies to assess the intra and inter-individual variation in the BCR repertoire over time in healthy individuals. We applied a standardized isotype-specific BCR repertoire deep sequencing protocol to a single highly sampled participant, and then evaluated the method in 9 further participants to comprehensively describe such variation. We assessed total repertoire metrics of mutation, diversity, VJ gene usage and isotype subclass usage as well as tracking specific BCR sequence clusters. There was good assay reproducibility (both in PCR amplification and biological replicates), but we detected striking fluctuations in the repertoire over time that we hypothesize may be due to subclinical immune activation. Repertoire properties were unique for each individual, which could partly be explained by a decrease in IgG2 with age, and genetic differences at the immunoglobulin locus. There was a small repertoire of public clusters (0.5, 0.3, and 1.4% of total IgA, IgG, and IgM clusters, respectively), which was enriched for expanded clusters containing sequences with suspected specificity toward antigens that should have been historically encountered by all participants through prior immunization or infection. We thus provide baseline BCR repertoire information that can be used to inform future study design, and aid in interpretation of results from these studies. Furthermore, our results indicate that BCR repertoire studies could be used to track changes in the public repertoire in and between populations that might relate to population immunity against infectious diseases, and identify the characteristics of inflammatory and immunological diseases.

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.