Search tips
Search criteria 


Logo of springeropenLink to Publisher's site
Metabolomics. 2017; 13(10): 122.
Published online 2017 September 4. doi:  10.1007/s11306-017-1264-1
PMCID: PMC5583274

A robust, single-injection method for targeted, broad-spectrum plasma metabolomics



Metabolomics is a powerful emerging technology for studying the systems biology and chemistry of health and disease. Current targeted methods are often limited by the number of analytes that can be measured, and/or require multiple injections.


We developed a single-injection, targeted broad-spectrum plasma metabolomic method on a SCIEX Qtrap 5500 LC-ESI-MS/MS platform. Analytical validation was conducted for the reproducibility, linearity, carryover and blood collection tube effects. The method was also clinically validated for its potential utility in the diagnosis of chronic fatigue syndrome (CFS) using a cohort of 22 males CFS and 18 age- and sex-matched controls.


Optimization of LC conditions and MS/MS parameters enabled the measurement of 610 key metabolites from 63 biochemical pathways and 95 stable isotope standards in a 45-minute HILIC method using a single injection without sacrificing sensitivity. The total imprecision (CVtotal) of peak area was 12% for both the control and CFS pools. The 8 metabolites selected in our previous study (PMID: 27573827) performed well in a clinical validation analysis even when the case and control samples were analyzed 1.5 years later on a different instrument by a different investigator, yielding a diagnostic accuracy of 95% (95% CI 85–100%) measured by the area under the ROC curve.


A reliable and reproducible, broad-spectrum, targeted metabolomic method was developed, capable of measuring over 600 metabolites in plasma in a single injection. The method might be a useful tool in helping the diagnosis of CFS or other complex diseases.

Electronic supplementary material

The online version of this article (doi:10.1007/s11306-017-1264-1) contains supplementary material, which is available to authorized users.

Keywords: Targeted metabolomics, LC-MS/MS, Broad-spectrum, Hydrophilic interaction chromatography, Chronic fatigue syndrome, Validation


Metabolomics has recently emerged as a useful tool with broad applications in clinical and translational research such as disease biomarker discovery, clinical drug safety evaluation and precision medicine (Patti et al. 2012; Becker et al. 2012). It can be performed using either untargeted or targeted approaches. Untargeted metabolomics generally requires high-resolution time-of-flight (TOF) instruments to generate precision accurate mass data for all the detectable ions (Spalding et al. 2016). Accurate metabolite identification in untargeted metabolomics is time-consuming and hindered by the requirement for extensive accurate mass databases and sophisticated MS deconvolution and metabolite identification steps. Even though thousands of features (molecular ions) are captured in untargeted metabolomics, only a limited number of metabolites, typically fewer than 0.2–5% of the features, can actually be identified and confirmed (Schymanski and Neumann 2013).

Targeted metabolomics offers several advantages over the untargeted approach. Generally, the targeted metabolomics is conducted using a triple quadrupole mass spectrometer (QqQ, also designated as MS/MS) in multiple reaction monitoring (MRM) mode, which is usually less expensive and more robust than TOF instruments (Ferrario et al. 2016). Since the MRM transitions are specifically optimized to the targeted metabolites using the authentic standards, data analysis is easier. Moreover, MRM-based targeted metabolomics is intrinsically more sensitive, selective, and quantitative, and has a broader dynamic range of 5–6 orders of magnitude than MS based on untargeted metabolomics in complex biological samples using TOF or QTOF instruments (Morin et al. 2013; Dillen et al. 2012).

MRM based targeted analysis by tandem mass spectrometry has become a routine diagnostic tool in clinical laboratories such as vitamin D analysis and the screening of inborn errors of metabolism for newborns (Farrell et al. 2012; la Marca 2014). However, the small number of metabolites measured by typical targeted methods (up to 300 metabolites for single-injection methods) limits the use of targeted metabolomics as a discovery toolbox (Yuan et al. 2012; Bajad et al. 2006). A single LC-MS/MS platform targeting a large scale of metabolites is highly desirable to fulfill different application requirements, reduce the cost and increase the productivity of metabolomics.

The analysis of complex samples requires efficient chromatographic separation of the metabolites prior to MS to reduce ion suppression and improve signal sensitivity. Hydrophilic interaction chromatography (HILIC) has proved to be an effective method for the separation of water-soluble metabolites that are only minimally retained in the reversed-phase liquid chromatography (RPLC) (Bajad et al. 2006; Ivanisevic et al. 2013). HILIC has also been optimized for the separation of phospholipids (Losito et al. 2013; Schwalbe-Herrmann et al. 2010).

Chronic fatigue syndrome (CFS) is a complex, multisystem disease that affects more than 2 million people in the United States (Prins et al. 2006). The symptoms of CFS include, but not limited to: post-exertional malaise, disabling fatigue, chronic pain, and nonrestorative sleep (Clayton 2015). No biochemical diagnostic lab test is currently available. Our previous study identified unique metabolic signatures in the plasma of patients with CFS using the targeted metabolomic approach (Naviaux et al. 2016). In this study, we described the full method with details on key analytical procedures such as method optimization, MRM transitions, linearity and carryover which would help other researchers in the field to adapt and apply this method to their own work. To confirm the inter-instrument and inter-tester reproducibility of our metabolomic method, we also reanalyzed 40 plasma samples from the previous CFS study on a different instrument by a different investigator.

Materials and methods

Chemicals and material

Unlabeled chemical standards were purchased from Sigma–Aldrich, Cayman Chemical, TCI America, Fisher Scientific, or Avanti Polar Lipids. HPLC grade (NH4)2CO3 (Catalog No. A651-500) and ACS grade NH4OH (Catalog No. AC423300250) were from ACROS Organics, Fisher Scientific. LC-MS grade water (Catalog No. LC365) and acetonitrile (Catalog No. LC015) were purchased from Honeywell Burdick & Jackson. Commercial stable isotopes were obtained from Cambridge Isotope. Blood collection tubes including Lithium heparin (Catalog No. 367,884), sodium heparin (Catalog No. 367,871), dipotassium (K2)-EDTA (Catalog No. 367,861) and serum separator tube (SSTs, Order No. 367,983) were from BD. Custom-synthesized uniformly 13C-labeled metabolite internal standards were produced via metabolic labeling by growing Escherichia coli NCM 3277, Caenorhabditis elegans N2, and Komagataella phaffii with 13C6-glucose and 13C-bicarbonate as sole carbon sources.

Metabolites extraction

All human samples were collected with signed informed consent in accordance with a human subject protocol approved by University of California, San Diego Institutional Review Board (IRB #14-0072). Plasma samples (90 µl) were mixed with 5 µl of commercial internal standards and 5 µl of custom-synthesized 13C labeled standards. The mixture was then extracted with 400 µl of cold (−20 °C) methanol: acetonitrile (50:50 MEOH:ACN). The extract was stored at −80 °C prior to LC-MS/MS. Detailed procedures were described in the online Supplemental Materials.

LC-MS/MS analysis

Details for the optimization of LC and MS/MS conditions were provided in the online Supplemental Materials. The optimal conditions are described below.

LC-MS/MS analysis was performed using a UFLC XR HPLC system (LC-20AD, Shimadzu) coupled to a Turbo V electrospray ionization source and a Qtrap 5500 mass spectrometer (SCIEX). The samples were separated in HILIC mode using a polymer-based NH2 4 µm, 250 × 2 mm column (Asahipak NH2P-40 2E, Shodex). The optimal LC conditions were as follows: Mobile phase A: 95% H2O with 20 mM (NH4)2CO3 and 5% ACN, pH 9.8. Mobile phase B: 100% ACN. The gradient was: 0–3.5 min 95% B, 3.6–8 min 85% B, 8.1–13 min 75% B, 14–30 min 0% B, 31–41 min 95% B, 41.1 min stop. The flow rate was 200 µl/min and the injection volume was 10 µl. Typical pump pressures ranged from 850 to 2800 psi over the course of the run.

The MS/MS detection was performed using electrospray ionization (ESI) and by scheduled multiple reaction monitoring (MRM). All 610 metabolites and 95 internal standards were targeted in a single injection using both negative and positive modes with rapid polarity switching (50 ms) and advanced scheduled MRM algorithm. The ESI source conditions were set as follows: electrospray voltage of −4500 V for negative mode and 5500 V for positive mode, source temperature of 500 °C, curtain gas of 30, ion source gas 1 and gas 2 of 35 psi, respectively. The mass transitions and compound-dependent parameters were optimized using the authentic standards and listed in Table S1. This table is formatted for easy compatibility with Analyst 1.6 software (SCIEX).

Analytical validation

A balanced precision experiment was performed to evaluate the reproducibility of the method using five replicates of the control and pooled CFS plasma analyzed on each of 5 days (Grant and Hoofnagle 2014). The mean intra-assay (CVintra), inter assay (CVinter) and total CVs for each pool were calculated. The linearity was evaluated using 5-point admixtures of control and CFS pools (100, 75, 50, 25, and 0% of the CFS pool in the mixture). Duplicates were performed for each point. The peak quality and carryover were estimated using the pooled adult human plasma spiked with stable isotope internal standards (Quality control samples, QC). Effect of blood collection tubes was evaluated using the samples collected by four types of tubes including Li-heparin, Na-heparin, K2-EDTA and serum separator tube (SSTs).

Details for the procedures of analytical validation were described in online Supplemental Materials.

Inter-instrument and inter-tester reproducibility

Clinical validation consisted of a cohort of 22 male CFS (53 ± 2.8 years old, mean ± SEM, range 21–67 years) and 18 male controls (53 ± 3.5 years old, mean ± SEM, range 23–69 years), who were enrolled in our previous study (Naviaux et al. 2016). These plasma samples were stored in −80 °C for about 1.5 years and reanalyzed on a different LC-MS/MS system by a different investigator. Metabolomic data was log 2 transformed and analyzed by partial least squares discriminant analysis (PLS-DA) in MetaboAnalyst ( Eight metabolites were selected for testing as a diagnostic classifier in the original study (Naviaux et al. 2016). Overfitting was minimized by selecting classifier metabolites to interrogate altered pathways using a combination of random forest analysis (Breiman 2001), forward selection, and backward elimination methods to identify a set of metabolites-a “biosignature”-that performed well as a diagnostic classifier by receiver operator characteristic (ROC) curve analysis (Szymanska et al. 2012; Xia et al. 2013). These same eight metabolites were then tested in the current study, after repeating the sample metabolomic analysis in this study, 1.5 years later. Classifier robustness was estimated by repeated double cross-validation (rdCV) and permutation testing 1000 times in MetaboAnalyst (Table S2). The correlation of peak area, z-scores and pathway impact scores between two analyses were also used to check the reproducibility of the method for the diagnosis of CFS.

Data analysis

The peak quality evaluation metrics including retention time (RT), peak area, peak width, peak width at 50% height, and signal to noise ratio (S/N) were generated in MultiQuant 3.0 (SCIEX) and log 2 transformed prior to statistical analysis. Coefficient of variation (CV), Pearson correlation, frequency distribution and linear regression analysis were conducted in Graphpad Prism 6.0. PLS-DA, heatmap, volcano plot analysis and ROC curve analysis were performed in metaboanalyst (


Targeted metabolites

Purified authentic standards were used to determine the MRM transitions, compound-specific parameters, and column retention times. A total of 610 endogenous metabolites were finally selected based on the availability of standards and reliable quantification on the LC columns (Fig. 1a, b). These metabolites cover a broad range of chemical classes and known to be essential in a variety of health and disease states. Table S1 lists the optimal MRM transition, source parameters, chemical class and the metabolic pathway of each metabolite.

Fig. 1
Total 705 metabolites were targeted on both negative and positive mode in a single run. a The number of targeted metabolites in the method. b Targeted endogenous metabolites with diverse chemical classes

The application of stable isotopes in metabolomics helps the correction of extraction loss, matrix effect, ion suppression and MS signal drift. They are also critical for inter-batch normalization for large-scale metabolomic studies. In this study, a mix of commercial stable isotopes and in-house prepared fully 13C-labeled metabolites were created and 95 stable isotope standards were targeted (Fig. 1a). In total, 705 metabolites were analyzed using the single injection method.

Method optimization

A detailed description of optimization steps and analysis is described in the Online Supplemental Results.


Aminopropyl columns from three different manufacturers: Luna NH2 (Phenomenex), inertsil NH2 (GL Science) and Asahipak NH2P-40 (Shodex) were evaluated. Based on our observation, we found the polymer-based column from Shodex to have a longer lifetime, typically over 1000 injections, and to be more stable at high pH than the silica based aminopropyl columns. Further optimization was conducted using a Shodex 250 × 2 mm aminopropyl column.

LC conditions

Mobile phase compositions, pH, and oven temperatures were carefully investigated to achieve the best separation (Fig. S1 and S2). As shown in Fig. 2a, the targeted 705 analytes were eluted across the entire gradient. In a given 2-min window, the maximum number of compounds scanned was only about 150 (Fig. 2b). Some representative examples are shown in Fig. 2c–f. Six different classes of membrane lipids were separated on the HILIC column including ceramides, glycosphingolipids, sphingomyelins, cardiolipins, phospholipids and bis(monoacylglycero)phosphates (Fig. 2c, d). The method also showed excellent separation of four types of glycerophospholipids including phosphatidylcholine (PC, RT: 4.58 min), phosphatidylethanolamine (PE, RT: 9.65 min), phosphatidylinositol A (PI, RT: 10.10 min), and phosphatidylserine (PS, RT: 12.64) (Fig. 2c, d). Baseline separation of several isobaric metabolites, including alanine and sarcosine (Fig. 2e) and methylmalonic acid and succinic acid (Fig. 2f) was achieved using this method.

Fig. 2
The chromatographic separation of 705 targeted metabolites on a shodex polymer-based aminopropyl column in HILIC mode in a single run. a The total ion chromatogram (TIC) of a quality control plasma extract analyzed using the developed method; b the frequency ...

MS/MS conditions

We optimized the source temperature and ion source gas using quality control (QC) plasma extract lots to determine the conditions suitable for our diverse classes of targeted metabolites. Optimal source parameters for our method were ion source gas at 35 psi, and source temperature set at 500 °C (Fig. S3). Traditional LC-MS/MS based metabolomics often requires two separate injections for negative and positive ions (Pluskal et al. 2010). In combination with the rapid polarity switching and scheduled MRM algorithm, we quantified all 705 MRM transitions in a single run without significantly reducing the sensitivity (Fig. S4).

Analytical validation

Peak quality evaluation

To characterize the quality of the peaks generated using the developed method, we performed the frequency distribution analysis of peak area, width, width at 50% of the peak height, points across the peak, and signal to noise ratio (S/N) for all metabolites detected in plasma. All the peak metrics showed a good fit to a lognormal distribution (Fig. 3). Peak width is an important factor that determines the actual dwell time spent on the peak and points across the peak in scheduled MRM algorithm. The average peak width of the targeted analytes was 35.6 s with 95% confidence interval (CI) between 10.1 and 124.9 s and the mean of peak width at 50% of peak height was 5.4 s (Fig. 3a, c). The average number of data points across the peak was 20, which ensures the accurate and reproducible quantification of the analytes (Fig. 3b). The 95% CI of signal to noise (S/N) for the targeted metabolites was between 8 and 1427 with the mean of 107 (Fig. 3d).

Fig. 3
Frequency distribution analysis of peak quality for the targeted metabolites in human plasma. a Histogram of peak width, b histogram of points across the peaks, c histogram of peak width at 50% of peak height and d histogram of signal to noise ratio (S/N). ...


Generating reproducible results is essential for the success of large-scale metabolomic studies that usually involve hundreds of samples and require multiple batches to complete the analysis of a large cohort of samples. The average CV for intra-batch chromatographic retention time (RT) variation was 0.11% (95% CI 0.02–0.34%) for the control pool and 0.12% (95% CI 0.02–0.39) for the CFS pool. The inter-batch RT variation was also very small, with average CVs of 0.14% (95% CI 0.02–0.44%) for the control pool and 0.15% (95% CI 0.02–0.54%) for the CFS pool. This is particularly important because the scheduled MRM algorithm relies on the stability of analyte retention time.

The peak area (AUCs) reproducibility was also excellent. The average CVintra and CVinter of peak area for the control pool was 7.4 and 9.5%, respectively, which was below the acceptance limit (20% RSD) of FDA’s guideline for bio-analytical method validation using LC-MS (Tiwari and Tiwari 2010) (Fig. 4a, c). Similarly, the average CVintra and CVinter of the peak area for the CFS pool was only 7.5 and 9.7% (Fig. 4b, c). The total variability was 12% for both control and the CFS pools.

Fig. 4
The peak area reproducibility of the developed targeted metabolomic method. a The intra-batch CV of the control pool for each of 5 days. b The intra-batch CV of the CFS pool for each of 5 days. c The inter-batch CV of the control and CFS pool ...


The linearity was evaluated by mixing the plasma of the healthy controls with CFS patients at 0:1, 1:3, 1:1, 3:1, and 1:0 ratios. Our method showed good linearity with regression coefficients >0.9 for the representative metabolites used for the diagnosis of CFS (Fig. 5a–f). This also indicated that the metabolic biomarkers identified using the developed metabolomic method have the capability for resolving meaningful differences between CFS patients and controls of as little as 20%. Eighty-nine percent (89%) of the detected metabolites had r value >0.90 (Fig. S5).

Fig. 5
Linearity evaluation. af Six representative metabolites selected for distinguishing male CFS patients from the controls. The linearity was evaluated by mixing the control pool with CFS pool at the ratios of 0:1, 1:3, 1:1, 3:1, and 1:0 ratios. ...

Effect of blood collection tubes

We investigated the differences in metabolomic profile between four different types of blood collection tubes including K2-EDTA, Na-heparin, Li-heparin and SSTs. Our results indicated that the anticoagulant has an impact on the detection of plasma metabolites (Fig.S6a). The most striking differences were the significant increase of phosphorylated nucleotides and a 10–50-fold decrease of nucleosides like adenosine in EDTA plasma compared with Li-heparin plasma (Fig.S6c). In contrast, the metabolic differences between Na-heparin and Li-heparin plasma were relatively small (Fig.S6d). In addition, there were dramatic metabolomic differences between serum and plasma with higher levels of eicosanoids, nucleotides and serotonin in serum samples (Fig. S6b and S6d).


The average carryover, expressed as the AUC percentage of the first blank observed after injection of a sample containing the upper limit of quantification for each analyte, was only 0.75%. Eighty-five percent of the analytes measured had carryovers of less than 2%, and none had more than 4% (Fig. S7).

Inter-instrument and inter-tester reproducibility

To test the inter-instrument and inter-tester reproducibility of our developed metabolomic method which is critical for analytical method transfer, we reanalyzed the 40 CFS and control plasma samples used in our previous study on a different instrument and by a different investigator. The AUCs of 335 metabolites in both control and CFS groups were well correlated between original study and the validation analysis (Fig. 6a, b). Z-score (absolute value >0.8) in the validation study was also well correlated with the previous values (r2 = 0.866, p < 0.0001, Fig. 6c). Similar to the original study, PLS-DA in the repeated analysis showed a distinct separation of male CFS from the control group (Fig. 6d). Pathway analysis revealed the significant correlation of pathway impact scores of disturbed pathways in CFS identified by both analyses (Fig. 6e). The 8 metabolite biomarkers identified in our previous study still performed well in distinguishing controls from CFS patients in the validation by AUROC analysis of 95% (95% CI 85–100%) (Fig. 6f).

Fig. 6
Inter-instrument and inter-tester validation of the developed targeted metabolomic method for its utility in the diagnosis of CFS. a AUCs of the control groups in the validation study were correlated with those in the original study (r2 = 0.879, ...


In this study, we developed, optimized, analytically and clinically validated a robust, broad-spectrum targeted metabolomic workflow based on HILIC-MS/MS (Fig.S8). The key features of our method can be summarized as follows:

Expanded coverage of the metabolome

Technical improvements in new QqQ instruments include faster scan speeds (as low as 2 ms dwell time for a single transition), rapid polarity switching (≤50 ms), high sensitivity, quantifiability, and large dynamic ranges of 5–6 orders of magnitude, making them ideal for targeted metabolomics. Previously published targeted metabolomic methods typically measure about 200–300 polar metabolites, or require several injections to accommodate positively and negatively charged analytes (Xiao et al. 2012; Zhou et al. 2016). In the present study, we dramatically expanded the metabolome coverage of targeted metabolomics. Using a single injection method, a total of 610 endogenous metabolites and 95 isotopes were targeted. These metabolites covered 63 major metabolic pathways involved in various critical physiological and pathological processes such as glycolysis, TCA cycle, purine and pyrimidine nucleosides and nucleotides, amino acid, fatty acid, phospholipid, and sphingolipid metabolism.

Powerful separation of diverse classes of metabolites

The greatest challenge in targeted metabolomics comes from the diversity of metabolites that differ greatly in their physical and chemical properties of size, hydrophobicity and charge (Dunn and Ellis 2005). In order to analyze a broader range of metabolites simultaneously, we have focused our efforts on the optimization of LC column separation performance. The separation of cell membrane lipids from the complex biological extract on a HILIC column can be challenging. Glycerophospholipids had been reported to be separated on HILIC columns (Cifkova et al. 2016; Zhu et al. 2012). Our method showed the separation of six different classes of membrane lipids (ceramides, glycosphingolipids, sphingomyelins, cardiolipins, glycerophospholipids bis(monoacylglycero)phosphates), and four subclasses of glycerophospholipids (PC, PI, PE and PS) on the Shodex column in HILIC mode.

The possibility for decreased sensitivity is always a concern when targeting a large number of metabolites in both negative and positive modes simultaneously. The good separation of diverse classes of molecules in our method enabled the quantification of all 705 metabolites in a single run without sacrificing sensitivity. Mean signal to noise ratio was 107, which indicated that the targeted metabolites could be accurately identified.

Our method allowed the baseline separation of some isomers of interest, which previously required chemical derivatization or special columns for accurate quantification. These included alanine and its isomer, sarcosine (Chung et al. 2015), and succinic acid and its isomer, methylmalonic acid (Carvalho and Kok 2008) (Fig. 2e, f). For isomeric species that could not be separated in our method, we measured the total of all molecular species such as the hexose monosaccharide pool (glucose, fructose, galactose, etc).

Enhanced reproducibility and stability

One of the common challenges with HILIC technique is chromatographic reproducibility (Cubbon et al. 2010). Retained salts from repeated injections of biological samples and the high pH in the mobile phase can potentially cause column degradation in a short period and result in peak broadening and sensitivity reduction (Contrepois et al. 2015). In this study, the polymer-based aminopropyl column from Shodex exhibited excellent intra-batch and long-term inter-batch reproducibility. In addition, the Shodex column was stable under the high pH condition for over 800–1000 injections without significant changes in peak shape, column pressure or sensitivity.

Another factor affecting data reproducibility is the number of points collected across the peak. Scanning too many metabolites at the same time can increase the cycle time and reduce the points across the peak, which in turn reduces peak area reproducibility. In our method, the gradient was carefully optimized so that even in the most intense scanning period in the gradient (RT 16–18 min), only about 150 metabolites were simultaneously targeted by the scheduled MRM algorithm. The average number of points across the peak were 20, which were enough for the generation of reproducible peak area each time.

Selection of sample collection tubes

Bias in sample collection is a frequently forgotten problem responsible for uncontrolled errors in metabolomic analysis and can have a major impact on the identification of biomarkers. Our results showed clear differences between K2-EDTA and Li-heparin plasma that related to the specific class of metabolites measured. The choice of tubes should therefore be made based on the aim of the study. If the study focuses on the analysis of phosphorylated nucleotides, EDTA tubes have some advantages, but come at the cost of several weaknesses. For more general, targeted, metabolomic analysis across many different classes of molecules, Li-heparin tubes are recommended. Our findings are in line with others who have recommended EDTA plasma for untargeted studies (Yin et al. 2013) and Li-heparin plasma for targeted studies (Dunn et al. 2011). Serum samples, which avoid the issue of chemical additives, but require a clotting step, may also have utility in some metabolomics studies. The coagulation cascade during the preparation of serum affected the levels of certain classes of metabolites such as eicosanoids, nucleotides and serotonin. When using archived samples from biobanks, critical attention should be paid to the type of collection tube that was used. Collection tube chemistry must be standardized for each metabolomic experimental series.

Inter-instrument and inter-tester validation for disease diagnosis

This study successfully replicated our recent findings of the unique chemical signature in CFS patients, which could be used for the diagnosis of CFS (Naviaux et al. 2016). Although the samples had been stored in −80 °C for about 1.5 years and a small number of metabolites might be changed because of degradation or ion suppression effects, the chemical signature and the diagnostic biomarkers remained unchanged. Good linearity was still achieved for these biomarkers in the 5-point mixing experiment. To validate the method independently, we repeated the analysis on this 40-sample study on a different Qtrap 5500 system by a different investigator. Even though nearly 95% of the absolute peak areas were different due to sensitivity and tuning differences between two instruments, the eight metabolites selected in our original study still performed with a 95% ROC curve accuracy in distinguishing controls from CFS patients (Fig. 6). And perhaps even more importantly, pathway analysis of the results validated the original study results, emphasizing the importance of sphingolipids and phospholipids in the biology of CFS in this cohort (Naviaux et al. 2016).

Analytical limitations

A mass resolution of ±0.7 da is a known limitation of QqQ mass spectrometers, which rely on characteristic fragmentation patterns to identify metabolites. Occasionally, isobaric molecules cannot be resolved in our method such as hexose sugars. The hexose pool was quantified. Peak splitting was observed for some of long-chain acylcarnitines such palmitoylcarnitine (C16) and myristoylcarnitine (C14) on the HILIC column. In addition, some neurotransmitters (norepinephrine, epinephrine and melatonin) were not well detected in our method. A reverse phase method might be an alternative solution for some of these metabolites.


We report an optimized method for targeted, broad-spectrum metabolomic analysis of 610 molecules from over 60 biochemical pathways in a single injection. This method performed well in rigorous analytical validation and was able to reproduce the metabolomic results from samples stored for over 1.5 years in a clinical validation experiment.

Electronic supplementary material

Below is the link to the electronic supplementary material.


The authors thank Jonathan Monk for computational and bioinformatics assistance. This work was supported by gifts from the UCSD Christini Fund, the Lennox Foundation, the Wright Family Foundation, the UCSD Mitochondrial Research Fund, and the Jane Botsford Johnson Foundation.

Author contributions

All authors confirmed they have contributed to the intellectual content of this paper and have met the following 3 requirements: (a) significant contributions to the conception and design, acquisition of data, or analysis and interpretation of data; (b) drafting or revising the article for intellectual content; and (c) final approval of the published article.


This work was supported by gifts from the UCSD Christini Fund, the Lennox Foundation, the Wright Family Foundation, the UCSD Mitochondrial Research Fund, and the Jane Botsford Johnson Foundation.

Compliance with ethical standards

Compliance with ethical standards

Conflict of interest

Upon manuscript submission, all authors completed the author disclosure form. Disclosures and/or potential conflicts of interest.

Consultant or advisory role

RKN is a scientific advisory board member for the Autism Research Institute (ARI) and the Open Medicine Foundation (OMF).

Role of sponsors

The funding organizations played no role in the design of the study, review and interpretation of data, and final approval of the manuscript.


Electronic supplementary material

The online version of this article (doi:10.1007/s11306-017-1264-1) contains supplementary material, which is available to authorized users.

Contributor Information

Kefeng Li, ude.dscu@ilk.

Robert K. Naviaux, ude.dscu@xuaivan.


  • Bajad SU, Lu W, Kimball EH, Yuan J, Peterson C, Rabinowitz JD. Separation and quantitation of water soluble cellular metabolites by hydrophilic interaction chromatography-tandem mass spectrometry. Journal of Chromatography. A. 2006;1125(1):76–88. doi: 10.1016/j.chroma.2006.05.019. [PubMed] [Cross Ref]
  • Becker S, Kortz L, Helmschrodt C, Thiery J, Ceglarek U. LC-MS-based metabolomics in the clinical laboratory. Journal of Chromatography. B, Analytical Technologies in the Biomedical and Life Sciences. 2012;883–884:68–75. doi: 10.1016/j.jchromb.2011.10.018. [PubMed] [Cross Ref]
  • Breiman L. Random forests. Machine Learning. 2001;45(1):5–32. doi: 10.1023/A:1010933404324. [Cross Ref]
  • Carvalho VM, Kok F. Determination of serum methylmalonic acid by alkylative extraction and liquid chromatography coupled to tandem mass spectrometry. Analytical Biochemistry. 2008;381(1):67–73. doi: 10.1016/j.ab.2008.06.023. [PubMed] [Cross Ref]
  • Chung TC, Li CT, Kou HS, Wu HL. High-performance liquid chromatographic analysis of sarcosine as a fluorescent levofloxacin derivative. Journal of Chromatographic Science. 2015;53(8):1310–1315. doi: 10.1093/chromsci/bmv010. [PubMed] [Cross Ref]
  • Cifkova E, Hajek R, Lisa M, Holcapek M. Hydrophilic interaction liquid chromatography-mass spectrometry of (lyso)phosphatidic acids, (lyso)phosphatidylserines and other lipid classes. Journal of Chromatography. A. 2016;1439:65–73. doi: 10.1016/j.chroma.2016.01.064. [PubMed] [Cross Ref]
  • Clayton EW. Beyond myalgic encephalomyelitis/chronic fatigue syndrome: An IOM report on redefining an illness. JAMA: The Journal of the American Medical Association. 2015;313(11):1101–1102. doi: 10.1001/jama.2015.1346. [PubMed] [Cross Ref]
  • Contrepois K, Jiang LH, Snyder M. Optimized analytical procedures for the untargeted metabolomic profiling of human urine and plasma by combining hydrophilic interaction (HILIC) and reverse-phase liquid chromatography (RPLC)-mass spectrometry. Molecular & Cellular Proteomics. 2015;14(6):1684–1695. doi: 10.1074/mcp.M114.046508. [PMC free article] [PubMed] [Cross Ref]
  • Cubbon S, Antonio C, Wilson J, Thomas-Oates J. Metabolomic applications of Hilic-Lc-Ms. Mass Spectrometry Reviews. 2010;29(5):671–684. doi: 10.1002/mas.20252. [PubMed] [Cross Ref]
  • Dillen L, Cools W, Vereyken L, Lorreyne W, Huybrechts T, de Vries R, et al. Comparison of triple quadrupole and high-resolution TOF-MS for quantification of peptides. Bioanalysis. 2012;4(5):565–579. doi: 10.4155/bio.12.3. [PubMed] [Cross Ref]
  • Dunn WB, Broadhurst D, Begley P, Zelena E, Francis-McIntyre S, Anderson N, et al. Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry. Nature Protocols. 2011;6(7):1060–1083. doi: 10.1038/nprot.2011.335. [PubMed] [Cross Ref]
  • Dunn WB, Ellis DI. Metabolomics: Current analytical platforms and methodologies. Trac-Trends in Analytical Chemistry. 2005;24(4):285–294. doi: 10.1016/j.trac.2004.11.021. [Cross Ref]
  • Farrell CJ, Martin S, McWhinney B, Straub I, Williams P, Herrmann M. State-of-the-art vitamin D assays: A comparison of automated immunoassays with liquid chromatography-tandem mass spectrometry methods. Clinical Chemistry. 2012;58(3):531–542. doi: 10.1373/clinchem.2011.172155. [PubMed] [Cross Ref]
  • Ferrario M, Cambiaghi A, Brunelli L, Giordano S, Caironi P, Guatteri L, et al. Mortality prediction in patients with severe septic shock: A pilot study using a target metabolomics approach. Scientific Reports. 2016;6:20391. doi: 10.1038/srep20391. [PMC free article] [PubMed] [Cross Ref]
  • Grant RP, Hoofnagle AN. From lost in translation to paradise found: Enabling protein biomarker method transfer by mass spectrometry. Clinical Chemistry. 2014;60(7):941–944. doi: 10.1373/clinchem.2014.224840. [PMC free article] [PubMed] [Cross Ref]
  • Ivanisevic J, Zhu ZJ, Plate L, Tautenhahn R, Chen S, O’Brien PJ, et al. Toward ‘omic scale metabolite profiling: A dual separation-mass spectrometry approach for coverage of lipid and central carbon metabolism. Analytical Chemistry. 2013;85(14):6876–6884. doi: 10.1021/ac401140h. [PMC free article] [PubMed] [Cross Ref]
  • la Marca G. Mass spectrometry in clinical chemistry: The case of newborn screening. Journal of Pharmaceutical and Biomedical Analysis. 2014;101:174–182. doi: 10.1016/j.jpba.2014.03.047. [PubMed] [Cross Ref]
  • Losito I, Patruno R, Conte E, Cataldi TR, Megli FM, Palmisano F. Phospholipidomics of human blood microparticles. Analytical Chemistry. 2013;85(13):6405–6413. doi: 10.1021/ac400829r. [PubMed] [Cross Ref]
  • Morin LP, Mess JN, Garofolo F. Large-molecule quantification: Sensitivity and selectivity head-to-head comparison of triple quadrupole with Q-TOF. Bioanalysis. 2013;5(10):1181–1193. doi: 10.4155/bio.13.87. [PubMed] [Cross Ref]
  • Naviaux RK, Naviaux JC, Li K, Bright AT, Alaynick WA, Wang L, et al. Metabolic features of chronic fatigue syndrome. Proceedings of theNational Academyof Sciences of the United States of America. 2016;113(37):E5472-5480. [PubMed]
  • Patti GJ, Yanes O, Siuzdak G. Innovation: Metabolomics: The apogee of the omics trilogy. Nature Reviews. Molecular Cell Biology. 2012;13(4):263–269. doi: 10.1038/nrm3314. [PMC free article] [PubMed] [Cross Ref]
  • Pluskal T, Nakamura T, Villar-Briones A, Yanagida M. Metabolic profiling of the fission yeast S. pombe: Quantification of compounds under different temperatures and genetic perturbation. Molecular Biosystems. 2010;6(1):182–198. doi: 10.1039/B908784B. [PubMed] [Cross Ref]
  • Prins JB, van der Meer JW, Bleijenberg G. Chronic fatigue syndrome. Lancet. 2006;367(9507):346–355. doi: 10.1016/S0140-6736(06)68073-2. [PubMed] [Cross Ref]
  • Schwalbe-Herrmann M, Willmann J, Leibfritz D. Separation of phospholipid classes by hydrophilic interaction chromatography detected by electrospray ionization mass spectrometry. Journal of Chromatography A. 2010;1217(32):5179–5183. doi: 10.1016/j.chroma.2010.05.014. [PubMed] [Cross Ref]
  • Schymanski EL, Neumann S. The Critical Assessment of small molecule identification (CASMI): Challenges and solutions. Metabolites. 2013;3(3):517–538. doi: 10.3390/metabo3030517. [PMC free article] [PubMed] [Cross Ref]
  • Spalding JL, Cho K, Mahieu NG, Nikolskiy I, Llufrio EM, Johnson SL, et al. Bar coding MS2 spectra for metabolite identification. Analytical Chemistry. 2016;88(5):2538–2542. doi: 10.1021/acs.analchem.5b04925. [PMC free article] [PubMed] [Cross Ref]
  • Szymanska E, Saccenti E, Smilde AK, Westerhuis JA. Double-check: Validation of diagnostic statistics for PLS-DA models in metabolomics studies. Metabolomics. 2012;8(Suppl 1):3–16. doi: 10.1007/s11306-011-0330-3. [PMC free article] [PubMed] [Cross Ref]
  • Tiwari G, Tiwari R. Bioanalytical method validation: An updated review. Pharm Methods. 2010;1(1):25–38. doi: 10.1016/S2229-4708(10)11004-8. [PMC free article] [PubMed] [Cross Ref]
  • Xia J, Broadhurst DI, Wilson M, Wishart DS. Translational biomarker discovery in clinical metabolomics: An introductory tutorial. Metabolomics. 2013;9(2):280–299. doi: 10.1007/s11306-012-0482-9. [PMC free article] [PubMed] [Cross Ref]
  • Xiao JF, Zhou B, Ressom HW. Metabolite identification and quantitation in LC-MS/MS-based metabolomics. Trends Analytical Chemistry. 2012;32:1–14. doi: 10.1016/j.trac.2011.08.009. [PMC free article] [PubMed] [Cross Ref]
  • Yin PY, Peter A, Franken H, Zhao XJ, Neukamm SS, Rosenbaum L, et al. Preanalytical aspects and sample quality assessment in metabolomics studies of human blood. Clinical Chemistry. 2013;59(5):833–845. doi: 10.1373/clinchem.2012.199257. [PubMed] [Cross Ref]
  • Yuan M, Breitkopf SB, Yang X, Asara JM. A positive/negative ion-switching, targeted mass spectrometry-based metabolomics platform for bodily fluids, cells, and fresh and fixed tissue. Nature Protocols. 2012;7(5):872–881. doi: 10.1038/nprot.2012.024. [PMC free article] [PubMed] [Cross Ref]
  • Zhou J, Liu H, Liu Y, Liu J, Zhao X, Yin Y. Development and evaluation of a parallel reaction monitoring strategy for large-scale targeted metabolomics quantification. Analytical Chemistry. 2016;88(8):4478–4486. doi: 10.1021/acs.analchem.6b00355. [PubMed] [Cross Ref]
  • Zhu C, Dane A, Spijksma G, Wang M, van der Greef J, Luo G, et al. An efficient hydrophilic interaction liquid chromatography separation of 7 phospholipid classes based on a diol column. Journal of Chromatography A. 2012;1220:26–34. doi: 10.1016/j.chroma.2011.11.034. [PubMed] [Cross Ref]

Articles from Springer Open Choice are provided here courtesy of Springer