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

No peptides and spectral but identified as a protein #513

Open
LiTing99-cutie opened this issue Nov 22, 2024 · 3 comments
Open

No peptides and spectral but identified as a protein #513

LiTing99-cutie opened this issue Nov 22, 2024 · 3 comments
Assignees

Comments

@LiTing99-cutie
Copy link

LiTing99-cutie commented Nov 22, 2024

Dear philosopher team,

After i ran a LFQ-MBR workflow using fragpipe, and inspected the result: i found a wierd row in the protein.tsv:sp|P27870|VAV_MOUSE P27870 VAV_MOUSE Vav1 845 Mus musculus Proto-oncogene vav 1:Experimental evidence at protein level 0.71 0.3871 0.8196 0 0 0 0 000 0 0,
since this protein has no peptides and spectra, why does it still identified as a protein?

Besides, i input 6 samples for running at a time, i notices the protein number would also be different from that of running one at a time. Which strategy should i used?

My platform and the philosopher version:

  • Linux
  • Version: v 5.1.1

Log file is in the attachments
log_2024-11-18_12-39-25.txt

@fcyu fcyu self-assigned this Nov 22, 2024
@fcyu
Copy link
Member

fcyu commented Nov 22, 2024

After i ran a LFQ-MBR workflow using fragpipe, and inspected the result: i found a wierd row in the protein.tsv:sp|P27870|VAV_MOUSE P27870 VAV_MOUSE Vav1 845 Mus musculus Proto-oncogene vav 1:Experimental evidence at protein level 0.71 0.3871 0.8196 0 0 0 0 000 0 0,
since this protein has no peptides and spectra, why does it still identified as a protein?

Can you share the file you mentioned? My guess is that the protein is from MBR.

Besides, i input 6 samples for running at a time, i notices the protein number would also be different from that of running one at a time. Which strategy should i used?

FDR control are different. Run 6 samples one at a time basically perform FDR filtering 6 times, but running 6 samples together perform one global FDR once. You should run 6 samples together.

Best,

Fengchao

@LiTing99-cutie
Copy link
Author

Of cource, i actually tested ionquant.mbr=0 and the row also existed, I uploaded the protein.txt in the attachment now.

Also when i used ionquant.mbr=0 and quantitation.run-label-free-quant=false, and the row turned into
sp|P27870|VAV_MOUSE P27870 VAV_MOUSE Vav1 845 Mus musculus Proto-oncogene vav 1:Experimental evidence at protein level 0.71 0.3871 0.8196 1 0 0 1 0 0 0 0 0, shows it has one total peptide and one spectral count. But when i grep VAV_MOUSE in the psm.tsv, ion.tsv and peptide.tsv, i got no result. It should be in the mapped protein colunm or something? I don't know why.

Also, i also noticed when all of the samples are inputted together, the razor peptide column of protein.tsv could be 0 sometimes. Is this normal? Without razor peptide but still identified as a protein?

protein.txt

@fcyu
Copy link
Member

fcyu commented Nov 23, 2024

Thanks for sharing the file. It seems to be a bug.....

Best,

Fengchao

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

No branches or pull requests

2 participants