Next Article in Journal
Fatty Acids and Membrane Lipidomics in Oncology: A Cross-Road of Nutritional, Signaling and Metabolic Pathways
Previous Article in Journal
Occurrence and Exposure Assessment of Mycotoxins in Ready-to-Eat Tree Nut Products through Ultra-High Performance Liquid Chromatography Coupled with High Resolution Q-Orbitrap Mass Spectrometry
Previous Article in Special Issue
MetPC: Metabolite Pipeline Consisting of Metabolite Identification and Biomarker Discovery Under the Control of Two-Dimensional FDR
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrating Two-Dimensional Gas and Liquid Chromatography-Mass Spectrometry for Untargeted Colorectal Cancer Metabolomics: A Proof-of-Principle Study

1
Department of Chemistry, University of Louisville, Louisville, KY 40292, USA
2
Department of Oncology, Wayne State University School of Medicine, Detroit, MI 48201, USA
3
Biostatistics Core, Karmanos Cancer Institute, Wayne State University, Detroit, MI 48201, USA
4
Department of Pathology, Wayne State University School of Medicine, Detroit, MI 48201, USA
*
Author to whom correspondence should be addressed.
Metabolites 2020, 10(9), 343; https://doi.org/10.3390/metabo10090343
Submission received: 1 August 2020 / Revised: 21 August 2020 / Accepted: 21 August 2020 / Published: 25 August 2020
(This article belongs to the Special Issue High-Throughput Metabolomics)

Abstract

:
Untargeted metabolomics is expected to lead to a better mechanistic understanding of diseases and thus applications of precision medicine and personalized intervention. To further increase metabolite coverage and achieve high accuracy of metabolite quantification, the present proof-of-principle study was to explore the applicability of integration of two-dimensional gas and liquid chromatography-mass spectrometry (GC × GC-MS and 2DLC-MS) platforms to characterizing circulating polar metabolome extracted from plasma collected from 29 individuals with colorectal cancer in comparison with 29 who remained cancer-free. After adjustment of multiple comparisons, 20 metabolites were found to be up-regulated and 8 metabolites were found to be down-regulated, which pointed to the dysregulation in energy metabolism and protein synthesis. While integrating the GC × GC-MS and 2DLC-MS data can dramatically increase the metabolite coverage, this study had a limitation in analyzing the non-polar metabolites. Given the small sample size, these results need to be validated with a larger sample size and with samples collected prior to diagnostic and treatment. Nevertheless, this proof-of-principle study demonstrates the potential applicability of integration of these advanced analytical platforms to improve discrimination between colorectal cancer cases and controls based on metabolite profiles in future studies.

1. Introduction

Untargeted metabolomics has been increasingly employed in the past decade for biomarker discoveries for various diseases/conditions as well as for metabolome-wide association studies [1,2]. This is expected to lead to a mechanistic understanding of diseases and thus applications of precision medicine and personalized intervention. Currently, multiple analytical techniques are available for untargeted metabolomics, and they are broadly grouped into gas chromatography (GC) or liquid chromatography (LC) mass spectrometry (MS)-based and nuclear magnetic resonance spectroscopy (NMR)-based studies [1,3,4]. Comparing with NMR, GC-MS and LC-MS are more sensitive and provide much-increased metabolite coverage. However, each of these two methods still has limited resolving power for analysis of complex samples such as metabolites in a biological sample, and therefore none of them is able to analyze all metabolites. Furthermore, it has been reported that metabolites detected by LC-MS and GC-MS are different and that the metabolite overlap and correlations between platforms are modest [5,6,7].
The majority of the untargeted metabolomics studies, thus far, were performed using GC-MS and/or LC-MS, where one-dimensional GC or LC was coupled with MS [8,9]. With the vast complexity and diverse chemical characteristics of metabolites as well as large concentration difference between metabolites in the metabolome, two-dimensional LC or GC coupled with MS have been developed and respectively used in untargeted metabolomics to achieve higher metabolite coverage [9,10]. To further increase metabolite coverage and achieve high accuracy of metabolite quantification, we have developed a method by analyzing the same metabolite sample on comprehensive two-dimensional gas chromatography-mass spectrometry (GC × GC-MS) and two-dimensional liquid chromatography-mass spectrometry (2DLC-MS). We then integrate their data for metabolite identification and quantification [7,11,12]. Such analysis does not only dramatically increase the metabolite coverage but also provide accurate metabolite quantification and reliable metabolic pathway assignment. These advanced techniques should be particularly useful in dealing with complex biological samples [13], such as plasma/serum, which contains not only human endogenous metabolites but also environmental metabolites, plant metabolites, pharmaceutical metabolites, and microbial metabolites in a wide range of concentrations [1]. Yet, the introduction of these advanced two-dimensional approaches to clinical samples has not progressed as rapidly as expected.
Altered energy metabolism, exampled by aerobic glycolysis, has been recognized recently as an emerging hallmark of cancer [14]. Furthermore, metabolic syndrome, consisting of a cluster of abnormal lipid, glucose, and energy metabolisms, predisposes individuals to various cancers [15,16,17]. Colorectal cancer (CRC) is the third most commonly diagnosed cancer in both genders and the second leading cause of cancer deaths in the USA [18]. Although metabolic syndrome is an established risk factor for CRC [15,16,17], there have been few tumor markers for CRC with good diagnostic, predictive, or prognostic values. Carcinoembryonic antigen (CEA) has been routinely used to detect CRC recurrence after surgery but has insufficient sensitivity to be used alone [19]. Thus, new sets of biomarkers would greatly facilitate CRC prevention and patient care. The objectives of the present proof-of-principal study was to explore the applicability of integration of GC × GC-MS and 2DLC-MS platforms to characterizing circulating metabolome in individuals who suffered CRC in comparison with those who remained cancer-free, both limited to non-smokers to minimize influences of tobacco metabolites.

2. Results

Our two-dimensional gas chromatography was configured in comprehensive mode, i.e., multiple fractions are collected from the first-dimension column with a modulation period of 2.00 s and each fraction was then subjected to the second-dimension column for further separation. Our 2DLC was configured in parallel mode [20]. Specifically, a sample was first delivered to two sample loops during sample loading. The two sample aliquots were then simultaneously injected onto a dual column setup, and parallel separations were performed on the HILIC and RPC columns, respectively. The eluates from the two columns were then merged and subjected to a mass spectrometer for further analysis. This strategy is simple yet effective for coupling HILIC and RPC to decrease analysis time and increasing throughput. Furthermore, the parallel 2DLC-MS configuration allows the use of two long columns and gradient time to increase separation power and does not suffer from the peak partition and solvent miscibility and solvent strength issues.
GC × GC-MS detected 52 to 126 annotated peaks (i.e., identified metabolites) from each sample with a median of 110, while 2DLC-MS (+) detected 116 to 145 annotated peaks with a median of 127 and 2DLC-MS (−) detected 91 to 137 annotated peaks with a median of 118. When the numbers of annotated peaks were compared by group, the case group had slightly higher numbers than or similar numbers with the control group, i.e., GC × GC-MS (median 112 vs.107), 2DLC-MS (+) (median 127 vs. 127), and 2DLC-MS (−) (median 118 vs.119), respectively, for the cases and controls. It should be noted that each annotated peak in GC × GC-MS refers to a unique metabolite, while two annotated peaks in 2DLC-MS (+) or 2DLC-MS (−) may represent the same metabolite due to parallel 2DLC configuration.
To study the metabolite identification and abundance changes between groups, we filtered out annotated peaks in the alignment table that presented in < 50% of the samples in each group. After this filtering, the GC × GC-MS dataset had 104 annotated peaks, 2DLC-MS (+) dataset had 112 annotated peaks, and 2DLC-MS (−) dataset had 106 annotated peaks. Figure 1 shows a Venn diagram among unique metabolites identified from multiple platforms after filtering. A total of 16 metabolites were commonly detected by all three platforms, while 19 metabolites by both GC × GC-MS and 2DLC-MS (+), 25 metabolites by both GC × GC-MS and 2DLC-MS (−), and 26 metabolites by both 2DLC-MS (+) and 2DLC-MS (−). Figure 2 shows the results of PLS-DA analysis using the merged data of the three platforms, which confirm that the samples are discriminated by group.
Table 1 shows the results of statistical significance tests, i.e., metabolites with significant abundance changes between cases and controls. After multiple comparison adjustment, 9 metabolites with significant abundance differences were detected by 2DLC-MS (+), and 10 metabolites detected by 2DLC-MS (−) and 9 metabolites detected by GC × GC-MS, with one overlapped metabolite malate between 2DLC-MS (−) and GC × GC-MS. Among these 28 metabolites with different relative abundances, 20 were up-regulated in the cases while 8 were down-regulated. Among the 20 up-regulated metabolites, more than 50% of fold changes were detected for 6 metabolites, including 3-hydroxybutyric acid (2.12), glutamic acid (2.00), beta-alanine (1.76), tricine (1.60), pyroglutamic acid (1.55), malate (1.51), and citric acid (1.51). Among the 8 down-regulated metabolites, the top 5 abundance decline in the cases was observed for propylene glycol (0.50), alanine (0.52), testosterone sulfate (0.54), pinolenic acid (0.56), and piperine (0.63).
It should be noted that 4 metabolites with significance abundance changes were actually detected in multiple platforms. For example, glutamic acid was detected by all three platforms. While its abundance has significant difference between the cases and the controls in the GC × GC-MS data with p = 0.00258 and fold change = 2.00, its abundance changes detected by 2DLC-MS (+) and 2DLC-MS (−) are not significant (p = 0.136 and fold change = 1.31 in 2DLC-MS (+); p = 0.257 and fold change = 1.19 in 2DLC-MS (−)). However, its abundance change has the same direction, i.e., all up-regulated in the cases. Manual data analysis shows that such discrepancy was mainly caused by the low instrument response (Figure 3), as we reported before [7,11].
Combining these three datasets, we further explored the metabolic pathway that was altered by disease status. A total of 13 pathways were detected with a p < 0.05. The top five pathways affected by disease history are aminoacyl-tRNA biosynthesis (6 hits), tricarboxylic acid (TCA) cycle (3 hits), glyoxylate and dicarboxylate metabolism (4 hits), beta-alanine metabolism (3 hits), and alanine, aspartate, and glutamate metabolism (4 hits) (Table 2 and Figure 4).
Figure 5 depicts the results of validating identified 28 potential metabolite biomarkers (Table 1) using the receiver operating characteristic (ROC) analysis. The metabolites selected for ROC analysis were 4-methyl-2-oxovaleric acid, aspartic acid, beta-alanine, citric acid, glycine, lysophosphatidic acid, neopentyl glycol, oleamide, propylene glycol, testosterone sulfate, and uridine. While 1SE and OPT were respectively used to select an optimal λ in LASSO logistic regression for the development of a prediction model to distinguish CRC cases from control groups, both tuning methods picked up the same set of potential biomarkers as the final list of markers. The area under the curves (AUC) of ROC is 0.954 for 1SE turned (Figure 5A) and 0.970 for OPT tuned LASSO logistic regression (Figure 5B).

3. Discussion

Following our previous study focused on volatile organic compounds (VOC) [21], the primary aim of this study was to apply multiple two-dimensional separation techniques with mass spectrometry for metabolite biomarker discovery from complex biological samples by untargeted metabolomics. In the past several years, a number of global metabolomic studies on CRC have been published. These studies utilized diverse biological specimens, including blood (plasma/serum) as the most common type [22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39], followed by tissue [27,39,40,41,42], feces [39,40,43], and urine [44,45], while others aimed specifically to detect volatile metabolites in breath and in other types of biospecimens [21,46]. To our knowledge, only two of these were prospective studies where biospecimens were collected before cancer development [26,34], aiming to identify etiological pathways. The rest and ours were based on samples collected after the cancer diagnosis, which was more focused on biomarker discovery. A minority of these studies employed more than one untargeted metabolomic platforms [24,26,31,32,34,43,45], but none of these incorporated multiple two-dimensional separations.
Because metabolic profiles are highly dependent on the choice of the matrix [27], the comparison across different types of biospecimens is not informative. Blood is the most widely used clinical specimens. Hashim et al. [23] recently reviewed the results of blood-based global metabolomic studies for CRC. The authors identified nine studies and found that the number of annotated metabolites in each study vastly varied from 14 [25] to 447 [34], not necessarily in a platform-dependent manner. Each study applied heterogeneous criteria for differential abundance between cases and controls, yielding the number of differential metabolites ranging from 5 [30] to 72 [31]. More recently published studies identified nearly or over 600 named metabolites [26,28], while another recent paper reported a vast number of differentially present metabolites, i.e., 693, in plasma of CRC cases compared with polyp patients [27]. It should be noted that the number of annotated metabolites is affected by many factors, including sample treatments, bioanalytical platforms used for analysis, and the methods of data analysis. Earlier studies with a larger number of annotated metabolites generally used an additional extraction and/or separation technique specifically for lipid metabolites, which are known to be less susceptible to degradation over long-term storage [47].
What was apparent in the aforementioned review and other recent papers on circulating CRC-associated metabolites is that very few metabolites have been detected consistently across the studies. Even if a metabolite was reported in more than one study, often the direction association was different [23]. For instance, there were no lipid metabolites commonly detected in the same direction in recent studies [26,28,29]. The most consistently observed in the same direction as well as by more than one platform were pyruvic acid (up-regulated) and tryptophan (down-regulated) [23]. However, in our study population, neither pyruvic acid nor tryptophan showed significant abundance difference between the cases and the controls, being the fold changes of these metabolites close to unity (0.94–0.99, data not shown).
On the other hand, some of the metabolites differentially present by CRC status in our study were noted by others. Specifically, increased relative abundances of glycine and glutamic acid have been reported in studies [25,30,33,35,40,44]. As in our study, three of these studies [35,40,41] detected increased abundances of both metabolites simultaneously, including two tissue-based studies [40,41] analyzing metabolites directly from CRC but using different analytical platforms (capillary electrophoresis (CE)-TOFMS, NMR, and GC-MS). This indicates that the source of these metabolites is more likely to be cancer itself, rather systemic reactions to the disease. The reported fold changes, if any, have been rather modest (<2.0), consistent with our results. In addition, an increased relative abundance of malate was reported by three serum-based studies with a similar fold change to ours (×1.3) [21,32,36]. Furthermore, increased abundances of uridine and isocitrate were respectively reported in two studies [43,48], and four other studies found a decreased relative abundance of leucine [25,29,30,36].
Based on metabolic pathway analysis, we found five biological pathways were affected by disease history. Among these, alanine, aspartate, glutamate, and β-alanine metabolism pathways are also reported to be dysregulated in both CRC tissue and blood by others [25,27,41,43]. These two metabolic pathways and the TCA cycle are known to be central to fulfill increased energy consumption and increased anabolic activities to sustain cellular growth and proliferation [49]. Glucose and glutamine are the main sources used to maintain active essential metabolic pathways, such as glycolysis and anaplerotic flux of the tricarboxylic acid (TCA) cycle. Dysregulation of the TCA cycle and glutamine addiction are characteristic features of glutamine metabolism in cancers [48]. Glutamate is formed through hydrolysis of glutamine, and it can be further converted to α-ketoglutarate and other TCA intermediates [40,49].
Up-regulation of β-alanine metabolism pathway has been noted in both CRC tissue and plasma by Wang et al. [27]. Besides CRC, increased β-alanine activity was positively associated with an aggressive type of estrogen receptor-negative breast cancer in both tumor tissues and established cell lines [50,51]. β-alanine is crucial for acetyl-CoA synthesis that is vital in TCA cycle and biomass production such as the synthesis of fatty acids, cholesterols, and acetylcholines, fostering cellular metabolism favorable for tumor cell growth [52]. Epithelial-mesenchymal transition induces cancer stem cell characteristics and promotes tumor invasiveness, which was associated with the metabolic reprogramming, including beta-alanine upregulation, and increased β-alanine levels in breast cancer tissue were predictive of overall survival [53]. It is also worth noting that dysregulation of β-alanine metabolism is a common feature of obesity, a known risk factor for CRC [54].
Glyoxylate and dicarboxylate metabolic pathway comprises an important subset of reactions involved in biosynthesis of carbohydrates from fatty acids or two-carbon precursors, which enter the system as acetyl-CoA and thus is closely related to several pathways discussed above [55]. Despite the lack of reports from CRC metabolic steadies, its dysregulation has been described in tumor tissue, peripheral blood, and urine from a variety of cancer-bearing humans and animals including gastric cancer, hepatocellular carcinoma, and myeloid leukemia [55,56,57]. This was confirmed in pathway analysis based on gene expression of aggressive breast cancer [58]. Thus, dysregulation of this pathway may represent a common metabolic profile of advanced cancer.
Our data also pointed to the dysregulation of the aminoacyl-tRNA biosynthesis pathway CRC patients. Aminoacyl-tRNA synthetases (ARSs) carry out the first step of protein translation by catalyzing the litigation of amino acids to their cognate tRNAs. Aberrant expression, mislocalization, and variant formation of ARSs have been recently observed in various cancer cells [59,60]. However, it has not been clear whether increased expression of ARSs in cancer tissue/cells results from the increased demand for protein synthesis in cancer cells or whether ARSs are drivers of cell transformation and growth. Overexpression of specific types of ARSs has also been noted. Particularly, an increase in the catalytic activity of methionyl-tRNA synthetase (MRS) that is required for translation initiation was reported in human colon cancer [61] and experimental overexpression of methionyl-tRNA in breast epithelial cell lines resulted in increased metabolic activity and cell proliferation [62]. Leucyl-tRNA synthetase overexpression has been linked to activation of mammalian target of rapamycin complex 1 (mTORC1)-signaling pathway that controls critical cellular processes, including cell growth, metabolism, and autophagy, and thus contributes to carcinogenesis [63]. In addition, leucyl-tRNA-derived small RNA was found essential for cell viability and tumor growth mediating ribosome biogenesis [64]. Furthermore, it has become clear that ARSs possess versatile non-canonical functions, interacting with diverse regulatory networks and fine-tuning post-transcriptional gene expression, which is relevant to various pathological conditions, including cancer [60,61,62,63].
While metabolomic profiling has been informative in identifying new therapeutic targets, the clinical utility of selected metabolomic markers has been tested in early detection of colorectal cancer [33,36]. Although most single metabolites showed poor to moderate values of area under the curves (AUC), high AUC values (>0.90) were achieved using a panel of selected metabolites in our and other studies. However, actual clinical translational values of this approach are still questionable. Metabolite abundance obtained from an untargeted metabolomic study is generally normalized relative values. To be used as a clinical test, cut-offs based on absolute concentrations are necessary. More importantly, little has been known about changes in these endogenous metabolites by physiological or pathological conditions other than colorectal cancer. It is possible that identified metabolites are merely markers for underlying obesity, inflammation, or aging [65,66,67]. Furthermore, since circulating levels of the metabolites are expected to be much lower than those in cancer tissue are and since these endogenous metabolites are not tissue-specific, high sensitivity will be required to detect changes in peripheral blood arising from the small early stage of cancer.
The strengths of current study include the controls derived from the same general population that gave rise to the cases, availability of the information of detailed personal characteristics, and integration of multiple advanced separation techniques in global metabolomics. We acknowledge limitations in our study. Being an early stage of the study, the sample size was small to obtain estimates with high precision. Besides, as biological samples used in this study were collected after cancer treatment/resection from the cases, the direct effect of the presence of the primary tumor was deemed to have faded, although underlying metabolic deviations that fostered cancer development were expected to remain. In fact, resection of CRC does not normalize metabolic shifts seen in breath VOCs [68] or in blood total metabolites [29]. Although the storage at −80 °C is considered appropriate for most metabolomic studies, metabolite degradation over the long-term storage is certainly possible [33]. However, it is reasonable to assume that the degradation would occur similarly over time between cases and controls. Yet, we acknowledge that the discriminating power may decline with longer storage due to the degradation of low abundance metabolites. Finally, the integration of lipid profiling, which was missing in this study, would provide a more comprehensive view of cancer-causing or cancer-induced metabolomic abnormalities

4. Materials and Methods

4.1. Materials

Sodium phosphate dibasic (Na2HPO4, purity 100.2%, Cat. No. S374-500), sodium phosphate monobasic (NaH2PO4, purity 99%, Cat. No. 389872500), and sodium hydroxide (NaOH, Cat. No. S399-500) were purchased from Fisher Scientific Inc. (Pittsburgh, PA, USA). Dichloromethane (CH2Cl2, purity ≥ 99.8%, Cat. No. 34856-1L), methoxyamine hydrochloride (purity ≥ 98%, Cat. No. 226904-5G), and a mixture containing C7–C26 n-alkanes (purity 95–99%, Cat. No. 49451U) were purchased from Sigma-Aldrich Corp. (St. Louis, MO, USA). N-methyl-N-[tert-butyldimethylsilyl trifluoroacetamide (MTBSTFA) with 1% tert-butyldimethylchlorosilane (TBDMCS) was purchased from Restek Corp. (Bellefonte, PA, USA). dH2O was prepared from a Millipore synergy system (Burlington, MA, USA).

4.2. Study Samples

The samples used in this study were plasma samples collected for a population-based CRC case-control study. It was conducted in Metropolitan Detroit during the period of 2003–2006 and designed to enroll frequency-matched cases and controls according to their gender, age, and county. The study was approved by the Institutional Review Board of Wayne State University. Details concerning the parent study were described elsewhere [69]. Blood was collected in ethylenediaminetetraacetic acid tubes at subjects’ home or other preferred locations and delivered to a processing laboratory with cooling gels where plasma samples were separated, aliquoted, and stored in −80 °C freezers. Dietary and other personal characteristics were obtained by the use of structured questionnaires. A total of 1421 subjects who were not smoking for at least preceding two years, either Caucasian or African American, donated blood samples in the parent study.
For this study, 30 cases and 30 controls were randomly selected from the original pool of the parent study, stratified by race and gender (50% white/black; 50% male/female), as well as no cancer targeted therapies (surgery, radiation, or chemotherapy). Because the freeze-thaw cycle affects the stability of metabolites, we used samples not previously thawed for this study. The mean age of both cases and control groups was equally 64 years old. Among selected, two misretrieved samples (one case and one control) were removed from the analysis, leaving 29 cases and 29 controls.

4.3. Sample Preparation

To extract polar metabolites from a plasma sample, 400 µL methanol was mixed with 100 µL plasma. After 2 min of vigorously vortex-mixing, the mixture was centrifuged at 4 °C, 14,000 rpm for 20 min. The supernatant was aliquoted into two 200 µL aliquots, and both aliquots were dried by speed vacuum at 4 °C to remove methanol followed by freeze-drying to remove water. Dried aliquots were stored at −80 °C until use. One aliquot was redissolved in 150 µL 20% acetonitrile for 2DLC-MS analysis, while the other aliquot was used for GC × GC-MS analysis after being derivatized by MTBSTFA.
For derivatization, the dried aliquot was dissolved in 30 μL methoxyamine hydrochloride (MH) solution (20 mg/mL in pyridine) and vigorously vortex-mixed for 2 min. Methoxylation was carried out at 60 °C for 1 h. Then, 30 μL MTBSTFA with 1% TBDMCS was added to the sample. Derivatization was carried out at 60 °C for 1 h. After centrifugation at 20 °C, 3000 rpm for 10 min, the supernatant was transferred into a GC vial. Methoxylation and derivatization were conducted right before GC × GC-MS analysis. For quality control, a pooled sample was prepared in the same way by mixing a small amount of each sample. The C7–C26 n-alkane mixture was dissolved in dichloromethane at a concentration of 20 µg/mL per compound and analyzed to convert retention time into retention index.

4.4. GC × GC-MS Analysis

GC × GC-MS analysis was performed on a LECO Pegasus 4D GC × GC-Time of Flight (TOF) system (LECO Corporation, Saint Joseph, MI, USA) consisting of a Gerstel MPS2 auto-sampler (GERSTEL Inc., Linthicum, MD, USA), an Agilent 6890 gas chromatography (Agilent Technologies, Santa Clara, CA, USA), a LECO two-stage cryogenic modulator, and a LECO Pegasus time-of-flight mass spectrometer with an electron ionization (EI) ion source. The first-dimension column was a J&W DB-5ms Ultra Inert column (60 m × 0.25 mm × 0.25 μm) and the second-dimension column was a J&W DB-17 column (1 m × 0.1 mm × 0.1 μm). Both columns were purchased from Agilent Technologies. The sample injection volume was 1 µL, with a split ratio of 10. The carrier gas was helium (99.999% purity) at 2.0 mL/min corrected constant flow via pressure ramps. The temperature of the inlet and the transfer line were all set to 280 °C. The primary oven temperature ramp for the first-dimension column was 60 °C for 0.5 min, then 5 °C/min to 270 °C, and hold for 15 min at 270 °C. The secondary oven temperature ramp for the second-dimension column was 10 °C warmer with the same temperature gradient as the primary oven. The modulator temperature was 15 °C warmer than the secondary oven, and modulation period was set to 2.00 s. The mass acquisition range was 29–800 Da at an acquisition rate of 200 spectra/s. The ion source temperature was 230 °C. The detector voltage was 1521 V with an electron energy of 70 eV. The acceleration voltage was turned on after a solvent delay of 645 s.

4.5. 2DLC-MS Analysis

2DLC-MS analysis was performed on a DIONEX UltiMate 3000 UHPLC system coupled with a Q Exactive HF Hybrid Quadrupole-Orbitrap mass spectrometer (Thermo Fisher Scientific, Inc., Langenselbold, Germany). A reversed-phase chromatography (RPC) column and a hydrophilic interaction chromatography (HILIC) column was configured in parallel mode in the LC system [21]. The RPC column was ACQUITY UPLC HSS T3 (150 × 2.1 mm, 1.8 μm, Waters Corporation, Milford, MA, USA) and the HILIC column was SeQuant ZIC-cHILIC (150 × 2.1 mm, 3 μm, Merck KGaA, Darmstadt, Germany) equipped with a guard column. The LC conditions were set as follows: mobile phase A was water with 0.1% formic acid for PRC and 10 mM ammonium acetate (pH adjusted to 3.25 with acetate) for HILIC. The mobile phase B was acetonitrile with 0.1% formic acid for both RPC and HILIC. The RPC gradient was 0 min, 5% B, hold for 5.0 min; 5.0 min to 6.1 min, 5% B to 15% B; 6.1 min to 10.0 min, 15% B to 60% B, hold for 2.0 min; 12.0 min to 14.0 min, 60% B to 100% B, hold for 13.0 min; 27.0 min to 27.1 min, 100% B to 5% B, hold for 12.9 min. The HILIC gradient was 0 to 5.0 min, 95% B to 35% B, hold for 1.0 min; 6.0 min to 6.1 min, 35% B to 5% B, hold for 16.9 min; 23.0 min to 23.1 min, 5% B to 95% B, hold for 16.9 min. The flow rate was 0.4 mL/min for RPC or 0.3 mL/min for HILIC. The column temperature was 40 °C for both RPC and HILIC.
Under the same LC conditions, all samples were firstly analyzed in full MS scan mode in a random order to obtain MS data for metabolite quantification. Samples from the same group were then combined as pool samples and analyzed in full MS/dd-MS2 scan mode under collision energy of 20, 40, and 60 eV, respectively. Both MS and MS/MS data were collected in negative and positive ionization modes using parameters recommended by the vendor (Supplemental Table S1).

4.6. Data Analysis

GC × GC-MS data were first analyzed by LECO’s instrument control software ChromaTOF (version 4.51.6.0) using parameters recommended by the vendor (Supplemental Table S2). Metabolites were initially assigned by mass spectral matching using the NIST/EPA/NIH mass spectral library (version 2.2) as reference. Any metabolite assignment with a similarity score below 500 (out of 1000) was discarded. The analysis results of each sample were exported as a peak list containing the information of top 10 metabolite assignments for each chromatographic peak. All peak lists were further analyzed by MetPP software for retention index (RI) matching, peak merging, and cross-sample peak list alignment [64,70]. The initially assigned metabolites were further filtered by RI matching with a threshold of p ≤ 0.001 (equivalent to an absolute RI window of ±59 RI units). To further reduce false positives, different chromatographic peaks with the same metabolite identification were removed except peaks that were verified by experiment information of authentic standards recorded in our in-house database.
The 2DLC-MS data were first converted into mzXML format by MSConvert software (Proteowizard, Palo Alto, CA, USA). The mzXML files were then subjected to MetSign software for spectrum deconvolution and cross-sample peak list alignment [71,72,73]. Metabolite identification was achieved using the 2DLC-MS/MS data of the pooled samples, as described in our previous work [7,11,12]. Briefly, MS/MS spectra of pool samples were first matched with the information of authentic standards of metabolites recorded in our in-house database, including parent ion m/z, retention time, and MS/MS spectrum similarity with thresholds of parent ion m/z variation < 5 ppm, retention time variation < 0.15 min, and MS/MS spectral similarity > 0.4. The experimental data of MS/MS spectra that do not have a match in the in-house database were further subjected to Compound Discoverer v2.0 (Thermo Fisher Scientific Inc., Langenselbold, Germany), where metabolite identification was achieved by parent ion matching with m/z variation < 5 ppm and MS/MS spectrum matching with a similarity threshold > 40 (out of 100).
After metabolite identification, data collected from three platforms, GC × GC-MS, 2DLC-MS (+), and 2DLC-MS (−) were merged and normalized by contrast-based normalization using Metsign software [71,72,73]. Partial least squares-discriminant analysis (PLS-DA) was performed to study the metabolic profile difference among groups. A pairwise two-tail t-test was carried out to determine whether a metabolite has a significant difference in abundance between two groups. Grubbs’ test was employed for outlier detection before the t-test, and the outliers were excluded for the t-test. After the t-test, the p-values were adjusted for multiple comparisons using permutation-based adjustments with 1000 bootstraps. A metabolite was considered to be significantly different in abundance between two groups if it was detected in more than 50% of samples in each group and had an adjusted p < 0.05.
We further validated identified potential metabolomic markers using the receiver operating characteristic (ROC) analysis. We used a penalized logistic regression for the development of a prediction model to distinguish CRC cases from control groups using the least absolute shrinkage and selection operator (LASSO) approach. The LASSO logistic regression has a tuning parameter. The lambda (λ), a tuning parameter of the LASSO logistic regression, was tuned using leave-one-out cross-validation. When selecting an optimal λ, we used two methods, one standard error (1SE) and minimum (OPT), based on the deviance (i.e., mean squared error). The method OPT selects the value of lambda with the smallest deviance, while the method 1SE chooses the value of lambda that minimizes the deviance plus one standard error.
Pathway analysis was performed using Metaboanalyst (version 4.0) [74]. Names of metabolites identified with significant abundance differences were uploaded and analyzed by the program. The data were matched against the KEGG pathway library for homo sapiens. The “Hypergeometric Test” and “Relative-betweenness Centrality” algorithms were, respectively, used for over-representation and pathway topology analysis.

5. Conclusions

This proof-of-principle study tested the integration of two two-dimensional GC- and LC-MS for metabolite biomarker discovery for colorectal cancer using archived human plasma specimens. After adjustment of multiple comparisons, 20 metabolites and 8 metabolites were found to be up-regulated and down-regulated, respectively. The metabolites with 50% or higher increase in cancer cases than controls included 3-hydroxybutyric acid, glutamic acid, beta-alanine, tricine, pyroglutamic acid, malate, and citric acid. These metabolites with significant abundance differences also pointed to the dysregulation in energy metabolism and protein synthesis. While integrating the GC × GC-MS and 2DLC-MS data could dramatically increase the metabolite coverage, we only analyzed the polar metabolites in this study. Given this limitation, caution should be excised not to over-interpret our results. Furthermore, due to the small sample size, these results need to be validated with a large sample size and with samples collected prior to diagnostic and treatment. Nevertheless, this proof-of-principle study demonstrates the potential applicability of integration of these advanced analytical platforms to improve discrimination between colorectal cancer cases and controls based on metabolite profiles in future studies.

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-1989/10/9/343/s1, Table S1: Mass spectrometer parameters for 2DLC-MS analysis. Table S2: ChromaTOF software parameters for GC × GC-MS analysis.

Author Contributions

Conceptualization, I.K., S.K. and X.Z.; methodology, I.K., F.Y., X.Y. and X.Z.; software, S.K., F.Y., X.Z. and X.Y.; validation, X.Z., X.Y. and F.Y.; statistical data analysis, S.K. and F.Y.; investigation, F.Y., S.K., X.Y., X.Z. and I.K.; resources, X.Z. and I.K.; data curation, F.Y., X.Z. and X.Y.; writing—original draft preparation, review and editing, F.Y., S.K., X.Y., X.Z. and I.K.; visualization, F.Y. and S.K.; supervision, I.K., S.K. and X.Z.; project administration, X.Z., S.K. and I.K.; funding acquisition, X.Z., I.K. and S.K. All authors have read and agreed to the published version of the manuscript.

Funding

The project was supported by grants R01-CA93817 from National Cancer Institute, Grant Boost Program from Wayne State University, 1P20GM113226, and S10OD020106 from NIH. The Biostatistics Core in the Karmanos Cancer Institute at Wayne State University is supported in part by the NIH Cancer Center Support Grant P30 CA022453.

Acknowledgments

We thank Barbara Rusin, Maria Samerson, and Ann Bankowski from Wayne State University for technical support and/or collection/selection of blood samples from cases and controls.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cui, L.; Lu, H.T.; Lee, Y.H. Challenges and emergent solutions for LC-MS/MS based untargeted metabolomics in diseases. Mass Spectrom. Rev. 2018, 37, 772–792. [Google Scholar] [CrossRef] [PubMed]
  2. Jones, D.P.; Park, Y.; Ziegler, T.R. Nutritional metabolomics: Progress in addressing complexity in diet and health. Annu. Rev. Nutr. 2012, 32, 183–202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Mishur, R.J.; Rea, S.L. Applications of mass spectrometry to metabolomics and metabonomics: Detection of biomarkers of aging and of age-related diseases. Mass Spectrom. Rev. 2012, 31, 70–95. [Google Scholar] [CrossRef] [PubMed]
  4. Nicholson, J.K.; Lindon, J.C. Systems biology—Metabonomics. Nature 2008, 455, 1054–1056. [Google Scholar] [CrossRef]
  5. Yu, B.; Zanetti, K.A.; Temprosa, M.; Albanes, D.; Appel, N.; Barrera, C.B.; Ben-Shlomo, Y.; Boerwinkle, E.; Casas, J.P.; Clish, C.; et al. The consortium of metabolomics studies (COMETS): Metabolomics in 47 prospective cohort studies. Am. J. Epidemiol. 2019, 188, 991–1012. [Google Scholar] [CrossRef] [Green Version]
  6. Suhre, K.; Meisinger, C.; Doring, A.; Altmaier, E.; Belcredi, P.; Gieger, C.; Chang, D.; Milburn, M.V.; Gall, W.E.; Weinberger, K.M.; et al. Metabolic footprint of diabetes: A multiplatform metabolomics study in an epidemiological setting. PLoS ONE 2010, 5, e13953. [Google Scholar] [CrossRef] [Green Version]
  7. Prodhan, M.A.I.; Shi, B.Y.; Song, M.; He, L.Q.; Yuan, F.; Yin, X.M.; Bohman, P.; McClain, C.J.; Zhang, X. Integrating comprehensive two-dimensional gas chromatography mass spectrometry and parallel two-dimensional liquid chromatography mass spectrometry for untargeted metabolomics. Analyst 2019, 144, 4331–4341. [Google Scholar] [CrossRef]
  8. Liesenfeld, D.B.; Habermann, N.; Owen, R.W.; Scalbert, A.; Ulrich, C.M. Review of mass spectrometry–based metabolomics in cancer research. Cancer Epidemiol. Biomark. Prev. 2013, 22, 2182–2201. [Google Scholar] [CrossRef] [Green Version]
  9. Navarro-Reig, M.; Jaumot, J.; Baglai, A.; Vivó-Truyols, G.; Schoenmakers, P.J.; Tauler, R. Untargeted comprehensive two-dimensional liquid chromatography coupled with high-resolution mass spectrometry analysis of rice metabolome using multivariate curve resolution. Anal. Chem. 2017, 89, 7675–7683. [Google Scholar] [CrossRef]
  10. Shi, X.; Wahlang, B.; Wei, X.L.; Yin, X.M.; Falkner, K.C.; Prough, R.A.; Kim, S.H.; Mueller, E.G.; McClain, C.J.; Cave, M.; et al. Metabolomic analysis of the effects of polychlorinated biphenyls in nonalcoholic fatty liver disease. J. Proteome Res. 2012, 11, 3805–3815. [Google Scholar] [CrossRef]
  11. He, L.Q.; Li, F.Y.; Yin, X.M.; Bohman, P.; Kim, S.; McClain, C.J.; Feng, W.K.; Zhang, X. Profiling of Polar Metabolites in mouse feces using four analytical platforms to study the effects Of cathelicidin-related antimicrobial peptide in alcoholic liver disease. J. Proteome Res. 2019, 18, 2875–2884. [Google Scholar] [CrossRef] [PubMed]
  12. Yuan, F.; Harder, J.; Ma, J.; Yin, X.; Zhang, X.; Kosiewicz, M.M. Using multiple analytical platforms to investigate the androgen depletion effects on fecal metabolites in a mouse model of systemic lupus erythematosus. J. Proteome Res. 2020, 19, 667–676. [Google Scholar] [CrossRef] [PubMed]
  13. Pirok, B.W.J.; Stoll, D.R.; Schoenmakers, P.J. Recent developments in two-dimensional liquid chromatography: Fundamental improvements for practical applications. Anal. Chem. 2019, 91, 240–263. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Hanahan, D.; Weinberg, R.A. Hallmarks of cancer: The next generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [Green Version]
  15. Bellastella, G.; Scappaticcio, L.; Esposito, K.; Giugliano, D.; Maiorino, M.I. Metabolic syndrome and cancer: “The common soil hypothesis”. Diabetes Res. Clin. Pract. 2018, 143, 389–397. [Google Scholar] [CrossRef]
  16. Micucci, C.; Valli, D.; Matacchione, G.; Catalano, A. Current perspectives between metabolic syndrome and cancer. Oncotarget 2016, 7, 38959–38972. [Google Scholar] [CrossRef] [Green Version]
  17. Luo, Z.J.; Saha, A.K.; Xiang, X.Q.; Ruderman, N.B. AMPK, the metabolic syndrome and cancer. Trends Pharmacol. Sci. 2005, 26, 69–76. [Google Scholar] [CrossRef]
  18. Howlader, N.; Noone, A.M.; Krapcho, M.; Miller, D.; Brest, A.; Yu, M.; Ruhl, J.; Tatalovich, Z.; Mariotto, A.; Lewis, D.R.; et al. SEER Cancer Statistics Review 1975–2016; National Cancer Institute: Bethesda, MD, USA, 2019. [Google Scholar]
  19. Nicholson, B.D.; Shinkins, B.; Pathiraja, I.; Roberts, N.W.; James, T.J.; Mallett, S.; Perera, R.; Primrose, J.N.; Mant, D. Blood CEA levels for detecting recurrent colorectal cancer. Cochrane Database Syst. Rev. 2015, 2015, CD011134. [Google Scholar] [CrossRef] [Green Version]
  20. Klavins, K.; Drexler, H.; Hann, S.; Koellensperger, G. Quantitative metabolite profiling utilizing parallel column analysis for simultaneous reversed-phase and hydrophilic interaction liquid chromatography separations combined with tandem mass spectrometry. Anal. Chem. 2014, 86, 4145–4150. [Google Scholar] [CrossRef]
  21. Kim, S.; Yin, X.; Prodhan, M.A.I.; Zhang, X.; Zhong, Z.; Kato, I. Global plasma profiling for colorectal cancer-associated volatile organic compounds: A proof-of-principle study. J. Chromatogr. Sci. 2019, 57, 385–396. [Google Scholar] [CrossRef]
  22. Chen, L.; Zhang, C.; Gui, Q.; Chen, Y.; Yang, Y. Ultraperformance liquid chromatography coupled with quadrupole timeofflight mass spectrometrybased metabolic profiling of human serum prior to and following radical resection of colorectal carcinoma. Mol. Med. Rep. 2015, 12, 6879–6886. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Hashim, N.A.A.; Ab-Rahim, S.; Suddin, L.S.; Saman, M.S.A.; Mazlan, M. Global serum metabolomics profiling of colorectal cancer. Mol. Clin. Oncol. 2019, 11, 3–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Long, Y.; Sanchez-Espiridion, B.; Lin, M.B.; White, L.; Mishra, L.; Raju, G.S.; Kopetz, S.; Eng, C.; Hildebrandt, M.A.T.; Chang, D.W.; et al. Global and targeted serum metabolic profiling of colorectal cancer progression. Cancer 2017, 123, 4066–4074. [Google Scholar] [CrossRef] [PubMed]
  25. Zamani, Z.; Arjmand, M.; Vahabi, F.; Eshaq Hosseini, S.M.; Fazeli, S.M.; Iravani, A.; Bayat, P.; Oghalayee, A.; Mehrabanfar, M.; Haj Hosseini, R.; et al. A metabolic study on colon cancer using 1h nuclear magnetic resonance spectroscopy. Biochem. Res. Int. 2014, 2014, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Shu, X.; Xiang, Y.B.; Rothman, N.; Yu, D.; Li, H.L.; Yang, G.; Cai, H.; Ma, X.; Lan, Q.; Gao, Y.T.; et al. Prospective study of blood metabolites associated with colorectal cancer risk. Int. J. Cancer 2018, 143, 527–534. [Google Scholar] [CrossRef] [Green Version]
  27. Wang, Z.Z.; Cui, B.B.; Zhang, F.; Yang, Y.; Shen, X.T.; Li, Z.; Zhao, W.W.; Zhang, Y.Y.; Deng, K.; Rong, Z.W.; et al. Development of a correlative strategy to discover colorectal tumor tissue derived metabolite biomarkers in plasma using untargeted metabolomics. Anal. Chem. 2019, 91, 2401–2408. [Google Scholar] [CrossRef]
  28. Zaimenko, I.; Jaeger, C.; Brenner, H.; Chang-Claude, J.; Hoffmeister, M.; Grotzinger, C.; Detjen, K.; Burock, S.; Schmitt, C.A.; Stein, U.; et al. Non-invasive metastasis prognosis from plasma metabolites in stage II colorectal cancer patients: The DACHS study. Int. J. Cancer 2019, 145, 221–231. [Google Scholar] [CrossRef]
  29. Geijsen, A.J.M.R.; Brezina, S.; Keski-Rahkonen, P.; Baierl, A.; Bachleitner-Hofmann, T.; Bergmann, M.M.; Boehm, J.; Brenner, H.; Chang-Claude, J.; van Duijnhoven, F.J.B.; et al. Plasma metabolites associated with colorectal cancer: A discovery-replication strategy. Int. J. Cancer 2019, 145, 1221–1231. [Google Scholar] [CrossRef] [Green Version]
  30. Deng, L.L.; Gu, H.W.; Zhu, J.J.; Gowda, G.A.N.; Djukovic, D.; Chiorean, E.G.; Raftery, D. Combining NMR and LC/MS using backward variable elimination: Metabolomics analysis of colorectal cancer, polyps and healthy controls. Anal. Chem. 2016, 88, 7975–7983. [Google Scholar] [CrossRef] [Green Version]
  31. Tan, B.B.; Qiu, Y.P.; Zou, X.; Chen, T.L.; Xie, G.X.; Cheng, Y.; Dong, T.T.; Zhao, L.J.; Feng, B.; Hu, X.F.; et al. Metabonomics Identifies Serum Metabolite Markers of Colorectal Cancer. J. Proteome Res. 2013, 12, 3000–3009. [Google Scholar] [CrossRef] [Green Version]
  32. Qiu, Y.P.; Cai, G.X.; Su, M.M.; Chen, T.L.; Zheng, X.J.; Xu, Y.; Ni, Y.; Zhao, A.H.; Xu, L.X.; Cai, S.J.; et al. Serum metabolite profiling of human colorectal cancer using GC-TOFMS and UPLC-QTOFMS. J. Proteome Res. 2009, 8, 4844–4850. [Google Scholar] [CrossRef] [PubMed]
  33. Nishiumi, S.; Kobayashi, T.; Ikeda, A.; Yoshie, T.; Kibi, M.; Izumi, Y.; Okuno, T.; Hayashi, N.; Kawano, S.; Takenawa, T.; et al. A novel serum metabolomics-based diagnostic approach for colorectal cancer. PLoS ONE 2012, 7, e40459. [Google Scholar] [CrossRef] [PubMed]
  34. Cross, A.J.; Moore, S.C.; Boca, S.; Huang, W.Y.; Xiong, X.Q.; Stolzenberg-Solomon, R.; Sinha, R.; Sampson, J.N. A prospective study of serum metabolites and colorectal cancer risk. Cancer 2014, 120, 3049–3057. [Google Scholar] [CrossRef] [PubMed]
  35. Uchiyama, K.; Yagi, N.; Mizushima, K.; Higashimura, Y.; Hirai, Y.; Okayama, T.; Yoshida, N.; Katada, K.; Kamada, K.; Handa, O.; et al. Serum metabolomics analysis for early detection of colorectal cancer. J. Gastroenterol. 2017, 52, 677–694. [Google Scholar] [CrossRef]
  36. Nishiumi, S.; Kobayashi, T.; Kawana, S.; Unno, Y.; Sakai, T.; Okamoto, K.; Yamada, Y.; Sudo, K.; Yamaji, T.; Saito, Y.; et al. Investigations in the possibility of early detection of colorectal cancer by gas chromatography/triple-quadrupole mass spectrometry. Oncotarget 2017, 8, 17115–17126. [Google Scholar] [CrossRef] [Green Version]
  37. Ma, Y.L.; Liu, W.J.; Peng, J.Y.; Huang, L.; Zhang, P.; Zhao, X.P.; Cheng, Y.Y.; Qin, H.L. A pilot study of gas chromatograph/mass spectrometry-based serum metabolic profiling of colorectal cancer after operation. Mol. Biol. Rep. 2010, 37, 1403–1411. [Google Scholar] [CrossRef]
  38. Ritchie, S.A.; Ahiahonu, P.W.K.; Jayasinghe, D.; Heath, D.; Liu, J.; Lu, Y.S.; Jin, W.; Kavianpour, A.; Yamazaki, Y.; Khan, A.M.; et al. Reduced levels of hydroxylated, polyunsaturated ultra long-chain fatty acids in the serum of colorectal cancer patients: Implications for early screening and detection. BMC Med. 2010, 8, 13. [Google Scholar] [CrossRef] [Green Version]
  39. Phua, L.C.; Chue, X.P.; Koh, P.K.; Cheah, P.Y.; Ho, H.K.; Chan, E.C.Y. Non-invasive fecal metabonomic detection of colorectal cancer. Cancer Biol. Ther. 2014, 15, 389–397. [Google Scholar] [CrossRef] [Green Version]
  40. Lin, Y.; Ma, C.C.; Bezabeh, T.; Wang, Z.N.; Liang, J.H.; Huang, Y.; Zhao, J.Y.; Liu, X.M.; Ye, W.; Tang, W.; et al. 1H NMR-based metabolomics reveal overlapping discriminatory metabolites and metabolic pathway disturbances between colorectal tumor tissues and fecal samples. Int. J. Cancer 2019, 145, 1679–1689. [Google Scholar] [CrossRef]
  41. Denkert, C.; Budczies, J.; Weichert, W.; Wohlgemuth, G.; Scholz, M.; Kind, T.; Niesporek, S.; Noske, A.; Buckendahl, A.; Dietel, M.; et al. Metabolite profiling of human colon carcinoma-deregulation of TCA cycle and amino acid turnover. Mol. Cancer 2008, 7, 72. [Google Scholar] [CrossRef] [Green Version]
  42. Chan, E.C.Y.; Koh, P.K.; Mal, M.; Cheah, P.Y.; Eu, K.W.; Backshall, A.; Cavill, R.; Nicholson, J.K.; Keun, H.C. Metabolic profiling of human colorectal cancer using high-resolution magic angle spinning nuclear magnetic resonance (HR-MAS NMR) spectroscopy and gas chromatography mass spectrometry (GC/MS). J. Proteome Res. 2009, 8, 352–361. [Google Scholar] [CrossRef] [PubMed]
  43. Brown, D.G.; Rao, S.; Weir, T.L.; O’Malia, J.; Bazan, M.; Brown, R.J.; Ryan, E.P. Metabolomics and metabolic pathway networks from human colorectal cancers, adjacent mucosa and stool. Cancer Metab. 2016, 4, 11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Qiu, Y.P.; Cai, G.X.; Su, M.M.; Chen, T.L.; Liu, Y.M.; Xu, Y.; Ni, Y.; Zhao, A.H.; Cai, S.J.; Xu, L.X.; et al. Urinary metabonomic study on colorectal cancer. J. Proteome Res. 2010, 9, 1627–1634. [Google Scholar] [CrossRef] [PubMed]
  45. Cheng, Y.; Xie, G.X.; Chen, T.L.; Qiu, Y.P.; Zou, X.; Zheng, M.H.; Tan, B.B.; Feng, B.; Dong, T.T.; He, P.A.; et al. Distinct urinary metabolic profile of human colorectal cancer. J. Proteome Res. 2012, 11, 1354–1363. [Google Scholar] [CrossRef]
  46. Di Lena, M.; Porcelli, F.; Altomare, D.F. Volatile organic compounds as new biomarkers for colorectal cancer: A review. Colorectal Dis. 2016, 18, 654–663. [Google Scholar] [CrossRef]
  47. Kyle, J.E.; Casey, C.P.; Stratton, K.G.; Zink, E.M.; Kim, Y.M.; Zheng, X.Y.; Monroe, M.E.; Weitz, K.K.; Bloodsworth, K.J.; Orton, D.J.; et al. Comparing identified and statistically significant lipids and polar metabolites in 15-year old serum and dried blood spot samples for longitudinal studies. Rapid Commun. Mass Spectrom. 2017, 31, 447–456. [Google Scholar] [CrossRef] [Green Version]
  48. Anne, L. (Ed.) Glutamine metabolism in cancer. In The Heterogeneity of Cancer Metabolism; Springer International Publishing: Cham, Swizerland, 2018; pp. 13–32. [Google Scholar]
  49. Amelio, I.; Cutruzzola, F.; Antonov, A.; Agostini, M.; Melino, G. Serine and glycine metabolism in cancer. Trends Biochem. Sci. 2014, 39, 191–198. [Google Scholar] [CrossRef]
  50. Budczies, J.; Brockmoller, S.F.; Muller, B.M.; Barupal, D.K.; Richter-Ehrenstein, C.; Kleine-Tebbe, A.; Griffin, J.L.; Oresic, M.; Dietel, M.; Denkert, C.; et al. Comparative metabolomics of estrogen receptor positive and estrogen receptor negative breast cancer: Alterations in glutamine and beta-alanine metabolism. J. Proteom. 2013, 94, 279–288. [Google Scholar] [CrossRef]
  51. Hutschenreuther, A.; Birkenmeier, G.; Bigl, M.; Krohn, K.; Birkemeyer, C. Glycerophosphoglycerol, beta-alanine and pantothenic acid as metabolic companions of glycolytic activity and cell migration in breast cancer cell lines. Metabolites 2013, 3, 1084–1101. [Google Scholar] [CrossRef] [Green Version]
  52. Dai, X.F.; Luo, Y.N.; Xu, Y.; Zhang, J.Y. Key indexes and the emerging tool for tumor microenvironment editing. Am. J. Cancer Res. 2019, 9, 1027–1042. [Google Scholar]
  53. Bhowmik, S.K.; Ramirez-Pena, E.; Arnold, J.M.; Putluri, V.; Sphyris, N.; Michailidis, G.; Putluri, N.; Ambs, S.; Sreekumar, A.; Mani, S.A. EMT-induced metabolite signature identifies poor clinical outcome. Oncotarget 2015, 6, 42651–42660. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Bellissimo, M.P.; Cai, Q.P.; Ziegler, T.R.; Liu, K.H.; Tran, P.H.; Voss, M.B.; Martin, G.S.; Jones, D.P.; Yu, T.W.; Alvarez, J.A. Plasma high-resolution metabolomics differentiates adults with normal weight obesity from lean individuals. Obesity 2019, 27, 1729–1737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Tsai, C.K.; Yeh, T.S.; Wu, R.C.; Lai, Y.C.; Chiang, M.H.; Lu, K.Y.; Hung, C.Y.; Ho, H.Y.; Cheng, M.L.; Lin, G.G. Metabolomic alterations and chromosomal instability status in gastric cancer. World J. Gastroenterol. 2018, 24, 3760–3769. [Google Scholar] [CrossRef] [PubMed]
  56. Cano, K.E.; Li, L.; Bhatia, S.; Bhatia, R.; Forman, S.J.; Chen, Y. NMR-based metabolomic analysis of the molecular pathogenesis of therapy-related myelodysplasia/acute myeloid leukemia. J. Proteome Res. 2011, 10, 2873–2881. [Google Scholar] [CrossRef] [Green Version]
  57. Li, F.; Li, X.; Miao, Y.; Shan, C.; Yuan, F.; Ma, C.; Wang, Q.; Chen, J.; Chen, Y. Uhplc–ms-based metabolomics analysis on mice bearing neoplasm (h22) for hispidulin. J. Pharm. Biomed. Anal. 2016, 125, 310–318. [Google Scholar] [CrossRef]
  58. Chen, W.Y.; Wu, F.; You, Z.Y.; Zhang, Z.M.; Guo, Y.L.; Zhong, L.X. Analyzing the differentially expressed genes and pathway cross-talk in aggressive breast cancer. J. Obstet. Gynaecol. Res. 2015, 41, 132–140. [Google Scholar] [CrossRef]
  59. Kim, S.; You, S.; Hwang, D. Aminoacyl-tRNA synthetases and tumorigenesis: More than housekeeping. Nat. Rev. Cancer 2011, 11, 708–718. [Google Scholar] [CrossRef]
  60. Kim, D.; Kwon, N.H.; Kim, S. Association of Aminoacyl-tRNA Synthetases with Cancer. In Aminoacyl-Trna Synthetases in Biology and Medicine; Springer: Dordrecht, The Netherlands, 2013; pp. 207–245. [Google Scholar] [CrossRef]
  61. Kushner, J.P.; Boll, D.; Quagliana, J.; Dickman, S. Elevated methionine-trna synthetase activity in human colon cancer. Proc. Soc. Exp. Biol. Med. 1976, 153, 273–276. [Google Scholar] [CrossRef]
  62. Pavon-Eternod, M.; Gomes, S.; Rosner, M.R.; Pan, T. Overexpression of initiator methionine tRNA leads to global reprogramming of tRNA expression and increased proliferation in human epithelial cells. RNA 2013, 19, 461–466. [Google Scholar] [CrossRef] [Green Version]
  63. Yao, P.; Fox, P.L. Aminoacyl-tRNA synthetases in medicine and disease. EMBO Mol. Med. 2013, 5, 332–343. [Google Scholar] [CrossRef]
  64. Wei, X.L.; Shi, X.; Koo, I.; Kim, S.; Schmidt, R.H.; Arteel, G.E.; Watson, W.H.; McClain, C.; Zhang, X. MetPP: A computational platform for comprehensive two-dimensional gas chromatography time-of-flight mass spectrometry-based metabolomics. Bioinformatics 2013, 29, 1786–1792. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Lustgarten, M.S.; Price, L.L.; Phillips, E.M.; Fielding, R.A. Serum Glycine Is Associated with Regional Body Fat and Insulin Resistance in Functionally-Limited Older Adults. PLoS ONE 2013, 8, e84034. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Palsson-McDermott, E.M.; O’Neill, L.A.J. The Warburg effect then and now: From cancer to inflammatory diseases. Bioessays 2013, 35, 965–973. [Google Scholar] [CrossRef] [PubMed]
  67. Salminen, A.; Kauppinen, A.; Hiltunen, M.; Kaarniranta, K. Krebs cycle intermediates regulate DNA and histone methylation: Epigenetic impact on the aging process. Ageing Res. Rev. 2014, 16, 45–65. [Google Scholar] [CrossRef]
  68. Altomare, D.F.; Di Lena, M.; Porcelli, F.; Travaglio, E.; Longobardi, F.; Tutino, M.; Depalma, N.; Tedesco, G.; Sardaro, A.; Memeo, R.; et al. Effects of Curative Colorectal Cancer Surgery on Exhaled Volatile Organic Compounds and Potential Implications in Clinical Follow-up. Ann. Surg. 2015, 262, 862–867. [Google Scholar] [CrossRef]
  69. Kato, I.; Land, S.; Majumdar, A.P.; Barnholtz-Sloan, J.; Severson, R.K. Functional polymorphisms to modulate luminal lipid exposure and risk of colorectal cancer. Cancer Epidemiol. 2010, 34, 291–297. [Google Scholar] [CrossRef] [Green Version]
  70. Koo, I.; Shi, X.; Kim, S.; Zhang, X. iMatch2: Compound identification using retention index for analysis of gas chromatography-mass spectrometry data. J. Chromatogr. A 2014, 1337, 202–210. [Google Scholar] [CrossRef] [Green Version]
  71. Wei, X.L.; Shi, X.; Kim, S.; Patrick, J.S.; Binkley, J.; Kong, M.Y.; McClain, C.; Zhang, X. Data Dependent Peak Model Based Spectrum Deconvolution for Analysis of High Resolution LC-MS Data. Anal. Chem. 2014, 86, 2156–2165. [Google Scholar] [CrossRef]
  72. Wei, X.L.; Sun, W.L.; Shi, X.; Koo, I.; Wang, B.; Zhang, J.; Yin, X.M.; Tang, Y.N.; Bogdanov, B.; Kim, S.; et al. MetSign: A Computational Platform for High-Resolution Mass Spectrometry-Based Metabolomics. Anal. Chem. 2011, 83, 7668–7675. [Google Scholar] [CrossRef] [Green Version]
  73. Wei, X.L.; Shi, X.; Kim, S.; Zhang, L.; Patrick, J.S.; Binkey, J.; McClain, C.; Zhang, X. Data Preprocessing Method for Liquid Chromatography-Mass Spectrometry Based Metabolomics. Anal. Chem. 2012, 84, 7963–7971. [Google Scholar] [CrossRef] [Green Version]
  74. Chong, J.; Soufan, O.; Li, C.; Caraus, I.; Li, S.; Bourque, G.; Wishart, D.S.; Xia, J. MetaboAnalyst 4.0: Towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018, 46, W486–W494. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Overlap of metabolites identified from 3 analytical platforms.
Figure 1. Overlap of metabolites identified from 3 analytical platforms.
Metabolites 10 00343 g001
Figure 2. Partial least squares-discriminant analysis (PLS-DA) results of data from 3 platforms.
Figure 2. Partial least squares-discriminant analysis (PLS-DA) results of data from 3 platforms.
Metabolites 10 00343 g002
Figure 3. Instrument response of glutamic acid detected by 3 platforms. (A) Three-dimensional chromatographic peak of glutamic acid in GC × GC-MS. (B) Extracted ion chromatogram of glutamic acid in the same sample detected by two-dimensional gas and liquid chromatography-mass spectrometry (2DLC-MS) (+). (C) Extracted ion chromatogram of glutamic acid in the same sample detected by 2DLC-MS (−).
Figure 3. Instrument response of glutamic acid detected by 3 platforms. (A) Three-dimensional chromatographic peak of glutamic acid in GC × GC-MS. (B) Extracted ion chromatogram of glutamic acid in the same sample detected by two-dimensional gas and liquid chromatography-mass spectrometry (2DLC-MS) (+). (C) Extracted ion chromatogram of glutamic acid in the same sample detected by 2DLC-MS (−).
Metabolites 10 00343 g003
Figure 4. Pathways affected by the metabolic differences between cases and controls. The color and size of circles represent the sizes of the p-values (the darker the smaller) and the pathway impact scores (the larger the bigger), respectively.
Figure 4. Pathways affected by the metabolic differences between cases and controls. The color and size of circles represent the sizes of the p-values (the darker the smaller) and the pathway impact scores (the larger the bigger), respectively.
Metabolites 10 00343 g004
Figure 5. Receiver operating characteristic (ROC) results of (A) one standard error (1SE) tuned and (B) minimum (OPT) tuned.
Figure 5. Receiver operating characteristic (ROC) results of (A) one standard error (1SE) tuned and (B) minimum (OPT) tuned.
Metabolites 10 00343 g005
Table 1. Metabolites with significant abundance change between cases and controls.
Table 1. Metabolites with significant abundance change between cases and controls.
Namep-ValueFold Change (Cases/Controls)Platform%FO *
Controls
n = 29
Cases
n = 29
3-Hydroxybutyric acid9.98 × 10−32.122DLC-MS (−)96.5596.55
4-Methyl-2-oxovaleric acid4.62 × 10−30.712DLC-MS (−)100.00100.00
Adonitol1.07 × 10−21.272DLC-MS (−)62.0762.07
Alanine4.98 × 10−30.522DLC-MS (−)100.00100.00
Arginine3.83 × 10−20.782DLC-MS (+)82.7682.76
Aspartic Acid5.09 × 10−30.742DLC-MS (+)79.3179.31
beta-Alanine4.15 × 10−21.76GC × GC-MS100.00100.00
Choline8.56 × 10−31.192DLC-MS (+)68.9772.41
Citric acid2.13 × 10−21.512DLC-MS (−)100.00100.00
Ethanolamine3.17 × 10−21.312DLC-MS (+)100.00100.00
Glutamic acid2.58 × 10−32.00GC × GC-MS68.9772.41
Glutaric acid4.86 × 10−21.162DLC-MS (−)96.5596.55
Glycine1.44 × 10−31.342DLC-MS (+)93.1089.66
Hydracrylic acid4.32 × 10−21.28GC × GC-MS100.00100.00
Lactic acid3.11 × 10−21.13GC × GC-MS100.00100.00
Lysine1.86 × 10−21.30GC × GC-MS79.3182.76
Lysophosphatidic acid4.57 × 10−21.232DLC-MS (+)100.00100.00
Malate2.09 × 10−21.372DLC-MS (−)100.00100.00
Malate4.91 × 10−21.51GC × GC-MS72.4172.41
Neopentyl glycol3.24 × 10−21.16GC × GC-MS100.00100.00
Oleamide3.40 × 10−21.112DLC-MS (+)72.4172.41
Pinolenic Acid2.26 × 10−20.562DLC-MS (+)72.4175.86
Piperine4.29 × 10−20.632DLC-MS (+)100.00100.00
Propylene glycol4.03 × 10−30.50GC × GC-MS96.55100.00
Pyroglutamic acid1.21 × 10−21.552DLC-MS (−)72.4172.41
Testosterone sulfate2.83 × 10−20.542DLC-MS (−)100.00100.00
Tricine1.81 × 10−21.60GC × GC-MS86.2182.76
Uridine2.48 × 10−21.252DLC-MS (−)96.5596.55
* FO: Percentage frequency of occurrence.
Table 2. Pathways affected by the metabolic differences between cases and controls.
Table 2. Pathways affected by the metabolic differences between cases and controls.
NameTotalHitsp-ValueImpact
Aminoacyl-tRNA biosynthesis4864.35 × 10−50
Alanine, aspartate, and glutamate metabolism2845.97 × 10−40.42068
Arginine biosynthesis1439.35 × 10−40.19289
Glyoxylate and dicarboxylate metabolism3241.01× 10−30.13757
beta-Alanine metabolism2133.19 × 10−30.39925
Glutathione metabolism2837.34 × 10−30.11548
Glycerophospholipid metabolism3631.48 × 10−20.07130
Butanoate metabolism1521.97 × 10−20
Histidine metabolism1622.23 × 10−20
Pantothenate and CoA biosynthesis1923.09 × 10−20.02143
Citrate cycle (TCA cycle)2023.40 × 10−20.13450
Pyruvate metabolism2224.06 × 10−20.03110
Propanoate metabolism2324.41 × 10−20.00000

Share and Cite

MDPI and ACS Style

Yuan, F.; Kim, S.; Yin, X.; Zhang, X.; Kato, I. Integrating Two-Dimensional Gas and Liquid Chromatography-Mass Spectrometry for Untargeted Colorectal Cancer Metabolomics: A Proof-of-Principle Study. Metabolites 2020, 10, 343. https://doi.org/10.3390/metabo10090343

AMA Style

Yuan F, Kim S, Yin X, Zhang X, Kato I. Integrating Two-Dimensional Gas and Liquid Chromatography-Mass Spectrometry for Untargeted Colorectal Cancer Metabolomics: A Proof-of-Principle Study. Metabolites. 2020; 10(9):343. https://doi.org/10.3390/metabo10090343

Chicago/Turabian Style

Yuan, Fang, Seongho Kim, Xinmin Yin, Xiang Zhang, and Ikuko Kato. 2020. "Integrating Two-Dimensional Gas and Liquid Chromatography-Mass Spectrometry for Untargeted Colorectal Cancer Metabolomics: A Proof-of-Principle Study" Metabolites 10, no. 9: 343. https://doi.org/10.3390/metabo10090343

APA Style

Yuan, F., Kim, S., Yin, X., Zhang, X., & Kato, I. (2020). Integrating Two-Dimensional Gas and Liquid Chromatography-Mass Spectrometry for Untargeted Colorectal Cancer Metabolomics: A Proof-of-Principle Study. Metabolites, 10(9), 343. https://doi.org/10.3390/metabo10090343

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