mzrtsim: Raw Data Simulation for Reproducible Gas/Liquid Chromatography Mass Spectrometry Based Non-targeted Metabolomics Data Analysis

Abstract

Reproducibility of data analysis is pivotal in the context of non-targeted metabolomics based on mass spectrometry coupled with chromatography. While various algorithms have been proposed for feature or peak extraction, their validation often revolves around a limited set of known compounds or standards. While data simulation is widely used in other omics studies, simulation are focused on the feature level, neglecting uncertainties inherent in the feature or peak extraction process for metabolomics mass spectrometry data. In this technique note, we introduce a R package called mzrtsim to simulate gas/liquid chromatography full scan raw data in the mzML format. Unlike simulations solely based on virtual features, our approach leverages experimental spectral data from the MassBank of North America (MoNA) and the human metabolome database (HMDB). We developed algorithms to simulate chromatographic peaks, accounting for the tailing factor. The results of our study demonstrate the potential of this tool for comparing established metabolomics software (e.g., XCMS, mzMine, and OpenMS) against a ground truth. We found that the investigated software introduced false positive peaks and/or loss compounds with less peaks. They also showed different sensitivity to tailing and leading peaks and we also found setting appropriate intensity cutoffs will not lose information for the following statistical analysis. This R package is free and available online (https://github.com/yufree/mzrtsim). TOC O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=103 SRC="FIGDIR/small/567024v3_ufig1.gif" ALT="Figure 1"> View larger version (23K): org.highwire.dtl.DTLVardef@2861forg.highwire.dtl.DTLVardef@1b36989org.highwire.dtl.DTLVardef@1d5c939org.highwire.dtl.DTLVardef@8ddce6_HPS_FORMAT_FIGEXP M_FIG C_FIG