Can the fragpipe analyze the dimethylation-peptides?

Question

Can the fragpipe analyze the dimethylation-peptides?

Closed this issue a month ago · 40 comments

Hi fcyu,
I know the DIANN can analyze the multiplexed DIA raw data (plexDIA). Can the Fragpipe do it?
If it can, how do I configure the parameters in the Linux workflow (attachment)?
Now I want to analyze the dimethyl labeling peptides from plexDIA (light：Δ0=28.0313，intermediate：Δ4=32.0564，heavy：Δ8=36.0757) using the Fragpipe in Linux. But I don't know how to set up the parameters.
Can you help me?

Answer 1 · 2024-06-04T15:14:46.000Z

Yes, FragPipe can. We have a "plex DIA" panel in the "Quant (DIA)" tab dedicated for that. Hover over the mouse and read the tooltip for more details.

Best,

Fengchao

Answer 2 · 2024-06-14T07:25:53.000Z

Hi Fengchao, I have set up the plexDIA parameter based on the "Quant (DIA)" tab dedicated for that for my dimethyl-peptides in the FragPipe in Linux(below screenshot1) . And the spectra library used in the analysis was constructed by AlphaPeptDeep. However, I run into a error when I run the script(screenshot2). I don’t know if I can correctly set up the light, medium or heavy paramerters. Can you help me? ***@***.*** ***@***.*** 发件人: Fengchao ***@***.***> 发送时间: 2024年6月11日 2:48 收件人: Nesvilab/FragPipe ***@***.***> 抄送: 宁小莲(Xiaolian Ning) ***@***.***>; Author ***@***.***> 主题: Re: [Nesvilab/FragPipe] Can the fragpipe analyze the dimethylation-peptides? (Issue #1609) 外部邮件/External Mail Closed #1609<#1609> as completed. ― Reply to this email directly, view it on GitHub<#1609 (comment)>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/AUSN4G7T3YSIUUMCNWF6WBTZGXYGXAVCNFSM6AAAAABIYC2LIOVHI2DSMVQWIX3LMV45UABCJFZXG5LFIV3GK3TUJZXXI2LGNFRWC5DJN5XDWMJTGEYDMMBXGQZDANI>. You are receiving this because you authored the thread.Message ID: ***@***.***>

Answer 3 · 2024-06-14T13:13:31.000Z

Could you upload your log file by replying to #1609

Thanks,

Fengchao

Answer 4 · 2024-06-15T06:29:27.000Z

Hi Fengchao,
Thanks your reply. The attachment is the log file.
log_2024-06-14_12-11-09.txt

Best

Answer 5 · 2024-06-15T14:18:36.000Z

You only enabled DIA-NN but it is the last step just for quantification. You need to run all steps starting from MSFragger. Here is a brief tutorial for the label-free DIA: https://fragpipe.nesvilab.org/docs/tutorial_DIA.html. The diaTOP_ABPP workflow is the one you should take a look at.

Best,

Fengchao

Answer 6 · 2024-06-16T12:22:25.000Z

Hi fucy,
I have built a spectral library, I think I just need to provide the library and DIA data file. In addition, I want to analyze the dimethyl-labeled DIA.
Best

Answer 7 · 2024-06-16T14:45:20.000Z

How did you build the library? Could you upload the log file?

It seems that there is Dimethyl string that is not recognized by DIA-NN, which makes sense since you didn't define it using the command line options. If you want the DIA-NN quant module to recognize Dimethyl, you need to type the additional command in the opt cmds box. Or, you could use our workflow to run from scratch.

Best,

Fengchao

Answer 8 · 2024-06-24T01:26:56.000Z

Hi fucy,
I used the DIANN in window to analyze the dimethyl-based DIA, it is ok by setting up modification in the command line options. but it doesn't work using the Fragpipe in linux.
The attachment is the log file to build the library using the AlphapeptDeep.
logfile.zip

Best

Answer 9 · 2024-06-24T01:46:24.000Z

I believe when you raw DIA-NN in Windows, you added additional command line options to let DIA-NN recognize your custom modifications in the spectral library. You should do the similar when running FragPipe in Linux.

Best,

Fengchao

Answer 10 · 2024-06-24T01:52:49.000Z

Hi fucy,
if I can add the custom modifications in the FragPipe in Linux, how I do it？ Can I copy the parameter and paste it in FragPipe? Please read the screenshots.

Best

Answer 11 · 2024-06-24T01:56:08.000Z

Yes, you should copy the DIA-NN's command to the cmd line opts box in the Quant (DIA) tab. And you should also uncheck plex DIA box, and remove all entries in the light, medium, and heavy boxes since you build the spectral library by your own.

Best,

Fengchao

Answer 12 · 2024-06-24T02:04:02.000Z

Thanks for your reply.
If I use your FragPipe workflow to run dimethyl-based DIA data from scratch, how I set up the parameters in Linux?
Now I want to built the library based on DDA fractions by FragPipe in Linux and analyze the DIA data.

Best

Answer 13 · 2024-06-24T02:11:42.000Z

Load the DIA_SpecLib_Quant.workflow in the workflow tab (remember to click "load").
Load both of your DDA and DIA data in the input LC-MS Files panel of the workflow tab.
Set your DDA files to the DDA data type, and your DIA files to the DIA data by using the set DDA and set DIA buttons.
Specify your fasta file in the database tab.
Add the dimethyl modifications to the variable modification panel in the MSFragger tab.
Enable the plex DIA checkbox, and configure the light and heavy labels in the plex DIA panel of the Quant (DIA) tab.
Go to the result tab, specify the path of the result folder, and click run.

Best,

Fengchao

Answer 14 · 2024-06-24T08:36:29.000Z

Hi fcyu,
I configured the parameters according to your recommends(read the screenshots belowed), are they correct?

Best

Answer 15 · 2024-06-24T12:28:54.000Z

I think the settings look good. If you are not sure, could send me the log file to take a look after FragPipe finishes.

Best,

Fengchao

Answer 16 · 2024-06-25T12:54:10.000Z

Hi fcyu,
Now I can run the FragPipe to analyze dimethyl-peptide DIA. But I have a little confuse that I configured the dimethyl as variable modification, but the log file reported that the 28.0313 was the fix modification(screenshot). why?

Best

Answer 17 · 2024-06-25T12:57:04.000Z

The library was bulit based on the dimethyl-peptides DDA rawdata using the FragPipe and all labels(light, intermedian, heavy) were configured as variable modification.

Answer 18 · 2024-06-25T12:57:15.000Z

Briefly speaking, because a better approach is setting the light as fixed, and the differences between the heavy/medium and the light as variable. FragPipe optimized that for you.

Check the messages a few lines below

Best,

Fengchao

Answer 19 · 2024-06-25T14:01:18.000Z

Thank you for your patient and detailed explanation. I understand now.

Best

Answer 20 · 2024-06-26T01:59:46.000Z

Hi fcyu,
I can run the Fragpipe to analyze the dimethyl-peptides DIA rawdata. I checked the results and found the label in N-terminal is not 'label-N-L(screenshot)', but is 'UniMod:1'. Why is it such?

Best

Answer 21 · 2024-06-26T02:02:33.000Z

UniMod:1 is another modification, Acetylation, which means that the peptide you were looking at did not have dimethyl. Could you check the psm.tsv and library.tsv to troubleshoot?

Thanks,

Fengchao

Answer 22 · 2024-06-26T02:25:30.000Z

The library was built using dimethyl-labeled DDA and the generated library have the dimethyl label.
Could it be that the N-terminus is acetylated, making it less reactive with dimethyl labeling?

Answer 23 · 2024-06-26T02:29:29.000Z

I apologize for the miscommunication. I meant you could check the same peptide ADQLTEEQIAEFK against the psm.tsv and library.tsv to see if it indeed have acetylation on N-term. What you showed is a different peptide.

Could it be that the N-terminus is acetylated, making it less reactive with dimethyl labeling?

Yes. That was why in the MSFragger settings, both N-term acetylation and dimethyl labeling were variable modifications (one or the other or none).

Best,

Fengchao

Answer 24 · 2024-06-26T02:58:02.000Z

Thank you for your detailed explanation. I have reviewed the labeling status of the same peptide in the psm.tsv file, as shown in the screenshot below:

Answer 25 · 2024-06-26T02:58:49.000Z

Thank you for your detailed explanation. I have reviewed the labeling status of the same peptide in the psm.tsv file, as shown in the screenshot below.

Answer 26 · 2024-06-26T02:59:14.000Z

The figure was not uploaded successfully.

Best,

Fengchao

Answer 27 · 2024-06-26T03:02:31.000Z

ningzaizai commented 22 days ago

Answer 28 · 2024-06-26T03:03:38.000Z

Thanks, so there are no issues in FragPipe. Good to know.

Best,

Fengchao

Answer 29 · 2024-06-26T03:08:06.000Z

Thank you very much for your patient and detailed explanation.

Answer 30 · 2024-06-26T06:49:01.000Z

Hi fcyu,
I now want to build a large DDA spectral library based on dimethylation for methylated DIA analysis. Should I only use dimethylation as the light label for building the spectral library in DDA, and then use Fragpipe for analysis where DIANN extracts ion chromatograms based on light, medium, and heavy label mass differences? Or do I need to label with light, medium, and heavy labels separately, collect DDA data for each, and then build the spectral library?

Answer 31 · 2024-06-26T12:23:11.000Z

Should I only use dimethylation as the light label for building the spectral library in DDA, and then use Fragpipe for analysis where DIANN extracts ion chromatograms based on light, medium, and heavy label mass differences?

Do you mean only generate light dimethyl labeled DDA data? Yes, it should work. But in the MSFragger search, since you load both DDA and DIA, you still need to specify light, medium, and heavy as variable modifications.

Or do I need to label with light, medium, and heavy labels separately, collect DDA data for each, and then build the spectral library?

Generate the DDA data with light, medium, and heavy labeling also works.

Best,

Fengchao

Answer 32 · 2024-06-26T13:49:31.000Z

I understand that I should use only the light labeling sample, then collect DDA fractions. When building spectra library using Fragpipe, MSFraggers should be set to use only the modifications corresponding to the light label (n28.0313 as variable modification, K28.0313 as fixed modification), and then generate a spectral library containing only the light label. In the plexDIA tab, set light, medium, and heavy labeling. Finally, DIANN generates corresponding spectral libraries for medium and heavy labels based on the mass differences between light and medium, and heavy labels. This way, I can simplify the experiment by using only the light label for labeling and acquiring DDA data.

Answer 33 · 2024-06-26T13:58:00.000Z

I understand that I should use only the light labeling sample, then collect DDA fractions.

If you want, you could use light, medium, and heavy samples to collect DDA fractions, respectively.

When building spectra library using Fragpipe, MSFraggers should be set to use only the modifications corresponding to the light label (n28.0313 as variable modification, K28.0313 as fixed modification), and then generate a spectral library containing only the light label.

No, you should still set all of light, medium, and heavy if you have loaded both DDA and DIA data, because there are three labels in DIA.

DIANN generates corresponding spectral libraries for medium and heavy labels based on the mass differences between light and medium, and heavy labels

It is not exactly how DIA-NN plexDIA works. DIA-NN only accept light-only spectral library. It did a "match-within-run" to detect the medium and heavy peptides.

I can simplify the experiment by using only the light label for labeling and acquiring DDA data.

You can only generate light DDA data. But remember to set all of light, medium, and heavy in the MSFragger search if you have loaded the DIA data.

Best,

Fengchao

Answer 34 · 2024-06-26T14:21:25.000Z

Thank you for your detailed reply!
I just tested the Fragpipe workflow by constructing spectra library using only the light label and only analyzing DIA data containing only the light label. However, strangely, although I only labeled with the light tag in DIA data, the results still contain information about medium and heavy labels

Answer 35 · 2024-06-26T14:28:45.000Z

Which file is the screenshot from? Do those medium and heavy peptides have non-zero intensity?

Thanks,

Fengchao

Answer 36 · 2024-06-26T14:55:15.000Z

The screenshot from the 'diann-output.tsv ' file.
Those medium and heavy peptides have intensity.

Answer 37 · 2024-06-26T14:56:57.000Z

Could you upload the file?

Thanks,

Fengchao

Answer 38 · 2024-06-27T00:22:47.000Z

Hi fcyu,
I have identified the issue: when I set the cmd-opts in the workflow file to --fixed-mod Dimethyl,28.0313,nK --channels Dimethyl,0,nK,0:0; Dimethyl,4,nK,4.0251:4.0251; Dimethyl,8,nK,8.0444:8.0444, even with only one label in the DIA data, the results show three labels. However, when I set only --fixed-mod Dimethyl,28.0313,nK --channels Dimethyl,4,nK,4.0251:4.0251, the corresponding labels are only Dimethyl-n-4 or Dimethyl-K-4 modifications.
Best

Answer 39 · 2024-06-27T00:39:50.000Z

The first command seems to be correct and the second one seems to be wrong, that was why you didn't see the light labels.

But, you should not add any command to DIA-NN's cmd opts if you enabled plex DIA, because FragPipe have taken care of that.

Check the DIA-NN commands in the log file and you will understand what I said.

Best,

Fengchao

Answer 40 · 2024-06-27T00:51:53.000Z

This analysis is based on the spectral library I built using DDA components, analyzing another set of DIA data. Therefore, this analysis is equivalent to having an existing spectral library. So I disabled plex DIA and configured the parameter in DIA-NN's cmd opts. Now, the result is ok.