Nesvilab / FragPipe

A cross-platform Graphical User Interface (GUI) for running MSFragger and Philosopher - powered pipeline for comprehensive analysis of shotgun proteomics data
http://fragpipe.nesvilab.org
Other
178 stars 37 forks source link

global normalization and interpreted output file #1225

Closed Ink-Mo closed 10 months ago

Ink-Mo commented 10 months ago

I am currently replicating the data analysis from the article "Integrated Proteogenomic Characterization of Clear Cell Renal Cell Carcinoma." I have encountered an issue where the normalization method mentioned in the article, Aij =(RCij/ MADi) * MAD0 +M0+log2(REFi) (REFik, was estimated using the weighted sum of the MS1 intensities of the top three most intense peptide ions,REFi = Mean(REFik,k=1,.,q))seems to differ from the global normalization method described in TMT-Integrator, which calculates yij'=(yij/m2)xm1+m0. I would like to confirm this.

Additionally, I have another question: If I have conducted a TMT-11plex experiment with all 11 labels applied to different samples, and I have imported the mzML data files into Fragpipe using the TMT10-bridge workflow, but at the end, I wish to exclude one channel for data analysis (removing one sample, leaving only 10 channels annotated), will this have any impact? What impact, if any, will it have on the final results for both a single TMT-11plex and two TMT-11plex experiments? My speculation is whether this would affect the global normalization process.

Lastly, I noticed that in the Quant module, "Use MS1 intensity“ and ”top3 ions" is selected by default. I would like to know if the TMT-report file does not report MS1 intensity value.

Thank you for your assistance.

anesvi commented 10 months ago

I am not sure I follow any of this. We used TMT-Integrator in the paper. It was an early version of TMT-Integrator, but nothing should be fundamentally different.

For the second question, if you have TMT 11 (that is, 10 samples + bridge sample, so 11 total), you should load TMT10-bridge workflow but change from 10 to 11. And specify the bridge sample ('reference' approach). The resulting table will have 10 samples (the bridge sample will not be there, since everything is calculated as ratio to bridge sample). You should not remove samples from the annotation.tsv file. If you do not want to use one of the samples, you should remove if from the final table. If you are concerned about global normalization, do not use it. We now use MD option (and almost never GN option). If you want to exclude the sample from GN, you can try annotating the sample you want to exclude as NA in the annotation.tsx. I think that way it will not be used.

"Lastly, I noticed that in the Quant module, "Use MS1 intensity“ and ”top3 ions" is selected by default. I would like to know if the TMT-report file does not report MS1 intensity value." Sorry I do not understand your question.

Alexey

Ink-Mo commented 10 months ago

I am not sure I follow any of this. We used TMT-Integrator in the paper. It was an early version of TMT-Integrator, but nothing should be fundamentally different.

For the second question, if you have TMT 11 (that is, 10 samples + bridge sample, so 11 total), you should load TMT10-bridge workflow but change from 10 to 11. And specify the bridge sample ('reference' approach). The resulting table will have 10 samples (the bridge sample will not be there, since everything is calculated as ratio to bridge sample). You should not remove samples from the annotation.tsv file. If you do not want to use one of the samples, you should remove if from the final table. If you are concerned about global normalization, do not use it. We now use MD option (and almost never GN option). If you want to exclude the sample from GN, you can try annotating the sample you want to exclude as NA in the annotation.tsx. I think that way it will not be used.

"Lastly, I noticed that in the Quant module, "Use MS1 intensity“ and ”top3 ions" is selected by default. I would like to know if the TMT-report file does not report MS1 intensity value." Sorry I do not understand your question.

Alexey

I am not sure I follow any of this. We used TMT-Integrator in the paper. It was an early version of TMT-Integrator, but nothing should be fundamentally different.

