reconstructChromPeakSpectra: Data independent acquisition (DIA): reconstruct MS2 spectra

View source: R/functions-XCMSnExp.R

reconstructChromPeakSpectraR Documentation

Data independent acquisition (DIA): reconstruct MS2 spectra

Description

Reconstructs MS2 spectra for each MS1 chromatographic peak (if possible) for data independent acquisition (DIA) data (such as SWATH). See the LC-MS/MS analysis vignette for more details and examples.

Usage

reconstructChromPeakSpectra(object, expandRt = 1, diffRt = 2,
  minCor = 0.8, intensity = "maxo",
  peakId = rownames(chromPeaks(object, msLevel = 1L)),
  BPPARAM = bpparam())

Arguments

object

XCMSnExp with identified chromatographic peaks.

expandRt

numeric(1) allowing to expand the retention time range for extracted ion chromatograms by a constant value (for the peak shape correlation).

diffRt

numeric(1) defining the maximal allowed difference between the retention time of the chromatographic peak (apex) and the retention times of MS2 chromatographic peaks (apex) to consider them as representing candidate fragments of the original ion.

minCor

numeric(1) defining the minimal required correlation coefficient for MS2 chromatographic peaks to be considered for MS2 spectrum reconstruction.

intensity

character(1) defining the column in the chromPeaks matrix that should be used for the intensities of the reconstructed spectra's peaks.

peakId

optional character vector with peak IDs (i.e. rownames of chromPeaks) of MS1 peaks for which MS2 spectra should be reconstructed. By default they are reconstructed for all MS1 chromatographic peaks.

BPPARAM

parallel processing setup. See bpparam() for more information.

Details

In detail, the function performs for each MS1 chromatographic peak:

  • Identify all MS2 chromatographic peaks from the isolation window containing the m/z of the ion (i.e. the MS1 chromatographic peak) with approximately the same retention time than the MS1 peak (accepted rt shift can be specified with the diffRt parameter).

  • Correlate the peak shapes of the candidate MS2 chromatographic peaks with the peak shape of the MS1 peak retaining only MS2 chromatographic peaks for which the correlation is > minCor.

  • Reconstruct the MS2 spectrum using the m/z of all above selected MS2 chromatographic peaks and their intensity (either "maxo" or "into"). Each MS2 chromatographic peak selected for an MS1 peak will thus represent one mass peak in the reconstructed spectrum.

The resulting Spectra object provides also the peak IDs of the MS2 chromatographic peaks for each spectrum as well as their correlation value.

Value

Spectra() with the reconstructed MS2 spectra for all MS1 peaks in object. Contains empty Spectrum2 objects for MS1 peaks for which reconstruction was not possible (either no MS2 signal was recorded or the correlation of the MS2 chromatographic peaks with the MS1 chromatographic peak was below threshold minCor. Spectra metadata columns "ms2_peak_id" and "ms2_peak_cor" (of type CharacterList() and NumericList() with length equal to the number of peaks per reconstructed MS2 spectrum) providing the IDs and the correlation of the MS2 chromatographic peaks from which the MS2 spectrum was reconstructed.

Author(s)

Johannes Rainer, Micheal Witting

See Also

findChromPeaksIsolationWindow() for the function to perform MS2 peak detection in DIA isolation windows and for examples.


xiaodfeng/DynamicXCMS documentation built on Aug. 6, 2023, 3:02 p.m.