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

Serum proteomics reveals high-affinity and convergent antibodies by tracking SARS-CoV-2 hybrid immunity to emerging variants of concern

The rapid spread of SARS-CoV-2 and its continuing impact on human health has prompted the need for effective and rapid development of monoclonal antibody therapeutics. In this study, we investigate polyclonal antibodies in serum and B cells from the whole blood of three donors with SARS-CoV-2 immunity to find high-affinity anti-SARS-CoV-2 antibodies to escape variants. Serum IgG antibodies were selected by their affinity to the receptor-binding domain (RBD) and non-RBD sites on the spike protein of Omicron subvariant B.1.1.529 from each donor. Antibodies were analyzed by bottom-up mass spectrometry, and matched to single- and bulk-cell sequenced repertoires for each donor. The antibodies observed in serum were recombinantly expressed, and characterized to assess domain binding, cross-reactivity between different variants, and capacity to inhibit RBD binding to host protein. Donors infected with early Omicron subvariants had serum antibodies with subnanomolar affinity to RBD that also showed binding activity to a newer Omicron subvariant BQ.1.1. The donors also showed a convergent immune response. Serum antibodies and other single- and bulk-cell sequences were similar to publicly reported anti-SARS-CoV-2 antibodies, and the characterized serum antibodies had the same variant-binding and neutralization profiles as their reported public sequences. The serum antibodies analyzed were a subset of anti-SARS-CoV-2 antibodies in the B cell repertoire, which demonstrates significant dynamics between the B cells and circulating antibodies in peripheral blood.

Bioinformatics in Russia: history and present-day landscape

Bioinformatics has become an interdisciplinary subject due to its universal role in molecular biology research. The current status of Russia's bioinformatics research in Russia is not known. Here, we review the history of bioinformatics in Russia, present the current landscape, and highlight future directions and challenges. Bioinformatics research in Russia is driven by four major industries: information technology, pharmaceuticals, biotechnology, and agriculture. Over the past three decades, despite a delayed start, the field has gained momentum, especially in protein and nucleic acid research. Dedicated and shared centers for genomics, proteomics, and bioinformatics are active in different regions of Russia. Present-day bioinformatics in Russia is characterized by research issues related to genetics, metagenomics, OMICs, medical informatics, computational biology, environmental informatics, and structural bioinformatics. Notable developments are in the fields of software (tools, algorithms, and pipelines), use of high computation power (e.g. by the Siberian Supercomputer Center), and large-scale sequencing projects (the sequencing of 100 000 human genomes). Government funding is increasing, policies are being changed, and a National Genomic Information Database is being established. An increased focus on eukaryotic genome sequencing, the development of a common place for developers and researchers to share tools and data, and the use of biological modeling, machine learning, and biostatistics are key areas for future focus. Universities and research institutes have started to implement bioinformatics modules. A critical mass of bioinformaticians is essential to catch up with the global pace in the discipline.

A monoclonal Trd chain supports the development of the complete set of functional γδ T cell lineages

The clonal selection theory describes key features of adaptive immune responses of B and T cells. For αβ T cells and B cells, antigen recognition and selection principles are known at a detailed molecular level. The precise role of the antigen receptor in γδ T cells remains less well understood. To better understand the role of the γδ T cell receptor (TCR), we generate an orthotopic TCRδ transgenic mouse model. We demonstrate a multi-layered functionality of γδ TCRs and diverse roles of CDR3δ-mediated selection during γδ T cell development. Whereas epithelial populations using Vγ5 or Vγ7 chains are almost unaffected in their biology in the presence of the transgenic TCRδ chain, pairing with Vγ1 positively selects γδ T cell subpopulations with distinct programs in several organs, thereby distorting the repertoire. In conclusion, our data support dictation of developmental tropism together with adaptive-like recognition principles in a single antigen receptor.

A scalable model for simulating multi-round antibody evolution and benchmarking of clonal tree reconstruction methods

Affinity maturation (AM) of B cells through somatic hypermutations (SHMs) enables the immune system to evolve to recognize diverse pathogens. The accumulation of SHMs leads to the formation of clonal lineages of antibody-secreting b cells that have evolved from a common naïve B cell. Advances in high-throughput sequencing have enabled deep scans of B cell receptor repertoires, paving the way for reconstructing clonal trees. However, it is not clear if clonal trees, which capture microevolutionary time scales, can be reconstructed using traditional phylogenetic reconstruction methods with adequate accuracy. In fact, several clonal tree reconstruction methods have been developed to fix supposed shortcomings of phylogenetic methods. Nevertheless, no consensus has been reached regarding the relative accuracy of these methods, partially because evaluation is challenging. Benchmarking the performance of existing methods and developing better methods would both benefit from realistic models of clonal lineage evolution specifically designed for emulating B cell evolution. In this paper, we propose a model for modeling B cell clonal lineage evolution and use this model to benchmark several existing clonal tree reconstruction methods. Our model, designed to be extensible, has several features: by evolving the clonal tree and sequences simultaneously, it allows modeling selective pressure due to changes in affinity binding; it enables scalable simulations of large numbers of cells; it enables several rounds of infection by an evolving pathogen; and, it models building of memory. In addition, we also suggest a set of metrics for comparing clonal trees and measuring their properties. Our results show that while maximum likelihood phylogenetic reconstruction methods can fail to capture key features of clonal tree expansion if applied naively, a simple post-processing of their results, where short branches are contracted, leads to inferences that are better than alternative methods.

New Views of Old Proteins: Clarifying the Enigmatic Proteome

All human diseases involve proteins, yet our current tools to characterize and quantify them are limited. To better elucidate proteins across space, time, and molecular composition, we provide a >10 years of projection for technologies to meet the challenges that protein biology presents. With a broad perspective, we discuss grand opportunities to transition the science of proteomics into a more propulsive enterprise. Extrapolating recent trends, we describe a next generation of approaches to define, quantify, and visualize the multiple dimensions of the proteome, thereby transforming our understanding and interactions with human disease in the coming decade.

Antibody repertoire sequencing analysis

High-throughput sequencing for B cell receptor (BCR) repertoire provides useful insights for the adaptive immune system. With the continuous development of the BCR-seq technology, many efforts have been made to develop methods for analyzing the ever-increasing BCR repertoire data. In this review, we comprehensively outline different BCR repertoire library preparation protocols and summarize three major steps of BCR-seq data analysis, i. e., V(D)J sequence annotation, clonal phylogenetic inference, and BCR repertoire profiling and mining. Different from other reviews in this field, we emphasize background intuition and the statistical principle of each method to help biologists better understand it. Finally, we discuss data mining problems for BCR-seq data and with a highlight on recently emerging multiple-sample analysis.

Gene prediction in the immunoglobulin loci

The V(D)J recombination process rearranges the variable (V), diversity (D), and joining (J) genes in the immunoglobulin loci to generate antibody repertoires. Annotation of these loci across various species and predicting the V, D, and J genes (IG genes) is critical for studies of the adaptive immune system. However, since the standard gene finding algorithms are not suitable for predicting IG genes, they have been semi-manually annotated in very few species. We developed the IGDetective algorithm for predicting IG genes and applied it to species with the assembled IG loci. IGDetective generated the first large collection of IG genes across many species and enabled their evolutionary analysis, including the analysis of the "bat IG diversity" hypothesis. This analysis revealed extremely conserved V genes in evolutionary distant species indicating that these genes may be subjected to the same selective pressure, e.g., pressure driven by common pathogens. IGDetective also revealed extremely diverged V genes and a new family of evolutionary conserved V genes in bats with unusual noncanonical cysteines. Moreover, in difference from all other previously reported antibodies, these cysteines are located within complementarity-determining regions. Since cysteines form disulfide bonds, we hypothesize that these cysteine-rich V genes might generate antibodies with noncanonical conformations and could potentially form a unique part of the immune repertoire in bats. We also analyzed the diversity landscape of the recombination signal sequences and revealed their features that trigger the high/low usage of the IG genes.

Variations in antibody repertoires correlate with vaccine responses

An important challenge in vaccine development is to figure out why a vaccine succeeds in some individuals and fails in others. Although antibody repertoires hold the key to answering this question, there have been very few personalized immunogenomics studies so far aimed at revealing how variations in immunoglobulin genes affect a vaccine response. We conducted an immunosequencing study of 204 calves vaccinated against bovine respiratory disease (BRD) with the goal to reveal variations in immunoglobulin genes and somatic hypermutations that impact the efficacy of vaccine response. Our study represents the largest longitudinal personalized immunogenomics study reported to date across all species, including humans. To analyze the generated data set, we developed an algorithm for identifying variations of the immunoglobulin genes (as well as frequent somatic hypermutations) that affect various features of the antibody repertoire and titers of neutralizing antibodies. In contrast to relatively short human antibodies, cattle have a large fraction of ultralong antibodies that have opened new therapeutic opportunities. Our study reveals that ultralong antibodies are a key component of the immune response against the costliest disease of beef cattle in North America. The detected variants of the cattle immunoglobulin genes, which are implicated in the success/failure of the BRD vaccine, have the potential to direct the selection of individual cattle for ongoing breeding programs.

A perspective toward mass spectrometry-based de novo sequencing of endogenous antibodies

A key step in therapeutic and endogenous humoral antibody characterization is identifying the amino acid sequence. So far, this task has been mainly tackled through sequencing of B-cell receptor (BCR) repertoires at the nucleotide level. Mass spectrometry (MS) has emerged as an alternative tool for obtaining sequence information directly at the – most relevant – protein level. Although several MS methods are now well established, analysis of recombinant and endogenous antibodies comes with a specific set of challenges, requiring approaches beyond the conventional proteomics workflows. Here, we review the challenges in MS-based sequencing of both recombinant as well as endogenous humoral antibodies and outline state-of-the-art methods attempting to overcome these obstacles. We highlight recent examples and discuss remaining challenges. We foresee a great future for these approaches making de novo antibody sequencing and discovery by MS-based techniques feasible, even for complex clinical samples from endogenous sources such as serum and other liquid biopsies.

Targeting B cells in the pre-phase of systemic autoimmunity globally interferes with autoimmune pathology

Systemic lupus erythematosus (SLE) is characterized by a loss of self-tolerance, systemic inflammation, and multi-organ damage. While a variety of therapeutic interventions are available, it has become clear that an early diagnosis and treatment may be key to achieve long lasting therapeutic responses and to limit irreversible organ damage. Loss of humoral tolerance including the appearance of self-reactive antibodies can be detected years before the actual onset of the clinical autoimmune disease, representing a potential early point of intervention. Not much is known, however, about how and to what extent this pre-phase of disease impacts the onset and development of subsequent autoimmunity. By targeting the B cell compartment in the pre-disease phase of a spontaneous mouse model of SLE we now show, that resetting the humoral immune system during the clinically unapparent phase of the disease globally alters immune homeostasis delaying the downstream development of systemic autoimmunity.

•

The clinically unapparent pre-phase of SLE impacts clinical disease

•

Autoreactive IgM antibodies represent a biomarker for early therapeutic intervention

•

Pre-phase B cells orchestrate clinical disease

•

Depleting pre-phase B cells diminishes disease pathology

Progress and challenges in mass spectrometry-based analysis of antibody repertoires

Humoral immunity is divided into the cellular B cell and protein-level antibody responses. High-throughput sequencing has advanced our understanding of both these fundamental aspects of B cell immunology as well as aspects pertaining to vaccine and therapeutics biotechnology. Although the protein-level serum and mucosal antibody repertoire make major contributions to humoral protection, the sequence composition and dynamics of antibody repertoires remain underexplored. This limits insight into important immunological and biotechnological parameters such as the number of antigen-specific antibodies, which are for example, relevant for pathogen neutralization, microbiota regulation, severity of autoimmunity, and therapeutic efficacy. High-resolution mass spectrometry (MS) has allowed initial insights into the antibody repertoire. We outline current challenges in MS-based sequence analysis of antibody repertoires and propose strategies for their resolution.

CD22 Controls Germinal Center B Cell Receptor Signaling, Which Influences Plasma Cell and Memory B Cell Output

Germinal center reactions are established during a thymus-dependent immune response. Germinal center (GC) B cells are rapidly proliferating and undergo somatic hypermutation in Ab genes. This results in the production of high-affinity Abs and establishment of long-lived memory cells. GC B cells show lower BCR-induced signaling when compared with naive B cells, but the functional relevance is not clear. CD22 is a member of the Siglec family and functions as an inhibitory coreceptor on B cells. Interestingly, GC B cells downregulate sialic acid forms that serve as high-affinity ligands for CD22, indicating a role for CD22 ligand binding during GC responses. We studied the role of CD22 in the GC with mixed bone marrow chimeric mice and found a disadvantage of CD22^-/- GC B cells during the GC reaction. Mechanistic investigations ruled out defects in dark zone/light zone distribution and affinity maturation. Rather, an increased rate of apoptosis in CD22^-/- GC B cells was responsible for the disadvantage, also leading to a lower GC output in plasma cells and memory B cells. CD22^-/- GC B cells showed a clearly increased calcium response upon BCR stimulation, which was almost absent in wild-type GC B cells. We conclude that the differential expression of the low-affinity cis CD22 ligands in the GC normally results in a strong attenuation of BCR signaling in GC B cells, probably due to higher CD22-BCR interactions. Therefore, attenuation of BCR signaling by CD22 is involved in GC output and B cell fate.

Cell-to-cell expression dispersion of B-cell surface proteins is linked to genetic variants in humans

Variability in gene expression across a population of homogeneous cells is known to influence various biological processes. In model organisms, natural genetic variants were found that modify expression dispersion (variability at a fixed mean) but very few studies have detected such effects in humans. Here, we analyzed single-cell expression of four proteins (CD23, CD55, CD63 and CD86) across cell lines derived from individuals of the Yoruba population. Using data from over 30 million cells, we found substantial inter-individual variation of dispersion. We demonstrate, via de novo cell line generation and subcloning experiments, that this variation exceeds the variation associated with cellular immortalization. We detected a genetic association between the expression dispersion of CD63 and the rs971 SNP. Our results show that human DNA variants can have inherently-probabilistic effects on gene expression. Such subtle genetic effects may participate to phenotypic variation and disease outcome.

Triqueneaux, Burny, Symmons et al. show association between gene expression noise and genotypes, using single-cell expression of four proteins across human-derived lymphoblastoid cell lines. This study suggests that very subtle regulatory effects of human DNA variants may contribute to phenotypic variation and disease outcome.

Differential Expression of IgM and IgD Discriminates Two Subpopulations of Human Circulating IgM+IgD+CD27+ B Cells That Differ Phenotypically, Functionally, and Genetically

The origin and function of blood IgM⁺IgD⁺CD27⁺ B cells is controversial, and they are considered a heterogeneous population. Previous staining of circulating B cells of healthy donors with rotavirus fluorescent virus-like particles allowed us to differentiate two subsets of IgM⁺IgD⁺CD27⁺: IgM^hi and IgM^lo B cells. Here, we confirmed this finding and compared the phenotype, transcriptome, in vitro function, and Ig gene repertoire of these two subsets. Eleven markers phenotypically discriminated both subsets (CD1c, CD69, IL21R, CD27, MTG, CD45RB, CD5, CD184, CD23, BAFFR, and CD38) with the IgM^hi phenotypically resembling previously reported marginal zone B cells and the IgM^lo resembling both naïve and memory B cells. Transcriptomic analysis showed that both subpopulations clustered close to germinal center-experienced IgM only B cells with a Principal Component Analysis, but differed in expression of 78 genes. Moreover, IgM^hi B cells expressed genes characteristic of previously reported marginal zone B cells. After stimulation with CpG and cytokines, significantly (p < 0.05) higher frequencies (62.5%) of IgM^hi B cells proliferated, compared with IgM^lo B cells (35.37%), and differentiated to antibody secreting cells (14.22% for IgM^hi and 7.19% for IgM^lo). IgM^hi B cells had significantly (p < 0.0007) higher frequencies of mutations in IGHV and IGKV regions, IgM^lo B cells had higher usage of IGHJ6 genes (p < 0.0001), and both subsets differed in their HCDR3 properties. IgM^hi B cells shared most of their shared IGH clonotypes with IgM only memory B cells, and IgM^lo B cells with IgM^hi B cells. These results support the notion that differential expression of IgM and IgD discriminates two subpopulations of human circulating IgM⁺IgD⁺CD27⁺ B cells, with the IgM^hi B cells having similarities with previously described marginal zone B cells that passed through germinal centers, and the IgM^lo B cells being the least differentiated amongst the IgM⁺CD27⁺ subsets.

Top-down Mass Spectrometry Analysis of Human Serum Autoantibody Antigen-Binding Fragments

Detecting autoimmune diseases at an early stage is crucial for effective treatment and disease management to slow disease progression and prevent irreversible organ damage. In many autoimmune diseases, disease-specific autoantibodies are produced by B cells in response to soluble autoantigens due to defects in B cell tolerance mechanisms. Autoantibodies accrue early in disease development, and several are so disease-specific they serve as classification criteria. In this study, we established a high-throughput, sensitive, intact serum autoantibody analysis platform based on the optimization of a one dimensional ultra-high-pressure liquid chromatography top-down mass spectrometry platform (1D UPLC-TDMS). This approach has been successfully applied to a 12 standard monoclonal antibody antigen-binding fragment (Fab) mixture, demonstrating the feasibility to separate and sequence intact antibodies with high sequence coverage and high sensitivity. We then applied the optimized platform to characterize total serum antibody Fabs in a systemic lupus erythematosus (SLE) patient sample and compared it to healthy control samples. From this analysis, we show that the SLE sample has many dominant antibody Fab-related mass features unlike the healthy controls. To our knowledge, this is the first top-down demonstration of serum autoantibody pool analysis. Our proposed approach holds great promise for discovering novel serum autoantibody biomarkers that are of interest for diagnosis, prognosis, and tolerance induction, as well as improving our understanding of pathogenic autoimmune processes.

Meta-Analysis for Correlating Structure of Bioactive Peptides in Foods of Animal Origin with Regard to Effect and Stability

Amino acid (AA) sequences of 807 bioactive peptides from foods of animal origin were examined in order to correlate peptide structure with activity (antihypertensive, antioxidative, immunomodulatory, antimicrobial, hypolipidemic, antithrombotic, and opioid) and stability in vivo. Food sources, such as milk, meat, eggs, and marine products, show different frequencies of bioactive peptides exhibiting specific effects. There is a correlation of peptide structure and effect, depending on type and position of AA. Opioid peptides contain a high percentage of aromatic AA residues, while antimicrobial peptides show an excess of positively charged AAs. AA residue position is significant, with those in the first and penultimate positions having the biggest effects on peptide activity. Peptides that have activity in vivo contain a high percentage (67%) of proline residues, but the positions of proline in the sequence depend on the length of the peptide. We also discuss the influence of processing on activity of these peptides, as well as methods for predicting release from the source protein and activity of peptides.

A Two-State Model of Tree Evolution and Its Applications to Alu Retrotransposition

Models of tree evolution have mostly focused on capturing the cladogenesis processes behind speciation. Processes that derive the evolution of genomic elements, such as repeats, are not necessarily captured by these existing models. In this article, we design a model of tree evolution that we call the dual-birth model, and we show how it can be useful in studying the evolution of short Alu repeats found in the human genome in abundance. The dual-birth model extends the traditional birth-only model to have two rates of propagation, one for active nodes that propagate often, and another for inactive nodes, that with a lower rate, activate and start propagating. Adjusting the ratio of the rates controls the expected tree balance. We present several theoretical results under the dual-birth model, introduce parameter estimation techniques, and study the properties of the model in simulations. We then use the dual-birth model to estimate the number of active Alu elements and their rates of propagation and activation in the human genome based on a large phylogenetic tree that we build from close to one million Alu sequences.

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.

Immunoglobulin gene analysis as a tool for investigating human immune responses

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

Analyzing Immunoglobulin Repertoires

Somatic assembly of T cell receptor and B cell receptor (BCR) genes produces a vast diversity of lymphocyte antigen recognition capacity. The advent of efficient high-throughput sequencing of lymphocyte antigen receptor genes has recently generated unprecedented opportunities for exploration of adaptive immune responses. With these opportunities have come significant challenges in understanding the analysis techniques that most accurately reflect underlying biological phenomena. In this regard, sample preparation and sequence analysis techniques, which have largely been borrowed and adapted from other fields, continue to evolve. Here, we review current methods and challenges of library preparation, sequencing and statistical analysis of lymphocyte receptor repertoire studies. We discuss the general steps in the process of immune repertoire generation including sample preparation, platforms available for sequencing, processing of sequencing data, measurable features of the immune repertoire, and the statistical tools that can be used for analysis and interpretation of the data. Because BCR analysis harbors additional complexities, such as immunoglobulin (Ig) (i.e., antibody) gene somatic hypermutation and class switch recombination, the emphasis of this review is on Ig/BCR sequence analysis.

Studying Antibody Repertoires with Next-Generation Sequencing

Next-generation sequencing is making it possible to study the antibody repertoire of an organism in unprecedented detail, and, by so doing, to characterize its behavior in the response to infection and in pathological conditions such as autoimmunity and cancer. The polymorphic nature of the repertoire poses unique challenges that rule out the use of many commonly used NGS methods and require tradeoffs to be made when considering experimental design.We outline the main contexts in which antibody repertoire analysis has been used, and summarize the key tools that are available. The humoral immune response to vaccination has been a particular focus of repertoire analyses, and we review the key conclusions and methods used in these studies.

The Antibody Repertoire of Colorectal Cancer

Immunotherapy is becoming increasingly important in the fight against cancers, using and manipulating the body's immune response to treat tumors. Understanding the immune repertoire-the collection of immunological proteins-of treated and untreated cells is possible at the genomic, but technically difficult at the protein level. Standard protein databases do not include the highly divergent sequences of somatic rearranged immunoglobulin genes, and may lead to miss identifications in a mass spectrometry search. We introduce a novel proteogenomic approach, AbScan, to identify these highly variable antibody peptides, by developing a customized antibody database construction method using RNA-seq reads aligned to immunoglobulin (Ig) genes.AbScan starts by filtering transcript (RNA-seq) reads that match the template for Ig genes. The retained reads are used to construct a repertoire graph using the "split" de Bruijn graph: a graph structure that improves on the standard de Bruijn graph to capture the high diversity of Ig genes in a compact manner. AbScan corrects for sequencing errors, and converts the graph to a format suitable for searching with MS/MS search tools. We used AbScan to create an antibody database from 90 RNA-seq colorectal tumor samples. Next, we used proteogenomic analysis to search MS/MS spectra of matched colorectal samples from the Clinical Proteomic Tumor Analysis Consortium (CPTAC) against the AbScan generated database. AbScan identified 1,940 distinct antibody peptides. Correlating with previously identified Single Amino-Acid Variants (SAAVs) in the tumor samples, we identified 163 pairs (antibody peptide, SAAV) with significant cooccurrence pattern in the 90 samples. The presence of coexpressed antibody and mutated peptides was correlated with survival time of the individuals. Our results suggest that AbScan (https://github.com/csw407/AbScan.git) is an effective tool for a proteomic exploration of the immune response in cancers.

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.

Reconstructing Antibody Repertoires from Error-Prone Immunosequencing Reads

Transforming error-prone immunosequencing datasets into Ab repertoires is a fundamental problem in immunogenomics, and a prerequisite for studies of immune responses. Although various repertoire reconstruction algorithms were released in the last 3 y, it remains unclear how to benchmark them and how to assess the accuracy of the reconstructed repertoires. We describe an accurate IgReC algorithm for constructing Ab repertoires from high-throughput immunosequencing datasets and a new framework for assessing the quality of reconstructed repertoires. Surprisingly, Ab repertoires constructed by IgReC from barcoded immunosequencing datasets in the blind mode (without using information about unique molecular identifiers) improved upon the repertoires constructed by the state-of-the-art tools that use barcoding. This finding suggests that IgReC may alleviate the need to generate repertoires using the barcoding technology (the workhorse of current immunogenomics efforts) because our computational approach to error correction of immunosequencing data is nearly as powerful as the experimental approach based on barcoding.

Computational Model Reveals Limited Correlation between Germinal Center B-Cell Subclone Abundancy and Affinity: Implications for Repertoire Sequencing

Immunoglobulin repertoire sequencing has successfully been applied to identify expanded antigen-activated B-cell clones that play a role in the pathogenesis of immune disorders. One challenge is the selection of the Ag-specific B cells from the measured repertoire for downstream analyses. A general feature of an immune response is the expansion of specific clones resulting in a set of subclones with common ancestry varying in abundance and in the number of acquired somatic mutations. The expanded subclones are expected to have BCR affinities for the Ag higher than the affinities of the naive B cells in the background population. For these reasons, several groups successfully proceeded or suggested selecting highly abundant subclones from the repertoire to obtain the Ag-specific B cells. Given the nature of affinity maturation one would expect that abundant subclones are of high affinity but since repertoire sequencing only provides information about abundancies, this can only be verified with additional experiments, which are very labor intensive. Moreover, this would also require knowledge of the Ag, which is often not available for clinical samples. Consequently, in general we do not know if the selected highly abundant subclone(s) are also the high(est) affinity subclones. Such knowledge would likely improve the selection of relevant subclones for further characterization and Ag screening. Therefore, to gain insight in the relation between subclone abundancy and affinity, we developed a computational model that simulates affinity maturation in a single GC while tracking individual subclones in terms of abundancy and affinity. We show that the model correctly captures the overall GC dynamics, and that the amount of expansion is qualitatively comparable to expansion observed from B cells isolated from human lymph nodes. Analysis of the fraction of high- and low-affinity subclones among the unexpanded and expanded subclones reveals a limited correlation between abundancy and affinity and shows that the low abundant subclones are of highest affinity. Thus, our model suggests that selecting highly abundant subclones from repertoire sequencing experiments would not always lead to the high(est) affinity B cells. Consequently, additional or alternative selection approaches need to be applied.

De Novo MS/MS Sequencing of Native Human Antibodies

One direct route for the discovery of therapeutic human monoclonal antibodies (mAbs) involves the isolation of peripheral B cells from survivors/sero-positive individuals after exposure to an infectious reagent or disease etiology, followed by single-cell sequencing or hybridoma generation. Peripheral B cells, however, are not always easy to obtain and represent only a small percentage of the total B-cell population across all bodily tissues. Although it has been demonstrated that tandem mass spectrometry (MS/MS) techniques can interrogate the full polyclonal antibody (pAb) response to an antigen in vivo, all current approaches identify MS/MS spectra against databases derived from genetic sequencing of B cells from the same patient. In this proof-of-concept study, we demonstrate the feasibility of a novel MS/MS antibody discovery approach in which only serum antibodies are required without the need for sequencing of genetic material. Peripheral pAbs from a cytomegalovirus-exposed individual were purified by glycoprotein B antigen affinity and de novo sequenced from MS/MS data. Purely MS-derived mAbs were then manufactured in mammalian cells to validate potency via antigen-binding ELISA. Interestingly, we found that these mAbs accounted for 1 to 2% of total donor IgG but were not detected in parallel sequencing of memory B cells from the same patient.

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.

VDJML: a file format with tools for capturing the results of inferring immune receptor rearrangements

Background

The genes that produce antibodies and the immune receptors expressed on lymphocytes are not germline encoded; rather, they are somatically generated in each developing lymphocyte by a process called V(D)J recombination, which assembles specific, independent gene segments into mature composite genes. The full set of composite genes in an individual at a single point in time is referred to as the immune repertoire. V(D)J recombination is the distinguishing feature of adaptive immunity and enables effective immune responses against an essentially infinite array of antigens. Characterization of immune repertoires is critical in both basic research and clinical contexts. Recent technological advances in repertoire profiling via high-throughput sequencing have resulted in an explosion of research activity in the field. This has been accompanied by a proliferation of software tools for analysis of repertoire sequencing data. Despite the widespread use of immune repertoire profiling and analysis software, there is currently no standardized format for output files from V(D)J analysis. Researchers utilize software such as IgBLAST and IMGT/High V-QUEST to perform V(D)J analysis and infer the structure of germline rearrangements. However, each of these software tools produces results in a different file format, and can annotate the same result using different labels. These differences make it challenging for users to perform additional downstream analyses.

Results

To help address this problem, we propose a standardized file format for representing V(D)J analysis results. The proposed format, VDJML, provides a common standardized format for different V(D)J analysis applications to facilitate downstream processing of the results in an application-agnostic manner. The VDJML file format specification is accompanied by a support library, written in C++ and Python, for reading and writing the VDJML file format.

Conclusions

The VDJML suite will allow users to streamline their V(D)J analysis and facilitate the sharing of scientific knowledge within the community. The VDJML suite and documentation are available from https://vdjserver.org/vdjml/. We welcome participation from the community in developing the file format standard, as well as code contributions.

B-cell repertoire dynamics after sequential hepatitis B vaccination and evidence for cross-reactive B-cell activation

BACKGROUND

A diverse B-cell repertoire is essential for recognition and response to infectious and vaccine antigens. High-throughput sequencing of B-cell receptor (BCR) genes can now be used to study the B-cell repertoire at great depth and may shed more light on B-cell responses than conventional immunological methods. Here, we use high-throughput BCR sequencing to provide novel insight into B-cell dynamics following a primary course of hepatitis B vaccination.

METHODS

Nine vaccine-naïve participants were administered three doses of hepatitis B vaccine (months 0, 1, and 2 or 7). High-throughput Illumina sequencing of the total BCR repertoire was combined with targeted sequencing of sorted vaccine antigen-enriched B cells to analyze the longitudinal response of both the total and vaccine-specific repertoire after each vaccine. ELISpot was used to determine vaccine-specific cell numbers following each vaccine.

RESULTS

Deconvoluting the vaccine-specific from total BCR repertoire showed that vaccine-specific sequence clusters comprised <0.1 % of total sequence clusters, and had certain stereotypic features. The vaccine-specific BCR sequence clusters were expanded after each of the three vaccine doses, despite no vaccine-specific B cells being detected by ELISpot after the first vaccine dose. These vaccine-specific BCR clusters detected after the first vaccine dose had distinct properties compared to those detected after subsequent doses; they were more mutated, present at low frequency even prior to vaccination, and appeared to be derived from more mature B cells.

CONCLUSIONS

These results demonstrate the high-sensitivity of our vaccine-specific BCR analysis approach and suggest an alternative view of the B-cell response to novel antigens. In the response to the first vaccine dose, many vaccine-specific BCR clusters appeared to largely derive from previously activated cross-reactive B cells that have low affinity for the vaccine antigen, and subsequent doses were required to yield higher affinity B cells.

Immunoglobulin Classification Using the Colored Antibody Graph

The somatic recombination of V, D, and J gene segments in B-cells introduces a great deal of diversity, and divergence from reference segments. Many recent studies of antibodies focus on the population of antibody transcripts that show which V, D, and J gene segments have been favored for a particular antigen, a repertoire. To properly describe the antibody repertoire, each antibody must be labeled by its constituting V, D, and J gene segment, a task made difficult by somatic recombination and hypermutation events. While previous approaches to repertoire analysis were based on sequential alignments, we describe a new de Bruijn graph-based algorithm to perform VDJ labeling and benchmark its performance.

Secreted autoantibody repertoires in Sjögren's syndrome and systemic lupus erythematosus: A proteomic approach

The structures of epitopes bound by autoantibodies against RNA-protein complexes have been well-defined over several decades, but little is known of the clonality, immunoglobulin (Ig) variable (V) gene usage and mutational status of the autoantibodies themselves at the level of the secreted (serum) proteome. A novel proteomic workflow is presented based on affinity purification of specific Igs from serum, high-resolution two-dimensional gel electrophoresis, and de novo and database-driven sequencing of V-region proteins by mass spectrometry. Analysis of anti-Ro52/Ro60/La proteomes in primary Sjögren's syndrome (SS) and anti-Sm and anti-ribosomal P proteomes in systemic lupus erythematosus (SLE) has revealed that these antibody responses are dominated by restricted sets of public (shared) clonotypes, consistent with common pathways of production across unrelated individuals. The discovery of shared sets of specific V-region peptides can be exploited for diagnostic biomarkers in targeted mass spectrometry platforms and for tracking and removal of pathogenic clones.