Peak Frequency Test

From MicrobeMS Wiki
Revision as of 17:40, 16 December 2024 by Laschp (talk | contribs) (→‎Output of a peak frequency test)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Introduction

Peak frequency tests are useful for biomarker screening in ensembles of mass spectra exhibiting certain degree of similarity. These tests require pre-processed spectra with valid peak tables as inputs and are suitable to systematically determine the number of peaks and their relative abundance throughout the spectral m/z intervals. For peak frequency tests it is required to compile a mass spectral data set and to perform pre-processing including peak picking using standardized parameters. Peak frequency tests include a systematic evaluation for peak abundance in specific m/z segments which width is defined by a program-wide parameter dx (sometimes denoted as ppm tolerance). The output of a peak frequency test is a table in which ranked absolute and relative peak frequency values are listed. In addition, the relative peak frequency is plotted for each m/z segment (see screenshot below).

Performing a peak frequency test

1. Load the mass spectral data files via the load spectra (Bruker data file format), import spectra from mzXML data, or the load MS multifile options of the File pulldown menu.
2. Perform spectral pre-processing and peak detection.
3. Note that spectra selected for peak frequency tests are required to contain valid peak tables. Spectra without associated peak table cannot be processed.
4. Select the respective mass spectra in the listbox at the top left corner (the listbox is labeled by MicrobeMS spectra ID`s). To select multiple spectra hold the <shift> key while selecting. Alternatively, frequency tests can be also carried out from labeled spectra, i.e. from spectra with a class assignment. Labeling, or class assignment, can be carried out by selecting the appropriate spectra and choosing class assignments --> class X from the Edit pulldown menu.
5. If necessary define the parameter dx (ppm tolerance) which has a default value of dx = 1000 (relative ppm units). The parameter dx defines the width of m/z segments in which spectra are divided when analyzing peak frequencies. Peaks found in the same m/z segment are considered identical while mass peaks in different segments are considered different peaks. An spectral segment at the position x has width of x * dx/10^6 and ranges from [x*(1-dx/(2*10^6))] to [x*(1+dx/(2*10^6))], for example from m/z 1999 to 2001 if x = 2000 (m/z) and dx = 1000.
6. When finished select peak frequency plots from the Analysis pulldown menu. Choose options from selection or from class X if the peak frequency test should involve selected spectra or spectra with an appropriate class labeling, respectively.

Output of a peak frequency test


Example of a peak frequency table taken from the log file of MicrobeMS:

Command line output of peak frequency tests
Plot of the peak frequency vs. m/z dependence