Reproducibility and Reuse of Adaptive Immune Receptor Repertoire Data

Vaccination induces broadly neutralizing antibody precursors to HIV gp41

A key barrier to the development of vaccines that induce broadly neutralizing antibodies (bnAbs) against human immunodeficiency virus (HIV) and other viruses of high antigenic diversity is the design of priming immunogens that induce rare bnAb-precursor B cells. The high neutralization breadth of the HIV bnAb 10E8 makes elicitation of 10E8-class bnAbs desirable; however, the recessed epitope within gp41 makes envelope trimers poor priming immunogens and requires that 10E8-class bnAbs possess a long heavy chain complementarity determining region 3 (HCDR3) with a specific binding motif. We developed germline-targeting epitope scaffolds with affinity for 10E8-class precursors and engineered nanoparticles for multivalent display. Scaffolds exhibited epitope structural mimicry and bound bnAb-precursor human naive B cells in ex vivo screens, protein nanoparticles induced bnAb-precursor responses in stringent mouse models and rhesus macaques, and mRNA-encoded nanoparticles triggered similar responses in mice. Thus, germline-targeting epitope scaffold nanoparticles can elicit rare bnAb-precursor B cells with predefined binding specificities and HCDR3 features.

Best practices for machine learning in antibody discovery and development

In the past 40 years, therapeutic antibody discovery and development have advanced considerably, with machine learning (ML) offering a promising way to speed up the process by reducing costs and the number of experiments required. Recent progress in ML-guided antibody design and development (D&D) has been hindered by the diversity of data sets and evaluation methods, which makes it difficult to conduct comparisons and assess utility. Establishing standards and guidelines will be crucial for the wider adoption of ML and the advancement of the field. This perspective critically reviews current practices, highlights common pitfalls and proposes method development and evaluation guidelines for various ML-based techniques in therapeutic antibody D&D. Addressing challenges across the ML process, best practices are recommended for each stage to enhance reproducibility and progress.

Neem leaf glycoprotein binding to Dectin-1 receptors on dendritic cell induces type-1 immunity through CARD9 mediated intracellular signal to NFκB

BACKGROUND

A water-soluble ingredient of mature leaves of the tropical mahogany 'Neem' (Azadirachta indica), was identified as glycoprotein, thus being named as 'Neem Leaf Glycoprotein' (NLGP). This non-toxic leaf-component regressed cancerous murine tumors (melanoma, carcinoma, sarcoma) recurrently in different experimental circumstances by boosting prime antitumor immune attributes. Such antitumor immunomodulation, aid cytotoxic T cell (Tc)-based annihilation of tumor cells. This study focused on identifying and characterizing the signaling gateway that initiate this systemic immunomodulation. In search of this gateway, antigen-presenting cells (APCs) were explored, which activate and induce the cytotoxic thrust in Tc cells.

METHODS

Six glycoprotein-binding C-type lectins found on APCs, namely, MBR, Dectin-1, Dectin-2, DC-SIGN, DEC205 and DNGR-1 were screened on bone marrow-derived dendritic cells from C57BL/6 J mice. Fluorescence microscopy, RT-PCR, flow cytometry and ELISA revealed Dectin-1 as the NLGP-binding receptor, followed by verifications through RNAi. Following detection of β-Glucans in NLGP, their interactions with Dectin-1 were explored in silico. Roles of second messengers and transcription factors in the downstream signal were studied by co-immunoprecipitation, western blotting, and chromatin-immunoprecipitation. Intracellularization of FITC-coupled NLGP was observed by processing confocal micrographs of DCs.

RESULTS

Considering extents of hindrance in NLGP-driven transcription rates of the cytokines IL-10 and IL-12p35 by receptor-neutralization, Dectin-1 receptors on dendritic cells were found to bind NLGP through the ligand's peripheral β-Glucan chains. The resulting signal phosphorylates PKCδ, forming a trimolecular complex of CARD9, Bcl10 and MALT1, which in turn activates the canonical NFκB-pathway of transcription-regulation. Consequently, the NFκB-heterodimer p65:p50 enhances Il12a transcription and the p50:p50 homodimer represses Il10 transcription, bringing about a cytokine-based systemic-bias towards type-1 immune environment. Further, NLGP gets engulfed within dendritic cells, possibly through endocytic activities of Dectin-1.

CONCLUSION

NLGP's binding to Dectin-1 receptors on murine dendritic cells, followed by the intracellular signal, lead to NFκB-mediated contrasting regulation of cytokine-transcriptions, initiating a pro-inflammatory immunopolarization, which amplifies further by the responding immune cells including Tc cells, alongside their enhanced cytotoxicity. These insights into the initiation of mammalian systemic immunomodulation by NLGP at cellular and molecular levels, may help uncovering its mode of action as a novel immunomodulator against human cancers, following clinical trials.

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.

Human immunoglobulin gene allelic variation impacts germline-targeting vaccine priming

Vaccine priming immunogens that activate germline precursors for broadly neutralizing antibodies (bnAbs) have promise for development of precision vaccines against major human pathogens. In a clinical trial of the eOD-GT8 60mer germline-targeting immunogen, higher frequencies of vaccine-induced VRC01-class bnAb-precursor B cells were observed in the high dose compared to the low dose group. Through immunoglobulin heavy chain variable (IGHV) genotyping, statistical modeling, quantification of IGHV1-2 allele usage and B cell frequencies in the naive repertoire for each trial participant, and antibody affinity analyses, we found that the difference between dose groups in VRC01-class response frequency was best explained by IGHV1-2 genotype rather than dose and was most likely due to differences in IGHV1-2 B cell frequencies for different genotypes. The results demonstrate the need to define population-level immunoglobulin allelic variations when designing germline-targeting immunogens and evaluating them in clinical trials.

The clinical applications of immunosequencing

Technological advances in high-throughput sequencing have opened the door for the interrogation of adaptive immune responses at unprecedented scale. It is now possible to determine the sequences of antibodies or T-cell receptors produced by individual B and T cells in a sample. This capability, termed immunosequencing, has transformed the study of both infectious and non-infectious diseases by allowing the tracking of dynamic changes in B and T cell clonal populations over time. This has improved our understanding of the pathology of cancers, autoimmune diseases, and infectious diseases. However, to date there has been only limited clinical adoption of the technology. Advances over the last decade and on the horizon that reduce costs and improve interpretability could enable widespread clinical use. Many clinical applications have been proposed and, while most are still undergoing research and development, some methods relying on immunosequencing data have been implemented, the most widespread of which is the detection of measurable residual disease. Here, we review the diagnostic, prognostic, and therapeutic applications of immunosequencing for both infectious and non-infectious diseases.

Dandelion uses the single-cell adaptive immune receptor repertoire to explore lymphocyte developmental origins

Assessment of single-cell gene expression (single-cell RNA sequencing) and adaptive immune receptor (AIR) sequencing (scVDJ-seq) has been invaluable in studying lymphocyte biology. Here we introduce Dandelion, a computational pipeline for scVDJ-seq analysis. It enables the application of standard V(D)J analysis workflows to single-cell datasets, delivering improved V(D)J contig annotation and the identification of nonproductive and partially spliced contigs. We devised a strategy to create an AIR feature space that can be used for both differential V(D)J usage analysis and pseudotime trajectory inference. The application of Dandelion improved the alignment of human thymic development trajectories of double-positive T cells to mature single-positive CD4/CD8 T cells, generating predictions of factors regulating lineage commitment. Dandelion analysis of other cell compartments provided insights into the origins of human B1 cells and ILC/NK cell development, illustrating the power of our approach. Dandelion is available at https://www.github.com/zktuong/dandelion .

Potential of TCR sequencing in graft-versus-host disease

Graft-versus-host disease (GvHD) remains one of the major complications following allogeneic haematopoietic stem cell transplantation (allo-HSCT). GvHD can occur in almost every tissue, with the skin, liver, and intestines being the mainly affected organs. T cells are implicated in initiating GvHD. T cells identify a broad range of antigens and mediate the immune response through receptors on their surfaces (T cell receptors, TCRs). The composition of TCRs within a T cell population defines the TCR repertoire of an individual, and this repertoire represents exposure to self and non-self proteins. Monitoring the changes in the TCR repertoire using TCR sequencing can provide an indication of the dynamics of a T cell population. Monitoring the frequency and specificities of specific TCR clonotypes longitudinally in different conditions and specimens (peripheral blood, GvHD-affected tissue samples) can provide insights into factors modulating immune reactions following allogeneic transplantation and will help to understand the underlying mechanisms mediating GvHD. This review provides insights into current studies of the TCR repertoire in GvHD and potential future clinical implications of TCR sequencing.

Vaccination induces HIV broadly neutralizing antibody precursors in humans

Broadly neutralizing antibodies (bnAbs) can protect against HIV infection but have not been induced by human vaccination. A key barrier to bnAb induction is vaccine priming of rare bnAb-precursor B cells. In a randomized, double-blind, placebo-controlled phase 1 clinical trial, the HIV vaccine-priming candidate eOD-GT8 60mer adjuvanted with AS01B had a favorable safety profile and induced VRC01-class bnAb precursors in 97% of vaccine recipients with median frequencies reaching 0.1% among immunoglobulin G B cells in blood. bnAb precursors shared properties with bnAbs and gained somatic hypermutation and affinity with the boost. The results establish clinical proof of concept for germline-targeting vaccine priming, support development of boosting regimens to induce bnAbs, and encourage application of the germline-targeting strategy to other targets in HIV and other pathogens.

AIRRscape: An interactive tool for exploring B-cell receptor repertoires and antibody responses

The sequencing of antibody repertoires of B-cells at increasing coverage and depth has led to the identification of vast numbers of immunoglobulin heavy and light chains. However, the size and complexity of these Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) datasets makes it difficult to perform exploratory analyses. To aid in data exploration, we have developed AIRRscape, an R Shiny-based interactive web browser application that enables B-cell receptor (BCR) and antibody feature discovery through comparisons among multiple repertoires. Using AIRR-seq data as input, AIRRscape starts by aggregating and sorting repertoires into interactive and explorable bins of germline V-gene, germline J-gene, and CDR3 length, providing a high-level view of the entire repertoire. Interesting subsets of repertoires can be quickly identified and selected, and then network topologies of CDR3 motifs can be generated for further exploration. Here we demonstrate AIRRscape using patient BCR repertoires and sequences of published monoclonal antibodies to investigate patterns of humoral immunity to three viral pathogens: SARS-CoV-2, HIV-1, and DENV (dengue virus). AIRRscape reveals convergent antibody sequences among datasets for all three pathogens, although HIV-1 antibody datasets display limited convergence and idiosyncratic responses. We have made AIRRscape available as a web-based Shiny application, along with code on GitHub to encourage its open development and use by immuno-informaticians, virologists, immunologists, vaccine developers, and other scientists that are interested in exploring and comparing multiple immune receptor repertoires.

Data Sharing and Reuse: A Method by the AIRR Community

High-throughput sequencing of adaptive immune receptor repertoires (AIRR, i.e., IG and TR ) has revolutionized the ability to study the adaptive immune response via large-scale experiments. Since 2009, AIRR sequencing (AIRR-seq) has been widely applied to survey the immune state of individuals (see "The AIRR Community Guide to Repertoire Analysis" chapter for details). One of the goals of the AIRR Community is to make the resulting AIRR-seq data FAIR (Findable, Accessible, Interoperable, and Reusable) (Wilkinson et al. Sci Data 3:1-9, 2016), with a primary goal of making it easy for the research community to reuse AIRR-seq data (Breden et al. Front Immunol 8:1418, 2017; Scott and Breden. Curr Opin Syst Biol 24:71-77, 2020). The basis for this is the MiAIRR data standard (Rubelt et al. Nat Immunol 18:1274-1278, 2017). For long-term preservation, it is recommended that researchers store their sequence read data in an INSDC repository. At the same time, the AIRR Community has established the AIRR Data Commons (Christley et al. Front Big Data 3:22, 2020), a distributed set of AIRR-compliant repositories that store the critically important annotated AIRR-seq data based on the MiAIRR standard, making the data findable, interoperable, and, because the data are annotated, more valuable in its reuse. Here, we build on the other AIRR Community chapters and illustrate how these principles and standards can be incorporated into AIRR-seq data analysis workflows. We discuss the importance of careful curation of metadata to ensure reproducibility and facilitate data sharing and reuse, and we illustrate how data can be shared via the AIRR Data Commons.

Adaptive Immune Receptor Repertoire (AIRR) Community Guide to TR and IG Gene Annotation

High-throughput sequencing of adaptive immune receptor repertoires (AIRR, i.e., IG and TR) has revolutionized the ability to carry out large-scale experiments to study the adaptive immune response. Since the method was first introduced in 2009, AIRR sequencing (AIRR-Seq) has been applied to survey the immune state of individuals, identify antigen-specific or immune-state-associated signatures of immune responses, study the development of the antibody immune response, and guide the development of vaccines and antibody therapies. Recent advancements in the technology include sequencing at the single-cell level and in parallel with gene expression, which allows the introduction of multi-omics approaches to understand in detail the adaptive immune response. Analyzing AIRR-seq data can prove challenging even with high-quality sequencing, in part due to the many steps involved and the need to parameterize each step. In this chapter, we outline key factors to consider when preprocessing raw AIRR-Seq data and annotating the genetic origins of the rearranged receptors. We also highlight a number of common difficulties with common AIRR-seq data processing and provide strategies to address them.

AIRR Community Guide to Planning and Performing AIRR-Seq Experiments

The development of high-throughput sequencing of adaptive immune receptor repertoires (AIRR-seq of IG and TR rearrangements) has provided a new frontier for in-depth analysis of the immune system. The last decade has witnessed an explosion in protocols, experimental methodologies, and computational tools. In this chapter, we discuss the major considerations in planning a successful AIRR-seq experiment together with basic strategies for controlling and evaluating the outcome of the experiment. Members of the AIRR Community have authored several chapters in this edition, which cover step-by-step instructions to successfully conduct, analyze, and share an AIRR-seq project.

High-throughput and single-cell T cell receptor sequencing technologies

T cells express T cell receptors (TCRs) composed of somatically recombined TCRα and TCRβ chains, which mediate recognition of major histocompatibility complex (MHC)-antigen complexes and drive the antigen-specific adaptive immune response to pathogens and cancer. The TCR repertoire in each individual is highly diverse, which allows for recognition of a wide array of foreign antigens, but also presents a challenge in analyzing this response using conventional methods. Recent studies have developed high-throughput sequencing technologies to identify TCR sequences, analyze their antigen specificities using experimental and computational tools, and pair TCRs with transcriptional and epigenetic cell state phenotypes in single cells. In this Review, we highlight these technological advances and describe how they have been applied to discover fundamental insights into T cell-mediated immunity.

Biological controls for standardization and interpretation of adaptive immune receptor repertoire profiling

Use of adaptive immune receptor repertoire sequencing (AIRR-seq) has become widespread, providing new insights into the immune system with potential broad clinical and diagnostic applications. However, like many high-throughput technologies, it comes with several problems, and the AIRR Community was established to understand and help solve them. We, the AIRR Community’s Biological Resources Working Group, have surveyed scientists about the need for standards and controls in generating and annotating AIRR-seq data. Here, we review the current status of AIRR-seq, provide the results of our survey, and based on them, offer recommendations for developing AIRR-seq standards and controls, including future work.

Opportunities and Challenges in Democratizing Immunology Datasets

The field of immunology is rapidly progressing toward a systems-level understanding of immunity to tackle complex infectious diseases, autoimmune conditions, cancer, and beyond. In the last couple of decades, advancements in data acquisition techniques have presented opportunities to explore untapped areas of immunological research. Broad initiatives are launched to disseminate the datasets siloed in the global, federated, or private repositories, facilitating interoperability across various research domains. Concurrently, the application of computational methods, such as network analysis, meta-analysis, and machine learning have propelled the field forward by providing insight into salient features that influence the immunological response, which was otherwise left unexplored. Here, we review the opportunities and challenges in democratizing datasets, repositories, and community-wide knowledge sharing tools. We present use cases for repurposing open-access immunology datasets with advanced machine learning applications and more.

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.

Benchmarking of T cell receptor repertoire profiling methods reveals large systematic biases

Monitoring the T cell receptor (TCR) repertoire in health and disease can provide key insights into adaptive immune responses, but the accuracy of current TCR sequencing (TCRseq) methods is unclear. In this study, we systematically compared the results of nine commercial and academic TCRseq methods, including six rapid amplification of complementary DNA ends (RACE)-polymerase chain reaction (PCR) and three multiplex-PCR approaches, when applied to the same T cell sample. We found marked differences in accuracy and intra- and inter-method reproducibility for T cell receptor α (TRA) and T cell receptor β (TRB) TCR chains. Most methods showed a lower ability to capture TRA than TRB diversity. Low RNA input generated non-representative repertoires. Results from the 5' RACE-PCR methods were consistent among themselves but differed from the RNA-based multiplex-PCR results. Using an in silico meta-repertoire generated from 108 replicates, we found that one genomic DNA-based method and two non-unique molecular identifier (UMI) RNA-based methods were more sensitive than UMI methods in detecting rare clonotypes, despite the better clonotype quantification accuracy of the latter.

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.

The adaptive immune receptor repertoire community as a model for FAIR stewardship of big immunology data

Systems biology involves network-oriented, computational approaches to modeling biological systems through analysis of big biological data. To contribute maximally to scientific progress, big biological data should be FAIR: findable, accessible, interoperable, and reusable. Here, we describe high-throughput sequencing data that characterize the vast diversity of B- and T-cell clones comprising the adaptive immune receptor repertoire (AIRR-seq data) and its contribution to our understanding of COVID-19 (coronavirus disease 19). We describe the accomplishments of the AIRR community, a grass-roots network of interdisciplinary laboratory scientists, bioinformaticians, and policy wonks, in creating and publishing standards, software and repositories for AIRR-seq data based on the FAIR principles.

immuneSIM: tunable multi-feature simulation of B- and T-cell receptor repertoires for immunoinformatics benchmarking

Summary

B- and T-cell receptor repertoires of the adaptive immune system have become a key target for diagnostics and therapeutics research. Consequently, there is a rapidly growing number of bioinformatics tools for immune repertoire analysis. Benchmarking of such tools is crucial for ensuring reproducible and generalizable computational analyses. Currently, however, it remains challenging to create standardized ground truth immune receptor repertoires for immunoinformatics tool benchmarking. Therefore, we developed immuneSIM, an R package that allows the simulation of native-like and aberrant synthetic full-length variable region immune receptor sequences by tuning the following immune receptor features: (i) species and chain type (BCR, TCR, single and paired), (ii) germline gene usage, (iii) occurrence of insertions and deletions, (iv) clonal abundance, (v) somatic hypermutation and (vi) sequence motifs. Each simulated sequence is annotated by the complete set of simulation events that contributed to its in silico generation. immuneSIM permits the benchmarking of key computational tools for immune receptor analysis, such as germline gene annotation, diversity and overlap estimation, sequence similarity, network architecture, clustering analysis and machine learning methods for motif detection.

Availability and implementation

The package is available via https://github.com/GreiffLab/immuneSIM and on CRAN at https://cran.r-project.org/web/packages/immuneSIM. The documentation is hosted at https://immuneSIM.readthedocs.io.

Contact

sai.reddy@ethz.ch or victor.greiff@medisin.uio.no

Supplementary information

Supplementary data are available at Bioinformatics online.

Benchmarking immunoinformatic tools for the analysis of antibody repertoire sequences

Summary

Antibody repertoires reveal insights into the biology of the adaptive immune system and empower diagnostics and therapeutics. There are currently multiple tools available for the annotation of antibody sequences. All downstream analyses such as choosing lead drug candidates depend on the correct annotation of these sequences; however, a thorough comparison of the performance of these tools has not been investigated. Here, we benchmark the performance of commonly used immunoinformatic tools, i.e. IMGT/HighV-QUEST, IgBLAST and MiXCR, in terms of reproducibility of annotation output, accuracy and speed using simulated and experimental high-throughput sequencing datasets.

We analyzed changes in IMGT reference germline database in the last 10 years in order to assess the reproducibility of the annotation output. We found that only 73/183 (40%) V, D and J human genes were shared between the reference germline sets used by the tools. We found that the annotation results differed between tools. In terms of alignment accuracy, MiXCR had the highest average frequency of gene mishits, 0.02 mishit frequency and IgBLAST the lowest, 0.004 mishit frequency. Reproducibility in the output of complementarity determining three regions (CDR3 amino acids) ranged from 4.3% to 77.6% with preprocessed data. In addition, run time of the tools was assessed: MiXCR was the fastest tool for number of sequences processed per unit of time. These results indicate that immunoinformatic analyses greatly depend on the choice of bioinformatics tool. Our results support informed decision-making to immunoinformaticians based on repertoire composition and sequencing platforms.

Availability and implementation

All tools utilized in the paper are free for academic use.

Supplementary information

Supplementary data are available at Bioinformatics online.

AncesTree: An interactive immunoglobulin lineage tree visualizer

High-throughput sequencing of human immunoglobulin genes allows analysis of antibody repertoires and the reconstruction of clonal lineage evolution. The study of antibodies (Abs) affinity maturation is of specific interest to understand the generation of Abs with high affinity or broadly neutralizing activities. Moreover, phylogenic analysis enables the identification of the key somatic mutations required to achieve optimal antigen binding. The Immcantation framework provides a start-to-finish set of analytical methods for high-throughput adaptive immune receptor repertoire sequencing (AIRR-Seq; Rep-Seq) data. Furthermore, Immcantation's Change-O package has developed IgPhyML, an algorithm designed to build specifically immunoglobulin (Ig) phylogenic trees. Meanwhile Phylip, an algorithm that has been originally developed for applications in ecology and macroevolution, can also be used for the phylogenic reconstruction of antibodies maturation pathway. To complement Ig lineages made by IgPhyML or Dnaml (Phylip), we developed AncesTree, a graphic user interface (GUI) that aims to give researchers the opportunity to interactively explore antibodies clonal evolution. AncesTree displays interactive immunoglobulins phylogenic tree, Ig related mutations and sequence alignments using additional information coming from specialized antibody tools. The GUI is a Java standalone application allowing interaction with Ig tree that can run under Windows, Linux and Mac OS.

RAbHIT: R Antibody Haplotype Inference Tool

SUMMARY

Antibody haplotype inference (chromosomal phasing) may have clinical implications for the identification of genetic predispositions to diseases. Yet, our knowledge of the genomic loci encoding for the variable regions of the antibody is only partial, mostly due to the challenge of aligning short reads from genome sequencing to these highly repetitive loci. A powerful approach to infer the content of these loci relies on analyzing repertoires of rearranged V(D)J sequences. We present here RAbHIT, an R Haplotype Antibody Inference Tool, that implements a novel algorithm to infer V(D)J haplotypes by adapting a Bayesian framework. RAbHIT offers inference of haplotype and gene deletions. It may be applied to sequences from naïve and non-naïve B-cells, sequenced by different library preparation protocols.

AVAILABILITY AND IMPLEMENTATION

RAbHIT is freely available for academic use from comprehensive R archive network (CRAN) (https://cran.r-project.org/web/packages/rabhit/) under CC BY-SA 4.0 license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The ADC API: A Web API for the Programmatic Query of the AIRR Data Commons

The Adaptive Immune Receptor Repertoire (AIRR) Community is a research-driven group that is establishing a clear set of community-accepted data and metadata standards; standards-based reference implementation tools; and policies and practices for infrastructure to support the deposit, curation, storage, and use of high-throughput sequencing data from B-cell and T-cell receptor repertoires (AIRR-seq data). The AIRR Data Commons is a distributed system of data repositories that utilizes a common data model, a common query language, and common interoperability formats for storage, query, and downloading of AIRR-seq data. Here is described the principal technical standards for the AIRR Data Commons consisting of the AIRR Data Model for repertoires and rearrangements, the AIRR Data Commons (ADC) API for programmatic query of data repositories, a reference implementation for ADC API services, and tools for querying and validating data repositories that support the ADC API. AIRR-seq data repositories can become part of the AIRR Data Commons by implementing the data model and API. The AIRR Data Commons allows AIRR-seq data to be reused for novel analyses and empowers researchers to discover new biological insights about the adaptive immune system.

Network Organization of Antibody Interactions in Sequence and Structure Space: the RADARS Model

Adaptive immunity in vertebrates is a complex self-organizing network of molecular interactions. While deep sequencing of the immune-receptor repertoire may reveal clonal relationships, functional interpretation of such data is hampered by the inherent limitations of converting sequence to structure to function. In this paper, a novel model of antibody interaction space and network, termed radial adjustment of system resolution, RAdial ADjustment of System Resolution (RADARS), is proposed. The model is based on the radial growth of interaction affinity of antibodies towards an infinity of directions in structure space, each direction corresponding to particular shapes of antigen epitopes. Levels of interaction affinity appear as free energy shells of the system, where hierarchical B-cell development and differentiation takes place. Equilibrium in this immunological thermodynamic system can be described by a power law distribution of antibody-free energies with an ideal network degree exponent of phi square, representing a scale-free fractal network of antibody interactions. Plasma cells are network hubs, memory B cells are nodes with intermediate degrees, and B1 cells function as nodes with minimal degree. Overall, the RADARS model implies that a finite number of antibody structures can interact with an infinite number of antigens by immunologically controlled adjustment of interaction energy distribution. Understanding quantitative network properties of the system should help the organization of sequence-derived predicted structural data.

Inter- and intraspecies comparison of phylogenetic fingerprints and sequence diversity of immunoglobulin variable genes

Protection and neutralization of a vast array of pathogens is accomplished by the tremendous diversity of the B cell receptor (BCR) repertoire. For jawed vertebrates, this diversity is initiated via the somatic recombination of immunoglobulin (Ig) germline elements. While it is clear that the number of these germline segments differs from species to species, the extent of cross-species sequence diversity remains largely uncharacterized. Here we use extensive computational and statistical methods to investigate the sequence diversity and evolutionary relationship between Ig variable (V), diversity (D), and joining (J) germline segments across nine commonly studied species ranging from zebrafish to human. Metrics such as guanine-cytosine (GC) content showed low redundancy across Ig germline genes within a given species. Other comparisons, including amino acid motifs, evolutionary selection, and sequence diversity, revealed species-specific properties. Additionally, we showed that the germline-encoded diversity differs across antibody (recombined V-D-J) repertoires of various B cell subsets. To facilitate future comparative immunogenomics analysis, we created VDJgermlines, an R package that contains the germline sequences from multiple species. Our study informs strategies for the humanization and engineering of therapeutic antibodies.

Structural diversity of B-cell receptor repertoires along the B-cell differentiation axis in humans and mice

Most current analysis tools for antibody next-generation sequencing data work with primary sequence descriptors, leaving accompanying structural information unharnessed. We have used novel rapid methods to structurally characterize the complementary-determining regions (CDRs) of more than 180 million human and mouse B-cell receptor (BCR) repertoire sequences. These structurally annotated CDRs provide unprecedented insights into both the structural predetermination and dynamics of the adaptive immune response. We show that B-cell types can be distinguished based solely on these structural properties. Antigen-unexperienced BCR repertoires use the highest number and diversity of CDR structures and these patterns of naïve repertoire paratope usage are highly conserved across subjects. In contrast, more differentiated B-cells are more personalized in terms of CDR structure usage. Our results establish the CDR structure differences in BCR repertoires and have applications for many fields including immunodiagnostics, phage display library generation, and “humanness” assessment of BCR repertoires from transgenic animals. The software tool for structural annotation of BCR repertoires, SAAB+, is available at https://github.com/oxpig/saab_plus.

B-cell receptors (BCR) are the major components of the adaptive immune system. These are immunoglobulin molecules that bind to foreign substances known as antigens. Each individual has a huge BCR repertoire, where each individual BCR has a specific binding site composed of the complementary-determining regions (CDRs) capable of recognising a specific antigen. Drug discovery and immunodiagnostics inspired by the adaptive immune system rely on our ability to accurately interrogate the structural diversity of the binding sites of the BCR repertoire. Here we report our novel rapid pipeline, SAAB+, which has enabled us to interrogate how the structure of the CDR changes in BCR repertoires along the B-cell differentiation axis. By analysing human and mouse BCR repertoires at an unprecedented scale, we observed species-specific structural predetermination and detected CDR dynamics across multiple stages of B-cell differentiation. We showed that naïve repertoires share the highest number and diversity of CDR structures, a pattern which was highly conserved in all B-cell donors. Our results suggest that increased B-cell differentiation is associated with a personalization of CDR structure usages. Finally, we established the differences in CDR usages between humans and mice, analysis with immediate relevance for BCR repertoire “humanness” assessment and rational immunotherapeutic engineering.

VDJbase: an adaptive immune receptor genotype and haplotype database

VDJbase is a publicly available database that offers easy searching of data describing the complete sets of gene sequences (genotypes and haplotypes) inferred from adaptive immune receptor repertoire sequencing datasets. VDJbase is designed to act as a resource that will allow the scientific community to explore the genetic variability of the immunoglobulin (Ig) and T cell receptor (TR) gene loci. It can also assist in the investigation of Ig- and TR-related genetic predispositions to diseases. Our database includes web-based query and online tools to assist in visualization and analysis of the genotype and haplotype data. It enables users to detect those alleles and genes that are significantly over-represented in a particular population, in terms of genotype, haplotype and gene expression. The database website can be freely accessed at https://www.vdjbase.org/, and no login is required. The data and code use creative common licenses and are freely downloadable from https://bitbucket.org/account/user/yaarilab/projects/GPHP.

The Diagnostic, Prognostic, and Therapeutic Potential of Adaptive Immune Receptor Repertoire Profiling in Cancer

Lymphocytes play a critical role in antitumor immune responses. They are directly targeted by some therapies, and the composition and spatial organization of intratumor T-cell populations is prognostic in some cancer types. A better understanding of lymphocyte population dynamics over the course of disease and in response to therapy is urgently needed to guide therapy decisions and to develop new therapy targets. Deep sequencing of the repertoire of antigen receptor-encoding genes expressed in a lymphocyte population has become a widely used approach for profiling the population's immune status. Lymphocyte antigen receptor repertoire deep sequencing data can be used to assess the clonal richness and diversity of lymphocyte populations; to track clone members over time, between tissues, and across lymphocyte subsets; to detect clonal expansion; and to detect the recruitment of new clones into a tissue. Repertoire sequencing is thus a critical complement to other methods of lymphocyte and immune profiling in cancer. This review describes the current state of knowledge based on repertoire sequencing studies conducted on human cancer patients, with a focus on studies of the T-cell receptor beta chain locus. The review then outlines important questions left unanswered and suggests future directions for the field.

The Immune Epitope Database and Analysis Resource Program 2003-2018: reflections and outlook

The Immune Epitope Database and Analysis Resource (IEDB) contains information related to antibodies and T cells across an expansive scope of research fields (infectious diseases, allergy, autoimmunity, and transplantation). Capture and representation of the data to reflect growing scientific standards and techniques have required continual refinement of our rigorous curation and query and reporting processes beginning with the automated classification of over 28 million PubMed abstracts, and resulting in easily searchable data from over 20,000 published manuscripts. Data related to MHC binding and elution, nonpeptidics, natural processing, receptors, and 3D structure is first captured through manual curation and subsequently maintained through recuration to reflect evolving scientific standards. Upon promotion to the free, public database, users can query and export records of specific relevance via the online web portal which undergoes iterative development to best enable efficient data access. In parallel, the companion Analysis Resource site hosts a variety of tools that assist in the bioinformatic analyses of epitopes and related structures, which can be applied to IEDB-derived and independent datasets alike. Available tools are classified into two categories: analysis and prediction. Analysis tools include epitope clustering, sequence conservancy, and more, while prediction tools cover T and B cell epitope binding, immunogenicity, and TCR/BCR structures. In addition to these tools, benchmarking servers which allow for unbiased performance comparison are also offered. In order to expand and support the user-base of both the database and Analysis Resource, the research team actively engages in community outreach through publication of ongoing work, conference attendance and presentations, hosting of user workshops, and the provision of online help. This review provides a description of the IEDB database infrastructure, curation and recuration processes, query and reporting capabilities, the Analysis Resource, and our Community Outreach efforts, including assessment of the impact of the IEDB across the research community.

Looking for therapeutic antibodies in next-generation sequencing repositories

Recently it has become possible to query the great diversity of natural antibody repertoires using next-generation sequencing (NGS). These methods are capable of producing millions of sequences in a single experiment. Here we compare clinical-stage therapeutic antibodies to the ~1b sequences from 60 independent sequencing studies in the Observed Antibody Space database, which includes antibody sequences from NGS analysis of immunoglobulin gene repertoires. Of 242 post-Phase 1 antibodies, we found 16 with sequence identity matches of 95% or better for both heavy and light chains. There are also 54 perfect matches to therapeutic CDR-H3 regions in the NGS outputs, suggesting a nontrivial amount of convergence between naturally observed sequences and those developed artificially. This has potential implications for both the legal protection of commercial antibodies and the discovery of antibody therapeutics.

Deep generative models for T cell receptor protein sequences

Probabilistic models of adaptive immune repertoire sequence distributions can be used to infer the expansion of immune cells in response to stimulus, differentiate genetic from environmental factors that determine repertoire sharing, and evaluate the suitability of various target immune sequences for stimulation via vaccination. Classically, these models are defined in terms of a probabilistic V(D)J recombination model which is sometimes combined with a selection model. In this paper we take a different approach, fitting variational autoencoder (VAE) models parameterized by deep neural networks to T cell receptor (TCR) repertoires. We show that simple VAE models can perform accurate cohort frequency estimation, learn the rules of VDJ recombination, and generalize well to unseen sequences. Further, we demonstrate that VAE-like models can distinguish between real sequences and sequences generated according to a recombination-selection model, and that many characteristics of VAE-generated sequences are similar to those of real sequences.

Recent progress in the analysis of αβT cell and B cell receptor repertoires

T cell receptors (TCRs) and B cell receptors (BCRs) are vertebrate evolution's best answer to the threat of microbial pathogens that can evolve much faster than ourselves. These antigen receptors are generated during T cell or B cell development by combinatorial rearrangement of germline genome V, D and J gene segments, and with junctional residues capable of enormous diversity. For decades the complexity of these receptor repertoires has limited their analysis, but advances in DNA sequencing technology and an array of complementary tools have now made their study much more tractable, filling a major gap in our ability to understand immunology as a system. Here, we summarize the recent approaches and discoveries that are enabling these advances, with some suggestions as to what may lie ahead.

Application of immune repertoire sequencing in cancer immunotherapy

With the prominent breakthrough in the field of tumor immunology, diverse cancer immunotherapies have attracted great attention in the last decade. The immune checkpoint inhibitors, adoptive cell therapies, and therapeutic cancer vaccines have already achieved impressive clinical success. However, the fact that only a small subset of patients with specific tumor types can benefit from these treatments limits the application of cancer immunotherapy. To seek out the molecular mechanisms behind this challenge and to select cancer precision medicine for different individuals, researchers apply the immune repertoire sequencing (IRS) to evaluate genetic responses of each patient to current immunotherapies. This review summarizes the technical advances and recent applications of IRS in cancer immunotherapy, indicates the limitations of this technique, and predicts future perspectives both in basic studies and clinical trials.

High-Quality Library Preparation for NGS-Based Immunoglobulin Germline Gene Inference and Repertoire Expression Analysis

Next generation sequencing (NGS) of immunoglobulin (Ig) repertoires (Rep-seq) enables examination of the adaptive immune system at an unprecedented level. Applications include studies of expressed repertoires, gene usage, somatic hypermutation levels, Ig lineage tracing and identification of genetic variation within the Ig loci through inference methods. All these applications require starting libraries that allow the generation of sequence data with low error rate and optimal representation of the expressed repertoire. Here, we provide detailed protocols for the production of libraries suitable for human Ig germline gene inference and Ig repertoire studies. Various parameters used in the process were tested in order to demonstrate factors that are critical to obtain high quality libraries. We demonstrate an improved 5'RACE technique that reduces the length constraints of Illumina MiSeq based Rep-seq analysis but allows for the acquisition of sequences upstream of Ig V genes, useful for primer design. We then describe a 5' multiplex method for library preparation, which yields full length V(D)J sequences suitable for genotype identification and novel gene inference. We provide comprehensive sets of primers targeting IGHV, IGKV, and IGLV genes. Using the optimized protocol, we produced IgM, IgG, IgK, and IgL libraries and analyzed them using the germline inference tool IgDiscover to identify expressed germline V alleles. This process additionally uncovered three IGHV, one IGKV, and six IGLV novel alleles in a single individual, which are absent from the IMGT reference database, highlighting the need for further study of Ig genetic variation. The library generation protocols presented here enable a robust means of analyzing expressed Ig repertoires, identifying novel alleles and producing individualized germline gene databases from humans.

Human adaptive immune receptor repertoire analysis-Past, present, and future

The genes encoding adaptive immune antigen receptors, namely the immunoglobulins expressed in membrane-bound or secreted forms by B cells, and the cell surface T cell receptors, are unique in human biology because they are generated by combinatorial rearrangement of the genomic DNA. The diversity of receptors so generated in populations of lymphocytes enables the human immune system to recognize antigens expressed by pathogens, but also underlies the pathological specificity of autoimmune diseases and the mistargeted immunity in allergies. Several recent technological developments, foremost among them the invention of high-throughput DNA sequencing instruments, have enabled much deeper and thorough evaluation of clones of human B cells and T cells and the antigen receptors they express during physiological and pathogenic immune responses. The evolutionary struggles between host adaptive immune responses and populations of pathogens are now open to greater scrutiny, elucidation of the underlying reasons for successful or failed immunity, and potential predictive modeling, than ever before. Here we give an overview of the foundations, recent progress, and future prospects in this dynamic area of research.

iReceptor: A platform for querying and analyzing antibody/B-cell and T-cell receptor repertoire data across federated repositories

Next-generation sequencing allows the characterization of the adaptive immune receptor repertoire (AIRR) in exquisite detail. These large-scale AIRR-seq data sets have rapidly become critical to vaccine development, understanding the immune response in autoimmune and infectious disease, and monitoring novel therapeutics against cancer. However, at present there is no easy way to compare these AIRR-seq data sets across studies and institutions. The ability to combine and compare information for different disease conditions will greatly enhance the value of AIRR-seq data for improving biomedical research and patient care. The iReceptor Data Integration Platform (gateway.ireceptor.org) provides one implementation of the AIRR Data Commons envisioned by the AIRR Community (airr-community.org), an initiative that is developing protocols to facilitate sharing and comparing AIRR-seq data. The iReceptor Scientific Gateway links distributed (federated) AIRR-seq repositories, allowing sequence searches or metadata queries across multiple studies at multiple institutions, returning sets of sequences fulfilling specific criteria. We present a review of the development of iReceptor, and how it fits in with the general trend toward sharing genomic and health data, and the development of standards for describing and reporting AIRR-seq data. Researchers interested in integrating their repositories of AIRR-seq data into the iReceptor Platform are invited to contact support@ireceptor.org.

Inferred Allelic Variants of Immunoglobulin Receptor Genes: A System for Their Evaluation, Documentation, and Naming

Immunoglobulins or antibodies are the main effector molecules of the B-cell lineage and are encoded by hundreds of variable (V), diversity (D), and joining (J) germline genes, which recombine to generate enormous IG diversity. Recently, high-throughput adaptive immune receptor repertoire sequencing (AIRR-seq) of recombined V-(D)-J genes has offered unprecedented insights into the dynamics of IG repertoires in health and disease. Faithful biological interpretation of AIRR-seq studies depends upon the annotation of raw AIRR-seq data, using reference germline gene databases to identify the germline genes within each rearrangement. Existing reference databases are incomplete, as shown by recent AIRR-seq studies that have inferred the existence of many previously unreported polymorphisms. Completing the documentation of genetic variation in germline gene databases is therefore of crucial importance. Lymphocyte receptor genes and alleles are currently assigned by the Immunoglobulins, T cell Receptors and Major Histocompatibility Nomenclature Subcommittee of the International Union of Immunological Societies (IUIS) and managed in IMGT®, the international ImMunoGeneTics information system® (IMGT). In 2017, the IMGT Group reached agreement with a group of AIRR-seq researchers on the principles of a streamlined process for identifying and naming inferred allelic sequences, for their incorporation into IMGT®. These researchers represented the AIRR Community, a network of over 300 researchers whose objective is to promote all aspects of immunoglobulin and T-cell receptor repertoire studies, including the standardization of experimental and computational aspects of AIRR-seq data generation and analysis. The Inferred Allele Review Committee (IARC) was established by the AIRR Community to devise policies, criteria, and procedures to perform this function. Formalized evaluations of novel inferred sequences have now begun and submissions are invited via a new dedicated portal (https://ogrdb.airr-community.org). Here, we summarize recommendations developed by the IARC-focusing, to begin with, on human IGHV genes-with the goal of facilitating the acceptance of inferred allelic variants of germline IGHV genes. We believe that this initiative will improve the quality of AIRR-seq studies by facilitating the description of human IG germline gene variation, and that in time, it will expand to the documentation of TR and IG genes in many vertebrate species.

Identification of Subject-Specific Immunoglobulin Alleles From Expressed Repertoire Sequencing Data

The adaptive immune receptor repertoire (AIRR) contains information on an individuals' immune past, present and potential in the form of the evolving sequences that encode the B cell receptor (BCR) repertoire. AIRR sequencing (AIRR-seq) studies rely on databases of known BCR germline variable (V), diversity (D), and joining (J) genes to detect somatic mutations in AIRR-seq data via comparison to the best-aligning database alleles. However, it has been shown that these databases are far from complete, leading to systematic misidentification of mutated positions in subsets of sample sequences. We previously presented TIgGER, a computational method to identify subject-specific V gene genotypes, including the presence of novel V gene alleles, directly from AIRR-seq data. However, the original algorithm was unable to detect alleles that differed by more than 5 single nucleotide polymorphisms (SNPs) from a database allele. Here we present and apply an improved version of the TIgGER algorithm which can detect alleles that differ by any number of SNPs from the nearest database allele, and can construct subject-specific genotypes with minimal prior information. TIgGER predictions are validated both computationally (using a leave-one-out strategy) and experimentally (using genomic sequencing), resulting in the addition of three new immunoglobulin heavy chain V (IGHV) gene alleles to the IMGT repertoire. Finally, we develop a Bayesian strategy to provide a confidence estimate associated with genotype calls. All together, these methods allow for much higher accuracy in germline allele assignment, an essential step in AIRR-seq studies.

Mosaic deletion patterns of the human antibody heavy chain gene locus shown by Bayesian haplotyping

Analysis of antibody repertoires by high-throughput sequencing is of major importance in understanding adaptive immune responses. Our knowledge of variations in the genomic loci encoding immunoglobulin genes is incomplete, resulting in conflicting VDJ gene assignments and biased genotype and haplotype inference. Haplotypes can be inferred using IGHJ6 heterozygosity, observed in one third of the people. Here, we propose a robust novel method for determining VDJ haplotypes by adapting a Bayesian framework. Our method extends haplotype inference to IGHD- and IGHV-based analysis, enabling inference of deletions and copy number variations in the entire population. To test this method, we generated a multi-individual data set of naive B-cell repertoires, and found allele usage bias, as well as a mosaic, tiled pattern of deleted IGHD and IGHV genes. The inferred haplotypes may have clinical implications for genetic disease predispositions. Our findings expand the knowledge that can be extracted from antibody repertoire sequencing data.

Using association rule mining and ontologies to generate metadata recommendations from multiple biomedical databases

Metadata-the machine-readable descriptions of the data-are increasingly seen as crucial for describing the vast array of biomedical datasets that are currently being deposited in public repositories. While most public repositories have firm requirements that metadata must accompany submitted datasets, the quality of those metadata is generally very poor. A key problem is that the typical metadata acquisition process is onerous and time consuming, with little interactive guidance or assistance provided to users. Secondary problems include the lack of validation and sparse use of standardized terms or ontologies when authoring metadata. There is a pressing need for improvements to the metadata acquisition process that will help users to enter metadata quickly and accurately. In this paper, we outline a recommendation system for metadata that aims to address this challenge. Our approach uses association rule mining to uncover hidden associations among metadata values and to represent them in the form of association rules. These rules are then used to present users with real-time recommendations when authoring metadata. The novelties of our method are that it is able to combine analyses of metadata from multiple repositories when generating recommendations and can enhance those recommendations by aligning them with ontology terms. We implemented our approach as a service integrated into the CEDAR Workbench metadata authoring platform, and evaluated it using metadata from two public biomedical repositories: US-based National Center for Biotechnology Information BioSample and European Bioinformatics Institute BioSamples. The results show that our approach is able to use analyses of previously entered metadata coupled with ontology-based mappings to present users with accurate recommendations when authoring metadata.

Antigen Receptor Sequence Reconstruction and Clonality Inference from scRNA-Seq Data

In this chapter, we describe TraCeR and BraCeR, our computational tools for reconstruction of paired full-length antigen receptor sequences and clonality inference from single-cell RNA-seq (scRNA-seq) data. In brief, TraCeR reconstructs T-cell receptor (TCR) sequences from scRNA-seq data by extracting sequencing reads derived from TCRs by aligning the reads from each cell against synthetic TCR sequences. TCR-derived reads are then assembled into full-length recombined TCR sequences. BraCeR builds on the TraCeR pipeline and accounts for somatic hypermutations (SHM) and isotype switching. Here we discuss experimental design, use of the tools, and interpretation of the results.

Expression Cloning of Antibodies from Single Human B Cells

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

Epitope Specific Antibodies and T Cell Receptors in the Immune Epitope Database

The Immune Epitope Database (IEDB) is a free public resource which catalogs experiments characterizing immune epitopes. To accommodate data from next generation repertoire sequencing experiments, we recently updated how we capture and query epitope specific antibodies and T cell receptors. Specifically, we are now storing partial receptor sequences sufficient to determine CDRs and VDJ gene usage which are commonly identified by repertoire sequencing. For previously captured full length receptor sequencing data, we have calculated the corresponding CDR sequences and gene usage information using IMGT numbering and VDJ gene nomenclature format. To integrate information from receptors defined at different levels of resolution, we grouped receptors based on their host species, receptor type and CDR3 sequence. As of August 2018, we have cataloged sequence information for more than 22,510 receptors in 18,292 receptor groups, shown to bind to more than 2,241 distinct epitopes. These data are accessible as full exports and through a new dedicated query interface. The later combines the new ability to search by receptor characteristics with previously existing capability to search by epitope characteristics such as the infectious agent the epitope is derived from, or the kind of immune response involved in its recognition. We expect that this comprehensive capture of epitope specific immune receptor information will provide new insights into receptor-epitope interactions, and facilitate the development of novel tools that help in the analysis of receptor repertoire data.

Predicting B cell receptor substitution profiles using public repertoire data

B cells develop high affinity receptors during the course of affinity maturation, a cyclic process of mutation and selection. At the end of affinity maturation, a number of cells sharing the same ancestor (i.e. in the same “clonal family”) are released from the germinal center; their amino acid frequency profile reflects the allowed and disallowed substitutions at each position. These clonal-family-specific frequency profiles, called “substitution profiles”, are useful for studying the course of affinity maturation as well as for antibody engineering purposes. However, most often only a single sequence is recovered from each clonal family in a sequencing experiment, making it impossible to construct a clonal-family-specific substitution profile. Given the public release of many high-quality large B cell receptor datasets, one may ask whether it is possible to use such data in a prediction model for clonal-family-specific substitution profiles. In this paper, we present the method “Substitution Profiles Using Related Families” (SPURF), a penalized tensor regression framework that integrates information from a rich assemblage of datasets to predict the clonal-family-specific substitution profile for any single input sequence. Using this framework, we show that substitution profiles from similar clonal families can be leveraged together with simulated substitution profiles and germline gene sequence information to improve prediction. We fit this model on a large public dataset and validate the robustness of our approach on two external datasets. Furthermore, we provide a command-line tool in an open-source software package (https://github.com/krdav/SPURF) implementing these ideas and providing easy prediction using our pre-fit models.

Antibody engineering can be greatly informed by knowledge about the underlying affinity maturation process. As such this can be probed by sequencing, but unfortunately, in practice often only one member of the clonal family is sequenced, making it difficult to determine a set of possible amino acid mutations that would retain the original antibody antigen binding affinity. We overcome this data sparsity by developing a statistical learning approach that leverages vast information about amino acid preferences available in public immune system repertoire data. We use a penalized regression approach to devise a flexible statistical model that integrates multiple sources of information into a coherent prediction framework and validate our prediction algorithm using subsampling and held out data.

Sequencing the peripheral blood B and T cell repertoire - Quantifying robustness and limitations

The adaptive immune response generates a large repertoire of T cells with T-cell receptors (TCR_alpha and TCR_beta) and B cells with immunoglobulins (Ig). The repertoire changes in response to antigen stimulation both through amplification of specific cells (clonal expansion) as well as somatic hypermutation of immunoglobulins. Alterations of the immune repertoire have been observed in response to acute disease, such as external pathogens, or chronic diseases, such as autoimmunity and cancer. Here we establish experimental and analytical protocols for quantifying the peripheral blood of healthy human individuals by profiling the immune repertoire for the Complementarity determining region 3 (CDR3) of the variable regions of TCR_beta (CDRβ3) and the IgG heavy chain (CDRH1, CDRH2, CDRH3). The results demonstrate that 40 ml of blood are sufficient to reliably capture the 10,000 most common TCR_beta and 1000 most common IgG and determine their relative frequency in the circulation. We conclude that by using an accessible sample size of human PBMC one is able to robustly monitor alterations in the immune repertoire.

AIRR Community Standardized Representations for Annotated Immune Repertoires

Increased interest in the immune system's involvement in pathophysiological phenomena coupled with decreased DNA sequencing costs have led to an explosion of antibody and T cell receptor sequencing data collectively termed "adaptive immune receptor repertoire sequencing" (AIRR-seq or Rep-Seq). The AIRR Community has been actively working to standardize protocols, metadata, formats, APIs, and other guidelines to promote open and reproducible studies of the immune repertoire. In this paper, we describe the work of the AIRR Community's Data Representation Working Group to develop standardized data representations for storing and sharing annotated antibody and T cell receptor data. Our file format emphasizes ease-of-use, accessibility, scalability to large data sets, and a commitment to open and transparent science. It is composed of a tab-delimited format with a specific schema. Several popular repertoire analysis tools and data repositories already utilize this AIRR-seq data format. We hope that others will follow suit in the interest of promoting interoperable standards.