Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

GTEx TMT sample mapping #59

Open
helenabea opened this issue Apr 20, 2024 · 1 comment
Open

GTEx TMT sample mapping #59

helenabea opened this issue Apr 20, 2024 · 1 comment

Comments

@helenabea
Copy link

Hello, I'm having a lot of difficulty understanding how to map pepquery output to TMT tag of GTEx so that I can correctly map peptide to sample_ID.

I know this issue was opened before, but I do not understand the answers, mainly because this is my first time using mass-spec data.

So basically, I have this peptide:

spectrum_title: SecondInstrument_Sample20_S2R23_061716_Fr04:25026:4
peptide: HQDWRSGEAEFPDEESDESSSTSETSEEEVQDR
modification: TMT_10-plex_of_peptide_N-term@0[229.1629]
exp_mass: 4042.6508257733763
pep_mass: 4042.68993351836
tol_ppm: 9.67369390842789
tol_da: 0.039107744983539305
isotope_error: 0.0
score: 17.858883971787446

I was able to track it to psm_rank.mgf file, which gives me the following information:

BEGIN IONS
TITLE=SecondInstrument_Sample20_S2R23_061716_Fr04:25026:4
PEPMASS=1011.669982910156 3234025.828125
RTINSECONDS=8040.620801
CHARGE=4+

However, I do not understand how I can map this to GTEx sample ID.
image
So basically, I do understand that SecondInstrument_Sample20 has samples from Esophagus - Muscularis, Nerve - Tibial, Spleen, Uterus, Liver, Heart - Left Ventricle, Skin - Sun Exposed (Lower leg). What I want is to correctly attribute peptide "HQDWRSGEAEFPDEESDESSSTSETSEEEVQDR" to the "Nerve - Tibial" sample, for example.

@wenbostar
Copy link
Collaborator

The GTEx dataset is generated using the SPS-MS3 workflow as described in this paper: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7575058/. In this workflow, peptide identification is based on MS2 spectra while quantification is based on MS3 spectra. You won't be able to get the quantification information from the spectra returned by PepQuery since the identification is based on MS2 spectra which don't contain quantification data. I'd like to suggest to use FragPipe for quantification. You could add the target peptides identified by PepQuery into a database then use FragPipe to quantify these peptides.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants