Nesvilab / FragPipe

A cross-platform proteomics data analysis suite
http://fragpipe.nesvilab.org
Other
205 stars 38 forks source link

Problem with 'PhilosopherFilter' exit code #579

Closed KR-15 closed 2 years ago

KR-15 commented 2 years ago

Hi,

I'm currently using FragPiPe with those specific versions:

Version info: FragPipe version 17.1 MSFragger version 3.4 Philosopher version 4.1.0 (build 1635640020)

I am using the LFQ-phospho workflow for basic phosphoenrichment experiment with 2 groups (8vs8) on timsTOF Pro data.

As the LFQ-phospho workflow only accepts mzML files as an input, I previously used MSConvert to transform .d raw files into mzML files and I think it goes well for this part (I used the tutorial to convert .d into mzML of this page https://fragpipe.nesvilab.org/docs/tutorial_convert.html). Then I run the LFQ-phospho workflow analysis without changing any default parameters (maybe only the "min ions" parameter) and I get an error message probably related with PhilosopherFilter as the end of the analysis:

INFO[18:06:06] Executing Filter v4.1.0
INFO[18:06:06] Processing peptide identification files
FATA[18:06:06] Cannot decode packed binary. expected element type but have Process 'PhilosopherFilter' finished, exit code: 1 Process returned non-zero exit code, stopping

You can find attached the log file of the analysis.

log_2022-01-11_18-06-07.txt

If you need other information from me do not hesitate.

Thank you

Best,

KR

fcyu commented 2 years ago

Seems to be an error from Philosopher:

PhilosopherFilter [Work dir: D:\Orbidata\IsSE_MaAU\Results_FragPipe_Phospho_mzML_PTMProphet\TRI_4]
C:\Users\Admin\Desktop\FragPipe-16.0\fragpipe\tools\philosopher\philosopher.exe filter --sequential --prot 0.01 --tag rev_ --pepxml D:\Orbidata\IsSE_MaAU\Results_FragPipe_Phospho_mzML_PTMProphet\TRI_4 --protxml D:\Orbidata\IsSE_MaAU\Results_FragPipe_Phospho_mzML_PTMProphet\combined.prot.xml --razorbin D:\Orbidata\IsSE_MaAU\Results_FragPipe_Phospho_mzML_PTMProphet\DMSO_1\.meta\razor.bin
INFO[18:06:06] Executing Filter  v4.1.0                     
INFO[18:06:06] Processing peptide identification files      
FATA[18:06:06] Cannot decode packed binary. expected element type <msms_pipeline_analysis> but have <msms_run_summary> 
Process 'PhilosopherFilter' finished, exit code: 1
Process returned non-zero exit code, stopping

Felipe @prvst , can you take a look when you have time?

Thanks,

Fengchao

prvst commented 2 years ago

Hi @KR-15, this error means that PTMProphet created a malformed XML file, with some missing tags. To avoid this you have to reduce the number of threads, and run it again.

KR-15 commented 2 years ago

Hi prvst,

Thank you for your reply.

I have relaunch the exact same analysis with only 4 CPUs and 2 CPUs instead of 8 CPUs (first try). It seems that the analysis runs further but there is still a new error message almost at the end of the analysis (see log files below). The error message is the same for both. Just to precise that there is a DIA-NN analysis that are currently running in the same computer but not all available threads are used for it.

For curiosity, why the number of threads used can impact the proper creation of a specific file (here xml) ?

Thank you for your help.

KR

log_2022-01-12_18-41-28_4_CPU.txt log_2022-01-13_10-44-29_2_CPU.txt

fcyu commented 2 years ago

Hi KR,

Could you please send the file 211220_IsSE_MaAU_nE_LFQ-P_CL_n1-8_F3_DMSO_2886.mzML to us? I can give you a link to upload the data if you need it.

Best,

Fengchao

KR-15 commented 2 years ago

Hi fcyu,

Yes of course, could you please provide a safe link for the upload ?

Best,

KR

fcyu commented 2 years ago

Sounds good! Here (https://www.dropbox.com/sh/i1zeuv8510hd9i9/AADXlXWgoAWbdFQm4UD0k4f9a?dl=0) is the link.

Thanks,

Fengchao

KR-15 commented 2 years ago

I had a problem uploading in the dropbox so I send a filesender link with de mzML file in your personal email address sorry.

Thanks for your help.

KR

fcyu commented 2 years ago

Hi KR,

The file is well received.

I perform an analysis using human database and the default LFQ-MBR workflow. But I cannot reproduce the issue.

I need all of the files in your DMSO_1 folder (corresponding to run 211220_IsSE_MaAU_nE_LFQ-P_CL_n1-8_F3_DMSO_2886.mzML), your fasta file, and fragger.params file to debug.

Thanks,

Fengchao

KR-15 commented 2 years ago

Hi,

The error occured with the LFQ-phospho worfklow so maybe it can come from this particular worfklow ?

Best,

KR

fcyu commented 2 years ago

Yes, I am using the LFQ plus phospho. I guess I need at least your fasta file and fragger.params to reproduce the error.

Best,

Fengchao

KR-15 commented 2 years ago

Hi Fengchao,

I have just sent you all the files that you may need to reproduce the error in a new filesender link.

Thanks for your time and help.

KR

fcyu commented 2 years ago

Hi KR,

The files are well received. I see that, this time, you also included a .d.zip file, which reminds me that your data is from timsTOF Pro. And you do specified IM-MS during the FragPipe analysis, but there is no ion mobility array in your mzML file. I think that is the cause of the error you were having. (When I was trying to reproducing the issue, I didn't specify IM-MS. Then IonQuant won't look for the ion mobility array. That is also why I could not reproduce the error).

Can you confirm that your data is from ddaPASEF with ion mobility module enabled. If so, you need to convert your .d according to the tutorial here (https://fragpipe.nesvilab.org/docs/tutorial_convert.html). I see that you were not using scanSumming filter, which you should.

Or, you can use use .d without any converting.

Best,

Fengchao

fcyu commented 2 years ago

Hi KR,

I think we have figured out the cause if the issue. Please read https://github.com/Nesvilab/FragPipe/issues/581#issuecomment-1013909453 for details.

Best,

Fengchao