Skip to main content
  • Methodology article
  • Open access
  • Published:

Overview of methodologies for T-cell receptor repertoire analysis

Abstract

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.

Background

T cell mediated antigen recognition depends on the interaction of the T-cell receptor (TCR) with the antigen-major histocompatibility complex (MHC) molecules (Fig. 1a). TCRs are highly diverse heterodimers, consisting of a combination of α and β chains (αβ TCR) expressed by the majority of T cells, or γδ chains (γδ TCR) expressed by T cells in peripheral blood (1–5%) and T cells found at mucosal sites [1]. Similar to immunoglobulins expressed by B cells – membrane bound immunoglobulins are often referred to as B-cell receptors (BCRs) – the TCR chains consist of a variable region, important for antigen recognition, and a constant region. The variable region of TCRα and δ chains is encoded by a number of variable (V) and joining (J) genes, while TCRβ and γ chains are additionally encoded by diversity (D) genes [2, 3]. During VDJ recombination, one random allele of each gene segment is recombined with the others to form a functional variable region (Fig. 1b). Recombination of the variable region with a constant gene segment results in a functional TCR chain transcript. Additionally, random nucleotides are added and/or deleted at the junction sites between the gene segments. This process leads to strong combinatorial (depending on which gene regions will recombine) and junctional diversity (which and how many nucleotides will be added/deleted), resulting in a large and highly variable TCR repertoire, which will ensure the identification of a plethora of antigens. Additional diversity is achieved by the pairing of α and β or γ and δ chains to form a functional TCR [4].

Fig. 1
figure 1

Interaction between an antigen presenting cell (APC) and a T cell, and V(D)J recombination. a Interaction between the antigen–major histocompatibility complex (MHC) and the αβ T-cell receptor (TCR). b V(D)J recombination: During T cell development, the loci that encode T-cell receptor α and β-chains are rearranged. For both loci, variable (V) and joining (J) gene segments, and an additional diversity (D) gene segment for the β-chain, are recombined to form the final rearranged TCR DNA sequence. This process also involves the deletion and insertion of nucleotides at the V-D, D-J and V-J junctions (not shown). Following transcription, the sequence between the recombined V (D)J regions and the gene encoding the constant (C) region is removed by splicing. The complementarity-determining region (CDR) 3 is encoded by the V (D) J junction, whereas the CDR1 and CDR2 loops are encoded within the germline V gene

Each TCR chain contains three hypervariable loops in its structure, termed complementarity determining regions (CDR1–3). CDR1 and 2 are encoded by V genes and are required for interaction of the TCR with the MHC complex. CDR3, however, is encoded by the junctional region between the V and J or D and J genes and is therefore highly variable. It plays an essential role in the interaction of the TCR with the peptide-MHC complex, as it is the region of the TCR in direct contact with the peptide antigen. For this reason, CDR3 is often used as the region of interest to determine T cell clonotypes, as it is highly unlikely that two T cells will express the same CDR3 nucleotide sequence, unless they have derived from the same clonally expanded T cell [2, 4].

The sum of all TCRs by the T cells of one individual is termed the TCR repertoire or TCR profile. The TCR repertoire can change greatly with the onset and progression of diseases, which is why scientists are becoming more and more interested in determining the immune repertoire status under different disease conditions, such as cancer, autoimmune, inflammatory and infectious diseases. For example, Muraro et al. used the TCR repertoire to analyse the effect of autologous stem cell transplantation on T cell populations in multiple sclerosis patients [5]. In cancer, cytotoxic T cells can kill tumour cells upon recognition of tumour specific antigens. Some studies have tried to identify specific T cell clonotypes involved in this process by analysing tumour infiltrating lymphocyte repertoires [6,7,8].

The main challenge while studying the immune repertoire is its diversity. VDJ recombination of the different TCR genes could theoretically generate between 1015 and 1020 TCR chains. Despite this, the actual diversity present in a human body is estimated at around 1013 different clonotypes [9], implying that the afore-described seemingly random TCR development is obviously not random at all and is subject to different constraints. Moreover, while there are TCRs that are common in the general population, recent high-resolution studies have shown that the majority of TCRs is rare (in analogy to common vs. rare genomic variants) [10, 11]. This is one of the reasons why precise methods are necessary to properly investigate complete individual immune-repertoires. In the past decades, different techniques were developed which enabled the study of the immune repertoire. Monoclonal antibodies allowed the analysis of specific V gene subgroups by fluorescence microscopy or flow cytometry, while quantitative polymerase chain reaction (PCR) strategies, in parallel with spectratyping techniques, were able to provide a rather low-resolution overview of the repertoire [12]. Despite these methods, for many years technical limitations made it difficult to create a comprehensive overview of real human TCR repertoires, until highly specific methods based on next generation sequencing (NGS) were developed, facilitating the parallel analysis of millions of TCR sequences. Nevertheless, it is still difficult to define a gold standard method, as every available method has its advantages and disadvantages.

In this article, we provide an overview of the currently available methodologies for TCR repertoire analysis, and we also describe the different aspects that a scientist should consider when choosing the appropriate method for the research question to be answered. We performed a small benchmark experiment that comprises some but not all available protocols for NGS-based immune repertoire analysis, and while our benchmark is by no means comprehensive and exhaustive, our results highlight some characteristics of the different methodologies and the approach in general that may serve as a guide for scientists that are interested and new in the field of immunogenetics.

Some publications focusing on immune repertoire profiling are already available and constitute an important source of information for any scientist interested in this research area [12,13,14,15,16,17,18,19]. However, direct comparisons and benchmarks of the most common methods are scarce.

Choosing the right starting material for TCR profiling

One of the most basic, yet important decisions a scientist should make when choosing a method for TCR analysis regards the starting material, i.e. whether to use genomic DNA (gDNA) or RNA. As discussed previously [14, 20], either starting material has advantages and disadvantages. The points in favour of gDNA are the higher stability and the presence of a single template per cell, which allows for better quantification of single TCR clones [21]. However, using gDNA does not provide any information on the expression level of the genes of interest and may lead to errors in the sequencing results due to introns, possible residuals of VDJ rearrangements and interfering priming sites found in the sample. When using RNA, quantification of single TCR clones is more challenging as a cell will contain multiple TCR transcripts. Many current methods, however, are designed for RNA as starting material, as the studied mRNA contains the final TCR products. Employing RNA potentially allows for sequencing of the entire J and V gene and it provides information about expression levels. Also, the quantity of starting material is a factor that needs to be taken into account. If only low quantities of starting material are available, this can be a limiting factor when selecting a suitable method, as some kits require a minimum input quantity and concentration of RNA or gDNA.

In general, due to the complexity of the target and the threat of batch effects that can affect the downstream data analysis, it is essential to ensure that the processing of all samples is as uniform as possible, for example by using the same concentration of starting material and trying to have a comparable number of reads for each sample [22].

TCR sequencing

In this article, we focused on the latest high-throughput sequencing (HTS) methods currently available for TCR repertoire profiling. The two options that we discuss are bulk sequencing of pooled immune cell populations or approaches allowing the analysis at the single cell level. We chose to concentrate on protocols for Illumina sequencing platforms, as this is the most widely established technology. Nonetheless, methods compatible with IonTorrent [23,24,25] and Roche 454 [26,27,28] exist.

While both methods for the analysis of single cells and cell populations (“bulk methods”) are available, we mainly focused on the latter, which are more commonly used to study TCR diversity and compare distinct repertoires in larger cohorts. The main disadvantage of bulk sequencing is that it can only provide information about the frequency of single TCR chains, but not their pairing. Single cell approaches are therefore becoming more and more important in immune repertoire studies, as they can accurately identify the pairs of the two TCR chains (αβ, γδ) at the cellular level, bringing repertoire analysis to a higher level of complexity. It is also the concrete chain pair that more accurately reflects the biological in vivo function. However, single cell sequencing approaches are currently more expensive, they cover often only a limited number of cells, as compared to bulk approaches, and, they require fresh material for the isolation and sorting of live cells, which is not always available especially when dealing with human diseased material. Some of the aspects we will consider in following paragraphs include the processing of RNA or gDNA samples prior to sequencing (library preparation) and choosing the adequate sequencing depth – the number of replicate reads necessary to efficiently detect the sequences of interest.

Bulk methods

There are different aspects to consider when choosing population-based TCR analysis methods. The methods we are going to compare differ in many aspects, from the type of starting material (gDNA or RNA) for the library preparation approach and the sequencing method. Some commercial companies offering immune repertoire analyses services are listed in Table 1. Our comparison mainly considers different library preparation methods, for which we will also discuss the different potential biases.

Table 1 Exemplary companies providing immune repertoire products and services

Choosing target sequences: Chains and CDR regions

Several companies offer library preparation and sequencing services for all TCR chains, but α chain and β chain remain the most common targets, as αβ T cells constitute the majority of the total T cell population [29]. Historically, the β chain was the main target studied due to its higher combinatorial potential compared to the α chain, which is due to the presence of the D gene component [14]. The β chain is also unique in each single cell, whereas it is possible that two α chains are expressed by the same cell, increasing the level of complexity [30]. γδ T-cell receptors are not widely studied, as γδ T cells only account for a small proportion of the total T cell population. The overall diversity of γδ TCRs is lower, compared to αβ TCRs, and there is an abundance bias based on which anatomical location is being analysed, as γδ T cells are found at higher frequency at mucosal sites. Therefore, they have been of less interest as peripheral blood samples are most widely studied [1]. PCR-based methods may amplify α and β chains simultaneously, but they are often separated and treated as two different samples in the last steps of the library preparation and during sequencing. This has been found to increase the precision and specificity of the outcome [31,32,33].

The CDR3 region is the preferential target of many TCR repertoire studies, due to its relevance for TCR-peptide interaction. To date, CDR1 and CDR2 have not attracted the same attention from the scientific community, because they do not directly interact with the antigen. However, CDR1 and CDR2 play an important role in making contact with the MHC molecule and thus influence the sensitivity and affinity of the TCR binding [34, 35]. Being aware of the sequence of the entire transcript, including CDR1 and CDR2, may be a great advantage for modelling the TCR structure and its binding properties. Not all methods are able to detect CDR1 and CDR2. This limitation applies especially to protocols using multiple primer sequences. Indeed, many allele-specific primers are designed in different positions of the V genes, often eliminating the possibility of sequencing outside CDR3.

Library preparation approach

We regard the library preparation approach as one of the key features to be considered when selecting a method. There are only a few techniques widely used for bulk analysis (Fig. 2). Most published methods are a variant of one of these approaches, the majority of which are PCR-based. Immune repertoire extrapolation from data generated through transcriptome sequencing is also feasible [36]; even though to date it has not been widely used and, given common transcriptome sequencing depths, this approach may be limited and reveal only a fraction of the TCR diversity as compared to a target-specific method.

Fig. 2
figure 2

Exemplary workflow of three principal methodologies for TCR library preparation. The figure depicts a simplified workflow of the library preparation procedure using multiplex PCR, targeted in-solution enrichment and 5’RACE-switch-oligo nested PCR. Multiplex PCR is suitable for both RNA and gDNA sequencing. Samples undergo cDNA synthesis and 1 or more PCR steps followed by adaptor ligation and sequencing. While the forward primers for cDNA synthesis are designed to cover all known V genes for both starting materials, the location and number of the reverse primers differs, due to introns in DNA. Target enrichment, also applicable to both gDNA and RNA, is preceded by a standard library preparation including fragmentation for gDNA or mRNA purification for RNA, followed by end-repairing, A-tailing and finally adaptor ligation. The enrichment of target sequences is then performed using RNA baits complementary to the sequence of interest. The RNA baits hybridize with molecules in the library, which are then retrieved using magnetic beads and can undergo further amplification before sequencing. Nested PCR based on the 5’RACE and switch-oligo approach (only for RNA) makes use of the incorporation of an adaptor molecule at the 5′ end of the cDNA during cDNA synthesis. The forward primer for a subsequent PCR is designed to bind to the 5′ adaptor sequence, while the reverse primer is designed to bind to the C-region of the transcript. Hence, only one primer pair is required to cover the complete spectrum of possible V genes. Subsequent nested PCRs performed in the same fashion may increase outcome specificity. Finally, adaptor ligation is performed. The procedures showed in this picture constitute only an example of the different available methods

Multiplex PCR

Given the diversity of the target, multiplex PCR approaches are among the most widespread, also offered as customizable service or kits by companies such as Adaptive Biotechnologies, BGI and iRepertoire [5, 31, 32, 37,38,39,40,41]. Primers for the J alleles or the constant region of the TCR α and β chains are used together with a mix of primers for all known V alleles. This results in a specific amplification of the TCR transcript across the CDR3 region. The multiplex approach can be used for both gDNA and RNA and the published protocols assure no cross-primer interference during amplification. However, this method cannot detect new V alleles variants due to the fixed set of primers used. Furthermore, multiplex PCR methods are subjected to amplification biases [42], which lead to better amplification of some alleles compared to others, thereby distorting the relative abundances of the resulting products. It is possible to correct for this kind of error by using a specific experimental design including adjustment of primer concentrations [43] and/or using molecular barcoding [44].

Target enrichment

A targeted enrichment method is available using e.g. Agilent’ RNA baits for capturing TCRs of αβ T cells. For starting material, gDNA or RNA is first processed with a standard sequencing library preparation kit (i.e. Illumina TruSeq or SureSelectXT from Agilent), followed by incubation of the samples with custom designed RNA baits. These RNA baits, which are complementary to the sequences of interest and tolerate a few different bases compared to the target, hybridize with the gDNA/cDNA target, allowing then for capturing it, and submitting the captured gDNA/cDNA to a further amplification step of the wanted sequences. This method, among other approaches, requires fewer PCR cycles and is thus less susceptible to PCR bias. Also α and β chains can be processed together, while it is suggested to separate the processing of the two chains for other methods in order to increase the quality of the outcome [45, 46].

5’RACE cDNA synthesis and nested PCR

For RNA samples, rapid amplification of 5′ complementary DNA ends (5’RACE) [47] employing the template-switch effect is becoming a gold standard for bulk TCR analysis [28, 48,49,50]. This method, marketed by Clonotech as “SMART” technology, relies on the terminal transferase activity of the reverse transcriptase enzyme, which incorporates additional nucleotides (usually dCTP) at the 3′ end of the cDNA molecule during the first strand synthesis reaction. A template-switch oligonucleotide containing an oligo(rG) sequence anchors to the non-template stretch of the first-strand cDNA, allowing the reverse transcriptase to switch templates and to continue replicating to the end of the oligonucleotide [51]. This enables the synthesis of cDNA strands containing the complete 5′ end of the mRNA, independent of the carried V allele, which enables capturing of all TCR variants present in the sample, provided that the integrity of the transcripts is conserved. Recently, Clonotech developed a commercial kit for TCR analysis using the afore-described template-switch technology. cDNA synthesis is carried out using primers against a small proportion of the target TCR mRNA transcript, the constant region. Consecutive PCRs may be carried out using a common adaptor as 5′ primer and constant region primers for the 3′ end. Ideally, these primers should be designed on nested sequences of the constant region in order to increase amplification specificity. Hence, only one primer set is required per reaction, avoiding the use of multiple primer sets and thus the associated amplification bias. The PCR products can be ligated to the appropriate sequencing adaptors and used for NGS sequencing.

General issues of TCR analysis

Despite the successful adoption and improvements of the mentioned approaches, every method based on PCR is still susceptible to a number of errors intrinsic to the particular technique, namely variable amplification efficiency due to differences in GC content, amplification stochasticity, template-switching and polymerase errors [52]. In addition, sequencing errors independent of the library preparation method used must always be taken in account. TCR sequencing is particularly vulnerable to sequencing errors, since a specific TCR may differ from another by only a single nucleotide. This fact makes it important to distinguish between PCR errors, sequencing errors and low frequency clonotypes. Thus, different techniques have been recently developed to overcome this issue. These methods comprise the usage of unique molecular identifiers (UMIs) introduced during cDNA synthesis to distinguish between single RNA molecules and minimize the impact of PCR amplification and sequencing errors [53], and specific algorithms for correction of this particular type of data [54].

Unique molecular identifiers

UMIs deserve special attention. The introduction of UMIs enabled the determination of the absolute count of RNA transcripts processed in a sample. UMI sequences, containing random nucleotides, are inserted into the template-switch oligonucleotide, which ligates with target molecules during cDNA synthesis, thereby uniquely barcoding every cDNA molecule in the sample with a different UMI. During data analysis, this allows retrieval of sequences originated from the same mRNA molecule after PCRs and sequencing. In a field like immune repertoire analysis, where target molecules may differ only by a single base, making the distinction between technical errors and biological differences is even more important; thus, using UMIs allows for a straightforward error correction. For these reasons their usage is increasing in immune repertoire profiling practice and herewith recommended [53, 55].

Sequencing platform and sequencing depth

Some sequencing platforms are more error prone than others, which means that a careful choice of the sequencing depth is needed in order to effectively manage the error rate, especially when not using molecular identifiers. A high sequencing depth allows for analysis of a more complete and complex repertoire [15], but deep sequencing is not always the best choice in immune repertoire analysis, depending on the purpose of the study. Disease-oriented analyses often look for highly expressed and clonally expanded TCRs. In this case, a superficial low-coverage screening of the immune repertoire may even be enough to catch the most common and expanded clonotypes in the sample. For this purpose, the Illumina MiSeq platform is commonly used, while the Illumina HiSeq is more often used for deep sequencing [14]. A general recommendation in order to achieve sufficient coverage for all sequences is to aim for at least 30,000 on-target reads, but ideally 100,000 reads per 10 ng of total starting RNA material (which relates to approximately 10,000 lymphocytes) should be performed [48].

Another issue that should be considered when selecting a sequencing platform is the diversity present in the library. Multiplex PCR and target enrichment based libraries have a higher diversity compared to 5’RACE-based libraries, which may all start with the same adaptor sequence at the 5′ end. Higher diversity makes sequencing with Illumina platforms easier. In case of low diversity samples, as 5′ RACE, some adjustments are necessary to improve the sequencing outcome. Including higher percentages of PhiX, which increases the diversity within the sequencing run, or the addition of random nucleotides in the used PCRs primers may overcome this issue [48].

Service vs. in-house methods

Samples processed by a company (see Table 1) undergo a standardized, robust and likely more reproducible workflow compared to a non-specialized laboratory. In addition, setting up a proper library preparation and bioinformatics data analysis pipeline is not trivial and having the analysis performed by experts in the field may help in saving time and effort. However, companies’ services are usually costlier than running the methods in-house if the instruments and personnel are readily available. Also, when thinking of data analysis, it is important to know prior to the experiment, what kind of data and format will be provided by the respective company. Some companies supply both raw sequencing data and even analysed data (e.g. BGI), or only raw data if the analysis is not covered by the service contract; others instead provide only the final analysed output and no raw data (e.g. Adaptive Biotechnologies). Ideally, both formats should be provided, in case one needs or wants to run additional analysis.

An additional option, when possible, may be to use a commercial kit for library preparation and then sequence in-house. This warrants complete oversight and control of sample processing and some companies provide a service of data analysis for self-performed sequencing (e.g. iRepertoire). Per-sample costs may be more expensive than for an in-house established protocol, but likely also less time consuming as the kits are standardized and contain thorough descriptions and advice for troubleshooting. “Open-source” protocols that are used in-house remain the more customizable option and enable for full control of every step of the process.

Data analysis

Over the past years different tools and strategies have been developed for immune repertoire analysis, of which some have been summarized in previous reviews [14, 15, 56,57,58]. Other methods, such as IMSEQ [59], TCRklass, iMonitor [60], LymAnalyzer [61] and RTRC [62], have however emerged since. A popular tool is MiXCR (previously MiTCR), developed by Bolotin et al. [33], which allows for a highly customizable analysis of both TCR and immunoglobulin sequences. This is the tool we chose for our analysis, as its parameters may be optimized for different data types, sources and desired outputs. Software specific for analysis of data containing UMIs are MIGEC [50] and pRESTO [63]. The tools listed above are mainly used for primary analysis, as the recovery of TCR sequences from raw data and successive clustering and annotation. LymAnalyzer additionally contains a feature for SNP calling and sequence mutation trees generation for IGs.

Further (secondary) data analysis of the immune repertoire classically involves the calculation of one or more diversity indices [64, 65]. Among the most widely used are the Shannon and the Simpson indices, as well as the Inverse-Simpson and the Gini indices. These differ for example in the consideration they give to factors as the species richness and the evenness of the dataset.

Another typical step of the analysis is the calculation of V and J gene usage in the different samples/datasets. The usage of different V and J genes is indeed not uniform. In literature, there are many examples of biased gene usage [10, 66]. A biased usage of specific genes may also be the result of alterations in the repertoire caused by diseases or other special conditions as organ transplantations.

Different tools for secondary TCR repertoire analysis and diversity estimation are available [67] and a list is available in the repertoire sequencing (Rep-seq) category of the Omic-tools community (https://omictools.com/rep-seq-category) [68]. Recent developments include VDJtools [69], which is capable of analysing outputs from the most common repertoire processing tools described above, and VDJviz [70], a webtool offering similar features as VDJtools. Another tool that provides TCR diversity measures and gene usage statistics computations is the R package “tcR”, which can be used to process the output files format of software as ImmunoSEQ [71, 72], IMSEQ, MiTCR, MiXCR, MIGEC and VDJtools [73]. Other approaches recently developed for estimation of TCR diversity are from Greiff et al., which creates a diversity profile using many diversity coefficients simultaneously [74], from Laydon et al., which introduces a new solution called DivE using rarefaction curves [75] and from Kaplinsky et al., which makes use of a maximum-likelihood based approach without assumptions of the complete repertoire clone richness [76].

When dealing with immune repertoire data, uniformity between samples is important. To this end, especially regarding diversity analyses, down- or re-sampling is a commonly used strategy to generate more comparable data. Similar data types may be easily encountered in the fields of ecology and metagenomics studies. Hence, data analysis packages intended for these disciplines may also be useful for immune repertoire analysis. An example for an interdisciplinary approach already used for TCR data is estimating the total diversity by using the “unseen species model” [77, 78]. The function for this model is provided for example by the “Vegan” R package, together with a series of common diversity measures and estimators [79]. More precise information on low and high complexity data analysis strategies are described in detail elsewhere [80,81,82,83].

Outlook: Single cell methods

Here, we want to briefly cite the main available approaches for single cell analysis of TCRs. There are two options commonly used for the analysis of the TCR repertoire from single cell suspensions: performing overall transcriptome sequencing and extrapolation of TCR information; or using approaches to specifically target the TCR transcript. To our knowledge, the most commonly used commercial workflow for TCR information extraction from transcriptome data is offered through the C1 system machine of the Fluidigm Corporation [84, 85].

Different methods have been established for specific targeting of the TCR transcript, of which some use multiplex primer sets [86]. For example, Han and colleagues used a multiplex PCR approach for TCR α and β chain analysis of FACS sorted cells [8]. Furthermore, they included additional non-interfering multiplex primers in their setup, enabling the parallel study of the expression level of phenotypic traits related to T cells, such as FOXP3, IL17A, TNF and others, thereby providing a more complete picture of the T cells of interest [8]. Another method, named “pairSEQ”, employs an experimental design that divides a sample into different subsets. It then uses combinatorics to evaluate unique TCR αβ chains in every subset [37].

Very recently, Wafergen Biosystems launched a new machine, the ICELL8 single-cell system. This machine may be used with a dedicated kit for TCR sequencing, based on the SMART technology of Takara Bio USA [87]. 10XGenomics, Inc. recently also launched a new dedicated kit for V(D)J analysis. A special note goes to the publicly available instrument and method published by McDaniel and colleagues, which can be used to process millions of cells and to potentially analyse all lymphocyte receptor chains (TCRαβ, TCRγδ, B cell heavy and light chains). This method provides instructions for the construction of a dedicated device, which makes use of an extended concept of the emulsion PCR technology [88], capturing individual molecules using primer-covered beads in droplets within an oil phase, and performing PCR reactions for each bead [89] [90]. This technology, applied to TCR analysis, was previously also published by Turchaninova et al. [91].

For a more detailed summary of single cell approaches to study the immune system, reviews on the topic have been published by Chattopadhyay et al. [92] and Proserpio et al. [93].

Results

We performed experiments using two of the methods presented above as the most common approaches for bulk immune repertoire sequencing, namely multiplex PCR and 5’RACE-based PCR, using total RNA as starting material. In addition, we compared our results with results provided by the BGI immune repertoire sequencing service, based on a multiplex PCR approach starting from gDNA.

Here, we present some considerations based on our own experiments and analyses, which may help to better understand the methods described until now.

Replicate correlations

iRepertoire® library preparation and 5’RACE-based PCR were performed in duplicates. Analysis of 5’RACE-based PCR results was performed in parallel for UMI corrected versus non-UMI corrected data. Previous studies have demonstrated the difficulty of detecting the entire TCR diversity of a sample, as it can vary consistently even between very close anatomical locations or between time points. We thus decided to use replicates of the same sample to better assess the stability of TCR diversity, especially when using superficial sequencing [11].

As mentioned before, quantity of starting material and sequencing depth play a major role in defining the extent of the captured TCR diversity. A popular recommendation is to use 100,000 reads to more than efficiently cover the diversity in 10 ng of total RNA [48]. We opted for a superficial sequencing approach (1 million reads for 500 ng total RNA), aiming to detect only the most abundant TCR clonotypes and to determine how many were identified in both duplicates, despite the use of a sequencing depth which likely does not cover the complete diversity present in the sample. For each method (iRepertoire, 5’RACE and 5’RACE UMI corrected) we compared the percentage of clonotypes shared between duplicates (Table 2), assessing their capacity of capturing TCR diversity.

Table 2 Percentages of CDR3 nucleotide sequences detected in both duplicates of the same method

When comparing the total detected repertoire, only 17–52% of observed clonotypes were shared between duplicates (Table 2). However, when we compared only the most abundant clonotypes of each duplicate (the 300, 100, 50 and 20 most abundant clones), we observed a drastic increase in the percentage of shared clonotypes (50–80%). The highest overlap was observed when comparing the most frequent 50 to 20 clonotypes. Thus, despite superficial sequencing depth, we could successfully detect the majority of the most abundant clonotypes in the samples. The increment in the percentage of shared species was even more marked for the UMI corrected data, demonstrating the usefulness of this strategy for error correction. Shared clonotype percentages were comparable between α and β chains.

Out of frame and stop codon containing sequences

During TCR data analysis it is common to encounter non-functional CDR3 sequences, which are out of frame or stop codon containing sequences. MiXCR labels these sequences with particular symbols, making it possible to exclude these sequences from the clonotype list.

We analysed the percentages of these non-functional sequences in the data obtained from the iRepertoire® kit or 5′ RACE based PCR and in the data provided by BGI, which we re-processed with MiXCR (Table 3). Our results imply that the percentage of both out of frame and stop codon containing sequences is higher in α chain results as compared to β chains, and that out of frame sequences are more commonly detected, as compared to stop codons. iRepertoire results appeared to contain the least non-functional CDR3 sequences. Samples analysed by BGI contained a significantly higher percentage of non-functional CDR3 sequences as compared to other methods, which was to be expected, as gDNA was used as starting material and the sequencing depth was significantly higher. As anticipated, the percentages of non-functional sequences decreased after correcting for UMIs. Percentages were comparable between the two patients we analysed.

Table 3 Percentages of detected CDR3 sequences that are out of frame or contain stop codons

Method comparison

Here, we compared results from three different methods performed with the same patient samples, namely BGI service (using gDNA), iRepertoire® kit and the in-house established 5’RACE-based PCR (using RNA). In addition, for 5’RACE data we compared both UMI-corrected and not corrected data. After clonotype grouping and export, all datasets were filtered by retaining only clonotypes detected with two or more counts. To make data sets with different sequencing depths comparable, we decided to assess only the most abundant clonotypes detected by each method. To set a threshold, we determined which methodology provided the lowest total number of detected TCR clonotypes. This applied to UMI-corrected 5′ RACE based PCR (1300–2400 species, depending on chain and patient), which was anticipated due to the strict filtering steps applied during correction. We therefore only considered the most abundant species detected by other methods above this threshold and excluded any sequences below it. Our aim was to determine if high frequency TCR clonotypes and their relative abundancies could consistently be captured using these different library preparation approaches.

When analysing the overlap of TCR clonotypes detected by the different methods, we found that less than 10% were captured by all four methods (Fig. 3, Additional file 1: Figure S1). These species were, however, detected at high abundance in all methods. The majority of clonotypes detected by one of the methods were uniquely detected by that particular method (up to 75%), at least among the highly abundant clones. As anticipated, the strongest overlap in species was observed between 5′ RACE-based PCR and UMI-corrected 5′ RACE-based PCR. The percentage of sequences commonly captured by three or all methods increased when considering only the highly abundant clonotypes (Additional file 1: Figure S2). This is in concordance with the observations made for replicate correlations (Table 2). However, clones common to all methods were 19–25% when comparing the top 100 clones and overlap not higher than 37%. These results demonstrate not only the diversity of the TCR repertoire, but also how clonotype abundance within the same sample may vary when using different techniques and correction methods. Indeed, clonotypes detected as highly abundant by only one method, may still be detected in others, but at low counts, resulting in inconsistent information about relative species frequencies. Using strict criteria of UMI correction lead to an involuntary loss of information, reflected by the fact that not UMI corrected clonotypes detected by multiple methods were lost after UMI correction. However, less stringent and conservative methods for UMI correction are available [50] and have recently been discussed by Smith, Heger and Sudbery [94].

Fig. 3
figure 3

Venn diagram showing the overlap between the most abundant TCR sequences detected by each method. The threshold was defined by the method which detected fewest TCR species (UMI-corrected 5’RACE-based PCR). The diagrams show the number and relative frequency of TCR sequences detected by only one up to all four methods. Sequences that were found only by one method may still be detected by other methods, but may not appear in the most abundant species and are thus not represented here. Data are shown for one patient for both (a) α and (b) β chains

Exemplary analysis of gene usage and diversity

As described in the data analysis paragraph, we assessed the V and J gene usage in our samples for α and β chain.

The analysis showed a method-dependent gene usage bias for both the V and the J genes of α and β chain (Fig. 4 and additional file 1: Figure S3, respectively). Relative differences in gene usage between the two patients seem to be conserved among different methods, in particular for the most used genes. For example, TRAV8–6 appears to be more used in patient 1 compared to patient 2 and vice versa for TRAV1–2. However, and most importantly, a method-dependent bias in the gene usage was observed. A possible explanation for such differences across methods could be the use of different primers during library preparation. We believe this is a critical observation and an aspect that should be considered prior to attempting to compare results derived from different technologies.

Fig. 4
figure 4

Alpha chain Variable (V) gene usage among methods and replicates. The heat map shows the gene usage proportion inside each sample for the V genes listed at the bottom on the figure. Each of the three parts of the heat map is representative of one of the methods. The samples are described by the patients “P1/2”, the method “iRepertoire/BGI/5’RACE/5’RACE + UMI” and the replicate “1 or 2” or “0” in case of BGI which does not include replicates. The figure was generated using the “ggplot2” R package

In addition to gene usage, we also performed an exemplary analysis to study the diversity among different samples and, in this particular case, methods. Results are shown in Fig. 5.

Fig. 5
figure 5

Dissimilarity dendrogram for alpha chain. The distance between patients, methods, and replicates, calculated using the Morisita index, is represented by the dendrograms. The samples are described by the patients “P1/2”, the method “iRepertoire/BGI/5’RACE/5’RACE + UMI” and the replicate “1 or 2”. The figure was generated using the “ggplot2” and the “ggdendro” R packages

As expected, samples from the same patient cluster together, as well as replicates from the same method. iRepertoire and 5’RACE PCR results seem to be the most distant from each other, while BGI results appears to cluster more with 5’RACE results, compared to iRepertoire (Fig. 5, Additional file 1: Figure S4).

Discussion

From our own experience, we learned that different methods are effective in TCR species identification, but to a different extent.

We prefer RNA over gDNA as starting material, because, even if more unstable, RNA contains the final TCR transcript, it does not add noise to the results due to incomplete VDJ-recombination products, and it allows for the use of the 5’RACE method. This is the approach we suggest for library preparation, for two main reasons: (1) it bypasses the amplification biases associated with multiplex PCR and (2) it is suitable for use with UMIs.

In general, we consider low abundance TCR counts as unstable, unless using particularly effective error correction strategies, such as UMIs. These error correction strategies need to be handled with care in order to use the right filtering criteria that fit the needs of the study. If the purposes of the project include analyses of low abundant species, we suggest a strict UMI filtering.

We did not observe a higher overlap with other methods while comparing UMI vs non UMI corrected data. Nevertheless, the concept behind UMI usage remains valid and we believe it should become mandatory in TCR analysis methodologies, in order to avoid PCR distortions, which currently increase the complexity of comparing samples and methods.

To avoid setting up the methods in-house, different kits are commercially available, with or without sequencing service (Table 1). When sequencing in-house, it is important to select the appropriate sequencing depth in relation to the quantity of starting material. More superficial sequencing approaches are cheaper, but may not be able to reflect the entire diversity and are not suitable for rare TCR analysis. However, superficial sequencing may suffice when only abundant clonotypes are of interest. Independent of which method is chosen, one should keep in mind that a uniform quantity and quality of starting material is essential for the positive outcome of any experiment.

Interestingly, we were able to detect some TCR species using low-coverage sequencing, which were not detected within the BGI data. We believe that this may be due to the difficulty of capturing the entire TCR diversity, even with deep sequencing, considering the extremely high variability of the target. Nevertheless, highly abundant species showed a significant overlap between methods. We think that deep sequencing is preferable for studies which aim to extensively analyse the repertoire of cohorts and populations, while low-coverage sequencing might be preferential for studies directed at the identification of already known or abundant clonotypes. Also, due to the high repertoire diversity, performing biological replicates when possible can help reinforce the analysis findings.

Conclusions

Table 4 summarizes the advantages and disadvantages of every method that we employed in our studies.

Table 4 Advantages and disadvantages of the tested techniques

It is important to consciously select a method by keeping in mind strengths and weaknesses of each approach, as well as the goal of the scientific project which the method will be applied to. It is also crucial to be careful when comparing relative abundances in data obtained with different methods, as they may be affected by method specific biases as PCR amplification efficiency discrepancies due to different primers.

Many of the mentioned methods may also be applicable to B-cell receptor analyses and as immune repertoire studies are increasing in number and complexity, we believe that an educated choice of methodology is one of the most important steps to achieve the desired results in the growing field of “immunogenetics”.

Methods

Sample selection

For benchmarking we used two explanted liver tissue samples from patients with primary sclerosing cholangitis (PSC), which were previously included in a study analysing PSC-related TCR repertoires [95]. PSC is a chronic inflammatory disease of the liver affecting the intra- and extra- hepatic bile ducts. PSC is characterized by high T cell infiltration; thus, it provides a suitable system to study T-cell receptor (TCR) repertoire signatures. To date, it is completely unknown what causes PSC and which is the triggering and/or driving antigen. Detecting disease-associated TCR signatures would be an important step towards the identification of the triggering/driving antigen. Therefore, systematic TCR profiling experiments in the organ of interest are of great importance to further understand the immunogenetics of PSC and of other immune-mediated diseases of unknown etiology.

gDNA and RNA were isolated simultaneously from whole fresh-frozen disease-affected tissue using the AllPrep DNA/RNA Mini Kit from Qiagen.

iRepertoire ®

As an exemplary method for RNA-based multiplex PCR we selected the iRprofile kit from iRepertoire Inc. We processed both samples in duplicates using 500 ng of total RNA. This kit is highly customizable and offers options for both T and B cell receptor sequencing, different receptor chains (αβ, γδ), gDNA or RNA, mouse or human, and sequencing platforms (Illumina, Roche 454). iRprofile contains separate reactions mixes for α and β chains, each uniquely barcoded. The protocol consists of two consecutive PCRs using multiplex primers specific for the V and J genes of the mentioned chains. The resulting products are then ready for sequencing. The sequencing data may then be sent to the company for demultiplexing. Basic data analysis is also performed and demultiplexed data may be requested as FASTA files.

5’RACE

The 5’RACE method we chose is an in-house adaptation of the protocol published by Mamedov et al. 2013 [48]. Briefly, the method entails 5’RACE-based cDNA synthesis using a 5′-template switch adapter containing 12 random nucleotides forming a UMI, followed by two consecutive nested PCRs. During the second PCR, Illumina adaptors are incorporated [96]. Consequently, custom sequencing primers were used in the following MiSeq run. α and β chains were amplified in the same reaction for cDNA synthesis and PCR 1, and they were separated and uniquely barcoded during PCR 2. Samples were processed in duplicates using 500 ng of total RNA.

BGI

5 μg of gDNA were sent to BGI for both α and β chain immune repertoire analysis (2.5 μg/chain) based on multiplex PCR. BGI performed basic data processing which included data filtering, removal of adapter contamination and low quality reads from raw reads and elimination of sequencing background. Alignment to V/D/J gene segments was carried out separately in IMGT database, and data were realigned for best results. BGI also carried out structural analysis which included CDR sequence and base composition, V/D/J recombination insertion and deletion. Data analysis included immune repertoire profiling and differential analysis of diversity between samples, differential expression analysis of clones between samples and differential expression analysis of clones between groups. In addition to complete data analysis all raw FASTQ files were provided by BGI.

Sequencing and data analysis

RNA sequencing was performed on an Illumina MiSeq 250PE. The sequencing platform used by BGI was Hiseq2000 100PE.

For data analysis, we used MiXCR (version 2.1.1), obtaining a ranked table of clonotypes including relative species abundances, nucleotide and amino acid CDR3 sequences and respective VDJ alleles as output. Sequences containing the same UMI were grouped under the same UMI signature. For each UMI, only the most abundant sequence was selected, while the others were considered PCR or sequencing errors (script used for UMI filtering is available as Additional file 2). UMI filtered results represent absolute mRNA transcript relative abundances in the original sample. BGI provided data obtained using an older version of the MiXCR software (MiTCR), which is why we reprocessed the raw data with the same version of the software we used for iRepertoire and 5’RACE data. Analysis parameters were optimized based on starting material.

Gene usage analysis was performed using the “geneusage” function of the “tcR” R package, while diversity analysis was performed using the “vegdist” function of the “Vegan” R package.

Abbreviations

5’RACE:

5′ rapid amplification of cDNA ends

BCR:

B-cell receptor

C:

Constant

CDR:

Complementary determining region

D:

Diversity

dCTP:

Deoxycytidine triphosphate

FACS:

Fluorescence-activated cell sorting

FOXP3:

Forkhead box P3

gDNA:

Genomic DNA

HTS:

High-throughput sequencing

IL17A:

Interleukin-17A

J:

Joining

MHC:

major histocompatibility complex

Mplex-PCR:

Multiplex PCR

mRNA:

messenger RNA

NGS:

Next generation sequencing

PCR:

Polymerase chain reaction

PE:

Paired-end

PSC:

Primary sclerosing cholangitis

TCR:

T-cell receptor

TNF:

Tumour necrosis factor

UMI:

Unique molecular identifier

V:

Variable

References

  1. Presti E Lo, Dieli F, Meraviglia S. Tumor-Infiltrating γδ T Lymphocytes: Pathogenic Role, Clinical Significance, and Differential Programing in the Tumor Microenvironment. Front Immunol. 2014;5. doi:10.3389/fimmu.2014.00607.

  2. Miles JJ, Douek DC. Price D a. Bias in the αβ T-cell repertoire: implications for disease pathogenesis and vaccination. Immunol Cell Biol. 2011;89:375–87. doi:10.1038/icb.2010.139.

    Article  CAS  Google Scholar 

  3. Burtrum DB, Kim S, Dudley EC, Hayday AC, Petrie HT. TCR gene recombination and alpha beta-gamma delta lineage divergence: productive TCR-beta rearrangement is neither exclusive nor preclusive of gamma delta cell development. J Immunol. 1996;157(10):4293-6.

  4. Turner SJ, Doherty PC, McCluskey J, Rossjohn J. Structural determinants of T-cell receptor bias in immunity. Nat Rev Immunol. 2006;6:883–94. doi:10.1038/nri1977.

    Article  CAS  Google Scholar 

  5. Muraro PA, Robins H, Malhotra S, Howell M, Phippard D, Desmarais C, et al. T cell repertoire following autologous stem cell transplantation for multiple sclerosis. J Clin Invest. 2014;124:1168–72. doi:10.1172/JCI71691.

    Article  CAS  Google Scholar 

  6. Bai X, Zhang Q, Wu S, Zhang X, Wang M, He F, et al. Characteristics of tumor infiltrating lymphocyte and circulating lymphocyte repertoires in pancreatic cancer by the sequencing of T cell receptors. Sci Rep. 2015;5:13664. doi:10.1038/srep13664.

    Article  Google Scholar 

  7. Gros A, Robbins PF, Yao X, Li YF, Turcotte S, Tran E, et al. PD-1 identifies the patient-specific CD8+ tumor-reactive repertoire infiltrating human tumors. J Clin Invest. 2014;124:2246–59. doi:10.1172/JCI73639.

    Article  CAS  Google Scholar 

  8. Han A, Glanville J, Hansmann L, Davis MM. Linking T-cell receptor sequence to functional phenotype at the single-cell level. Nat Biotechnol. 2014;32:684–92. doi:10.1038/nbt.2938.

    Article  CAS  Google Scholar 

  9. Laydon DJ, Bangham CRM, Asquith B. Estimating T-cell repertoire diversity: limitations of classical estimators and a new approach. Philos Trans R Soc Lond Ser B Biol Sci. 2015;370:20140291. doi:10.1098/rstb.2014.0291.

    Article  Google Scholar 

  10. Freeman JD, Warren RL, Webb JR, Nelson BH, Holt RA. Profiling the T-cell receptor beta-chain repertoire by massively parallel sequencing. Genome Res. 2009;19:1817–24. doi:10.1101/gr.092924.109.

    Article  CAS  Google Scholar 

  11. Warren RL, Freeman JD, Zeng T, Choe G, Munro S, Moore R, et al. Exhaustive T-cell repertoire sequencing of human peripheral blood samples reveals signatures of antigen selection and a directly measured repertoire size of at least 1 million clonotypes. Genome Res. 2011;21:790–7. doi:10.1101/gr.115428.110.

    Article  CAS  Google Scholar 

  12. Six A, Mariotti-Ferrandiz ME, Chaara W, Magadan S, Pham H-P, Lefranc M-P, et al. The past, present, and future of immune repertoire biology - the rise of next-generation repertoire analysis. Front Immunol. 2013;4:413. doi:10.3389/fimmu.2013.00413.

    Article  Google Scholar 

  13. Nikolich-Zugich J, Slifka MK, Messaoudi I. The many important facets of T-cell repertoire diversity. Nat Rev Immunol. 2004;4:123–32.

    Article  CAS  Google Scholar 

  14. Woodsworth DJ, Castellarin M, Holt R a. Sequence analysis of T-cell repertoires in health and disease. Genome Med. 2013;5:98. doi:10.1186/gm502.

    Article  Google Scholar 

  15. Calis JJA, Rosenberg BR. Characterizing immune repertoires by high throughput sequencing: strategies and applications. Trends Immunol. 2014;35:581–90.

    Article  CAS  Google Scholar 

  16. Benichou J, Ben-Hamo R, Louzoun Y, Efroni S. Rep-Seq: uncovering the immunological repertoire through next-generation sequencing. Immunology. 2012;135:183–91. doi:10.1111/j.1365-2567.2011.03527.x.

    Article  CAS  Google Scholar 

  17. Han Y, Li H, Guan Y, Huang J. Immune repertoire: a potential biomarker and therapeutic for hepatocellular carcinoma. Cancer Lett. 2016;379:206–12. doi:10.1016/j.canlet.2015.06.022.

    Article  CAS  Google Scholar 

  18. Baum PD, Venturi V, Price D a. Wrestling with the repertoire: the promise and perils of next generation sequencing for antigen receptors. Eur J Immunol. 2012;42:2834–9. doi:10.1002/eji.201242999.

    Article  CAS  Google Scholar 

  19. Hou X-L, Wang L, Ding Y-L, Xie Q, Diao H-Y. Current status and recent advances of next generation sequencing techniques in immunological repertoire. Genes Immun. 2016;17:153–64. doi:10.1038/gene.2016.9.

    Article  CAS  Google Scholar 

  20. Six A, Mariotti-Ferrandiz ME, Chaara W, Magadan S, Pham HP, Lefranc MP, et al. The past, present, and future of immune repertoire biology - the rise of next-generation repertoire analysis. Front Immunol. 2013;4 NOV:1–16.

    Google Scholar 

  21. Dziubianau M, Hecht J, Kuchenbecker L, Sattler A, Stervbo U, Rödelsperger C, et al. TCR repertoire analysis by next generation sequencing allows complex differential diagnosis of T cell-related pathology. Am J Transplant. 2013;13:2842–54.

    Article  CAS  Google Scholar 

  22. Mamedov IZ, Britanova OV, Zvyagin IV, Turchaninova MA, Bolotin DA, Putintseva EV, et al. Preparing unbiased T-cell receptor and antibody cDNA libraries for the deep next generation sequencing profiling. Front Immunol. 2013;4:456. doi:10.3389/fimmu.2013.00456.

    Article  Google Scholar 

  23. Tamura K, Hazama S, Yamaguchi R, Imoto S, Takenouchi H, Inoue Y, et al. Characterization of the T cell repertoire by deep T cell receptor sequencing in tissues and blood from patients with advanced colorectal cancer. Oncol Lett. 2016;11:3643–9. doi:10.3892/ol.2016.4465.

    Google Scholar 

  24. Yew PY, Alachkar H, Yamaguchi R, Kiyotani K, Fang H, Yap KL, et al. Quantitative characterization of T-cell repertoire in allogeneic hematopoietic stem cell transplant recipients. Bone Marrow Transplant. 2015;50:1227–34. doi:10.1038/bmt.2015.133.

    Article  CAS  Google Scholar 

  25. Fang H, Yamaguchi R, Liu X, Daigo Y, Yew PY, Tanikawa C, et al. Quantitative T cell repertoire analysis by deep cDNA sequencing of T cell receptor α and β chains using next-generation sequencing (NGS). Oncoimmunology. 2014;3:e968467. doi:10.4161/21624011.2014.968467.

    Article  Google Scholar 

  26. Al-Hussaini M, Rettig MP, Ritchey JK, Karpova D, Uy GL, Eissenberg LG, et al. Targeting CD123 in acute myeloid leukemia using a T-cell-directed dual-affinity retargeting platform. Blood. 2016;127:122–31. doi:10.1182/blood-2014-05-575704.

    Article  CAS  Google Scholar 

  27. O’Connell AE, Volpi S, Dobbs K, Fiorini C, Tsitsikov E, de Boer H, et al. Next generation sequencing reveals skewing of the T and B cell receptor repertoires in patients with wiskott-Aldrich syndrome. Front Immunol. 2014;5:340. doi:10.3389/fimmu.2014.00340.

    Google Scholar 

  28. Kitaura K, Shini T, Matsutani T, Suzuki R. A new high-throughput sequencing method for determining diversity and similarity of T cell receptor (TCR) α and β repertoires and identifying potential new invariant TCR α chains. BMC Immunol. 2016;17:38. doi:10.1186/s12865-016-0177-5.

    Article  Google Scholar 

  29. Holtmeier W, Kabelitz D. γδ T cells link innate and adaptive immune responses focus on human Vγ9/Vδ2 and Vδ1 T cells. Chem Immunol Allergy Basel. 2005;86:151–83.

    Article  CAS  Google Scholar 

  30. Padovan E, Casorati G, Dellabona P, Meyer S, Brockhaus M, Lanzavecchia A. Expression of two T cell receptor alpha chains: dual receptor T cells. Science. 1993;262:422–4. doi:10.1126/science.8211163.

    Article  CAS  Google Scholar 

  31. Zeng G, Huang Y, Huang Y, Lyu Z, Lesniak D, Randhawa P. Antigen-specificity of T-cell infiltrates in biopsies with T-cell mediated rejection and BK Polyomavirus Viremia: analysis by next generation sequencing. Am J Transplant. 2016; doi:10.1111/ajt.13911.

  32. Sims JS, Grinshpun B, Feng Y, Ung TH, Neira JA, Samanamud JL, et al. Diversity and divergence of the glioma-infiltrating T-cell receptor repertoire. Proc Natl Acad Sci. 2016;113:E3529–37. doi:10.1073/pnas.1601012113.

    Article  CAS  Google Scholar 

  33. Bolotin D A, Poslavsky S, Mitrophanov I, Shugay M, Mamedov IZ, Putintseva E V, et al. MiXCR: software for comprehensive adaptive immunity profiling. Nat Methods 2015;12:380–381. doi:10.1038/nmeth.3364.

  34. Lynch JN, Donermeyer DL, Weber KS, Kranz DM, Allen PM. Subtle changes in TCRα CDR1 profoundly increase the sensitivity of CD4 T cells. Mol Immunol. 2013;53:283–94. doi:10.1016/j.molimm.2012.08.020.

    Article  CAS  Google Scholar 

  35. Birnbaum ME, Mendoza JL, Sethi DK, Dong S, Glanville J, Dobbins J, et al. Deconstructing the peptide-MHC specificity of t cell recognition. Cell. 2014;157:1073–87.

    Article  CAS  Google Scholar 

  36. Brown SD, Raeburn LA, Holt RA, Woodsworth D, Castellarin M, Holt R, et al. Profiling tissue-resident T cell repertoires by RNA sequencing. Genome Med. 2015;7:125. doi:10.1186/s13073-015-0248-x.

    Article  CAS  Google Scholar 

  37. Howie B, Sherwood AM, Berkebile AD, Berka J, Emerson RO, Williamson DW, et al. High-throughput pairing of T cell receptor and sequences. Sci Transl Med. 2015;7:301ra131. doi:10.1126/scitranslmed.aac5624.

    Article  Google Scholar 

  38. Lord J, Chen J, Thirlby RC, Sherwood AM, Carlson CS. T-cell receptor sequencing reveals the Clonal diversity and overlap of colonic Effector and FOXP3+ T cells in ulcerative colitis. Inflamm Bowel Dis. 2015;21:19–30. doi:10.1097/MIB.0000000000000242.

    Article  Google Scholar 

  39. Madi A, Shifrut E, Reich-Zeliger S, Gal H, Best K, Ndifon W, et al. T-cell receptor repertoires share a restricted set of public and abundant CDR3 sequences that are associated with self-related immunity. Genome Res. 2014; doi:10.1101/gr.170753.113.

  40. Rossetti M, Spreafico R, Consolaro A, Leong JY, Chua C, Massa M, et al. TCR repertoire sequencing identifies synovial Treg cell clonotypes in the bloodstream during active inflammation in human arthritis. Ann Rheum Dis. 2016;:annrheumdis-2015-208992. doi:10.1136/annrheumdis-2015-208992.

  41. Ndifon W, Gal H, Shifrut E, Aharoni R, Yissachar N, Waysbort N, et al. Chromatin conformation governs T-cell receptor Jβ gene segment usage. Proc Natl Acad Sci U S A. 2012;109:15865–70. doi:10.1073/pnas.1203916109.

    Article  CAS  Google Scholar 

  42. Okino ST, Kong M, Sarras H, Wang Y. Evaluation of bias associated with high-multiplex, target-specific pre-amplification. Biomol Detect Quantif. 2016;6:13–21.

    Article  CAS  Google Scholar 

  43. Carlson CS, Emerson RO, Sherwood AM, Desmarais C, Chung M-W, Parsons JM, et al. Using synthetic templates to design an unbiased multiplex PCR assay. Nat Commun. 2013;4:20216–21. doi:10.1038/ncomms3680.

    Article  Google Scholar 

  44. Peng Q, Vijaya Satya R, Lewis M, Randad P, Wang Y, Mortazavi A, et al. Reducing amplification artifacts in high multiplex amplicon sequencing by using molecular barcodes. BMC Genomics. 2015;16:589. doi:10.1186/s12864-015-1806-8.

    Article  Google Scholar 

  45. Linnemann C, Heemskerk B, Kvistborg P, Kluin RJC, Bolotin D a, Chen X, et al. High-throughput identification of antigen-specific TCRs by TCR gene capture. Nat Med. 2013;19:1534–41. doi:10.1038/nm.3359.

    Article  CAS  Google Scholar 

  46. Linnemann C, Mezzadra R, Schumacher TNM. TCR repertoires of intratumoral T-cell subsets. Immunol Rev. 2014;257:72–82.

    Article  CAS  Google Scholar 

  47. Rapid amplification of 5′ complementary DNA ends (5′ RACE). Nat Methods. 2005;2:629–30. doi:10.1038/nmeth0805-629.

  48. Mamedov IZ, Britanova O V., Zvyagin I V., Turchaninova M a., Bolotin D a., Putintseva E V., et al. Preparing unbiased T-cell receptor and antibody cDNA libraries for the deep next generation sequencing profiling. Front Immunol. 2013;4:456.

  49. Matz M, Shagin D, Bogdanova E, Britanova O, Lukyanov S, Diatchenko L, et al. Amplification of cDNA ends based on template-switching effect and step-out PCR. Nucleic Acids Res. 1999;27:1558–60.

    Article  CAS  Google Scholar 

  50. Shugay M, Britanova OV, Merzlyak EM, Turchaninova MA, Mamedov IZ, Tuganbaev TR, et al. Towards error-free profiling of immune repertoires. Nat Methods. 2014;11:653–5. doi:10.1038/nmeth.2960.

    Article  CAS  Google Scholar 

  51. Zhu YY, Machleder EM, Chenchik A, Li R, Siebert PD. Reverse transcriptase template switching: a SMARTtm approach for full-length cDNA library construction. BioTechniques. 2001;30:892–7.

    CAS  Google Scholar 

  52. Kebschull JM, Zador AM. Sources of PCR-induced distortions in high-throughput sequencing data sets. Nucleic Acids Res. 2015;43:e143. doi:10.1093/nar/gkv717.

    Article  Google Scholar 

  53. Kivioja T, Vähärautio A, Karlsson K, Bonke M, Enge M, Linnarsson S, et al. Counting absolute numbers of molecules using unique molecular identifiers. Nat Methods. 2011;9:72–4. doi:10.1038/nmeth.1778.

    Article  Google Scholar 

  54. Bolotin D a, Mamedov IZ, Britanova OV, Zvyagin IV, Shagin D, Ustyugova SV, et al. Next generation sequencing for TCR repertoire profiling: platform-specific features and correction algorithms. Eur J Immunol. 2012;42:3073–83.

    Article  CAS  Google Scholar 

  55. Egorov ES, Merzlyak EM, Shelenkov AA, Britanova OV, Sharonov GV, Staroverov DB, et al. Quantitative profiling of immune repertoires for minor lymphocyte counts using unique molecular identifiers. J Immunol. 2015;194:6155–63. doi:10.4049/jimmunol.1500215.

    Article  CAS  Google Scholar 

  56. Kidd B a, Peters L a, Schadt EE, Dudley JT. Unifying immunology with informatics and multiscale biology. Nat Immunol. 2014;15:118–27. doi:10.1038/ni.2787.

    Article  CAS  Google Scholar 

  57. Li S, Lefranc M-P, Miles JJ, Alamyar E, Giudicelli V, Duroux P, et al. IMGT/HighV QUEST paradigm for T cell receptor IMGT clonotype diversity and next generation repertoire immunoprofiling. Nat Commun. 2013;4 doi:10.1038/ncomms3333.

  58. Thomas N, Heather J, Ndifon W, Shawe-Taylor J, Chain B. Decombinator: a tool for fast, efficient gene assignment in T-cell receptor sequences using a finite state machine. Bioinformatics. 2013;29:542–50.

    Article  CAS  Google Scholar 

  59. Kuchenbecker L, Nienen M, Hecht J, Neumann a. U, Babel N, Reinert K, et al. IMSEQ--a fast and error aware approach to immunogenetic sequence analysis. Bioinformatics. 2015; May:1–9. doi:10.1093/bioinformatics/btv309.

  60. Zhang W, Du Y, Su Z, Wang C, Zeng X, Zhang R, et al. Imonitor: a robust pipeline for TCR and BCR repertoire analysis. Genetics. 2015;201:459–72.

    Article  CAS  Google Scholar 

  61. Yu Y, Ceredig R, Seoighe C. LymAnalyzer: a tool for comprehensive analysis of next generation sequencing data of T cell receptors and immunoglobulins. Nucleic Acids Res. 2016;44:e31. doi:10.1093/nar/gkv1016.

    Article  Google Scholar 

  62. Gerritsen B, Pandit A, Andeweg AC, de Boer RJ. RTCR: a pipeline for complete and accurate recovery of T cell repertoires from high throughput sequencing data. Bioinformatics. 2016;32:3098–106. doi:10.1093/bioinformatics/btw339.

    Article  Google Scholar 

  63. Vander Heiden JA, Yaari G, Uduman M, Stern JNH, O’Connor KC, Hafler DA, et al. pRESTO: a toolkit for processing high-throughput sequencing raw reads of lymphocyte receptor repertoires. Bioinformatics. 2014;30:1930–2. doi:10.1093/bioinformatics/btu138.

    Article  CAS  Google Scholar 

  64. Morris EK, Caruso T, Buscot F, Fischer M, Hancock C, Maier TS, et al. Choosing and using diversity indices: insights for ecological applications from the German biodiversity Exploratories. Ecol Evol. 2014;4:3514–24. doi:10.1002/ece3.1155.

    Article  Google Scholar 

  65. Barwell LJ, Isaac NJB, Kunin WE. Measuring β-diversity with species abundance data. J Anim Ecol. 2015;84:1112–22. doi:10.1111/1365-2656.12362.

    Article  Google Scholar 

  66. Miles JJ, Douek DC, Price DA. Bias in the ab T-cell repertoire: implications for disease pathogenesis and vaccination. Immunol Cell Biol. 2011;89:375–87.

    Article  CAS  Google Scholar 

  67. Venturi V, Kedzierska K, Turner SJ, Doherty PC, Davenport MP. Methods for comparing the diversity of samples of the T cell receptor repertoire. J Immunol Methods. 2007;321:182–95.

    Article  CAS  Google Scholar 

  68. Henry VJ, Bandrowski AE, Pepin A-S, Gonzalez BJ, Desfeux A. OMICtools: an informative directory for multi-omic data analysis. Database (Oxford). 2014;2014:bau069. doi:10.1093/database/bau069.

  69. Shugay M, Bagaev D V., Turchaninova M a., Bolotin D a., Britanova O V., Putintseva E V., et al. VDJtools: unifying post-analysis of T cell receptor repertoires. PLoS Comput Biol 2015;11:e1004503. doi:10.1371/journal.pcbi.1004503.

  70. Bagaev DV, Zvyagin IV, Putintseva EV, Izraelson M, Britanova OV, Chudakov DM, et al. VDJviz: a versatile browser for immunogenomics data. BMC Genomics. 2016;17:453. doi:10.1186/s12864-016-2799-7.

    Article  Google Scholar 

  71. Robins HS, Campregher PV, Srivastava SK, Wacher A, Turtle CJ, Kahsai O, et al. Comprehensive assessment of T-cell receptor beta-chain diversity in alphabeta T cells. Blood. 2009;114:4099–107. doi:10.1182/blood-2009-04-217604.

    Article  CAS  Google Scholar 

  72. Klinger M, Pepin F, Wilkins J, Asbury T, Wittkop T, Zheng J, et al. Multiplex Identification of Antigen-Specific T Cell Receptors Using a Combination of Immune Assays and Immune Receptor Sequencing. PLoS One. 2015;10:e0141561. doi:10.1371/journal.pone.0141561.

  73. Nazarov V, Pogorelyy M, Komech E, Zvyagin I, Bolotin D, Shugay M, et al. tcR: an R package for T cell receptor repertoire advanced data analysis. BMC Bioinformatics. 2015;16:175. doi:10.1186/s12859-015-0613-1.

    Article  Google Scholar 

  74. Greiff V, Bhat P, Cook SC, Menzel U, Kang W, Reddy ST. A bioinformatic framework for immune repertoire diversity profiling enables detection of immunological status. Genome Med. 2015;7:49. doi:10.1186/s13073-015-0169-8.

    Article  Google Scholar 

  75. Laydon DJ, Bangham CRM, Asquith B. Estimating T-cell repertoire diversity: limitations of classical estimators and a new approach. Phil Trans R Soc B. 2015;370, 20140291 doi:10.1098/rstb.2014.0291.

  76. Kaplinsky J, Arnaout R, Segers-Nolten IM, Loken MR, Antia R. Robust estimates of overall immune-repertoire diversity from high-throughput measurements on samples. Nat Commun. 2016;7:11881. doi:10.1038/ncomms11881.

    Article  CAS  Google Scholar 

  77. Britanova OV, Putintseva EV, Shugay M, Merzlyak EM, Turchaninova M a, Staroverov DB, et al. Age-related decrease in TCR repertoire diversity measured with deep and normalized sequence profiling. J Immunol. 2014;192:2689–98. doi:10.4049/jimmunol.1302064.

    Article  CAS  Google Scholar 

  78. EFRON B, THISTED R. Estimating the number of unseen species: how many words did Shakespeare know? Biometrika. 1976;63:435–47. doi:10.1093/biomet/63.3.435.

    Google Scholar 

  79. VEGAN DP. A package of R functions for community ecology. J Veg Sci. 2003;14:927–30. doi:10.1111/j.1654-1103.2003.tb02228.x.

    Article  Google Scholar 

  80. Heather JM, Ismail M, Oakes T, Chain B. High-throughput sequencing of the T-cell receptor repertoire: pitfalls and opportunities. Brief Bioinform. 2017;13:bbw138. doi:10.1093/bib/bbw138.

    Article  Google Scholar 

  81. Venturi V, Kedzierska K, Turner SJ, Doherty PC, Davenport MP. Methods for comparing the diversity of samples of the T cell receptor repertoire. J Immunol Methods. 2007;321:182–95. doi:10.1016/j.jim.2007.01.019.

    Article  CAS  Google Scholar 

  82. Mora T, Walczak AM. Quantifying lymphocyte receptor diversity. bioRxiv. 2016. http://www.biorxiv.org/content/early/2016/04/03/046870.

  83. Michal Or-Guil., Carmen Molina-Paris., Rob J. De Boer., Miles Davenport., Ramit Mehr., Veronika Zarnitsyna. Immune System Modeling and Analysis. Frontiers Media SA; 2015. https://books.google.de/books?id=fIxkCgAAQBAJ&dq=immune+diversity+indices&hl=it&source=gbs_navlinks_s.

  84. Redmond D, Poran A, Elemento O, Han A, Glanville J, Hansmann L, et al. Single-cell TCRseq: paired recovery of entire T-cell alpha and beta chain transcripts in T-cell receptors from single-cell RNAseq. Genome Med. 2016;8:80. doi:10.1186/s13073-016-0335-7.

    Article  Google Scholar 

  85. Stubbington MJT, Lönnberg T, Proserpio V, Clare S, Speak AO, Dougan G, et al. T cell fate and clonality inference from single-cell transcriptomes. Nat Methods. 2016;13:329–32. doi:10.1038/nmeth.3800.

    Article  Google Scholar 

  86. Kim S-M, Bhonsle L, Besgen P, Nickel J, Backes A, Held K, et al. Analysis of the paired TCR α- and β-chains of single human T cells. PLoS One. 2012;7:e37338. doi:10.1371/journal.pone.0037338.

    Article  CAS  Google Scholar 

  87. Taylor S, Yasuyama N, Farmer A. A SMARTer approach to profiling the human T-cell receptor repertoire. J Immunol. 2016;196(1 Supplement):209.5.

    Google Scholar 

  88. Williams R, Peisajovich SG, Miller OJ, Magdassi S, Tawfik DS, Griffiths AD. Amplification of complex gene libraries by emulsion PCR. Nat Methods. 2006;3:545–50. doi:10.1038/nmeth896.

    Article  CAS  Google Scholar 

  89. McDaniel JR, DeKosky BJ, Tanno H, Ellington AD, Georgiou G. Ultra-high-throughput sequencing of the immune receptor repertoire from millions of lymphocytes. Nat Protoc. 2016;11:429–42. doi:10.1038/nprot.2016.024.

    Article  CAS  Google Scholar 

  90. Best K, Oakes T, Heather JM, Shawe-Taylor J, Chain B, Polz MF, et al. Computational analysis of stochastic heterogeneity in PCR amplification efficiency revealed by single molecule barcoding. Sci Rep. 2015;5:14629. doi:10.1038/srep14629.

    Article  CAS  Google Scholar 

  91. Turchaninova MA, Britanova OV, Bolotin DA, Shugay M, Putintseva EV, Staroverov DB, et al. Pairing of T-cell receptor chains via emulsion PCR. Eur J Immunol. 2013;43:2507–15. doi:10.1002/eji.201343453.

    Article  CAS  Google Scholar 

  92. Chattopadhyay PK, Gierahn TM, Roederer M, Love JC. Single-cell technologies for monitoring immune systems. Nat Immunol. 2014;15:128–35. doi:10.1038/ni.2796.

    Article  CAS  Google Scholar 

  93. Proserpio V, Mahata B. Single-cell technologies to study the immune system. Immunology. 2016;147:133–40. doi:10.1111/imm.12553.

    Article  CAS  Google Scholar 

  94. Smith T, Heger A, Sudbery I. UMI-tools: Modeling sequencing errors in unique molecular identifiers to improve quantification accuracy. Genome Res. 2017;27:491–9.

    Article  CAS  Google Scholar 

  95. Liaskou E, Henriksen EKK, Holm K, Kaveh F, Hamm D, Fear J, et al. High-throughput T-cell receptor sequencing across chronic liver diseases reveals distinct disease-associated repertoires. Hepatology. 2016;63:1608–19. doi:10.1002/hep.28116.

    Article  CAS  Google Scholar 

  96. Caporaso JG, Lauber CL, Walters W A, Berg-Lyons D, Lozupone C A, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci U S A. 2011;108 Suppl:4516–22.

Download references

Acknowledgements

Andre Franke, Marie Dowds and Elisa Rosati receive infrastructure support from the DFG Cluster of Excellence No. 306 „Inflammation at Interfaces” and the DFG Research Training Group (RTG) No. 1743 „Genes, Environment and Inflammation”.

Funding

This work was founded by the Norwegian PSC Research Center (NoPSC) and the DFG Cluster of Excellence No. 306 “Inflammation at Interfaces”.

Availability of data and materials

The datasets analysed during the current study are available in the European Nucleotide Archive (ENA), with study accession number PRJEB21352 (http://www.ebi.ac.uk/ena/data/view/PRJEB21352).

Author information

Authors and Affiliations

Authors

Contributions

ER performed the RNA-based benchmarking experiments, data analyses and drafted the first version of the manuscript. EL and THK provided samples for the benchmarking and data. EKH helped to design manuscript Figs. MD, THK and AF designed and supervised the work. All authors critically read the manuscript and approved the final version.

Corresponding author

Correspondence to Andre Franke.

Ethics declarations

Ethics approval and consent to participate

All samples were obtained with Local Research and Ethics approval and informed patient consent (LREC references: 2003/242 South Birmingham REC, renewed 2012; and 06/Q2702/61 Black Country REC).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional files

Additional file 1:

Supplementary material. Supplementary figures mentioned in the main manuscript and their titles and legends. Figure S1. Venn diagram showing the overlap between the most abundant TCR sequences detected by each method. The threshold was defined by the method which detected fewest TCR species (UMI-corrected 5’RACE-based PCR). The diagrams show the number and relative frequency of TCR sequences detected by only one up to all four methods. Sequences which were found only by one method may still be detected by other methods, but may not appear in the most abundant species and are thus not represented here. Data are shown for two patients and for both α and β chains. Figure S2. Venn diagram showing the overlap between the top 300 most abundant TCR sequences detected by each method. For every technique (BGI, 5’RACE PCR, 5’RACE PCR UMI-corrected, iRepertoire kit) the 300 most abundant clonotypes were considered. The diagrams show how many sequences were found to be present also in the most abundant 300 clonotypes of other methods. Sequences which are shown to be found only by one method may still be detected by other methods, but don’t appear in the 300 most abundant species of these. Data are shown for two patients and for both α and β chains. Figure S3. Variable (V) and Joining (J) genes usage among methods and replicates. The heat map shows the gene usage proportion inside each sample for the Variable and Joining region alleles listed at the bottom on the figure. Each of the three parts of the heat map is representative of one of the methods. The samples are described by the patients “P1/2”, the method “iRepertoire/BGI/5’RACE/5’RACE + UMI” and the replicate “1 or 2” or “0” in case of BGI which doesn’t have replicates. The figure was generated using the “ggplot2” R package. a) Beta chain V genes. b) Alpha chain J genes. c) Beta chain J genes. Figure S4. Dissimilarity dendrogram for beta chain. The distance between patients, methods, and replicates, calculated using the Morisita index, is represented by the dendrograms. The samples are described by the patients “P1/2”, the method “iRepertoire/BGI/5’RACE/5’RACE + UMI” and the replicate “1 or 2”. iRepertoire replicate 2 for beta chain did not satisfy the data analysis quality criteria and was therefore excluded. The figure was generated using the “ggplot2” and the “ggdendro” R packages. (DOCX 3501 kb)

Additional file 2:

UMI filtering script (python). The script used to filter the unique molecular identifiers. (PY 4 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rosati, E., Dowds, C.M., Liaskou, E. et al. Overview of methodologies for T-cell receptor repertoire analysis. BMC Biotechnol 17, 61 (2017). https://doi.org/10.1186/s12896-017-0379-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12896-017-0379-9

Keywords