annotate_gcr: Add positional annotation to MSstats group comparison result

Description Usage Arguments Details Value Examples

Description

This function takes a peptide-level group comparison result as output by MSstats::groupComparison and adds positional annotation to each peptide feature.

Usage

1
annotate_gcr(name, MSstats.GC, label)

Arguments

MSstats.GC

MSstats group comparison result as generated by groupComparison()

label

Label (character) of the comparision

name

file name (character) of a fasta file in working directory

Details

This function maps all peptide features to a reference proteome using their stripped aa sequence and perfect string matching. In addition, it determines the N-terminal modification status of each feature.

Important prerequisites:

This functions expects the results of a MSstats peptide-level significance analysis as input, since positional proteomics is a completely peptide centric approach.

Peptide-level significance analysis can be enforced easily by copying the content of PeptideSequence to ProteinName prior to calling MSstats::dataProcess().

Value

This function returns a table (data.frame) in long format. For none proteotypic peptides this function returns multiple rows, each listing the peptide in a different positional context. This ensures that all possible peptide origins can be considered in the downstream analysis.

Examples

1
annotate_gcr("uniprot-proteome_UP000000589.fasta", sample.data, "control vs. GluC")

tobiasko/posprot documentation built on May 26, 2019, 5:33 a.m.