Abstract
Background
Comprehensive twodimensional gas chromatography timeofflight mass spectrometry (GCxGC/TOFMS) has been used for metabolite profiling in metabolomics. However, there is still much experimental variation to be controlled including both withinexperiment and betweenexperiment variation. For efficient analysis, an ideal peak alignment method to deal with such variations is in great need.
Results
Using experimental data of a mixture of metabolite standards, we demonstrated that our method has better performance than other existing method which is not modelbased. We then applied our method to the data generated from the plasma of a rat, which also demonstrates good performance of our model.
Conclusions
We developed a modelbased peak alignment method to process both homogeneous and heterogeneous experimental data. The unique feature of our method is the only modelbased peak alignment method coupled with metabolite identification in an unified framework. Through the comparison with other existing method, we demonstrated that our method has better performance. Data are available at http://stage.louisville.edu/faculty/x0zhan17/software/softwaredevelopment/mspa webcite. The R source codes are available at http://www.biostat.iupui.edu/~ChangyuShen/CodesPeakAlignment.zip webcite.
Trial Registration
2136949528613691
Background
Comprehensive twodimensional gas chromatography timeofflight mass spectrometry (GCxGC/TOFMS) has been employed in analysis of complex samples in metabolomics studies, especially in cancer [13]. However, experiment output still suffers from substantial variations, challenging the data interpretation in these studies. For this reason, the methodological development at the data analysis stage has become a crucial research area, which is still at its infancy. Throughout the paper, we mean experiment by technical run, which is performed by machine.
In practice, there are two types of experiment variations: variation within an experiment and variation between experiments. The need to control variation within an experiment has been reported by many researchers [46]. Theoretically, all instrument signals generated by one metabolite species should be reported as a single peak within an experiment. However, in reality, multiple peak entries can occur due to the abnormality in peak shape, high sensitivity of the peak detection algorithm, and experimental cause such as short modulation time [46]. To reduce such variation, peak merging should be done before proceeding to any subsequent analysis. On the other hand, variation between experimental runs is induced by factors such as difference in experiment configuration and runtorun variations. Betweenexperiment variation usually is much higher in magnitude than withinexperiment variation. Specifically, retention time (RT) depends heavily on experimental setup by the nature of experiment. For analysis purpose, different retention times of the same metabolite between experiments should be adjusted for further analysis. The process of such an alignment is usually referred to as peak alignment. Typically, a peak alignment process includes two steps: peak matching where the identity of peaks from each experiment is matched, and RT adjustment that is based on the results of the peak matching.
Several studies addressed the alignment issue of metabolomic profiling from GCxGC/TOFMS experiment [49]. From a methodological perspective, we can classify them into three categories. In the first generation methods [79], alignment implies RT adjustment, which is solely based on data of the retention time without the input of metabolite identification. For example, algorithms for aligning local region of interest were introduced in Fraga et al. and Mispelaar et al. [7,8]. And then, algorithm to align entire chromatogram in twodimensional GC was suggested by Pierce et al. [9]. However, the limitation of those methods is that aligning metabolites by using twodimensional retention times only may produce false alignment because some metabolites with similar chemical functional groups have similar retention times in both GC dimensions [5]. For this reason, the second generation methods have been developed exploiting two different types of information: closeness in twodimensional retention times and spectrum similarity [46]. Three wellknown methods of this generation are MSort [4], DISCO [5], and mSPA [6]. Since DISCO is a modified version of MSort, they are similar in many respects. Thus, our focus is on the difference between mSPA and the other two. First of all, the process of peak alignment in mSPA includes metabolite matching only, without reference to RT adjustment while the other two address both. Second, Kim et al. [6] defined and used a mixture similarity score, weighted average of RT distance and spectra similarity for peak matching while the other two used both information sequentially in each step. Third, as a spectrum similarity measure, Kim et al. [6] used dot product and the other two used Pearson's correlation coefficient. More comparison among these three methods can be found in [6].
Our method, as a third generation method, is unique in that it is modelbased approach. Compared to the second generation methods, our method is different in many respects. Compared to mSPA, our method considers rank distance of retention time instead of Euclidean distance which is used in mSPA. As a spectrum similarity measure, we use cosine score, angle between two spectra while mSPA uses dot product which is the cosine value of the angle. Also, our method covers both homogeneous and heterogeneous data while mSPA can handle homogeneous data only. For clarity, when we get data under the exactly same experiment configuration, we call it homogeneous. Otherwise, we call it heterogeneous. Most of all, our method uses posterior probability for metabolite matching based on an empirical Bayes model. The mSPA, however, defines an ad hoc likelihood function and maximises the function. Compared with DISCO, there are some aspects in common: both methods (1) can be applied to homogeneous and heterogeneous data, (2) address both peak matching and RT adjustment. On the other hand, they differ in four key ways: (1) our method does not need any RT transformation at the pairwise peak matching stage, (2) we use posterior probability as a matching confidence, (3) we use latticewise method for RT adjustment, not peakwise, (4) as a similarity measure, we use mixture similarity score with cosine score involved, but DISCO use Pearson's correlation coefficient.
Since DISCO can be applied to heterogeneous as well as homogeneous data, we compare our method with DISCO. In what follows, we provide a brief description of the model. Then we demonstrate the performance of our method with a mixture of standard compounds and a rat plasma data.
Results
Experiment datasets
Two different types of experiments are analyzed in this study: a mixture of standard compounds (Experiment I) and a rat plasma (Experiment II). We have three sets of homogeneous data from Experiment I corresponding to three different temperature gradients, respectively: dataset1 (5°C/min) with 10 replicates, dataset2 (7°C/min) with 2 replicates and dataset3 (10°C/min) with 4 replicates. To produce a heterogeneous dataset by using three homogeneous datasets available, we selected one technical replicate from each dataset and combined them, which is called dataset4. Thus, we have 10, 2, 4, 3 technical replicates in each dataset from Experiment I. From Experiment II, we have 5 homogeneous technical replicates but, no heterogeneous output. More details of experiments are given in Additional file 1.
Additional file 1. supplementary materials. This file include formula derivation and some results including tables and plots.
Format: PDF Size: 650KB Download file
This file can be viewed with: Adobe Acrobat Reader
Overview of algorithm
To help understanding of our results, we summarize our algorithm briefly because the order of results follows that of our algorithm. After peak merging, we select two experiment outputs and calculate matching confidence of peaks in the form of posterior probability through the empirical Bayes model. Based on these matching confidence (Equation 8 in Methods section), we select metabolite pairs with high matching confidence by applying cutoff value to the posterior probability. We then continue the same pairwise process for all other experiment outputs and generate representative landmark peaks. Given landmark peaks, we adjust RT of all peaks with respect to these peaks.
Peak merging
Peak merging is performed based on the result obtained by ChromaTOF software. In the case that multiple peaks exist, we select the peak with maximum peak area and remove the others [6]. The number of compounds before and after peak merging is summarized in Tables 1 and 2.
Table 1. Summary of Experiment I: number of compounds after/before peak merging
Table 2. Summary of Experiment II: number of compounds after/before peak merging
Landmark peaks
Here we choose threshold value (h = 40) which is used to calculate a_{j, }b_{j }and b_{j}* in layer 2 of our model (see Methods section). Also, we use weight (w = 0.1) for mixture score and apply cutoff value of 0.9 to matching confidence, posterior probability of correct match for landmark peak selection.
In Experiment I, 11, 40, 28 and 24 landmark peaks were selected for each dataset, respectively. In Experiment II, 31 landmark peaks were selected.
Peak alignment results
As an efficient way to illustrate alignment results in each dimension of RT (say marginal view), we consider kernel density estimate (KDE) along with normal kernel, which can be considered as a continuous version of histogram. Each KDE plot was made by using retention times of peaks to show the density of retention time. The brief introduction of KDE is given in Additional file 1 (see Section 14). KDE plots before/after RT adjustment for homogeneous data (dataset3) are given in Figure 1. Kernel density estimates corresponding to the first retention time before (left)/after (right) RT adjustment are provided in the top row and those corresponding to the second retention time are provided in the bottom row. Based on the KDE plots for homogeneous data, it is clear that after RT adjustment modes from 4 curves are well overlapped and are sharper. In other words, there are more densities around the mode and more distinct hills and valleys after adjustment, implying that peaks are well aligned.
Figure 1. Experiment I: KDE plots before/after peak alignment for homogeneous experiment (dataset3). Top: kernel density estimate (KDE) corresponding to the first retention time before (left)/after (right) peak alignment. Bottom: kernel density estimate of the second retention time before (left)/after (right) peak alignment.
Similarly, KDE plots for heterogeneous case (dataset4) before/after RT adjustment are given in Figure 2. As expected, there are more runtorun variations in heterogeneous data than homogeneous ones and after RT adjustment we can see much more dramatic change in heterogeneous case. It is clear that the different RT ranges (presented by different colors in Figure 2) are well adjusted after peak alignment.
Figure 2. Experiment I: KDE plots before/after peak alignment for heterogeneous experiment (dataset4). Top: kernel density estimate (KDE) corresponding to the first retention time before (left)/after (right) peak alignment. Bottom: kernel density estimate of the second retention time before (left)/after (right) peak alignment.
The four KDE plots for more complicated biological sample (Experiment II) are given in Figure 3. We see the same situation as seen in standard mixture data, i.e., sharper and distinct hills and valleys after alignment.
Figure 3. Experiment II: KDE plots before/after peak alignment for rat plasma data. Top: kernel density estimate (KDE) of the first retention time before (left)/after (right) peak alignment. Bottom: kernel density estimate of the second retention time before (left)/after (right) peak alignment.
For two dimensional view on alignment results (say joint view), scatter plots of RT after RT adjustment corresponding to KDE plots in Figures 2 and 3 are given in Figure 4. We can see after alignment that peaks are superimposed very well, implying that peaks with similar retention times are well aligned. More results for other datasets are provided in Additional file 1.
Figure 4. Two dimensional view on alignment results. Scatter plots after peak alignment for dataset4 (left) and rat plasma (right).
Performance comparison
For comparison, we select an existing peak alignment method which is able to align both homogeneous and heterogeneous data, DISCO [5]. As a comparison measure, we consider F1 score, with higher value implying better performance. F1 score is calculated under the assumption that ChromaTOF identification results are gold standard, harmonic mean of positive predictive value (PPV) and sensitivity [6,10]. For our method, we considered four different weights (w = 0.1, 0.2, 0.3, and 0.5) and 5 different cutoff values (0.6, 0.7, 0.8, 0.9, and 0.95) for posterior probability calculation. For each combination of w and cutoff values, we calculated F1 score for each paired technical runs within each dataset. For instance, given a pair of experimental outputs, we get 20 F1 scores, i.e., each score is calculated with each parameter combination, respectively. Among them, we selected the best F1 score. Then, average of such best F1 values obtained from all pairs was calculated. Once it is done, we repeat the same calculation for each dataset: dataset1,..., dataset4 and rat plasma data. Results corresponding to dataset1, dataset4 and rat plasma are summarized in Table 3. More results are provided in Additional file 1.
Table 3. Averaged best F1 score: our method v.s. DISCO
Based on the results, it is clear that the performance of our method is better than DISCO except homogeneous case from Experiment I. As seen in the results, our method has better performance for complex data. More precisely, the difference in performance is getting bigger as the complexity of the data increases.
We investigated the relationship between alignment results by gold standard and our method. For the purpose of comparison, we assumed that the identification by ChromaTOF is correct and then aligned the peaks based on their assigned names, resulting in the alignment by gold standard (GS). However, there might exist false positives in the aligned peak list of the gold standard if ChromaTOF assigned a wrong name to a compound. Actually, it is known that the accuracy of ChromaTOF identification is about 80%. To examine the concordance in the results of both methods, three sets of experimental pairs with the best F1 score (i.e., each pair come from each of three different datasets respectively) were selected. In Experiment I, R_{8}^{5 }and R_{9}^{5 }(w = 0.2 and cutoff = 0.8) were selected for homogeneous case (F1 = 0.93) and R_{1}^{5 }and R_{1}^{7 }(w = 0.1 and cutoff = 0.6) were selected for heterogeneous case (F1 = 0.86). Numbers in parenthesis in Figure 5 represent the number of compound pairs matched by each method in Experiment I. In the case of homogeneous data (left panel in Figure 5), our method found 71 peak pairs and 67 pairs of them had the same compound name identified by ChromaTOF. On the other hand, for heterogeneous case (right panel in Figure 5), our method found 68 matching pairs and 60 of them were verified to have same compound name by ChromaTOF.
Figure 5. Experiment I: Venn diagram corresponding to best F1 measure for homogeneous(left), heterogeneous(right) data. R_{8}^{5 }and R_{9}^{5 }(w = 0.2 and cutoff = 0.8) were selected for homogeneous case and R_{1}^{5 }and R_{1}^{7 }(w = 0.1 and cutoff = 0.6) were selected for heterogeneous case. The best F1 scores for each case are 0.93 and 0.86, respectively. GS stands for gold standard.
Similarly, Venn diagrams corresponding to the best F1 score for Experiment II (rat plasma) are given in Figure 6. D_{3 }and D_{4 }(w = 0.1 and cutoff = 0.8) were selected (F1 = 0.62). Our method found 337 peak pairs and 181 pairs of them had the same compound name identified by ChromaTOF. More results for other datasets are provided in Additional file 1.
Figure 6. Experiment II: Venn diagram corresponding to best F1 score for rat plasma experiment. D_{3 }and D_{4 }(w = 0.1 and cutoff = 0.8) were used. GS stands for gold standard.
Manual validation
To investigate the discordance in the results of both methods, we manually checked the alignment results for homogeneous data from Experiment I (left panel in Figure 5). Since 67 common peak pairs were aligned by both methods, our focus was on the other parts. The four peaks aligned by our method, which have different names given by ChromaTOF, were further analyzed based on raw image data. It is clear that one of them is correctly aligned. However, the other three might be incorrect. In addition to that, six peak pairs which had the same compound name by ChromaTOF, but not aligned by our method were also manually examined. Among them, ChromaTOF made wrong decision for three of them while the other three are correct.
To further investigate the discordant results, we selected two peaks (CAS 193395 and 629925) in target. They were aligned by both method but, aligned to different compounds in library (denoted by * in Additional file 1, Table S8). Those peak pairs aligned by GS were not the best peak pairs by our method, but the second best. For example, a compound with CAS(193395) in target peak list was aligned to the compound with CAS(191242) in reference peak list by our method, not to the compound with CAS(193395) because peak pair (193395 and 191242) had a similarity score of 9.67 while peak pair (193395 and 193395) had a similarity score of 9.69 (i.e., in our method, a pair with the smaller similarity score is selected). Similarly, even though the compound (CAS 629925) was assigned to different compound in library, the difference in score was not substantial. More details about the four and six peaks aligned by each method only are summarized in Additional file 1 (see Section 4.5).
Discussion
Compared to the existing peak alignment method DISCO that can also be applied to heterogeneous as well as homogeneous data, our method has many unique properties. First, our method is a modelbased approach. Second, unlike the DISCO, our method does not need any type of data transformation such as zscore transformation at any stage of the process. For example, in case of heterogeneous data, DISCO first transforms retention times using zscores in order to reduce the retention time variation among different GCxGCMS experiments. In other words, data transformation converts heterogeneous data into (pseudo) homogeneous data. Although a welldefined transformation will definitely improve the peak alignment (especially, for heterogeneous data) through the reduction of false positives caused by retention time variation, it is not easy to find an optimal transformation function that can handle all kinds of variation. Therefore, requirement of transformation is usually considered as a downside and it is avoided. Third, once representative landmark peaks are obtained, we make grid net using their retention times and adjust retention times of all peaks with respect to the grid net sequentially in both dimensions (all peaks simultaneously in each dimension).
With a mixture of metabolite standards and a rat plasma data, it is shown that our method has good performance in terms of F1 score. The F1 score requires gold standard (GS) and we used ChromaTOF identification results as a tentative gold standard. In the case of mixture of metabolite standards, even though more variations are involved in heterogeneous data, F1 score is always good regardless of data type. On the other hand, F1 score for complicated data from rat plasma is not that much high even though it is homogeneous. Compared to standard mixture data, although we see some decrease in F1 value for complicated data, the value from our method is still higher than that of the other existing method.
We compared our matching results with the tentative GS. Through the comparison, even though there is some discordance, we see high level of concordance in matching results by both methods, resulting in high F1 score for our method.
Conclusion
In this paper, we developed a modelbased peak alignment method to handle experiment variations, which can be applied to both homogeneous and heterogeneous experiments. Our method utilises a part of the output of ChromaTOF software as input data. The workflow of the method consists of two steps: pairwise peak matching and retention time adjustment. Due to the use of landmark peak lists composed of peak pairs with high matching confidence, our approach produces good quality of peak alignment.
In the peak alignment context, the excellent performance of our method at the data processing stage will have an enormous positive effect on subsequent analysis. For example, even though experiments are performed under the different experiment configurations or even at different times, the data aligned by our method can be used as input for further analysis: for example, time course metabolomic data analysis. Thus, the area to which our method can be applied might be extended to metabolite biomarker finding and metabolite clustering. Furthermore, as a future study, we we will study the relationship between peak alignment and peak identification to improve the accuracy of both preprocessing.
Methods
Our empirical Bayes model for pairwise peak matching between technical runs utilises the same structural hierarchy as the model constructed for metabolite identifications in [11]. Here we briefly review the model. Suppose that we have two experiment outputs and consider one of them as reference and the other as target in the context of peak alignment.
Model Review
We consider a hierarchical statistical model with four layers. All layers together address the process of our algorithm. Here is a brief overview of our hierarchical model: (i) we first check if a compound in library is in sample, (ii) depending on the information given in (i), we check if the compound is matched to any compound in sample, (iii) we then check if the match is correct because our matching using similarity score is not always correct, (iv) we finally estimate the distribution of similarity score.
Layer 1: We consider the marginal probability that each metabolite in the reference is present in target:
where N is the number of the peaks in the reference.
Layer 2: Given the Y_{j }information, we consider the conditional probability of metabolite j being matched to some target metabolite. According to the value of Y_{j}, two different conditional probabilities are considered: P[Z_{j }= 1Y_{j }= 0] and P[Z_{j }= 1Y_{j }= 1]. Note that even though a metabolite j does not exist in target (Y_{j }= 0), there is some chance for the metabolite to be claimed as present (P[Z_{j }= 1Y_{j }= 0] > 0). For the case Y_{j }= 0, we consider the following model:
Where and η_{0 }is an unknown constant. η_{0 }is a kind of auxiliary parameter which is not of our interest because it is related to metabolites with no matching neighbor. The b_{j }is defined using the metabolite in reference:
where is a mixture similarity score between peaks q and k in the reference, C is the set of peaks in the reference, and I(·) is the indicator function.
Similarly, we consider the following model for the case Y_{j }= 1:
Where and η_{1 }is an unknown constant, which is not of our interest. The b_{j}* is defined using the metabolite in reference:
where b_{j}* includes metabolite j itself as a neighbor to account for the fact that Y_{j }= 1.
Layer 3: For reference metabolites matched to at least one target metabolite, we consider conditional probability of W_{jl }given Y_{j }and Z_{j}, the correctness of those matches. For those matches of metabolite j with Y_{j }= 1 and Z_{j }= 1, we consider the following model:
Since is between 0 and 1, this implies that our matching is not always correct even though metabolite j is true positive.
Layer 4: we use a mixture model to characterize the distribution of the mixture similarity scores:
where f is the mixture of densities f_{T }and f_{F }that are the distributions of the scores of the correct matches and incorrect matches, respectively, and ϕ_{T }and ϕ_{F }are corresponding parameters.
Rationale behind the model: the rationale behind the use of a logistic function in layer 2 results from logistic regression. In other words, we investigated the relationship between Z and corresponding competition scores by logistic regression. Then, we found that quadratic function is statistically significant (see Section 13 in Additional file 1). Note that score function (f) consists of two score density functions: f = πf_{T }+ (1π)f_{F. }According to the distribution of observed scores, the specification of the score functions is decided. Here we assume normality. According to the distribution of observed scores, either f_{T }or f_{F }could be assumed a normal mixture. The parameter vector to be estimated is More details about model description can be found in [11].
Matching Confidence
The matching confidence of metabolite j in reference to a target metabolite can be calculated as the posterior probability of W_{jl}:
where is the estimated parameter vector. Note that this matching confidence plays a key role in the first step of peak alignment procedure.
Since we treat Y and W as the unobserved variables, we employ ExpectationMaximization (EM) algorithm to handle such latent variables [12]. More details about EM are provided in Additional file 1.
Mixture similarity score
Mixture similarity score is defined as:
where S(A, B) is mixture similarity score between two peaks A and B. Note that w is weight (0≤w≤1), D is rank distance based on retention time, and C is cosine score, angle between two peaks in high dimensional space (see Additional file 1 for details). Clearly, there is unit imbalance between RT distance and spectrum similarity. To balance them, we rescale rank distance (D) and angle (C) in Equation (9). Since considering rank of RT as a measure of closeness in RT reduces false positive rate [5], we take into account the elution order of RT in Equation (9). When calculating similarity score, we prefer small value of w(0 ≤ w ≤ 0.5) in order for spectrum similarity to play a important role in similarity score.
Peak alignment algorithm
As aforementioned, our algorithm for peak alignment consists of two steps: pairwise peak matching and retention time adjustment. The process of our algorithm is summarized in Figure 7.
Figure 7. Schematic representation of peak alignment algorithm. Five key steps of peak alignment algorithm.
Pairwise peak matching
Given two experiment outputs, we calculate peak matching confidence through empirical Bayes model. Once the same calculation is done for all pairwise outputs, we then select peaks connecting through all outputs, which are called landmark peaks. As an illustration, suppose that we have 3 experiment outputs: C1, C2, and C3. In Figure 8, "o" presents metabolites within each experimental output and a connection line presents metabolite pairs with matching confidence greater than prespecified cutoff value. The landmark peaks denoted by * are selected because those peaks exist in all outputs. We define representative landmark peaks as an average of those retention times (RT) across experiments.
Figure 8. An illustrative example of landmark peak selection. 3 experiment outputs: C1, C2, and C3. Six peaks denoted by * are selected as landmark peaks. o presents compound within each output and  presents matching pairs with matching confidence greater than prespecified cutoff value.
Retention time adjustment
Given representative landmark peaks, we align all peaks in each experiment output with respect to them sequentially in both dimensions (especially, latticewise, not peakwise). For example, a target peak to be aligned (t^{m }denoted by • in Figure 9) is moved to * after RT adjustment. Xaxis presents the first dimensional retention time and Yaxis presents the second dimensional retention time. The dotted grid presents target retention time (t_{L}^{1 }and t_{H}^{1}: the closest target RT1 below and above t^{m}, t_{L}^{2 }and t_{H}^{2}: the closest target RT2 below and above t^{m}). The solid grid presents reference retention time (r_{L}^{1 }and r_{H}^{1}: the corresponding reference RT1 aligned to target RT1, r_{L}^{2 }and r_{H}^{2}: the corresponding reference RT2 aligned to target RT2). The mathematical formula for the first dimensional retention time of new aligned peak (*) is given:
Figure 9. A graphical representation of retention time adjustment. A target peak (t^{m }denoted by ·) to be aligned is moved to new position denoted by (*) after peak alignment. Xaxis presents the first retention time and Yaxis presents the second retention time. The dotted grid presents target retention time (t_{L}^{1 }and t_{H}^{1}: the closest target RT1 below and above t^{m}, t_{L}^{2 }and t_{H}^{2}: the closet target RT2 below and above t^{m}). The solid grid presents reference retention time (r_{L}^{1 }and r_{H}^{1}: the corresponding reference RT1 aligned to target RT1, r_{L}^{2 }and r_{H}^{2}: the corresponding reference RT2 aligned to target RT2).
Where Similarly, new value for the second dimensional retention time can be obtained. In summary, retention time adjustment process consists of two sequential steps: the first dimension RT adjustment and the second dimension RT adjustment. RTs for all peaks in each dimension are aligned simultaneously with respect to the reference lattice, which is constructed based on representative landmark peaks. More details are provided in Additional file 1.
Authors' contributions
JJ and CS designed and formulated the statistical model. JJ developed the programs to implement the model. XZ and XS designed the two experiments. XS conducted the experiments. JJ and SK conceived the study and conducted comparison study. All authors read and approved the final manuscript.
Acknowledgement
This work was supported by the National Institutes of Health [1R01GM087735 to J.J., X.Z., and C.S.] and an Indiana University Showalter Research Trust Funding Award to J.J and C.S. and the Department of Defense grant [BC030400 to J.J. and C.S.] and the Department of Energy grant [DEEM0000197 to S.K.]
References

Kind T, Tolstikov V, Fiehn O, Weiss RH: A comprehensive urinary metabolomic approach for identifying kidney cancer.
Analytical Biochemistry 2007, 363:185195. PubMed Abstract  Publisher Full Text

Sreekumar A, Poisson LM, Rajendiran T, Khan AP, Cao Q, Yu J, Laxman B, Mehra R, Lonigro RJ, Li Y, Nyati MK, Ahsam A, KalyanaSundaram S, Han B, Cao X, Byun J, Omenn GS, Ghosh D, Pennathur S, Alexander DC, Berger A, Shester JR, Wei JT, Varambally S, Beecher C, Chinnaiyan AM: Metabolomic profiles delineate potential role for sarcosine in prostate cancer progression.
Nature 2009, 457:910915. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Hu JD, Tang HQ, Zhang Q, Fan J, Hong J, Gu JZ, Chen JL: Prediction of gastric cancer metastasis through urinary metabolomic investigation using GC/MS.
World Journal of Gastroenterololy 2011, 17:727734. Publisher Full Text

Oh C, Huang X, Regnier FE, Buck C, Zhang X: Comprehensive twodimensional gas chromatography/timeofflight mass spectrometry peak sorting algorithm.
Journal of Chromatography 2008, 1179:205215. PubMed Abstract  Publisher Full Text

