Metabolic and gut microbial characterization of obesity-prone mice

Feb 22, 2019 - Obesity is characterized with high heterogeneity due to genetic abnormality, energy imbalance and/or gut dysbiosis. Obesity-prone (OP) ...
0 downloads 0 Views 889KB Size
Subscriber access provided by Macquarie University

Article

Metabolic and gut microbial characterization of obesity-prone mice under high-fat diet Yu Gu, Can Liu, Ningning Zheng, Wei Jia, Weidong Zhang, and Houkai Li J. Proteome Res., Just Accepted Manuscript • DOI: 10.1021/acs.jproteome.8b00945 • Publication Date (Web): 22 Feb 2019 Downloaded from http://pubs.acs.org on February 23, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Metabolic and gut microbial characterization of obesity-prone mice under high-fat diet Yu Gu1, Can Liu2, Ningning Zheng1, Wei Jia3,4, Weidong Zhang1,5*, Houkai Li1* 1

Institute of Interdisciplinary Integrative Medicine Research, Shanghai University of

Traditional Chinese Medicine, Shanghai, China 2

Department of Biochemistry and Molecular Biology, Bengbu Medical College, Anhui

Province, China 3

Shanghai Key Laboratory of Diabetes Mellitus and Center for Translational Medicine,

Shanghai Jiao Tong University Affiliated Sixth People's Hospital, Shanghai, China 4 University 5

of Hawaii Cancer Center, Honolulu, Hawaii, USA

Department of Phytochemistry, College of Pharmacy, Second Military Medical

University, Shanghai, China

*Corresponding

author:

Weidong Zhang: Address: No. 1200 Cai Lun Road, Pudong New District, Shanghai, China; No.

325

Guo

He

Road,

Yangpu

District,

Shanghai,

China.

E-mail:

[email protected].

Houkai Li: Address: No. 1200 Cai Lun Road, Pudong New District, Shanghai, China. Email: [email protected].

1

ACS Paragon Plus Environment

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ABSTRACT Obesity is characterized with high heterogeneity due to genetic abnormality, energy imbalance and/or gut dysbiosis. Obesity-prone (OP) and -resistant (OR) phenotypes are frequently observed in rodents even under high-fat diet (HFD). However, the underlying mechanisms are largely unknown. Male C57BL/6J mice were fed with chow or HFD for 8 weeks. OP and OR mice were defined based on body weight gain, and an integrated serum metabolic and gut microbial profiling was performed by gas chromatography-mass spectroscopy-based metabolomic and pyrosequencing of 16S rDNA of cecum contents. Sixty differential metabolites were identified in comparisons among Con, OP and OR groups, in which 27 were OP-related. These differential metabolites are mainly involved in glycolysis, lipids and amino acids metabolism, and TCA cycle. Meanwhile, OP mice had distinct profile in gut microbiota compared to OR or Con mice, which showed reduced ratio of Firmicutes to Bacteroidetes and increased Proteobacteria. Moreover, gut microbial alteration of OP mice was correlated with the changes of the key serum metabolites. OPenriched Parasutterella from Proteobacteria phylum correlated to most of metabolites, suggesting it was essential in obesity. OP mice are distinct in metabolic and gut microbial profiles, and OP-related metabolites and bacteria are of significance for understanding obesity development.

KEYWORDS Metabolic profiling; Gut microbiota; Obesity-prone; 16S rDNA; Obesity-resistant

2

ACS Paragon Plus Environment

Page 2 of 33

Page 3 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

1 INTRODUCTION Obesity is the critical risk factor for metabolic diseases including type 2 diabetes, cardiovascular disease, metabolic syndrome (encompassing hypertension, dyslipidemia, and insulin resistance), and non-alcoholic fatty liver disease

1, 2.

The etiology of obesity

involves complicated interactions between genetic and environmental factors such as dietary nutrients and physical activity 3. In addition, numerous investigations have indicated that gut dysbiosis is intricately involved in the development of obesity

4, 5

and

obesity-related metabolic diseases 6, 7. Meanwhile, obesity is characterized with substantial heterogeneity either in humans or rodent models 8, 9. For instance, many maternal twins even with highly identical genetics exhibit significant differences in body weight and variable susceptibility to metabolic diseases

9-11.

Similarly, the extent of body weight gain of rodents with same genetic

background usually vary greatly when fed with identical high-fat diet (HFD). Animals that are sensitive to HFD-induced obesity are defined as obesity-prone (OP) phenotype, while those insensitive ones are defined as obesity-resistant (OR)

12.

Substantial evidence has

revealed that the development of OP or OR phenotype is related with the difference in neural regulatory pathways 13-15. On the other hand, our previous investigation has observed significant differences in both hepatic gene expression and metabolic profiles of liver, serum, and urine samples between OP and OR rats by using combined transcriptomics and metabolomics. We identify some gut microbiota-derived metabolites such as phydroxyphenylacetic acid that are of significant difference between OP and OR phenotypes, suggesting that the development of OP phenotype might be associated with alteration of gut microbiota 14. Consistently, a recent study shows that bacteria of Bifidobacterium genus 3

ACS Paragon Plus Environment

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

from Actinobacteria phylum were significantly enriched in OP mice compared to that of OR ones

16.

Therefore, OP phenotype is associated with either the metabolic or gut

microbial alteration, however, little is known about the crosstalk between host metabolism and gut microbiota in regard to the differences between OP and OR phenotypes. In our current study, both OP and OR mice were obtained based on the extent of body weight gain among a group of mice under HFD feeding. Then, the serum metabolic and gut microbial profiling was performed by using GC/MS-based untargeted metabolomics and 16S rDNA sequencing on the bacterial genomic DNA from cecum contents, respectively. Our results showed distinct separation in either serum metabolic or gut microbial profile between OP and OR mice. Moreover, OP-related differential metabolites were identified which mainly included fatty acids, amino acids and their metabolites, formic acid, fructose, phosphoethanolamine, glycolic acid and etc, as well as the involved metabolic pathways. Meanwhile, gut microbial profiling revealed that the compositional alteration of gut bacteria was not only correlated with the body weight phenotypes, but also showed high correlation with some differential metabolites which highlights the crosstalk between host metabolism and gut microbiota in mediating the different phenotypes.

2 EXPERIMENTAL PROCEDURES 2.1 Animal treatment and sample collection Six-week-old male C57BL/6J mice (n=22) were purchased from Shanghai Laboratory Animal (Shanghai, China), and housed in a regulated barrier system facility at 23°C–24°C with 60%±10% relative humidity and a 12-h light/dark cycle. All of the mice were fed a standard chow diet during an initial 1-week acclimation period. Then, 6 mice were 4

ACS Paragon Plus Environment

Page 4 of 33

Page 5 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

randomly selected for control group with continued standard chow diet (average initial body weight 21.5 ± 0.5 g), while the rest 16 mice were changed into HFD feeding (average initial body weight 21.0 ± 1.0 g). The detailed composition and percentage of both standard and HFD diet (Nantong Trophy Experimental Animal Feed Co., Ltd) was shown in Supplemental Table S1. After 8 weeks feeding, mice in the HFD group were ranked base on the extent of their body weight gain. Six mice with the highest or lowest body weight gain (about 1/3 in total) were defined as OP or OR mice respectively according to previous publication

12.

Finally, mice were sacrificed after anesthesia with 10% chloral hydrate

intraperitoneally for collecting liver and white adipose tissue, blood and cecum contents after an overnight fasting (16 hours). Serum samples were obtained by centrifuging blood with 4000 rpm at 4°C. All the animal experiments were admitted by Animal Experiment of Shanghai University of Traditional Chinese Medicine (Shanghai, China), and the protocol was approved by the institutional Animal Ethics Committee.

2.2 Lipids measurement Liver samples were subjected to isopropanol extraction before triglycerides (TG) and total cholesterol (TC) measurement according to the previous method 17. 50-mg liver tissue was homogenized in 1-mL chloroform: methanol (2: 1) before agitation at 200 rpm for 2 hours at room temperature. Then 20% volume of 0.9% NaCl (w/v) was added and subsequently centrifugated at 2000 rpm for 10 min. The lower solvent layer was mixed with 1% TritonX100 before being dried under N2. The sample was dissolved with 100-μL isopropanol. Finally, the liver TG and TC contents were measured with the TG or TC assay kits (jiancheng Bioengineer Institute, Nanjing, China). Consistently, the serum TG and TC 5

ACS Paragon Plus Environment

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

contents were measured with the same method.

2.3 Serum sample preparation for metabolic profiling Serum metabolic profiling was performed by using GC/MS-based metabolomic approach as previously described. Briefly, 20-μL aliquot of each serum sample was added with 10-μL internal standard (50 μg/mL phenylalanine-13C9-15N and dulcitol). Then, 80μL cold methanol-chloroform (3:1, v/v) was added into each sample, mixed for 60 s, and centrifuged at 14000 g for 15 min at 4 °C. 80-μL aliquot of supernatant was collected and transferred to the GC vial for dryness under N2. The dried samples were dissolved in 30 μL of methoxylamine hydrochloride in pyridine (20 mg/mL), and incubated at 70 °C for 60 min. Samples were derivatized by adding 30-μL BSTFA (1% TMCS) and incubating at 70 °C for 60 min. Quality control (QC) sample was prepared by mixing the same volume of serum from each sample, and prepared in parallel.

2.4 GC/MS analysis The derivatized samples were analyzed in an Agilent 7890A gas chromatography system coupled to an Agilent 5975C MSD system (Agilent Technologies, Santa Clara, CA, USA) and an Agilent J&W HP-5MS fused-silica capillary column (30 m × 0.25 mm × 0.25 μm). The analysis was performed using helium as the carrier gas (constant flow rate = 1 mL/min). An injection volume was 1 μL. Solvent delay time was 5.5 min, and collision energy was 70 eV. The initial oven was heated to 70 °C for 2 min, ramped to 160 °C at 6 °C/min, to 240 °C at 10 °C/min, and then increased at a rate of 20 °C/min to a final temperature of 300 °C for 6 min. The injector, transfer line, and ion source (electron impact) 6

ACS Paragon Plus Environment

Page 6 of 33

Page 7 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

were heated to 250 °C, 260 °C, and 230 °C, respectively. Mass data were acquired in fullscan mode (m/z 50−600). All the samples were analyzed in a random sequence.

2.5 Data processing and statistical analysis The raw GC/MS data were preprocessed according to the described method including peak picking, alignment, deconvolution, and further processing 18. The detailed structure identification of metabolites was performed with the following method. AMDIS software was used to deconvolute mass spectra, and the purified mass spectra were automatically matched with an in-housed standard library including retention time and mass spectra, Golm Metabolome Database, and Agilent Fiehn GC/MS Metabolomics RTL Library. Then the normalized identified metabolite data were imported into the SIMCA software, version 14.1 (Umetrics, Umeå, Sweden), which was used to perform the principal component analysis (PCA) and partial least-squares discriminant analysis (PLS-DA). The metabolites were filtered by variable influence on projection (VIP) selection using values obtained from the PLS-DA model and the filtering conditions VIP > 1 and p < 0.01, with p-values determined using two-tailed Student’s t-tests.

2.6 Bacterial genomic DNA extraction of cecum contents To perform 16S rDNA sequencing, we first carefully pooled the cecum contents from the 6 mice of each group (Con, OP and OR). Then, the pooled samples were divided into 3 replicates for bacterial genomic DNA extraction using the QIAamp fast DNA Stool Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol 19. Briefly, 200 mg of cecum contents were added with 1 mL inhibit EX buffer and vibrated. After heating 7

ACS Paragon Plus Environment

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

the suspension at 70 ℃

Page 8 of 33

and vortex, sample was centrifuged for 1 min. Then the

supernatant was added into the 1.5 mL microcentrifuge tube with Proteinase K, and mixed thoroughly with 200 μL buffer AL before incubation at 70 ℃ for 10 min. 200 μL of ethanol was added to lysate and mix. 600 μL of lysate was applied to QIAmap spin column and then centrifuged. QIAmap spin column was placed in a new 2 mL collection rube. Then 500 μL butter AW1 and AW2 was added into spin column, respectively, with repeating above operation. Subsequently, 200 μL buffer ATE was pipetted onto QIAamp membrane to obtain purified DNA.

2.7 16S rDNA pyrosequencing Purified bacterial genomic DNA was used as a template for the amplification and DNA sequencing of the V3 region of 16S rRNA gene. The used primers were 338F (5’CCTACGGGAGGCAGCAG-3’) and 518R (5’- ATTACCGCGGCTGCTGG-3’). The PCR amplification with bar-coded primers, pyrosequencing of PCR amplicons, and quality control of raw data were performed as described previously 20.

2.8 Bioinformatics and multivariate statistical analysis The acquired valid and representative sequences of each sample were compared Greengenes database using the nearest alignment space termination algorithm

21,

and

constructed a neighbor-joining tree with ARB 22. Operational taxonomic units (OTUs) were delineated at 97% similarity level with Mothur software. The representative sequence of each OTU was selected with the most abundance and subjected to RDP classifier for taxonomical assignment with a bootstrap cutoff of 60% 8

ACS Paragon Plus Environment

23, 24.

Weighted Fast UniFrac

Page 9 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

principal coordinate analysis (PCoA), was performed with the phylogenetic tree constructed by each OTU generated with QIIME 25. The comparisons of bacteria relative abundance at phylum level between groups were performed with student’s t test after the normality analysis of all data. The top 30 genera were analyzed in R 3.2.4 with pheatmap packages.

