Home  · Classes  · Annotated Classes  · Modules  · Members  · Namespaces  · Related Pages
IsobaricAnalyzer

Extracts and normalizes isobaric labeling information from an MS experiment.

pot. predecessor tools $ \longrightarrow $ IsobaricAnalyzer $ \longrightarrow $ pot. successor tools
FileConverter IDMapper
FileFilter

Extract the isobaric reporter ion intensities (currently iTRAQ 4plex and 8plex and TMT 6plex are supported) from raw MS2 data, does isotope corrections and stores the resulting quantitation as consensusXML, where each consensus centroid corresponds to one isobaric MS2 scan (e.g., HCD). The position of the centroid is the precursor position, its sub-elements are the channels (thus having m/z's of 113-121 for iTRAQ 126-131 for TMT respectively).

Isotope correction is done using non-negative least squares (NNLS), i.e.,

Minimize ||Ax - b||, subject to x >= 0, where b is the vector of observed reporter intensities (with 'contaminating' isotope species), A is a correction matrix (as supplied by the manufacturer AB Sciex) and x is the desired vector of corrected (real) reporter intensities. Other software solves this problem using an inverse matrix multiplication, but this can yield entries in x which are negative. In a real sample, this solution cannot possibly be true, so usually negative values (= negative reporter intensities) are set to 0. However, a negative result usually means, that noise was not accounted for thus we use NNLS to get a non-negative solution, without the need to truncate negative values. In (the usual) case that inverse matrix multiplication yields only positive values, our NNLS will give the exact same optimal solution.

The correction matrices can be found (and changed) in the INI file. However, these matrices for both 4plex and 8plex iTRAQ are now stable, and every kit delivered should have the same isotope correction values. Thus, there should be no need to change them, but feel free to compare the values in the INI file with your kit's Certificate. For TMT 6plex the values have to be adapted for each kit.

After this quantitation step, you might want to annotate the consensus elements with the respective identifications, obtained from an identification pipeline. Note that quantification is solely on peptide level at this stage. In order to obtain protein quantifications, you can try TextExporter to obtain a simple text format which you can feed to other software tools (e.g., R), or you can try ProteinQuantifier.

The command line parameters of this tool are:

IsobaricAnalyzer -- Calculates isobaric quantitative values for peptides
Version: 1.11.1 Nov 14 2013, 11:18:15, Revision: 11976

Usage:
  IsobaricAnalyzer <options>

This tool has algoritm parameters that are not shown here! Please check the ini file for a detailed descripti
on or use the --helphelp option.

Options (mandatory options marked with '*'):
  -type <mode>       Isobaric Quantitation method used in the experiment. (default: 'itraq4plex' valid: 'itra
                     q4plex', 'itraq8plex', 'tmt6plex')
  -in <file>*        Input raw/picked data file  (valid formats: 'mzML')
  -out <file>*       Output consensusXML file with quantitative information (valid formats: 'consensusXML')
                     
Common TOPP options:
  -ini <file>        Use the given TOPP INI file
  -threads <n>       Sets the number of threads allowed to be used by the TOPP tool (default: '1')
  -write_ini <file>  Writes the default configuration file
  -id_pool <file>    ID pool file to DocumentID's for all generated output files. Disabled by default. (Set 
                     to 'main' to use /Users/aiche/dev/openms/openms-1.11.1/share/OpenMS/IDPool/IDPool.txt)
  --help             Shows options
  --helphelp         Shows all options (including advanced)

The following configuration subsections are valid:
 - extraction       Parameters for the channel extraction.
 - itraq4plex       Algorithm parameters for itraq4plex
 - itraq8plex       Algorithm parameters for itraq8plex
 - quantification   Parameters for the peptide quantification.
 - tmt6plex         Algorithm parameters for tmt6plex

You can write an example INI file using the '-write_ini' option.
Documentation of subsection parameters can be found in the doxygen documentation or the INIFileEditor.
Have a look at the OpenMS documentation for more information.

INI file documentation of this tool:

Legend:
required parameter
advanced parameter
+IsobaricAnalyzerCalculates isobaric quantitative values for peptides
version1.11.1 Version of the tool that generated this parameters file.
++1Instance '1' section for 'IsobaricAnalyzer'
typeitraq4plex Isobaric Quantitation method used in the experiment.itraq4plex,itraq8plex,tmt6plex
in input raw/picked data file input file*.mzML
out output consensusXML file with quantitative informationoutput file*.consensusXML
log Name of log file (created only when specified)
debug0 Sets the debug level
threads1 Sets the number of threads allowed to be used by the TOPP tool
no_progressfalse Disables progress logging to command linetrue,false
id_pool ID pool file to DocumentID's for all generated output files. Disabled by default. (Set to 'main' to use /Users/aiche/dev/openms/openms-1.11.1/share/OpenMS/IDPool/IDPool.txt)
testfalse Enables the test mode (needed for internal use only)true,false
+++extractionParameters for the channel extraction.
select_activationHigh-energy collision-induced dissociation Operate only on MSn scans where any of its precursors features a certain activation method (e.g., usually HCD for iTRAQ). Set to empty string if you want to disable filtering.Collision-induced dissociation,Post-source decay,Plasma desorption,Surface-induced dissociation,Blackbody infrared radiative dissociation,Electron capture dissociation,Infrared multiphoton dissociation,Sustained off-resonance irradiation,High-energy collision-induced dissociation,Low-energy collision-induced dissociation,Photodissociation,Electron transfer dissociation,
reporter_mass_shift0.1 Allowed shift (left to right) in Da from the expected position.1e-08:0.5
min_precursor_intensity1 Minimum intensity of the precursor to be extracted. MS/MS scans having a precursor with a lower intensity will not be considered for quantitation.0:∞
keep_unannotated_precursortrue Flag if precursor with missing intensity value or missing precursor spectrum should be included or not.true,false
min_reporter_intensity0 Minimum intenesity of the individual reporter ions to be used extracted.0:∞
discard_low_intensity_quantificationsfalse Remove all reporter intensities if a single reporter is below the threshold given in min_reporter_intensity.true,false
min_precursor_purity0 Minimum fraction of the total intensity in the isolation window of the precursor spectrum attributable to the selected precursor.0:1
precursor_isotope_deviation0.02 Maximum allowed deviation between theoretical and observed isotopic peaks of the precursor peak in the isolation window to be counted as part of the precursor.0:∞
+++itraq4plexAlgorithm parameters for itraq4plex
channel_114_description Description for the content of the 114 channel.
channel_115_description Description for the content of the 115 channel.
channel_116_description Description for the content of the 116 channel.
channel_117_description Description for the content of the 117 channel.
reference_channel114 Number of the reference channel (114-117).114:117
correction_matrix[0.0/1.0/5.9/0.2, 0.0/2.0/5.6/0.1, 0.0/3.0/4.5/0.1, 0.1/4.0/3.5/0.1] Override default values (see Documentation); use the following format: <-2Da>/<-1Da>/<+1Da>/<+2Da> ; e.g. '0/0.3/4/0' , '0.1/0.3/3/0.2'
+++itraq8plexAlgorithm parameters for itraq8plex
channel_113_description Description for the content of the 113 channel.
channel_114_description Description for the content of the 114 channel.
channel_115_description Description for the content of the 115 channel.
channel_116_description Description for the content of the 116 channel.
channel_117_description Description for the content of the 117 channel.
channel_118_description Description for the content of the 118 channel.
channel_119_description Description for the content of the 119 channel.
channel_121_description Description for the content of the 121 channel.
reference_channel113 Number of the reference channel (113-121). Please note that 120 is not valid.113:121
correction_matrix[0.00/0.00/6.89/0.22, 0.00/0.94/5.90/0.16, 0.00/1.88/4.90/0.10, 0.00/2.82/3.90/0.07, 0.06/3.77/2.99/0.00, 0.09/4.71/1.88/0.00, 0.14/5.66/0.87/0.00, 0.27/7.44/0.18/0.00] Override default values (see Documentation); use the following format: <-2Da>/<-1Da>/<+1Da>/<+2Da> ; e.g. '0/0.3/4/0' , '0.1/0.3/3/0.2'
+++quantificationParameters for the peptide quantification.
isotope_correctiontrue Enable isotope correction (highly recommended). Note that you need to provide a correct isotope correction matrix otherwise the tool will fail or produce invalid results.true,false
normalizationfalse Enable normalization of channel intensities with respect to the reference channel. The normalization is done by using the Median of Ratios (every channel / Reference). Also the ratio of medians (from any channel and reference) is provided as control measure!true,false
+++tmt6plexAlgorithm parameters for tmt6plex
channel_126_description Description for the content of the 126 channel.
channel_127_description Description for the content of the 127 channel.
channel_128_description Description for the content of the 128 channel.
channel_129_description Description for the content of the 129 channel.
channel_130_description Description for the content of the 130 channel.
channel_131_description Description for the content of the 131 channel.
reference_channel126 Number of the reference channel (126-131).126:131
correction_matrix[0.0/0.0/0.0/0.0, 0.0/0.0/0.0/0.0, 0.0/0.0/0.0/0.0, 0.0/0.0/0.0/0.0, 0.0/0.0/0.0/0.0, 0.0/0.0/0.0/0.0] Override default values (see Documentation); use the following format: <-2Da>/<-1Da>/<+1Da>/<+2Da> ; e.g. '0/0.3/4/0' , '0.1/0.3/3/0.2'

OpenMS / TOPP release 1.11.1 Documentation generated on Thu Nov 14 2013 11:19:24 using doxygen 1.8.5