VDJviz: a versatile browser for immunogenomics data
© The Author(s). 2016
Received: 1 December 2015
Accepted: 25 May 2016
Published: 13 June 2016
The repertoire of T- and B-cell receptor sequences encodes the antigen specificity of adaptive immunity system, determines its present state and guides its ability to mount effective response against encountered antigens in future. High throughput sequencing of immune repertoires (Rep-Seq) is a promising technique that allows to profile millions of antigen receptors of an individual in a single experiment. While a substantial number of tools for mapping and assembling Rep-Seq data were published recently, the field still lacks an intuitive and flexible tool that can be used by researchers with little or no computational background for in-depth analysis of immune repertoire profiles.
Here we report VDJviz, a web tool that can be used to browse, analyze and perform quality control of Rep-Seq results generated by various pre-processing software. On a set of real data examples we show that VDJviz can be used to explore key repertoire characteristics such as spectratype, repertoire clonality, V-(D)-J recombination patterns and to identify shared clonotypes. We also demonstrate the utility of VDJviz in detection of critical Rep-Seq biases such as artificial repertoire diversity and cross-sample contamination.
VDJviz is a versatile and lightweight tool that can be easily employed by biologists, immunologists and immunogeneticists for routine analysis and quality control of Rep-Seq data. The software is freely available for non-commercial purposes, and can be downloaded from: https://github.com/antigenomics/vdjviz.
KeywordsHigh-throughput sequencing Immunology Browser T-cell B-cell Repertoire sequencing
A diverse repertoire of T- and B-cell antigen receptors is a critical component of host defense system in vertebrates called adaptive immunity which ensures readiness and ability to detect and mount an effective response against the great variety of encountered pathogens. T- and B-cell receptor repertoire is formed by genomic rearrangement of Variable (V), Diversity (D) and Joining (J) segment loci in a process called V-(D)-J recombination . Each of the resulting segment junctions carry complementarity determining region 3 (CDR3) that plays a key role in antigen recognition and largely defines the specificity of T-cells and immunoglobulins . Recent advances in molecular methods and high-throughput sequencing allow to profile antigen receptor repertoires using a technique called Rep-Seq . Raw receptor sequences produced by Rep-Seq can be processed by one of the existing bioinformatics software tools (http://omictools.com/rep-seq-c424-p1.html) to map V-(D)-J junctions and extract CDR3 regions, thus forming a set of clonotypes - unique combinations of V, D and J segments, and CDR3 sequence. Those mappings are then assembled to estimate individual clonotype frequencies that reflect clonal expansions caused by antigen recognition, peripheral selection  and convergent V-(D)-J recombination processes . Resulting datasets are inherently complex due to extremely high diversity of T- and B-cell receptor sequences and a plethora of physiological factors that shape the repertoire structure .
A parser that supports output from 6 commonly used Rep-Seq processing software: MiTCR, MIGEC, MiXCR, IgBlast, IMGT HighVQuest and ImmunoSEQ platform; as well as an internal concise tab-delimited format.
An intuitive clonotype table viewer with V-(D)-J markup that can be used to navigate through the entire sample and perform complex searches.
Comprehensive single sample analysis modules calculating basic repertoire statistics and providing interactive spectratype, V/J segment usage and clonality plots.
Extended multi-sample analysis that includes clonotype tracking and sample intersection with a flexible set of clonotype matching rules, repertoire diversity comparison using rarefaction and simple side-by-side comparison of single sample analysis results.
Export of analysis results and dataset sharing.
VDJviz is a web based GUI application that uses VDJtools API (https://github.com/mikessh/vdjtools) as a back-end. The software utilizes Play framework (https://www.playframework.com/) for running the server instance and state-of-art web graphics libraries such as D3js (http://d3js.org/) for visualization. The reason for choosing Play framework is its stability and ease of deployment, while D3js allows us to create complex interactive plots. The browser is lightweight and uses around 4GM RAM to host several users analyzing 25 samples of up to 10,000 clonotypes, which is the upload limit for the demo version available online. This limit can be removed for local installations to allow browsing large samples using better hardware. In most cases users can also down-sample clonotype abundance tables (http://vdjtools-doc.readthedocs.org/en/latest/preprocess.html#downsample) to view large samples with commodity hardware.
Results and discussion
To the best of our knowledge, in contrast to the rich set browsers available in the field of genomics (e.g. Refs. [16–18]), the only published software that falls in immune repertoire browser category so far is IMEX . Existing unpublished solutions for immune repertoire browsing include VDJserver (https://vdjserver.org/), Vidjil browser (http://www.vidjil.org/#browser) and ImmunoSEQ Analyzer (https://clients.adaptivebiotech.com/). In this section we will first compare the functionality of VDJviz with aforementioned web tools and demonstrate VDJviz features on the set of relevant examples further in the text.
IMEX is a closed-source GUI software that allows computing basic repertoire statistics, analyzing V-D-J segment usage, performing diversity estimation and provides some options for comparing clonotype tables. The software is implemented using .NET technology and can natively run on Windows. Running it on Unix-based systems requires setting up the Mono Framework (http://www.mono-project.com/). There are several general limitations of IMEX compared to VDJviz. First, IMEX limits its analysis to datasets produced by IMGT High-V/Quest software while VDJviz allows both IMGT High-V/Quest input and input generated by other software tools. IMGT High-V/Quest is frequently used by immunologists, however the current upload batch size of 0.5 mln reads and variable submission/processing times makes it unfeasible for analysis of large datasets containing tens of millions of reads. Next, currently IMEX is limited to TRB and IGH loci while no such limitation exists in VDJviz. However the most important issue with IMEX is the way it estimates one of the key immune repertoire parameters, repertoire diversity . IMEX fits the a × (1-exp(−b × n)) + k × n function, where n is sampling depth, a is real number of clonotypes and k is the error rate, using an optimization algorithm to the rarefaction curve obtained by random re-sampling. This empirical model can produce spurious results in some common settings and cannot reliably distinguish rare clonotypes and errors. For example, let us consider a highly diverse and uniform repertoire (say, naive T-cells) and note that corresponding Rep-Seq data can have negligible error rate if produced using high-fidelity protocols [12, 21]. The rarefaction curve in error-free setting is a linear function of sample size . On the other hand, the optimal parameters for model used in IMEX can be selected as k = 1 and a = 0/b = any or b = 0/a = any in this setting, thus either rendering all clonotypes as erroneous or providing an arbitrary number of clonotypes in a sample that depends on the seed of the random number generator used by the optimization algorithm. VDJviz, on the other hand implements a robust and commonly used rarefaction algorithm  leaving the choice of error correction strategy up to the user.
Vidjill browser is an extension of recently published Vidjil Rep-Seq processing software . The major difference between Vidjil browser and VDJviz lies in the repertoire browsing implementation and repertoire analysis features. Vidjil browser operates with V-D-J signatures of clonotype clusters and implements a graphical clonotype tracking interface with an aim to facilitate clonotype tracking for MRD detection and monitoring. VDJviz, on the other hand, lists individual clonotypes in tabular format and all the relevant information such as V,D and J segments and the CDR3 region sequence, which allows to directly browse the clonal composition of sample and perform clonotype table searches using pattern-matching and filters. Notably, VDJviz implements some of the commonly used analysis modes such as diversity estimation and spectratyping that are not present in Vidjil browser. VDJviz also implements basic clonotype tracking functionality in its cross-sample intersection and clonotype search modules. VDJviz doesn’t limit clonotype tracking to samples coming from the same donor, allowing to match clonotypes based on CDR3 amino acid sequences and therefore allows exploring clonotypes shared by several different donors.
VDJserver software, being in beta version, incorporates V-D-J mapping engine and requires to upload raw sequencing data, which can be both considered as a benefit and a limitation comparing to VDJviz that accepts processed data in multiple formats. While doing data processing on server side facilitates analysis for data produced using common library preparation protocols, it is unfeasible to implement a general algorithm that covers all possible customizations of those protocols and complex cases such as multiplexing and unique molecular identifier tagging . The output provided by VDJtools includes segment usage chart and V-D-J mapping statistics, while clonotype tables are only available as a downloadable plain text file, which is far less than the functionality provided by VDJviz, Vidjil browser and IMEX.
ImmunoSEQ analyzer is a commercial software and supports only customer data produced by corresponding commercial assay. ImmunoSEQ has a rich feature set, some of which are not present in VDJviz, namely a variety scatterplots for sample comparison, immunoglobulin somatic hypermutation and edit distance analysis. VDJviz, on the other hand, offers more options for diversity estimation including rarefaction analysis and clonality plot, clonotype-level detalization for sample intersection and a powerful clonotype search engine. Clonotype search algorithms of VDJviz are also more flexible: various filters such as segment filter can be used in combination, user can search for CDR3 sequence patterns and several clonotype matching modes are supported, for example amino acid-not-nucleotide matching that can be used to filter cross-sample contaminations.
Below we present six example cases that demonstrate the usability of VDJviz for common immune repertoire analysis tasks, in-depth browsing of repertoire clonal composition and detection of Rep-Seq artifacts. Data for reproducing all the examples presented here is available in the “examples” folder of VDJviz source code repository, all figures in this paper are screenshots of VDJviz browser interface.
Example 1: spectratyping
Spectratype can also be used to spot out-of-frame clonotypes and a thorough look at Fig. 1 reveals the abundance of out-of-frame clonotypes in 6 yo sample (small bars corresponding to CDR3 lengths that are not a multiple of 3). Summary report provided by VDJviz shows that the total abundance of out-of-frame clonotypes is ~2 times more for the young donor compared to aged one (P < 0.0001, Fisher’s exact test). Out-of-frame TCR sequences are extremely useful for studying V-(D)-J recombination mechanics [27–31] as they are not subject to thymic selection. However, they are relatively rare in mRNA-based samples due to nonsense-mediated mRNA decay. The result shown on Fig. 1 suggests that a deeper sampling is required to detect a sufficient number of out-of-frame clonotypes for repertoires having a high fraction of expanded clonotypes.
Example 2: variable and joining segment usage
Example 3: clonality analysis
Example 4: rarefaction and error correction
First, we tested quality-based filtering (without using UMI information) by removing all clonotypes that have at least one low-quality base (less than Phred quality 20 and 35) in their CDR3 sequence. Such filtering has a relatively small effect on the observed diversity which is more than 10 times higher than the value expected from UMI-corrected results. Notably, using only clonotypes that were found in both samples results in observed diversity that is still ~5 times larger than the expected value, confirming the previous observations that errors that result in artificial diversity are highly reproducible . This indicates that using replicate-based error correction to investigate repertoire diversity  is a strategy that should be applied with a great caution. Using frequency-based error correction for quality-filtered sample, namely requiring more than 1:20 abundance ratio difference for merging clonotypes that differ by a single mismatch, showed the best result, yet the observed diversity was still ~60 % higher than the expected value obtained using UMI-correction.
The accuracy and pitfalls of quality- and frequency-based error filtering strategies were previously characterized using a synthetic dataset . As for the comprehensive characterization of the accuracy of UMI-based techniques that are considered as gold-standard in present example, the reader should refer to recently published studies [12, 38].
Example 5: errors and contamination in repertoire sequencing data
Next, we have browsed samples from recently published minimal residual disease (MRD) study  to address the issue of cross-sample contamination. In lymphomas, MRD can be monitored by tracking the malignant clonotype sequence in post-treatment immune repertoire. Cross-sample contamination, however, can be a serious issue in this case: contamination with the malignant clonotype which is usually highly abundant can lead to false positive MRD detection. We have checked for the cancer clonotype sequence reported for patient PT-2 in post-treatment samples of other patients (Fig. 5b). Notably, repertoires of 6 out of 42 patients appear to contain exactly the same sequence. This can be hardly explained by coincidence, as no other nucleotide variants were found for the dominant clonotype’s amino acid sequence ruling out convergent recombination. High number of added N-nucleotides in V-D-J junction also supports the fact that CDR3 nucleotide sequence matching in 6 samples simply by chance is highly improbable (P < 10−62). Notably, the most abundant contamination is present at the level of 24 reads per 100,000 in this example case. Thus, while the method is extremely sensitive and can detect MRD at a level of 1 read per 100,000, such contaminations can severely dampen method’s precision.
Example 6: public clonotypes
While the examples demonstrated here mostly deal with TCR beta sequences, VDJviz can also handle TCR alpha (see Additional file 1: Figure S1), gamma and delta sequences, as well as immunoglobulin sequences (see Additional file 2: Figure S2), albeit with no support for hypermutations in CDR1,2 and framework regions. Overall the examples presented here demonstrate that the analysis modes provided by VDJviz are highly informative and can be used both for explorative analysis and for quality control. The latter is crucial as a multitude of biases can arise due to complexity of Rep-Seq data. While those biases can be dealt with using corresponding techniques or removed manually, their extent should be routinely checked every time an analysis of Rep-Seq data is performed.
While VDJviz web tool can be extended in many ways by adding new analysis types, the most important challenge is to implement intuitive interface for visualizing somatic hypermutations in B-cell repertoires  and novel paired-chain Rep-Seq data [45–47].
As we have demonstrated, VDJviz allows to have a grasp of immune repertoire structure for samples of interest in several clicks and can be easily used by immunologists and biologists with little computational knowledge. VDJviz is not limited to a single library preparation protocol or Rep-Seq processing software including highly popular IMGT HighVQuest  and ImmunoSEQ platforms (http://www.adaptivebiotech.com/immunoseq). VDJviz allows great flexibility and can be easily installed as a local server, therefore we believe that in perspective it can become a handy tool-of-choice for immunologists routinely working with immune repertoire data.
Availability and requirements
• Both standalone and online demo VDJviz versions can be found at https://github.com/antigenomics/vdjviz.
• Operating system(s): platform independent.
• Other requirements: Java 1.8, Play Framework.
• License: free for non-profit and academic use.
This work was supported by the Russian Science Foundation project №14-14-00533. Shugay M and Putintseva EV are supported by individual fellowship mol-a-dk RFBR grants 16-34-60179 and 16-34-60178. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Library preparations was carried out in part using equipment provided by the Shemyakin–Ovchinnikov Institute of Bioorganic Chemistry Core Facility.
DVB and MS developed software. IVZ, OVB, MI, and EVP prepared the cDNA libraries and performed software testing. MS worked on the manuscript. DMC worked on the manuscript and supervised the work. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis 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.
- Janeway CA. Immunobiology: the immune system in health and disease. 8th ed. New York: Garland Science; 2012.Google Scholar
- Xu JL, Davis MM. Diversity in the CDR3 region of V(H) is sufficient for most antibody specificities. Immunity. 2000;13(1):37–45.View ArticlePubMedGoogle Scholar
- Benichou J, Ben-Hamo R, Louzoun Y, Efroni S. Rep-Seq: uncovering the immunological repertoire through next-generation sequencing. Immunology. 2012;135(3):183–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Rocha B, von Boehmer H. Peripheral selection of the T cell repertoire. Science. 1991;251(4998):1225–8.View ArticlePubMedGoogle Scholar
- Quigley MF, Greenaway HY, Venturi V, Lindsay R, Quinn KM, Seder RA, et al. Convergent recombination shapes the clonotypic landscape of the naive T-cell repertoire. Proc Natl Acad Sci U S A. 2010;107(45):19414–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Nikolich-Zugich J, Slifka MK, Messaoudi I. The many important facets of T-cell repertoire diversity. Nat Rev Immunol. 2004;4(2):123–32.View ArticlePubMedGoogle Scholar
- Georgiou G, Ippolito GC, Beausang J, Busse CE, Wardemann H, Quake SR. The promise and challenge of high-throughput sequencing of the antibody repertoire. Nat Biotechnol. 2014;32(2):158–68.View ArticlePubMedPubMed CentralGoogle Scholar
- Brochet X, Lefranc MP, Giudicelli V. IMGT/V-QUEST: the highly customized and integrated system for IG and TR standardized V-J and V-D-J sequence analysis. Nucleic Acids Res. 2008;36(Web Server issue):W503–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Bolotin DA, Shugay M, Mamedov IZ, Putintseva EV, Turchaninova MA, Zvyagin IV, et al. MiTCR: software for T-cell receptor sequencing data analysis. Nat Methods. 2013;10(9):813–4.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- Ye J, Ma N, Madden TL, Ostell JM. IgBLAST: an immunoglobulin variable domain sequence analysis tool. Nucleic Acids Res. 2013;41(Web Server issue):W34–40.View ArticlePubMedPubMed CentralGoogle Scholar
- 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(6):653–5.View ArticlePubMedGoogle Scholar
- Bolotin DA, Poslavsky S, Mitrophanov I, Shugay M, Mamedov IZ, Putintseva EV, et al. MiXCR: software for comprehensive adaptive immunity profiling. Nat Methods. 2015;12(5):380–1.View ArticlePubMedGoogle Scholar
- 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(4):e31.View ArticlePubMedGoogle Scholar
- Kuchenbecker L, Nienen M, Hecht J, Neumann AU, Babel N, Reinert K, et al. IMSEQ-a fast and error aware approach to immunogenetic sequence analysis. Bioinformatics. 2015;31(18):2963–71.View ArticlePubMedGoogle Scholar
- Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29(1):24–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Nicol JW, Helt GA, Blanchard Jr SG, Raja A, Loraine AE. The Integrated Genome Browser: free software for distribution and exploration of genome-scale datasets. Bioinformatics. 2009;25(20):2730–1.View ArticlePubMedPubMed CentralGoogle Scholar
- Meyer F, Paarmann D, D’Souza M, Olson R, Glass EM, Kubal M, et al. The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics. 2008;9:386.View ArticlePubMedPubMed CentralGoogle Scholar
- Schaller S, Weinberger J, Jimenez-Heredia R, Danzer M, Oberbauer R, Gabriel C, et al. ImmunExplorer (IMEX): a software framework for diversity and clonality analyses of immunoglobulins and T cell receptors on the basis of IMGT/HighV-QUEST preprocessed NGS data. BMC Bioinformatics. 2015;16:252.View ArticlePubMedPubMed CentralGoogle Scholar
- 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:413.View ArticlePubMedPubMed CentralGoogle Scholar
- Vander Heiden JA, Yaari G, Uduman M, Stern JN, O’Connor KC, Hafler DA, et al. pRESTO: a toolkit for processing high-throughput sequencing raw reads of lymphocyte receptor repertoires. Bioinformatics. 2014;30(13):1930–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Colwell RK, Chao A, Gotelli NJ, Lin S, Mao CX, Chazdon RL, et al. Models and estimators linking individual-based and sample-based rarefaction, extrapolation and comparison of assemblages. J Plant Ecol. 2012;5(1):3–21.View ArticleGoogle Scholar
- Giraud M, Salson M, Duez M, Villenet C, Quief S, Caillault A, et al. Fast multiclonal clusterization of V(D)J recombinations from high-throughput sequencing. BMC Genomics. 2014;15:409.View ArticlePubMedPubMed CentralGoogle Scholar
- Britanova OV, Putintseva EV, Shugay M, Merzlyak EM, Turchaninova MA, Staroverov DB, et al. Age-related decrease in TCR repertoire diversity measured with deep and normalized sequence profiling. J Immunol. 2014;192(6):2689–98.View ArticlePubMedGoogle Scholar
- Kivioja T, Vaharautio A, Karlsson K, Bonke M, Enge M, Linnarsson S, et al. Counting absolute numbers of molecules using unique molecular identifiers. Nat Methods. 2012;9(1):72–4.View ArticleGoogle Scholar
- Yager EJ, Ahmed M, Lanzer K, Randall TD, Woodland DL, Blackman MA. Age-associated decline in T cell repertoire diversity leads to holes in the repertoire and impaired immunity to influenza virus. J Exp Med. 2008;205(3):711–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Murugan A, Mora T, Walczak AM, Callan Jr CG. Statistical inference of the generation probability of T-cell receptors from sequence repertoires. Proc Natl Acad Sci U S A. 2012;109(40):16161–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Putintseva EV, Britanova OV, Staroverov DB, Merzlyak EM, Turchaninova MA, Shugay M, et al. Mother and child t cell receptor repertoires: deep profiling study. Front Immunol. 2013;4:463.View ArticlePubMedPubMed CentralGoogle Scholar
- 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;24(10):1603–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Zvyagin IV, Pogorelyy MV, Ivanova ME, Komech EA, Shugay M, Bolotin DA, et al. Distinctive properties of identical twins’ TCR repertoires revealed by high-throughput sequencing. Proc Natl Acad Sci U S A. 2014;111(16):5980–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Rubelt F, Bolen CR, McGuire HM, Heiden JA, Gadala-Maria D, Levin M, et al. Individual heritable differences result in unique cell lymphocyte receptor repertoires of naive and antigen-experienced cells. Nat Commun. 2016;7:11112.View ArticlePubMedGoogle Scholar
- Britanova OV, Bochkova AG, Staroverov DB, Fedorenko DA, Bolotin DA, Mamedov IZ, et al. First autologous hematopoietic SCT for ankylosing spondylitis: a case report and clues to understanding the therapy. Bone Marrow Transplant. 2012;47(11):1479–81.View ArticlePubMedGoogle Scholar
- Mamedov IZ, Britanova OV, Bolotin DA, Chkalina AV, Staroverov DB, Zvyagin IV, et al. Quantitative tracking of T cell clones after haematopoietic stem cell transplantation. EMBO Mol Med. 2011;3(4):201–7.View ArticlePubMedPubMed CentralGoogle Scholar
- 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(3):1168–72.View ArticlePubMedPubMed CentralGoogle Scholar
- Chao A, Chiu C, Hsieh TC, Davis T, Nipperess DA, Faith DP. Rarefaction and extrapolation of phylogenetic diversity. Methods Ecol Evol. 2015;6(4):380–8.View ArticleGoogle Scholar
- 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(12):6155–63.View ArticlePubMedGoogle Scholar
- Qi Q, Liu Y, Cheng Y, Glanville J, Zhang D, Lee JY, et al. Diversity and clonal selection in the human T-cell repertoire. Proc Natl Acad Sci U S A. 2014;111(36):13139–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Khan TA, Friedensohn S, de Vries AR, Straszewski J, Ruscheweyh HJ, Reddy ST. Accurate and predictive antibody repertoire profiling by molecular amplification fingerprinting. Science Advances. 2016;2(3), e1501371.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu D, Sherwood A, Fromm JR, Winter SS, Dunsmore KP, Loh ML, et al. High-throughput sequencing detects minimal residual disease in acute T lymphoblastic leukemia. Sci Transl Med. 2012;4(134):134ra163.View ArticleGoogle Scholar
- Venturi V, Price DA, Douek DC, Davenport MP. The molecular basis for public T-cell responses? Nat Rev Immunol. 2008;8(3):231–8.View ArticlePubMedGoogle Scholar
- 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(10):1817–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang C, Sanders CM, Yang Q, Schroeder Jr HW, Wang E, Babrzadeh F, et al. High throughput sequencing reveals a complex pattern of dynamic interrelationships among human T cell subsets. Proc Natl Acad Sci U S A. 2010;107(4):1518–23.View ArticlePubMedPubMed CentralGoogle Scholar
- 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. Scientific Reports. 2015;5:13664.View ArticlePubMedPubMed CentralGoogle Scholar
- Shugay M, Bolotin DA, Putintseva EV, Pogorelyy MV, Mamedov IZ, Chudakov DM. Huge overlap of individual TCR beta repertoires. Front Immunol. 2013;4:466.View ArticlePubMedPubMed CentralGoogle Scholar
- 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(9):2507–15.View ArticlePubMedGoogle Scholar
- DeKosky BJ, Ippolito GC, Deschner RP, Lavinder JJ, Wine Y, Rawlings BM, et al. High-throughput sequencing of the paired human immunoglobulin heavy and light chain repertoire. Nat Biotechnol. 2013;31(2):166–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Howie B, Sherwood AM, Berkebile AD, Berka J, Emerson RO, Williamson DW, et al. High-throughput pairing of T cell receptor alpha and beta sequences. Sci Transl Med. 2015;7(301):301ra131.View ArticlePubMedGoogle Scholar
- Feng Y, van der Veeken J, Shugay M, Putintseva EV, Osmanbeyoglu HU, Dikiy S, et al. A mechanism for expansion of regulatory T-cell repertoire and its role in self-tolerance. Nature. 2015;528(7580):132–6.PubMedPubMed CentralGoogle Scholar