2.9 Correlation analysis between bacterial taxonomy and differential metabolites Correlation analysis between bacterial taxonomy at genus level and differential metabolites was performed with Spearman’s correlation coefficient (SCC). The significant correlation between metabolites and microbial genera was performed with the criteria of both p < 0.05 and SCC > 0.7 and < -0.7, which was visualized in a cross-correlation maps including positive (red line) or negative (blue line) relationship.

3 RESULTS 3.1 Characterization of OP phenotype in HFD feeding mice After 8-week HFD or chow diet feeding, the phenotypic characteristics, body weight gain, organ weight and serum biochemistry were examined and shown in Figure 1. First, we observed that mice in HFD group showed different responses towards HFD feeding, and therefore OP and OR mice were determined on the basis of their body weight gain (Fig. 1A). Consistently, OP mice showed heavier WAT mass and WAT index compared to OR or Con mice (Fig. 1B-C), as well as the liver weight (Fig. 1D) and liver TG (Fig. 1F). The 5 (Fig. 1G). In addition, the levels of both serum TG and TC were significantly higher in OP mice compared to Con group, but of no difference between OP and OR mice (Fig. 1H-I), 9

ACS Paragon Plus Environment

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

as well as the comparable levels of liver TC among groups (Fig. 1G). Accordingly, the OP mice were of different characters compared to either OR or Con mice.

Fig. 1. Phenotypes of the control (Con, fed standard chow diet), obesity-prone (OP, fed high-fat diet), and obesity-resistant (OR, fed high-fat diet) mice (n = 6 per group). (A) Body weight gain. (B) WAT (white adipose tissue) weight. (C) WAT index (the weight of the WAT accounts for the body weight). (D) Liver weight. (E) Liver index (the weight of the liver accounts for the body weight). (F) Liver TG (G) Liver TC. (H) Serum TG. (I) Serum TC. The p value was calculated by two-tailed Student’s t-test. All data are presented as mean ± SEM.

3.2 Distinct metabolic changes in OP mice To elucidate the metabolic character of OP mice, the GC/MS-based untargeted metabolic profiling was performed on serum samples. First, 182 endogenous metabolites 10

ACS Paragon Plus Environment

Page 10 of 33

Page 11 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

were determined in serum samples. Then, an unsupervised PCA was performed to visualize the general differences among samples. The PCA scatter plot showed that samples from the Con group were distinctly separated from the rest, whereas samples of OP mice also showed separation trend with OR mice (Fig. 2A). The supervised PLS-DA model showed clear separation among groups (Fig. 2B). Meanwhile, models of PCA and PLS-DA between either OP or OR and Con groups were conducted respectively (Supplemental Figure S1A-F). The evaluation data for the quality of PCA and PLS-DA models were summarized in Supplemental Table S2.

Fig. 2. Analyses of serum metabolites from Con, OR and OP group based on GC/MS (n = 6 per group). The SIMCA-P-derived (A) PCA and (B) PLS-DA plots among the Con, OR, and OP mice. Differential metabolites based on PLS-DA results using VIP>1 and P1 at multivariate statistical analysis or p1.2.

3.7 Correlation analysis between obesity-related phenotypes, bacterial genera and differential metabolites Given the fact that the development of obesity phenotype is associated with both the imbalance of metabolism and gut dysbiosis 36, we therefore are interested in correlating the 27 differential metabolites (clusters of I, III and V) with the 14 bacterial genera that were obviously altered in OP group compared to either Con or OR groups by spearman’s correlation analysis. The differential metabolites in cluster I (12 metabolites), and III (6 metabolites) were acquired by comparing OP with either Con or OR mice, while cluster V (9 metabolites) represented the commonly altered differential metabolites between OP and Con/OR groups. The correlation analysis showed that genera, Anaerotruncus, Butyricimonas, Parasutterella, Butyricicoccus, Rikenella, Xylanibacter, were positively correlated the differential metabolites from both cluster I and V, as well as the obesityrelated

phenotypes.

On

the

contrary,

genera,

19

ACS Paragon Plus Environment

Alistipes,

Oscillibacter,

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 33

Lachnospiracea_incertae_sedis, Bacteroides and Wolinella, were negatively correlated with the metabolites from both cluster I and V, as well as the obesity-related phenotypes. The 14 selected genera showed minor to moderate correlation with the 6 differential metabolites in cluster III, except for the negative correlation between Alistipes with the 6 metabolites (Supplemental Figure S3). Then, a further cross-correlation map was conducted with the criteria of SCC > 0.7 or < -0.7 and p < 0.05, in which 19 differential metabolites and 11 bacterial genera were included (Fig. 5). Parasutterella genus (belonging to Proteobacteria phylum) were positively correlated with metabolites such as 2-hydroxyisovaleric acid, linoleic acid, glycerol-2-phosphate, glyceryl monooleate, uridine, glutamic acid, stearic acid, adenine, methylphosphate, cis-5,8,11-eicosatrienoic acid, whereas Wolinella genus (belonging to Proteobacteria phylum) were negatively correlated with metabolites such as 2hydroxyisovaleric acid, pantothenic acid, cis-5,8,11-eicosatrienoic acid, and DHA, suggesting that bacteria within these two types of genera might perform different roles in OP mice. Six genera within Bacteroidetes phylum were included in this correlation analysis such

as

Bacteroides,

Butyricimonas,

Alistipes,

Rikenella,

Xylanibacter

and

Lachnospiracea_incertae_sedis. Interestingly, only a small number of metabolites were correlated with these 6 genera in Bacteroidetes phylum. For example, Bacterroides were negatively 2-hydroxyisovaleric acid. Butyricimonas genus were positively correlated with alanine, but negatively correlated with hypotaurine. Lachnospiracea_incertae_sedis genus were negatively correlated with glutamic acid. In addition, Alistipes genus was negatively correlated with glycolic acid and 2-hydroxyadipic acid. Besides, 3 genera within Firmicutes were included such as Oscillibacter, Butyricicoccus and Anaerotruncus. 20

ACS Paragon Plus Environment

Page 21 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Oscillibacter genus was negatively correlated with 4 metabolites including 4hydroxylbutyric acid, pantothenic acid, methylphosphate, and cis-5,8,11-eicosatrienoic acid. Butyricicoccus genus was positively correlated with 5 metabolites including palmitic acid, linoleic acid, oleic acid, stearic acid and glyceryl monooleate. Finally, only DHA was showed positive correlation with bacteria in Anaerotruncus genus.

Fig. 5. Spearman correlations between differential metabolites and genera in OP mice. Correlation linkage chart with positive correlations indicated by red lines and negative correlations indicated by blue lines. Spearman’s correlation coefficient less than -0.7 or more than 0.7 with p < 0.05 was selected. The selected genera were shown in cluster with its phylum and family information.

