Exploring QuiXoT features

.Here we describe some of the features of QuiXoT you can use in everyday's work analysing quantitative proteomics experiments. If this is the first time you run the program, you might be interested first in checking the unit tests for QuiXoT.

Analysis 1: 18O, HR

We will check some features using the 18O experiment using high resolution spectrometry, which we can see in the Test 1 of the unit test. If you didn't generate the QuiXML file (including quantification and statistics) following those steps, you can use the file VSMC_QuiXML_bs_quant_stats.xml in the VSMC_result folder (remember to move this file together with the binStack if you want to see the spectra; or just use the QuiXML alone if you do not need them).

1.1 Getting started

Open QuiXoT.exe, drag and drop anywhere on its window the abovementioned QuiXML file, and select the 18, HR, SEQUEST strategy. You should see this window:

1.2 Managing spectra

We might want to start checking how are the spectra (if you did not copy the binStack, you can skip this and go to generating graphs). Click on button spectrum, and then click on any row in the datagrid. For the fourth row (with FirstScan = 8502) you should see something like this

The experimental spectrum is drawn blue colour, while the red colour is the theoretical prediction (taking into account the peptide sequence, the isotope distribution, the labelling, and the labelling efficiency). The horizontal "lids" of the quantified peaks depict the tolerance used to consider if the theoretical peak matched (or not) the experimental one (note that there is a minimum "lid size" to visualise it). We can enlarge it to see better:

Or select the first four quantified peaks with the mouse (click and drag horizontally):

In this view we can see that apparently no co-eluted peaks are present.

Note that on the top of the screen we have an indicator of which is the precursor mass of the peptide that has been matched:

Let's look at a specific spectrum. To make it faster, we will filter the spectrum with FirstScan = 3017 (write the filter in the filter field and then click on the filter button:

The open the spectrum, enlarge

And then select the third quantified (red) peak:

(Note the warning on the top: the precursor is outside this window.)

This is a clear example of a co-eluted peak, so we can label the corresponding peptide to filter it out of the statistics later. We can go to the Label4 column, and directly write down s_coelution (or any other tag we prefer, as far as we can filter it easily later; starting the tags related to scans as s_, those related to peptides as p_, and those related to proteins as q_ is a good practice to distinguish between filters):

We can right-click on the spectrum to either zoom out or export data to a text file that can be treated by another software.

1.3 Checking information

Bottom left panels

Close the spectrum-window, and remove the filter in the filter field (delete and click on filter button).

You can see that in the bottom left corner there is a lot of information about every quantification. You can click with the mouse on any row, and go down with the arrow button: this information shows the data of each scan-peptide-protein of every row:

There are two panels, one for the quantitative information of spectra, peptides and proteins, and another with the identification information (prior to QuiXoT analysis). The details of this are explained in DataGrid information in QuiXoT.

Central bottom bars

Next to these panels there are three bars with important information for 18O experiments:

It is easier to explain what they are showing by an example. For example, if we select the spectrum with FirstScan = 1492,

we see it looks like the non-labelled species (first sample) is slightly more abundant than the labelled one (second sample). However, QuiXoT analysis shows that in this spectrum, the labelling efficiency was f = 0.8989. This means that there is a 10% of peptides coming from the second sample which failed to be labelled, and hence are added up to the non-labelled peptides. Correcting this, QuiXoT assigns a A = 100,501, and B = 101,294, which leads to a ratio of 0.992 and a log2Ratio = -0.113 (so, actually the second sample is slightly more abundant than the first one).

A more extreme case is shown for FirstScan = 3041:

also here looks like the two samples are almost equally abundant. However, it seems that the identified peptide is very poorly labelled, with labelling efficiency f = 0.3925. This means that most of the peptide in the second sample has not been labelled, adding itself up to the first, non-labelled sample. The QuiXoT analysis shows that abundance in the first sample is A = 4775, and for the second is B = 28,740, so that the log2Ratio = -2.59.

This is shown in the bottom-central bars:

The meaning of these bars is:

• the first bar (bluish green) shows the corrected abundance of the first sample,
• the second bar diplays the corrected abundance of the second sample, split in three colours:

• purple for the amount that has been fully labelled, i.e., with two 18O (both available 16O oxygen atoms in the carboxylic group have been replaced by 18O), so the labelled peptides have been added in the spectrum counts with a 4 Dalton separation
• yellow for the amount that has been partially labelled, i.e., labelled with only one 18O (only one of the two available 16O oxygen atoms in the carboxylic group has been replaced by 18O), so the labelled peptides have been added in the spectrum with a 2 Dalton separation (being superimposed to two of the non-labelled isotopologues)
• red, for the amount that has been fully non-labelled, i.e. despite the peptides come from the second sample, they are superimposed to the peptides from the first sample

• the third bar shows labelling efficiency, f, which is a ratio between 0 and 1; when f > 0.6, it is green colour, and when it is <= 0.6 it is red.

Bottom right panel

In the bottom right panel we have information about the statistical analysis:

Here we can see:

• the weight calibration constant K; we see we have used here the recommended value for these experiments (K = 40),
• the variances for each of the three levels taken into consideration:

• sigma2s for the variance of the spectrum level,
• sigma2p for the variance of the peptide level,
• sigma2q for the variance of the protein level.

This panel includes also controls for statistics:

• var conf, to calculate the limits of confidence of the variances,
• change values, to manually modify the variances and the weight calibration constant,
• var calc, to calculate automatically the variances,
• stats, to perform the satistics.

Top right panels

There are two small panels on the top right area:

The sample data panel, as it can be seen, has information about the number of scans, peptides and proteins in the analysis.

The load status panel is for advanced users, and allows reloading specific files. The files for isotopes, amino acids and QuantitationMethods refer to XML configuration files present in the conf folder of QuiXoT.

Finally, there are three buttons in the middle of the main window:

They are especially useful in order for the display of spectra, peptides and proteins. By default all the spectra are shown, so in every row, peptides and proteins (and their associated information, such as Xp, Wp, etc) appear repeated as many times as needed. However, this might not be desirable when we want to focus on peptides or proteins. Pressing the peptides or proteins uses scan or peptide indices to select only the first scan of each peptide or the first peptide of each protein,

1.4 Handling graphs

An interesting feature is the way graphs can be handled. We will start representing the spectrum weight (Ws) vs the log2Ratio (Xs). First, click on spectrum, and request the corresponding graph (don't forget to click on the scans button if you clicked on peptides or proteins in the previous subsection):

The resulting graph should be:

You might want to change some things. For example, centring the axis, or using a different notation for the axis. In this case, you can right-click on the graph, chose select bounds, and change the values of limits, X tick format and max # of X ticks to:

Then, we should obtain this:

Now, let's assume we want to know to which spectra are associated the four datapoints in the middle left. We can select the corresponding area with the mouse to zoom it:

The data are automatically filtered in the datagrid, so we can examine them,