High-throughput sequencing of the T-cell receptor repertoire: pitfalls and opportunities

Deciphering Abdominal Aortic Diseases Through T-Cell Clonal Repertoire of Perivascular Adipose Tissue

BACKGROUND

Recent studies suggest that immune-mediated inflammation of perivascular adipose tissue of abdominal aortic aneurysms (AAAs) contributes to disease development and progression. Whether the perivascular adipose tissue of AAA is characterized by a specific adaptive immune signature remains unknown.

METHODS AND RESULTS

To investigate this hypothesis, we sequenced the T-cell receptor β-chain in the perivascular adipose tissue of patients with AAA and compared it with patients with aortic occlusive disease, who share the former anatomical site of the lesion and risk factors but differ in pathogenic mechanisms. Our results demonstrate that patients with AAA have a lower repertoire diversity than those with aortic occlusive disease and significant differences in variable/joining gene segment usage. Furthermore, we identified a set of 7 public T-cell receptor β-chain clonotypes that distinguished AAA and aortic occlusive disease with very high accuracy. We also found that the T-cell receptor β-chain repertoire differentially characterizes small and large AAAs (aortic diameter<55 mm and ≥55 mm, respectively).

CONCLUSIONS

This work supports the hypothesis that T cell-mediated immunity is fundamental in AAA pathogenesis and opens up new clinical perspectives.

T-cell immunity induced and reshaped by an anti-HPV immuno-oncotherapeutic lentiviral vector

We recently developed an immuno-oncotherapy against human papillomavirus (HPV)-induced tumors based on a lentiviral vector encoding the Early E6 and E7 oncoproteins of HPV16 and HPV18 genotypes, namely "Lenti-HPV-07". The robust and long-lasting anti-tumor efficacy of Lenti-HPV-07 is dependent on CD8+ T-cell induction and remodeling of the tumor microenvironment. Here, we first established that anti-vector immunity induced by Lenti-HPV-07 prime has no impact on the efficacy of a homologous boost to amplify anti-HPV T-cell immunity. To longitudinally monitor the evolution of the T-cell repertoire generated after the prime, homologous or heterologous boost with Lenti-HPV-07, we tracked T-cell clonotypes by deep sequencing of T-Cell Receptor (TCR) variable β and α chain mRNA, applied to whole peripheral blood cells (PBL) and a T cell population specific of an immunodominant E7HPV16 epitope. We observed a hyper-expansion of clonotypes post prime, accompanied by increased frequencies of HPV-07-specific T cells. Additionally, there was a notable diversification of clonotypes post boost in whole PBL, but not in the E7HPV16-specific T cells. We then demonstrated that the effector functions of such Lenti-HPV-07-induced T cells synergize with anti-checkpoint inhibitory treatments by systemic administration of anti-TIM3 or anti-NKG2A monoclonal antibodies. While Lenti-HPV-07 is about to enter a Phase I/IIa clinical trial, these results will help better elucidate its mode of action in immunotherapy against established HPV-mediated malignancies.

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.

Are homeostatic mechanisms aiding the reconstitution of the T-cell pool during lymphopenia in humans?

A timely recovery of T-cell numbers following haematopoietic stem-cell transplantation (HSCT) is essential for preventing complications, such as increased risk of infection and disease relapse. In analogy to the occurrence of lymphopenia-induced proliferation in mice, T-cell dynamics in humans are thought to be homeostatically regulated in a cell density-dependent manner. The idea is that T cells divide faster and/or live longer when T-cell numbers are low, thereby helping the reconstitution of the T-cell pool. T-cell reconstitution after HSCT is, however, known to occur notoriously slowly. In fact, the evidence for the existence of homeostatic mechanisms in humans is quite ambiguous, since lymphopenia is often associated with infectious complications and immune activation, which confound the study of homeostatic regulation. This calls into question whether homeostatic mechanisms aid the reconstitution of the T-cell pool during lymphopenia in humans. Here we review the changes in T-cell dynamics in different situations of T-cell deficiency in humans, including the early development of the immune system after birth, healthy ageing, HIV infection, thymectomy and hematopoietic stem cell transplantation (HSCT). We discuss to what extent these changes in T-cell dynamics are a side-effect of increased immune activation during lymphopenia, and to what extent they truly reflect homeostatic mechanisms.

Advances in antibody discovery from human BCR repertoires

Antibodies make up an important and growing class of compounds used for the diagnosis or treatment of disease. While traditional antibody discovery utilized immunization of animals to generate lead compounds, technological innovations have made it possible to search for antibodies targeting a given antigen within the repertoires of B cells in humans. Here we group these innovations into four broad categories: cell sorting allows the collection of cells enriched in specificity to one or more antigens; BCR sequencing can be performed on bulk mRNA, genomic DNA or on paired (heavy-light) mRNA; BCR repertoire analysis generally involves clustering BCRs into specificity groups or more in-depth modeling of antibody-antigen interactions, such as antibody-specific epitope predictions; validation of antibody-antigen interactions requires expression of antibodies, followed by antigen binding assays or epitope mapping. Together with innovations in Deep learning these technologies will contribute to the future discovery of diagnostic and therapeutic antibodies directly from humans.

NoisET: Noise Learning and Expansion Detection of T-Cell Receptors

High-throughput sequencing of T- and B-cell receptors makes it possible to track immune repertoires across time, in different tissues, in acute and chronic diseases and in healthy individuals. However, quantitative comparison between repertoires is confounded by variability in the read count of each receptor clonotype due to sampling, library preparation, and expression noise. We review methods for accounting for both biological and experimental noise and present an easy-to-use python package NoisET that implements and generalizes a previously developed Bayesian method. It can be used to learn experimental noise models for repertoire sequencing from replicates, and to detect responding clones following a stimulus. We test the package on different repertoire sequencing technologies and data sets. We review how such approaches have been used to identify responding clonotypes in vaccination and disease data. Availability: NoisET is freely available to use with source code at github.com/statbiophys/NoisET.

A hierarchy of selection pressures determines the organization of the T cell receptor repertoire

We systematically examine the receptor repertoire in T cell subsets in young, adult, and LCMV-infected mice. Somatic recombination generates diversity, resulting in the limited overlap between nucleotide sequences of different repertoires even within the same individual. However, statistical features of the repertoire, quantified by the V gene and CDR3 k-mer frequency distributions, are highly conserved. A hierarchy of immunological processes drives the evolution of this structure. Intra-thymic divergence of CD4+ and CD8+ lineages imposes subtle but dominant differences observed across repertoires of all subpopulations in both young and adult mice. Differentiation from naive through memory to effector phenotype imposes an additional gradient of repertoire diversification, which is further influenced by age in a complex and lineage-dependent manner. The distinct repertoire of CD4+ regulatory T cells is more similar to naive cells in young mice and to effectors in adults. Finally, we describe divergent (naive and memory) and convergent (CD8+ effector) evolution of the repertoire following acute infection with LCMV. This study presents a quantitative framework that captures the structure of the repertoire in terms of its fundamental statistical properties and describes how this structure evolves as individual T cells differentiate, migrate and mature in response to antigen exposure.

Haematopoietic Stem Cell Transplantation Results in Extensive Remodelling of the Clonal T Cell Repertoire in Multiple Sclerosis

Autologous haematopoietic stem cell transplantation (AHSCT) is a vital therapeutic option for patients with highly active multiple sclerosis (MS). Rates of remission suggest AHSCT is the most effective form of immunotherapy in controlling the disease. Despite an evolving understanding of the biology of immune reconstitution following AHSCT, the mechanism by which AHSCT enables sustained disease remission beyond the period of lymphopenia remains to be elucidated. Auto-reactive T cells are considered central to MS pathogenesis. Here, we analyse T cell reconstitution for 36 months following AHSCT in a cohort of highly active MS patients. Through longitudinal analysis of sorted naïve and memory T cell clones, we establish that AHSCT induces profound changes in the dominant T cell landscape of both CD4+ and CD8+ memory T cell clones. Lymphopenia induced homeostatic proliferation is followed by clonal attrition; with only 19% of dominant CD4 (p <0.025) and 13% of dominant CD8 (p <0.005) clones from the pre-transplant repertoire detected at 36 months. Recovery of a thymically-derived CD4 naïve T cell repertoire occurs at 12 months and is ongoing at 36 months, however diversity of the naïve populations is not increased from baseline suggesting the principal mechanism of durable remission from MS after AHSCT relates to depletion of putative auto-reactive clones. In a cohort of MS patients expressing the MS risk allele HLA DRB1*15:01, public clones are probed as potential biomarkers of disease. AHSCT appears to induce sustained periods of disease remission with dynamic changes in the clonal T cell repertoire out to 36 months post-transplant.

The importance of taking ART appropriately in children and adolescents with HIV-1 to reach the highest capacity of immune function later in life

Current antiretroviral therapy (ART) guidelines recommend treating all children with HIV-1 infection. This has changed from the broader use of ART to treat children to improve morbidity and minimise mortality. However, prior to current recommendations, not everyone with HIV-1 received timely treatment. What happens to the paediatric immune system when HIV-1 replication is not appropriately supressed remains unclear. 11 samples from adolescents with HIV-1 on ART and uninfected controls in the UK, aged 12-25 years, were examined; overall, adolescents with CD4⁺ counts > 500/μl and a viral load < 50 copies/ml were compared with adolescents with CD4⁺ counts < 500/μl and a viral load > 50 copies/ml at time of sampling. Measurements of thymic output were combined with high throughput next generation sequencing and bioinformatics to systematically organize CD4⁺ and CD8⁺ T cell receptor (TCR) repertoires. TCR repertoire diversity, clonal expansions, TCR sequence sharing, and formation of TCR clusters in HIV-1 infected adolescents with successful HIV-1 suppression were compared to adolescents with ineffective HIV-1 suppression. Thymic output and CD4⁺ T cell numbers were decreased in HIV-1 infected adolescents with poor HIV-1 suppression. A strong homeostatic TCR response, driven by the decreased CD4⁺ T cell compartment and reduced thymic output was observed in the virally uncontrolled HIV-1-infected adolescents. Formation of abundant robust TCR clusters and structurally related TCRs were found in the adolescents with effective HIV-1 suppression. Numerous CD4⁺ T cell numbers in the virally controlled adolescents emphasize the importance of high thymic output and formation of robust TCR clusters in the maintenance of HIV-1 suppression. While the profound capacity for immune recovery in children may allow better opportunity to deal with immunological stress, when ART is taken appropriately, this study demonstrates new insights into the unique paediatric immune system and the immunological changes when HIV-1 replication is ongoing.

Cyclophosphamide enhances the antitumor potency of GITR engagement by increasing oligoclonal cytotoxic T cell fitness

Only a subset of cancer patients responds to checkpoint blockade inhibition in the clinic. Strategies to overcome resistance are promising areas of investigation. Targeting glucocorticoid-induced tumor necrosis factor receptor-related protein (GITR) has shown efficacy in preclinical models, but GITR engagement is ineffective in controlling advanced, poorly immunogenic tumors, such as B16 melanoma, and has not yielded benefit in clinical trials. The alkylating agent cyclophosphamide (CTX) depletes regulatory T cells (Tregs), expands tumor-specific effector T cells (Teffs) via homeostatic proliferation, and induces immunogenic cell death. GITR agonism has an inhibitory effect on Tregs and activates Teffs. We therefore hypothesized that CTX and GITR agonism would promote effective antitumor immunity. Here we show that the combination of CTX and GITR agonism controlled tumor growth in clinically relevant mouse models. Mechanistically, we show that the combination therapy caused tumor cell death, clonal expansion of highly active CD8+ T cells, and depletion of Tregs by activation-induced cell death. Control of tumor growth was associated with the presence of an expanded population of highly activated, tumor-infiltrating, oligoclonal CD8+ T cells that led to a diminished TCR repertoire. Our studies show that the combination of CTX and GITR agonism is a rational chemoimmunotherapeutic approach that warrants further clinical investigation.

Features of repertoire diversity and gene expression in human cytotoxic T cells following allogeneic hematopoietic cell transplantation

Cytomegalovirus reactivation is still a critical concern following allogeneic hematopoietic cell transplantation, and cellular immune reconstitution of cytomegalovirus-specific cytotoxic T-cells is necessary for the long-term control of cytomegalovirus reactivation after allogeneic hematopoietic cell transplantation. Here we show the features of repertoire diversity and the gene expression profile of HLA-A24 cytomegalovirus-specific cytotoxic T-cells in actual recipients according to the cytomegalovirus reactivation pattern. A skewed preference for BV7 genes and sequential “G” amino acids motif is observed in complementarity-determining region-3 of T cell receptor-β. Increased binding scores are observed in T-cell clones with complementarity-determining region-3 of T cell receptor-β with a “(G)GG” motif. Single-cell RNA-sequence analyses demonstrate the homogenous distribution of the gene expression profile in individual cytomegalovirus-specific cytotoxic T-cells within each recipient. On the other hand, bulk RNA-sequence analyses reveal that gene expression profiles among patients are different according to the cytomegalovirus reactivation pattern, and are associated with cytokine production or cell division. These methods and results can help us to better understand immune reconstitution following hematopoietic cell transplantation, leading to future studies on the clinical application of adoptive T-cell therapies.

Cytomegalovirus reactivation is an important concern after allogeneic stem cell transplantation (allo-HCT) or organ transplantation. Here, Hideki Nakasone et al. investigate changes in repertoire diversity and gene expression among clinically-transferred T cells to improve our understanding of immune reconstitution following allo-HCT.

Tools for fundamental analysis functions of TCR repertoires: a systematic comparison

The full set of T cell receptors (TCRs) in an individual is known as his or her TCR repertoire. Defining TCR repertoires under physiological conditions and in response to a disease or vaccine may lead to a better understanding of adaptive immunity and thus has great biological and clinical value. In the past decade, several high-throughput sequencing-based tools have been developed to assign TCRs to germline genes and to extract complementarity-determining region 3 (CDR3) sequences using different algorithms. Although these tools claim to be able to perform the full range of fundamental TCR repertoire analyses, there is no clear consensus of which tool is best suited to particular projects. Here, we present a systematic analysis of 12 available TCR repertoire analysis tools using simulated data, with an emphasis on fundamental analysis functions. Our results shed light on the detailed functions of TCR repertoire analysis tools and may therefore help researchers in the field to choose the right tools for their particular experimental design.

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.

Plasticity of the Immune System in Children Following Treatment Interruption in HIV-1 Infection

It is intriguing that, unlike adults with HIV-1, children with HIV-1 reach a greater CD4+ T cell recovery following planned treatment cessation. The reasons for the better outcomes in children remain unknown but may be related to increased thymic output and diversity of T cell receptor repertoires. HIV-1 infected children from the PENTA 11 trial tolerated planned treatment interruption without adverse long-term clinical, virological, or immunological consequences, once antiretroviral therapy was re-introduced. This contrasts to treatment interruption trials of HIV-1 infected adults, who had rapid changes in T cells and slow recovery when antiretroviral therapy was restarted. How children can develop such effective immune responses to planned treatment interruption may be critical for future studies. PENTA 11 was a randomized, phase II trial of planned treatment interruptions in HIV-1-infected children (ISRCTN 36694210). In this sub-study, eight patients in long-term follow-up were chosen with CD4+ count>500/ml, viral load <50c/ml at baseline: four patients on treatment interruption and four on continuous treatment. Together with measurements of thymic output, we used high-throughput next generation sequencing and bioinformatics to systematically organize memory CD8+ and naïve CD4+ T cell receptors according to diversity, clonal expansions, sequence sharing, antigen specificity, and T cell receptor similarities following treatment interruption compared to continuous treatment. We observed an increase in thymic output following treatment interruption compared to continuous treatment. This was accompanied by an increase in T cell receptor clonal expansions, increased T cell receptor sharing, and higher sequence similarities between patients, suggesting a more focused T cell receptor repertoire. The low numbers of patients included is a limitation and the data should be interpreted with caution. Nonetheless, the high levels of thymic output and the high diversity of the T cell receptor repertoire in children may be sufficient to reconstitute the T cell immune repertoire and reverse the impact of interruption of antiretroviral therapy. Importantly, the effective T cell receptor repertoires following treatment interruption may inform novel therapeutic strategies in children infected with HIV-1.

T cell receptor revision and immune repertoire changes in autoimmune diseases

Autoimmune disease (AID) is a condition in which the immune system breaks down and starts to attack the body. Some common AIDs include systemic lupus erythematosus, rheumatoid arthritis, type 1 diabetes mellitus and so forth. The changes in T-cell receptor (TCR) repertoire have been found in several autoimmune diseases, and may be responsible for the breakdown of peripheral immune tolerance. In this review, we discussed the processes of TCR revision in peripheral immune environment, the changes in TCR repertoire that occurred in various AIDs, and the specifically expanded T cell clones. We hope our discussion can provide insights for the future studies, helping with the discovery of disease biomarkers and expanding the strategies of immune-targeted therapy. HighlightsRestricted TCR repertoire and biased TCR-usage are found in a variety of AIDs.TCR repertoire shows tissue specificity in a variety of AID diseases.The relationship between TCR repertoire diversity and disease activity is still controversial in AIDs.Dominant TCR clonotypes may help to discover new disease biomarkers and expand the strategies of immune-targeted therapy.

Autoencoder based local T cell repertoire density can be used to classify samples and T cell receptors

Recent advances in T cell repertoire (TCR) sequencing allow for the characterization of repertoire properties, as well as the frequency and sharing of specific TCR. However, there is no efficient measure for the local density of a given TCR. TCRs are often described either through their Complementary Determining region 3 (CDR3) sequences, or theirV/J usage, or their clone size. We here show that the local repertoire density can be estimated using a combined representation of these components through distance conserving autoencoders and Kernel Density Estimates (KDE). We present ELATE-an Encoder-based LocAl Tcr dEnsity and show that the resulting density of a sample can be used as a novel measure to study repertoire properties. The cross-density between two samples can be used as a similarity matrix to fully characterize samples from the same host. Finally, the same projection in combination with machine learning algorithms can be used to predict TCR-peptide binding through the local density of known TCRs binding a specific target.

Generative models of T-cell receptor sequences

T-cell receptors (TCR) are key proteins of the adaptive immune system, generated randomly in each individual, whose diversity underlies our ability to recognize infections and malignancies. Modeling the distribution of TCR sequences is of key importance for immunology and medical applications. Here, we compare two inference methods trained on high-throughput sequencing data: a knowledge-guided approach, which accounts for the details of sequence generation, supplemented by a physics-inspired model of selection; and a knowledge-free variational autoencoder based on deep artificial neural networks. We show that the knowledge-guided model outperforms the deep network approach at predicting TCR probabilities, while being more interpretable, at a lower computational cost.

Selection influences naive CD8+ TCR-β repertoire sharing

Within each individual, the adaptive immune system generates a repertoire of cells expressing receptors capable of recognizing diverse potential pathogens. The theoretical diversity of the T-cell receptor (TCR) repertoire exceeds the actual size of the T-cell population in an individual by several orders of magnitude - making the observation of identical TCRs in different individuals extremely improbable if all receptors were equally likely. Despite this disparity between the theoretical and the realized diversity of the repertoire, these 'public' receptor sequences have been identified in autoimmune, cancer and pathogen interaction contexts. Biased generation processes explain the presence of public TCRs in the naive repertoire, but do not adequately explain the different abundances of these public TCRs. We investigate and characterize the distribution of genomic TCR-β sequences of naive CD8+ T cells from three genetically identical mice, comparing non-productive (non-functional sequences) and productive sequences. We find public TCR-β sequences at higher abundances compared with unshared sequences in the productive, but not in the non-productive, repertoire. We show that neutral processes such as recombination biases, codon degeneracy and generation probability do not fully account for these differences, and conclude that thymic or peripheral selection plays an important role in increasing the abundances of public TCR-β sequences.

CDR3 and V genes show distinct reconstitution patterns in T cell repertoire post-allogeneic bone marrow transplantation

Restoration of T cell repertoire diversity after allogeneic bone marrow transplantation (allo-BMT) is crucial for immune recovery. T cell diversity is produced by rearrangements of germline gene segments (V (D) and J) of the T cell receptor (TCR) α and β chains, and selection induced by binding of TCRs to MHC-peptide complexes. Multiple measures were proposed for this diversity. We here focus on the V-gene usage and the CDR3 sequences of the beta chain. We compared multiple T cell repertoires to follow T cell repertoire changes post-allo-BMT in HLA-matched related donor and recipient pairs. Our analyses of the differences between donor and recipient complementarity determining region 3 (CDR3) beta composition and V-gene profile show that the CDR3 sequence composition does not change during restoration, implying its dependence on the HLA typing. In contrast, V-gene usage followed a time-dependent pattern, initially following the donor profile and then shifting back to the recipients' profile. The final long-term repertoire was more similar to that of the recipient's original one than the donor's; some recipients converged within months, while others took multiple years. Based on the results of our analyses, we propose that donor-recipient V-gene distribution differences may serve as clinical biomarkers for monitoring immune recovery.

CD4+ T Follicular Helper Cells in Human Tonsils and Blood Are Clonally Convergent but Divergent from Non-Tfh CD4+ Cells

T follicular helper (Tfh) cells are fundamental for B cell selection and antibody maturation in germinal centers. Circulating Tfh (cTfh) cells constitute a minor proportion of the CD4⁺ T cells in peripheral blood, but their clonotypic relationship to Tfh populations resident in lymph nodes and the extent to which they differ from non-Tfh CD4⁺ cells have been unclear. Using donor-matched blood and tonsil samples, we investigate T cell receptor (TCR) sharing between tonsillar Tfh cells and peripheral Tfh and non-Tfh cell populations. TCR transcript sequencing reveals considerable clonal overlap between peripheral and tonsillar Tfh cell subsets as well as a clear distinction between Tfh and non-Tfh cells. Furthermore, influenza-specific cTfh cell clones derived from blood can be found in the repertoire of tonsillar Tfh cells. Therefore, human blood samples can be used to gain insight into the specificity of Tfh responses occurring in lymphoid tissues, provided that cTfh subsets are studied.

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.

Predicting Cross-Reactivity and Antigen Specificity of T Cell Receptors

Adaptive immune recognition is mediated by specific interactions between heterodimeric T cell receptors (TCRs) and their cognate peptide-MHC (pMHC) ligands, and the methods to accurately predict TCR:pMHC interaction would have profound clinical, therapeutic and pharmaceutical applications. Herein, we review recent developments in predicting cross-reactivity and antigen specificity of TCR recognition. We discuss current experimental and computational approaches to investigate cross-reactivity and antigen-specificity of TCRs and highlight how integrating kinetic, biophysical and structural features may offer valuable insights in modeling immunogenicity. We further underscore the close inter-relationship of these two interconnected notions and the need to investigate each in the light of the other for a better understanding of T cell responsiveness for the effective clinical applications.

Characteristics of TCR Repertoire Associated With Successful Immune Checkpoint Therapy Responses

Immunotherapies have revolutionized cancer treatment. In particular, immune checkpoint therapy (ICT) leads to durable responses in some patients with some cancers. However, the majority of treated patients do not respond. Understanding immune mechanisms that underlie responsiveness to ICT will help identify predictive biomarkers of response and develop treatments to convert non-responding patients to responding ones. ICT primarily acts at the level of adaptive immunity. The specificity of adaptive immune cells, such as T and B cells, is determined by antigen-specific receptors. T cell repertoires can be comprehensively profiled by high-throughput sequencing at the bulk and single-cell level. T cell receptor (TCR) sequencing allows for sensitive tracking of dynamic changes in antigen-specific T cells at the clonal level, giving unprecedented insight into the mechanisms by which ICT alters T cell responses. Here, we review how the repertoire influences response to ICT and conversely how ICT affects repertoire diversity. We will also explore how changes to the repertoire in different anatomical locations can better correlate and perhaps predict treatment outcome. We discuss the advantages and limitations of current metrics used to characterize and represent TCR repertoire diversity. Discovery of predictive biomarkers could lie in novel analysis approaches, such as network analysis of amino acids similarities between TCR sequences. Single-cell sequencing is a breakthrough technology that can link phenotype with specificity, identifying T cell clones that are crucial for successful ICT. The field of immuno-sequencing is rapidly developing and cross-disciplinary efforts are required to maximize the analysis, application, and validation of sequencing data. Unravelling the dynamic behavior of the TCR repertoire during ICT will be highly valuable for tracking and understanding anti-tumor immunity, biomarker discovery, and ultimately for the development of novel strategies to improve patient outcomes.

The Whole Body as the System in Systems Immunology

The human immune system is comprised of a diverse and interactive network of specialized cells localized in diverse tissues throughout the body, where they mediate protection against pathogens and environmental insults while maintaining tissue homeostasis. Although much of our understanding of human immunology has derived from studies of peripheral blood, recent work utilizing human tissue resources and innovative computational methods have employed a whole-body, systems-based approach, revealing tremendous complexity and heterogeneity of the immune system within individuals and across the population. In this review, we discuss how tissue localization, developmental and age-associated changes, and conditions of health and disease shape the immune response, as well as how improved understanding of interindividual and tissue-specific immunity can be leveraged for developing targeted therapeutics.

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.

Inferring the immune response from repertoire sequencing

High-throughput sequencing of B- and T-cell receptors makes it possible to track immune repertoires across time, in different tissues, and in acute and chronic diseases or in healthy individuals. However, quantitative comparison between repertoires is confounded by variability in the read count of each receptor clonotype due to sampling, library preparation, and expression noise. Here, we present a general Bayesian approach to disentangle repertoire variations from these stochastic effects. Using replicate experiments, we first show how to learn the natural variability of read counts by inferring the distributions of clone sizes as well as an explicit noise model relating true frequencies of clones to their read count. We then use that null model as a baseline to infer a model of clonal expansion from two repertoire time points taken before and after an immune challenge. Applying our approach to yellow fever vaccination as a model of acute infection in humans, we identify candidate clones participating in the response.

Massively parallel interrogation and mining of natively paired human TCRαβ repertoires

T cells engineered to express antigen-specific T cell receptors (TCRs) are potent therapies for viral infections and cancer. However, efficient identification of clinical candidate TCRs is complicated by the size and complexity of T cell repertoires and the challenges of working with primary T cells. Here, we present a high-throughput method to identify TCRs with high functional avidity from diverse human T cell repertoires. The approach uses massively parallel microfluidics to generate libraries of natively paired, full-length TCRαβ clones, from millions of primary T cells, which are then expressed in Jurkat cells. The TCRαβ-Jurkat libraries enable repeated screening and panning for antigen-reactive TCRs using peptide:MHC binding and cellular activation. We captured >2.9 million natively paired TCRαβ clonotypes from six healthy human donors and identified rare (<0.001% frequency) viral antigen–reactive TCRs. We also mined a tumor-infiltrating lymphocyte (TIL) sample from a melanoma patient and identified several tumor-specific TCRs, which, after expression in primary T cells, led to tumor cell killing.

High-throughput sequencing of CD4+ T cell repertoire reveals disease-specific signatures in IgG4-related disease

Background

CD4⁺ T cells play critical roles in the pathogenesis of IgG4-related disease (IgG4-RD). The aim of this study was to investigate the TCR repertoire of peripheral blood CD4⁺ T cells in IgG4-RD.

Methods

The peripheral blood was collected from six healthy controls and eight IgG4-RD patients. TCR β-chain libraries of CD4⁺ T cells were constructed by 5′-rapid amplification of cDNA ends (5′-RACE) and sequenced by Illumina Miseq platform. The relative similarity of TCR repertoires between samples was evaluated according to the total frequencies of shared clonotypes (metric F), correlation of frequencies of shared clonotypes (metric R), and total number of shared clonotypes (metric D).

Results

The clonal expansion and diversity of CD4⁺ T cell repertoire were comparable between healthy controls and IgG4-RD patients, while the proportion of expanded and coding degenerated clones, as an indicator of antigen-driven clonal expansion, was significantly higher in IgG4-RD patients. There was no significant difference in TRBV and TRBJ gene usage between healthy controls and IgG4-RD patients. The complementarity determining region 3 (CDR3) length distribution was skewed towards longer fragments in IgG4-RD. Visualization of relative similarity of TCR repertoires by multi-dimensional scaling analysis showed that TCR repertoires of IgG4-RD patients were separated from that of healthy controls in F and D metrics. We identified 11 IgG4-RD-specific CDR3 amino acid sequences that were expanded in at least 2 IgG4-RD patients, while not detected in healthy controls. According to TCR clonotype networks constructed by connecting all the CDR3 sequences with a Levenshtein distance of 1, 3 IgG4-RD-specific clusters were identified. We annotated the TCR sequences with known antigen specificity according to McPAS-TCR database and found that the frequencies of TCR sequences associated with each disease or immune function were comparable between healthy controls and IgG4-RD patients.

Conclusion

According to our study of CD4⁺ T cells from eight IgG4-RD patients, TCR repertoires of IgG4-RD patients were different from that of healthy controls in the proportion of expanded and coding degenerated clones and CDR3 length distribution. In addition, IgG4-RD-specific TCR sequences and clusters were identified in our study.

VDJdb: a curated database of T-cell receptor sequences with known antigen specificity

The ability to decode antigen specificities encapsulated in the sequences of rearranged T-cell receptor (TCR) genes is critical for our understanding of the adaptive immune system and promises significant advances in the field of translational medicine. Recent developments in high-throughput sequencing methods (immune repertoire sequencing technology, or RepSeq) and single-cell RNA sequencing technology have allowed us to obtain huge numbers of TCR sequences from donor samples and link them to T-cell phenotypes. However, our ability to annotate these TCR sequences still lags behind, owing to the enormous diversity of the TCR repertoire and the scarcity of available data on T-cell specificities. In this paper, we present VDJdb, a database that stores and aggregates the results of published T-cell specificity assays and provides a universal platform that couples antigen specificities with TCR sequences. We demonstrate that VDJdb is a versatile instrument for the annotation of TCR repertoire data, enabling a concatenated view of antigen-specific TCR sequence motifs. VDJdb can be accessed at https://vdjdb.cdr3.net and https://github.com/antigenomics/vdjdb-db.

Single T Cell Sequencing Demonstrates the Functional Role of αβ TCR Pairing in Cell Lineage and Antigen Specificity

Although structural studies of individual T cell receptors (TCRs) have revealed important roles for both the α and β chain in directing MHC and antigen recognition, repertoire-level immunogenomic analyses have historically examined the β chain alone. To determine the amount of useful information about TCR repertoire function encoded within αβ pairings, we analyzed paired TCR sequences from nearly 100,000 unique CD4+ and CD8+ T cells captured using two different high-throughput, single-cell sequencing approaches. Our results demonstrate little overlap in the healthy CD4+ and CD8+ repertoires, with shared TCR sequences possessing significantly shorter CDR3 sequences corresponding to higher generation probabilities. We further utilized tools from information theory and machine learning to show that while α and β chains are only weakly associated with lineage, αβ pairings appear to synergistically drive TCR-MHC interactions. Vαβ gene pairings were found to be the TCR feature most informative of T cell lineage, supporting the existence of germline-encoded paired αβ TCR-MHC interaction motifs. Finally, annotating our TCR pairs using a database of sequences with known antigen specificities, we demonstrate that approximately a third of the T cells possess α and β chains that each recognize different known antigens, suggesting that αβ pairing is critical for the accurate inference of repertoire functionality. Together, these findings provide biological insight into the functional implications of αβ pairing and highlight the utility of single-cell sequencing in immunogenomics.

RNase H-dependent PCR-enabled T-cell receptor sequencing for highly specific and efficient targeted sequencing of T-cell receptor mRNA for single-cell and repertoire analysis

RNase H-dependent PCR-enabled T-cell receptor sequencing (rhTCRseq) can be used to determine paired alpha/beta T-cell receptor (TCR) clonotypes in single cells or perform alpha and beta TCR repertoire analysis in bulk RNA samples. With the enhanced specificity of RNase H-dependent PCR (rhPCR), it achieves TCR-specific amplification and addition of dual-index barcodes in a single PCR step. For single cells, the protocol includes sorting of single cells into plates, generation of cDNA libraries, a TCR-specific amplification step, a second PCR on pooled sample to generate a sequencing library, and sequencing. In the bulk method, sorting and cDNA library steps are replaced with a reverse-transcriptase (RT) reaction that adds a unique molecular identifier (UMI) to each cDNA molecule to improve the accuracy of repertoire-frequency measurements. Compared to other methods for TCR sequencing, rhTCRseq has a streamlined workflow and the ability to analyze single cells in 384-well plates. Compared to TCR reconstruction from single-cell transcriptome sequencing data, it improves the success rate for obtaining paired alpha/beta information and ensures recovery of complete complementarity-determining region 3 (CDR3) sequences, a prerequisite for cloning/expression of discovered TCRs. Although it has lower throughput than droplet-based methods, rhTCRseq is well-suited to analysis of small sorted populations, especially when analysis of 96 or 384 single cells is sufficient to identify predominant T-cell clones. For single cells, sorting typically requires 2-4 h and can be performed days, or even months, before library construction and data processing, which takes ~4 d; the bulk RNA protocol takes ~3 d.

Systematic comparative study of computational methods for T-cell receptor sequencing data analysis

High-throughput sequencing technologies have exposed the possibilities for the in-depth evaluation of T-cell receptor (TCR) repertoires. These studies are highly relevant to gain insights into human adaptive immunity and to decipher the composition and diversity of antigen receptors in physiological and disease conditions. The major objective of TCR sequencing data analysis is the identification of V, D and J gene segments, complementarity-determining region 3 (CDR3) sequence extraction and clonality analysis. With the advancement in sequencing technologies, new TCR analysis approaches and programs have been developed. However, there is still a deficit of systematic comparative studies to assist in the selection of an optimal analysis approach. Here, we present a detailed comparison of 10 state-of-the-art TCR analysis tools on samples with different complexities by taking into account many aspects such as clonotype detection [unique V(D)J combination], CDR3 identification or accuracy in error correction. We used our in silico and experimental data sets with known clonalities enabling the identification of potential tool biases. We also established a new strategy, named clonal plane, which allows quantifying and comparing the clonality of multiple samples. Our results provide new insights into the effect of method selection on analysis results, and it will assist users in the selection of an appropriate analysis method.

High-Throughput Sequence Analysis of Peripheral T-Cell Lymphomas Indicates Subtype-Specific Viral Gene Expression Patterns and Immune Cell Microenvironments

