Next Article in Journal
Potentially Treatable Disorder Diagnosed Post Mortem by Exome Analysis in a Boy with Respiratory Distress
Next Article in Special Issue
Transcriptome and Gene Ontology (GO) Enrichment Analysis Reveals Genes Involved in Biotin Metabolism That Affect l-Lysine Production in Corynebacterium glutamicum
Previous Article in Journal
Tamoxifen Treatment of Breast Cancer Cells: Impact on Hedgehog/GLI1 Signaling
Previous Article in Special Issue
Multiple Evolutionary Selections Involved in Synonymous Codon Usages in the Streptococcus agalactiae Genome
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolomic Profiling of Plasma from Melioidosis Patients Using UHPLC-QTOF MS Reveals Novel Biomarkers for Diagnosis

1
State Key Laboratory of Emerging Infectious Diseases, The University of Hong Kong, Pokfulam, Hong Kong, China
2
Research Centre of Infection and Immunology, The University of Hong Kong, Pokfulam, Hong Kong, China
3
Carol Yu Centre for Infection, The University of Hong Kong, Pokfulam, Hong Kong, China
4
Department of Microbiology, The University of Hong Kong, Pokfulam, Hong Kong, China
5
Department of Medicine, The University of Hong Kong, Pokfulam, Hong Kong, China
6
Department of Pathology, The University of Hong Kong, Pokfulam, Hong Kong, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2016, 17(3), 307; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms17030307
Submission received: 8 January 2016 / Revised: 15 February 2016 / Accepted: 22 February 2016 / Published: 27 February 2016
(This article belongs to the Special Issue Microbial Genomics and Metabolomics)

Abstract

:
To identify potential biomarkers for improving diagnosis of melioidosis, we compared plasma metabolome profiles of melioidosis patients compared to patients with other bacteremia and controls without active infection, using ultra-high-performance liquid chromatography-electrospray ionization-quadruple time-of-flight mass spectrometry. Principal component analysis (PCA) showed that the metabolomic profiles of melioidosis patients are distinguishable from bacteremia patients and controls. Using multivariate and univariate analysis, 12 significant metabolites from four lipid classes, acylcarnitine (n = 6), lysophosphatidylethanolamine (LysoPE) (n = 3), sphingomyelins (SM) (n = 2) and phosphatidylcholine (PC) (n = 1), with significantly higher levels in melioidosis patients than bacteremia patients and controls, were identified. Ten of the 12 metabolites showed area-under-receiver operating characteristic curve (AUC) >0.80 when compared both between melioidosis and bacteremia patients, and between melioidosis patients and controls. SM(d18:2/16:0) possessed the largest AUC when compared, both between melioidosis and bacteremia patients (AUC 0.998, sensitivity 100% and specificity 91.7%), and between melioidosis patients and controls (AUC 1.000, sensitivity 96.7% and specificity 100%). Our results indicate that metabolome profiling might serve as a promising approach for diagnosis of melioidosis using patient plasma, with SM(d18:2/16:0) representing a potential biomarker. Since the 12 metabolites were related to various pathways for energy and lipid metabolism, further studies may reveal their possible role in the pathogenesis and host response in melioidosis.

Graphical Abstract

1. Introduction

Melioidosis is a disease caused by the highly pathogenic gram-negative bacterium, Burkholderia pseudomallei (B. pseudomallei). The disease is often serious and potentially fatal, most commonly manifested as severe community-acquired pneumonia and sepsis. Although the disease is mainly endemic in Southeast Asia and Northern Australia, melioidosis has been increasingly reported in areas of the Asia-Pacific region, including India [1,2], Mauritius [3], South, Central and North America [4,5,6], and West and East Africa [7,8], which may suggest an expanding geographical distribution. B. pseudomallei is known to be a natural saprophyte, and therefore melioidosis is believed to be acquired through contact with contaminated soil and water in the environment [9,10]. Illness can be presented as an acute, subacute, or chronic process, with an incubation period of up to 26 years [11]. The disease manifestations can range from subclinical infection, localized abscesses, and pneumonia to fulminant sepsis, leading to a mortality rate of up to 19% [12]. Besides human, melioidosis also affects various animals in endemic areas [10,13]. Treatment of melioidosis is often difficult, as B. pseudomallei is usually resistant to multiple antibiotics and prolonged antibiotics are required to prevent relapse [14,15]. Moreover, diagnostic and therapeutic resources in endemic areas are often limited, which have hindered efforts to improve treatment outcomes.
Diagnosis of melioidosis can be difficult because of several reasons. First, B. pseudomallei may not be readily isolated from clinical specimens. Second, even if it is successfully isolated, commercial bacterial identification systems often cannot differentiate B. pseudomallei from closely related species such as Burkholderia thailandensis and members of Burkholderia cepacia complex [16]. Therefore, new molecular techniques are often required for more accurate species identification [14,17,18,19,20,21,22,23]. Despite these new technologies, the diagnostic problems associated with culture-negative cases remain unresolved. Although different serological tests have been developed to help diagnose culture-negative melioidosis, their clinical usefulness is limited by the low sensitivities and specificities [24,25]. The availability of alternative techniques for improved diagnosis of melioidosis is thus eagerly awaited, and such techniques should be able to differentiate between melioidosis and infections caused by common Gram-negative bacteria including the closely related Burkholderia species.
Metabolomics is a new research platform for systematic studies of small molecules of a specific system such as a cell, tissue or an organism. The metabolic profiles between such systems can then be compared, thus allowing the identification of specific metabolite markers. The technique has been used to characterize various infectious diseases or pathogens [26,27,28,29,30,31,32,33,34,35,36,37,38]. By exploring the metabolomes of culture supernatant, we have identified specific biomarkers that are produced by a unique thiamine degradation pathway in B. pseudomallei [32]. We have also recently reported the use of metabolomics to identify novel biomarkers in plasma of tuberculosis patients, which may be useful for diagnosis [37]. Despite being an important pathogen, no studies have reported the use of metabolomics to explore specific biomarkers in plasma of melioidosis patients.
We hypothesize that there are specific biomarkers that may be detected in plasma of melioidosis patients. To identify potential biomarkers for the non-invasive diagnosis of melioidosis, we applied the metabolomics technology for metabolite profiling of plasma samples from melioidosis patients, using ultra-high-performance liquid chromatography-electrospray ionization-quadrupole time-of-flight mass spectrometry (UHPLC-ESI-QTOFMS). Multi- and univariate statistical analyses of the metabolome data were used to identify specific metabolites that are present in significantly higher levels in plasma of melioidosis patients than in plasma of patients with other bacteremia or controls without infections. The diagnostic performances of the identified biomarkers were evaluated using receiver operating characteristic curve (ROC) analysis. In this pilot study, untargeted metabolomics on plasma sample were conducted with the aim to explore potential diagnostic biomarkers and biological pathways involved in host–B. pseudomallei interaction.

2. Results

2.1. Metabolomic Profiling of Plasma Samples from Melioidosis Patients, Patients with Other Bacteremia and Controls

The metabolomes of 76 plasma samples from the three groups (22 samples from newly-diagnosed melioidosis patients, 24 samples from patients with other bacteremia and 30 samples from controls without active infection) were compared using UHPLC-QTOFMS [38]. The metabolites could be separated well using the UPLC-MS method with sub-micron particle size, 1.7 µm packing. The base peak chromatographic profiles showed stable retention time for all peaks without observable drift, supporting the stability and reliability of the accurate-mass QTOF system and metabolomic profiling data.

2.2. Omics-Based Statistical and Bioinformatic Analysis for Identification of Biomarkers

A total of 2424 molecular features were obtained by XCMS package [39] and subjected to MetaboAnalyst 3.0 software [40] for statistical analysis. For multivariate analysis, principal component analysis (PCA) revealed that the three groups were clustered separately, with 50.0% of the total variance among the three groups represented by the first two principal components (PCs), where principal component 1 (PC1) and PC2 explained 36.5% and 13.5% of the variance, respectively (Figure 1). In particular, the melioidosis group could be distinguished from the bacteremia and control groups based on the first two PCs, with clear separation along PC2 dimension.
Univariate analysis using one-way analysis of variance (ANOVA) identified 764 statistically significant features with variable importance in the projection (VIP) score >1 and p < 0.05 when compared both between melioidosis patients and patients with other bacteremia, and between melioidosis patients and controls. Further volcano plot analysis revealed 131 significant features with fold-change (FC) >1.5 and p < 0.05 by Student’s t-test when compared both between melioidosis patients and patients with other bacteremia (Figure 2A), and between melioidosis patients and controls (Figure 2B), which were subjected to univariate ROC analysis.
A total of 12 metabolites with area-under-receiver operating characteristic curve (AUC) >0.80 when compared between melioidosis patients and patients with other bacteremia, or between melioidosis patients and controls, were identified. They were identified as metabolites belonging to four lipid classes, acylcarnitine (six metabolites), lysophosphatidylethanolamine (LysoPE) (three metabolites), sphingomyelins (SM) (two metabolites) and phosphatidylcholine (PC) (one metabolite), using LC-MS/MS analyses by elution order, MS/MS fragmentation and predicted molecular formula (Table 1 and Table 2, Figure 3). The identification and assignment of lipid classes were based on the fingerprint fragment and specific neutral losses. Specifically, the presence of m/z 60 (trimethylamine ion, C3H10N+) and m/z 85 (C4H5O2+) for acylcarnitine, m/z 44 (C2H5N+) and neutral loss of phosphorylethanolamine for LysoPE, m/z 184 (Phosphocholine ion, C5H15NO4P+) and neutral loss of water for SM, and m/z 60 (trimethylamine ion, C3H10N+) and m/z 184 (phosphocholine ion, C5H15NO4P+) for PC, were key features employed for their identification. In addition, the identities of six biomarkers, including decanoylcarnitine, decanoylcarnitine, and l-octanoylcarnitine LysoPE(16:0/0:0), LysoPE(18:0/0:0) and PC(16:0/16:0), were confirmed by matching the retention time (RT) and MS/MS fragmentation patterns of authentic chemical standards, where available. Chemical standards for the other six metabolites were not available for comparison. The details of the fragments in each MS/MS spectrum for each identified metabolite are shown in Figure 3. The 12 metabolites were related to various pathways for energy and lipid metabolism, including tricarboxylic acid (TCA) cycle, fatty acid β-oxidation, fatty acid de novo synthesis, linoleic acid, α-linoleic acid, arachidonic acid, phospholipid and sphingolipid metabolism.

2.3. Diagnostic Performance of Metabolites

The AUC, sensitivity and specificity for ROC curves calculated for the 12 metabolites at optimal cutoffs are summarized in Table 2. Box-whisker plots revealed that they all exhibited significantly higher levels in plasma samples of melioidosis patients than in samples of patients with other bacteremia and controls (p < 0.01 by Student’s t-test) (Figure 4). Among the 12 biomarkers, 10 showed AUC > 0.80 when compared both between melioidosis and bacteremia patients, and between melioidosis patients and controls. SM(d18:2/16:0) possessed the largest AUC when compared both between melioidosis patients and patients with other bacteremia (AUC 0.998, sensitivity 100% and specificity 91.7%), and between melioidosis patients and controls (AUC 1.000, sensitivity 96.7% and specificity 100%) (Figure 5).

3. Discussion

Using metabolomics approach, we identified 12 novel biomarkers in plasma of melioidosis patients with significantly higher levels than in plasma of patients with other bacteremia and controls without active infection. In this study, samples from patients with other bacteremia were included because these patients may present with sepsis mimicking melioidosis where B. pseudomallei bacteremia can occur. Therefore, the present biomarkers may help differentiate melioidosis from bacteremia caused by common bacterial species such as Escherichia coli. In particular, SM(d18:2/16:0) represents the most promising biomarker for melioidosis, with AUC of 0.998, sensitivity of 100.0% and specificity of 91.7% when compared to other bacteremia, and AUC of 1.000, sensitivity of 96.7% and specificity of 100.0% when compared to controls without active infection. However, the present study is limited by the small number of patients with melioidosis included, which is partly due to clinical difficulties in the diagnosis of melioidosis and our relatively low disease prevalence when compared to other Southeast Asian countries such as Thailand. Further studies with inclusion of more cases from endemic regions are required to validate the diagnostic potential of the present biomarkers.
The high plasma concentrations of six medium- to long-chain (C13 to C19) acylcarnitines in melioidosis patients may reflect changes in fatty acid (FA) β-oxidation during infection. Acylcarnitines are synthesized from acyl-CoAs with transfer of hydroxyl group of carnitine, and are transported from the intermembraneous space into the mitochondrial matrix for FA β-oxidation [41]. In mice infected with B. pseudomallei, transcriptomics studies have revealed changes in transcript levels of various enzymes involved in FA β-oxidation [42]. Interestingly, higher levels of medium-chain (C5 to C10) acylcarnitines have also been found in plasma of patients with systemic inflammatory response syndrome (SIRS) due to severe sepsis/septic shock than those with SIRS due to non-infective causes [43]. Since accumulation of medium-chain acyl-CoAs in the mitochondria is toxic for the cell, the synthesis of acylcarnitines may have a protective effect during infection [43]. Further studies may help elucidate the role of these medium- to long-chain acylcarnitines in the host response to melioidosis.
The high levels of three LysoPE, LysoPE (16:0/0:0), LysoPE (0:0/18:0) and LysoPE (18:0/0:0), in plasma of melioidosis patients may be the result of changes in phospholipid metabolism or cellular damage from systemic infection. LysoPE, a constituent of cell membranes, is derived from the hydrolysis of PE, which is catalyzed by phospholipase A2 (PLA2) [44]. It has been shown that LysoPEs can stimulate invariant natural killer T cell activation through self-antigenicity, suggesting a possible role in innate immunity during infection [45]. In addition, LysoPE was shown to enhance the ingestion activity of macrophage on IgG-coated target cells via Fc receptors [46]. It remains to be determined if LysoPE may be involved in the innate immunity against B. pseudomallei. Our findings may also be in line with previous observations of higher plasma levels of type II PLA2 in sepsis patients, which significantly correlated with TNF-α, IL-6 and IL-8 levels [47]. On the other hand, elevated plasma levels of LysoPE (16:0/0:0) have been observed in rats with induced-liver injuries as a result of massive destruction of cell membranes [48,49]. Therefore, it is also possible that the high LysoPE levels may be the consequence of severe organ damage during systemic melioidosis.
The high levels of SM(d16:1/16:0) and SM(d18:2/16:0) may reflect changes in sphingolipid metabolism during melioidosis. SMs, structural components of cell membranes, are synthesized by the transfer of PC to ceramides by sphingomyelin synthase and degraded back to ceramides by sphingomyelinase [50]. While these SMs are likely produced by the host, it is interesting to note that the genome of B. pseudomallei also possessed genes homologous to hemolytic phospholipase C (PlcH) of Pseudomonas aeruginosa with sphingomyelin synthase and sphingomyelinase activity [50,51,52]. Further studies are required to better understand the possible role of sphingomyelin metabolism in the pathogenesis of melioidosis.
The PC(16:0/16:0) detected in plasma of melioidosis patients may be the result of innate immune response in the lungs. PC(16:0/16:0) or dipalmitoylphosphatidylcholine (DPCC) is a major component of pulmonary surfactant [53,54]. Melioidosis is believed to be acquired through inhalation of contaminated aerosals. B. pseudomallei has been found to induce pro-inflammatory cytokines from macrophages and alveolar type II pneumocytes (ATII) cells [55], the latter being responsible for secretion of surfactants. We speculate that the production of PC(16:0/16:0) may be upregulated in ATII cells when in contact with B. pseudomallei, which may be secreted in plasma leading to the elevated PC(16:0/16:0) levels. However, comparison with plasma of patients with other causes of pneumonia would be important to determine if PC(16:0/16:0) is specific to melioidosis or may represent a general biomarker for pneumonia.

4. Materials and Methods

4.1. Patient and Control Samples

Clinical samples were collected from patients hospitalized in Queen Mary Hospital, Hong Kong. A total of 22 plasma samples from five patients with newly-diagnosed melioidosis, 24 plasma samples from 24 patients with bacteremia caused by other bacterial species and 30 controls without active infections were included for UHPLC-QTOFMS analysis. Plasma samples from melioidosis patients were collected before commencement of antibiotic treatment. The diagnosis of melioidosis was made according to compatible clinical features, and with either isolation of B. pseudomallei from clinical samples and/or positive antibodies against B. pseudomallei as determined by enzyme-linked immunosorbent assay as described previously [14]. All five patients with melioidosis were immunocompromised with underlying diseases, including HIV, autoimmune vasculitis, diabetes mellitus, lymphoma and acute myeloid leukemia. B. pseudomallei was isolated from the blood cultures of three patients, while the other two cases were diagnosed by positive antibodies against B. pseudomallei. Plasma samples were collected at admission from patients with other bacteremia were used. Causative agents of other bacteremia included Aeromonas caviae (n = 1), Bacteroides (n = 1), E. coli (n = 17), Klebisella pnuemoniae (n = 1), Prevotella species (n = 1), Proteus mirabilis (n = 2), Streptococcus mitis (n = 1), and Streptococcus pneumoniae (n = 2) (two patients had two different bacterial isolates from the same blood culture). Controls included patients with no clinical evidence of active infection. This study has been approved by the Institutional Review Board, the University of Hong Kong/Hospital Authority of Hong Kong West Cluster under reference number UW 13-265.

4.2. Chemicals and Reagents

LC-MS grade water, methanol and acetonitrile were purchased from J.T. Baker (Center Valley, PA, USA). High-performance liquid chromatography (HPLC)-grade ethanol and acetone were purchased from Merck & Co. (Kenilworth, NJ, USA). Formic Acid was of American Chemical Society (ACS) reagent grade from Sigma-Aldrich (Saint Louis, MO, USA). Decanoylcarnitine, LysoPE(16:0/0:0), LysoPE(18:0/0:0) and PC(16:0/16:0) was purchased from Avanti Polar Lipid (Alabaster, AL, USA). Decanoylcarnitine and l-octanoylcarnitine were purchased from Sigma-Aldrich (Saint Louis, MO, USA).

4.3. Sample Preparation

Blood samples were collected in heparin bottles, transferred immediately to the laboratory, and centrifuged at 3000 rpm at 4 °C for 10 min to obtain the plasma fractions. For metabolomics analysis, 100 µL of plasma was thawed at 4 °C and plasma proteins were precipitated with 400 µL of methanol/ethanol/acetone mixture at a ratio of 1:1:1 (v/v/v). The sample extract was vigorously vortexed for 1 min, and centrifuged at 14,000 rpm at 4 °C for 10 min. The supernatant was collected for UHPLC-ESI-QTOFMS analysis. All specimens were immediately kept at −80 °C until analysis and stored within one week. The thawed specimens were analyzed within 48 h in a random manner to prevent the batch effect.

4.4. Untargeted Metabolomics Profiling of Patient Plasma Using UHPLC-ESI-QTOFMS

The metabolomic profiling of plasma supernatants was performed as describe previously with modifications [37], using Agilent 1290 Infinity UHPLC (Agilent Technologies, Waldbronn, Germany) coupled with Agilent 6540 UHD Accurate-Mass QTOF system (Agilent Technologies, Santa Clara, CA, USA) accompanied with a MassHunter Workstation software for QTOF (version B.03.01 for Data Acquisition, Agilent Technologies, Santa Clara, CA, USA). Waters Acquity UPLC BEH C18 column (2.1 × 100 mm, 1.7 µm) (Waters, Milford, MA, USA) was used for the separation with the injection volume of 5 µL. The column and autosampler temperature were maintained at 45 and 10 °C, respectively. The separation was performed at a flow rate of 0.4 mL/min under a gradient program in which mobile phase A was composed of LC-MS grade water containing 0.1% formic acid (v/v) and mobile phase B was composed of acetonitrile. The gradient program was applied as follows: t = 0 min, 5% B; t = 0.5 min, 5% B; t = 7 min, 48% B; t = 20 min, 78% B; t = 27 min, 80% B; t = 31 min, 99.5% B; t = 36.5 min, 99.5% B; t = 36.51 min, 5% B. The stop time was 40 min. The ESI mass spectra were acquired in both positive and negative ion modes using Agilent Jet Stream ESI source (Agilent Technologies, Santa, CA, USA) with capillary voltages at +3800 and −3500 V, respectively. Other source conditions were kept constant in all the experiments as follow: gas temperature was kept constant at 300 °C, drying gas (nitrogen) was set at the rate of 7 L/min, and the pressure of nebulizer gas (nitrogen) was 40 psi. The sheath gas was kept at a flow rate of 10 L/min at a temperature of 350 °C. The voltages of the Fragmentor, Skimmer 1, and OctopoleRFPeak were 135, 50 and 500 V, respectively. The mass data were collected between m/z 80 and 1700 at the acquisition rate of 2 scans per second. Two reference masses at m/z 121.0509 (protonated molecular ion of C5H4N4) and m/z 922.0098 (protonated molecular ion of C18H18O6N3P3F24) for positive mode, and m/z 119.0363 (deprotonated molecular ion of C5H4N4) and m/z 966.0007 (formate adduct of C18H18O6N3P3F24) for negative mode were used as constant mass correction during LC-MS run. Product ion scanning experiments were conducted using ultra-high purity N2 as collision energy with same parameters set in MS acquisition, the collision energy (CE) was set at 10, 20 or 40 eV to generate the best quality of MS/MS spectra for the putative identification and structural elucidation of the significant metabolites.

4.5. Data Processing and Statistical Analysis

Multivariate and univariate analysis was performed to identify molecular features that discriminate melioidosis patients from patients with other bacteremia patients and controls as described previously [37]. Multivariate analysis was performed on a total of 76 LC-MS data of plasma samples from three groups (22, 24 and 30 samples from melioidosis patients, bacteremia with other bacteremia and controls without active infection respectively). The raw LC-MS data were converted into mzData format using Agilent MassHunter Qualitative Analysis software (version B.05.00, Agilent Technologies, USA) and subsequently processed using open-source XCMS package [39] operating in R [56], which adopted different peak detection and alignment as well as data filtering with centWave algorithms. Data was further processed with normalization, scaling, filtering and statistical analysis using MetaboAnalyst 3.0 [40]. The data were mean-centered, square root scaled and normalized such that the sum of squares for each chromatogram equaled on for statistical analysis [57]. Insignificant features between melioidosis patients and patients with bacteremia or controls were filtered out using both uni- and multivariate analyses. For multivariate analysis, PCA was performed for unsupervised analysis on all LC-MS features using MetaboAnalyst 3.0.
For univariate analysis, statistical significance of features was determined among melioidosis patients, patients with other bacteremia and controls using one-way ANOVA with Tukey’s post-hoc test. p < 0.05 was considered to be statistically significant. Significant features with FC >1.5 by volcano plots and p < 0.05 by Student’s t-test between melioidosis patients and patients with bacteremia, and between melioidosis patients and controls were identified. Common significant features were subject to univariate ROC analysis using web-based ROCCET [58]. The classical ROC curve analysis was performed and AUC was calculated by Monte Carlo Cross Validation (MCCV) using sub-sampling. In addition, the optimal cutoffs for the given metabolite were computed to obtain the sensitivity, specificity, and confidence intervals at different cut-offs for the evaluation of the recognition and prediction ability with respect to each variable. Significant features with AUC ≥0.8 obtained from either comparison between melioidosis patients and patients with other bacteremia, or between melioidosis patients and controls were identified. Box-whisker plots were generated and p values were calculated by the Student’s t-test using Analyse-it software (Analyse-it Software, Leeds, UK). Multivariate ROC curves were further generated using ROCCET. The procedures were repeated multiple times to calculate the performance and confidence interval of the model using support vector machines (SVM). The predicted class probabilities for each sample were evaluated with the best classifier (based on AUC) with confusion matrix.