Wang B, Fang A, Heim J, Bogdanov B, Pugh S, Libardoni M, Zhang X: ISCO: distance and spectrum correlation optimization alignment for twodimensional gas chromatography timeofflight mass spectrometrybased metabolomics.
Anal Chem 2010, 82:50695081. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Kim S, Fang A, Wang B, Jeong J, Zhang X: An optimal peak alignment for comprehensive twodimensional gas chromatography mass spectrometry using mixture similarity measure.
Bioinformatics 2011, 27:16601666. PubMed Abstract  Publisher Full Text

Fraga CG, Prazen BJ, Synovec RE: Objective data alignment and chemometric analysis of comprehensive twodimensional separations with runtorun peak shifting on both dimensions.

Mispelaar VG, Tas AC, Smilde AK, Schoenmakers PJ, van Asten AC: Quantitative analysis of target components by comprehensive twodimensional gas chromatography.
Journal of Chematography A 2003, 1019:1529. Publisher Full Text

Pierce KM, Wood LF, Wright BW, Synovec RE: A comprehensive twodimensional retention time alignment algorithm to enhance chemometric analysis of comprehensive twodimensional separation data.
Analytical Chemistry 2005, 77:77357743. PubMed Abstract  Publisher Full Text

Horai H, Arita M, Kanaya S, Nihei Y, Ikeda T, Suwa K, Ojima Y, Tanaka K, Tanaka S, Aoshima K, Oda Y, Kakazu Y, Kusano M, Tohge T, Matsuda F, Sawada Y, Hirai YM, Nakanishi H, Ikeda K, Akimoto N, Maoka T, Takahashi H, Ara T, Sakurai N, Suzuki H, Shibata D, Neumann S, Lida T, Tanaka K, Funatsu K, Matsuura F, Soga T, Taguchi R, Saito K, Nishioka T: MassBank: a public repository for sharing mass spectral data for life sciences.
J. Mass Spectrometry 2010, 45:703714. Publisher Full Text

Jeong J, Shi X, Zhang X, Kim S, Shen C: An empirical Bayes model using a competition score for metabolite identification in gas chromatography mass spectrometry.
BMC Bioinformatics 2011, 12:392. PubMed Abstract  BioMed Central Full Text  PubMed Central Full Text

Dempster AP, Laird NM, Rubin DB: Maximum likelihood from incomplete data via the EM algorithm.