 Research
 Open Access
 Published:
A coherent mathematical characterization of isotope trace extraction, isotopic envelope extraction, and LCMS correspondence
BMC Bioinformatics volume 16, Article number: S1 (2015)
Abstract
Background
Liquid chromatographymass spectrometry is a popular technique for highthroughput protein, lipid, and metabolite comparative analysis. Such statistical comparison of millions of data points requires the generation of an interrun correspondence. Though many techniques for generating this correspondence exist, few if any, address certain wellknown runtorun LCMS behaviors such as elution order swaps, unbounded retention time swaps, missing data, and significant differences in abundance. Moreover, not all extant correspondence methods leverage the rich discriminating information offered by isotope envelope extraction informed by isotope trace extraction. To date, no attempt has been made to create a formal generalization of extant algorithms for these problems.
Results
By enumerating extant objective functions for these problems, we elucidate discrepancies between known LCMS data behavior and extant approaches. We propose novel objective functions that more closely model known LCMS behavior.
Conclusions
Through instantiating the proposed objective functions in the form of novel algorithms, practitioners can more accurately capture the known behavior of isotope traces, isotopic envelopes, and replicate LCMS data, ultimately providing for improved quantitative accuracy.
Background
Liquid chromatographymass spectrometry (LCMS) is a popular technique for elucidating the composition of liquid samples. Data processing considerations are essential to accurately determine the identity of molecules (analytes such as lipids or peptides) contained in the sample (a process called identification), as well as their quantity in sample (a process called quantification).
Information about sample quantity is captured directly in survey scans, or MS (aka MS1) data. Fragmentation spectra of one or more analytes constitute MS/MS (or MS2) data, and this information is typically used to corroborate or ascertain the identity of a molecule. Partitioning/clustering MS1 signal from complex samples and mapping the signal to other analyses (correspondence) is challenging. Some quantification strategies bypass these challenges by using information derived directly or indirectly from MS/MS data. These methods include spectral counting [1] and isobaric tags for relative and absolute quantification (iTRAQ) [2]. Though these methods have been successful, the amount of quantifiable signal embedded in MS1 data is estimated to far exceed what is currently available by MS/MS [3]; however, most MS1 data remains unused by current software. Hence, improving methods for partitioning and mapping MS1 signal stands to significantly (˜10 fold) increase the sensitivity of a typical labelfree or isotopelabeling MSomics experiment, both for experiments currently being run and for past experiments where raw data is still available.
Subdivision of raw mass spectrometer output data into smaller signal partitions attributed to specific analytes in the sample is critical prior to achieving analyte identification and quantification. The larger partition unit, called an isotopic envelope trace, is the signal pattern generated by each analyte/charge combination (see Figure 1). Because mass spectrometers can only detect charged analytes, the sample must be subjected to an ionization method, which imputes a charge on each detected analyte. Since multiple instances of each component exist in the sample, and since each instance is charged independently, there exist in each output the signals of multiple analytes, each with (potentially) multiple charge states. These create a distinct signalthe isotopic envelope tracefor the total signal detected for each analyte/charge state combination. Each isotopic envelope trace is composed of a series of isotope traces, which are manifestations of the fact that each analyte is composed of chemically similar compounds that differ in the weight of certain isotopes (such as ^{12}C vs ^{13}C). At each charge state, each molecular variant of the analyte is detected at a particular m/z offset, creating one isotope trace per molecular variant/chargestate/analyte combination.
Mass spectrometry data, in its raw form, is not ideal for isotope trace extraction or subsequent processing. After internally accumulating signal over discrete time slices, the mass spectrometer outputs raw data condensed into the form of many narrow profiles wherever signal is present. Conversion to centroid mode integrates the abundance of each of these profiles into a single tuple called a centroid. This is considered a routine conversion for which ample software is readily available. We adopt the typical convention of using centroid data.
Despite the ubiquity of LCMS experiments, to the best of our knowledge, no concise, complete description of the LCMS isotope trace and isotopic envelope extraction problems exists. Here, we describe constructs for isotope traces and isotopic envelopes, as well as formally describe the relationship of centroids, isotope traces and isotopic envelopes. In this context, we review extant objective functions for isotope trace extraction, isotopic envelope extraction, and correspondence. Finally, we propose novel objective functions for each of these tasks that address shortcomings in current approaches.
Results and discussion
Isotope trace extraction
The most important data processing step in a typical quantitative LCMS pipeline is isotope trace extraction [4]. Clustering centroids into isotope traces is a nontrivial problem due to the many sources of noise affecting centroid mass and abundance. Sources of noise affecting centroids include chemistry effects due to chromatography, abundance inaccuracy due to ionization efficiencies, m/z deviation due to machine calibration, occlusion/adulteration of lowabundance signal due to dynamic range limitations, and compounded inaccuracies in masstocharge ratio (m/z) and abundance due to centroid construction. Of course, these complications are propagated from the clustering of isotope traces to the clustering of isotopic envelopes to the identification of crossexperiment correspondence.
A centroid is denoted as c = (µ, τ, α) where µ, τ, α are values for m/z, retention time (RT), and abundance, respectively. A single MS run produces a set of centroids C={\left\{{c}_{i}\right\}}_{i=0}^{n}, where n can readily reach into the millions.
An isotope trace F ⊂ C is defined as a set of centroids: F={\left\{{c}_{i}\right\}}_{i=0}^{m}, with each set F constrained so that all members of a given isotope trace F are within a distance threshold θ from other centroids in their neighborhood ϒ (see Figure 2):
where θ is a function of centroid m/z, RT, and abundance, δ_{ F } is a distance function based on m/z, RT, and abundance, and ϒ is a neighborhood demarcated by m/z, RT, and abundance. Additionally, the slope of a (abundanceweighted) linear regressor estimate for an isotopic trace is very nearly infinite (in the m/z, RTplane). One way to formalize this is to use a weighted, inverse variant of the TheilSen estimator as follows (see Figure 3):
where c^{α} is the abundance of centroid c and c^{µ} is the m/z of centroid c.
Note that the behavior of isotope traces are dependent on all three MS dimensions although many common approaches to isotope trace extraction ignore one or more of these dimensions. For example, most proprietary MS software uses hard m/z bins for isotope trace extraction.
Extant objective functions
The prominent algorithms for isotope trace extraction include centWave [5], MatchedFilter [5], centroidPicker [6], massifquant [7], and MaxQuant [8].
MatchedFilter operates on the simplifying assumptions that 1) isotope traces are completely contained within preprocessed hard m/z bins and 2) the shapes of all isotope traces in a run can be fit to the same shape. MatchedFilter minimizes the error of a Gaussian fit over prospective isotope traces, by attempting to find the set of isotope traces \mathcal{F}, a scaling factor b_{ F } , and mean retention time F^{t} for each isotope trace that minimizes the summed abundance error over all isotope traces. Note the use of a single, global variance σ, an average RT width for all F ∈\mathcal{F}:
The centWave algorithm extracts isotope traces that fit a scaled and translated Ricker wavelet ζ (commonly called a Mexican hat function). The fit is calculated as a convolution between the shape function and the signal intensity (abundance), so the goal is to maximize the objective function:
where
with isotope tracespecific scaling parameter bF and translation parameter t_{ F } chosen to maximize the convolutional fit over isotope trace F .
The algorithm centroidPicker uses heuristic operations on a neighborhood graph to separated the data into connected components. It connects an undirected graph G = (C, N) of centroids where the edges N are constrained such that:
for some intensity threshold θ and centroid distance function δ_{ c }, resulting in G being composed of one or more connected components, each considered one isotope trace. Thus, \mathcal{F}=\left\{{F}_{i}\forall {c}_{k}\in {F}_{i},{\exists}_{cl\in {F}_{i}}\left\{{c}_{l}\in \Upsilon \left({c}_{k}\right)\right\}\right\}, where the neighborhood function ϒ (c) returns the set of nodes connected to c (and is symmetric because G is undirected).
The objective functions for massifquant and MaxQuant define \mathcal{F} as the set of all F formed by iterating over values of time t, and adding c if c^{τ} = t and \left{c}^{\mu}{c}_{*}^{\mu}\right<\in, where c_{∗} ∈ F and {c}^{\tau}{c}_{*}^{\tau}\le {c}^{\tau}{c}_{j}^{\tau} for all c_{ j } ∈ F. For massifquant, ∈ is prescribed by a Kalman filter induced from the variance in c^{µ} and c^{α} for all c_{ j } ∈ F such that {c}_{j}^{\tau}<t, with the added constraint that c^{τ} be unique in F . MaxQuant defines ∈ simply as a distance threshold of 7 ppm m/z.
Proposed objective functions
We define F^{µ}, the m/z of isotope trace F, given by the weighted m/z of its component centroids:
and using it propose an alternative objective function for isotope trace extraction:
where, again, centroid clustering \mathcal{F} and retention time means F^{t} are chosen to minimize the Gaussian fit error; however, rather than using a single global variance in the RT dimension, each isotope trace F has a local variance σ_{ F }; in addition, the scaling factors have become timedependent scalar functions b_{ F }(·). The second Gaussian factor, parameterized by mean F^{µ} and variance function h(·), models the m/z width of the isotope trace, which is a function of the abundance α. Isotope traces splay at low abundance and narrow at high abundance; thus, both the variance h(·) and the scaling factors a_{ F }(·) are modeled as functions dependent on the abundance α. Note that while variance is traceindependent (depending only on abundance), each isotope trace has its own scaling function (which in turn is dependent on abundance).
Alleviating current limitations in isotopic trace extraction
Current objective functions for isotopic trace extraction fail to capture isotopic trace behavior formalized in this section: namely, a pattern of centroids forming a generally tight distribution through time around a specific m/z, with variation occurring as a factor of abundance, with normal abundance traces splaying at the beginning and end of elution, and lower abundance traces displaying high m/z variance in general. Moreover, isotope traces are skewed in time, with sharp onset of intensity followed by a postpeak long tail. The shape of traces is almost never strictly Gaussian (or even symmetric), as chromatography almost always deviates from the Gaussian in heading (which is more steep) and in tailing (which is less steep). Our objective functions account for each of these behaviors.
Isotopic envelope extraction
The LCMS clustering problem is defined as a twostep partitioning problem. In the first step, isotope trace extraction, we require a partition ϕ of the set of all centroids C into the set of isotope traces \mathcal{F}, \varphi \left(C\right)={\left\{{F}_{i}\right\}}_{i=1}^{r}=\mathcal{F}with the properties:
In other words, 1) all centroids are assigned to an isotope trace; 2) isotope traces can't share centroids. Because any sensor's detection of a physical system will deviate somewhat from the true physical system, we can expect MS detections to contain extraneous centroids. However, all signal ought to be accounted for (even if some identified "traces" eventually are identified as noise) and, in a platonic model, ought to be assigned to an isotope trace.
In the second step, isotopic envelope extraction, we require a partition ψ of the set of isotope traces \mathcal{F} into the set of isotopic envelopes \epsilon ,\phantom{\rule{2.77695pt}{0ex}}\psi \left(\mathcal{F}\right)={\left\{{E}_{i}\right\}}_{i=1}^{p}=\epsilonwith the property
The choice of partitions φ and ψ is guided by a set of distance functions Δ that define distances between centroids, isotope traces, isotopic envelopes, etc. and objective functions λ_{ F } and λ_{ E } that describe "good" isotope traces and isotopic envelopes, respectively. The choice of distance and objective functions, along with choice of optimization procedure, characterizes an algorithmic approach for solving this clustering problem. A defining general property of isotopic envelopes, however, is the regular spacing between component isotope traces. In addition, for virtually all molecules from biological sources we expect that if there is an isotope with index j and an isotope with index j + 2, then there exists an isotope with index j + 1.
An isotopic envelope E is the set of isotope traces F_{ i } that are produced by a given analyte/charge state combination: E={\left\{{F}_{i}\right\}}_{i=0}^{q} subject to the constraint that the m/z difference between each consecutive (assuming an ordering of centroids from least mass to greatest mass) isotope trace in E must be equivalent to \frac{k}{{z}_{E}}+\in, where k is the mass of a neutron, z_{ E } is the integer charge of E and ∈ is a noise tolerance parameter. That is, assuming an indexing function {\iota}^{\mathit{\mu}}:\epsilon \times \mathcal{N}\to \mathcal{F} that returns the ith least massive isotope trace in an isotopic envelope:
The m/z m of the jth isotope trace in E must be roughly equivalent to
where \stackrel{\u0303}{m} is the uncharged molecular weight of the ion.
Every isotope trace consists of signal from at least one isotopic envelope, and, in the case of overlapping isotopic envelopes, an isotope trace may be composed of signal from more than one isotopic envelope.
Extant objective functions
FeatureFinder [9] is an isotopic envelope extraction algorithm in OpenMS that searches directly for E. Although the details are not completely clear, it appears that the algorithm attempts to minimize
where the G_{ E } compute a comparison between the (µ, τ, α) values for a centroid and the expected centroid values obtained from a heuristic isotopic envelope shape. Note that isotopic trace extraction is ignored.
MSInspect [10], another approach to isotopic envelope extraction, groups all coeluting signals and compares them to a simulated envelope calculated from a Poisson distribution parameterized by m/z, with the goal being to minimize the KL divergence between the Poisson distribution and the "distribution" of abundance in an instantaneous profile of the envelope at time τ :
where the notation c ∈^{τ} F means that c ∈ F at time τ, E is the maximal intensity (instantaneous) isotopic envelope (at time τ), \widehat{P}\left(\cdot \right) is the ratio of the intensity of isotope trace F (at time τ) to the total intensity of all isotope traces F ∈ E (at time τ), and P_{ m }(·) is the value of the Poisson distribution at c^{µ}.
Proposed objective functions
We propose an alternative objective function for isotopic envelope extraction:
where β is a relative importance weighting coefficient. The first term computes the deviation of member isotope traces from the expected chargebased m/z intervalwe want the isotope traces in envelope E to fit expected m/z spacing:
The second term computes the deviation in elution time of member isotope traceswe want all the isotope traces in isotopic envelope E to coelute within a small time window:
where F^{τ} could be defined analogously to Equation 7, could be the maximum intensity for isotopic trace F or could be some other reasonable definition for isotopic trace elution time.
We want to optimize ε and the z_{ E } so that λ_{ E } is minimized; that is, we want to find chargestate/isotopicenvelope pairs such that the errors in expected m/z and coelution time are minimized.
The isotopic envelope extraction segment of the MaxQuant [8] algorithm is one of the possible instantiations of this objective function, though many possibilities exist for how to set the allowable m/z and RT error and how to generate the prerequisite list of isotope traces.
Alleviating current limitations in isotopic envelope extraction
Isotopic envelopes are rich with data: the expectation of contiguous isotope traces with a uniform m/z charge gap, and similar maximal abundance across all isotope traces. Accounting for this behavior is not possible without adopting an isotope tracecentric approach to data extraction. Reliance upon maximal elution time alonean approach that is susceptible to conflation with overlapping envelopes in complex samplesis not a sensitive approach in envelopes of lower abundance, where maximal elution times are not pronounced. Moreover, by first finding the isotope traces, the exact m/z of each isotope trace can be calculated using a weighted average, alleviating the need for larger than theoretically justified isotope trace gaps, which will not be sensitive in complex samples with overlapping isotopic envelopes. Instead, the proposed objective functions leverage a precise and reliable m/z charge gap and adjacency of isotope traces along with maximal elution times, using all the information in the data.
Correspondence
The final objective of almost every MS experiment is the differential analysis of more than one MS run. This comparison allows the identification of significant quantity and component differences, useful for applications such as drug design, disease treatment, biological processes research and chemical forensics. Correspondence yields a mapping between isotopic envelopes in different runs (see Figure 4), a prerequisite for differential analysis.
The combination of noise from within one run (enumerated above) and noise from run to runmost notable in retention time shifts, where an isotopic envelope appears at a different retention time or with a compressed or stretched RT length compared to another runmake LCMS correspondence nontrivial.
The correspondence mapping should again optimize an objective function which, in turn, characterizes an algorithm choice for solving the correspondence problem.
Extant objective functions
According to a recent review on LCMS correspondence algorithms [11], all extant approaches use either centroid data or a reduction of isotopic envelope traces into a single centroid. Of the almost sixty algorithms reviewed there, nearly all use the same objective functionfinding a family of onetoone partial functions χ_{ r } : ε_{ r } → ε_{ ∗ } (a different function for each experimental run r), where ε_{ ∗ } is the set of envelopes from a reference run, that minimizes global RT and m/z distance between isotopic envelopes (in any of their reduced forms, according to the authors):
where δ()^{τ,µ} is a distance function defined over RT and m/z.
The continuous profile model (CPM) [12] uses a different objective function, and thus is free from the reference requirement that most other algorithms have, allowing for a symmetric solution (one that is not dependent on the choice of a reference run). Additionally, the mapping is somewhat more localized than that of most correspondence algorithms. CPM minimizes the log likelihood of differences between a hidden Markov model mτ of the RT of a latent run and observed runs:
where D is the set of observed runs.
Proposed objective functions
In contrast to existing LCMS correspondence objective functions, the objective functions suggested here use the entire isotopic envelope. This allows greater discrimination by using isotope trace quantity and spacing to match isotopic envelopes from different runs. This extra discrimination is essential given the amount of RT variance and (to a lesser degree) m/z variance present in the data.
Let R be a set of runs, each of which has an associated set of isotopic envelopes {\epsilon}_{r}={\left\{{E}_{i}^{r}\right\}}_{i=1}^{pr},1\le r\le \leftR\right and let \stackrel{\u0303}{\epsilon}={\cup}_{r}{\epsilon}_{r}. We seek to find a binary equivalence relation ρ that induces a set of correspondence classes over \stackrel{\u0303}{\epsilon} that is reflexive (an envelope corresponds with itself), symmetric (if envelope E_{1} from run 1 corresponds with envelop E_{2} from run 2, then E_{2} also corresponds with E_{1}) and transitive (if envelope E_{1} from run 1 corresponds with envelope E_{2} from run 2 and envelope E_{2} corresponds with envelope E_{3} from run 3, then E_{1} corresponds with E_{3}); and if \rho \left({E}_{i}^{r},{E}_{j}^{s}\right)=\text{TRUE}, then for k ≠ i, \rho \left({E}_{k}^{r},{E}_{j}^{s}\right)=\text{FALSE} and for k ≠ j, \rho \left({E}_{i}^{r},{E}_{k}^{s}\right)=\text{FALSE} (an envelope from one run may have 0 or 1 matches from any other run; note that due to reflexivity, this also means that two nonidentical envelopes from the same run never correspond).
This relation should minimize

The difference in charge state between corresponding isotopic envelopes, {\delta}_{charge}.

The difference in m/z between isotope traces in corresponding isotopic envelopes, {\delta}_{m{z}_{it}}.

The difference in elution duration between isotope traces in corresponding isotopic envelopes, {\delta}_{dur}.

The difference in isotope abundance ratios between corresponding isotopic envelopes, {\delta}_{ratio}.

The difference in m/z between corresponding isotopic envelopes, {\delta}_{m{z}_{ie}}.

The number of singleton correspondence classes, {\delta}_{orphan}.

The difference in retention time between corresponding isotopic envelopes, {\delta}_{rt}.
An objective function incorporating all of these variables can take many forms, with perhaps the simplest generalization being a weighted linear combination, with weighting coefficients ω allowing relative prioritization:
with the summation over ρ(E_{1}, E_{2}) meaning a summation taken over all pairs of envelopes E_{1}, {E}_{2}\in \stackrel{\u0303}{\epsilon} for which ρ(E_{1}, E_{2}) = TRUE. Given the weighting coefficients ω, the most desirable correspondence would be that induced by the relation ρ* that minimizes λ_{ corr } (see Figure 4),
Alleviating current limitations in correspondence
Recently, several ubiquitous shortcomings were identified in a review of over 50 LCMS correspondence algorithms [11]. The most significant of these shortcomings was the fact that all current LCMS correspondence algorithms make model assumptions that fail to capture common behavior. In other words, each algorithm is constructed in such a way that the algorithm is guaranteed to get the wrong answer under certain conditions that are common to real LCMS data. The behaviors discussed included the ideas that:

Not all analytes appear in all replicates.

Elution order can swap.

Shifts occur in m/z as well as in RT.
Some correspondence methods reduce isotopic envelopes to a single point representation. This deprives the method of a rich source of distinguishing data found in full isotopic envelopesthe expectation of contiguous isotope traces with a uniform m/z charge gap, number of isotope traces, and relative abundance ratio of isotope traces. Similarly, most correspondence algorithms conduct an initial RT alignment, where signals (almost always muchreduced from the full isotopic envelope, and rarely built up from isotope traces to isotopic envelopes) are shifted up or down in RT (preserving original order) in order to most closely match a reference run. This is invariably followed by direct matching. The problem is that the initial warping is a lossy procedure that adulterates the original RT time, which would be useful to probabilistically ascertaining the closest corresponding isotopic envelope.
The proposed objective function does not force matches between runs, as it is very common for species to either not be present or fall below the signaltonoise ratio in differential studies. Instead, the proposed objective function leverages the full breadth of isotope envelope information, allowing a rigorous direct comparison of candidate correspondences based on all available data to select the most likely correspondence (in the sense of minimizing error), or no correspondence at all if that is the most likely case given the data.
Conclusions
We present a concise attempt to formalize LCMS data clustering problems, describing the constructs of isotope traces and isotopic envelopes and their relational structure. We provide a review of current approaches to isotope trace extraction and LCMS correspondence, and propose novel objective functions for both tasks that address shortcomings in current methods.
References
 1.
Choi H, Fermin D, Nesvizhskii AI: Significance analysis of spectral count data in labelfree shotgun proteomics. Mol Cell Proteomics. 2008, 7 (12): 23732385. 10.1074/mcp.M800203MCP200.
 2.
Wiese S, Reidegeld KA, Meyer HE, Warscheid B: Protein labeling by iTRAQ: a new tool for quantitative mass spectrometry in proteome research. Proteomics. 2007, 7 (3): 340350. 10.1002/pmic.200600422.
 3.
Michalski A, Cox J, Mann M: More than 100,000 Detectable Peptide Species Elute in Single Shotgun Proteomics Runs but the Marjority is Inaccessible to DataDependent LCMS/MS. Journal of Proteome Research. 2011, 10: 17851793. 10.1021/pr101060v.
 4.
Cappadona S, Baker PR, Cutillas PR, Heck AJ, van Breukelen B: Current challenges in software solutions for mass spectrometrybased quantitative proteomics. Amino Acids. 2012, 43 (3): 10871108. 10.1007/s0072601212898.
 5.
Tautenhahn R, Bottcher C, Neumann S: Highly sensitive feature detection for high resolution LC/MS. BMC Bioinformatics. 2008, 9 (1): 50410.1186/147121059504.
 6.
Pluskal T, Castillo S, VillarBriones A, Oresic M: MZmine 2: Modular framework for processing, visualizing, and analyzing mass spectrometrybased molecular profile data. BMC Bioinformatics. 2010, 11 (1): 39510.1186/1471210511395.
 7.
Conley CJ, Smith R, Torgrip RJ, Taylor RM, Tautenhahn R, Prince JT: Massifquant: opensource Kalman filter based XCMS isotope trace feature detection. Bioinformatics. 2014, 359
 8.
Cox J, Mann M: MaxQuant enables high peptide identification rates, individualized ppbrange mass accuracies and proteomewide protein quantification. Nature Biotechnology. 2008, 26 (12): 13671372. 10.1038/nbt.1511.
 9.
Weisser H, Nahnsen S, Grossmann J, Nilse L, Quandt A, Brauer H, Sturm M, Kenar E, Kohlbacher O, Aebersold R, et al: An automated pipeline for highthroughput labelfree quantitative proteomics. Journal of Proteome Research. 2013
 10.
Bellew M, Coram M, Fitzgibbon M, Igra M, Randolph T, Wang P, May D, Eng J, Fang R, Lin C, et al: A suite of algorithms for the comprehensive analysis of complex protein mixtures using highresolution LCMS. Bioinformatics. 2006, 22 (15): 19021909. 10.1093/bioinformatics/btl276.
 11.
Smith R, Ventura D, Prince JT: LCMS Alignment in Theory and Practice: A Comprehensive Algorithmic Review. Briefings in Bioinformatics. 2013
 12.
Listgarten J, Neal RM, Roweis ST, Wong P, Emili A: Difference detection in LCMS data for protein biomarker discovery. Bioinformatics. 2007, 23 (2): 198204. 10.1093/bioinformatics/btl553.
Author information
Affiliations
Corresponding author
Additional information
Competing interests and declarations
The authors declare that they have no competing interests. The publication costs for this article were funded by the University of Montana Office of Research and Sponsored Programs.
This article has been published as part of BMC Bioinformatics Volume 16 Supplement 7, 2015: Selected articles from The 11th Annual Biotechnology and Bioinformatics Symposium (BIOT2014): Bioinformatics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/16/S7.
Authors' contributions
RS, JTP and DV all contributed in writing this manuscript.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Smith, R., Prince, J.T. & Ventura, D. A coherent mathematical characterization of isotope trace extraction, isotopic envelope extraction, and LCMS correspondence. BMC Bioinformatics 16, S1 (2015). https://doi.org/10.1186/1471210516S7S1
Published:
DOI: https://doi.org/10.1186/1471210516S7S1
Keywords
 Proteomics
 Correspondence
 Alignment
 Lipidomics