For the second question, if you have TMT 11 (that is, 10 samples + bridge sample, so 11 total), you should load TMT10-bridge workflow but change from 10 to 11. And specify the bridge sample ('reference' approach). The resulting table will have 10 samples (the bridge sample will not be there, since everything is calculated as ratio to bridge sample). You should not remove samples from the annotation.tsv file. If you do not want to use one of the samples, you should remove if from the final table. If you are concerned about global normalization, do not use it. We now use MD option (and almost never GN option). If you want to exclude the sample from GN, you can try annotating the sample you want to exclude as NA in the annotation.tsx. I think that way it will not be used.

"Lastly, I noticed that in the Quant module, "Use MS1 intensity“ and ”top3 ions" is selected by default. I would like to know if the TMT-report file does not report MS1 intensity value." Sorry I do not understand your question.

Alexey

I apologize for the confusion. In Fragpipe's Quant(Isobaric) module, the ratio to abundance conversion is default to have "Use MS1 intensity" and "Top 3 ions" checked. I would like to know where in the TMT report file or which file I can find the MS1 intensity values.

Thinks

anesvi commented 10 months ago

PSM.tsv have intensity column. Those intensities are used to estimate the protein abundance in the reference sample, and convert from ratio to abundance values ( shown in abundance_… files)

Get Outlook for iOShttps://aka.ms/o0ukef


From: Ink-Mo @.> Sent: Friday, August 25, 2023 10:06:08 PM To: Nesvilab/FragPipe @.> Cc: Nesvizhskii, Alexey @.>; Assign @.> Subject: Re: [Nesvilab/FragPipe] global normalization and interpreted output file (Issue #1225)

External Email - Use Caution

I am not sure I follow any of this. We used TMT-Integrator in the paper. It was an early version of TMT-Integrator, but nothing should be fundamentally different.

For the second question, if you have TMT 11 (that is, 10 samples + bridge sample, so 11 total), you should load TMT10-bridge workflow but change from 10 to 11. And specify the bridge sample ('reference' approach). The resulting table will have 10 samples (the bridge sample will not be there, since everything is calculated as ratio to bridge sample). You should not remove samples from the annotation.tsv file. If you do not want to use one of the samples, you should remove if from the final table. If you are concerned about global normalization, do not use it. We now use MD option (and almost never GN option). If you want to exclude the sample from GN, you can try annotating the sample you want to exclude as NA in the annotation.tsx. I think that way it will not be used.

"Lastly, I noticed that in the Quant module, "Use MS1 intensity“ and ”top3 ions" is selected by default. I would like to know if the TMT-report file does not report MS1 intensity value." Sorry I do not understand your question.

Alexey

I am not sure I follow any of this. We used TMT-Integrator in the paper. It was an early version of TMT-Integrator, but nothing should be fundamentally different.

For the second question, if you have TMT 11 (that is, 10 samples + bridge sample, so 11 total), you should load TMT10-bridge workflow but change from 10 to 11. And specify the bridge sample ('reference' approach). The resulting table will have 10 samples (the bridge sample will not be there, since everything is calculated as ratio to bridge sample). You should not remove samples from the annotation.tsv file. If you do not want to use one of the samples, you should remove if from the final table. If you are concerned about global normalization, do not use it. We now use MD option (and almost never GN option). If you want to exclude the sample from GN, you can try annotating the sample you want to exclude as NA in the annotation.tsx. I think that way it will not be used.

"Lastly, I noticed that in the Quant module, "Use MS1 intensity“ and ”top3 ions" is selected by default. I would like to know if the TMT-report file does not report MS1 intensity value." Sorry I do not understand your question.

Alexey

I apologize for the confusion. In Fragpipe's Quant(Isobaric) module, the ratio to abundance conversion is default to have "Use MS1 intensity" and "Top 3 ions" checked. I would like to know where in the TMT report file or which file I can find the MS1 intensity values.

Thinks

— Reply to this email directly, view it on GitHubhttps://github.com/Nesvilab/FragPipe/issues/1225#issuecomment-1694117206, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AIIMM6ZYHX2OK4J6N5I5UZTXXFKZBANCNFSM6AAAAAA36DLOVE. You are receiving this because you were assigned.Message ID: @.***>


Electronic Mail is not secure, may not be read every day, and should not be used for urgent or sensitive issues