- Open Access
Heat-treated high-fat diet modifies gut microbiota and metabolic markers in apoe−/− mice
Nutrition & Metabolism volume 13, Article number: 22 (2016)
High-fat diet has been known to have adverse effects on metabolic markers, as well as the gut microbiota. However, the effect of heat processing of high-fat diet, which leads to formations of advanced glycation end products (AGEs) has not been clearly distinguished from the effect of unheated fat. This study compared the effect of high-fat diet with heat-treated high-fat diet on adiposity, atherosclerosis and gut microbiota composition in the caecum of apoe −/− mice.
Male apoe −/− mice were fed either low-fat (LF) control diet, high-fat (40 E% saturated fat, HF) control diet, or heat-treated high-fat (200 °C for 10 min, HT) diet, for 8 weeks. The plasma samples were used in the analysis of Nε-carboxy-methyl-lysine (CML) and Nε-carboxy-ethyl-lysine (CEL). The heart samples were analysed for atherosclerotic plaques, and the DNA from caecum was extracted and analysed for microbiota composition using 16S rRNA gene sequencing on a Miseq instrument. Additionally, the functions of microbial communities were also predicted based on the bacterial 16S rRNA gene sequence using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt).
Here we found that HT modifies gut microbiota composition and host adiposity. Prediction of bacterial gene functions based on 16S rRNA gene sequence revealed that HF increased bacterial genera enriched in lipid metabolism genes, while HT did not. Plasma CML and CEL increased 1.7 and 2.5 times, respectively, in mice fed HT as compared to mice fed HF. Despite lower adiposity, mice fed HT maintained atherosclerosis and displayed enlarged spleens.
The results suggested that heat processing of high-fat diet modifies the substrates reaching the lower gut of apoe −/− mice, resulting in different effects on gut microbiota composition. AGEs seem to maintain the effect on atherosclerosis, despite lower adiposity, and causing enlarged spleens, which possibly reflect elevated levels of inflammation in the body.
Modern diet involves heat processing of food, which leads to formation of Maillard reaction products (MRPs). These products are responsible for the aroma, colour, flavour, and texture as well as deterioration of nutritional value in the processed food and adverse health implications [1, 2]. Recently, more attention has been paid to the health implication of Maillard reaction products, especially those that lead to irreversible modifications in amino acids, peptides and proteins; advanced glycation end products (AGEs) [3–6]. AGEs are suggested to be associated to risk markers for inflammation mediated pathological conditions; such as cardiovascular disease, Alzheimer’s disease and diabetes-linked-complications [5, 7–12]. The mechanism for AGEs induced toxicity can be mediated through binding to AGE receptors (RAGE) and inducing pro-inflammatory cascade [8, 13, 14], and consequently associated pathological conditions. Cross link formation affecting the structure of the protein may also affect the digestibility and thus modified proteins may reach the lower gut where they can be fermented by gut microbiota .
The gut microbiota is intimately involved in numerous aspects of host physiology. The overall balance in the gut microbiota composition is an important key factor ensuring normal host functions . Several studies have revealed that presence or absence of gut microbiota or some specific groups of bacteria have been found to contribute to development of many diseases such as type-2 diabetes [17, 18], artherosclerosis [19, 20] and systemic inflammatory response syndrome . The composition of the gut microbiota is known to be affected by dietary components reaching the lower gut. So far, the role of process induced structure alteration in proteins on the digestibility and possible metabolism by gut microbiota is an area that is largely unexplored.
Atherosclerosis is one form of chronic inflammatory response to accumulation of fat in blood vessels. The apolipoprotein E knock out (apoe −/−) mouse is an animal model used for assessing effects of nutritional agents on lipid metabolism, inflammation and atherosclerosis . Apoe−/− mice have an impaired ability for clearance of plasma lipoprotein, leading to the development of atherosclerosis in a short time. Several studies have been focusing on the pathological effects of high-fat diet. High intake of dietary fat has been known to induce several dynamic metabolic alterations, especially atherosclerosis, as well as changes in gut microbiota composition . However, the fact that intake of fat is usually accompanied by heat processing , the specific effect of high fat intake versus the effect of AGEs has not been clearly distinguished. In the present study, we aimed to compare the effect of high-fat diet (40 E% saturated fat) with heat-treated (200 °C for 10 min) high-fat diet on adiposity, atherosclerosis and gut microbiota composition in the caecum of apoe −/− mice. Moreover, the study compared the possible effect of intake of high-fat diet on the above mentioned end points by comparing with intake of low-fat diet. Additionally, the functions of microbial communities were also predicted based on the bacterial 16S rRNA gene sequence using a recently developed software, Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) .
Male apoe −/− mice (Scanbur AB, Karlslunde, Denmark), 6 weeks of age, were adapted to the environment at the animal facility for 2 weeks before starting the experiment. At the age of 8 weeks, the mice were randomly divided into three weight-matched groups (n = 10, five mice/cage). The mice were fed low-fat (LF) control diet, high-fat (HF) control diet, or heat-treated high-fat (HT) diet. A picture of HF and HT is available as Additional file 1: Figure S1. The mice were weighed every week. After 8 weeks, the mice were anesthetized with isofluorane (Abbott Scandinavia AB, Solna, Sweden) and terminated by heart puncture. The liver, spleen and epididymal fat pads were weighed. Blood plasma, heart, and caecum were collected and frozen at −40 °C until further analyses. The plasma samples were used in the analysis of CML and CEL. The heart samples were analysed for atherosclerotic plaques, and the DNA from caecum was extracted and analysed for microbiota composition using 16S rRNA gene sequencing on a Miseq instrument.
High performance liquid chromatography coupled to tandem mass spectrometry for determination of CML and CEL
Plasma samples (50 μL) were hydrolysed for 12 h at 110 °C using 6 M HCl, together with isotope labelled d4-CML and d4-CEL as internal standard (Larodan Fine Chemicals AB, Malmö Sweden), and analysed using high-pressure liquid chromatography mass spectrometry (HPLC-MS/MS). The chromatographic separation of CML and CEL in the hydrolysed samples was performed using Accela UHPLC pump with auto-injector. Detection was performed using a LTQ VelosPro Orbitrap massspectrometer (Thermo Scientific, Waltham, USA) run in positive electrospray ionization ion trap tandem mass spectrometry (MS/MS) mode, detecting two selective-reaction monitoring (SRM) transitions for each analyte, and internal standard. The Xcalibur software (Thermo Scientific) was used both for data acquisition and evaluation. Solid phase extraction, chromatographic parameters, ion source parameter and the SRM transitions are the same as described in Tareke E. et al. .
Quantification of atherosclerosis
Frozen 10 μm sections of the aortic root region of the heart were prepared using a cryostat (Leica CM 1950, Leica Biosystems, Nusslock, Germany) and stained with Oil Red O (Histolab, Gothenburg, Sweden) and haematoxylin (Mayer’s HTX, Histolab, Gothenburg, Sweden). Amount of atherosclerotic plaques covering the vessel area was quantified by a blinded observer using BioPixQ 2.0 (Biopix software, Gothenburg, Sweden). Three sections/mouse, with well-oriented root regions and all three valve cusps present, were evaluated.
Caecal microbiota composition
A total of 22 caecal samples (LF (n = 4), HF (n = 8), HT (n = 9)) were used in the microbiota analysis. Caecal tissue and content were thawed on ice and DNA was extracted using the QIAamp DNA Stool Mini Kit (Qiagen), with an addition of a bead beating step. Sterile glass beads (1 mm) were added in combination with stool lysis buffer and cell disruption was performed for 2 × 2 min at 25 Hz using a TissueLyser (Qiagen), followed by a heating step at 95 °C for 5 min. After lysis, DNA-damaging substances and PCR inhibitors were removed using InhibitEX tablet (provided with the kit) and the DNA was purified on QIAamp Mini spin columns. A normalized input of 5 ng/μL of the DNA was used in PCR reactions where 16S rRNA genes were amplified prior to sequencing.
The V1-V3 regions of 16S rRNA genes were amplified using a limited cycle PCR with forward and reverse primers containing Illumina adapter sequences and dual-index barcodes used for tagging each sample, primer sequences are listed in Table 1. Paired-end sequencing with a read length of 2×300 bp was performed on a Miseq Instrument (Illumina, Inc, San Diego, California) using a Miseq v3 reagent kit (Illumina, Inc, San Diego, California). Sequences were analysed using Quantitative Insights into Microbial Ecology (QIIME), as previously described . The Operational Taxonomic Units (OTUs) were picked using a closed-reference OTU picking with UCLUST  algorithm against the Greengenes  database (version 13_8), pre-clustered at 97 % identify. After quality filtering, a total of 2,280,567 reads were included for downstream analyses and an average of 103,662 sequences (SD: 18,962) were assigned to each sample (ranging from 67,817 to 145,015 sequences). The relative abundances of bacterial taxa at phylum and genus levels were thereafter compared between groups. To correct for sampling depth differences, 67,817 reads/sample were randomly selected and utilized for further calculation of alpha-diversity and weighted and unweighted Unifrac as well as correlations between the gut microbiota and different metabolic markers.
Prediction of bacterial metagenomes using PICRUSt
The bacterial metagenomes were reconstructed based on the bacterial 16S rRNA gene sequence using the open-source software, PICRUSt . The OTU table generated in QIIME from the 16S sequencing data at 67,817 rarefied sequences per sample was used as an input. The copy number per OTU was normalized before the metagenome was predicted using Kyoto Encyclopedia of Genes and Genomes (KEGG) database . The output from metagenome prediction was an annotated table of predicted gene family counts for each sample, where gene families were grouped by KEGG Orthology (KO) identifiers.
One-way analysis of variance (ANOVA) was used to calculate the differences in metabolic markers in different groups of mice using Graphpad Prism 6. Differences in within-community richness (α-diversity) were calculated in QIIME using a non-parametric t-test and the P-values were corrected for multiple comparisons using False Discovery Rate (FDR) correction . Differences in community composition among groups of samples (β-diversity) were analyzed using the non-parametric analysis of similarity (ANOSIM)  statistical test in QIIME on both unweighted- and weighted Unifrac phylogenetic metrics. Graphpad Prism 6 was also used to identify significant differences in taxonomic distributions at phylum and genus levels between the groups of mice using two-way ANOVA and the Holm-Sidak method of correction for multiple comparisons. All the differences were considered to be significant at P < 0.05. Moreover, linear discriminant analysis (LDA) effect size (LEfSe) analysis  was performed to identify differentially abundant bacterial taxa (from phylum to species level) as well as the predicted metagenomes in each group of mice. A strict version of LEfSe was used, the discovered bacterial taxa and predicted genes with LDA score higher than two were considered to be enriched in the respective group as compared to all other groups. Partial Least Squares (PLS) Regression analysis was performed using SIMCA14 (Umetrics, Umeå, Sweden) and PLS loading score scatter plots were plotted to illustrate the associations between the parameters. Pearson’s correlation was calculated for each pairwise combination of bacterial genera versus relative spleen weight, plasma CEL and CML, and atherosclerosis plaque using Graphpad Prism 6 and the P-values were then corrected by Benjamini-Hochberg procedure for multiple comparisons [30, 33].
Considering no difference in food consumption of the mice in different groups (Fig. 1a), mice fed HF had significantly higher body weight at the end point as compared to mice fed LF (P < 0.001, Fig. 1b). The liver weight did not differ between groups (Fig. 1c), while the epididymal fat pad weight was significantly lower in the HT and LF groups compared to the HF group (Fig. 1d). Relative spleen weight was significantly higher in HT mice as compared to HF mice (P < 0.05, Fig. 1e). Atherosclerotic plaque size (% of vessel area) was higher in HF and HT as compared to LF (both P < 0.05, Fig. 1f). Analysis of Advanced glycation end products (CML and CEL levels) in plasma of mice revealed significant increases in levels of CML and CEL in plasma of HT mice as compared to LF and HF (all P < 0.05, Fig. 1g). CML increased 1.7 times in HT mice as compared to LF and HF, while CEL increased 3.5 and 2.5 times in HT as compared to LF and HF, respectively. There was no significant difference in plasma CML and CEL levels between HF and LF mice.
Microbial diversity indices
Estimation of α-diversity of bacterial 16S rRNA gene at the sequence number of 67,817 sequences/sample (Fig. 2a) demonstrated that HF had a tendency for decreased α-diversity as compared to LF. Moreover, HT lowered the α-diversity even more as the observed OTUs and phylogenetic diversity whole tree (PD whole-tree) indices showed significant differences between HT and LF (P < 0.05), although there was no significant difference with Chao1 index. ANOSIM analyses of both unweighted and weighted Unifrac revealed significant differences between the communities of the mice in the three treatment groups (P < 0.01 and P < 0.05, respectively). Principle coordinate analysis (PCoA) plots for unweighted and weighted Unifrac are shown in Fig. 2b and c, respectively.
Phylum-level taxonomic distributions
Relative abundance at phylum-level (Fig. 2d) revealed that Firmicutes was found to be the most dominant phylum in all groups and was followed by Bacteroidetes. The level of Firmicutes was significantly higher in HT than LF and HF (both P < 0.01), whereas the level of Bacteroidetes was significantly lower in HT than LF and HF (P < 0.001 and P < 0.01, respectively).
Genus-level taxonomic distributions
The most abundant genera presented in the mice (more than 5 % relative abundance) are shown in Fig. 2e. As compared to LF, HF significantly increased Mucispirillum (P < 0.05) and unclassified family of Clostridiales (P < 0.05) and reduced unclassified genus of S24-7 (P < 0.001), while similar effects were observed in HT with significant increase in the relative abundance of unclassified family of Clostridiales (P < 0.001) and significant decrease in the relative abundance of unclassified genus of S24-7 (P < 0.01). Moreover, a significant decrease of an unclassified genus of Rikenellaceae (P < 0.001) was also observed in HT, as compared to LF. In addition, HT was found to reduce the relative abundance of an unclassified genus of Rikenellaceae (P < 0.001) as well as to increase the relative abundance of Allobaculum (P < 0.001) and unclassified family of Clostridiales (P < 0.001) as compared to HF.
Differences in bacterial composition between groups of treatment using LEfSe
The strict version of LEfSe (all against all) detected 24 bacterial clades showing statistically significant differences between different groups of treatment with LDA scores higher than two (Fig. 3a). A Cladogram of these bacterial taxa is shown in Additional file 2: Figure S2. The most enriched bacteria in LF were Sutterella (from phylum to genus), Anaeroplasma (from phylum to genus), Adlercreutzia and Lactobacillus. Mucispirillum schaedleri (from phylum to species) and Lactococcus (from family to genus) were found to be the most enriched bacteria in HF while Allobaculum (from order to genus) was found to be enriched in HT.
Prediction of functional structure of genes based on 16S gene information using PICRUSt
By applying PICRUSt to the 16S sequencing data, a total of 328 genes were predicted to be involved in different functions (data not shown). Assigning the strict version (all against all) of LEfSe to the predicted genes from PICRUSt, detected a total of 11, 9 and 5 genes enriched in LF, HF and HT, respectively, with LDA scores higher than two (Fig. 3b). Of these, 7, 8 and 4 genes, enriched in LF, HF and HT, respectively, were involved in metabolism. Predicted metabolism genes enriched in LF encoded for metabolism of two amino acid groups; arginine and proline (KO00330), glycine, serine and threonine (KO00260), energy metabolism, particularly for methane metabolism (KO00680), enzyme families, especially for peptidases, and cofactors and vitamins, especially for vitamin B6 (KO00750). Predicted genes encoding for lipid metabolism, especially for lipid biosynthesis proteins (KO01004) and glycerophospholipid (KO00564) were enriched in HF. Moreover, metabolism of pyruvate (KO00620), lysine biosynthesis (KO00300), beta-alanine (KO00410), as well as pantothenate and CoA biosynthesis (KO00770) were also enriched in HF. These genes enriched in HF were shown to be involved in a cascade of the fatty acid biosynthesis KEGG pathway. Predicted genes encoding for carbohydrate metabolism; fructose and mannose metabolism (KO00051) and pentose and glucuronate interconversions (KO00040), energy metabolism encoding for photosynthesis proteins and an unclassified metabolism (no KO number) were enriched in HT.
Additionally, genes enriched in LF were involved in genetic information, especially for replication and repair of chromosome, ribosome biogenesis and an unclassified translation protein. A gene enriched in HF was involved in cellular processes, especially for cell growth and death and a gene enriched in HT was involved in an unclassified function (Fig. 3b).
Gut microbiota and metabolic biomarkers
Loading and score scatter (Fig. 4) PLS plots illustrated associations between the gut microbiota and different biomarkers as well as revealed separate clusters of the mice by different treatments (Fig. 4, small panel). The results supported the LEfSe results (Fig. 3a and Additional file 2: Figure S2) where those bacteria enriched in different groups of treatment were located in the corresponding area of the PLS plot (Fig. 4). Sutterella, Anaeroplasma, Adlercreutzia and Lactobacillus which were found to be enriched in LF were located in the lower component of the PLS plot. Mucispirillum and Lactococcus enriched in HF were located in the upper right component, while Allobaculum was located in the upper left component of the PLS plot. Correlations between the gut microbiota and different biomarkers were also analyzed. Relative spleen weight was found to be positively correlated with Dehalobacterium (P < 0.01), Clostridium (P < 0.05), Prevotella (P < 0.05), unclassified family of Clostridiales (P < 0.05) and negatively correlated with unclassified genus of Rikenellaceae (P < 0.05) (Figs. 4 and 5). Plasma CEL and CML levels were located in the right component of and close to the relative spleen weight on the PLS plot. These suggested positive correlations between plasma CML and CEL versus relative spleen weight, and thus versus those bacterial genera that correlated with spleen weight, although the correlations were not statistically significant after correction for multiple corrections. No correlation was found between the gut microbiota and amount of atherosclerotic plaques.
The findings from this study revealed that heat processing of fat led to changes in its effects on metabolic markers and the gut microbiota in apoe −/− mice as compared to unheated fat. Mice fed HT displayed elevated circulating levels of CML and CEL as well as enlarged spleens, to a higher extent than mice fed HF. Enlargement of the spleen has previously been observed even in healthy mice with intra-peritoneal administration of RAGE and this has been suggested to be exerted through stimulation of adaptive immune system in response to inflammation . Interestingly, the amounts of atherosclerotic plaque were similar in HF and HT, which were both higher than LF. The amount of atherosclerotic plaques was anticipated to be higher in mice fed HT as heat processing of fat would add up the effect of AGEs to the effect of fat on atherosclerosis. The lower epididymal fat-pad weight observed in HT mice, to the same level as in LF mice, could possibly provide an explanation. Chemical modifications of fat by heat-treatment might abate the effect of unheated fat, leaving only the effect of AGEs on atherosclerosis. Thus, HT mice reached similar levels of atherosclerosis as HF mice despite lower adiposity, which may suggest an additive effect of AGEs on atherosclerosis in the HT mice. Heat-treated LF has not been included in this study as the aim was only to compare the effects of fat and heat-treatment of HF (HT), using LF as a baseline control. However, heat-treatment of LF diet may provide more insight into the role of fat and AGEs on atherogenesis.
The results from 16S sequencing of caecal samples showed that HF and HT had an impact on the composition of the gut microbiota in the mice. Both HF and HT lowered the diversity of the gut microbiota, although the change was only significant with HT. The most dominant phyla in all groups were Firmicutes and Bacteroidetes, whose relative abundance were significantly altered with HT, as well as with HF but to a lesser extent. An imbalanced ratio of Firmicutes and Bacteroidetes has been found to be associated with risk factors of obesity in both animals and humans. High dimensional comparison using LEfSe analysis showed that both HF and HT mice differed in their gut microbiota composition and functions, as compared to each other and to LF mice. The LF group had a high diversity of bacterial genera from different phyla. Although most groups of bacteria found enriched in LF are not well known regarding their roles for the host, the predicted gene functions revealed increases of genes involved in a broad range of metabolism pathways including amino acids, enzymes, as well as cofactors and vitamins. Mucispirillum found enriched in the mice fed HF, is one of the putative mucin degraders  and has been found in a study of Belzer et al., to be involved in the onset of symptomatic colitis in mice . Interestingly, the lactic acid bacteria, Lactococcus, found enriched also in the mice fed HF, has been found in a study of Parks et al., to have positive correlation with body fat percentage gain in mice fed high-fat, high-sucrose diet . Although genes enriched in HF mice were encoded for several metabolism pathways, it is remarkable that most of these were involved in a cascade of a fatty acid synthesis pathway. Some bacteria have been known to synthesize glycerophospholipid on the inner leaflet of the inner membrane and different lipid biosynthetic mechanisms have previously been defined . These bacterial genera could therefore generate even more fat to the host’s energy metabolism. In addition, the gut microbiota composition has also been shown to play a role in how much energy is harvested from the diet [39–41], and indeed, mice fed HF had significantly higher epididymal fat pad weight, as well as slightly higher body weight than mice fed HT. Interestingly, Allobaculum which has previously been found to be decreased in mice fed high-fat diet [37, 42], became enriched in HT mice in this study. Although these studies showed somewhat similar 16S gene results as observed in this current study, it is worth mentioning that some of them were performed on animal faeces as opposed to caecum and these two locations in the gut may harbour different bacterial species. Together with the finding that no bacterial genes encoding for lipid metabolism were enriched in mice fed HT, we hypothesize that heat processing of fat changed the chemical properties of the fat, allowing different substrates to reach the lower gut. However, the exact mechanism of how heat-treated fat appears to be less accessible to the gut microbiota than unheated fat needs further investigation.
In conclusion, intake of high-fat diet led to increased adiposity, decreased microbiota diversity and altered gut microbiota composition. Chemical modifications of high-fat diet by heat processing seemed to reduce the substrates available for members of the gut microbiota having lipid metabolism genes, resulting in lower adiposity. However, the overall alterations of the gut microbiota are still towards those associated with risk factors of obesity. Moreover, AGEs formed in HT could be a possible cause of spleen enlargement, which possibly reflect elevated levels of inflammation in the body. Lower intake of heat-processed food may need to be taken into consideration when developing a dietary preventative approach of the metabolic syndrome.
The study was approved by the local ethical review committee for animal experiments in Lund, Sweden (approval number M-295-12).
Advanced Glycation End products
analysis of similarity
analysis of variance
apolipoprotein E knock out
high-pressure liquid chromatography mass spectrometry
Kyoto Encyclopedia of Genes and Genomes
Kyoto Encyclopedia of Genes and Genomes Orthology
linear discriminant analysis
linear discriminant analysis effect size
Maillard reaction products
tandem mass spectrometry
operational taxonomic unit
principle coordinate analysis
polymerase chain reaction
Phylogenetic Investigation of Communities by Reconstruction of Unobserved States
Partial Least Squares
Quantitative Insights into Microbial Ecology
ribosomal ribonucleic acid
Finot PA. Historical perspective of the Maillard reaction in food science. Ann N Y Acad Sci. 2005;1043:1–8.
Newton AE, Fairbanks AJ, Golding M, Andrewes P, Gerrard JA. The role of the Maillard reaction in the formation of flavour compounds in dairy products--not only a deleterious reaction but also a rich source of flavour compounds. Food Funct. 2012;3(12):1231–41.
Chuyen NV. Toxicity of the AGEs generated from the Maillard reaction: on the relationship of food-AGEs and biological-AGEs. Mol Nutr Food Res. 2006;50(12):1140–9.
Nemet I, Varga-Defterdarovic L, Turk Z. Methylglyoxal in food and living organisms. Mol Nutr Food Res. 2006;50(12):1105–17.
Vistoli G, De Maddis D, Cipak A, Zarkovic N, Carini M, Aldini G. Advanced glycoxidation and lipoxidation end products (AGEs and ALEs): an overview of their mechanisms of formation. Free Radic Res. 2013;47 Suppl 1:3–27.
Vlassara H, Striker GE. AGE restriction in diabetes mellitus: a paradigm shift. Nat Rev Endocrinol. 2011;7(9):526–39.
Alexiou P, Chatzopoulou M, Pegklidou K, Demopoulos VJ. RAGE: a multi-ligand receptor unveiling novel insights in health and disease. Curr Med Chem. 2010;17(21):2232–52.
Chuah YK, Basir R, Talib H, Tie TH, Nordin N. Receptor for advanced glycation end products and its involvement in inflammatory diseases. Int J Inflam. 2013;2013:403460.
Hegab Z, Gibbons S, Neyses L, Mamas MA. Role of advanced glycation end products in cardiovascular disease. World J Cardiol. 2012;4(4):90–102.
Pasinetti GM, Zhao Z, Qin W, Ho L, Shrishailam Y, Macgrogan D, Ressmann W, Humala N, Liu X, Romero C, et al. Caloric intake and Alzheimer's disease. Experimental approaches and therapeutic implications. Interdiscip Top Gerontol. 2007;35:159–75.
Baynes JW, Thorpe SR. Glycoxidation and lipoxidation in atherogenesis. Free Radic Biol Med. 2000;28(12):1708–16.
Miyazawa T, Nakagawa K, Shimasaki S, Nagai R. Lipid glycation and protein glycation in diabetes and atherosclerosis. Amino Acids. 2012;42(4):1163–70.
Poulsen MW, Hedegaard RV, Andersen JM, de Courten B, Bugel S, Nielsen J, Skibsted LH, Dragsted LO. Advanced glycation endproducts in food and their effects on health. Food Chem Toxicol. 2013;60:10–37.
Rojas A, Delgado-Lopez F, Gonzalez I, Perez-Castro R, Romero J, Rojas I. The receptor for advanced glycation end-products: a complex signaling scenario for a promiscuous receptor. Cell Signal. 2013;25(3):609–14.
Mills DJS, Tuohy KM, Booth J, Buck M, Crabbe MJC, Gibson GR, Ames JM. Dietary glycated protein modulates the colonic microbiota towards a more detrimental composition in ulcerative colitis patients and non-ulcerative colitis subjects. J Appl Microbiol. 2008;105(3):706–14.
Sekirov I, Russell SL, Antunes LC, Finlay BB. Gut microbiota in health and disease. Physiol Rev. 2010;90(3):859–904.
Cani PD, Bibiloni R, Knauf C, Waget A, Neyrinck AM, Delzenne NM, Burcelin R. Changes in gut microbiota control metabolic endotoxemia-induced inflammation in high-fat diet-induced obesity and diabetes in mice. Diabetes. 2008;57(6):1470–81.
Cani PD, Neyrinck AM, Fava F, Knauf C, Burcelin RG, Tuohy KM, Gibson GR, Delzenne NM. Selective increases of bifidobacteria in gut microflora improve high-fat-diet-induced diabetes in mice through a mechanism associated with endotoxaemia. Diabetologia. 2007;50(11):2374–83.
Bjorkbacka H, Kunjathoor VV, Moore KJ, Koehn S, Ordija CM, Lee MA, Means T, Halmen K, Luster AD, Golenbock DT, et al. Reduced atherosclerosis in MyD88-null mice links elevated serum cholesterol levels to activation of innate immunity signaling pathways. Nat Med. 2004;10(4):416–21.
Bukowska H, Pieczul-Mroz J, Jastrzebska M, Chelstowski K, Naruszewicz M. Decrease in fibrinogen and LDL-cholesterol levels upon supplementation of diet with Lactobacillus plantarum in subjects with moderately elevated cholesterol. Atherosclerosis. 1998;137(2):437–8.
Shimizu K, Ogura H, Goto M, Asahara T, Nomoto K, Morotomi M, Yoshiya K, Matsushima A, Sumi Y, Kuwagata Y, et al. Altered gut flora and environment in patients with severe SIRS. J Trauma. 2006;60(1):126–33.
Kolovou G, Anagnostopoulou K, Mikhailidis DP, Cokkinos DV. Apolipoprotein E knockout models. Curr Pharm Des. 2008;14(4):338–51.
Manzel A, Muller DN, Hafler DA, Erdman SE, Linker RA, Kleinewietfeld M. Role of "Western diet" in inflammatory autoimmune diseases. Curr Allergy Asthma Rep. 2014;14(1):404.
Langille MGI, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Vega Thurber RL, Knight R, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotech. 2013;31(9):814–21.
Tareke E, Forslund A, Lindh CH, Fahlgren C, Östman E. Isotope dilution ESI-LC-MS/MS for quantification of free and total Nε-(1-Carboxymethyl)-l-Lysine and free Nε-(1-Carboxyethyl)-l-Lysine: Comparison of total Nε-(1-Carboxymethyl)-l-Lysine levels measured with new method to ELISA assay in gruel samples. Food Chem. 2013;141(4):4253–9.
Caporaso J, Kuczynski J, Stombaugh J, Bittinger K, Bushman F, Costello E, Fierer N, Peña A, Goodrich J, Gordon J, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1.
DeSantis T, Hugenholtz P, Larsen N, Rojas M, Brodie E, Keller K, Huber T, Dalevi D, Hu P, Andersen G. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–72.
Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40(Database issue):D109–14.
BaY HY. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300.
Clarke KR. Non-parametric multivariate analyses of changes in community structure. Aust J Ecol. 1993;18(1):117–43.
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60.
McDonald JH. Multiple comparisons: Controlling the false discovery rate: Benjamini–Hochberg procedure. In: Handbook of Biological Statistics. Baltimore, Maryland: Sparkly House Publishing; 2014. p. 254–60.
Brisslert M, Amu S, Pullerits R. Intra-peritoneal sRAGE treatment induces alterations in cellular distribution of CD19+, CD3+ and Mac-1+ cells in lymphoid organs and peritoneal cavity. Cell Tissue Res. 2013;351(1):139–48.
Robertson BR, O'Rourke JL, Neilan BA, Vandamme P, On SL, Fox JG, Lee A. Mucispirillum schaedleri gen. nov., sp. nov., a spiral-shaped bacterium colonizing the mucus layer of the gastrointestinal tract of laboratory rodents. Int J Syst Evol Microbiol. 2005;55(Pt 3):1199–204.
Belzer C, Gerber GK, Roeselers G, Delaney M, DuBois A, Liu Q, Belavusava V, Yeliseyev V, Houseman A, Onderdonk A, et al. Dynamics of the Microbiota in Response to Host Infection. PLoS One. 2014;9(7), e95534.
Parks Brian W, Nam E, Org E, Kostem E, Norheim F, Hui Simon T, Pan C, Civelek M, Rau Christoph D, Bennett Brian J, et al. Genetic Control of Obesity and Gut Microbiota Composition in Response to High-Fat, High-Sucrose Diet in Mice. Cell Metab. 2013;17(1):141–52.
Parsons JB, Rock CO. Bacterial lipids: Metabolism and membrane homeostasis. Prog Lipid Res. 2013;52(3):249–76.
Velagapudi VR, Hezaveh R, Reigstad CS, Gopalacharyulu P, Yetukuri L, Islam S, Felin J, Perkins R, Borén J, Orešič M, et al. The gut microbiota modulates host energy and lipid metabolism in mice. J Lipid Res. 2010;51(5):1101–12.
Bäckhed F, Ding H, Wang T, Hooper LV, Koh GY, Nagy A, Semenkovich CF, Gordon JI. The gut microbiota as an environmental factor that regulates fat storage. Proc Natl Acad Sci U S A. 2004;101(44):15718–23.
Ley RE, Bäckhed F, Turnbaugh P, Lozupone CA, Knight RD, Gordon JI. Obesity alters gut microbial ecology. Proc Natl Acad Sci U S A. 2005;102(31):11070–5.
Ravussin Y, Koren O, Spor A, LeDuc C, Gutman R, Stombaugh J, Knight R, Ley RE, Leibel RL. Responses of Gut Microbiota to Diet Composition and Weight Loss in Lean and Obese Mice. Obesity (Silver Spring, Md). 2012;20(4). doi:10.1038/oby.2011.1111.
This research was funded by the Antidiabetic Food Center (AFC), a VINN Excellence Center funded by VINNOVA at Lund University, the Direktör Albert Påhlsson Foundation and the Per Håkansson Foundation. Chrysoula Matziouridou is greatly acknowledged for performing the atherosclerosis assessment.
The authors declare that they have no competing interests.
ET designed the study, FF performed the animal study, NM performed gut microbiota analyses and ET performed AGEs analysis. All authors took part in writing the paper. All authors read and approved the final manuscript.
Unheated (left) and heat-treated (200 °C for 10 min, right) high-fat diets. (PDF 241 kb)
Cladogram plot showing bacterial taxa with LDA scores higher than two using strict version of LEfSe. Different colors represent the most abundant taxa in different groups of mice (yellow indicating LF, red indicating HF and blue indicating HT). Circles represent phylogenetic levels from phylum to genus. The sizes of the circle are proportional to the taxon’s abundance. (PDF 436 kb)