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

Human IgMhiCD300a+ B Cells Are Circulating Marginal Zone Memory B Cells That Respond to Pneumococcal Polysaccharides and Their Frequency Is Decreased in People Living with HIV

CD300a is differentially expressed among B cell subsets, although its expression in immunoglobulin (Ig)M+ B cells is not well known. We identified a B cell subset expressing CD300a and high levels of IgM (IgMhiCD300a+). The results showed that IgMhiCD300a+ B cells were CD10-CD27+CD25+IgDloCD21hiCD23-CD38loCD1chi, suggesting that they are circulating marginal zone (MZ) IgM memory B cells. Regarding the immunoglobulin repertoire, IgMhiCD300a+ B cells exhibited a higher mutation rate and usage of the IgH-VDJ genes than the IgM+CD300a- counterpart. Moreover, the shorter complementarity-determining region 3 (CDR3) amino acid (AA) length from IgMhiCD300a+ B cells together with the predicted antigen experience repertoire indicates that this B cell subset has a memory phenotype. IgM memory B cells are important in T cell-independent responses. Accordingly, we demonstrate that this particular subset secretes higher amounts of IgM after stimulation with pneumococcal polysaccharides or a toll-like receptor 9 (TLR9) agonist than IgM+CD300a- cells. Finally, the frequency of IgMhiCD300a+ B cells was lower in people living with HIV-1 (PLWH) and it was inversely correlated with the years with HIV infection. Altogether, these data help to identify a memory B cell subset that contributes to T cell-independent responses to pneumococcal infections and may explain the increase in severe pneumococcal infections and the impaired responses to pneumococcal vaccination in PLWH.

Age-dependent nasal immune responses in non-hospitalized bronchiolitis children

Bronchiolitis in children is associated with significant rates of morbidity and mortality. Many studies have been performed using samples from hospitalized bronchiolitis patients, but little is known about the immunological responses from infants suffering from mild/moderate bronchiolitis that do not require hospitalization. We have studied a collection of nasal lavage fluid (NLF) samples from outpatient bronchiolitis children as a novel strategy to unravel local humoral and cellular responses, which are not fully characterized. The children were age-stratified in three groups, two of them (GI under 2-months, GII between 2-4 months) presenting a first episode of bronchiolitis, and GIII (between 4 months and 2 years) with recurrent respiratory infections. Here we show that elevated levels of pro-inflammatory cytokines (IL1β, IL6, TNFα, IL18, IL23), regulatory cytokines (IL10, IL17A) and IFNγ were found in the three bronchiolitis cohorts. However, little or no change was observed for IL33 and MCP1, at difference to previous results from bronchiolitis hospitalized patients. Furthermore, our results show a tendency to IL1β, IL6, IL18 and TNFα increased levels in children with mild pattern of symptom severity and in those in which non RSV respiratory virus were detected compared to RSV+ samples. By contrast, no such differences were found based on gender distribution. Bronchiolitis NLFs contained more IgM, IgG1, IgG3 IgG4 and IgA than NLF from their age-matched healthy controls. NLF from bronchiolitis children predominantly contained neutrophils, and also low frequency of monocytes and few CD4+ and CD8+ T cells. NLF from infants older than 4-months contained more intermediate monocytes and B cell subsets, including naïve and memory cells. BCR repertoire analysis of NLF samples showed a biased VH1 usage in IgM repertoires, with low levels of somatic hypermutation. Strikingly, algorithmic studies of the mutation profiles, denoted antigenic selection on IgA-NLF repertoires. Our results support the use of NLF samples to analyze immune responses and may have therapeutic implications.

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

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

Immune Repertoire Analysis on High-Performance Computing Using VDJServer V1: A Method by the AIRR Community

AIRR-seq data sets are usually large and require specialized analysis methods and software tools. A typical Illumina MiSeq sequencing run generates 20-30 million 2 × 300 bp paired-end sequence reads, which roughly corresponds to 15 GB of sequence data to be processed. Other platforms like NextSeq, which is useful in projects where the full V gene is not needed, create about 400 million 2 × 150 bp paired-end reads. Because of the size of the data sets, the analysis can be computationally expensive, particularly the early analysis steps like preprocessing and gene annotation that process the majority of the sequence data. A standard desktop PC may take 3-5 days of constant processing for a single MiSeq run, so dedicated high-performance computational resources may be required.VDJServer provides free access to high-performance computing (HPC) at the Texas Advanced Computing Center (TACC) through a graphical user interface (Christley et al. Front Immunol 9:976, 2018). VDJServer is a cloud-based analysis portal for immune repertoire sequence data that provides access to a suite of tools for a complete analysis workflow, including modules for preprocessing and quality control of sequence reads, V(D)J gene assignment, repertoire characterization, and repertoire comparison. Furthermore, VDJServer has parallelized execution for tools such as IgBLAST, so more compute resources are utilized as the size of the input data grows. Analysis that takes days on a desktop PC might take only a few hours on VDJServer. VDJServer is a free, publicly available, and open-source licensed resource. Here, we describe the workflow for performing immune repertoire analysis on VDJServer's high-performance computing.

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.

The methods and advances of adaptive immune receptors repertoire sequencing

The adaptive immune response is a powerful tool, capable of recognizing, binding to, and neutralizing a vast number of internal and external threats via T or B lymphatic receptors with widespread sets of antigen specificities. The emergence of high-throughput sequencing technology and bioinformatics provides opportunities for research in the fields of life sciences and medicine. The analysis and annotation for immune repertoire data can reveal biologically meaningful information, including immune prediction, target antigens, and effective evaluation. Continuous improvements of the immunological repertoire sequencing methods and analysis tools will help to minimize the experimental and calculation errors and realize the immunological information to meet the clinical requirements. That said, the clinical application of adaptive immune repertoire sequencing requires appropriate experimental methods and standard analytical tools. At the population cell level, we can acquire the overview of cell groups, but the information about a single cell is not obtained accurately. The information that is ignored may be crucial for understanding the heterogeneity of each cell, gene expression and drug response. The combination of high-throughput sequencing and single-cell technology allows us to obtain single-cell information with low-cost and high-throughput. In this review, we summarized the current methods and progress in this area.

The Identity Card of T Cells-Clinical Utility of T-cell Receptor Repertoire Analysis in Transplantation

There is a clear medical need to change the current strategy of "one-size-fits-all" immunosuppression for controlling transplant rejection to precision medicine and targeted immune intervention. As T cells play a key role in both undesired graft rejection and protection, a better understanding of the fate and function of both alloreactive graft-deteriorating T cells and those protecting to infections is required. The T-cell receptor (TCR) is the individual identity card of each T cell clone and can help to follow single specificities. In this context, tracking of lymphocytes with certain specificity in blood and tissue in clinical follow up is of especial importance. After overcoming technical limitations of the past, novel molecular technologies opened new avenues of diagnostics. Using advantages of next generation sequencing, a method was established for T-cell tracing by detection of variable TCR region as identifiers of individual lymphocyte clones. The current review describes principles of laboratory and computational methods of TCR repertoire analysis, and gives an overview on applications for the basic understanding of transplant biology and immune monitoring. The review also delineates methodological pitfalls and challenges. With the outlook on prediction of antigens in immune-mediated processes including those of unknown causative pathogens, monitoring the fate and function of individual T cell clones, and the adoptive transfer of protective effector or regulatory T cells, this review highlights the current and future capability of TCR repertoire analysis.

The Pipeline Repertoire for Ig-Seq Analysis

With the advent of high-throughput sequencing of immunoglobulin genes (Ig-Seq), the understanding of antibody repertoires and their dynamics among individuals and populations has become an exciting area of research. There is an increasing number of computational tools that aid in every step of the immune repertoire characterization. However, since not all tools function identically, every pipeline has its unique rationale and capabilities, creating a rich blend of useful features that may appear intimidating for newcomer laboratories with the desire to plunge into immune repertoire analysis to expand and improve their research; hence, all pipeline strengths and differences may not seem evident. In this review we provide a practical and organized list of the current set of computational tools, focusing on their most attractive features and differences in order to carry out the characterization of antibody repertoires so that the reader better decides a strategic approach for the experimental design, and computational pathways for the analyses of immune repertoires.

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.

Efficient Selection of Antibodies Reactive to Homologous Epitopes on Human and Mouse Hepatocyte Growth Factors by Next-Generation Sequencing-Based Analysis of the B Cell Repertoire

: YYB-101 is a humanized rabbit anti-human hepatocyte growth factor (HGF)-neutralizing antibody currently in clinical trial. To test the effect of HGF neutralization with antibody on anti-cancer T cell immunity, we generated surrogate antibodies that are reactive to the mouse homologue of the epitope targeted by YYB-101. First, we immunized a chicken with human HGF and monitored changes in the B cell repertoire by next-generation sequencing (NGS). We then extracted the VH gene repertoire from the NGS data, clustered it into components by sequence homology, and classified the components by the change in the number of unique VH sequences and the frequencies of the VH sequences within each component following immunization. Those changes should accompany the preferential proliferation and somatic hypermutation or gene conversion of B cells encoding HGF-reactive antibodies. One component showed significant increases in the number and frequencies of unique VH sequences and harbored genes encoding antibodies that were reactive to human HGF and competitive with YYB-101 for HGF binding. Some of the antibodies also reacted to mouse HGF. The selected VH sequences shared 98.3% identity and 98.9% amino acid similarity. It is therefore likely that the antibodies encoded by them all react to the epitope targeted by YYB-101.

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

Background

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

Results

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

Conclusion

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

Reproducibility and Reuse of Adaptive Immune Receptor Repertoire Data

High-throughput sequencing (HTS) of immunoglobulin (B-cell receptor, antibody) and T-cell receptor repertoires has increased dramatically since the technique was introduced in 2009 (1–3). This experimental approach explores the maturation of the adaptive immune system and its response to antigens, pathogens, and disease conditions in exquisite detail. It holds significant promise for diagnostic and therapy-guiding applications. New technology often spreads rapidly, sometimes more rapidly than the understanding of how to make the products of that technology reliable, reproducible, or usable by others. As complex technologies have developed, scientific communities have come together to adopt common standards, protocols, and policies for generating and sharing data sets, such as the MIAME protocols developed for microarray experiments. The Adaptive Immune Receptor Repertoire (AIRR) Community formed in 2015 to address similar issues for HTS data of immune repertoires. The purpose of this perspective is to provide an overview of the AIRR Community’s founding principles and present the progress that the AIRR Community has made in developing standards of practice and data sharing protocols. Finally, and most important, we invite all interested parties to join this effort to facilitate sharing and use of these powerful data sets (join@airr-community.org).

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.

T Cell Receptor Profiling in Type 1 Diabetes

PURPOSE OF REVIEW

The genetic susceptibility and dominant protection for type 1 diabetes (T1D) associated with human leukocyte antigen (HLA) haplotypes, along with minor risk variants, have long been thought to shape the T cell receptor (TCR) repertoire and eventual phenotype of autoreactive T cells that mediate β-cell destruction. While autoantibodies provide robust markers of disease progression, early studies tracking autoreactive T cells largely failed to achieve clinical utility.

RECENT FINDINGS

Advances in acquisition of pancreata and islets from T1D organ donors have facilitated studies of T cells isolated from the target tissues. Immunosequencing of TCR α/β-chain complementarity determining regions, along with transcriptional profiling, offers the potential to transform biomarker discovery. Herein, we review recent studies characterizing the autoreactive TCR signature in T1D, emerging technologies, and the challenges and opportunities associated with tracking TCR molecular profiles during the natural history of T1D.

Reproducibility and Reuse of Adaptive Immune Receptor Repertoire Data

High-throughput sequencing (HTS) of immunoglobulin (B-cell receptor, antibody) and T-cell receptor repertoires has increased dramatically since the technique was introduced in 2009 (1-3). This experimental approach explores the maturation of the adaptive immune system and its response to antigens, pathogens, and disease conditions in exquisite detail. It holds significant promise for diagnostic and therapy-guiding applications. New technology often spreads rapidly, sometimes more rapidly than the understanding of how to make the products of that technology reliable, reproducible, or usable by others. As complex technologies have developed, scientific communities have come together to adopt common standards, protocols, and policies for generating and sharing data sets, such as the MIAME protocols developed for microarray experiments. The Adaptive Immune Receptor Repertoire (AIRR) Community formed in 2015 to address similar issues for HTS data of immune repertoires. The purpose of this perspective is to provide an overview of the AIRR Community's founding principles and present the progress that the AIRR Community has made in developing standards of practice and data sharing protocols. Finally, and most important, we invite all interested parties to join this effort to facilitate sharing and use of these powerful data sets (join@airr-community.org).