OpenMS
|
Identifies precursor mass traces and tries to correlate them with fragment ion mass traces in SWATH maps.
This algorithm will try to correlate the masstraces to find co-eluting traces and cluster them.
This program looks at mass traces in a precursor MS1 map and tries to correlate them with features found in the corresponding MS2 map based on their elution profile. It uses
It does a separate correlation analysis on the MS1 and the MS2 map, both produces a set of pseudo spectra. In a second (optional) step, the MS2 pseudo spectra are correlated with the MS1 traces and the most likely precursor is assigned to the pseudo spectrum.
It is based on the following papers: ETISEQ – an algorithm for automated elution time ion sequencing of concurrently fragmented peptides for mass spectrometry-based proteomics BMC Bioinformatics 2009, 10:244 doi:10.1186/1471-2105-10-244 ; http://www.biomedcentral.com/1471-2105/10/244 they use FFT to correlate and then use lag of at least 1 scan and pearson correlation of 0.7 to assign precursors to product ions If one fragment matches to multiple precursors, it is assigned to all of them. If it doesn't match any, it is assigned to all
The command line parameters of this tool are:
ClusterMassTracesByPrecursor -- Correlate precursor masstraces with fragment ion masstraces in SWATH maps based on their elution profile. Full documentation: http://www.openms.de/doxygen/release/3.2.0/html/TOPP_ClusterMassTracesByPrecursor.html Version: 3.2.0 Sep 18 2024, 16:00:56, Revision: e231942 To cite OpenMS: + Pfeuffer, J., Bielow, C., Wein, S. et al.. OpenMS 3 enables reproducible analysis of large-scale mass spec trometry data. Nat Methods (2024). doi:10.1038/s41592-024-02197-7. Usage: ClusterMassTracesByPrecursor <options> Options (mandatory options marked with '*'): -in_ms1 <file>* MS1 mass traces (valid formats: 'consensusXML') -in_swath <file>* MS2 / SWATH mass traces (valid formats: 'consensusXML') -out <file>* Output file (valid formats: 'mzML') -assign_unassigned_to_all Assign unassigned MS2 fragments to all precursors (only for ms1_centrif) -min_pearson_correlation <double> Minimal pearson correlation score to match elution profiles to each othe r. (default: '0.7') -max_lag <number> Maximal lag (e.g. by how many spectra the peak may be shifted at most). This parameter will depend on your chromatographic setup but a number between 1 and 3 is usually sensible. (default: '1') -min_nr_ions <number> Minimal number of ions to report a spectrum. (default: '3') -max_rt_apex_difference <double> Maximal difference of the apex in retention time (in seconds). This is a hard parameter, all profiles further away will not be considered at all. (default: '5.0') -swath_lower <double> Swath lower isolation window (default: '0.0') -swath_upper <double> Swath upper isolation window (default: '0.0') 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)
INI file documentation of this tool: