- Methodology article
- Open access
- Published:

# Correcting for the effects of natural abundance in stable isotope resolved metabolomics experiments involving ultra-high resolution mass spectrometry

*BMC Bioinformatics*
**volumeÂ 11**, ArticleÂ number:Â 139 (2010)

## Abstract

### Background

Stable isotope tracing with ultra-high resolution Fourier transform-ion cyclotron resonance-mass spectrometry (FT-ICR-MS) can provide simultaneous determination of hundreds to thousands of metabolite isotopologue species without the need for chromatographic separation. Therefore, this experimental metabolomics methodology may allow the tracing of metabolic pathways starting from stable-isotope-enriched precursors, which can improve our mechanistic understanding of cellular metabolism. However, contributions to the observed intensities arising from the stable isotope's natural abundance must be subtracted (deisotoped) from the raw isotopologue peaks before interpretation. Previously posed deisotoping problems are sidestepped due to the isotopic resolution and identification of individual isotopologue peaks. This peak resolution and identification come from the very high mass resolution and accuracy of FT-ICR-MS and present an analytically solvable deisotoping problem, even in the context of stable-isotope enrichment.

### Results

We present both a computationally feasible analytical solution and an algorithm to this newly posed deisotoping problem, which both work with any amount of ^{13}C or ^{15}N stable-isotope enrichment. We demonstrate this algorithm and correct for the effects of ^{13}C natural abundance on a set of raw isotopologue intensities for a specific phosphatidylcholine lipid metabolite derived from a ^{13}C-tracing experiment.

### Conclusions

Correction for the effects of ^{13}C natural abundance on a set of raw isotopologue intensities is computationally feasible when the raw isotopologues are isotopically resolved and identified. Such correction makes qualitative interpretation of stable isotope tracing easier and is required before attempting a more rigorous quantitative interpretation of the isotopologue data. The presented implementation is very robust with increasing metabolite size. Error analysis of the algorithm will be straightforward due to low relative error from the implementation itself. Furthermore, the algorithm may serve as an independent quality control measure for a set of observed isotopologue intensities.

## Background

Application of mass spectrometry to stable isotope tracing experiments for the elucidation of glucose dates back to at least the early 1980's [1, 2]. The general scheme for these experiments is to supply a labeled precursor such as uniformly-labeled ^{13}C glucose ([U-^{13}C]-glucose) to a bacterial culture, tissue culture, or a whole multicellular organism and then extract a set of cellular or excreted metabolites for analysis [3, 4]. For identified metabolites, specific patterns of isotopologues are usually observed, which are then interpreted within the context of known cellular metabolic pathways [3â€“5]. Recently, we applied this technique to elucidate specific aspects of lipid metabolism [6].

The ultra-high resolution capability of Fourier transform-ion cyclotron resonance-mass spectrometry (FT-ICR-MS) makes it possibility to identify simultaneously hundreds, if not thousands, of metabolites from crude cell extracts without the need for chromatographic separation [6]. The better than 1 ppm mass accuracy of state-of-the-art FT-ICR-MS is often high enough to provide mass-to-charge ratios (m/z) down to the 3^{rd} and 4^{th} decimal place for metabolites less than a few thousand Daltons. This is accurate enough to distinguish relativistic mass differences between expected isotopes of CHONPS elements and unambiguously determine the isotope-specific molecular formula of an individual peak. Furthermore, the FT-ICR-MS's high mass resolution allows for the direct detection or deconvolution of individual isotopologues or mass-equivalent sets of isotopomers for a given metabolite.

Isotopologue identification and quantification of thousands of metabolites in these metabolomic experiments can provide a wealth of data for modeling the flux through metabolic networks. But before isotopologue intensity data can be properly interpreted, the contributions from isotopic natural abundance must be factored out (deisotoped). This is a computationally expensive and analytically intractable problem for data from lower mass resolution spectrometers where individual isotopically-resolved isotopologues cannot be distinguished [7]. In these instances, numerical methods have been employed to approximate and subtract the contributions from isotopic natural abundance [4, 7â€“9]. Some of these calculations are aimed at a different deisotoping problem, namely identifying the related isotopologues and calculating the monoisotopic mass from its isotopic mass distribution [10, 11]. Fortuitously, with the isotope-resolved isotopologue peaks from FT-ICR-MS histograms, we can pose a similar but distinct problem that allows for the derivation of a computationally tractable analytical solution. In addition, isotopologues derived from the same molecule (or very similar set of molecules) neatly handle peak intensity referencing issues by providing a natural internal reference.

## Results

### Derivation of the analytical solution

Equation 1 represents the relative distribution of carbon isotopologues from natural abundance only, as a sum of multinomial coefficients multiplied by the intensity of I_{M+0}, the theoretically untainted ^{12}C monoisotopic peak. The terms being summed are similar in form to those presented in Snider, 2007. I_{M+i;NA} is the expected intensity of the i^{th} isotopologue peak representing i additional nucleons. NAx_{C} is the fractional natural abundance of the ^{X}C isotope. C_{Max} is the number of carbons in the molecule. The multinomial coefficients, derived from the multinomial theorem with 3 variables represent the number of possible isotopomers of identical mass for a molecule with C_{Max} carbons given 3 isotopes of carbon: ^{12}C, ^{13}C, and ^{14}C.

Isotopologue peaks containing ^{14}C are typically not observed, since the isotope is very rare. Moreover, due to the very high mass resolution in FT-ICR-MS histograms, isotopologue peaks representing molecules comprised exclusively of the major isotope of CHONPS elements (expected elements for biological systems) along with ^{13}C, are completely resolved/deconvoluted and identified. Thus, we can ignore the contributions from ^{14}C and from minor isotopes of all other elements excluding carbon. This simplifies the calculation to a *single term with a binomial coefficient* (binomial term) shown in Equation 2, where NA13_{C} â‰ˆ 0.01109. The binomial coefficient represents the number of possible isotopomers of identical mass for a molecule with C_{Max} carbons given only 2 isotopes of carbon: ^{12}C and ^{13}C.

At natural abundance, each peak, I_{M+i;NA}, is directly related to the theoretically untainted ^{12}C monoisotopic peak, I_{M+0}, that has a fractional intensity of 1 when dividing by the sum of isotopologue intensities. However, once ^{13}C is incorporated into the molecule from a labeling source, the calculation of the contributions from natural abundance becomes more complex [8, 9]. The effects of ^{13}C natural abundance now depend on the amount of ^{13}C label already present. With each ^{12}C/^{13}C isotopologue resolved in the FT-ICR-MS histogram, we can use a series of binomial terms to accurately describe and correct for ^{13}C natural abundance. Equation 3 shows the basic form of these binomial terms as B_{C}(n, k) where k represents the total number of ^{13}C carbons present, n represents the number of ^{13}C carbons due to incorporation from a labeling source, and k-n is the number of ^{13}C carbons due to natural abundance. The binomial coefficient in Equation 3 enumerates the number of ways that k-n ^{13}C carbons can be incorporated into the molecule when n carbons are already labeled with ^{13}C. Equation 4 shows the first series needed in the correction, B_{C}sum(n) which represents the fraction of I_{M+i} intensity that is converted to other isotopologues due to the effects of natural abundance.

Equation 5 shows the full correction as the original isotopologue intensity minus natural abundance contributions based on lower mass untainted isotopologue intensities. Division by the fractional intensity, 1 - B_{C}sum(i), compensates for natural abundance effects that lower the intensity of the given isotopologue. As illustrated in Table 1, Equation 5 must be applied in a sequential fashion starting with i = 0, since the results of each step are needed in subsequent steps. In other words, the natural abundance corrected intensities of isotopologues with lower ^{13}C incorporation from labeling are needed to calculate the natural abundance correction of isotopologues with higher ^{13}C incorporation from labeling.

Since ^{15}N incorporation can be distinguished from ^{13}C incorporation due to the very high mass resolution in FT-ICR-MS histograms, it takes only a trivial conversion of Equations 3, 4, and 5 to handle labeling in ^{14}N/^{15}N isotopologues. We simply replace N_{Max} for C_{Max} and NA15_{N} for NA13_{C}. However, handling all of the mixed ^{14}N/^{15}N/^{12}C/^{13}C isotopologues that arise from simultaneous ^{13}C and ^{15}N labeling requires a series of two binomial terms multiplied together as shown in Equations 6 and 7. Given the peaks are isotopically resolved, there are CMax * NMax separate observable isotopologues, whose intensities are represented by I_{M+i, j;NA}. The multiplied binomial terms, B_{C}(x, i) * B_{N}(y, j), describe the combined effects from both carbon and nitrogen natural abundance.

A version of each equation in larger fonts is available in Additional file 1.

### Implementation of the algorithm

We implemented Equations 3, 4, and 5 as an iterative algorithm in the Perl programming language [Additional file 2]. Iteration allows the algorithm to partially compensate for missing (zero intensity) isotopologues. The algorithm (Figure 1) starts with C_{Max} and the observed ^{12}C/^{13}C isotopologue intensities contaminated by contributions from ^{13}C natural abundance. Based on C_{Max}, the algorithm precalculates the binomial coefficients needed in later steps using Equations 3 and 4. During each iteration, the algorithm performs three steps. In step 1, the algorithm calculates the set of uncontaminated ^{12}C/^{13}C isotopologue intensities using Equation 5 and the observed intensities supplemented with calculated contaminated intensities for missing isotopologues. From Equation 5, it is apparent that this must be done in ascending mass order starting with I_{M+0}. Sometimes, small negative uncontaminated intensities arise from errors in the observed intensities. These negative intensities are flattened to zero, since they have no basis in reality. Next, the algorithm renormalizes the uncontaminated intensities based on the sum of observed intensities. This is required since missing isotopologues were supplemented with calculated values and negative intensities are flattened to zero. In step 2, the algorithm calculates the set of contaminated intensities based on the uncontaminated set by solving for I_{M+i;NA} in Equation 5. In step 3, the algorithm calculates the absolute difference between observed and calculated contaminated intensities. If this difference decreases, the algorithm performs another iteration until no more improvement is seen. Finally, the algorithm prints the results and ends.

### Testing the implementation

We created several sets of simulated isotopologues (test sets) with varying levels of ^{13}C-labeling and added the expected contributions (contamination) from ^{13}C natural abundance by solving for I_{M+i;NA} in Equation 5. We then tested the implementation with these test sets. Figure 2 shows the results for three of these test sets of a hypothetical metabolite with 20 carbon atoms. The ^{13}C natural abundance contaminated intensities are in red and corrected intensities in green. The red bars in Figure 2A represents the expected observed isotopologue intensities when no ^{13}C-labeling is present. This naturally collapses into a single green ^{12}C monoisotopic peak with correction. Figure 2B shows the contaminated and corrected isotopologue intensities when equal amounts of ^{13}C-labeling for 8, 10, and 12 carbons are present. There is a tapering phenomenon observed in the contaminated intensities due to the fact that the number of carbons affecting the intensities decreases with increasing amounts of ^{13}C-labeling. Equation 3 captures this phenomenon within its binomial coefficient where it is further demonstrated in Figure 2C with natural abundance having no effect on a metabolite with 100% ^{13}C-labeling.

The implementation is also quite efficient even in an interpreted programming language like Perl. 10,000 repetitions of this algorithm for all 3 simulated test sets took only 17 seconds on a single core of an Intel T7200 Core 2 Duo mobile processor with 2GB of RAM and running release 5.3 of the RedHat Enterprise Linux operating system. The implementation is also very accurate. Given the perfect data in these three simulated test sets, the largest error was 4.12 Ã— 10^{-16} seen in the I_{M+1} corrected intensity for the test set representing no ^{13}C-labeling (Figure 2A). Furthermore, the implementation appears quite robust since the relative error actually decreases as the number of carbons (C_{Max}) increases. At a C_{Max} = 100, the relative error is 6.77 Ã— 10^{-17}. This implementation does have some numerical limitations, for example, the C_{Max} must be less than 270 carbons due to all numerical quantities being represented as double precision (64 bit) floating point numbers. However, this limitation is easily overcome by using higher precision floating point numbers.

### Application to phosphatidylcholine 34:1 observed isotopologue intensities

Figure 3 shows the two sets of ^{12}C/^{13}C isotopologue intensities for phosphatidylcholine 34:1 (34 carbons in 2 fatty acid chains with only 1 double bond), with ^{13}C natural abundance contaminated intensities in red and corrected intensities in green. The algorithm converged within 8 iterations to produce the corrected intensity results. In comparing the contaminated and corrected intensities, the most significant changes are seen in isotopologues 0-4 and 16-20. The drastic drop in I_{M+1}, I_{M+2}, and I_{M+4} isotopologues make the incorporation of ^{13}C-labeled glycerol much clearer. Also, the drop in I_{M+16}, I_{M+18}, and I_{M+20} isotopologues supports the expected incorporation of ^{13}C-labeled acetyl groups in the fatty acid chain biosynthesis.

## Discussion and Conclusions

Overall, correcting for the effects of natural abundance makes interpretation of isotopologue intensities from stable isotope tracing experiments easier within the context of cellular metabolism. Such a correction is required before using more quantitative methods of interpretation. Since the relative error is virtually zero with perfectly simulated data and the algorithm is very robust with increasing C_{Max}, the accuracy of this correction is really only limited by the error in the isotopologue intensities themselves. Thus, the propagation of data error through this algorithm should be straightforward to analyze and quantify. Moreover, from Equation 5 it is evident that effects from natural abundance significantly link together groups of observed isotopologue intensities. This difference between calculated and observed intensities should be highly sensitive to the error in a set of isotopologue intensities. Therefore, this difference should be usable as an independent check on the quality of the observed set of isotopologue intensities. Such a quality control check would be especially useful when it is not possible or practical to repeat experiments or to determine whether additional experiments are necessary.

## Methods

### Cell Culture and FT-ICR-MS

We separated glycerophospholipids from crude cell extracts derived from MCF7-LCC2 cells in tissue culture after 24 hours of labeling with uniformly labeled ^{13}C-glucose. We analyzed the sample on a hybrid linear ion trap 7T FT-ICR mass spectrometer (Finnigan LTQ FT, Thermo Electron, Bremen, Germany) equipped with a TriVersa NanoMate ion source (Advion BioSciences, Ithaca, NY) as described elsewhere [6].

## Funding

Supported in part by NSF EPSCoR grant #EPS-0447479.

## References

Wolfe R: Tracers in Metabolic Research: Radioisotope and Stable Isotope.

*Mass Spectrometry Methods*1984, 9: 287.Fan T, Lane A, Higashi R, Farag M, Gao H, Bousamra M, Miller D: Altered regulation of metabolic pathways in human lung cancer discerned by 13 C stable isotope-resolved metabolomics (SIRM).

*Molecular Cancer*2009, 8: 41. 10.1186/1476-4598-8-41Lane A, Fan T, Higashi R: Isotopomer-based metabolomic analysis by NMR and mass spectrometry.

*Methods in cell biology*2008, 84: 541. full_textFischer E, Sauer U: Metabolic flux profiling of Escherichia coli mutants in central carbon metabolism using GC-MS.

*European Journal of Biochemistry*2003, 270(5):880â€“891. 10.1046/j.1432-1033.2003.03448.xFan T, Lane A: Structure-based profiling of metabolites and isotopomers by NMR.

*Progress in Nuclear Magnetic Resonance Spectroscopy*2008, 52(2â€“3):69â€“117. 10.1016/j.pnmrs.2007.03.002Lane A, Fan T, Xie Z, Moseley H, Higashi R: Isotopomer analysis of lipid biosynthesis by high resolution mass spectrometry and NMR.

*Analytica Chimica Acta*2009, 651(2):201â€“208. 10.1016/j.aca.2009.08.032Snider R: Efficient calculation of exact mass isotopic distributions.

*Journal of the American Society for Mass Spectrometry*2007, 18(8):1511â€“1515. 10.1016/j.jasms.2007.05.016Van Winden W, Wittmann C, Heinzle E, Heijnen J: Correcting mass isotopomer distributions for naturally occurring isotopes.

*Biotechnology and bioengineering*2002, 80(4):477â€“479. 10.1002/bit.10393Zhang X, Hines W, Adamec J, Asara J, Naylor S, Regnier F: An automated method for the analysis of stable isotope labeling data in proteomics.

*Journal of the American Society for Mass Spectrometry*2005, 16(7):1181â€“1191. 10.1016/j.jasms.2005.03.016Wehofsky M, Hoffmann R: Automated deconvolution and deisotoping of electrospray mass spectra.

*Journal of Mass Spectrometry*2002, 37(2):223â€“229. 10.1002/jms.278Alecio M, Ferrige A, Pannell L, Ray R: Comparison of Methods for Detecting and Deisotoping Weak, High Charge Signals in Data.

*2006: American Society for Mass Spectrometry*2006.

## Acknowledgements

I thank Drs. T. W-M. Fan. A.N. Lane and R.M. Higashi for support and helpful discussion.

## Author information

### Authors and Affiliations

### Corresponding author

## Additional information

### Authors' contributions

The author derived the analytical solution, implemented the algorithm, tested the implementation, applied the algorithm to the lipid metabolite experimental data, and wrote the manuscript.

## Electronic supplementary material

### 12859_2009_3596_MOESM2_ESM.PL

Additional file 2: **Perl program implementing the algorithm**. Perl program implementing the algorithm displayed in a ASCII text file. (PL 6 KB)

## Authorsâ€™ original submitted files for images

Below are the links to the authorsâ€™ original submitted files for images.

## 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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## About this article

### Cite this article

Moseley, H.N. Correcting for the effects of natural abundance in stable isotope resolved metabolomics experiments involving ultra-high resolution mass spectrometry.
*BMC Bioinformatics* **11**, 139 (2010). https://doi.org/10.1186/1471-2105-11-139

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/1471-2105-11-139

## Comments

View archived comments (1)