| OpenMS 2.5.0

FeatureFinderCentroided

The feature detection application for quantitation (centroided).

pot. predecessor tools | FeatureFinderCentroided | pot. successor tools
PeakPickerWavelet | FeatureLinkerUnlabeled (or another feature grouping tool)
SeedListGenerator | MapAlignerPoseClustering (or another alignment tool)

Reference:
Weisser et al.: An automated pipeline for high-throughput label-free quantitative proteomics (J. Proteome Res., 2013, PMID: 23391308).

This module identifies "features" in a LC/MS map. By feature, we understand a peptide in a MS sample that reveals a characteristic isotope distribution. The algorithm computes positions in rt and m/z dimension and a charge estimate of each peptide.

The algorithm identifies pronounced regions of the data around so-called `seeds`. In the next step, we iteratively fit a model of the isotope profile and the retention time to these data points. Data points with a low probability under this model are removed from the feature region. The intensity of the feature is then given by the sum of the data points included in its regions.

How to find suitable parameters and details of the different algorithms implemented are described in the TOPP tutorial.

Specialized tools are available for some experimental techniques: IsobaricAnalyzer.

The command line parameters of this tool are:

    FeatureFinderCentroided -- Detects two-dimensional features in LC-MS data.

    Full documentation: http://www.openms.de/documentation/TOPP_FeatureFinderCentroided.html

    Version: 2.5.0-GIT-NOTFOUND-GIT-NOTFOUND Sep 10 2020, 12:20:39, Revision: GIT-NOTFOUND

    To cite OpenMS:

      Rost HL, Sachsenberg T, Aiche S, Bielow C et al.. OpenMS: a flexible open-source software platform for mass spectrometry data analysis. Nat Meth. 2016; 13, 9: 741-748. doi:10.1038/nmeth.3959.

    Usage:

      FeatureFinderCentroided <options>

    This tool has algorithm parameters that are not shown here! Please check the ini file for a detailed descript

    ion or use the --helphelp option.

    Options (mandatory options marked with '*'):

      -in <file>*        Input file (valid formats: 'mzML')

      -out <file>*       Output file (valid formats: 'featureXML')

      -seeds <file>      User specified seed list (valid formats: 'featureXML')

    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

      --help             Shows options

      --helphelp         Shows all options (including advanced)

    The following configuration subsections are valid:

     - algorithm   Algorithm section

    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.

    For more information, please consult the online documentation for this tool:

      - http://www.openms.de/documentation/TOPP_FeatureFinderCentroided.html

INI file documentation of this tool:

Legend:
required parameter
advanced parameter
+FeatureFinderCentroided Detects two-dimensional features in LC-MS data.
version 2.5.0-GIT-NOTFOUND-GIT-NOTFOUND Version of the tool that generated this parameters file.
++1 Instance '1' section for 'FeatureFinderCentroided'
in input file input file*.mzML
out output file output file*.featureXML
seeds User specified seed list input file*.featureXML
out_mzq Optional output file of MzQuantML.output file*.mzq
log Name of log file (created only when specified)
debug 0 Sets the debug level
threads 1 Sets the number of threads allowed to be used by the TOPP tool
no_progress false Disables progress logging to command line true,false
force false Overwrite tool specific checks.true,false
test false Enables the test mode (needed for internal use only)true,false
+++algorithm Algorithm section
debug false When debug mode is activated, several files with intermediate results are written to the folder 'debug' (do not use in parallel mode).true,false
++++intensity Settings for the calculation of a score indicating if a peak's intensity is significant in the local environment (between 0 and 1)
bins 10 Number of bins per dimension (RT and m/z). The higher this value, the more local the intensity significance score is.
This parameter should be decreased, if the algorithm is used on small regions of a map.1:∞
++++mass_trace Settings for the calculation of a score indicating if a peak is part of a mass trace (between 0 and 1).
mz_tolerance 0.03 Tolerated m/z deviation of peaks belonging to the same mass trace.
It should be larger than the m/z resolution of the instrument.
This value must be smaller than that 1/charge_high!0.0:∞
min_spectra 10 Number of spectra that have to show a similar peak mass in a mass trace.1:∞
max_missing 1 Number of consecutive spectra where a high mass deviation or missing peak is acceptable.
This parameter should be well below 'min_spectra'!0:∞
slope_bound 0.1 The maximum slope of mass trace intensities when extending from the highest peak.
This parameter is important to seperate overlapping elution peaks.
It should be increased if feature elution profiles fluctuate a lot.0.0:∞
++++isotopic_pattern Settings for the calculation of a score indicating if a peak is part of a isotopic pattern (between 0 and 1).
charge_low 1 Lowest charge to search for.1:∞
charge_high 4 Highest charge to search for.1:∞
mz_tolerance 0.03 Tolerated m/z deviation from the theoretical isotopic pattern.
It should be larger than the m/z resolution of the instrument.
This value must be smaller than that 1/charge_high!0.0:∞
intensity_percentage 10.0 Isotopic peaks that contribute more than this percentage to the overall isotope pattern intensity must be present.0.0:100.0
intensity_percentage_optional 0.1 Isotopic peaks that contribute more than this percentage to the overall isotope pattern intensity can be missing.0.0:100.0
optional_fit_improvement 2.0 Minimal percental improvement of isotope fit to allow leaving out an optional peak.0.0:100.0
mass_window_width 25.0 Window width in Dalton for precalculation of estimated isotope distributions.1.0:200.0
abundance_12C 98.930000000000007 Rel. abundance of the light carbon. Modify if labeled.0.0:100.0
abundance_14N 99.632000000000005 Rel. abundance of the light nitrogen. Modify if labeled.0.0:100.0
++++seed Settings that determine which peaks are considered a seed
min_score 0.8 Minimum seed score a peak has to reach to be used as seed.
The seed score is the geometric mean of intensity score, mass trace score and isotope pattern score.
If your features show a large deviation from the averagene isotope distribution or from an gaussian elution profile, lower this score.0.0:1.0
++++fit Settings for the model fitting
max_iterations 500 Maximum number of iterations of the fit.1:∞
++++feature Settings for the features (intensity, quality assessment, ...)
min_score 0.7 Feature score threshold for a feature to be reported.
The feature score is the geometric mean of the average relative deviation and the correlation between the model and the observed peaks.0.0:1.0
min_isotope_fit 0.8 Minimum isotope fit of the feature before model fitting.0.0:1.0
min_trace_score 0.5 Trace score threshold.
Traces below this threshold are removed after the model fitting.
This parameter is important for features that overlap in m/z dimension.0.0:1.0
min_rt_span 0.333 Minimum RT span in relation to extended area that has to remain after model fitting.0.0:1.0
max_rt_span 2.5 Maximum RT span in relation to extended area that the model is allowed to have.0.5:∞
rt_shape symmetric Choose model used for RT profile fitting. If set to symmetric a gauss shape is used, in case of asymmetric an EGH shape is used.symmetric,asymmetric
max_intersection 0.35 Maximum allowed intersection of features.0.0:1.0
reported_mz monoisotopic The mass type that is reported for features.
'maximum' returns the m/z value of the highest mass trace.
'average' returns the intensity-weighted average m/z value of all contained peaks.
'monoisotopic' returns the monoisotopic m/z value derived from the fitted isotope model.maximum,average,monoisotopic
++++user-seed Settings for user-specified seeds.
rt_tolerance 5.0 Allowed RT deviation of seeds from the user-specified seed position.0.0:∞
mz_tolerance 1.1 Allowed m/z deviation of seeds from the user-specified seed position.0.0:∞
min_score 0.5 Overwrites 'seed:min_score' for user-specified seeds. The cutoff is typically a bit lower in this case.0.0:1.0
++++debug
pseudo_rt_shift 500.0 Pseudo RT shift used when .1.0:∞

For the parameters of the algorithm section see the algorithms documentation:
centroided
In the following table you can find example values of the most important parameters for different instrument types.
These parameters are not valid for all instruments of that type, but can be used as a starting point for finding suitable parameters.

'centroided' algorithm:

| Q-TOF | LTQ Orbitrap
intensity:bins | 10 | 10
mass_trace:mz_tolerance | 0.02 | 0.004
isotopic_pattern:mz_tolerance | 0.04 | 0.005

For the centroided algorithm centroided data is needed. In order to create centroided data from profile data use the PeakPickerWavelet.