4.6. Metabolite Identification

Features with significant differences were selected for product ion scanning (PIS) experiments. MS/MS spectra for the potential biomarkers and commercially available reference standards, including decanoylcarnitine, decanoylcarnitine, l-octanoylcarnitine, LysoPE(16:0/0:0), LysoPE(18:0/0:0) and PC(16:0/16:0), were processed using Agilent MassHunter Qualitative Analysis software (version B.05.00, Agilent Technologies, USA) to generate potential molecular formula based on the accurate mass and isotopic pattern recognitions of parent and fragment ions. All putative identities were confirmed by matching with entries in the METLIN database [59], Human Metabolome Database (HMDB) [60], MassBank [61], LipidMaps [62], KEGG (Kyoto Encyclopedia of Genes and Genomes) database [63] using exact molecular weights, nitrogen rule or MS/MS fragmentation pattern data and literature search. Efforts were made to distinguish metabolites from the other isobaric compounds whenever possible by its elution order and virtue of difference in fragmentation pattern corresponding to its structural characteristics. Putative identities of six biomarkers were confirmed by comparing their chromatographic RT and MS/MS spectra with those obtained from commercially available standards of decanoylcarnitine, decanoylcarnitine, l-octanoylcarnitine, LysoPE(16:0/0:0), LysoPE(18:0/0:0) and PC(16:0/16:0). The collision energy was set at 20 eV for generating the MS/MS spectra.

5. Conclusions

In this study, we compared the metabolome profiles of plasma of patients with melioidosis to those of patients with other bacteremia and controls without active infections. Twelve significant metabolites with significantly higher levels were identified in melioidosis patients than bacteremia patients and controls. These 12 metabolites, including l-octanoylcarnitine, decanoylcarnitine, dodecanoylcarnitine, LysoPE(16:0/0:0), LysoPE(18:0/0:0), PC(16:0/16:0), LysoPE(0:0/18:0), l-hexanoylcarnitine, SM(d16:1/16:0), 2-decenoylcarnitine, SM(d18:2/16:0) and trans-2-dodecenoylcarnitine, are involved in various pathways for energy and lipid metabolism. The present study demonstrates the potential of metabolomics in identifying novel biomarkers in studying infectious diseases. Further studies may reveal the potential of these metabolites as diagnostics biomarkers for melioidosis and their possible role in the pathogenesis and host response in melioidosis.

Acknowledgments

This work is partly supported by the Hong Kong Special Administrative RegionResearch Fund for the Control of Infectious Diseases (Commissioned Study HK-09-01-10) of the Food and Health Bureau; University Development Fund and Strategic Research Fund, The University of Hong Kong; and donation from Eunice Lam on emerging infectious disease research.

Author Contributions

Susanna K. P. Lau and Patrick C. Y. Woo conceived and designed the experiments; Kim-Chung Lee, George C. S. Lo, Vanessa S. Y. Ding, Wang-Ngai Chow, Tony Y. H. Ke and Shirly O. T. Curreem performed the experiments; Susanna K. P. Lau, Kim-Chung Lee, George C. S. Lo, Vanessa S. Y. Ding, Wang-Ngai Chow, Tony Y. H. Ke and Shirly O. T. Curreem analyzed the data; Susanna K. P. Lau, Kelvin K. W. To, Deborah T. Y. Ho, Siddharth Sridhar, Sally C. Y. Wong, Jasper F. W. Chan, Ivan F. N. Hung, Kong-Hong Sze, Ching-Wan Lam, Kwok-Yung Yuen and Patrick C. Y. Woo contributed reagents/materials/analysis tools; and Susanna K. P. Lau, Kim-Chung Lee, George C. S. Lo, Vanessa S. Y. Ding, Wang-Ngai Chow and Patrick C. Y. Woo wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ACSAmerican Chemical Society
ANOVAanalysis of variance
ATIIalveolar type II pneumocytes
AUCarea-under-receiver operating characteristic curve
CEcollision energy
DPCCdipalmitoylphosphatidylcholine
FAfatty acid
FCfold-change
IL-6interleukin-6
IL-8interleukin-8
LysoPElysophosphatidylethanolamine
MCCVMonte Carlo Cross Validation
PCphosphatidylcholine
PCAprincipal component analysis
PISproduct ion scanning
PLA2phospholipase A2
PlcHphospholipase C
ROCreceiver operating characteristic curve
RTretention time
SMsphingomyelins
SIRSsystemic inflammatory response syndrome
SVMsupport vector machines
TCAtricarboxylic acid
TNF-αTumor necrosis factor-α
UHPLC-ESI-QTOFMSultra-high-performance liquid chromatography-electrospray ionization-quadruple time-of-flight-mass spectrometry
VIPvariable importance in the projection

References

  1. Anuradha, K.; Meena, A.K.; Lakshmi, V. Isolation of Burkholderia pseudomallei from a case of septicaemia—A case report. Indian J. Med. Microbiol. 2003, 21, 129–132. [Google Scholar] [PubMed]
  2. Saravu, K.; Mukhopadhyay, C.; Vishwanath, S.; Valsalan, R.; Docherla, M.; Vandana, K.E.; Shastry, B.A.; Bairy, I.; Rao, S.P. Melioidosis in southern India: Epidemiological and clinical profile. Southeast. Asian J. Trop. Med. Public Health 2010, 41, 401–409. [Google Scholar] [PubMed]
  3. Issack, M.I.; Bundhun, C.D.; Gokhool, H. Melioidosis in Mauritius. Emerg. Infect. Dis. 2005, 11, 139–140. [Google Scholar] [CrossRef] [PubMed]
  4. Inglis, T.J.; Rolim, D.B.; Sousa Ade, Q. Melioidosis in the Americas. Am. J. Trop Med. Hyg. 2006, 75, 947–954. [Google Scholar] [PubMed]
  5. Rolim, D.B.; Vilar, D.C.; Sousa, A.Q.; Miralles, I.S.; de Oliveira, D.C.; Harnett, G.; O’Reilly, L.; Howard, K.; Sampson, I.; Inglis, T.J. Melioidosis, northeastern Brazil. Emerg. Infect. Dis 2005, 11, 1458–1460. [Google Scholar] [CrossRef] [PubMed]
  6. Stewart, T.; Engelthaler, D.M.; Blaney, D.D.; Tuanyok, A.; Wangsness, E.; Smith, T.L.; Pearson, T.; Komatsu, K.K.; Keim, P.; Currie, B.J.; et al. Epidemiology and investigation of melioidosis, southern Arizona. Emerg. Infect. Dis. 2011, 17, 1286–1288. [Google Scholar] [CrossRef] [PubMed]
  7. Salam, A.P.; Khan, N.; Malnick, H.; Kenna, D.T.; Dance, D.A.; Klein, J.L. Melioidosis acquired by traveler to Nigeria. Emerg. Infect. Dis 2011, 17, 1296–1298. [Google Scholar] [CrossRef] [PubMed]
  8. Cuadros, J.; Gil, H.; Miguel, J.D.; Marabe, G.; Gomez-Herruz, T.A.; Lobo, B.; Marcos, R.; Anda, P. Case report: Melioidosis imported from west Africa to Europe. Am. J. Trop. Med. Hyg. 2011, 85, 282–284. [Google Scholar] [CrossRef] [PubMed]
  9. Currie, B.J.; Mayo, M.; Anstey, N.M.; Donohoe, P.; Haase, A.; Kemp, D.J. A cluster of melioidosis cases from an endemic region is clonal and is linked to the water supply using molecular typing of Burkholderia pseudomallei isolates. Am. J. Trop. Med. Hyg. 2001, 65, 177–179. [Google Scholar] [PubMed]
  10. Lau, S.K.; Chan, S.Y.; Curreem, S.O.; Hui, S.W.; Lau, C.C.; Lee, P.; Ho, C.C.; Martelli, P.; Woo, P.C. Burkholderia pseudomallei in soil samples from an oceanarium in Hong Kong detected using a sensitive PCR assay. Emerg. Microbes Infect. 2014, 3, e69. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Mays, E.E.; Ricketts, E.A. Melioidosis: Recrudescence associated with bronchogenic carcinoma twenty-six years following initial geographic exposure. Chest 1975, 68, 261–263. [Google Scholar] [CrossRef] [PubMed]
  12. Currie, B.J.; Fisher, D.A.; Howard, D.M.; Burrow, J.N.; Lo, D.; Selva-Nayagam, S.; Anstey, N.M.; Huffam, S.E.; Snelling, P.L.; Marks, P.J.; et al. Endemic melioidosis in tropical northern Australia: A 10-year prospective study and review of the literature. Clin. Infect. Dis. 2000, 31, 981–986. [Google Scholar] [CrossRef] [PubMed]
  13. Dance, D.A. Melioidosis. Curr. Opin. Infect. Dis. 2002, 15, 127–132. [Google Scholar] [CrossRef] [PubMed]
  14. Woo, P.C.; Lau, S.K.; Woo, G.K.; Fung, A.M.; Ngan, A.H.; Hui, W.T.; Yuen, K.Y. Seronegative bacteremic melioidosis caused by Burkholderia pseudomallei with ambiguous biochemical profile: Clinical importance of accurate identification by 16s rRNA gene and groEL gene sequencing. J. Clin. Microbiol. 2003, 41, 3973–3977. [Google Scholar] [CrossRef] [PubMed]
  15. Cheng, A.C.; Currie, B.J. Melioidosis: Epidemiology, pathophysiology, and management. Clin. Microbiol. Rev. 2005, 18, 383–416. [Google Scholar] [CrossRef] [PubMed]
  16. Weissert, C.; Dollenmaier, G.; Rafeiner, P.; Riehm, J.; Schultze, D. Burkholderia pseudomallei misidentified by automated system. Emerg. Infect. Dis. 2009, 15, 1799–1801. [Google Scholar] [CrossRef] [PubMed]
  17. Godoy, D.; Randle, G.; Simpson, A.J.; Aanensen, D.M.; Pitt, T.L.; Kinoshita, R.; Spratt, B.G. Multilocus sequence typing and evolutionary relationships among the causative agents of melioidosis and glanders, Burkholderia pseudomallei and Burkholderia mallei. J. Clin. Microbiol. 2003, 41, 2068–2079. [Google Scholar] [CrossRef] [PubMed]
  18. Lau, S.K.; Tang, B.S.; Curreem, S.O.; Chan, T.M.; Martelli, P.; Tse, C.W.; Wu, A.K.; Yuen, K.Y.; Woo, P.C. Matrix-assisted laser desorption ionization-time of flight mass spectrometry for rapid identification of Burkholderia pseudomallei: Importance of expanding databases with pathogens endemic to different localities. J. Clin. Microbiol. 2012, 50, 3142–3143. [Google Scholar] [CrossRef] [PubMed]
  19. Inglis, T.J.; Healy, P.E.; Fremlin, L.J.; Golledge, C.L. Use of matrix-assisted laser desorption/ionization time-of-flight mass spectrometry analysis for rapid confirmation of Burkholderia pseudomallei in septicemic melioidosis. Am. J. Trop Med. Hyg. 2012, 86, 1039–1042. [Google Scholar] [CrossRef] [PubMed]
  20. Payne, G.W.; Vandamme, P.; Morgan, S.H.; Lipuma, J.J.; Coenye, T.; Weightman, A.J.; Jones, T.H.; Mahenthiralingam, E. Development of a recA gene-based identification approach for the entire Burkholderia genus. Appl. Environ. Microbiol. 2005, 71, 3917–3927. [Google Scholar] [CrossRef] [PubMed]
  21. Schmoock, G.; Ehricht, R.; Melzer, F.; Rassbach, A.; Scholz, H.C.; Neubauer, H.; Sachse, K.; Mota, R.A.; Saqib, M.; Elschner, M. DNA microarray-based detection and identification of Burkholderia mallei, Burkholderia pseudomallei and Burkholderia spp. Mol. Cell. Probes 2009, 23, 178–187. [Google Scholar] [CrossRef] [PubMed]
  22. Woo, P.C.; Woo, G.K.; Lau, S.K.; Wong, S.S.; Yuen, K. Single gene target bacterial identification. groEL gene sequencing for discriminating clinical isolates of Burkholderia pseudomallei and Burkholderia thailandensis. Diagn. Microbiol. Infect. Dis. 2002, 44, 143–149. [Google Scholar] [CrossRef]
  23. Ho, C.C.; Lau, C.C.; Martelli, P.; Chan, S.Y.; Tse, C.W.; Wu, A.K.; Yuen, K.Y.; Lau, S.K.; Woo, P.C. Novel pan-genomic analysis approach in target selection for multiplex PCR identification and detection of Burkholderia pseudomallei, Burkholderia thailandensis, and Burkholderia cepacia complex species: A proof-of-concept study. J. Clin. Microbiol. 2011, 49, 814–821. [Google Scholar] [CrossRef] [PubMed]
  24. Puthucheary, S.D.; Anuar, A.S.; Tee, T.S. Burkholderia thailandensis whole cell antigen cross-reacts with B. pseudomallei antibodies from patients with melioidosis in an immunofluorescent assay. Southeast. Asian J. Trop. Med. Public Health 2010, 41, 395–400. [Google Scholar] [PubMed]
  25. Limmathurotsakul, D.; Chantratita, N.; Teerawattanasook, N.; Piriyagitpaiboon, K.; Thanwisai, A.; Wuthiekanun, V.; Day, N.P.; Cooper, B.; Peacock, S.J. Enzyme-linked immunosorbent assay for the diagnosis of melioidosis: Better than we thought. Clin. Infect. Dis. 2011, 52, 1024–1028. [Google Scholar] [CrossRef] [PubMed]
  26. Olivier, I.; Loots, D.T. A metabolomics approach to characterise and identify various Mycobacterium species. J. Microbiol. Methods 2012, 88, 419–426. [Google Scholar] [CrossRef] [PubMed]
  27. Olivier, I.; Loots, D.T. An overview of tuberculosis treatments and diagnostics. What role could metabolomics play? J. Cell. Tissue Res. 2011, 11, 2655–2671. [Google Scholar]
  28. Meissner-Roloff, R.J.; Koekemoer, G.; Warren, R.M.; Loots, D.T. A metabolomics investigation of a hyper- and hypo-virulent phenotype of Beijing lineage m. Tuberculosis. Metabolomics 2012, 8, 1194–1203. [Google Scholar] [CrossRef]
  29. Tam, E.W.; Chen, J.H.; Lau, E.C.; Ngan, A.H.; Fung, K.S.; Lee, K.C.; Lam, C.W.; Yuen, K.Y.; Lau, S.K.; Woo, P.C. Misidentification of Aspergillus nomius and Aspergillus tamarii as Aspergillus flavus: Characterization by internal transcribed spacer, β-tubulin, and calmodulin gene sequencing, metabolic fingerprinting, and matrix-assisted laser desorption ionization-time of flight mass spectrometry. J. Clin. Microbiol. 2014, 52, 1153–1160. [Google Scholar] [PubMed]
  30. To, K.K.; Fung, A.M.; Teng, J.L.; Curreem, S.O.; Lee, K.C.; Yuen, K.Y.; Lam, C.W.; Lau, S.K.; Woo, P.C. Characterization of a Tsukamurella pseudo-outbreak by phenotypic tests, 16s rRNA sequencing, pulsed-field gel electrophoresis, and metabolic footprinting. J. Clin. Microbiol. 2013, 51, 334–338. [Google Scholar] [CrossRef] [PubMed]
  31. Woo, P.C.; Lam, C.W.; Tam, E.W.; Leung, C.K.; Wong, S.S.; Lau, S.K.; Yuen, K.Y. First discovery of two polyketide synthase genes for mitorubrinic acid and mitorubrinol yellow pigment biosynthesis and implications in virulence of Penicillium marneffei. PLoS Negl. Trop. Dis. 2012, 6, e1871. [Google Scholar] [CrossRef] [PubMed]
  32. Lau, S.K.; Lam, C.W.; Curreem, S.O.; Lee, K.C.; Chow, W.N.; Lau, C.C.; Sridhar, S.; Wong, S.C.; Martelli, P.; Hui, S.W.; et al. Metabolomic profiling of Burkholderia pseudomallei using uhplc-esi-q-tof-ms reveals specific biomarkers including 4-methyl-5-thiazoleethanol and unique thiamine degradation pathway. Cell. Biosci. 2015, 5, 26. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Lee, K.C.; Tam, E.W.; Lo, K.C.; Tsang, A.K.; Lau, C.C.; To, K.K.; Chan, J.F.; Lam, C.W.; Yuen, K.Y.; Lau, S.K.; et al. Metabolomics analysis reveals specific novel tetrapeptide and potential anti-inflammatory metabolites in pathogenic Aspergillus species. Int. J. Mol. Sci. 2015, 16, 13850–13867. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Lau, S.K.; Lam, C.W.; Curreem, S.O.; Lee, K.C.; Lau, C.C.; Chow, W.N.; Ngan, A.H.; To, K.K.; Chan, J.F.; Hung, I.F.; et al. Identification of specific metabolites in culture supernatant of Mycobacterium tuberculosis using metabolomics: Exploration of potential biomarkers. Emerg. Microbes Infect. 2015, 4, e6. [Google Scholar] [CrossRef] [PubMed]
  35. Madigan, C.A.; Cheng, T.Y.; Layre, E.; Young, D.C.; McConnell, M.J.; Debono, C.A.; Murry, J.P.; Wei, J.R.; Barry, C.E., 3rd; Rodriguez, G.M.; et al. Lipidomic discovery of deoxysiderophores reveals a revised mycobactin biosynthesis pathway in Mycobacterium tuberculosis. Proc. Natl. Acad. Sci. USA 2012, 109, 1257–1262. [Google Scholar] [CrossRef] [PubMed]
  36. Layre, E.; Lee, H.J.; Young, D.C.; Martinot, A.J.; Buter, J.; Minnaard, A.J.; Annand, J.W.; Fortune, S.M.; Snider, B.B.; Matsunaga, I.; et al. Molecular profiling of Mycobacterium tuberculosis identifies tuberculosinyl nucleoside products of the virulence-associated enzyme rv3378c. Proc. Natl. Acad. Sci. USA 2014, 111, 2978–2983. [Google Scholar] [CrossRef] [PubMed]
  37. Lau, S.K.; Lee, K.C.; Curreem, S.O.; Chow, W.N.; To, K.K.; Hung, I.F.; Ho, D.T.; Sridhar, S.; Li, I.W.S.; Wu, A.L.; et al. Plasma metabolome changes in patients with tuberculosis using untargeted uhplc-qtof ms method reveal novel biomarkers for diagnosis. J. Clin. Microbiol. 2015, 53, 3750–3759. [Google Scholar] [PubMed]
  38. To, K.K.; Lee, K.C.; Wong, S.S.; Lo, K.C.; Lui, Y.M.; Jahan, A.S.; Wu, A.L.; Ke, Y.H.; Law, C.Y.; Sze, K.H.; et al. Lipid mediators of inflammation as novel plasma biomarkers to identify patients with bacteremia. J. Infect. 2015, 70, 433–444. [Google Scholar] [CrossRef] [PubMed]
  39. XCMS. Available online: http://www.bioconductor.org/packages/2.8/bioc/html/xcms.html (accessed on 30 November 2015).
  40. MetaboAnalyst 3.0—A Comprehersive Tool Suite Dor Metabolomic Data Analysis. Available online: http://www.metaboanalyst.ca/ (accessed on 1 December 2015).
  41. Indiveri, C.; Iacobazzi, V.; Tonazzi, A.; Giangregorio, N.; Infantino, V.; Convertini, P.; Console, L.; Palmieri, F. The mitochondrial carnitine/acylcarnitine carrier: Function, structure and physiopathology. Mol. Asp. Med. 2011, 32, 223–233. [Google Scholar] [CrossRef] [PubMed]
  42. Chin, C.Y.; Monack, D.M.; Nathan, S. Genome wide transcriptome profiling of a murine acute melioidosis model reveals new insights into how Burkholderia pseudomallei overcomes host innate immunity. BMC Genom. 2010, 11, 672. [Google Scholar] [CrossRef] [PubMed]
  43. Schmerler, D.; Neugebauer, S.; Ludewig, K.; Bremer-Streck, S.; Brunkhorst, F.M.; Kiehntopf, M. Targeted metabolomics for discrimination of systemic inflammatory disorders in critically ill patients. J. Lipid Res. 2012, 53, 1369–1375. [Google Scholar] [CrossRef] [PubMed]
  44. Tepperman, B.L.; Soper, B.D. The role of phospholipase A2 in calcium-ionophore-mediated injury to rat gastric mucosal cells. Dig. Dis. Sci. 1999, 44, 494–502. [Google Scholar] [CrossRef] [PubMed]
  45. Facciotti, F.; Ramanjaneyulu, G.S.; Lepore, M.; Sansano, S.; Cavallari, M.; Kistowska, M.; Forss-Petter, S.; Ni, G.; Colone, A.; Singhal, A.; et al. Peroxisome-derived lipids are self antigens that stimulate invariant natural killer T cells in the thymus. Nat. Immunol. 2012, 13, 474–480. [Google Scholar] [CrossRef] [PubMed]
  46. Yamamoto, N.; Ngwenya, B.Z. Activation of mouse peritoneal macrophages by lysophospholipids and ether derivatives of neutral lipids and phospholipids. Cancer Res. 1987, 47, 2008–2013. [Google Scholar] [PubMed]
  47. Endo, S.; Inada, K.; Nakae, H.; Takakuwa, T.; Yamada, Y.; Suzuki, T.; Taniguchi, S.; Yoshida, M.; Ogawa, M.; Teraoka, H. Plasma levels of type II phospholipase A2 and cytokines in patients with sepsis. Res. Commun. Mol. Pathol. Pharmacol. 1995, 90, 413–421. [Google Scholar] [PubMed]
  48. Gonzalez, E.; van Liempd, S.; Conde-Vancells, J.; Gutierrez-de Juan, V.; Perez-Cormenzana, M.; Mayo, R.; Berisa, A.; Alonso, C.; Marquez, C.A.; Barr, J.; et al. Serum uplc-ms/ms metabolic profiling in an experimental model for acute-liver injury reveals potential biomarkers for hepatotoxicity. Metabolomics 2012, 8, 997–1011. [Google Scholar] [CrossRef] [PubMed]
  49. Yang, J.; Wang, H.; Xu, W.; Hao, D.; Du, L.; Zhao, X.; Sun, C. Metabolomic analysis of rat plasma following chronic low-dose exposure to dichlorvos. Hum. Exp. Toxicol. 2013, 32, 196–205. [Google Scholar] [CrossRef] [PubMed]
  50. Luberto, C.; Stonehouse, M.J.; Collins, E.A.; Marchesini, N.; El-Bawab, S.; Vasil, A.I.; Vasil, M.L.; Hannun, Y.A. Purification, characterization, and identification of a sphingomyelin synthase from Pseudomonas aeruginosa. PlcH is a multifunctional enzyme. J. Biol. Chem. 2003, 278, 32733–32743. [Google Scholar] [CrossRef] [PubMed]
  51. Korbsrisate, S.; Tomaras, A.P.; Damnin, S.; Ckumdee, J.; Srinon, V.; Lengwehasatit, I.; Vasil, M.L.; Suparak, S. Characterization of two distinct phospholipase c enzymes from Burkholderia pseudomallei. Microbiology 2007, 153, 1907–1915. [Google Scholar] [CrossRef] [PubMed]
  52. Tuanyok, A.; Tom, M.; Dunbar, J.; Woods, D.E. Genome-wide expression analysis of Burkholderia pseudomallei infection in a hamster model of acute melioidosis. Infect. Immun. 2006, 74, 5465–5476. [Google Scholar] [CrossRef] [PubMed]
  53. Theodorou, I.G.; Botelho, D.; Schwander, S.; Zhang, J.; Chung, K.F.; Tetley, T.D.; Shaffer, M.S.; Gow, A.; Ryan, M.P.; Porter, A.E. Static and dynamic microscopy of the chemical stability and aggregation state of silver nanowires in components of murine pulmonary surfactant. Environ. Sci. Technol. 2015, 49, 8048–8056. [Google Scholar] [CrossRef] [PubMed]
  54. Creuwels, L.A.; van Golde, L.M.; Haagsman, H.P. The pulmonary surfactant system: Biochemical and clinical aspects. Lung 1997, 175, 1–39. [Google Scholar] [CrossRef] [PubMed]
  55. Lu, R.; Popov, V.; Patel, J.; Eaves-Pyles, T. Burkholderia mallei and Burkholderia pseudomallei stimulate differential inflammatory responses from human alveolar type II cells (ATII) and macrophages. Front. Cell. Infect. Microbiol. 2012, 2, 165. [Google Scholar] [CrossRef] [PubMed]
  56. The R Project for Statistical Computing. Available online: https://www.r-project.org/ (accessed on 30 November 2015).
  57. Basanta, M.; Jarvis, R.M.; Xu, Y.; Blackburn, G.; Tal-Singer, R.; Woodcock, A.; Singh, D.; Goodacre, R.; Thomas, C.L.; Fowler, S.J. Non-invasive metabolomic analysis of breath using differential mobility spectrometry in patients with chronic obstructive pulmonary disease and healthy smokers. Analyst 2010, 135, 315–320. [Google Scholar] [CrossRef] [PubMed]
  58. Xia, J.; Broadhurst, D.I.; Wilson, M.; Wishart, D.S. Translational biomarker discovery in clinical metabolomics: an introductory tutorial. Metabolomics 2013, 9, 280–299. [Google Scholar] [CrossRef] [PubMed]
  59. Scripps Center for Meatbolomics – METLIN. Available online: http://metlin.scripps.edu/ (accessed on 1 December 2015).
  60. The Human Metabolome Database. Available online: http://www.hmdb.ca/ (accessed on 1 December 2015).
  61. Massbank. Available online: http://www.massbank.jp/ (accessed on 1 December 2015).
  62. LIPID Metabolites and Pathways Strategy (LIPID MAPS®) – Lipidomics Gateway. Available online: http://www.lipidmaps.org/ (accessed on 1 December 2015).
  63. KEGG: Kyoto Encyclopedia of Genes and Genomes. Available online: http://www.genome.jp/kegg/ (accessed on 1 December 2015).
