Testing quantms with large TMT dataset

· toanphung's blog


The proteomics dataset in picture for this analysis is from Target ALS. First attempt to analyze the dataset which include 224 plasma raw files and 198 csf raw files separately.

SDRF file creation #

Due to the large amount of files and since its is a TMT experiment too, a script was used to make the sdrf file

The columns used for the file

 1COLS = [
 2    'source name',
 3    'characteristics[organism]',
 4    'characteristics[organism part]',
 5    'characteristics[age]',
 6    'characteristics[sex]',
 7    'characteristics[cell type]',
 8    'characteristics[cell line]',
 9    'characteristics[disease]',
10    'characteristics[depletion]',
11    'characteristics[developmental stage]',
12    'characteristics[individual]',
13    'characteristics[pooled sample]',
14    'characteristics[biological replicate]',
15    'assay name',
16    'technology type',
17    'comment[technical replicate]',
18    'comment[data file]',
19    'comment[file uri]',
20    'comment[fraction identifier]',
21    'comment[label]',
22    'comment[instrument]',
23    'comment[proteomics data acquisition method]',
24    'comment[dissociation method]',
25    'comment[modification parameters]',
26    'comment[modification parameters]',
27    'comment[modification parameters]',
28    'comment[modification parameters]',
29    'comment[cleavage agent details]',
30    'comment[precursor mass tolerance]',
31    'comment[fragment mass tolerance]',
32    'comment[sdrf version]',
33    'comment[sdrf template]',
34    'factor value[disease]',
35]

Modification columns used the following

1MOD_PARAMS = [
2    'NT=TMTpro;AC=UNIMOD:2016;PP=Any N-term;MT=Fixed',
3    'NT=TMTpro;AC=UNIMOD:2016;TA=K;MT=Fixed',
4    'NT=Carbamidomethyl;AC=UNIMOD:4;TA=C;MT=Fixed',
5    'NT=Oxidation;AC=UNIMOD:35;TA=M;MT=Variable',
6]

Beyond that, the following are the other acquisition constants

1INSTRUMENT   = 'NT=Orbitrap Eclipse;AC=MS:1003028'
2ACQ_METHOD   = 'Data-Dependent Acquisition'
3DISSOCIATION = 'AC=MS:1000422;NT=HCD'
4PREC_TOL     = '10 ppm'
5FRAG_TOL     = '0.02 Da'

In factor value column Norm is used for pooled rows.

The database utilized was the current bersion at the time of the uniprot human swissprot fasta database in combination with a potential contaminant database.

Analysis run #

I ran into a problem due to the usage of TMT35plex in the original experiment with the additional deuterium channels and it seemed at the time that quantms currently support maximum TMT18plex and PRIDE does not have the ontology terms for the new deuterium channels. To make it work, I stripped the sdrf tables of any row with deuterium channels.

last updated: