Profile-Based LC-MS data alignment--a Bayesian approach.
A Bayesian alignment model (BAM) is proposed for alignment of liquid chromatography-mass spectrometry (LC-MS) data. BAM belongs to the category of profile-based approaches, which are composed of two major components: a prototype function and a set of mapping functions. Appropriate estimation of these functions is crucial for good alignment results. BAM uses Markov chain Monte Carlo (MCMC) methods to draw inference on the model parameters and improves on existing MCMC-based alignment methods through 1) the implementation of an efficient MCMC sampler and 2) an adaptive selection of knots. A block Metropolis-Hastings algorithm that mitigates the problem of the MCMC sampler getting stuck at local modes of the posterior distribution is used for the update of the mapping function coefficients. In addition, a stochastic search variable selection (SSVS) methodology is used to determine the number and positions of knots. We applied BAM to a simulated data set, an LC-MS proteomic data set, and two LC-MS metabolomic data sets, and compared its performance with the Bayesian hierarchical curve registration (BHCR) model, the dynamic time-warping (DTW) model, and the continuous profile model (CPM). The advantage of applying appropriate profile-based retention time correction prior to performing a feature-based approach is also demonstrated through the metabolomic data sets.
Department of Electrical and Computer Engineering,Virginia Tech, Washington, DC 20057, USA
Tsung-Heng
TH
Tsai
Mahlet G
MG
Tadesse
Yue
Y
Wang
Habtom W
HW
Ressom
D000465:Algorithms; D001499:Bayes Theorem; D001798:Blood Proteins / Q000032:analysis / Q000737:chemistry; D002853:Chromatography, Liquid / Q000379:methods*; D019295:Computational Biology / Q000379:methods*; D003198:Computer Simulation; D030562:Databases, Protein; D006801:Humans; D008390:Markov Chains; D013058:Mass Spectrometry / Q000379:methods*; D055442:Metabolome; D009010:Monte Carlo Method
D001798:Blood Proteins