Certain peripheral T-cell lymphomas (PTCLs) have been associated with viral infection, particularly infection with Epstein-Barr virus (EBV). However, a comprehensive virome analysis across PTCLs has not previously been reported. Here we utilized published whole-transcriptome RNA sequencing (RNA-seq) data sets from seven different PTCL studies and new RNA-seq data from our laboratory to screen for virus association, to analyze viral gene expression, and to assess B- and T-cell receptor diversity paradigms across PTCL subtypes. In addition to identifying EBV in angioimmunoblastic T-cell lymphoma (AITL) and extranodal NK/T-cell lymphoma (ENKTL), two PTCL subtypes with well-established EBV associations, we also detected EBV in several cases of anaplastic large-cell lymphoma (ALCL), and we found evidence of infection by the oncogenic viruses Kaposi's sarcoma-associated herpesvirus and human T-cell leukemia virus type 1 in isolated PTCL cases. In AITLs, EBV gene expression analysis showed expression of immediate early, early, and late lytic genes, suggesting either low-level lytic gene expression or productive infection in a subset of EBV-infected B-lymphocyte stromal cells. Deconvolution of immune cell subpopulations demonstrated a greater B-cell signal in AITLs than in other PTCL subtypes, consistent with a larger role for B-cell support in the pathogenesis of AITL. Reconstructed T-cell receptor (TCR) and B-cell receptor (BCR) repertoires demonstrated increased BCR diversity in AITLs, consistent with a possible EBV-driven polyclonal response. These findings indicate potential alternative roles for EBV in PTCLs, in addition to the canonical oncogenic mechanisms associated with EBV latent infection. Our findings also suggest the involvement of other viruses in PTCL pathogenesis and demonstrate immunological alterations associated with these cancers.IMPORTANCE In this study, we utilized next-generation sequencing data from 7 different studies of peripheral T-cell lymphoma (PTCL) patient samples to globally assess viral associations, provide insights into the contributions of EBV gene expression to the tumor phenotype, and assess the unique roles of EBV in modulating the immune cell tumor microenvironment. These studies revealed potential roles for EBV replication genes in some PTCL subtypes, the possible role of additional human tumor viruses in rare cases of PTCLs, and a role for EBV in providing a unique immune microenvironmental niche in one subtype of PTCLs. Together, these studies provide new insights into the understudied role of tumor viruses in PTCLs.

Characterization of T-Cell Receptor Repertoire in Patients with Rheumatoid Arthritis Receiving Biologic Therapies

Rheumatoid arthritis (RA) is a systematic autoimmune disease, predominantly causing chronic polyarticular inflammation and joint injury of patients. For the treatment of RA, biologic disease-modifying antirheumatic drugs (bDMARDs) have been used to reduce inflammation and to interfere with disease progression through targeting and mediating the immune system. Although the therapeutic effects of bDMARDs in RA patients have been widely reported, whether these drugs also play important roles in T-cell repertoire status is still unclear. We therefore designed the study to identify the role of T-cell repertoire profiles in RA patients with different types of bDMARD treatments. A high-throughput sequencing approach was applied to profile the T-cell receptor beta chain (TCRB) repertoire of circulating T lymphocytes in eight patients given adalimumab (anti-TNF-α) with/without the following use of either rituximab (anti-CD20) or tocilizumab (anti-IL6R). We subsequently analyzed discrepancies in the clonal diversity and CDR3 length distribution as well as usages of the V and J genes of TCRB repertoire and interrogated the association between repertoire diversity and disease activities followed by the treatment of bDMARDs in these RA patients. All groups of patients showed well-controlled DAS28 scores (<2.6) after different treatment regimens of drugs and displayed no significant statistical differences in repertoire diversity, distribution of CDR3 lengths, and usage of V and J genes of TCRB. Nonetheless, a trend between overall TCRB repertoire diversity and disease activity scores in all bDMARD-treated RA patients was observed. Additionally, age was found to be associated with repertoire diversity in RA patients treated with bDMARDs. Through the profiling of the TCR repertoire in RA patients receiving different biologic medications, our study indicated an inverse tendency between TCR repertoire diversity and disease activity after biologic treatment in RA patients.

An Economical, Quantitative, and Robust Protocol for High-Throughput T Cell Receptor Sequencing from Tumor or Blood

The T cell receptor repertoire provides a window to the cellular adaptive immune response within a tumor, and has the potential to identify specific and personalized biomarkers for tracking host responses during cancer therapy, including immunotherapy. We describe a protocol for amplifying, sequencing, and analyzing T cell receptors which is economical, robust, sensitive, and versatile. The key experimental step is the ligation of a single-stranded oligonucleotide to the 3' end of the T cell receptor cDNA, which allows easy amplification of all possible rearrangements using only a single set of primers per locus, while simultaneously introducing a unique molecular identifier to label each starting cDNA molecule. After sequencing, this molecular identifier can be used to correct both sequence errors and the effects of differential PCR amplification efficiency, thus producing a more accurate measure of the true T cell receptor frequency within the sample. This method has been applied to the analysis of unfractionated human tumor lysates, subpopulations of tumor-infiltrating lymphocytes, and peripheral blood samples from patients with a variety of solid tumors.

Vedolizumab is associated with changes in innate rather than adaptive immunity in patients with inflammatory bowel disease

OBJECTIVE

Vedolizumab, a monoclonal antibody directed against the integrin heterodimer α4β7, is approved for the treatment of Crohn's disease and ulcerative colitis. The efficacy of vedolizumab has been suggested to result from inhibition of intestinal T cell trafficking although human data to support this conclusion are scarce. We therefore performed a comprehensive analysis of vedolizumab-induced alterations in mucosal and systemic immunity in patients with inflammatory bowel disease (IBD), using anti-inflammatory therapy with the TNFα antibody infliximab as control.

DESIGN

Immunophenotyping, immunohistochemistry, T cell receptor profiling and RNA sequencing were performed using blood and colonic biopsies from patients with IBD before and during treatment with vedolizumab (n=18) or, as control, the anti-TNFα antibody infliximab (n=20). Leucocyte trafficking in vivo was assessed using single photon emission computed tomography and endomicroscopy.

RESULTS

Vedolizumab was not associated with alterations in the abundance or phenotype of lamina propria T cells and did not affect the mucosal T cell repertoire or leucocyte trafficking in vivo. Surprisingly, however, α4β7 antibody treatment was associated with substantial effects on innate immunity including changes in macrophage populations and pronounced alterations in the expression of molecules involved in microbial sensing, chemoattraction and regulation of the innate effector response. These effects were specific to vedolizumab, not observed in response to the TNFα antibody infliximab, and associated with inhibition of intestinal inflammation.

CONCLUSION

Our findings suggest that modulation of innate immunity contributes to the therapeutic efficacy of vedolizumab in IBD.

TRIAL REGISTRATION NUMBER

NCT02694588.

Evaluation of the capacities of mouse TCR profiling from short read RNA-seq data

Profiling T cell receptor (TCR) repertoire via short read transcriptome sequencing (RNA-Seq) has a unique advantage of probing simultaneously TCRs and the genome-wide RNA expression of other genes. However, compared to targeted amplicon approaches, the shorter read length is more prone to mapping error. In addition, only a small percentage of the genome-wide reads may cover the TCR loci and thus the repertoire could be significantly under-sampled. Although this approach has been applied in a few studies, the utility of transcriptome sequencing in probing TCR repertoires has not been evaluated extensively. Here we present a systematic assessment of RNA-Seq in TCR profiling. We evaluate the power of both Fluidigm C1 full-length single cell RNA-Seq and bulk RNA-Seq in characterizing the repertoires of different diversities under either naïve conditions or after immunogenic challenges. Standard read length and sequencing coverage were employed so that the evaluation was conducted in accord with the current RNA-Seq practices. Despite high sequencing depth in bulk RNA-Seq, we encountered difficulty quantifying TCRs with low transcript abundance (<1%). Nevertheless, top enriched TCRs with an abundance of 1–3% or higher can be faithfully detected and quantified. When top TCR sequences are of interest and transcriptome sequencing is available, it is worthwhile to conduct a TCR profiling using the RNA-Seq data.

Advanced model systems and tools for basic and translational human immunology

There are fundamental differences between humans and the animals we typically use to study the immune system. We have learned much from genetically manipulated and inbred animal models, but instances in which these findings have been successfully translated to human immunity have been rare. Embracing the genetic and environmental diversity of humans can tell us about the fundamental biology of immune cell types and the elasticity of the immune system. Although people are much more immunologically diverse than conventionally housed animal models, tools and technologies are now available that permit high-throughput analysis of human samples, including both blood and tissues, which will give us deep insights into human immunity in health and disease. As we gain a more detailed picture of the human immune system, we can build more sophisticated models to better reflect this complexity, both enabling the discovery of new immunological mechanisms and facilitating translation into the clinic.

The expanding role of systems immunology in decoding the T cell receptor repertoire

T cells play a crucial role in the immune system's defense against many infectious diseases, including persistent infections for which no effective vaccines currently exist. The T cell component of the adaptive immune system is highly complex involving a constantly evolving landscape of various inter-related T cell populations. These T cell populations are characterized by their phenotypic and functional properties as well as the collection, or repertoire, of T cell receptors (TCR) that mediate T cell recognition of antigenic peptides derived from pathogens. Understanding the various processes and factors that impact the development and evolution of the broader T cell repertoire available to recognize and respond to pathogens and the characteristics of antigen-experienced T cell repertoires associated with effective immune control of pathogens is critical to the rational design of T cell-based vaccines and therapies. In this article we discuss, using examples of recent research, the promise that systems immunology approaches, involving quantitative analysis and mathematical and computational modeling of immunological data, hold for decoding the complex TCR repertoire system in the current era of advancing technologies.

Exploring the pre-immune landscape of antigen-specific T cells

Background

Adaptive immune responses to newly encountered pathogens depend on the mobilization of antigen-specific clonotypes from a vastly diverse pool of naive T cells. Using recent advances in immune repertoire sequencing technologies, models of the immune receptor rearrangement process, and a database of annotated T cell receptor (TCR) sequences with known specificities, we explored the baseline frequencies of T cells specific for defined human leukocyte antigen (HLA) class I-restricted epitopes in healthy individuals.

Methods

We used a database of TCR sequences with known antigen specificities and a probabilistic TCR rearrangement model to estimate the baseline frequencies of TCRs specific to distinct antigens epitopespecificT-cells. We verified our estimates using a publicly available collection of TCR repertoires from healthy individuals. We also interrogated a database of immunogenic and non-immunogenic peptides is used to link baseline T-cell frequencies with epitope immunogenicity.

Results

Our findings revealed a high degree of variability in the prevalence of T cells specific for different antigens that could be explained by the physicochemical properties of the corresponding HLA class I-bound peptides. The occurrence of certain rearrangements was influenced by ancestry and HLA class I restriction, and umbilical cord blood samples contained higher frequencies of common pathogen-specific TCRs. We also identified a quantitative link between specific T cell frequencies and the immunogenicity of cognate epitopes presented by defined HLA class I molecules.

Conclusions

Our results suggest that the population frequencies of specific T cells are strikingly non-uniform across epitopes that are known to elicit immune responses. This inference leads to a new definition of epitope immunogenicity based on specific TCR frequencies, which can be estimated with a high degree of accuracy in silico, thereby providing a novel framework to integrate computational and experimental genomics with basic and translational research efforts in the field of T cell immunology.

Electronic supplementary material

The online version of this article (10.1186/s13073-018-0577-7) contains supplementary material, which is available to authorized users.

Recasting Human Vδ1 Lymphocytes in an Adaptive Role

γδ T cells are unconventional lymphocytes commonly described as 'innate-like' in function, which can respond in both a T cell receptor (TCR)-independent and also major histocompatibility complex (MHC)-unrestricted TCR-dependent manner. While the relative importance of TCR recognition had remained unclear, recent studies revealed that human Vδ1 T cells display unexpected parallels with adaptive αβ T cells. Vδ1 T cells undergo profound and highly focussed clonal expansion from an initially diverse and private TCR repertoire, most likely in response to specific immune challenges. Concomitantly, they differentiate from a Vδ1 T cell naïve (Tnaïve) to a Vδ1 T cell effector (Teffector) phenotype, marked by the downregulation of lymphoid homing receptors and upregulation of peripheral homing receptors and effector markers. This suggests that an adaptive paradigm applies to Vδ1 T cells, likely involving TCR-dependent but MHC-unrestricted responses to microbial and non-microbial challenges.

RepSeq Data Representativeness and Robustness Assessment by Shannon Entropy

High-throughput sequencing (HTS) has the potential to decipher the diversity of T cell repertoires and their dynamics during immune responses. Applied to T cell subsets such as T effector and T regulatory cells, it should help identify novel biomarkers of diseases. However, given the extreme diversity of TCR repertoires, understanding how the sequencing conditions, including cell numbers, biological and technical sampling and sequencing depth, impact the experimental outcome is critical to proper use of these data. Here, we assessed the representativeness and robustness of TCR repertoire diversity assessment according to experimental conditions. By comparative analyses of experimental datasets and computer simulations, we found that (i) for small samples, the number of clonotypes recovered is often higher than the number of cells per sample, even after removing the singletons; (ii) high-sequencing depth for small samples alters the clonotype distributions, which can be corrected by filtering the datasets using Shannon entropy as a threshold; and (iii) a single sequencing run at high depth does not ensure a good coverage of the clonotype richness in highly polyclonal populations, which can be better covered using multiple sequencing. Altogether, our results warrant better understanding and awareness of the limitation of TCR diversity analyses by HTS and justify the development of novel computational tools for improved modeling of the highly complex nature of TCR repertoires.

Computational Strategies for Dissecting the High-Dimensional Complexity of Adaptive Immune Repertoires

The adaptive immune system recognizes antigens via an immense array of antigen-binding antibodies and T-cell receptors, the immune repertoire. The interrogation of immune repertoires is of high relevance for understanding the adaptive immune response in disease and infection (e.g., autoimmunity, cancer, HIV). Adaptive immune receptor repertoire sequencing (AIRR-seq) has driven the quantitative and molecular-level profiling of immune repertoires, thereby revealing the high-dimensional complexity of the immune receptor sequence landscape. Several methods for the computational and statistical analysis of large-scale AIRR-seq data have been developed to resolve immune repertoire complexity and to understand the dynamics of adaptive immunity. Here, we review the current research on (i) diversity, (ii) clustering and network, (iii) phylogenetic, and (iv) machine learning methods applied to dissect, quantify, and compare the architecture, evolution, and specificity of immune repertoires. We summarize outstanding questions in computational immunology and propose future directions for systems immunology toward coupling AIRR-seq with the computational discovery of immunotherapeutics, vaccines, and immunodiagnostics.

Comparative analysis of murine T-cell receptor repertoires

For understanding the rules and laws of adaptive immunity, high-throughput profiling of T-cell receptor (TCR) repertoires becomes a powerful tool. The structure of TCR repertoires is instructive even before the antigen specificity of each particular receptor becomes available. It embodies information about the thymic and peripheral selection of T cells; the readiness of an adaptive immunity to withstand new challenges; the character, magnitude and memory of immune responses; and the aetiological and functional proximity of T-cell subsets. Here, we describe our current analytical approaches for the comparative analysis of murine TCR repertoires, and show several examples of how these approaches can be applied for particular experimental settings. We analyse the efficiency of different metrics used for estimation of repertoire diversity, repertoire overlap, V-gene and J-gene segments usage similarity, and amino acid composition of CDR3. We discuss basic differences of these metrics and their advantages and limitations in different experimental models, and we provide guidelines for choosing an efficient way to lead a comparative analysis of TCR repertoires. Applied to the various known and newly developed mouse models, such analysis should allow us to disentangle multiple sophisticated puzzles in adaptive immunity.

Deep sequencing of the T cell receptor β repertoire reveals signature patterns and clonal drift in atherosclerotic plaques and patients

The T cell receptor (TCR) β repertoire directly reflects the status of T cell function. Meanwhile, the immune/inflammatory responses regulated by T cells are the critical determinants of atherosclerosis development. However, due to technical limitations, the composition and molecular characteristics of the TCR repertoire in atherosclerotic patients have not been fully elucidated. In the present study, we use powerful immune repertoire sequencing technology to study this issue. Results show that the utilization of V and/or J genes and the diversity of TCRβ repertoire in atherosclerotic plaques are significantly reduced compared to those in the peripheral blood of normal subjects and atherosclerotic patients. The frequencies of the common T cell clones with certain lengths of the complement determining region 3 regions are notably different among all groups. The high-frequency common clones are also increased in the atherosclerotic plaques compared to that in the other two groups. The expansion of several T cell clonotypes (V29-1J2-1, V20-1J1-6, V6-3J2-7 and V11-2J2-2) is validated in atherosclerotic patients. In short, this study reveals that the diversity of TCR β repertoire significantly decreases in atherosclerotic plaques, probably because of the reduced utilization of VJ genes and marked expansion of some T cell subclones. It provides the basis for understanding the roles of T lymphocytes in the pathogenesis of atherosclerosis.

Quantitative Characterization of the T Cell Receptor Repertoire of Naïve and Memory Subsets Using an Integrated Experimental and Computational Pipeline Which Is Robust, Economical, and Versatile

The T cell receptor (TCR) repertoire can provide a personalized biomarker for infectious and non-infectious diseases. We describe a protocol for amplifying, sequencing, and analyzing TCRs which is robust, sensitive, and versatile. The key experimental step is ligation of a single-stranded oligonucleotide to the 3' end of the TCR cDNA. This allows amplification of all possible rearrangements using a single set of primers per locus. It also introduces a unique molecular identifier to label each starting cDNA molecule. This molecular identifier is used to correct for sequence errors and for effects of differential PCR amplification efficiency, thus producing more accurate measures of the true TCR frequency within the sample. This integrated experimental and computational pipeline is applied to the analysis of human memory and naive subpopulations, and results in consistent measures of diversity and inequality. After error correction, the distribution of TCR sequence abundance in all subpopulations followed a power law over a wide range of values. The power law exponent differed between naïve and memory populations, but was consistent between individuals. The integrated experimental and analysis pipeline we describe is appropriate to studies of T cell responses in a broad range of physiological and pathological contexts.

Systems immunology: just getting started

Systems-biology approaches in immunology take various forms, but here we review strategies for measuring a broad swath of immunological functions as a means of discovering previously unknown relationships and phenomena and as a powerful way of understanding the immune system as a whole. This approach has rejuvenated the field of vaccine development and has fostered hope that new ways will be found to combat infectious diseases that have proven refractory to classical approaches. Systems immunology also presents an important new strategy for understanding human immunity directly, taking advantage of the many ways the immune system of humans can be manipulated.

Overview of methodologies for T-cell receptor repertoire analysis

Background

The T-cell receptor (TCR), located on the surface of T cells, is responsible for the recognition of the antigen-major histocompatibility complex, leading to the initiation of an inflammatory response. Analysing the TCR repertoire may help to gain a better understanding of the immune system features and of the aetiology and progression of diseases, in particular those with unknown antigenic triggers. The extreme diversity of the TCR repertoire represents a major analytical challenge; this has led to the development of specialized methods which aim to characterize the TCR repertoire in-depth. Currently, next generation sequencing based technologies are most widely employed for the high-throughput analysis of the immune cell repertoire.

Results

Here, we report on the latest methodological advancements in the field by describing and comparing the available tools; from the choice of the starting material and library preparation method, to the sequencing technologies and data analysis.

Finally, we provide a practical example and our own experience by reporting some exemplary results from a small internal benchmark study, where current approaches from the literature and the market are employed and compared.

Conclusions

Several valid methods for clonotype identification and TCR repertoire analysis exist, however, a gold standard method for the field has not yet been identified. Depending on the purpose of the scientific study, some approaches may be more suitable than others. Finally, due to possible method specific biases, scientists must be careful when comparing results obtained using different methods.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-017-0379-9) contains supplementary material, which is available to authorized users.

Immune dysregulation in immunodeficiency disorders: The role of T-cell receptor sequencing

Immune dysregulation is a prominent feature of primary immunodeficiency disorders, which commonly manifested as autoimmunity, cytopenias and inflammatory bowel disease. In partial T-cell immunodeficiency disorders, it has been proposed that the imbalance between effector and regulatory T-cells drives the breakdown of peripheral tolerance. While there is no robust test for immune dysregulation, the T-cell receptor repertoire is used as a surrogate marker, and has been shown to be perturbed in a number of immunodeficiency disorders featuring immune dysregulation including Omenn's Syndrome, Wiskott-Aldrich Syndrome, and common variable immunodeficiency. This review discusses how recent advances in TCR next-generation sequencing and bioinformatics have led to the in-depth characterization of CDR3 sequences and an exponential growth in examinable parameters. Specifically, we highlight the use of junctional diversity as a means to differentiate intrinsic T-cell defects from secondary causes of repertoire perturbation in primary immunodeficiency disorders. However, key questions, such as the identity of antigenic targets for large, expanded T-cell clonotypes, remain unanswered despite the fact that such clones are likely to play a pathogenic role in driving immune dysregulation and autoimmunity. Finally, we discuss a number of emerging technologies such as in silico reconstruction, high-throughput pairwise αβ sequencing and single-cell RNAseq that offer the potential to define the antigenic epitope and function of a given T-cell, thereby enhancing our understanding in this field.