Next Article in Journal
Early Postnatal Genistein Administration Affects Mice Metabolism and Reproduction in a Sexually Dimorphic Way
Previous Article in Journal
Hypoxia Performance Curve: Assess a Whole-Organism Metabolic Shift from a Maximum Aerobic Capacity towards a Glycolytic Capacity in Fish
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Azoxymethane Alters the Plasma Metabolome to a Greater Extent in Mice Fed a High-Fat Diet Compared to an AIN-93 Diet

1
Grand Forks Human Nutrition Research Center, Agricultural Research Service, United States Department of Agriculture, Grand Forks, ND 58203, USA
2
Department of Surgery and University of Kansas Cancer Center, Kansas City, KS 66160, USA
3
School of Public Health and Health Sciences, University of Massachusetts, Amherst, MA 01003, USA
*
Author to whom correspondence should be addressed.
Submission received: 7 May 2021 / Revised: 1 July 2021 / Accepted: 2 July 2021 / Published: 9 July 2021
(This article belongs to the Special Issue Colon Cancer Metabolism)

Abstract

:
Consumption of a high-fat diet (HFD) links obesity to colon cancer in humans. Our data show that a HFD (45% energy fat versus 16% energy fat in an AIN-93 diet (AIN)) promotes azoxymethane (AOM)-induced colonic aberrant crypt foci (ACF) formation in a mouse cancer model. However, the underlying metabolic basis remains to be determined. In the present study, we hypothesize that AOM treatment results in different plasma metabolomic responses in diet-induced obese mice. An untargeted metabolomic analysis was performed on the plasma samples by gas chromatography time-of-flight mass spectrometry (GC-TOF-MS). We found that 53 of 144 identified metabolites were different between the 4 groups of mice (AIN, AIN + AOM, HFD, HFD + AOM), and sparse partial least-squares discriminant analysis showed a separation between the HFD and HFD + AOM groups but not the AIN and AIN + AOM groups. Moreover, the concentrations of dihydrocholesterol and cholesterol were inversely associated with AOM-induced colonic ACF formation. Functional pathway analyses indicated that diets and AOM-induced colonic ACF modulated five metabolic pathways. Collectively, in addition to differential plasma metabolomic responses, AOM treatment decreases dihydrocholesterol and cholesterol levels and alters the composition of plasma metabolome to a greater extent in mice fed a HFD compared to the AIN.

1. Introduction

Consumption of a high-fat diet (HFD) links obesity to colon cancer in humans [1,2]. The global obesity epidemic has, in part, been attributed to the adoption of Western lifestyle practices, including increased consumption of high-energy diets [1,2,3]. Diet-induced obesity is now established as a risk factor for cancer, and overweight and obesity affect two-thirds of Americans and an estimated 2.3 billion people worldwide [4]. In mice, consumption of a HFD can lead to the accumulation of excess body fat that is associated with adipose tissue dysfunction and a chronic state of low-grade inflammation known to promote tumor development [5,6]. However, there are few mechanistic studies on early colonic tumorigenesis concerning the underlying metabolic regulation in the context of HFD-induced obesity.
We recently reported that in a mouse model, a HFD promotes colonic aberrant crypt foci (ACF, putative preneoplastic lesions) formation accompanied by increased systemic levels of proinflammatory cytokines [7]. We therefore posit that determining plasma metabolomic responses may lead to a greater understanding of the metabolic regulation concerning a HFD and colonic ACF formation. As small-molecule metabolites exhibit an organisms’ physiological condition at a given moment, these metabolites represent a dynamic situation in response to biochemical and pathological changes [8,9]. Thus, untargeted metabolite profiling (e.g., metabolomics) is an effective approach that aids in determining these underlying biochemical changes [10]. It is known that a HFD promotes AOM-induced ACF formation in a mouse model, and differential metabolites were identified in the serum of colon cancer patients compared to the healthy controls [7,11,12]. Therefore, we hypothesize that AOM treatment may cause different plasma metabolomic responses in AOM mouse colon cancer models fed with the HFD compared to the AIN.

2. Results

2.1. Identified Metabolites and Their Group Separation

The plasma samples were collected at the end of the study (Week 14) [7]. We identified 144 metabolites (Table S1) from 555 discrete signals detected in the plasma by using GC-TOFMS. There were 53 of 144 metabolites significantly different between the 4 groups of mice (AIN, AIN + AOM; HFD, HFD + AOM), and the relative values for these 53 metabolites compared to the AIN group are shown in Table 1.
To examine the diet-AOM effect on metabolome, sparse partial least-squares discriminant analysis (sPLS-DA) was used to visualize the metabolite group separation and individual sample variability (Figure 1A). We found that there was a separation between the HFD and HFD + AOM groups but not the AIN and AIN + AOM groups.
The 10 major metabolites influencing separation along Component 1 were (1) beta-sitosterol (>120% increase), dihydrocholesterol (>89% increase), alpha-tocopherol (>94% increase), methanol phosphate (>18% increase), cholesterol (>32% increase), arachidonic acid (>6% increase), and nicotinamide (>58% increase), which were the highest concentrations in the HFD group compared to the other three groups; (2) the concentrations of citric acid, isocitric acid and myristic acid were decreased by (>23%) in the HFD group compared to the other three groups (Figure 1B).
The 10 major metabolites influencing separation along Component 2 were: (1) oxoproline (>29% increase), oleic acid (>72% increase), 3-ureidopropionate (>67% increase), 2-aminobutyric acid (>58% increase), ornithine (>65% increase), and uracil (>25% increase) were the highest concentration in the HFD + AOM group compared to the otherthree3 groups; (2) palmitoleic acid (>30% decrease), alpha-ketoglutarate (>28% decrease), malic acid (>40% decrease), and citrulline (>30% decrease) exhibited a lower concentration in the HFD and HFD + AOM groups compared with the AIN or AIN + AOM group (Figure 1C).

2.2. Diet and AOM Interaction

To visualize the effect of diet and AOM interaction on the abundance of these 53 metabolites, a two-way ANOVA (including interaction tests) was performed. In Figure 2, there were 25 metabolites (green circle) that exhibited overall diet x AOM interactions in which (1) the concentration of 12 metabolites (the overlapped green and red circles) differed due to different diets while the concentration of 18 metabolites (the overlapped green and blue circles) differed because of AOM treatment; (2) the concentration of 10 metabolites (the overlapped green, red and blue circles) differed due to both diet and AOM treatment.
These 25 interactive metabolites included (1) the concentrations of alpha-tocopherol, dihydrocholesterol and nicotinamide were increased (> 58%) in the HFD group compared to the other 3 groups; in contrast, the concentrations of citric acid and myristic acid were decreased at least 23% in the HFD group compared to the other 3 groups (Table 1). These 5 metabolites comprise 50% of component 1 (Figure 1B); (2) the concentrations of oxoproline, 2-aminobutyric acid, uracil and ornithine were increased (> 29%) in the HFD + AOM group compared to the other 3 groups while palmitoleic acid was decreased (> 80%) in the (HFD or HFD + AOM) group compared to the (LFD or LFD + AOM) group (Table 1). These 5 metabolites consist 50% of component 2 (Figure 1C).

2.3. Metabolic Pathways of Altered Metabolites

To examine the biological significance of these altered 53 metabolites, we determined the potential pathways using the metabolite set enrichment analysis. These pathways (Figure 3, Table 2) included (1) the citric acid cycle: oxoglutaric acid, succinic acid, isocitric acid, cis-aconitic acid, citric acid, pyruvic acid, fumaric acid; (2) arginine biosynthesis: citrulline, aspartic acid, ornithine, oxoglutaric acid, fumaric acid; (3) aminoacyl-tRNA biosynthesis: phenylalanine, glycine, aspartic acid, serine, methionine, lysine, threonine, tyrosine; (4) alanine, aspartate, and glutamate metabolism: aspartic acid, citric acid, fumaric acid, pyruvic acid, succinic acid, oxoglutaric acid; and (5) glyoxylate and dicarboxylate metabolism: cis-aconitic acid, citric acid, serine, glycine, isocitric acid, pyruvic acid (Figure 3, Table 1 and Table 2).

2.4. Correlation between Colonic ACF and Dihydrocholesterol/Cholesterol

A scatterplot analysis was used to determine the correlation between the number of ACF and the concentration of dihydrocholesterol/cholesterol. The AOM-induced colonic ACF data were taken from our previous study [7], and dihydrocholesterol/cholesterol data were from the current study (Table 1 and Table S1). The analysis (Figure 4) showed that there was an inverse association between the ACF number and dihydrocholesterol/cholesterol concentrations.

3. Discussion

Epidemiological and experimental data suggest that obesity increases colon cancer risk [2,3,6,7,8,9]. Our previous data demonstrate that a HFD promotes colonic ACF formation with increased levels of circulating proinflammatory cytokines in a mouse model [7]. In the present study, we found that 53 of 144 identified metabolites were different between the 4 experimental groups (Table 1, Figure S1). These data (1) confirm the hypothesis that AOM treatment causes different plasma metabolomic responses in AOM mouse colon cancer model fed with the HFD compared to the AIN; (2) support the notion that metabolomic fingerprints in the plasma may be used as potential markers for identifying changes in the connection between diet and colon cancer [13].
The separation between the HFD and HFD + AOM treatment groups but not the AIN and AIN + AOM groups (Figure 1A) indicates that AOM-induced colonic ACF formation altered the plasma metabolome to a greater extent in mice fed a HFD compared to the AIN. Compared to Component 2 of the sPLS-DA, Component 1 played a greater role (Figure 1B,C) in separating treatment groups [14,15]. The concentrations of alpha-tocopherol, nicotinamide, and beta-sitosterol were higher in the HFD group compared to the AIN group (Table 1, Figure 1B), which may be due to the increased dietary intakes (e.g., the alpha-tocopherol and beta-sitosterol in corn-oil) and homeostatic adaptation during the HFD consumption (e.g., nicotinamide) [7,16]. This observation is initially counterintuitive given the fact that alpha-tocopherol, beta-sitosterol, and nicotinamide are potential anticancer compounds [17,18,19]. For example, alpha-tocopherol exhibits superior antioxidant and anti-inflammatory properties via downregulation of transcription factor NF-kB activation [17]; beta-sitosterol and nicotinamide modulate multiple cell signaling pathways, including cell proliferation, angiogenesis, metastasis, inflammation, genomic stability, and immune-response [19,20,21]. However, compared to the AIN, the proinflammatory impact of HFD-induced obesity may outweigh the beneficial effects of these compounds (e.g., alpha-tocopherol); subsequently, AOM treatment induces a higher number of colonic ACF in the HFD group compared to the AIN group [7].
As colonic ACF formation is critically dependent on AOM bioactivation by the gut bacteria without affecting the daily food intake [7,22,23], the higher colonic ACF number may reduce the uptake of alpha-tocopherol, beta-sitosterol, and nicotinamide to a greater extent in the HFD + AOM group compared to the AIN + AOM group [7,24]. Thus, a low uptake efficacy of these beneficial compounds in the HFD + AOM group might further promote the vicious cycle of AOM-induced colonic ACF formation.
There are epidemiological reports on an inverse relationship between plasma cholesterol levels and certain colon cancer risk populations (e.g., during the 10 years preceding the cancer) [25,26,27]. In line with these epidemiological data [25,26,27], the concentrations of dihydrocholesterol and cholesterol in the AIN + AOM and HFD + AOM groups were decreased by 62% and 21% and 75% and 39% when compared to the AIN and HFD groups, respectively (Table 1, Figure 1B). The connection between declining serum cholesterol concentrations and colon cancer [25,26,27] may involve several factors. First, epithelial cell mutations occur during colonic tumorigenesis cause a decrease in intestinal cholesterol uptake [24,28]. Second, to meet the increasing demand for cell proliferation, cancer cells significantly enhance cholesterol absorption and synthesis [29,30]. Third, an intestinal cachexia may be a secondary result of metabolic and nutritional change in advanced tumorigenesis [27].
Dihydrocholesterol biosynthesis occurs through two major pathways: (1) unabsorbed cholesterol in the colon is biohydrogenated to dihydrocholesterol via bacterial enzymes [31,32], and (2) dihydrocholesterol is synthesized via a pathway with 7 alpha-hydroxylated C27-steroids as substrates in the liver [33]. Our data demonstrate that the levels of both plasma dihydrocholesterol (62 to 75% decrease) and cholesterol (21 to 39% decrease) are sensitive markers, which is inversely correlated to AOM-induced colonic aberrant crypt formation in the AIN and HFD groups (Figure 4) [7]. Because cholesterol is a precursor of dihydrocholesterol in the biosynthetic pathway [31,32,33], and the diet and AOM exhibit interactive effects on dihydrocholesterol but not cholesterol (Figure 2), these molecular events may account for a greater decrease of dihydrocholesterol concentrations compared to that of total plasma cholesterol (Table 1). As dihydrocholesterol is widely used as a surrogate marker of cholesterol synthesis and absorption [34], our data suggest that plasma dihydrocholesterol may be a more sensitive marker than cholesterol (62 to 75% vs. 21 to 39% decrease) for epidemiological studies [25,26,27] on colon cancer, although future human studies are needed to verify this finding.
Further functional analysis demonstrates that 53 plasma metabolites were involved in five major metabolic pathways (Figure 3, Table 2). First, the citrate cycle is a critical metabolic pathway that utilizes glucose, amino acids, and fatty acids. Emerging evidence demonstrates that cancer cells heavily rely on the citrate cycle for energy production and macromolecule synthesis [35]. Consequently, an altered citrate cycle pathway was detected in this study (Figure 3, Table 2), which is consistent with the serum metabolomic data from colon cancer patients [36,37]. Further studies are warranted to determine the connection between human colon cancer and plasma metabolomic profiles of the citrate cycle pathway in the context of obesity. Second, arginine biosynthesis and aminoacyl-tRNA biosynthesis are two pathways, which are involved in endothelial cell migration and angiogenesis in tumorigenesis [38,39]. However, there are scant data examining plasma metabolomic profiling in this regard. A change of these two pathways (Figure 3, Table 2) provides new insights into the plasma metabolome. Third, (a) alanine, aspartate, and glutamate metabolism and (b) glyoxylate and dicarboxylate metabolism are the other two pathways with no immediate functional connections with colon cancer. These two pathways are important for cellular homeostatic regulation because alanine, aspartate, and glutamate are derived from intermediates of central metabolism such as the citrate cycle pathway [40,41], a change in these two pathways (Figure 3, Table 2) may be the secondary results from the citrate cycle pathway.
In this report, although our data on differential plasma metabolomic responses are interesting, there are a few limitations. For example, our plasma metabolomic data show only a disease-associated relationship but not a causal relationship, and these data are from a single-time-point experiment. Because the cause of colon cancer is involved multiple metabolic pathways (collective effects), it is currently difficult for us to experimentally validate one or a few disease-causal metabolites. However, in the future, we will be able to shed some light on disease-causal metabolites if more metabolomic data/studies (with multiple experimental time points) are available.

4. Materials and Methods

4.1. Animals, Diets, and AOM Treatment

Three- to four-week-old male C57BL/6 mice (Harlan, Madison, WI, USA) were individually housed in PlexiglasTM ventilated cages within a pathogen-free facility that maintained a 12 h light/dark cycle and a temperature of 22 ± 1 °C. Mice were given free access to food and deionized water. This study was approved by the Animal Care and Use Committee of the Grand Forks Human Nutrition Research Center (Protocol Code: HZ13M2, 2013), and animals were maintained in accordance with NIH guidelines for the care and use of laboratory animals. The study design, diet composition, and preparation have been previously reported [7]. Briefly, C57BL/6 mice were fed either an AIN or HFD (n = 25/group). On Week 3, within a given diet group, mice received either weekly intraperitoneal injection of the colon carcinogen, AOM (Sigma, St. Louis, MO, USA) (n = 15/group), at a concentration of 8 mg/kg body weight [42] or phosphate-buffered saline (PBS, pH = 7.4) carrier solution (n = 10/group) for 4 weeks. At the termination of the experiment, mice were fasted for 6 h and then euthanized with a mixture of ketamine and xylazine (100 mg/kg body weight). Plasma samples were collected at the end of the study (Week 14) and stored at −80 °C for metabolomic analysis.

4.2. Plasma Metabolomics

Metabolomic analysis was performed at the West Coast Metabolomics Center (University of California, Davis Genomic Center, Davis, CA, USA) [14,15]. Plasma samples were extracted and derivatized by silylation and methyloximation and analyzed by GC-TOF-MS for untargeted metabolomics. Data were processed at the West Coast Metabolomics Center using the BinBase database [43]. Metabolite quantifier ion peak heights were normalized to the sum intensities of all known compounds and used for the follow-up statistical analyses (Table S1).

4.3. Statistical and Bioinformatic Analysis

Bioinformatic analysis: to avoid a skewed distribution, obtained (peak intensity) data were normalized by Log10 transformation with an autoscaling method which is highly recommended for most peak intensity data [44,45]. In addition, the suitability (our data distribution) of the above-normalized method was the best when compared to other available normalized methods (e.g., Pareto scaling) using the MetaboAnalyst software (version 5.0, McGill University, Sainte Anne de Bellevue, Quebec, Canada) [44,45]. The metabolite group separation and functional pathways were analyzed by sPLS-DA, heatmap clustering, and metabolite set enrichment analysis (MSEA) using the MetaboAnalyst software [44,45], respectively. The effects of diet (AIN or HFD) and AOM treatment (with or without) on the relative abundance of plasma metabolites were analyzed using a two-way analysis of variance (ANOVA) corrected by a false discovery rate (FDR) of 0.05 and Tukey’s contrasts for post hoc comparisons. Results are given as mean ± standard deviation (SD). JMP V14 (SAS Institute Inc., Cary, NC, USA) and MetaboAnalyst software (version 4.0, McGill University, Sainte Anne de Bellevue, QC, Canada) were used for all statistical analyses [44,45].

5. Conclusions

Our data demonstrated that AOM, a commonly used carcinogen in a chemical-induced colon cancer model, altered the plasma metabolome to a greater extent in mice fed the HFD compared to the AIN. These plasma metabolites were involved in five major metabolic pathways included (1) citrate cycle, (2) Arg biosynthesis, (3) aminoacyl-tRNA biosynthesis, (4) Ala, Asp, and Glu metabolism, and (5) glyoxylate and dicarboxylate metabolism. Our newly identified (inverse) association between plasma dihydrocholesterol levels and AOM-induced colonic ACF formation in a mouse model warrants future exploration as one of the unique plasma “multiple analytes” at the quest of human colon cancer biomarker search.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/metabo11070448/s1. Figure S1: Hierarchical clustering heatmap analysis of 53 differential metabolites between experimental groups. Each column represents the metabolite abundance of one animal plasma sample. Without AOM treatment (control), n = 10/group (AIN or HFD group). With AOM, treatment n = 15/group (AIN + AOM or HFD + AOM group). Table S1: Metabolites identified in plasma from C57BL/6 mice fed the AIN93G (AIN) or the high-fat diet (HFD) and with or without AOM treatment.

Author Contributions

Conceptualization, H.Z.; methodology, H.Z.; software, H.Z.; validation, H.Z.; formal analysis, H.Z., S.U., Z.L. and M.R.B.; data curation, H.Z.; writing—original draft preparation, H.Z., S.U., Z.L. and M.R.B.; writing—review and editing, H.Z., S.U., Z.L. and M.R.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the US Department of Agriculture, Agricultural Research Service, and Research Project 3062-51000-050-00D.

Institutional Review Board Statement

The animal experiments were approved by the Animal Care and Use Committee of the Grand Forks Human Nutrition Research Center (Protocol Code: HZ13M2, 2013), and animals were maintained in accordance with NIH guidelines for the care and use of laboratory animals.

Informed Consent Statement

Not applicable.

Data Availability Statement

All plasma metabolomics data used in this work are available as a Supplementary excel file (Table S1).

Acknowledgments

We greatly appreciate Kay Keehr and James Lindlauf for the technical support.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AC, aberrant crypt; ACF, aberrant crypt foci; AIN, a modified AIN93G formulation; Ala, alanine; ANOVA, analysis of variance; AOM, azoxymethane; Arg, arginine; Asp, aspartate; FDR, false discovery rate; Glutamate, Glu; sPLS-DA, sparse partial least-squares discriminant analysis; 2D, two-dimensional; GC-TOFMS; time-of-flight mass spectrometry; KEGG, Kyoto encyclopedia of genes and genomes; TCA cycle, citrate cycle.

References

  1. Zeng, H.; Lazarova, D.L. Obesity-related colon cancer: Dietary factors and their mechanisms of anticancer action. Clin. Exp. Pharmacol. Physiol. 2012, 39, 161–167. [Google Scholar] [CrossRef]
  2. Roslan, N.H.; Makpol, S.; Yusof, Y.A.M. A Review on Dietary Intervention in Obesity Associated Colon Cancer. Asian Pac. J. Cancer Prev. 2019, 20, 1309–1319. [Google Scholar] [CrossRef] [Green Version]
  3. O’Neill, A.M.; Burrington, C.M.; Gillaspie, E.A.; Lynch, D.T.; Horsman, M.; Greene, M.W. High-fat Western diet–induced obesity contributes to increased tumor growth in mouse models of human colon cancer. Nutr. Res. 2016, 36, 1325–1334. [Google Scholar] [CrossRef] [PubMed]
  4. Siegel, R.L.; Miller, K.D.; Fedewa, S.A.; Ahnen, D.J.; Meester, R.G.S.; Barzi, A.; Jemal, A. Colorectal cancer statistics, 2017. CA Cancer J. Clin. 2017, 67, 177–193. [Google Scholar] [CrossRef] [PubMed]
  5. Tarasiuk, A.; Mosińska, P.; Fichna, J. The mechanisms linking obesity to colon cancer: An overview. Obes. Res. Clin. Pract. 2018, 12, 251–259. [Google Scholar] [CrossRef]
  6. Penrose, H.M.; Heller, S.; Cable, C.; Nakhoul, H.; Baddoo, M.; Flemington, E.; Crawford, S.E.; Savkovic, S.D. High-fat diet induced leptin and Wnt expression: RNA-sequencing and pathway analysis of mouse colonic tissue and tumors. Carcinogenesis 2017, 38, 302–311. [Google Scholar] [CrossRef] [Green Version]
  7. Zeng, H.; Ishaq, S.; Liu, Z.; Bukowski, M. Colonic aberrant crypt formation accompanies an increase of opportunistic pathogenic bacteria in C57BL/6 mice fed a high-fat diet. J. Nutr. Biochem. 2018, 54, 18–27. [Google Scholar] [CrossRef] [Green Version]
  8. Chen, H.-H.; Tseng, Y.J.; Wang, S.-Y.; Tsai, Y.-S.; Chang, C.-S.; Kuo, T.-C.; Yao, W.-J.; Shieh, C.-C.; Wu, C.-H.; Kuo, P.-H. The metabolome profiling and pathway analysis in metabolic healthy and abnormal obesity. Int. J. Obes. 2015, 39, 1241–1248. [Google Scholar] [CrossRef] [PubMed]
  9. Palviainen, M.; Laukkanen, K.; Tavukcuoglu, Z.; Velagapudi, V.; Kärkkäinen, O.; Hanhineva, K.; Auriola, S.; Ranki, A.; Siljander, P. Cancer Alters the Metabolic Fingerprint of Extracellular Vesicles. Cancers 2020, 12, 3292. [Google Scholar] [CrossRef]
  10. Zhang, A.; Sun, H.; Wang, X. Power of metabolomics in biomarker discovery and mining mechanisms of obesity. Obes. Rev. 2012, 14, 344–349. [Google Scholar] [CrossRef]
  11. Ried, J.S.; Shin, S.-Y.; Krumsiek, J.; Illig, T.; Theis, F.J.; Spector, T.D.; Adamski, J.; Wichmann, H.-E.; Strauch, K.; Soranzo, N.; et al. Novel genetic associations with serum level metabolites identified by phenotype set enrichment analyses. Hum. Mol. Genet. 2014, 23, 5847–5857. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Guijas, C.; Montenegro-Burke, J.R.; Warth, B.; Spilker, M.; Siuzdak, G. Metabolomics activity screening for identifying metabolites that modulate phenotype. Nat. Biotechnol. 2018, 36, 316–320. [Google Scholar] [CrossRef] [PubMed]
  13. Zarei, I.; Oppel, R.C.; Borresen, E.C.; Brown, R.J.; Ryan, E.P. Modulation of plasma and urine metabolome in colorectal cancer survivors consuming rice bran. Integr. Food Nutr. Metab. 2019, 6, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Fiehn, O.; Garvey, W.T.; Newman, J.; Lok, K.H.; Hoppel, C.L.; Adams, S. Plasma Metabolomic Profiles Reflective of Glucose Homeostasis in Non-Diabetic and Type 2 Diabetic Obese African-American Women. PLoS ONE 2010, 5, e15234. [Google Scholar] [CrossRef] [Green Version]
  15. Piccolo, B.D.; Keim, N.L.; Fiehn, O.; Adams, S.H.; Van Loan, M.D.; Newman, J.W. Habitual physical activity and plasma metabolomic patterns distinguish individuals with low vs. high weight loss during controlled energy restriction. J. Nutr. 2015, 145, 681–690. [Google Scholar] [CrossRef] [Green Version]
  16. Hong, S.; Moreno-Navarrete, J.M.; Wei, X.; Kikukawa, Y.; Tzameli, I.; Prasad, D.; Lee, Y.; Asara, J.M.; Fernández-Real, J.M.; Maratos-Flier, E.; et al. Nicotinamide N-methyltransferase regulates hepatic nutrient metabolism through Sirt1 protein stabilization. Nat. Med. 2015, 21, 887–894. [Google Scholar] [CrossRef] [Green Version]
  17. Peh, H.Y.; Tan, W.S.D.; Liao, W.; Wong, W.F. Vitamin E therapy beyond cancer: Tocopherol versus tocotrienol. Pharmacol. Ther. 2016, 162, 152–169. [Google Scholar] [CrossRef]
  18. Nikas, I.P.; Paschou, S.A.; Ryu, H.S. The Role of Nicotinamide in Cancer Chemoprevention and Therapy. Biomolecules 2020, 10, 477. [Google Scholar] [CrossRef] [Green Version]
  19. Bin Sayeed, M.S.; Ameen, S.S. Beta-Sitosterol: A Promising but Orphan Nutraceutical to Fight against Cancer. Nutr. Cancer 2015, 67, 1214–1220. [Google Scholar] [CrossRef] [PubMed]
  20. Buqué, A.; Bloy, N.; Kroemer, G.; Galluzzi, L. Possible mechanisms of cancer prevention by nicotinamide. Br. J. Pharmacol. 2021, 178, 2034–2040. [Google Scholar] [CrossRef] [PubMed]
  21. Fania, L.; Mazzanti, C.; Campione, E.; Candi, E.; Abeni, D.; Dellambra, E. Role of Nicotinamide in Genomic Stability and Skin Cancer Chemoprevention. Int. J. Mol. Sci. 2019, 20, 5946. [Google Scholar] [CrossRef] [Green Version]
  22. Megaraj, V.; Ding, X.; Fang, C.; Kovalchuk, N.; Zhu, Y.; Zhang, Q.-Y. Role of Hepatic and Intestinal P450 Enzymes in the Metabolic Activation of the Colon Carcinogen Azoxymethane in Mice. Chem. Res. Toxicol. 2014, 27, 656–662. [Google Scholar] [CrossRef] [PubMed]
  23. Takada, H.; Hirooka, T.; Hiramatsu, Y.; Yamamoto, M. Effect of beta-glucuronidase inhibitor on azoxymethane-induced colonic carcinogenesis in rats. Cancer Res. 1982, 42, 331–334. [Google Scholar]
  24. Yu, B.; Peng, X.-H.; Wang, L.-Y.; Wang, A.-B.; Su, Y.-Y.; Chen, J.-H.; Zhang, X.-W.; Zhao, D.-Z.; Wang, H.; Pang, D.-X.; et al. Abnormality of intestinal cholesterol absorption in ApcMin/+ mice with colon cancer cachexia. Int. J. Clin. Exp. Pathol. 2019, 12, 759–767. [Google Scholar]
  25. Broitman, S.; Cerda, S.; Wilkinson, J. Cholesterol metabolism and colon cancer. Prog. Food Nutr. Sci. 1993, 17, 1–40. [Google Scholar]
  26. Winawer, S.J.; Flehinger, B.J.; Buchalter, J.; Herbert, E.; Shike, M. Declining serum cholesterol levels prior to diagnosis of colon cancer. A time-trend, case-control study. JAMA 1990, 263, 2083–2085. [Google Scholar] [CrossRef]
  27. Chung, Y.W.; Han, D.S.; Park, Y.K.; Son, B.K.; Paik, C.H.; Lee, H.L.; Jeon, Y.C.; Sohn, J.H. Association of obesity, serum glucose and lipids with the risk of advanced colorectal adenoma and cancer: A case-control study in Korea. Dig. Liver Dis. 2006, 38, 668–672. [Google Scholar] [CrossRef]
  28. Hamada, K.; Umemoto, A.; Kajikawa, A.; Tanaka, M.; Seraj, M.J.; Nakayama, M.; Kubota, A.; Monden, Y. Mucosa-specific DNA adducts in human small intestine: A comparison with the colon. Carcinogenesis 1994, 15, 2677–2680. [Google Scholar] [CrossRef]
  29. Zhang, X.; Zhao, X.W.; Liu, D.B.; Han, C.Z.; Du, L.L.; Jing, J.X.; Wang, Y. Lipid levels in serum and cancerous tissues of colorectal cancer patients. World J. Gastroenterol. 2014, 20, 8646–8652. [Google Scholar] [CrossRef] [PubMed]
  30. Miyamoto, S.; Narita, T.; Komiya, M.; Fujii, G.; Hamoya, T.; Nakanishi, R.; Tamura, S.; Kurokawa, Y.; Takahashi, M.; Mutoh, M. Novel screening system revealed that intracellular cholesterol trafficking can be a good target for colon cancer prevention. Sci. Rep. 2019, 9, 6192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Myant, N. The Biology of Cholesterol and Related Steroids; Heinemann Medical Books: London, UK, 1981; p. 910. [Google Scholar]
  32. Wang, X.; Guan, L.; Zhao, Y.; Lei, L.; Liu, Y.; Ma, K.Y.; Wang, L.; Man, S.W.; Wang, J.; Huang, Y.; et al. Plasma cholesterol-lowering activity of dietary dihydrocholesterol in hypercholesterolemia hamsters. Atherosclerosis 2015, 242, 77–86. [Google Scholar] [CrossRef] [PubMed]
  33. Skrede, S.; Björkhem, I.; Buchmann, M.S.; Hopen, G.; Fausa, O. A novel pathway for biosynthesis of cholestanol with 7 alpha-hydroxylated C27-steroids as intermediates, and its importance for the accumulation of cholestanol in cerebrotendinous xanthomatosis. J. Clin. Investig. 1985, 75, 448–455. [Google Scholar] [CrossRef]
  34. Lütjohann, D.; Björkhem, I.; Friedrichs, S.; Kerksiek, A.; Lövgren-Sandblom, A.; Geilenkeuser, W.-J.; Ahrends, R.; Andrade, I.; Ansorena, D.; Astiasarán, I.; et al. First international descriptive and interventional survey for cholesterol and non-cholesterol sterol determination by gas- and liquid-chromatography-Urgent need for harmonisation of analytical methods. J. Steroid Biochem. Mol. Biol. 2019, 190, 115–125. [Google Scholar] [CrossRef]
  35. Anderson, N.M.; Mucka, P.; Kern, J.G.; Feng, H. The emerging role and targetability of the TCA cycle in cancer metabolism. Protein Cell 2018, 9, 216–237. [Google Scholar] [CrossRef]
  36. Cho, S.; Song, N.; Choi, J.-Y.; Shin, A. Effect of Citric Acid Cycle Genetic Variants and Their Interactions with Obesity, Physical Activity and Energy Intake on the Risk of Colorectal Cancer: Results from a Nested Case-Control Study in the UK Biobank. Cancers 2020, 12, 2939. [Google Scholar] [CrossRef]
  37. Tan, B.; Qiu, Y.; Zou, X.; Chen, T.; Xie, G.; Cheng, Y.; Dong, T.; Zhao, L.; Feng, B.; Hu, X.; et al. Metabonomics Identifies Serum Metabolite Markers of Colorectal Cancer. J. Proteome Res. 2013, 12, 3000–3009. [Google Scholar] [CrossRef] [Green Version]
  38. Cui, H.; Yang, X.; Wang, Z.; Li, G.; Li, L.; Huo, S.; Zhang, B.; He, R.; Chen, K.; Xu, B.; et al. Tetrahydropalmatine triggers angiogenesis via regulation of arginine biosynthesis. Pharmacol. Res. 2021, 163, 105242. [Google Scholar] [CrossRef]
  39. Mirando, A.C.; Abdi, K.; Wo, P.; Lounsbury, K.M. Assessing the effects of threonyl-tRNA synthetase on angiogenesis-related responses. Methods 2017, 113, 132–138. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Sun, R.; Gu, J.; Chang, X.; Liu, F.; Liang, Y.; Yang, X.; Liang, L.; Tang, D. Metabonomics study on orthotopic transplantion mice model of colon cancer treated with Astragalus membranaceus-Curcuma wenyujin in different proportions via UPLC-Q-TOF/MS. J. Pharm. Biomed. Anal. 2021, 193, 113708. [Google Scholar] [CrossRef] [PubMed]
  41. Reitzer, L. Biosynthesis of Glutamate, Aspartate, Asparagine, L -Alanine, and D -Alanine. EcoSal Plus 2004, 1. [Google Scholar] [CrossRef] [PubMed]
  42. Padidar, S.; Farquharson, A.J.; Williams, L.M.; Kearney, R.; Arthur, J.R.; Drew, J.E. High-Fat Diet Alters Gene Expression in the Liver and Colon: Links to Increased Development of Aberrant Crypt Foci. Dig. Dis. Sci. 2012, 57, 1866–1874. [Google Scholar] [CrossRef] [PubMed]
  43. Lai, Z.; Tsugawa, H.; Wohlgemuth, G.; Mehta, S.; Mueller, M.; Zheng, Y.; Ogiwara, A.; Meissen, J.; Showalter, M.; Takeuchi, K.; et al. Identifying metabolites by integrating metabolome databases with mass spectrometry cheminformatics. Nat. Methods 2018, 15, 53–56. [Google Scholar] [CrossRef] [PubMed]
  44. Chong, J.; Wishart, D.S.; Xia, J. Using MetaboAnalyst 4.0 for Comprehensive and Integrative Metabolomics Data Analysis. Curr. Protoc. Bioinform. 2019, 68, e86. [Google Scholar] [CrossRef] [PubMed]
  45. van den Berg, R.A.; Hoefsloot, H.C.J.; Westerhuis, J.A.; Smilde, A.K.; Van Der Werf, M.J. Centering, scaling, and transformations: Improving the biological information content of metabolomics data. BMC Genom. 2006, 7, 142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Two-dimensional (2D) sPLS-DA of the 4 experimental groups (A) and loading plots of 10 metabolites that are most significant in group separation among the four groups for Component 1 (B) and Component 2 (C). Mice without AOM treatment (control), n = 10/group (AIN or HFD group). Mice with AOM treatment (AIN + AOM or HFD + AOM group), n = 15/group.
Figure 1. Two-dimensional (2D) sPLS-DA of the 4 experimental groups (A) and loading plots of 10 metabolites that are most significant in group separation among the four groups for Component 1 (B) and Component 2 (C). Mice without AOM treatment (control), n = 10/group (AIN or HFD group). Mice with AOM treatment (AIN + AOM or HFD + AOM group), n = 15/group.
Metabolites 11 00448 g001
Figure 2. The interaction effect of diet and AOM on 53 plasma metabolites. Plot of significant interacting metabolites (green circle, 25 of 53 metabolites) by diet vs. AOM interaction, and two-way ANOVA between diet (red circle) and AOM (blue circle) with an FDR-adjusted p < 0.05. Mice without AOM treatment (control) n = 10/group (AIN or HFD group). Mice with AOM treatment, n = 15/group (AIN + AOM or HFD + AOM group).
Figure 2. The interaction effect of diet and AOM on 53 plasma metabolites. Plot of significant interacting metabolites (green circle, 25 of 53 metabolites) by diet vs. AOM interaction, and two-way ANOVA between diet (red circle) and AOM (blue circle) with an FDR-adjusted p < 0.05. Mice without AOM treatment (control) n = 10/group (AIN or HFD group). Mice with AOM treatment, n = 15/group (AIN + AOM or HFD + AOM group).
Metabolites 11 00448 g002
Figure 3. Pathway analysis of 53 differential metabolites and five major metabolic pathways *: (1) citrate cycle; (2) Arg biosynthesis; (3) aminoacyl-tRNA biosynthesis; (4) Ala, Asp, and Glu metabolism; (5) glyoxylate and dicarboxylate metabolism. * Only the pathways with statistical significance are highlighted.
Figure 3. Pathway analysis of 53 differential metabolites and five major metabolic pathways *: (1) citrate cycle; (2) Arg biosynthesis; (3) aminoacyl-tRNA biosynthesis; (4) Ala, Asp, and Glu metabolism; (5) glyoxylate and dicarboxylate metabolism. * Only the pathways with statistical significance are highlighted.
Metabolites 11 00448 g003
Figure 4. Scatterplot analysis of the correlation between (A) the number of ACF and the relative concentrations (quantifier ion peak heights) of dihydrocholesterol; (B) the number of ACF and the relative concentrations (quantifier ion peak heights) of cholesterol. With n = 50 each set, both (A,B) show an inverse association with a significant Spearman correlation (p < 0.0001), respectively.
Figure 4. Scatterplot analysis of the correlation between (A) the number of ACF and the relative concentrations (quantifier ion peak heights) of dihydrocholesterol; (B) the number of ACF and the relative concentrations (quantifier ion peak heights) of cholesterol. With n = 50 each set, both (A,B) show an inverse association with a significant Spearman correlation (p < 0.0001), respectively.
Metabolites 11 00448 g004
Table 1. The effect of HFD feeding and AOM treatment on the abundance of plasma metabolites.
Table 1. The effect of HFD feeding and AOM treatment on the abundance of plasma metabolites.
MetabolitesAINAIN + AOMHFD HFD + AOM
Palmitoleic acid1.00 ± 0.40 a1.11 ± 0.27 a0.19 ± 0.05 b0.30 ± 0.11 b
Myristic acid1.00 ± 0.18 a1.02 ± 0.10 a0.57 ± 0.06 c0.74 ± 0.12 b
Oleic acid1.00 ± 1.46 bc0.73 ± 1.32 c1.99 ± 0.42 ab2.71 ± 0.95 a
Malic acid1.00 ± 0.26 a1.23 ± 0.42 a0.63 ± 0.08 b0.60 ± 0.12 b
Beta-sitosterol1.00 ± 0.35 b0.46 ± 0.11 c2.22 ± 0.65 a0.64 ± 0.10 bc
Oxoproline1.00 ± 0.14 b1.08 ± 0.15 b1.13 ± 0.17 b1.42 ± 0.21 a
Alpha-tocopherol 1.00 ± 0.23 b0.67 ± 0.10 c1.94 ± 0.46 a0.84 ± 0.17 bc
Alpha-ketoglutarate1.00 ± 0.22 a1.15 ± 0.34 a0.72 ± 0.09 b0.63 ± 0.13 b
Fumaric acid1.00 ± 0.14 a1.17 ± 0.32 a0.76 ± 0.05 b0.75 ± 0.11 b
Citric acid1.00 ± 0.24 ab1.08 ± 0.14 a0.68 ± 0.07 c0.88 ± 0.10 b
3-hydroxybutyric acid1.00 ± 0.64 c1.32 ± 0.62 bc2.08 ± 0.50 a1.80 ± 0.57 ab
Saccharic acid1.00 ± 0.18 c1.14 ± 0.21 bc1.49 ± 0.36 a1.33 ± 0.21 ab
Citrulline1.00 ± 0.21 a0.87 ± 0.17 ab0.77 ± 0.13 bc0.64 ± 0.12 c
2-hydroxyglutaric acid1.00 ± 0.13 ab1.12 ± 0.17 a0.86 ± 0.11 b0.87 ± 0.14 b
2-aminobutyric acid1.00 ± 0.37 b0.87 ± 0.36 b0.88 ± 0.24 b1.58 ± 0.54 a
Palmitic acid1.00 ± 0.15 ab1.09 ± 0.14 a0.86 ± 0.08 b0.88 ± 0.13 b
3-ureidopropionate1.00 ± 0.29 b1.35 ± 0.49 b1.19 ± 0.47 b2.02 ± 0.55 a
Myo-inositol1.00 ± 0.24 b0.97 ± 0.13 b1.37 ± 0.32 a1.09 ± 0.27 b
Dihydrocholesterol1.00 ± 0.26 b0.48 ± 0.12 c1.89 ± 0.55 a0.48 ± 0.12 c
Glucuronic acid1.00 ± 0.33 c1.98 ± 0.51 ab1.32 ± 0.37 bc2.48 ± 0.96 a
Lysine1.00 ± 0.37 a0.77 ± 0.22 ab0.57 ± 0.25 b0.58 ± 0.15 b
Cholesterol1.00 ± 0.14 b0.79 ± 0.14 c1.32 ± 0.22 a0.81 ± 0.12 c
Isothreonic acid1.00 ± 0.09 b1.07 ± 0.08 b1.22 ± 0.20 a1.03 ± 0.06 b
Uracil1.00 ± 0.33 b1.04 ± 0.15 b0.95 ± 0.17 b1.29 ± 0.23 a
Ornithine1.00 ± 0.27 b1.39 ± 0.62 b1.05 ± 0.46 b2.04 ± 0.33 a
Succinic acid1.00 ± 0.21 ab1.27 ± 0.47 a0.87 ± 0.11 b0.87 ± 0.24 b
N-acetyl-d-tryptophan1.00 ± 0.22 a1.00 ± 0.26 a0.68 ± 0.22 b0.87 ± 0.18 ab
Pyruvic acid1.00 ± 0.45 ab1.15 ± 0.65 a0.58 ± 0.18 b0.67 ± 0.34 b
2,3-dihydroxybutanoic acid 1.00 ± 0.19 b1.11 ± 0.12 ab1.00 ± 0.12 b1.29 ± 0.24 a
Linoleic acid1.00 ± 0.26 b1.10 ± 0.26 ab0.99 ± 0.20 b1.41 ± 0.50 a
Beta-alanine1.00 ± 0.32 a0.97 ± 0.33 a0.60 ± 0.21 b0.86 ± 0.23 ab
Maleimide1.00 ± 0.13 b0.88 ± 0.14 b1.26 ± 0.14 a0.85 ± 0.14 b
Isocitric acid1.00 ± 0.17 a1.15 ± 0.23 a0.77 ± 0.11 b1.08 ± 0.14 a
Tyrosine1.00 ± 0.19 a1.12 ± 0.20 a0.75 ± 0.22 b1.06 ± 0.14 a
Xylitol1.00 ± 0.54 b1.58 ± 0.42 a0.90 ± 0.23 b2.03 ± 0.65 a
Nicotinamide1.00 ± 0.26 b0.92 ± 0.20 bc1.58 ± 0.24 a0.76 ± 0.15 c
Putrescine1.00 ± 0.34 a0.77 ± 0.25 ab0.70 ± 0.10 b0.70 ± 0.12 b
Threonic acid1.00 ± 0.15 b1.16 ± 0.24 ab1.39 ± 0.22 a1.01 ± 0.24 b
Methanolphosphate1.00 ± 0.12 b0.69 ± 0.11 c1.18 ± 0.15 a0.65 ± 0.09 c
Threonine1.00 ± 0.23 a0.76 ± 0.13 b0.74 ± 0.16 b0.75 ± 0.12 b
Serine1.00 ± 0.22 ab0.82 ± 0.14 b0.80 ± 0.21 b1.02 ± 0.14 a
Glycine1.00 ± 0.32 a0.82 ± 0.12 ab0.78 ± 0.09 b0.80 ± 0.12 b
Salicylic acid1.00 ± 0.16 ab1.15 ± 0.21 ab0.95 ± 0.16 b1.20 ± 0.19 a
Erythritol1.00 ± 0.13 b1.13 ± 0.07 a1.03 ± 0.10 b1.08 ± 0.06 ab
Sophorose1.00 ± 0.32 ab0.88 ± 0.29 b1.28 ± 0.44 a0.78 ± 0.18 b
Thymidine1.00 ± 0.21 ab1.03 ± 0.24 a0.79 ± 0.13 b1.20 ± 0.20 a
Aminomalonate1.00 ± 0.34 a0.76 ± 0.14 bc0.93 ± 0.25 ab0.67 ± 0.11 c
Arachidonic acid1.00 ± 0.11 a0.71 ± 0.14 b1.06 ± 0.16 a0.60 ± 0.10 b
Aconitic acid1.00 ± 0.35 ab1.09 ± 0.32 a0.75 ± 0.30 b1.14 ± 0.25 a
Phenylalanine1.00 ± 0.19 ab0.97 ± 0.24 ab0.81 ± 0.31 b1.13 ± 0.20 ab
Aspartic acid1.00 ± 0.33 ab1.01 ± 0.28 ab0.72 ± 0.20 b1.13 ± 0.27 a
Methionine1.00 ± 0.24 ab0.94 ± 0.26 ab0.73 ± 0.30 b1.08 ± 0.23 a
Glycerol-alpha-phosphate1.00 ± 0.29 a0.79 ± 0.08 b0.89 ± 0.14 ab0.79 ± 0.10 b
Values are means ± SDs, n = 10/group (AIN or HFD group) and n = 15/group (AIN + AOM or HFD + AOM group). Data from the HFD, AIN + AOM, and HFD + AOM groups were converted to fold changes compared to the AIN group. For a given metabolite, if two values share at least one letter, then the difference between them is not statistically significant. However, if they do not have a letter in common, then the difference between them is statistically significant, p < 0.05 adjusted by the FDR method.
Table 2. The major metabolic pathways * involved in 53 differential metabolites (between experimental groups).
Table 2. The major metabolic pathways * involved in 53 differential metabolites (between experimental groups).
KEGG PathwayNumber of Metabolites Identifiedp *Impact **
Citrate cycle (TCA cycle)7<0.00010.35
Arg biosynthesis5<0.0030.29
Aminoacyl-tRNA biosynthesis8<0.0050.17
Ala, Asp, and Glu metabolism6<0.010.27
Glyoxylate and dicarboxylate metabolism6<0.030.20
* Only the pathways with statistical significance are listed, and p-values are obtained by the over-representation analysis and adjusted by Holm and FDR methods. ** Impact is the pathway impact score obtained by the pathway topology analysis.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zeng, H.; Umar, S.; Liu, Z.; Bukowski, M.R. Azoxymethane Alters the Plasma Metabolome to a Greater Extent in Mice Fed a High-Fat Diet Compared to an AIN-93 Diet. Metabolites 2021, 11, 448. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo11070448

AMA Style

Zeng H, Umar S, Liu Z, Bukowski MR. Azoxymethane Alters the Plasma Metabolome to a Greater Extent in Mice Fed a High-Fat Diet Compared to an AIN-93 Diet. Metabolites. 2021; 11(7):448. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo11070448

Chicago/Turabian Style

Zeng, Huawei, Shahid Umar, Zhenhua Liu, and Michael R. Bukowski. 2021. "Azoxymethane Alters the Plasma Metabolome to a Greater Extent in Mice Fed a High-Fat Diet Compared to an AIN-93 Diet" Metabolites 11, no. 7: 448. https://0-doi-org.brum.beds.ac.uk/10.3390/metabo11070448

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