VDJdb in 2019: database extension, new analysis infrastructure and a T-cell receptor motif compendium (2024)

  • Journal List
  • Nucleic Acids Res
  • v.48(D1); 2020 Jan 8
  • PMC6943061

As a library, NLM provides access to scientific literature. Inclusion in an NLM database does not imply endorsem*nt of, or agreement with, the contents by NLM or the National Institutes of Health.
Learn more: PMC Disclaimer | PMC Copyright Notice

VDJdb in 2019: database extension, new analysis infrastructure and a T-cell receptor motif compendium (1)

Link to Publisher's site

Nucleic Acids Res. 2020 Jan 8; 48(D1): D1057–D1062.

Published online 2019 Oct 7. doi:10.1093/nar/gkz874

PMCID: PMC6943061

PMID: 31588507

Dmitry V Bagaev,1,2 Renske M A Vroomans,3,4 Jerome Samir,5,6 Ulrik Stervbo,7,8 Cristina Rius,9 Garry Dolton,9 Alexander Greenshields-Watson,9 Meriem Attaf,9 Evgeny S Egorov,2 Ivan V Zvyagin,1,2 Nina Babel,7,8 David K Cole,9,10 Andrew J Godkin,9 Andrew K Sewell,9 Can Kesmir,11 Dmitriy M Chudakov,1,2,12 Fabio Luciani,5,6 and Mikhail Shugay1,2,12

Author information Article notes Copyright and License information PMC Disclaimer

Associated Data

Data Availability Statement

Abstract

Here, we report an update of the VDJdb database with a substantial increase in the number of T-cell receptor (TCR) sequences and their cognate antigens. The update further provides a new database infrastructure featuring two additional analysis modes that facilitate database querying and real-world data analysis. The increased yield of TCR specificity identification methods and the overall increase in the number of studies in the field has allowed us to expand the database more than 5-fold. Furthermore, several new analysis methods are included. For example, batch annotation of TCR repertoire sequencing samples allows for annotating large datasets on-line. Using recently developed bioinformatic methods for TCR motif mining, we have built a reduced set of high-quality TCR motifs that can be used for both training TCR specificity predictors and matching against TCRs of interest. These additions enhance the versatility of the VDJdb in the task of exploring T-cell antigen specificities. The database is available at https://vdjdb.cdr3.net.

INTRODUCTION

Knowing the exact antigen specificity of a given T-cell is key to solving numerous problems of both basic and applied immunology research: from discovering the specificity profile of TCR repertoire sequencing samples (1,2), to finding associations between autoimmunity and foreign mimics of self-antigens (3), and designing of personalized tumor immunotherapies (4). The field of molecular methods designed for studying antigen-specific T-cells is developing at a high pace: novel methodologies based on single-cell T-cell sequencing allow simultaneous detection of TCR sequence, T-cell phenotype and a vast array of antigen specificities (5). Resulting data, however, still need to be properly quality-controlled, and organized in the form of a database that is both comprehensive and easy to query.

After the first version of VDJdb (6) and a pathology-associated TCR database (McPAS-TCR (7)) were published, a commonly used iEDB database that describes immunogenic antigens was also modified to add metadata related to TCR and B-cell receptor sequences (8), highlighting the overall demand for such data in the field. A number of methods for TCR specificity prediction were also reported recently, many of which rely on VDJdb data for training and validating classifiers (9–12). The latter demonstrates the overall potential of the VDJdb database for developing better bioinformatic methods for TCR sequence analysis.

In this 2019 update,we focused on both accumulating the large amount of data generated by recent studies and providing an interface facilitating web-based analysis of adaptive immune receptor repertoire sequencing (AIRR-Seq, (13)) datasets. Given the large amount of data currently stored in VDJdb, we provided a reduced dataset of high-quality motifs that facilitates identification of TCR residues critical for recognition of certain antigens.

MATERIALS AND METHODS

Data acquisition and processing

All data acquired from published studies was manually parsed into VDJdb format according to VDJdb guidelines (https://github.com/antigenomics/vdjdb-db/blob/master/README.md) and quality-controlled both manually and using previously reported automated scripts (6). Summary statistics for VDJdb records were computed using an R notebook provided at https://github.com/antigenomics/vdjdb-db/blob/master/summary/vdjdb_summary.Rmd. Datasets from 10X genomics were downloaded from https://support.10xgenomics.com/single-cell-vdj/datasets (‘Application Note - A New Way of Exploring Immunity’ section, datasets ‘CD8+ T cells of Healthy Donor’ 1–4, available under the Creative Commons Attribution license) and processed using in-house scripts (available at https://bitbucket.org/kirbyvisp/10x-tcr/src/master/) that perform stringent filtering on the CDR3 length and composition, and tetramer read counts, yielding over 20 000 unique antigen-specific receptors. We’ve also performed several rounds of manual proofreading for the database and fixed a number of typos, mostly related to ambiguous segment naming (e.g. cases when V segments named according to Arden nomenclature were imported as IMGT segment names).

Updated VDJdb web browser implementation details

Since 2017, we upgraded VDJdb web server to run on the latest Play framework (v2.7.2, https://www.playframework.com) with Akka HTTP server to improve the overall performance. We have also fully re-implemented frontend using Angular (https://angular.io) to provide a faster and more responsive interface. Importantly, we have implemented a fully documented REST API that can be used to query the database and can be found at https://vdjdb-web.readthedocs.io/en/latest/api.html. We also facilitated local VDJdb web server installation by providing a Docker image available at https://cloud.docker.com/u/bvdmitri/repository/docker/bvdmitri/vdjdb-web.

TCR motif database

We have used the TCRNET implementation (14) in VDJtools (15) to identify TCR nodes in VDJdb TCR similarity network that have more neighbors than expected by chance, allowing for a single amino acid substitution in the CDR3 region. Only epitopes assigned to at least 30 distinct TCR amino acid sequences were considered. Selected nodes and their first neighbors were left in the TCR similarity network and sets of hom*ologous TCR sequences (motifs) were defined for each epitope as connected components of the resulting graph. Position weight matrices (PWMs) for CDR3 amino acid sequences of inferred motifs were constructed using connected components of the graph. PWM normalization was performed by using the probability in a control set as the information measure, where control set is a set TCR sequences having the same V/J genes and CDR3 length coming from a pool of healthy donor samples. Details of this procedure are summarized in an R markdown notebook available at https://github.com/antigenomics/vdjdb-motifs.

RESULTS

Timeline of data accumulation and perspectives

VDJdb database is substantially expanded compared to the previous report (6): since the establishment of the database a total of 155 published studies were processed and added resulting in 61 049 TCR specificity records in 2019 compared to only 5491 in 2017. We used the publication dates of papers added to VDJdb to calculate the rate of accumulation of TCR specificity knowledge (Figure ​(Figure1).1). The figure clearly shows that the number of records grows slowly prior to 2017, while a rapid growth in the number of records occurred in the last two years. The latter can be explained by the establishment of AIRR-Seq techniques producing very large amounts of TCR sequences (13) as a method of choice for performing the readout of TCR specificity assays.

Open in a separate window

Figure 1.

Summary statistics of VDJdb records according to July 2019 database release. Plots show the cumulative number of unique TCR sequences (up to unique V/J gene and CDR3 amino acid), antigens, MHC alleles and publications added to the database arranged by the publication date of corresponding papers. Colored lines represent records that only have TCR alpha (TRA, blue) or beta (TRB, green) chain or both chains (‘paired’ records, red). Note that as several VDJdb records can link a single TCR sequence to different metadata (e.g. another study and donor, or distinct epitope in case of cross-reactivity), the total number of unique TCR sequences (n = 42 211) is less than the total number of VDJdb records (n = 61 049).

Notably, while the number of unique epitopes with known specific TCRs increased, it is still dwarfed by the number of unique TCRs in the database, as there are currently only 212 epitopes in VDJdb. This highlights one of the most important factors currently limiting our ability to study and predict TCR specificity for a diverse range of antigens. We hope that in future this can be resolved using technologies allowing simultaneous testing for several antigen specificities (5).

Prior to 2019, there was a substantial lack of paired TCR records as most high-throughput studies focused on TCR beta chain only. However, with the advent of the droplet based single-cell sequencing techniques and the subsequent AIRR-seq focused improvements, the number of paired TCR alpha and beta records increased a lot. For example, one of the most recent studies added to VDJdb, a paired-chain dataset from 10X Genomics compendium (see Materials and Methods section), features 40 unique epitopes, more than any other high-throughput dataset produced so far.

Online batch analysis of AIRR-Seq samples

Due to a growing number of immunological studies that choose AIRR-Seq technology to survey adaptive immunity, one of the main tasks for the current VDJdb update was to provide the ability to query a large TCR sequence set (e.g. whole blood or tissue specific T cell repertoire) against the entire database. We extended the functionality of the VDJdb browser and included batch upload and query options for AIRR-Seq datasets (Annotation tab of the web interface). Batch query supports several commonly used data formats (IMGT/HighV-QUEST, MiXCR, ImmunoSEQ, etc), borrowing functionality from the VDJviz browser (16).

Upon uploading a sample, TCR sequences are aligned against VDJdb records and users are provided with a list of VDJdb hits containing all necessary information regarding alignment and metadata (Figure ​(Figure2A).2A). Users can specify parameters of the TCR sequence hom*ology search and/or pre-filter the VDJdb database to search in the specific data subset. Charts with summary statistics are supplied together with the analysis results and can be customized using different normalization options (Figure ​(Figure2B).2B). The results of repertoire annotation can be easily exported as tab-delimited tables and used for downstream analysis.

Open in a separate window

Figure 2.

Batch query interface of the VDJdb web browser. (A) TCR sequence annotations provided for an example AIRR-Seq sample with default matching criteria. The sample (HIP02877) represents an individual carrying a HLA-B*35 allele and is taken from the Emerson etal.'s study (17). Note a prominent EBV-specific clonal expansion restricted to this allele at the top of the annotations list. (B) Summary statistics charts comparing HIP02877 (‘B35+.txt’, left) to HIP13994 (‘CMV+.txt’, right) sample representing a CMV+ individual from the same study.

Evaluating and exploring TCR motifs

Similarity analysis of TCR sequences specific for a certain antigen reveals a complex network structure containing both large interconnected hubs of highly similar TCRs and a diverse set of unique TCR variants that cannot be co-clustered (6). Reducing the complexity of the VDJdb TCR network can both greatly speed-up database queries for large AIRR-Seq datasets and make it easier to visualize and interpret TCR binding motifs. The latter can be achieved by applying a recently developed de novo TCR motif discovery algorithm that allows to distinguish hubs of similar TCRs that arise due to antigen-specific enrichment from hubs that arise simply due to biases in the V(D)J rearrangement process (14).

Upon selecting such hubs in the VDJdb TCR network with a stringent criterion that only allows for a single substitution in the CDR3 region (see Material and Methods section) we obtained a set of 501 motifs specific for 40 epitopes for human and mouse TCR alpha and beta chains. Additionally, we prepared a set of normalized motif PWMs that control for amino acid biases for a given V/J combination and CDR3 length. This normalization removes germline-encoded residues that are not unique to a given motif: the conserved Cys and Phe/Trp and other flanking CDR3 residues that rarely interact with an epitope. Users can navigate the TCR motifs using the VDJdb motif browser (Motifs tab of the web interface) as shown in Figure ​Figure33 and query both full-length or partial CDR3 sequences against the database. Note that we provide motif PWMs to illustrate epitope-specific TCR patterns and do not report their accuracy for TCR specificity classification as this is beyond the scope of the present work.

Open in a separate window

Figure 3.

VDJdb motif browser interface. (A) Navigation tab showing the tree of available epitope motifs that can be selected to view PWMs of CDR3 amino acid sequences. The top plot shows the most abundant TCR beta chain motif for the A*02:GIL Influenza epitope. Motifs normalized for V(D)J rearrangement background are shown. (B) An example of CDR3 sequence query with matching amino acids highlighted. CDR3 sequence (CAEDNNARLMF) of the TCR alpha chain from the 3O4L PDB structure (TCR bound to A*02:GLC EBV epitope) was used as a query.

DISCUSSION

During the last two years we have greatly increased the number of entries in the VDJdb, and made substantial improvements to the database structure and web interface. Currently, VDJdb represents the largest open data set of TCRs with characterized specificity. We believe that the comprehensiveness of the VDJdb make it an attractive database for the benchmarking of TCR specificity prediction algorithms and other basic studies investigating TCR:peptide:MHC interactions. The future development of the database will primarily followtwo major directions. First, we aim at extending the set of existing records with putative antigen-specific TCRs using bioinformatic methods for de novo prediction of cognate TCR sequences. The latter is critical to provide a reasonable coverage for the overall diversity of an AIRR-Seq sample that can reach ∼106 TCR variants. Second, we aim at modelling TCR:peptide:MHC structures based on existing templates and VDJdb records in order to provide a basis for investigating the complexity of interactions involving residues of the TCR heterodimer and bring more insight into TCR recognition of antigens. We believe that the recent expansion of the VDJdb database both in terms of the underlying dataset and functionality, will make it a favorable resource for tasks related to annotation of T-cell AIRR-seq samples and basic studies of TCR specificity.

DATA AVAILABILITY

VDJdb database is available at https://vdjdb.cdr3.net.

ACKNOWLEDGEMENTS

We would like to thank all users of the VDJdb database for their valuable feedback.

Author contributions: D.V.B. developed VDJdb web server frontend. D.V.B. and M.S. developed VDJdb web server backend. R.V., J.S., U.S., C.R.R., G.D., A.G.W., M.A., E.S.E., I.V.Z., N.B., D.K.C., A.J.G., A.K.S., D.M.C., C.K., F.L. and M.S. processed, submitted and curated datasets for the VDJdb database. D.V.B., R.V., C.K., D.K.C., J.S., F.L., I.V.Z., D.M.C., U.S. and M.S. wrote the manuscript. M.S. supervised the study.

FUNDING

Russian Science Foundation [17-15-01495]. Funding for open access charge:Russian Science Foundation.

Conflict of interest statement. None declared.

REFERENCES

1. Pogorelyy M.V., Minervina A.A., Shugay M., Chudakov D.M., Lebedev Y.B., Mora T., Walczak A.M.. Detecting T cell receptors involved in immune responses from single repertoire snapshots. PLOS Biol. 2019; 17:e3000314. [PMC free article] [PubMed] [Google Scholar]

2. DeWitt W.S., Smith A., Schoch G., Hansen J.A., Matsen F.A., Bradley P.. Human T cell receptor occurrence patterns encode immune history, genetic background, and receptor specificity. eLife. 2018; 7:e38358. [PMC free article] [PubMed] [Google Scholar]

3. Rojas M., Restrepo-Jiménez P., Monsalve D.M., Pacheco Y., Acosta-Ampudia Y., Ramírez-Santana C., Leung P.S.C., Ansari A.A., Gershwin M.E., Anaya J.-M.. Molecular mimicry and autoimmunity. J. Autoimmun. 2018; 95:100–123. [PubMed] [Google Scholar]

4. Jiang N., Schonnesen A.A., Ma K.-Y.. Ushering in integrated T Cell repertoire profiling in cancer. Trends Cancer. 2019; 5:85–94. [PMC free article] [PubMed] [Google Scholar]

5. Bentzen A.K., Marquard A.M., Lyngaa R., Saini S.K., Ramskov S., Donia M., Such L., Furness A.J.S., McGranahan N., Rosenthal R. et al... Large-scale detection of antigen-specific T cells using peptide-MHC-I multimers labeled with DNA barcodes. Nat. Biotechnol. 2016; 34:1037–1045. [PubMed] [Google Scholar]

6. Shugay M., Bagaev D.V., Zvyagin I.V., Vroomans R.M., Crawford J.C., Dolton G., Komech E.A., Sycheva A.L., Koneva A.E., Egorov E.S. et al... VDJdb: a curated database of T-cell receptor sequences with known antigen specificity. Nucleic Acids Res. 2018; 46:D419–D427. [PMC free article] [PubMed] [Google Scholar]

7. Tickotsky N., Sagiv T., Prilusky J., Shifrut E., Friedman N.. McPAS-TCR: a manually curated catalogue of pathology-associated T cell receptor sequences. Bioinforma. Oxf. Engl. 2017; 33:2924–2929. [PubMed] [Google Scholar]

8. Mahajan S., Vita R., Shackelford D., Lane J., Schulten V., Zarebski L., Jespersen M.C., Marcatili P., Nielsen M., Sette A. et al... Epitope specific antibodies and T Cell receptors in the immune epitope database. Front. Immunol. 2018; 9:2688. [PMC free article] [PubMed] [Google Scholar]

9. Meysman P., DeNeuter N., Gielis S., BuiThi D., Ogunjimi B., Laukens K.. On the viability of unsupervised T-cell receptor sequence clustering for epitope preference. Bioinformatics. 2019; 35:1461–1468. [PubMed] [Google Scholar]

10. Jurtz V.I., Jessen L.E., Bentzen A.K., Jespersen M.C., Mahajan S., Vita R., Jensen K.K., Marcatili P., Hadrup S.R., Peters B. et al... NetTCR: sequence-based prediction of TCR binding to peptide-MHC complexes using convolutional neural networks. 2018; 03 October 2018, preprint: not peer reviewed 10.1101/433706. [CrossRef]

11. Gielis S., Moris P., Neuter N.D., Bittremieux W., Ogunjimi B., Laukens K., Meysman P.. TCRex: a webtool for the prediction of T-cell receptor sequence epitope specificity. 2018; 22 July 2018, preprint: not peer reviewed 10.1101/373472. [CrossRef]

12. Bradley P., Thomas P.G.. Using T Cell receptor repertoires to understand the principles of adaptive immune recognition. Annu. Rev. Immunol. 2019; 37:547–570. [PubMed] [Google Scholar]

13. Rubelt F., Busse C.E., Bukhari S.A.C., Bürckert J.-P., Mariotti-Ferrandiz E., Cowell L.G., Watson C.T., Marthandan N., Faison W.J., Hershberg U. et al... Adaptive Immune Receptor Repertoire Community recommendations for sharing immune-repertoire sequencing data. Nat. Immunol. 2017; 18:1274–1278. [PMC free article] [PubMed] [Google Scholar]

14. Ritvo P.-G., Saadawi A., Barennes P., Quiniou V., Chaara W., Soufi K.E., Bonnet B., Six A., Shugay M., Mariotti-Ferrandiz E. et al... High-resolution repertoire analysis reveals a major bystander activation of Tfh and Tfr cells. Proc. Natl. Acad. Sci. U.S.A. 2018; 115:9604–9609. [PMC free article] [PubMed] [Google Scholar]

15. Bagaev D.V., Zvyagin I.V., Putintseva E.V., Izraelson M., Britanova O.V., Chudakov D.M., Shugay M.. VDJviz: a versatile browser for immunogenomics data. BMC Genomics. 2016; 17:453. [PMC free article] [PubMed] [Google Scholar]

16. Pogorelyy M.V., Shugay M.. A framework for annotation of antigen specificities in high-throughput T-cell repertoire sequencing studies. Front. Immunol. 2019; 10:2159. [PMC free article] [PubMed] [Google Scholar]

17. Emerson R.O., DeWitt W.S., Vignali M., Gravley J., Hu J.K., Osborne E.J., Desmarais C., Klinger M., Carlson C.S., Hansen J.A. et al... Immunosequencing identifies signatures of cytomegalovirus exposure history and HLA-mediated effects on the T cell repertoire. Nat. Genet. 2017; 49:659–665. [PubMed] [Google Scholar]

Articles from Nucleic Acids Research are provided here courtesy of Oxford University Press

VDJdb in 2019: database extension, new analysis infrastructure and a T-cell receptor motif compendium (2024)

References

Top Articles
The Best Items Against Mordekaiser in League of Legends
67 Free Spirit Quotes for Those Who Live Differently
U Arizona Phonebook
Proko Funeral Home Obituary
Can I order a cake online from Stater Bros? - Chef's Resource
Discount Tire Parkdale Mall
Did Madelyn Cline Get Plastic Surgery? The Young Actress Looks Way Too Gorgeous These Days! - Celebs In-depth
Bhm Urban Dictionary
Pnc.comnc
Craigslist Farm And Garden Farmington Nm
Best Restaurants In The Financial District Sf
Best Boxing Gyms Near Me
Termos de Uso e Direitos Autorais - Responsabilidade Social
Framingham Line Mbta
Google Maps VS Waze : quelle est la meilleure application GPS ?
Bible Gateway passage: Romans 12 - New American Standard Bible 1995
Volvo All Engines (2010 Emissions) Fault Code: SPN 97 FMI 3 PID: Water In Fuel Indicator - EcuForce Truck
5 Ways to Get Overplucked Eyebrows to Grow Back
Acris Transfer Tax Calculator
Craigslist Philadelphia Craigslist: Your Guide to the City's Online Marketplace - First Republic Craigslist
i - Imaginary unit - MATLAB
Doodle Cricket Unblocked
Ncaa Predictions Cbs
Little Known Tools To Find Lenovo Parts - Upgrades And Options
Love and Other Alien Experiences
Thrillers won by Texas A&M have defined Southwest Classic at AT&T Stadium over past decade
Full Body Massage Near Me - Find Full Body Massage Places on Booksy.com! [US]
Quizwiz - Ace Your Homework & Exams, Now With ChatGPT AI
Page 2679 – Christianity Today
The Heartbreaking Story Of The Kind Vet Who Was Shot
Www Maurices Comenity
Collision Masters Fairbanks
Craigslist Las Vegas Rvs For Sale By Owner
Disinfection Services Lawton OK | Disinfection Services Cleaning Company
Trulia Home Sales
Susan Semetsis Obituary
Siena Weersverwachting
Bee And Willow Bar Cart
Used Gas Scooters For Sale Near Me For Adults
Where to Watch and Stream 'The Blackening': Showtimes, Digital, and Blu-ray
Mason County Mapsifter
Frank 26 Forum
Sodium Carbonate Sds Fisher
Chris and Rich Robinson swore never to speak again. But for the Black Crowes, rock heals all wounds
London Ups Store
Basis Phoenix Primary Calendar
Razor Edge Gotti Pitbull Price
popsike.com - very rare GENE VINCENT EP Made in France vg+
QPR 2-0 Hull City: Chris Willock and Ilias Chair net in victory
First United Bank Seminole Tx
Latest Posts
Article information

Author: Nathanial Hackett

Last Updated:

Views: 6529

Rating: 4.1 / 5 (52 voted)

Reviews: 91% of readers found this page helpful

Author information

Name: Nathanial Hackett

Birthday: 1997-10-09

Address: Apt. 935 264 Abshire Canyon, South Nerissachester, NM 01800

Phone: +9752624861224

Job: Forward Technology Assistant

Hobby: Listening to music, Shopping, Vacation, Baton twirling, Flower arranging, Blacksmithing, Do it yourself

Introduction: My name is Nathanial Hackett, I am a lovely, curious, smiling, lively, thoughtful, courageous, lively person who loves writing and wants to share my knowledge and understanding with you.