- Linux OS and the g++ compiler
- Python 3 (to generate certain output)
Please download TCRMatch from the latest release and compile from source by running:
make- 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 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 |
- 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
"Towards a Mathematical Foundation of Immunology and Amino Acid Chains" arXiv:1205.6031v2