Protein Quantification on Skyline

support
Protein Quantification on Skyline cataoliv  2016-07-19 08:18
 
Hello,

My name is Catarina and I am starting to work with SRM.
I have a doubçt regarding the protein quantification method used during Skyline analysis. I analyzed three biological replicates and three technical of each one in two conditions (Control vs Interaction) . If I perform a group comparison of Control and Interaction I obtain a table with a Fold change result and the adjusted P-value.
My question is which test is performed to obtain the p-value, once I have a couple of peptides that present fold changes in different directions and it still gives me a significant p-value, which does not fit with the raw results.
Hope to hearing from you soon,

Best,
 
 
Nick Shulman responded:  2016-07-19 11:36
The p-values come from a T-test.
A low p-value means that the observed variation among the replicates within a group was smaller than the variation between the two groups.
 
cataoliv responded:  2016-07-20 04:39
Thank you for your prompt response.
For the comparison of two conditions does the program use the paired t-test? Another question is if the t-test is assuming technical replicates and biological replicates or if it uses all the values all together indepently of being technical or biological replicates?
I do not know if I made myself clear, but when we do the t-test in excel we obtain different p-values, and this is why we have some doubts.

Thank you.
 
Nick Shulman responded:  2016-07-20 06:18
If you have specified something for the "Identity Annotation" in the group comparison, then the replicates get grouped together based on that annotation, and the values for each get averaged together before being used in the t-test.

If you are calculating the fold change "per peptide" (and not "per protein"), then you can use the document grid to see the normalized area values that are used in the t-test. Just make sure that the normalization method on:
Settings > Peptide Settings > Quantification
is the same as the normalization method specified for your group comparison.

In the Document Grid, you can customize the view and find the "Normalized Area" column (there's a search button on the toolbar in the View Editor-- it looks like binoculars). The "Normalized Area" column is under "Peptide Results/Quantification".

In order to calculate the fold change, Skyline takes the logarithm (base 2) of the normalized areas for each of the replicates. If there are technical replicates, then those log areas are averaged together for each group of technical replicates.

Then, a linear regression is performed where the x-value is 0 or 1 depending on whether the replicate was in the control group, and the y-value is the log area. Skyline raises 2 to the power of that slope, and that is the fold change for that peptide.

No, it is not a paired t-test.
 
cataoliv responded:  2016-07-20 06:34

Thank you.

We were thinking on doing a paired t-test because we wanted to compare our two conditions (control vs interaction) of each biological replicate. Is there any option to do it in the skyline?
 
Nick Shulman responded:  2016-07-20 06:52
You should use MSstats. MSstats was written in R and has access to more complicated statistics than what we have implemented in Skyline.

Here's the MSstats tutorial:
https://skyline.gs.washington.edu/labkey/wiki/home/software/Skyline/page.view?name=tutorial_grouped

I know that in MSstats if the same "BioReplicate" annotation value appears in the case and control group, then MSstats assumes that you are doing a time course experiment, and you get something like a paired t-test.