4 DISCUSSION In our current study, we first produced OP and OR phenotype in male C57BL/6J mice under 8-weeks HFD feeding, which is consistent with previous report 37. The gut microbial 21

ACS Paragon Plus Environment

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

and metabolic profiling revealed that OP mice were distinguished from either OR or Con mice. The 27 differential metabolites between OP and Con/OR mice were mainly involved in pathways of FA metabolism, TCA cycle, urea cycle and oxidative stress. Moreover, the correlation analysis showed that 11 genera were positively or negatively correlated with these OP-related differential metabolites, as well as obesity phenotypes including body weight gain, WAT weight and index, liver weight and TG levels. Obesity results from the interaction between genetic and environmental factors 38, and numerous studies have revealed the critical roles of energy imbalance 9, 39, as well as neural dysregulation in determining obesity formation

15, 40.

Moreover, recent advances in gut

microbiota have highlighted that gut dysbiosis is the common basis for various diseases such as obesity and obese-related metabolic diseases3, 16. Our previous investigation shows that OP and OR rats have significantly different metabolism under HFD feeding, in which some gut microbial-related metabolites are detected 14. In this study, we demonstrated that OP mice had very different metabolic profile compared to that of either Con or OR mice. Although the body weight gain was not statistically significant between OP and Con groups, OP mice showed significant increases in WAT weight (index), serum and liver TG, and serum TC, as well as higher extent of body weight gain than Con mice. Therefore, the OPrelated metabolic changes between OP and Con groups are mainly associated with these obese phenotypes, in addition to probable dietary impact. Subsequent analysis identified 21 differential metabolites between OP and OR mice, in which 6 of them were commonly altered between Con and OR mice and were excluded for further analysis. There were 6 differential metabolites in cluster III including arachidonic acid, 2-hydroxyadipic acid, hypotaurine, erithronic acid, fructose and glycolic acid, which 22

ACS Paragon Plus Environment

Page 22 of 33

Page 23 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

are of higher concentration in OP mice than OR. Since OP and OR mice were fed with identical HFD, the acquired differential metabolites between OP and OR mice could be attributed to their metabolic differences, instead of dietary impact. Meanwhile, 9 differential metabolites stood out that were significantly up-regulated in OP mice compared to either Con or OR mice including DHA, cis-5,8,11-eicosatrienoic acid, methylphosphate, adenine, stearic acid, phosphoethanolamine, pantothenic acid, glutamic acid, and fumaric acid. Previous studies have shown that n-6 PUFA can stimulate adipogenesis by upregulating peroxisome proliferator-activated receptor γ (PPARγ) activity and sterol regulatory element-binding protein-1c (SREBP-1c) expression resulting to lipogenesis and TG storage in adipose tissue 41, 42, whereas opposite effects are observed by n-3 PUFA such as ARA 42, 43. Therefore, our current results suggested disordered fatty acid metabolism in OP mice. In addition to the OP-related metabolic characters, our results showed that OP mice had distinct feature in gut microbiota compared to either Con or OR mice. Previously, a large number of publications have confirmed that gut dysbiosis plays critical roles in obesity development either in humans or animal models

36, 44.

Although the composition of gut

microbiota is easily affected by diets, our current results revealed that the compositional alteration of gut microbiota in OP mice were diets-independent because of the observed significant differences in gut microbiota between OP and OR/Con mice. The decreased ratio of Firmicutes to Bacteroidetes is predominantly observed in obese subjects 35. In our current study, we found that OP mice showed decreased ratio of Firmicutes to Bacteroidetes compared to either OR or Con mice. Nevertheless, contradictory results are also reported in other studies

45, 46.

Proteobacteria is a major phylum of gram-negative 23

ACS Paragon Plus Environment

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

bacteria, which consists a large number of pathogens with an outer membrane composed of lipopolysaccharides (LPS). Obesity is also characterized as endoxemia because of the release of LPS into circulation

47.

Accordingly, bacteria in Proteobacteria phylum are

usually enriched in obesity and metabolic diseases 20. Our current results showed that the relative abundance of Proteobacteria phylum mildly increased in OP mice. Moreover, bacteria of Parasutterella genus were obviously enriched in OP mice compared to both OR and Con group, as well as some unclassified genera within Proteobacteria phylum. Consistently, relative higher abundance of Parasutterella has also been reported in OP rats 48.

Although the metabolic alteration and gut microbial contribution to obesity formation have been well-investigated, the correlation between the altered metabolism and gut dysbiosis in obese phenotype has been rarely addressed, let alone between OP and OR phenotypes. In our current study, the correlation analysis was evaluated between the obeserelated phenotypes, 14 genera of OP-related bacteria and 27 differential metabolites from cluster I, III and V. Interestingly, the relative abundance of Alistipes genus from Bacteroidetes phylum was negatively correlated with most of the metabolites and obeserelated phenotypes. Since Alistipes genus belongs to acetic acid-producing genera 49, and decreased levels of SCFAs are frequently observed in obese subjects 50, our current results suggested that the formation of OP phenotype might be associated with the abundance of SCFAs-producing bacteria. Butyricicoccus genus are SCFAs-producing bacteria, which are usually reduced in obese-related disorders and enriched by various interventions against obesity

51.

In our current study, we found that the relative abundance of Butyricicoccus

genus was moderately enriched in OP mice. Moreover, the positively correlated 24

ACS Paragon Plus Environment

Page 24 of 33

Page 25 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

metabolites with Butyricicoccus genus were mainly middle chain fatty acids such as palmitic acid, linoleic acid, oleic acid, and stearic acid. These types of fatty acids usually play vicious roles such as aggravating insulin resistance and inflammation in metabolic disorders, excluding oleic acid

41, 52-54.

Therefore, our current results suggested that the

physiological functions of bacteria in Butyricicoccus genus were much more complicated than expected. Meanwhile, the relative abundance of Butyricimonas genus was reduced in OP mice compared to OR mice. Butyricimonas genus are also butyric acid-producing bacteria belong to Bacteroidetes phylum 55, which are usually reduced in obesity and obeserelated metabolic disorders

56,

and enriched by various anti-obesity treatments such as

pomegranate extract accompanied with body weight reduction and gut microbiome recovery

57.

Interestingly, the increased Butyricimonas genus were positively correlated

with alanine, but negatively with hypotaurine in OP mice. Hypotaurine is an intermediate in the biosynthesis of taurine from cysteine. Recent study reveals that infants born to obese mothers have enriched microbiome associated with a number of metabolic pathways such as taurine and hypotaurine pathway, suggesting that the susceptibility to obesity formation of these infants may be associated with the altered gut microbiome and taurine metabolism pathway 31. In addition, hypotaurine is one of the endogenous neurotransmitters of glycine receptor, and inhibition of glycine receptor-1 in the dorsal vagal complex can reduce food intake and body weight gain in rats

58.

Our current results, on one hand, consistently

revealed that increased hypotaurine might be involved in OP formation. On the other hand, the negative correlation between hypotaurine and butyric acid-producing bacteria in Butyricimonas genus suggested that bacteria in Butyricimonas genus might have profound ways in affecting host metabolism, in addition to butyric acid production. 25

ACS Paragon Plus Environment

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Proteobacteria phylum are gram-negative and mainly pathogens. We observed reduced abundance of Wolinella and enriched Parasutterella genera in OP mice. The role of bacteria of Wolinella genus in obesity is poorly known yet. Nevertheless, it is noticeable that the relative abundance of Parasutterella genus was significantly enriched in OP mice, which was also positively correlated with a variety of metabolites and obesity-related phenotypes. The increased abundance of Parasutterella genus has been observed in obesity-related conditions. For example, Li R et al investigated the antibiotic exposureinduced body weight gain and gut microbiome alteration. Their results indicate that administration of florfenicol or azithromycin lead to obviously increased adipogenesis and body weight gain, which are accompanied with increases in abundance of Alistipes, Desulfovibrio, Parasutterella and Rikenella genera

59.

Parasutterella genus has been observed in OP rats

Meanwhile, higher abundance of

48.

Accordingly, these observations

suggest that enriched bacteria in Parasutterella genus might play important role in obesity formation. In addition, the abundance of Parasutterella genus was positively correlated with 10 metabolites in our current study, and most of the metabolites are associated with obesity phenotypes. For example, glutamic acid takes part in obesity formation and even toxic to body 27, 28. Moreover, recent study reveals that higher levels of serum glutamic acid in obese subjects, which is due to reduced abundance of B. thetaiotaomicron, a glutamatefermenting commensal. Interestingly, B. thetaiotaomicron administration not only decreases plasma glutamate concentration, but also alleviates diet-induced obesity in mice 60.

Accordingly, the positive correlation between serum glutamate and Parasutterella

genus in our current report suggested that bacteria from Parasutterella genus might be also involved in elevating serum levels of glutamate, which is worthy further investigation. 26

ACS Paragon Plus Environment

Page 26 of 33

Page 27 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

Altogether, our current report highlights that OP mice are characterized with distinct serum metabolic and gut microbial profiles. The altered metabolic pathways are mainly involved glycolysis, lipid metabolism, amino acid metabolism, and TCA cycle. Moreover, gut microbial alteration of OP mice is correlated with the changes of serum metabolites which provides novel information in regard to the gut microbial contribution to obesity phenotypes. Our current study provides novel evidence of some OP-related bacteria genera and potential biomarkers for OP phenotype. Future studies are warranted to reveal the underlying mechanisms that links between metabolic and gut microbial profiles.

ASSOCIATED CONTENT Supporting Information Figure S1. The SIMCA-P-derived PCA and PLS-DA plots for the GC/MS spectral analysis of serum samples (n = 6 per group). Figure S2. Alpha diversity of the gut microbiota in the among Con, OR and OP groups (n = 3 per group). Figure S3. Association map of the threetiered analysis integrating the phenotype, genera and differential metabolites. Table S1: The different substances and percentage between standard diet and HFD diet. Table S2: Summary of the GC/MS data sets used in PCA and PLS-DA modeling.

AUTHOR INFORMATION

Corresponding Author Weidong Zhang: Address: No. 1200 Cai Lun Road, Pudong New District, Shanghai, China; No. 325 Guo He Road, Yangpu District, Shanghai, China. E-mail: [email protected]. Houkai Li: Address: No. 1200 Cai Lun Road, Pudong New District, Shanghai, China. 27

ACS Paragon Plus Environment

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

E-mail: [email protected]. ORCID Houkai Li: 0000-0003-2846-7895 Conflict of Interest The authors declare no competing financial interest. Author Contributions Yu Gu was responsible for data analysis, figures and manuscript writing. Ningning Zheng and Can Liu conducted animal experiment. Wei Jia helped in project design. Weidong Zhang and Houkai Li designed the study and manuscript revision.

ACKNOWLEDGMENTS This work was funded by the National Natural Science Foundation of China (No. 81873059 & 81673662), National Key Research and Development Program of China (No. 2017YFC1700200), the Program for Professor of Special Appointment (Eastern Scholar) and Shuguang Scholar (16SG36) at Shanghai Institutions of Higher Learning from Shanghai Municipal Education Commission.

ABBREVIATIONS OP: obesity-prone; OR: obesity-resistant; HFD: high-fat diet; WAT: adipose tissue; TG: triglycerides; TC: total cholesterol; SCC: Spearman’s correlation coefficient; FAs: fatty acids; SFAs: long chain fatty acids; PUFAs: polyunsaturated fatty acids; DHA: cis4,7,10,13,16,19-docosahexaenoic acid; ARA: arachidonic acid; LA: linoleic acid; MUFA: monounsaturated fatty acid; GSH: glutathione; ROS: reactive oxygen species; LPS: 28

ACS Paragon Plus Environment

Page 28 of 33

Page 29 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

lipopolysaccharide.

REFERENCES 1. Kalarchian, M. A.; Marcus, M. D.; Levine, M. D.; Courcoulas, A. P.; Pilkonis, P. A.; Ringham, R. M.; Soulakova, J. N.; Weissfeld, L. A.; Rofey, D. L., Psychiatric disorders among bariatric surgery candidates: Relationship to obesity and functional health status. American Journal of Psychiatry 2007, 164, (2), 328-334. 2. Ghazi, A. A.; Amirbaigloo, A., Hypocalcemia and osteomalacia after bariatric surgery. Clin Cases Miner Bone Metab 2017, 14, (2), 227-229. 3. Cox, A. J.; West, N. P.; Cripps, A. W., Obesity, inflammation, and the gut microbiota. Lancet Diabetes Endocrinol 2015, 3, (3), 207-15. 4. Choquet, H.; Meyre, D., Genomic insights into early-onset obesity. Genome Medicine 2010, 2. 5. Vinolo, M. A.; Rodrigues, H. G.; Nachbar, R. T.; Curi, R., Regulation of inflammation by short chain fatty acids. Nutrients 2011, 3, (10), 858-76. 6. Kang, Y.; Cai, Y., The development of probiotics therapy to obesity: a therapy that has gained considerable momentum. Hormones (Athens) 2018, 17, (2), 141-151. 7. Shen, F.; Zheng, R. D.; Sun, X. Q.; Ding, W. J.; Wang, X. Y.; Fan, J. G., Gut microbiota dysbiosis in patients with non-alcoholic fatty liver disease. Hepatobiliary Pancreat Dis Int 2017, 16, (4), 375-381. 8. Kral, T. V., Eating Behaviors and Weight Development in Obesity-Prone Children and the Importance of the Research of Albert J. Stunkard. Curr Obes Rep 2016, 5, (1), 126-31. 9. Poret, J. M.; Souza-Smith, F.; Marcell, S. J.; Gaudet, D. A.; Tzeng, T. H.; Braymer, H. D.; Harrison-Bernard, L. M.; Primeaux, S. D., High fat diet consumption differentially affects adipose tissue inflammation and adipocyte size in obesity-prone and obesity-resistant rats. Int J Obes (Lond) 2018, 42, (3), 535-541. 10. Sun, J.; Qiao, Y.; Qi, C.; Jiang, W.; Xiao, H.; Shi, Y. H.; Le, G. W., High-fat-diet induced obesity is associated with decreased antiinflammatory Lactobacillus reuteri sensitive to oxidative stress in mouse Peyer's patches. Nutrition 2016, 32, (2), 265-272. 11. Choi, J. Y.; McGregor, R. A.; Kwon, E. Y.; Kim, Y. J.; Han, Y.; Park, J. H. Y.; Lee, K. W.; Kim, S. J.; Kim, J.; Yun, J. W.; Choi, M. S., The metabolic response to a high-fat diet reveals obesityprone and -resistant phenotypes in mice with distinct mRNA-seq transcriptome profiles. International Journal of Obesity 2016, 40, (9), 1452-1460. 12. Zhang, X. Y.; Shen, W.; Liu, D. Z.; Wang, D. H., Diversity of Thermogenic Capacity Predicts Divergent Obesity Susceptibility in a Wild Rodent. Obesity (Silver Spring) 2018, 26, (1), 111-118. 13. Diane, A.; Vine, D. F.; Russell, J. C.; Heth, C. D.; Pierce, W. D.; Proctor, S. D., Interrelationship of CB1R and OBR pathways in regulation of metabolic, neuroendocrine, and behavioral responses to food restriction and voluntary wheel running. Journal of Applied Physiology 2014, 117, (2), 97104. 14. Li, H.; Xie, Z.; Lin, J.; Song, H.; Wang, Q.; Wang, K.; Su, M.; Qiu, Y.; Zhao, T.; Song, K.; Wang, X.; Zhou, M.; Liu, P.; Zhao, G.; Zhang, Q.; Jia, W., Transcriptomic and metabonomic profiling of obesity-prone and obesity-resistant rats under high fat diet. J Proteome Res 2008, 7, (11), 477583. 15. Derman, R. C.; Ferrario, C. R., Enhanced incentive motivation in obesity-prone rats is mediated by NAc core CP-AMPARs. Neuropharmacology 2018, 131, 326-336. 16. Zhang, P.; Yu, Y.; Qin, Y.; Zhou, Y.; Tang, R.; Wang, Q.; Li, X.; Wang, H.; Weston-Green, K.; Huang, X. F.; Zheng, K., Alterations to the microbiota-colon-brain axis in high-fat-diet-induced obese mice compared to diet-resistant mice. J Nutr Biochem 2018, 65, 54-65. 17. Folch, J.; Lees, M.; Sloane Stanley, G. H., A simple method for the isolation and purification of total lipides from animal tissues. J Biol Chem 1957, 226, (1), 497-509. 18. Gao, X.; Pujos-Guillot, E.; Sebedio, J. L., Development of a quantitative metabolomic approach to study clinical human fecal water metabolome based on trimethylsilylation derivatization 29

ACS Paragon Plus Environment

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

and GC/MS analysis. Anal Chem 2010, 82, (15), 6447-56. 19. McOrist, A. L.; Jackson, M.; Bird, A. R., A comparison of five methods for extraction of bacterial DNA from human faecal samples. J Microbiol Methods 2002, 50, (2), 131-9. 20. Zhang, X.; Zhao, Y.; Zhang, M.; Pang, X.; Xu, J.; Kang, C.; Li, M.; Zhang, C.; Zhang, Z.; Zhang, Y.; Li, X.; Ning, G.; Zhao, L., Structural changes of gut microbiota during berberine-mediated prevention of obesity and insulin resistance in high-fat diet-fed rats. PLoS One 2012, 7, (8), e42529. 21. DeSantis, T. Z.; Hugenholtz, P.; Larsen, N.; Rojas, M.; Brodie, E. L.; Keller, K.; Huber, T.; Dalevi, D.; Hu, P.; Andersen, G. L., Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol 2006, 72, (7), 5069-72. 22. Ludwig, W.; Strunk, O.; Westram, R.; Richter, L.; Meier, H.; Yadhukumar; Buchner, A.; Lai, T.; Steppi, S.; Jobb, G.; Forster, W.; Brettske, I.; Gerber, S.; Ginhart, A. W.; Gross, O.; Grumann, S.; Hermann, S.; Jost, R.; Konig, A.; Liss, T.; Lussmann, R.; May, M.; Nonhoff, B.; Reichel, B.; Strehlow, R.; Stamatakis, A.; Stuckmann, N.; Vilbig, A.; Lenke, M.; Ludwig, T.; Bode, A.; Schleifer, K. H., ARB: a software environment for sequence data. Nucleic Acids Res 2004, 32, (4), 1363-71. 23. Schloss, P. D.; Westcott, S. L.; Ryabin, T.; Hall, J. R.; Hartmann, M.; Hollister, E. B.; Lesniewski, R. A.; Oakley, B. B.; Parks, D. H.; Robinson, C. J.; Sahl, J. W.; Stres, B.; Thallinger, G. G.; Van Horn, D. J.; Weber, C. F., Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol 2009, 75, (23), 7537-41. 24. Cole, J. R.; Wang, Q.; Cardenas, E.; Fish, J.; Chai, B.; Farris, R. J.; Kulam-Syed-Mohideen, A. S.; McGarrell, D. M.; Marsh, T.; Garrity, G. M.; Tiedje, J. M., The Ribosomal Database Project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res 2009, 37, (Database issue), D141-5. 25. Caporaso, J. G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F. D.; Costello, E. K.; Fierer, N.; Pena, A. G.; Goodrich, J. K.; Gordon, J. I.; Huttley, G. A.; Kelley, S. T.; Knights, D.; Koenig, J. E.; Ley, R. E.; Lozupone, C. A.; McDonald, D.; Muegge, B. D.; Pirrung, M.; Reeder, J.; Sevinsky, J. R.; Turnbaugh, P. J.; Walters, W. A.; Widmann, J.; Yatsunenko, T.; Zaneveld, J.; Knight, R., QIIME allows analysis of high-throughput community sequencing data. Nat Methods 2010, 7, (5), 335-6. 26. Hammad, S. S.; Jones, P. J., Dietary Fatty Acid Composition Modulates Obesity and Interacts with Obesity-Related Genes. Lipids 2017, 52, (10), 803-822. 27. Miranda, R. A.; da Silva Franco, C. C.; de Oliveira, J. C.; Barella, L. F.; Tofolo, L. P.; Ribeiro, T. A.; Pavanello, A.; da Conceicao, E. P.; Torrezan, R.; Armitage, J.; Lisboa, P. C.; de Moura, E. G.; de Freitas Mathias, P. C.; Vieira, E., Cross-fostering reduces obesity induced by early exposure to monosodium glutamate in male rats. Endocrine 2017, 55, (1), 101-112. 28. Chakraborty, S. P., Patho-physiological and Toxicological Aspects of Monosodium Glutamate. Toxicol Mech Methods 2018, 1-35. 29. Hernandez, L. F.; Buchwald, P.; Abdulreda, M. H., Effect of Arginase-1 Inhibition on the Incidence of Autoimmune Diabetes in NOD Mice. Curr Res Diabetes Obes J 2018, 5, (3), 1-13. 30. Guevara-Cruz, M.; Vargas-Morales, J. M.; Mendez-Garcia, A. L.; Lopez-Barradas, A. M.; Granados-Portillo, O.; Ordaz-Nava, G.; Rocha-Viggiano, A. K.; Gutierrez-Leyte, C. A.; MedinaCerda, E.; Rosado, J. L.; Morales, J. C.; Torres, N.; Tovar, A. R.; Noriega, L. G., Amino acid profiles of young adults differ by sex, body mass index and insulin resistance. Nutr Metab Cardiovasc Dis 2018, 28, (4), 393-401. 31. Cerdo, T.; Ruiz, A.; Jauregui, R.; Azaryah, H.; Torres-Espinola, F. J.; Garcia-Valdes, L.; Segura, M. T.; Suarez, A.; Campoy, C., Maternal obesity is associated with gut microbial metabolic potential in offspring during infancy. Journal of Physiology and Biochemistry 2018, 74, (1), 159169. 32. Dou, L.; Jourde-Chiche, N.; Faure, V.; Cerini, C.; Berland, Y.; Dignat-George, F.; Brunet, P., The uremic solute indoxyl sulfate induces oxidative stress in endothelial cells. J Thromb Haemost 2007, 5, (6), 1302-8. 33. Sripetchwandee, J.; Chattipakorn, N.; Chattipakorn, S. C., Links Between Obesity-Induced Brain Insulin Resistance, Brain Mitochondrial Dysfunction, and Dementia. Front Endocrinol (Lausanne) 2018, 9, 496. 34. Pereira, R. M.; Botezelli, J. D.; da Cruz Rodrigues, K. C.; Mekary, R. A.; Cintra, D. E.; Pauli, J. R.; da Silva, A. S. R.; Ropelle, E. R.; de Moura, L. P., Fructose Consumption in the Development 30

ACS Paragon Plus Environment

Page 30 of 33

Page 31 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

of Obesity and the Effects of Different Protocols of Physical Exercise on the Hepatic Metabolism. Nutrients 2017, 9, (4). 35. Martins, L. M. S.; Perez, M. M.; Pereira, C. A.; Costa, F. R. C.; Dias, M. S.; Tostes, R. C.; Ramos, S. G.; de Zoete, M. R.; Ryffel, B.; Silva, J. S.; Carlos, D., Interleukin-23 promotes intestinal T helper type17 immunity and ameliorates obesity-associated metabolic syndrome in a murine high-fat diet model. Immunology 2018, 154, (8), 624-636. 36. Kho, Z. Y.; Lal, S. K., The Human Gut Microbiome - A Potential Controller of Wellness and Disease. Front Microbiol 2018, 9, 1835. 37. Nam, M.; Choi, M. S.; Jung, S.; Jung, Y.; Choi, J. Y.; Ryu, D. H.; Hwang, G. S., Lipidomic Profiling of Liver Tissue from Obesity-Prone and Obesity-Resistant Mice Fed a High Fat Diet. Sci Rep 2015, 5, 16984. 38. Schrempft, S.; van Jaarsveld, C. H. M.; Fisher, A.; Herle, M.; Smith, A. D.; Fildes, A.; Llewellyn, C. H., Variation in the Heritability of Child Body Mass Index by Obesogenic Home Environment. JAMA Pediatr 2018, 172, (12), 1153-1160. 39. Lee, J. R.; Tapia, M. A.; Nelson, J. R.; Moore, J. M.; Gereau, G. B.; Childs, T. E.; Vieira-Potter, V. J.; Booth, F. W.; Will, M. J., Sex dependent effects of physical activity on diet preference in rats selectively bred for high or low levels of voluntary wheel running. Behav Brain Res 2018, 359, 95103. 40. Kreutzer, C.; Peters, S.; Schulte, D. M.; Fangmann, D.; Turk, K.; Wolff, S.; van Eimeren, T.; Ahrens, M.; Beckmann, J.; Schafmayer, C.; Becker, T.; Kerby, T.; Rohr, A.; Riedel, C.; Heinsen, F. A.; Degenhardt, F.; Franke, A.; Rosenstiel, P.; Zubek, N.; Henning, C.; Freitag-Wolf, S.; Dempfle, A.; Psilopanagioti, A.; Petrou-Papadaki, H.; Lenk, L.; Jansen, O.; Schreiber, S.; Laudes, M., Hypothalamic Inflammation in Human Obesity Is Mediated by Environmental and Genetic Factors. Diabetes 2017, 66, (9), 2407-2415. 41. Naughton, S. S.; Mathai, M. L.; Hryciw, D. H.; McAinch, A. J., Linoleic acid and the pathogenesis of obesity. Prostaglandins Other Lipid Mediat 2016, 125, 90-9. 42. Mejia-Barradas, C. M.; Del-Rio-Navarro, B. E.; Dominguez-Lopez, A.; Campos-Rodriguez, R.; Martinez-Godinez, M.; Rojas-Hernandez, S.; Lara-Padilla, E.; Abarca-Rojano, E.; Miliar-Garcia, A., The consumption of n-3 polyunsaturated fatty acids differentially modulates gene expression of peroxisome proliferator-activated receptor alpha and gamma and hypoxia-inducible factor 1 alpha in subcutaneous adipose tissue of obese adolescents. Endocrine 2014, 45, (1), 98-105. 43. Muhlhausler, B. S.; Cook-Johnson, R.; James, M.; Miljkovic, D.; Duthoit, E.; Gibson, R., Opposing effects of omega-3 and omega-6 long chain polyunsaturated Fatty acids on the expression of lipogenic genes in omental and retroperitoneal adipose depots in the rat. J Nutr Metab 2010, 2010, 1-9. 44. Moreno-Indias, I.; Cardona, F.; Tinahones, F. J.; Queipo-Ortuno, M. I., Impact of the gut microbiota on the development of obesity and type 2 diabetes mellitus. Front Microbiol 2014, 5, 190. 45. Riva, A.; Borgo, F.; Lassandro, C.; Verduci, E.; Morace, G.; Borghi, E.; Berry, D., Pediatric obesity is associated with an altered gut microbiota and discordant shifts in Firmicutes populations. Environ Microbiol 2017, 19, (1), 95-105. 46. Walters, W. A.; Xu, Z.; Knight, R., Meta-analyses of human gut microbes associated with obesity and IBD. FEBS Lett 2014, 588, (22), 4223-33. 47. Cani, P. D.; Neyrinck, A. M.; Fava, F.; Knauf, C.; Burcelin, R. G.; Tuohy, K. M.; Gibson, G. R.; Delzenne, N. M., Selective increases of bifidobacteria in gut microflora improve high-fat-dietinduced diabetes in mice through a mechanism associated with endotoxaemia. Diabetologia 2007, 50, (11), 2374-83. 48. Duca, F. A.; Sakar, Y.; Lepage, P.; Devime, F.; Langelier, B.; Dore, J.; Covasa, M., Replication of obesity and associated signaling pathways through transfer of microbiota from obese-prone rats. Diabetes 2014, 63, (5), 1624-36. 49. Yin, J.; Li, Y.; Han, H.; Chen, S.; Gao, J.; Liu, G.; Wu, X.; Deng, J.; Yu, Q.; Huang, X.; Fang, R.; Li, T.; Reiter, R. J.; Zhang, D.; Zhu, C.; Zhu, G.; Ren, W.; Yin, Y., Melatonin reprogramming of gut microbiota improves lipid dysmetabolism in high-fat diet-fed mice. J Pineal Res 2018, 65, (4), e12524. 50. Jakobsdottir, G.; Xu, J.; Molin, G.; Ahrne, S.; Nyman, M., High-fat diet reduces the formation of butyrate, but increases succinate, inflammation, liver fat and cholesterol in rats, while dietary 31

