Peak width estimation in mass spectra

Patent metadata

Field	Value
Publication number	US-11721534-B2
Application number	US-202117186650-A
Country	US
Kind code	B2
Filing date	Feb 26, 2021
Priority date	Jul 10, 2020
Publication date	Aug 8, 2023
Grant date	Aug 8, 2023

Abstract

The disclosure relates to a method for processing a mass spectrum, comprising: (i) providing the mass spectrum which contains a plurality of data pairs, each data pair being representative of a mass value or mass-related value on a mass scale or mass-related scale and an abundance value or abundance-related value associated with the respective mass value or mass-related value, (ii) selecting a plurality of intervals on the mass scale or mass-related scale, each interval containing a multitude of the said data pairs, (iii) for each interval, applying a first mathematical-statistical analysis to the said data pairs contained in the respective interval in order to derive an interval-specific peak width, and (iv) using the said interval-specific peak widths to determine an estimated peak width for each mass value or mass-related value on the mass scale or mass-related scale.

First claim

The invention claimed is:

1. A method for processing a mass spectrum, which contains a first group of signals including actual ionic mass peaks of analytical interest and a second group of signals including omnipresent or ubiquitous noise, such as chemical or electronic noise, and improving a discrimination of the first group of signals from the second group of signals in the mass spectrum, comprising:
providing the mass spectrum which contains a plurality of data pairs, each data pair being representative of a mass value or mass-related value on a mass scale or mass-related scale and an abundance value or abundance-related value associated with the respective mass value or mass-related value, the mass spectrum being acquired from a sample by a mass spectrometry instrument chosen from among the group including: linear time-of-flight mass analyzer, reflector time-of-flight mass analyzer, orthogonal time-of-flight mass analyzer, Fourier Transform Ion Cyclotron Resonance (FT ICR) mass analyzer, and mass analyzer of the Kingdon type,
selecting a plurality of intervals on the mass scale or mass-related scale, each interval containing a multitude of the said data pairs,
for each interval, applying a first mathematical-statistical analysis to the said data pairs contained in the respective interval in order to derive an interval-specific peak width,
using the said interval-specific peak widths to determine an estimated peak width for each mass value or mass-related value on the mass scale or mass-related scale, and - complementing the mass spectrum with the estimated peak widths and subjecting it to a second mathematical-statistical analysis, wherein the second mathematical-statistical analysis is a peak picking or peak detection algorithm used to differentiate between the first group of signals and the second group of signals in the mass spectrum, the first group of signals being chosen from among the group including: lipids, metabolites, glycans, peptides, and proteins.

2. The method of claim 1 , wherein the derivation of the interval-specific peak widths includes computing for each interval an interval-specific degree of reliability, and wherein the determination of the estimated peak widths uses the interval-specific peak widths weighted as a function of the respective degrees of reliability.

3. The method of claim 1 , wherein the estimated peak widths are determined by at least one of interpolating between the interval-specific peak widths, extrapolating from the interval-specific peak widths, fitting a curve to the interval-specific peak widths, and performing a regression analysis to derive a mathematical relation mapping a mass value or mass-related value to an estimated peak width.

4. The method of claim 3 , wherein the estimated peak widths are determined by linear or spline interpolation.

5. The method of claim 3 , wherein the estimated peak widths are determined by fitting of a constant, linear, or higher degree polynomial.

6. The method of claim 1 , wherein the first mathematical-statistical analysis includes computing an autocorrelation of the mass spectrum within the respective interval, wherein the autocorrelation is the mass spectrum's correlation with a shifted version of itself.

7. The method of claim 6 , wherein the interval-specific peak width is determined as a full width at half maximum (FWHM) by computing the maximum (R max ) and minimum (R min ) of the autocorrelation and finding the position where the autocorrelation crosses a level given by their arithmetic mean, (R max +R min )/2.

8. The method of claim 6 , wherein the computation of interval-specific degrees of reliability includes computing for each interval a signal power of the mass spectrum and a total variation of the autocorrelation.

9. The method of claim 6 , wherein the autocorrelation is computed for a lag ranging between 0 and 1 atomic mass units or Daltons on the mass scale or mass-related scale.

10. The method of claim 1 , wherein the values on the mass scale or mass-related scale are expressed in atomic mass units or Daltons.

11. The method of claim 10 , wherein the length of each interval is between 50 and 1000 atomic mass units or Daltons.

12. The method of claim 1 , wherein, prior to the first mathematical-statistical analysis, the mass spectrum within each interval is resampled to equidistant values on the mass scale or mass-related scale.

13. The method of claim 12 , wherein a constant spacing between the equidistant mass values or mass-related values is chosen for each interval such that the resampled mass spectrum within the respective interval consists of the same number of data pairs as the original mass spectrum.

14. The method of claim 12 , wherein a constant spacing between the equidistant mass values or mass-related values is larger than or equal to one thousandth of an atomic mass unit or one milli-Dalton.

15. The method of claim 1 , wherein the plurality of intervals are chosen to be contiguous or discontiguous on the mass scale or mass-related scale.

16. The method of claim 1 , wherein the method is carried out during a mass spectrometry imaging experiment.

17. The method of claim 1 , wherein providing the mass spectrum comprises aggregating a plurality of individual mass spectra acquired by the same mass spectrometer over a limited period of time.

18. The method of claim 17 , wherein the plurality of individual mass spectra is acquired from a sample comprising one of tissue, a tissue section, and an extraction from tissue.

19. The method of claim 17 , wherein the aggregating includes averaging.

20. A mass spectrometer having an operating system which is adapted and configured to execute a method according to claim 1 .

21. A method for processing a mass spectrum, which contains a first group of signals including actual ionic mass peaks of analytical interest and a second group of signals including omnipresent or ubiquitous noise, such as chemical or electronic noise, and improving a discrimination of the first group of signals from the second group of signals in the mass spectrum, comprising:
providing the mass spectrum which contains a plurality of data pairs, each data pair being representative of a mass value or mass-related value on a mass scale or mass-related scale and an abundance value or abundance-related value associated with the respective mass value or mass-related value, the mass spectrum being acquired from a sample by a mass spectrometry instrument chosen from among the group including: linear time-of-flight mass analyzer, reflector time-of-flight mass analyzer, orthogonal time-of-flight mass analyzer, Fourier Transform Ion Cyclotron Resonance (FT ICR) mass analyzer, and mass analyzer of the Kingdon type,
selecting a plurality of intervals on the mass scale or mass-related scale, each interval containing a multitude of the said data pairs,
for each interval, applying a first mathematical-statistical analysis to the said data pairs contained in the respective interval in order to derive an interval-specific peak width,
using the said interval-specific peak widths to determine an estimated peak width for a plurality of selected mass values or mass-related values on the mass scale or mass-related scale, and
complementing the mass spectrum with the estimated peak widths and subjecting it to a second mathematical-statistical analysis, wherein the second mathematical-statistical analysis is a peak picking or peak detection algorithm used to differentiate between the first group of signals and the

Assignees

• Bruker Daltonik Gmbh

Inventors

• Boskamp Tobias

Classifications

• H01J49/0036 Primary

  Step by step routines describing the handling of the data generated during a measurement · CPC title

• G01N33/4833 Primary

  of solid biological material, e.g. tissue samples, cell cultures (tissue in vivo A61B5/00; cell suspensions G01N33/48735) · CPC title

• G01N27/62 Primary

  by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode · CPC title

• G06F17/18

  for evaluating statistical data {, e.g. average values, frequency distributions, probability functions, regression analysis (forecasting specially adapted for a specific administrative, business or logistic context G06Q10/04)} · CPC title

• G06F17/15

  Correlation function computation {including computation of convolution operations (arithmetic circuits for sum of products per se, e.g. multiply-accumulators G06F7/5443; digital filters, e.g. FIR, IIR, adaptive filters H03H17/00)} · CPC title