Figure 1. Principal component analysis (PCA) score plot in positive mode based on human plasma of 22 melioidosis, 24 bacteremia and 30 controls without active infection. The PCA score plots showed that samples from melioidosis patients, bacteremia patients and controls without active infection were clustered separately.
Figure 1. Principal component analysis (PCA) score plot in positive mode based on human plasma of 22 melioidosis, 24 bacteremia and 30 controls without active infection. The PCA score plots showed that samples from melioidosis patients, bacteremia patients and controls without active infection were clustered separately.
Ijms 17 00307 g001
Figure 2. Volcano Plot, using fold-change (FC) >1.5 and p-values cut-off <0.05 using Student’s t-test. The statistical analyses were performed for comparison between melioidosis and bacteremia patients in (A) positive mode as well as melioidosis patients and controls without active infection in (B) positive mode.
Figure 2. Volcano Plot, using fold-change (FC) >1.5 and p-values cut-off <0.05 using Student’s t-test. The statistical analyses were performed for comparison between melioidosis and bacteremia patients in (A) positive mode as well as melioidosis patients and controls without active infection in (B) positive mode.
Ijms 17 00307 g002
Figure 3. MS/MS mass spectra and predicted structures with expected fragmentation profiles of the 12 biomarkers in melioidosis patient plasma: (A) l-octanoylcarnitine; (B) decanoylcarnitine; (C) dodecanoylcarnitine; (D) lysophosphatidylethanolamine (LysoPE)(16:0/0:0); (E) LysoPE(18:0/0:0); (F) phosphatidylcholine PC(16:0/16:0); (G) LysoPE(0:0/18:0) (H) l-hexanoylcarnitine; (I) sphingomyelins SM(d16:1/16:0); (J) 2-decenoylcarnitine; (K) SM(d18:2/16:0); and (L) trans-2-dodecenoylcarnitine with or without comparison to commercially available standards.
Figure 3. MS/MS mass spectra and predicted structures with expected fragmentation profiles of the 12 biomarkers in melioidosis patient plasma: (A) l-octanoylcarnitine; (B) decanoylcarnitine; (C) dodecanoylcarnitine; (D) lysophosphatidylethanolamine (LysoPE)(16:0/0:0); (E) LysoPE(18:0/0:0); (F) phosphatidylcholine PC(16:0/16:0); (G) LysoPE(0:0/18:0) (H) l-hexanoylcarnitine; (I) sphingomyelins SM(d16:1/16:0); (J) 2-decenoylcarnitine; (K) SM(d18:2/16:0); and (L) trans-2-dodecenoylcarnitine with or without comparison to commercially available standards.
Ijms 17 00307 g003aIjms 17 00307 g003bIjms 17 00307 g003cIjms 17 00307 g003dIjms 17 00307 g003e
Figure 4. Box-and-whiskers plots representing relative abundance of: (A) l-hexanoylcarnitine; (B) l-octanoylcarnitine; (C) 2-decenoylcarnitine; (D) decanoylcarnitine; (E) trans-2-dodecenoylcarnitine; (F) dodecanoylcarnitine; (G) LysoPE(16:0/0:0); (H) LysoPE(0:0/18:0); (I) LysoPE(18:0/0:0); (J) SM(d16:1/16:0); (K) SM(d18:2/16:0); and (L) PC(16:0/16:0) in plasma of melioidosis patients, bacteremia patients and controls without active infections. The relative abundance of each metabolite in plasma of melioidosis patients was significantly higher than the other two groups using Student’s t-test (p-value < 0.01).
Figure 4. Box-and-whiskers plots representing relative abundance of: (A) l-hexanoylcarnitine; (B) l-octanoylcarnitine; (C) 2-decenoylcarnitine; (D) decanoylcarnitine; (E) trans-2-dodecenoylcarnitine; (F) dodecanoylcarnitine; (G) LysoPE(16:0/0:0); (H) LysoPE(0:0/18:0); (I) LysoPE(18:0/0:0); (J) SM(d16:1/16:0); (K) SM(d18:2/16:0); and (L) PC(16:0/16:0) in plasma of melioidosis patients, bacteremia patients and controls without active infections. The relative abundance of each metabolite in plasma of melioidosis patients was significantly higher than the other two groups using Student’s t-test (p-value < 0.01).
Ijms 17 00307 g004
Figure 5. Receiver operating characteristic (ROC) curves for SM(d18:2/16:0) with Area Under the Curves (AUC) with 95% confidence interval, when compared (A) between melioidosis and other bacteremia patients; and (B) between melioidosis patients and controls without active infection.
Figure 5. Receiver operating characteristic (ROC) curves for SM(d18:2/16:0) with Area Under the Curves (AUC) with 95% confidence interval, when compared (A) between melioidosis and other bacteremia patients; and (B) between melioidosis patients and controls without active infection.
Ijms 17 00307 g005
Table 1. Plasma metabolites with higher levels in melioidosis patients compared to bacterimia patients and controls.
Table 1. Plasma metabolites with higher levels in melioidosis patients compared to bacterimia patients and controls.
CompoundExperimental Mass, m/zIonRetention Time (min)MS/MS Fragment MassesElemental CompositionMetabolite Class
l-Hexanoylcarnitine260.1842[M + H]+4.4760.0808, 85.0285, 99.0803, 144.1018, 201.1117C13H25NO4acylcarnitine
l-Octanoylcarnitine288.2157[M + H]+6.7860.0805, 85.0283, 127.1110, 144.1019, 229.1438C15H29NO4acylcarnitine
2-Decenoylcarnitine314.2326[M + H]+7.8460.0806, 85.0281, 144.1015, 153.1257, 255.1591C17H31NO4acylcarnitine
Decanoylcarnitine316.2476[M + H]+8.6260.0806, 85.0284, 144.1017, 155.1424, 257.1748C17H33NO4acylcarnitine
Trans-2-dodecenoylcarnitine342.2636[M + H]+9.3860.0807, 85.0283, 144.1019, 181.1584, 283.1880C19H35NO4acylcarnitine
Dodecanoylcarnitine344.2775[M + H]+10.2060.0806, 85.0284, 144.1008, 183.1735, 285.2088C19H37NO4acylcarnitine
LysoPE(16:0/0:0)454.2934[M + H]+14.6444.0496, 62.0598, 216.0642, 239.2361,
257.2530, 313.2740, 393.2423, 436.2774
C21H44NO7Plysophosphatidylethanolamine
LysoPE(0:0/18:0)482.3244[M + H]+16.9444.0494, 216.0628, 267.2644, 285.2747, 341.3060C23H48NO7Plysophosphatidylethanolamine
LysoPE(18:0/0:0)482.3251[M + H]+17.6144.0497, 62.0600, 216.0618, 267.2672,
285.2777, 341.3058, 421.2718, 464.3112
C23H48NO7Plysophosphatidylethanolamine
SM(d16:1/16:0)675.5444[M + H]+27.8860.0808, 104.1072, 184.0735, 236.2355C37H75N2O6Psphingomyelins
SM(d18:2/16:0)701.5605[M + H]+28.6460.0802, 104.1068, 184.0736, 262.2575, 683.5484C39H77N2O6Psphingomyelins
PC(16:0/16:0)734.5616[M + H]+30.5860.0801, 104.1060, 184.0727, 478.3251, 496.3353C40H80NO8Pphosphatidylcholine
Table 2. The area-under-receiver operating characteristic curve (AUC), sensitivity and specificity for receiver operating characteristic (ROC) curves calculated at optimal cutoff as well as p-value and Fold-change for the twelve significant metabolites.
Table 2. The area-under-receiver operating characteristic curve (AUC), sensitivity and specificity for receiver operating characteristic (ROC) curves calculated at optimal cutoff as well as p-value and Fold-change for the twelve significant metabolites.
Significant MetabolitesMelioidosis vs. BacteremiaMelioidosis vs. Control without Active Infections
AUC a95% CI b Sensitivity (%)Specificity (%)p-Value c Fold-ChangeAUC a95% CI bSensitivity (%)Specificity (%)p-Value cFold-Change
l-Hexanoylcarnitine0.6650.497–0.80581.850.01.12 × 10−22.32 ↑d0.8490.728–0.95983.381.81.10 × 10−33.33 ↑d
l-Octanoylcarnitine0.8560.710–0.96486.479.26.34 × 10−43.49 ↑0.8270.693–0.94283.386.42.90 × 10−32.76 ↑
2-Decenoylcarnitine0.8390.700–0.93777.379.24.77 × 10−33.63 ↑0.8290.682–0.95290.081.81.90 × 10−33.56 ↑
Decanoylcarnitine0.8500.728–0.97277.387.51.46 × 10−43.53 ↑0.8210.680–0.94983.377.31.45 × 10−32.66 ↑
Trans-2-dodecenoylcarnitine0.8500.736–0.96577.379.21.20 × 10−52.40 ↑0.8860.776–0.96590.072.75.43 × 10−82.88 ↑
Dodecanoylcarnitine0.8220.700–0.93572.791.73.10 × 10−52.46 ↑0.7410.583–0.87776.772.79.60 × 10−31.66 ↑
LysoPE(16:0/0:0)0.9790.947–1.00090.995.81.51 × 10−105.20 ↑0.8120.672–0.93373.381.82.56 × 10−62.01 ↑
LysoPE(0:0/18:0)0.9940.982–1.00095.5100.06.08 × 10−107.51 ↑0.8190.658–0.92690.077.31.96 × 10−62.23 ↑
LysoPE(18:0/0:0)0.9980.993–1.000100.095.81.77 × 10−116.09 ↑0.8560.725–0.95890.077.31.41 × 10−72.16 ↑
SM(d16:1/16:0)0.9680.927–1.00090.991.77.19 × 10−103.41 ↑0.8840.783–0.98383.377.32.67 × 10−82.16 ↑
SM(d18:2/16:0)0.9980.993–1.000100.091.71.88 × 10−123.32 ↑1.0001.000–1.00096.7100.01.28 × 10−132.65 ↑
PC(16:0/16:0)0.8350.695–0.97677.387.57.01 × 10−59.64 ↑0.8700.724–0.98993.381.82.77 × 10−621.72 ↑
a AUC = area-under-receiver operating characteristic curve; b CI = confidence interval; c All p-values were calculated using Student’s t-test; d↑ = Higher level comparing melioidosis to the respective groups.

Share and Cite

MDPI and ACS Style

Lau, S.K.P.; Lee, K.-C.; Lo, G.C.S.; Ding, V.S.Y.; Chow, W.-N.; Ke, T.Y.H.; Curreem, S.O.T.; To, K.K.W.; Ho, D.T.Y.; Sridhar, S.; et al. Metabolomic Profiling of Plasma from Melioidosis Patients Using UHPLC-QTOF MS Reveals Novel Biomarkers for Diagnosis. Int. J. Mol. Sci. 2016, 17, 307. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms17030307

AMA Style

Lau SKP, Lee K-C, Lo GCS, Ding VSY, Chow W-N, Ke TYH, Curreem SOT, To KKW, Ho DTY, Sridhar S, et al. Metabolomic Profiling of Plasma from Melioidosis Patients Using UHPLC-QTOF MS Reveals Novel Biomarkers for Diagnosis. International Journal of Molecular Sciences. 2016; 17(3):307. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms17030307

Chicago/Turabian Style

Lau, Susanna K. P., Kim-Chung Lee, George C. S. Lo, Vanessa S. Y. Ding, Wang-Ngai Chow, Tony Y. H. Ke, Shirly O. T. Curreem, Kelvin K. W. To, Deborah T. Y. Ho, Siddharth Sridhar, and et al. 2016. "Metabolomic Profiling of Plasma from Melioidosis Patients Using UHPLC-QTOF MS Reveals Novel Biomarkers for Diagnosis" International Journal of Molecular Sciences 17, no. 3: 307. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms17030307

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop