Skip to content

Latest commit

 

History

History
85 lines (68 loc) · 6.05 KB

README.md

File metadata and controls

85 lines (68 loc) · 6.05 KB
Raw

TCRMatch

Maintainer: Will Chronister, Austin Crinklaw

Please open any issues on Github

Requirements:

  • Linux OS and the g++ compiler
  • Python 3 (to generate certain output)

Installation:

Please download TCRMatch from the latest release and compile from source by running:

make

Usage:

Input

  • A text file containing newline separated CDR3beta amino acid sequences
ASSLRWGDEQY
SADRVGNTLY
ASSFGREQY
ASSQDQVQDTQY
ASSPGQGPYEQY
ATSSGTGAYEQY
ASEAGGSYEQY
ASGDRGADTQY
ASSQAGAYEQY

Alternatively, you may upload files using the AIRR tsv format specified here

To generate the initial scores, please run

./tcrmatch -i input_file -t num_threads [-s threshold] [-d /path/to/database] [-a] > output_file
  • -i specifies the path to the input file
  • -a specifies input file is in AIRR TSV format
  • -t specifies the number of threads to use if operating on more than 1 core
  • -s optional parameter to specify the threshold (default is .97, in alignment with manuscript)
  • -d optional parameter to specify where the database is located, point this to a file of newline separated CDR3b sequences to test your own private set

To update the IEDB_data.tsv file located in the data folder, navigate to the scripts directory and simply run

./update.sh

This creates an archive folder in the top directory "data_archive" and moves the current IEDB_data.tsv file to said directory. This file will be renamed to the current day_month_year_IEDB_data.tsv, and a new version will be downloaded to data/IEDB_data.tsv.

Output

  • Output file has 7 columns in TSV format.
  • First column is the user provided input sequence.
  • Second column contains a CDR3beta sequence from the IEDB scoring higher than the threshold
  • The third column contains the score based on the kernel similarity metric (ranging from threshold to .99)
  • The fourth column contains the epitopes associated with the matched IEDB sequence. If there are multiple, they will be comma separated.
  • The fifth column contains the receptor groups (these can be used to reference the sequences back to the IEDB. If the receptor group is 12843, you can find all the receptors associated by visiting https://iedb.org/receptor/12843)
  • The sixth column contains the name of the specific antigen for the matched receptor
  • The seventh column contains the name of the source organism for the antigen

The last two columns are not always available depending on the data that was curated.

*Note that due to the trimming of IEDB CDR3b sequences, there may be multiple receptor groups that map to the same trimmed sequence. In these cases, the receptor groups will be comma separated and each link to their own unique receptor on the IEDB.

input_sequence match_sequence score epitopes receptor_group antigen source_organism
ASSSRSSYEQY ASSNRSSYEQY 0.97 GILGFVFTL 16861 Matrix protein 1 Influenza A virus
ASSSRSSYEQY ASSARSSYEQY 0.98 GILGFVFTL,RAKFKQLL,GLCTLVAML 53258 Matrix protein 1,Trans-activator protein BZLF1,Transcriptional regulator IE63 homolog Influenza A virus,Human herpesvirus 4 (Epstein Barr virus),Human herpesvirus 4 (Epstein Barr virus)
ASSSRSSYEQY ASSSRASYEQY 0.98 GILGFVFTL,GILGFVFTL 55749,35208 Matrix protein 1,Matrix protein 1 Influenza A virus,Influenza A virus
ASSSRSSYEQY ASSERSSYEQY 0.97 AMKRHGLDNYREYSLGN 21200
ASSSRSSYEQY ASSKRSSYEQY 0.97 GILGFVFTL 41132 Matrix protein 1 Influenza A virus
ASSSRSSYEQY ASSRSSYEQY 0.98 FVCNLLLLFVTVYSHLLLV,GILGFVFTL 120740,13575 ORF3a protein [Severe acute respiratory syndrome coronavirus 2],Matrix protein 1 SARS-CoV2,Influenza A virus
ASSSRSSYEQY ASSQRSSYEQY 0.97 GILGFVFTL,GILGFVFTL 35204,54706 Matrix protein 1,Matrix protein 1 Influenza A virus,Influenza A virus

How it works:

  • TCRMatch implements a similarity metric as defined in arXiv:1205.6031v2
  • Each input sequence is compared to all TCR CDR3beta sequences characterized in the IEDB
  • Sequences that have a score greater than the threshold (.97 by default) are returned as matches

References:

Chronister, William D. et al. "Tcrmatch: Predicting T-Cell Receptor Specificity Based On Sequence Similarity To Previously Characterized Receptors". Frontiers In Immunology, vol 12, 2021. Frontiers Media SA, doi:10.3389/fimmu.2021.640725.

"Towards a Mathematical Foundation of Immunology and Amino Acid Chains" arXiv:1205.6031v2