ACS Paragon Plus Environment

Journal of Proteome Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

fibre counteracts these effects. PLoS One 2013, 8, (11), e80476. 51. Pataky, Z.; Genton, L.; Spahr, L.; Lazarevic, V.; Terraz, S.; Gaia, N.; Rubbia-Brandt, L.; Golay, A.; Schrenzel, J.; Pichard, C., Impact of Hypocaloric Hyperproteic Diet on Gut Microbiota in Overweight or Obese Patients with Nonalcoholic Fatty Liver Disease: A Pilot Study. Dig Dis Sci 2016, 61, (9), 2721-31. 52. Sergi, D.; Morris, A. C.; Kahn, D. E.; McLean, F. H.; Hay, E. A.; Kubitz, P.; MacKenzie, A.; Martinoli, M. G.; Drew, J. E.; Williams, L. M., Palmitic acid triggers inflammatory responses in N42 cultured hypothalamic cells partially via ceramide synthesis but not via TLR4. Nutr Neurosci 2018, 1-14. 53. Palomer, X.; Pizarro-Delgado, J.; Barroso, E.; Vazquez-Carrera, M., Palmitic and Oleic Acid: The Yin and Yang of Fatty Acids in Type 2 Diabetes Mellitus. Trends Endocrinol Metab 2018, 29, (3), 178-190. 54. Miao, H.; Chen, L.; Hao, L.; Zhang, X.; Chen, Y.; Ruan, Z.; Liang, H., Stearic acid induces proinflammatory cytokine production partly through activation of lactate-HIF1alpha pathway in chondrocytes. Sci Rep 2015, 5, 13092. 55. Shin, J. H.; Sim, M.; Lee, J. Y.; Shin, D. M., Lifestyle and geographic insights into the distinct gut microbiota in elderly women from two different geographic locations. J Physiol Anthropol 2016, 35, (1), 31. 56. Garcia-Mantrana, I.; Selma-Royo, M.; Alcantara, C.; Collado, M. C., Shifts on Gut Microbiota Associated to Mediterranean Diet Adherence and Specific Dietary Intakes on General Adult Population. Front Microbiol 2018, 9, 890. 57. Gonzalez-Sarrias, A.; Romo-Vaquero, M.; Garcia-Villalba, R.; Cortes-Martin, A.; Selma, M. V.; Espin, J. C., The Endotoxemia Marker Lipopolysaccharide-Binding Protein is Reduced in Overweight-Obese Subjects Consuming Pomegranate Extract by Modulating the Gut Microbiota: A Randomized Clinical Trial. Molecular Nutrition & Food Research 2018, 62, (11), 1-32. 58. Yue, J. T.; Abraham, M. A.; Bauer, P. V.; LaPierre, M. P.; Wang, P.; Duca, F. A.; Filippi, B. M.; Chan, O.; Lam, T. K., Inhibition of glycine transporter-1 in the dorsal vagal complex improves metabolic homeostasis in diabetes and obesity. Nat Commun 2016, 7, 13501. 59. Li, R.; Wang, H.; Shi, Q.; Wang, N.; Zhang, Z.; Xiong, C.; Liu, J.; Chen, Y.; Jiang, L.; Jiang, Q., Effects of oral florfenicol and azithromycin on gut microbiota and adipogenesis in mice. PLoS One 2017, 12, (7), e0181690. 60. Liu, R. X.; Hong, J.; Xu, X. Q.; Feng, Q.; Zhang, D. Y.; Gu, Y. Y.; Shi, J.; Zhao, S. Q.; Liu, W.; Wang, X. K.; Xia, H. H.; Liu, Z. P.; Cui, B.; Liang, P. W.; Xi, L. Q.; Jin, J. B.; Ying, X. Y.; Wang, X. L.; Zhao, X. J.; Li, W. Y.; Jia, H. J.; Lan, Z.; Li, F. Y.; Wang, R.; Sun, Y. K.; Yang, M. L.; Shen, Y. X.; Jie, Z. Y.; Li, J. H.; Chen, X. M.; Zhong, H. Z.; Xie, H. L.; Zhang, Y. F.; Gu, W. Q.; Deng, X. X.; Shen, B. Y.; Xu, X.; Yang, H. M.; Xu, G. W.; Bi, Y. F.; Lai, S. H.; Wang, J.; Qi, L.; Madsen, L.; Wang, J. Q.; Ning, G.; Kristiansen, K.; Wang, W. Q., Gut microbiome and serum metabolome alterations in obesity and after weight-loss intervention. Nature Medicine 2017, 23, (7), 859-868.

32

ACS Paragon Plus Environment

Page 32 of 33

Page 33 of 33 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Proteome Research

For TOC Only

33

ACS Paragon Plus Environment