Advances and applications of immune receptor sequencing in systems immunology

Abalign: a comprehensive multiple sequence alignment platform for B-cell receptor immune repertoires

The utilization of high-throughput sequencing (HTS) for B-cell receptor (BCR) immune repertoire analysis has become widespread in the fields of adaptive immunity and antibody drug development. However, the sheer volume of sequences generated by these experiments presents a challenge in data processing. Specifically, multiple sequence alignment (MSA), a critical aspect of BCR analysis, remains inadequate for handling massive BCR sequencing data and lacks the ability to provide immunoglobulin-specific information. To address this gap, we introduce Abalign, a standalone program specifically designed for ultrafast MSA of BCR/antibody sequences. Benchmark tests demonstrate that Abalign achieves comparable or even better accuracy than state-of-the-art MSA tools, and shows remarkable advantages in terms of speed and memory consumption, reducing the time required for high-throughput analysis from weeks to hours. In addition to its alignment capabilities, Abalign offers a broad range of BCR analysis features, including extracting BCRs, constructing lineage trees, assigning VJ genes, analyzing clonotypes, profiling mutations, and comparing BCR immune repertoires. With its user-friendly graphic interface, Abalign can be easily run on personal computers instead of computing clusters. Overall, Abalign is an easy-to-use and effective tool that enables researchers to analyze massive BCR/antibody sequences, leading to new discoveries in the field of immunoinformatics. The software is freely available at http://cao.labshare.cn/abalign/.

The influences of microbial colonisation and germ-free status on the chicken TCRβ repertoire

Microbial colonisation is paramount to the normal development of the immune system, particularly at mucosal sites. However, the relationships between the microbiome and the adaptive immune repertoire have mostly been explored in rodents and humans. Here, we report a high-throughput sequencing analysis of the chicken TCRβ repertoire and the influences of microbial colonisation on tissue-resident TCRβ+ cells. The results reveal that the microbiome is an important driver of TCRβ diversity in both intestinal tissues and the bursa of Fabricius, but not in the spleen. Of note, public TCRβ sequences (shared across individuals) make a substantial contribution to the repertoire. Additionally, different tissues exhibit biases in terms of their V family and J gene usage, and these effects were influenced by the gut-associated microbiome. TCRβ clonal expansions were identified in both colonised and germ-free birds, but differences between the groups were indicative of an influence of the microbiota. Together, these findings provide an insight into the avian adaptive immune system and the influence of the microbiota on the TCRβ repertoire.

T-cell receptor variable region usage in Chagas disease: A systematic review of experimental and human studies

T cells recognize their ligand, the peptide major histocompatibility complex (MHC), via the T-cell receptor (TCR), which is composed of covalently linked α and β or γ and δ chains. This recognition is critical for T-cell ontogeny and controls the selection, activation, and function of T lymphocytes. Specific TCR αβ variable regions have been associated with immunopathogenesis of Chagas disease. Here, we present a systematic review that compiles experimental in vivo and human data regarding the preferential expression of variable alpha (Vα) and variable beta (Vβ) chain regions in Trypanosoma cruzi infection. The original studies indexed in PubMed/Medline, Scopus, and Web of Science databases were screened according to the PRISMA strategy. The analysis showed that expression of TCR Vα subfamilies were evaluated in one human study, and, unlike TCR Vβ, TCR Vα presented a more restricted usage. Despite the great variability in the usage of TCR Vβ regions in human Chagas disease, a down-regulation of TCR Vβ5 expression by T cells from patients in the acute phase of the disease was shown. Opposingly, this TCR region was found overly expressed in CD4+ T cells from chronic Chagas patients. It was also demonstrated that murine Vβ9+ T cells derived from nonlymphoid organs of T. cruzi-infected animals had a modulatory profile, while splenic Vβ9+ T cells produced inflammatory cytokines, indicating that although they display the same TCR Vβ region usage, these cells are functionally distinct. Despite the limitations of few papers and year of publication of the studies, compiling the data derived from them reveals that further investigation of TCR usage will point to their potential role in protective or pathogenic responses, as biomarkers of disease progression, and in the search for dominant peptides potentially useful for the development of vaccines or therapies.

Chagas disease is a neglected tropical disease, caused by infection with Trypanosoma cruzi. Differential expression of certain T-cell receptor (TCR) variable regions has been associated with the immunopathogenesis of Chagas disease. Here, we present a systematic review that compiled experimental in vivo and human data regarding the preferential expression of TCR alpha and beta chain variable regions in Chagas disease. The original studies indexed in the PubMed/Medline, Scopus, and Web of Science databases were screened according to the PRISMA strategy. Despite the great variability in the use of TCR Vβ in T. cruzi infection, the outcomes indicate that there is a down-regulation of TCR Vβ5 expression in T cells from patients in the acute phase of Chagas disease. However, this region is preferentially expressed by CD4+ T cells from chronic Chagas patients. Additionally, it has been demonstrated that murine Vβ9+ T cells derived from nonlymphoid organs displayed a modulatory profile, while splenic Vβ9+ T cells produced inflammatory cytokines, indicating that although they express the same TCR Vβ region, these cells are functionally distinct. Information on TCR expression, specificity and function have critical impact on vaccine design.

Different B cell subpopulations show distinct patterns in their IgH repertoire metrics

Several human B cell subpopulations are recognised in the peripheral blood, which play distinct roles in the humoral immune response. These cells undergo developmental and maturational changes involving VDJ recombination, somatic hypermutation and class switch recombination, altogether shaping their immunoglobulin heavy chain (IgH) repertoire. Here, we sequenced the IgH repertoire of naïve, marginal zone, switched and plasma cells from 10 healthy adults along with matched unsorted and in silico separated CD19⁺ bulk B cells. Using advanced bioinformatic analysis and machine learning, we show that sorted B cell subpopulations are characterised by distinct repertoire characteristics on both the individual sequence and the repertoire level. Sorted subpopulations shared similar repertoire characteristics with their corresponding in silico separated subsets. Furthermore, certain IgH repertoire characteristics correlated with the position of the constant region on the IgH locus. Overall, this study provides unprecedented insight over mechanisms of B cell repertoire control in peripherally circulating B cell subpopulations.

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.

Population variability in the generation and selection of T-cell repertoires

The diversity of T-cell receptor (TCR) repertoires is achieved by a combination of two intrinsically stochastic steps: random receptor generation by VDJ recombination, and selection based on the recognition of random self-peptides presented on the major histocompatibility complex. These processes lead to a large receptor variability within and between individuals. However, the characterization of the variability is hampered by the limited size of the sampled repertoires. We introduce a new software tool SONIA to facilitate inference of individual-specific computational models for the generation and selection of the TCR beta chain (TRB) from sequenced repertoires of 651 individuals, separating and quantifying the variability of the two processes of generation and selection in the population. We find not only that most of the variability is driven by the VDJ generation process, but there is a large degree of consistency between individuals with the inter-individual variance of repertoires being about ∼2% of the intra-individual variance. Known viral-specific TCRs follow the same generation and selection statistics as all TCRs.

OLGA: fast computation of generation probabilities of B- and T-cell receptor amino acid sequences and motifs

MOTIVATION

High-throughput sequencing of large immune repertoires has enabled the development of methods to predict the probability of generation by V(D)J recombination of T- and B-cell receptors of any specific nucleotide sequence. These generation probabilities are very non-homogeneous, ranging over 20 orders of magnitude in real repertoires. Since the function of a receptor really depends on its protein sequence, it is important to be able to predict this probability of generation at the amino acid level. However, brute-force summation over all the nucleotide sequences with the correct amino acid translation is computationally intractable. The purpose of this paper is to present a solution to this problem.

RESULTS

We use dynamic programming to construct an efficient and flexible algorithm, called OLGA (Optimized Likelihood estimate of immunoGlobulin Amino-acid sequences), for calculating the probability of generating a given CDR3 amino acid sequence or motif, with or without V/J restriction, as a result of V(D)J recombination in B or T cells. We apply it to databases of epitope-specific T-cell receptors to evaluate the probability that a typical human subject will possess T cells responsive to specific disease-associated epitopes. The model prediction shows an excellent agreement with published data. We suggest that OLGA may be a useful tool to guide vaccine design.

AVAILABILITY AND IMPLEMENTATION

Source code is available at https://github.com/zsethna/OLGA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

B-cell receptor repertoire sequencing in patients with primary immunodeficiency: a review

The advent of next-generation sequencing (NGS) now allows a detailed assessment of the adaptive immune system in health and disease. In particular, high-throughput B-cell receptor (BCR) repertoire sequencing provides detailed information about the functionality and abnormalities of the B-cell system. However, it is mostly unknown how the BCR repertoire is altered in the context of primary immunodeficiencies (PID) and whether findings are consistent throughout phenotypes and genotypes. We have performed an extensive literature search of the published work on BCR repertoire sequencing in PID patients, including several forms of predominantly antibody disorders and combined immunodeficiencies. It is somewhat surprising that BCR repertoires, even from severe clinical phenotypes, often show only mild abnormalities and that diversity or immunoglobulin gene segment usage is generally preserved to some extent. Despite the great variety of wet laboratory and analytical methods that were used in the different studies, several findings are common to most investigated PIDs, such as the increased usage of gene segments that are associated with self-reactivity. These findings suggest that BCR repertoire characteristics may be used to assess the functionality of the B-cell compartment irrespective of the underlying defect. With the use of NGS approaches, there is now the opportunity to apply BCR repertoire sequencing to multiple patients and explore the PID BCR repertoire in more detail. Ultimately, using BCR repertoire sequencing in translational research could aid the management of PID patients by improving diagnosis, estimating functionality of the immune system and improving assessment of prognosis.

Informatics for cancer immunotherapy

The rapid development of immunomodulatory cancer therapies has led to a concurrent increase in the application of informatics techniques to the analysis of tumors, the tumor microenvironment, and measures of systemic immunity. In this review, the use of tumors to gather genetic and expression data will first be explored. Next, techniques to assess tumor immunity are reviewed, including HLA status, predicted neoantigens, immune microenvironment deconvolution, and T-cell receptor sequencing. Attempts to integrate these data are in early stages of development and are discussed in this review. Finally, we review the application of these informatics strategies to therapy development, with a focus on vaccines, adoptive cell transfer, and checkpoint blockade therapies.