Original Research ARTICLE
Interactions between Obesity Status and Dietary Intake of Monounsaturated and Polyunsaturated Oils on Human Gut Microbiome Profiles in the Canola Oil Multicenter Intervention Trial (COMIT)
- 1Department of Human Nutritional Sciences, University of Manitoba, Winnipeg, MB, Canada
- 2Richardson Centre for Functional Foods and Nutraceuticals, University of Manitoba, Winnipeg, MB, Canada
- 3Gut Microbiome Laboratory, Department of Animal Science, University of Manitoba, Winnipeg, MB, Canada
- 4Department of Medical Microbiology and Infectious Diseases, University of Manitoba, Winnipeg, MB, Canada
Long-term dietary fatty acid intake is believed to induce changes in the human gut microbiome which might be associated with human health or obesity status; however, considerable debate remains regarding the most favorable ratios of fatty acids to optimize these processes. The objective of this sub-study of a double-blinded randomized crossover clinical study, the canola oil multi-center intervention trial, was to investigate effects of five different novel oil blends fed for 30 days each on the intestinal microbiota in 25 volunteers with risk of metabolic syndrome. The 60 g treatments included three MUFA-rich diets: (1) conventional canola oil (Canola); (2) DHA-enriched high oleic canola oil (CanolaDHA); (3) high oleic canola oil (CanolaOleic); and two PUFA-rich diets: (4) a blend of corn/saﬄower oil (25:75) (CornSaff); and (5) a blend of flax/saﬄower oil (60:40) (FlaxSaff). Stool samples were collected at the end of each period. DNA was extracted and amplified for 16S rRNA gene pyrosequencing. A total of 17 phyla and 187 genera were identified. While five novel oil treatments failed to alter bacterial phyla composition, obese participants resulted in a higher proportion of Firmicutes to Bacteroidetes than overweight or normal weight groups (P = 0.01). Similarly at the genus level, overall bacterial distribution was highly associated with subjects’ body mass index (BMI). Treatment effects were observed between MUFA- and PUFA-rich diets, with the three MUFA diets elevating Parabacteroides, Prevotella, Turicibacter, and Enterobacteriaceae’s populations, while the two PUFA-rich diets favored the higher abundance of Isobaculum. High MUFA content feedings also resulted in an increase of Parabacteroides and a decrease of Isobaculum in obese, but not overweight subjects. Data suggest that BMI is a predominant factor in characterization of human gut microbiota profile, and that MUFA-rich and PUFA-rich diets impact the composition of gut microbiota at lower taxonomical levels mainly in obese subjects.
The prevalence of metabolic syndrome (MetS) is dramatically increasing worldwide (Mente et al., 2009), with major attention directed on examining effects of dietary fat types, and particularly the optimal ratio of saturated fatty acid (SFA), monounsaturated fatty acid (MUFA), and polyunsaturated fatty acid (PUFA). High amounts of n-6 PUFA, mainly linoleic acid (LA) replacing the SFA contents found in present Western diets, have been shown to prevent cardiovascular diseases (Kim et al., 2013). Other types of PUFA including n-3 PUFA-rich diets also resulted in weight loss in humans (Madsen et al., 2005). The Mediterranean diet has long been associated with decreased prevalence of obesity and diabetes, and increased longevity (Bergouignan et al., 2009), partly because of the favorable effects of its MUFA content from olive oil, mainly oleic acid (OA). Therefore, although current dietary recommendations suggest reductions in intakes of SFA, they fall short of providing the optimal ratios of unsaturated fatty acids needed to prevent MetS development. Novel modified oils, designed as sources of blended MUFA and PUFA contents, are expected to significantly improve risk factors of MetS, including central obesity, body mass index (BMI), blood pressure, impaired glucose tolerance, lipid profile, age and lifestyle (Alberti et al., 2005). Recent evidence suggest that the composition of gut microbiota may play a role in metabolism and adiposity, implicated in MetS risk (Ley, 2010).
Each individual houses a specific and diverse composition of gut microbiome, which potentially impacts their health. Differences in the gut microbiota patterns have been observed between normal weight and obese animals (Ley et al., 2005; Turnbaugh et al., 2008; Myles et al., 2014). In humans, the gut microbiota has also been reported to impact the physiological state of obesity (Turnbaugh et al., 2009a). Although the links between the composition of gut microbiota and specific conditions associated with obesity are not clearly understood, an increased ratio of two dominant bacterial phyla, Firmicutes to Bacteroidetes appears to associate with obese-type humans and further correlate to obesity and MetS development (Ley et al., 2006). Therefore, manipulation of gut microbiome composition in obese subjects may have potential therapeutic implications for maintaining human health.
The composition of the gut microbiota can change dramatically in response to long-term dietary intake as nutrients obtained through foods are essential to these bacteria. A recent study comparing the gut microbiota of children from Burkina Faso and Italy revealed that the ratio of two dominant phyla, Bacteroidetes and Firmicutes, were entirely different between the two groups (De Filippo et al., 2010). Here, differences in dietary patterns presumably contributed to the shifts in microbiota, which in turn influenced overall host physiology and metabolism. Studies on traditional Western diets rich in SFA but low in fiber content demonstrated that the diets rich in SFA possessed stronger effects than unsaturated fatty acids in shifting gut microbiota profiles toward those in obese individuals (de Wit et al., 2012; Liu et al., 2012). On the other hand, very few studies to date have investigated the impacts of novel vegetable oils rich in MUFA or PUFA on the gut microbiota in a longer-term human clinical trial.
Therefore, the cause-and-effect relationships between any shift of the gut microbiota, the obese state itself, and the diet remain unknown (Conterno et al., 2011).
Here, we hypothesized that healthy dietary oil treatments with diverse fatty acid compositions, particularly those rich in MUFA or PUFA, may influence the bacterial composition in human gut. Objectives of present study were to investigate (1) whether significant differences in the gut microbiota composition (as represented by the feces Krause and Khafipour, 2011) could be produced in response to MUFA or PUFA oil blends; (2) how the obese status influences the gut bacterial communities during the interventions; and (3) whether the change in gut microbiota composition is correlated with the change in biomarkers of MetS.
Materials and Methods
The Canola Oil Multi-center Intervention Trial (COMIT) was conducted at Richardson Centre for Functional Foods and Nutraceuticals (RCFFN) at University of Manitoba, the Institute of Nutraceuticals and Functional Foods (INAF) at Laval University, and the Department of Nutritional Sciences at Pennsylvania State University (PSU) between September 2010 and March 2012. The intervention studies were reviewed and approved by institutional ethics boards in respective participating universities. Written consents were obtained from all subjects as prescribed by Research Ethics Boards at all clinical centers. The protocol of the present sub-study of COMIT was approved by Biomedical Research Ethics Board (BREB) at University of Manitoba. The clinical trial was registered with clinicaltrials.gov (NCT01351012).
The COMIT study was designed as a random, controlled, double-blind, crossover clinical trial on volunteers at risk of MetS, as previously described (Senanayake et al., 2014). Adult men and women with at least one of the following cardiovascular risk factors were recruited for the study: waist circumference ≥94 cm for men and ≥80 cm for women, triglyceride (TG): ≥ 1.7 mmol/L, HDL cholesterol (HDL-C): <1 mmol/L (men) or <1.3 mmol/L (women), blood pressure: ≥130 mmHg (systolic) and/or ≥ 85 mmHg (diastolic) and glucose: ≥ 5.5 mmol/L. Participants with thyroid disease, diabetes mellitus, kidney disease, liver disease, current smokers, or people drinking more than two alcoholic beverages per week were excluded from the study. During the trial, volunteers who took any medication known to affect endothelial function during the trial were released from protocol, but those who were taking blood pressure medication with a constant dose were included. All participants provided written consent before clinical trial started. The endpoint data for serum lipid variables [Total Cholesterol (TC), HDL-C, LDL Cholesterol (LDL-C), TG] for the 25 subjects included in the microbiome study were extracted from previously reported COMIT data (Jones et al., 2014), reanalyzed (Supplementary Figure S7) and reinterpreted in the context of their correlations with microbial changes evaluated in this research.
The five-phase randomized full-feeding study design provided subjects with a consistent, individual weight-maintaining diet with a fixed 7-day rotation isocaloric menu of three meals and two snacks a day, including 50% carbohydrate, 15% protein and 35% fat of total energy of 3000 Kcal per day. A daily intake of 60 g dietary oils was equally distributed to two identical sizes of beverage shakes at breakfast and supper. Five oil treatments included: (1) canola oil [Canola; 63% MUFA, 20% LA, 10% α-linolenic acid (ALA)]; (2) DHA enriched canola-oil [CanolaDHA; 64% MUFA, 13% LA, 6% docosahexaenoic acid (DHA)]; (3) high OA canola oil (CanolaOleic; 72% MUFA, 15% LA, 2% ALA); (4) a blend of corn oil/saﬄower oil (CornSaff; 18% MUFA, 69% LA); and (5) a blend of flax oil/saﬄower oil (FlaxSaff; 18% MUFA, 38% LA, 32% ALA) (Table 1). All five diets were low in SFA, and SFA was replaced with different combinations of unsaturated fatty acids. Three canola-based diets were rich in MUFA while two saﬄower oil blends were high in n-3 PUFA or n-6 PUFA. To compare the impacts of MUFA-rich diets (Canola, CanolaDHA, and CanolaOleic) with those of PUFA-rich diets (CornSaff and FlaxSaff) on the gut microbiota, a contrast analysis was performed. Comparisons of subgroups between CanolaDHA and CanolaOleic, and CornSaff and FlaxSaff were also conducted. All the meals and shakes were prepared fresh in the metabolic kitchen. Participants were instructed to visit the research site and consume at least one meal containing one dose of treatment beverage (usually breakfast) under the supervision of the clinical staff in the cafeteria on weekdays, while the remaining meals and weekend proportions were packed for off-site consumption. No outside food was allowed to consume during the intervention period. Compliance was assessed by clinical staff by the rate of completion of meals provided under supervision as well as by the presence of food not consumed in their returned meal bags packed for off-site consumption. Both participants and study coordinators were blinded to the treatments. Each treatment phase lasted 30 days in duration and separated with 4 weeks washout periods.
TABLE 1. Fatty acid composition of five dietary oil treatments1 (g) consumed at 60 g/d (35% of energy intake based on 3000 kcal/d).
Twenty-five participants from RCFFN (Winnipeg, MB, Canada) arm of the COMIT trial agreed to provide 5 g stool samples at the end of each intervention phase for the present study (Supplementary Tables S1 and S2). Subjects were instructed to collect their own samples from only one bowel movement in the privacy of their home the day before they visited RCFFN on either day 29 or day 30 of their endpoint visits. Subjects were directed to store fecal specimens immediately after collection at -20°C freezer until collector tubes were handed over to clinical coordinators at RCFFN. Fecal specimens subsequently labeled and stored at -80°C until further analysis.
The frozen fecal specimens in collector tubes were wrapped in individual zipper bags and placed in a cooler with ice packs at RCFFN. After transport to the Gut Microbiome and Large Animal Biosecurity Laboratories at Department of Animal Science, University of Manitoba, fecal samples were thoroughly thawed and homogenized at room temperature. Genomic DNA was then extracted from approximately 250 mg fecal mass using ZR Fecal DNA Kit (D6010, Zymo Research Corp., Orange, CA, USA), which included a bead-beating step to lyse microbial cells, according to the manufacturer’s protocol. To match the concentration requirement for pyrosequencing, DNA quantity was determined under a Beckman DU/800 spectrophotometer (Beckman Coulter, Inc., Fullerton, CA, USA). Genomic DNA was normalized to achieve a concentration of 20 ng/μL and quality checked by PCR amplification of the 16S ribosomal RNA (rRNA) gene, 27f (5′-GAAGAGTTTGATCATGGCTCAG-3′) and 342r (5′-CTGCTGCCTCCCGTAG-3′) (Sepehri et al., 2007; Khafipour et al., 2009). Amplicons were verified by agarose gel electrophoresis.
Approximately 200 ng aliquots of high quality, inhibitor-free genomic DNA were sent to the Research and Testing Laboratories (Lubbock, TX, USA) for bacterial rRNA gene tag-encoded GS FLX-Titanium amplicon pyrosequencing (Dowd et al., 2008). Briefly, a mixture of Hot Start, HotStar high-fidelity Taq polymerases, and Titanium reagents were used to perform a one-step PCR (35 cycles). Pyrosequencing was performed using a 454 GS FLX-Titanium Sequencing System (454 Life Sciences, Roche Company, Branford, CT, USA). The primers 28f (GAGTTTGATCMTGGCTCAG) and 519r (GTNTTACNGCGGCKGCTG) were utilized to target the variable regions V1–V3 of the bacterial 16S rRNA genes. The sequencing data are uploaded into the Sequence Read Archive of NCBI1 and are accessible through accession number SRR2959981.
Sequence Classification and Diversity Analysis
Sequencing data, which are categorical data, were edited, transformed and classified as described previously (Li et al., 2012). In general, all failed sequence reads, low-quality sequence ends, tags, and non-bacterial ribosome sequences and chimeras were removed from the dataset. Using software mothur (version 1.30.2) (Schloss et al., 2009), the second round of sequence quality control and assignments of operational taxonomic unit (OTU) were performed. All sequences less than 200 bp or sequences having one or more ambiguous base or containing a homopolymer length more than 8 bp were excluded from the data set. High quality 16S rRNA gene sequences were identified and aligned using the comprehensive rRNA database Silva (Pruesse et al., 2007) to reduce the noise from pyrosequencing data. The remaining sequence position columns and sequences were used to build a distance matrix with a distance threshold of 0.15. Using the furthest neighbor algorithm with a cutoff of 95% similarity, these sequences were clustered to OTU.
Within community diversity (α-diversity) was conducted based on OTU counts to evaluate the biodiversity of the bacterial population in the fecal samples at the genus level. Richness indices, Chao1 and abundance based coverage estimation (ACE), Shannon index and Simpson diversity were calculated to estimate the number of OTU that were present in each individual sample. Appropriate statistical considerations were used to exclude samples with extremely high or low diversity as outliers. A rarefaction curves, which allows the calculation of sequence richness, were also performed for five treatment groups and BMI groups using mothur (Schloss et al., 2009) based on randomly re-sampling the pool of sequence numbers.
To evaluate the differences in community composition among different treatments and BMI populations, β-diversity was measured by calculating the Bray-Curtis distance (Clarke et al., 2006). Principal coordinate analysis (PCoA) was applied on resulting distance matrices to generate two-dimensional plots using PRIMER v6 software (Warwick and Clarke, 2006). Permutational multivariate analysis of variance (PERMANOVA) was performed to assess significant differences of β-diversity among the treatments and BMI groups on the basis of distance measures using permutation methods (Anderson, 2005).
Based on the outcomes from mothur, the phyla and genera in relatively low abundance (under 0.1% of community) were removed from further analysis. SAS version 9.2 (SAS Institute Inc., Cary, NC, USA) was used for data analysis. Normality of residuals for α-diversity indices was tested by Kolmogorov–Smirnov tests and histograms were used to test the Gaussian nature of the dataset. The significance of differences between treatments was analyzed using mixed model with treatment followed by Tukey adjustment.
Normal distribution was tested by Kolmogorov–Smirnov tests and histograms, and Box-Cox power transformation macro was used if necessary when analyzing statistical differences among treatments or contrasts at phylum and genus levels. Depending on the distribution of the residuals, MIXED or GLIMMIX models in SAS with mixed-effect analysis of variance with treatment as fixed effect, and subject as a random effect followed by Tukey adjustment, were conducted to estimate the significant differences. Differences between treatments were considered significant at P < 0.05 while trends were observed at P < 0.1.
The non-parametric Spearman correlation matrix was generated and Spearman’s ρ and P-values were calculated to investigate the correlations between fecal microbiome and serum lipid profiles using the PAST software (PAleontological Statistics; Hammer et al., 2001). The results were visualized over heatmap using core plot package of R (R Core Team, 2013). Taxa with relative abundance greater than 0.1% of community were included in the analysis. Strong and weak correlations were presented by dark and light colors, respectively. The Spearman’s Rho denotes whether the correlation between the taxa of interest and the selected parameter was positive (closer to 1, blue squares) or negative (closer to -1, red squares).
Partial least square discriminant analysis (PLS-DA) was performed on bacterial genera to identify the effects of dietary oil treatments on the bacterial community, according to our previously described analytical method with modifications (Li et al., 2012). For this analysis, data were scaled using Unit Variance in SIMCA-P+ 13.0 (Umetrics, Umea, Sweden). Cross-validation and permutation testing were conducted to determine significant PLS components in the optimal model. To avoid over parameterization of the model, a cut-off level was set for variable influence on projection (VIP) value in each model to improve R2 and Q2 with better prediction of the bacterial distribution. Significant differences across treatments or BMI groups were expressed in scatter and score plots, according to the PLS regression coefficients. Each individual genus 95% confidence interval determined the positive or negative correlations with dietary treatments or BMI groups.
Sample Assessment by Pyrosequencing
After the 5-phase clinical trial, a total of 66 samples were collected from 25 subjects. Missing samples resulted from participants either forgetting to collect samples or failing to provide samples at expected sampling dates (Supplementary Tables S1). In total, 209,490 sequences were first generated. The minimum, median and maximum lengths of sequences were 250, 393, and 698 bp, respectively. Screening, filtering, and pre-clustering processes resulted in 142,887 sequences.
Alpha and β-Diversity Analyses
Bacterial richness and diversity in individual samples under different oil treatments and BMI were calculated (Table 2). Oil treatments and BMI had no significant impact on Chao1, ACE, Shannon and Simpson indices of α-diversity (Figure 1A and Table 2). However, the rarefaction curves showed higher richness and diversity in overweight and obese subjects compared to the normal weight participants (Figure 1B). Similarity or differences in the gut microbiota among different oil treatments and BMI (β-diversity) were compared using PCoA and PERMANOVA analyses of Bray-Curtis distances (Supplementary Figures S1–S6). Comparisons included MUFA vs. PUFA (P = 0.84), normal weight vs. overweight (P = 0.11), normal weight vs. obese (P < 0.01), overweight vs. obese (P = 0.03), MUFA vs. PUFA within overweight group (P = 0.57), and MUFA vs. PUFA within obese group (P = 0.99). The β-diversity did not change among individual oil treatments.
FIGURE 1. (A,B) Rarefaction analysis for the observed species. The rarefaction curve is generated using Chao1 richness estimator according to oil treatments (A) or subjects’ body mass index (BMI; B). Samples have been rarified at an even depth of 21,000 sequences per oil treatment or BMI group.
Overall Gut Microbiota Composition
A total of 15 phyla were identified in all the samples from 99.8% of all sequences (Supplementary Tables S3), of which four phyla were considered as abundant within the community (larger than 1%), including Firmicutes, Bacteroidetes, Actinobacteria, and Proteobacteria (Table 3). The abundance of Aquificae, Fusobacteria, and Verrucomicrobia were between 0.1 and 1% of community. The remaining eight phyla, which were in low abundance (under 0.1%) included Acidobacteria, Chrysiogenetes, Cyanobacteria, Deinococcus-Thermus, Nitrospira, Synergistetes, TM7, and Tenericutes. Overall, the phylum distribution did not fluctuate across the five oil treatments or among MUFA vs. PUFA groups. However, BMI status had an impact on the microbiota composition at the phylum level with obese group having greater proportion of Firmicutes (P = 0.02) compared to the combined group of normal and overweight subjects (data not shown). The average ratio of Bacteroidetes to Firmicutes was 0.15 across five diets with no difference among dietary interventions (data not shown).
At the genus level, 187 genera were determined using 16S rRNA gene pyrosequencing. Majority of sequences were identified at the genus (g.) level while some sequences could only be classified up to phylum (p.), class (c.), order (o.), or family (f.) levels. Supplementary Tables S4 and S5 present the distribution of taxa with a relative abundance >0.1% and <0.1% of community, respectively, among diets. The abundance of g. Faecalibacterium tended to differ across treatments (P = 0.06) where CanolaOleic feeding was associated with the highest level of Faecalibacterium and CanoalDHA feeding with the lowest level. Data also showed a tendency for higher proportion of g. Parabacteroides (P = 0.09) but lower proportion of g. Isobaculum (P = 0.08) when the three MUFA-rich diets compared to two PUFA-rich diets. Genus Blautia was observed to favor the CanolaOleic feeding compared to the CanolaDHA diet (P = 0.09). No differences were observed between n-6 PUFA-rich treatment CornSaff and n-3 PUFA-rich diet FlaxSaff.
The PLS-DA analysis (given for the cutoff value of 1.0) showed further differences at the lower taxonomical levels. Genera Parabacteroides, Prevotella, Turicibacter, and f. Enterobacteriaceae were correlated to MUFA-rich diets, while g. Isobaculum was correlated to PUFA-rich diets (R2 = 0.43, Q2 = 0.07; Figure 2). Comparison between CanolaDHA and CanolaOleic indicated that CanolaDHA was correlated to f. Lachnospiraceae and p. Firmicutes, whereas CanolaOleic was associated with g. Faecalibacterium and Coprobacillus (R2= 0.78, Q2= 0.45; Figure 3). The comparison between two PUFA-rich diets showed CornSaff treatment had an impact on the populations of g. Eggerthella, Slackia, Soehngenia, Anaerostipes, Robinsoniella, Phascolarctobacterium, while FlaxSaff failed to have such significant impact (R2= 0.67, Q2= 0.22; data not shown).
FIGURE 2. Partial least square discriminant analysis (PLS-DA) loading plot based on the relative abundances of bacterial genera in the gut microbiota and their association with MUFA-rich or PUFA-rich diets. The presenting genera are chosen at VIP cutoff of 1.0 and colored according to their corresponding phyla. The size of circles is indicative of taxa abundance. Taxa closer to MUFA or PUFA indicate positive association with either treatment group. The R2 (=0.43) and Q2 (=0.07) estimates were calculated based on three components. While majority of sequences were identified at the genus level (G), some could only be affiliated to phylum (P), class (C), order (O), or family (F) levels.
FIGURE 3. Partial least square discriminant analysis loading plot based on the relative abundances of bacterial genera in the gut microbiota and their association with CanolaDHA or CanolaOleic diets. The presenting genera are chosen at VIP cutoff of 1.0 and colored according to their corresponding phyla. The size of circles is indicative of taxa abundance. Taxa closer to CanolaDHA or CanolaOleic diets indicate positive association with either treatment group. The R2 (=0.78) and Q2 (=0.45) estimates were calculated based on three components. While majority of sequences were identified at the genus level (G), some could only be affiliated to phylum (P), class (C), order (O), or family (F) levels.
Gut Microbiota Profiles Within Individual Obese Status
At the phylum level, a significantly higher proportion of Firmicutes was observed in obese group (BMI > 30) compared to the combined group of normal weight (BMI < 25) and overweight subjects (BMI between 25 and 30; P = 0.02; data not shown). At the genus level, PLS-DA analysis confirmed a significant difference in the composition of bacteria among three BMI phenotype groups (R2 = 0.60, Q2 = 0.32; Figure 4).
FIGURE 4. Partial least square discriminant analysis score scatter plot of the gut microbiota comparing subjects with different body mass index (BMI). The distribution of fecal samples (n = 66) indicated a significant difference in the composition of bacterial genera among obese (BMI > 30), overweight (BMI between 25 and 30), and normal weight (BMI < 25) groups. The R2 (=0.60) and Q2 (=0.32) estimates were calculated based on three components.
Since the number of subjects in the normal weight category was insufficient, the analysis on contrasts and subgroups only focused on two representative populations, overweight and obese groups. In overweight people, g. Streptococcus, Tepidimicrobium, Robinsoniella, and Turicibacter were correlated to MUFA-rich diets, while g. Coriobacterium and Mogibacterium were associated with PUFA-rich diets (R2 = 0.69, Q2= 0.26; Figure 5). In the comparison between CanolaDHA and CanolaOleic, g. Adlercreutizia, Coriobacterium, Alistipes, and Robinsoniella showed strong correlations with CanolaDHA while g. Lactobacillus was associated with CanolaOleic (R2 = 0.90, Q2= 0.60; data not shown). In the comparison between the two PUFA-rich diets, CornSaff was highly associated with g. Adlercreutizia while FlaxSaff was correlated with g. Collinsella, Barnesiella, Streptococcus, Roseburia, Coprobacillus, and family Peptostreptococcaceae (R2 = 0.98, Q2= 0.74; data not shown).
FIGURE 5. Partial least square discriminant analysis loading plot based on the relative abundances of bacterial genera in the gut microbiota and their association with MUFA-rich or PUFA-rich diets in overweight population (BMI between 25 and 30). The presenting genera are chosen at VIP cutoff of 1.0 and colored according to their corresponding phyla. The size of circles is indicative of taxa abundance. Taxa closer to MUFA or PUFA indicate positive association with either treatment group. The R2 (=0.69) and Q2 (=0.26) estimates were calculated based on three components. While majority of sequences were identified at the genus level (G), some could only be affiliated to phylum (P) class (C), order (O), or family (F) levels.
In obese subjects, g. Parabacteroides, Prevotella, Flexithrix, Fusibacter, f. Enterobacteriaceae, and p. Firmicutes were correlated to MUFA-rich diets, but no specific taxa was associated with PUFA-rich diets (R2 = 0.66, Q2= -0.20; Figure 6). In comparison of CanolaDHA and CanolaOleic, only g. Parasutterlla correlated with CanolaDHA (R2 = 0.91, Q2= 0.29; data not shown). Between the two PUFA-rich diets, g. Collinsella, Hydrogenobaculum, and Parabacteroides were impacted by the CornSaff, while g. Clostridium was correlated to the FlaxSaff diet (R2 = 0.98, Q2= 0.63; data not shown).
FIGURE 6. Partial least square discriminant analysis (PLS-DA) loading plot based on the relative abundances of bacterial genera in the gut microbiota and their association with MUFA-rich or PUFA-rich diets in obese population (BMI > 30). The presenting genera are chosen at VIP cutoff of 1.0 and colored according to their corresponding phyla. The size of circles is indicative of taxa abundance. Taxa closer to MUFA or PUFA indicate positive association with either treatment group. The R2 (=0.66) and Q2 (=-0.20) estimates were calculated based on three components. While majority of sequences were identified at the genus level (G), some could only be affiliated to phylum (P) class (C), order (O), or family (F) levels.
Correlations between Serum Lipid Variables and Gut Microbiota Profiles
The endpoint data for serum lipid variables of the 25 subjects included in this sub-study were extracted from previously reported COMIT data (Jones et al., 2014) and reanalyzed. All changes in lipid profiles followed the same trends of the entire study population (Supplementary Figure S7). The correlations between lipid profiles and bacterial phyla are presented in Supplementary Tables S6. Across treatments, serum TG was negatively correlated with Aquificae (Spearman’s ρ = -0.27, P = 0.02) but positively correlated with Cyanobacteria (Spearman’s ρ = 0.24, P = 0.05). In contrast, LDL-C was positively correlated with p. Proteobacteria (Spearman’s ρ = 0.28, P = 0.01) and similar tendency was observed for TC (Spearman’s ρ = 0.21, P = 0.08). HDL-C, however, was positively correlated with Verrucomicrobia (Spearman’s ρ = 0.21, P = 0.08). In CanolaDHA treatment, TC levels trended to negatively correlate with Bacteroidetes (Spearman’s ρ = -0.51, P = 0.06) and positively correlated with Firmicutes (Spearman’s ρ = 0.55, P = 0.04), resulting in negative correlation trend with the ratio of Bacteroidetes to Firmicutes (Spearman’s Rho = -0.51, P = 0.06). In CornSaff treatment, TC levels were positively correlated with Bacteroidetes (Spearman’s ρ = 0.64, P = 0.02) and Bacteroidetes to Firmicutes ratio (Spearman’s ρ = -0.65, P = 0.02). The correlations between serum lipid variables and bacterial composition at the lower taxonomical levels are presented in Figure 7.
FIGURE 7. Correlations between fecal microbiota compositions and serum lipid profiles. The non-parameteric Spearman correlation coefficient and P-values were presented in a heatmap format. Strong and weak correlations are presented by dark and light colors, respectively. The intensity bar shows positive and negative correlation from blue to red, respectively. ∗P < 0.1; ∗∗P < 0.05; ∗∗∗P < 0.01. F.B, ratio of Firmicutes to Bacteroidetes.
We investigated the impact of five different healthy unsaturated fatty acid oil treatments with a dietary fat content of less than 35% of energy intake on gut microbiota profiles of volunteers with risks of MetS. Our experimental diets included three MUFA-rich (Canola, CanolaDHA, and CanolaOleic) and two PUFA-rich (CornSaff and FlaxSaff) diets and were designed to maintain the body weight during the treatment periods. Major findings on post-treatment outcomes indicate that diet, metabolic profiles, obesity state, and individual responses are all contributing at different degrees to the shifts in the gut microbiota compositions. Overall, the impact of our healthy unsaturated fatty acid oil treatments on the microbiota composition was relatively minor at the phylum level and mainly associated with microbiota shifts at the lower taxonomical levels, whereas BMI contributed to a significant shift at the phylum level with greater Firmicutes and less Bacteroidetes in obese group compared to the combined group of normal and overweight subjects.
Several studies have investigated the effects of high fat diets mainly from SFA on the gut microbiota composition [8, 11, 14, 30, 31, 32]. Turnbaugh et al. (2008) showed that diet-induced obesity can quickly develop with altered microbiota profiles in mice fed prototypic Western diets with high SFA content. In a study comparing different types of high fat diets on the profile of gut bacteria in a mouse model, Liu et al. (2012) observed that consumption of a SFA-rich diet resulted in a significant decrease in the abundance of p. Bacteroidetes compared to either n-3 PUFA-rich or n-6 PUFA-rich diets. Similarly, a clinical study by Ley et al. (2006) described that the relative proportion of p. Bacteroidetes decreased in obese people compared with lean subjects, and this proportion also increased when the obese individual lost weight due to either carbohydrate-restricted or fat-restricted diets. Nielsen et al. (2007) reported that supplementation of fish oil for 1 month in human infants influenced the intestinal microbiota composition suggesting that dietary n-3 fatty acids may influence the bacterial adhesion to the intestinal mucosa, which in turn associates to the immune response and fat accumulation. Furthermore, Andersen et al. (2011) reported a significant increase in p. Bacteroidetes following the n-6 PUFA treatment. In the present study we did not have a control low fat or a control high fat treatment and the proportion of dietary fat content in all five treatments was less than 35% of total energy intake, which is the recommended level of fat intake in North America to maintain body weight (Grundy et al., 2002). As a result, our treatments were not considered high fat diet and we did not expect to see major shifts at the phylum level in the bacterial composition among n-3 PUFA (FlaxSaff), n-6 PUFA (CornSaff) and DHA-rich (CanolaDHA) treatments, or the combination of MUFA-rich vs. PUFA-rich diets. As expected the differences observed were at the lower taxonomical levels, for example MUFA-rich diets increased the populations of g. Prevotella and Parabacteroides within p. Bacteroidetes, and their impact was more pronounced in the obese people compared to normal and overweight subjects.
The individuals’ BMI has a major impact on gut microbiota. A study by Ley et al. (2005) used obese mice with ob/ob, lean ob/+, and wild-type siblings, and their ob/+ mothers to investigate the impact of obese genotypes and phenotypes on the gut microbiota. Results from a 16S rRNA gene sequence observation showed that ob/ob mice had a decreased abundance of Bacteroidetes and an increased proportion of Firmicutes. Authors also reported that the ratio of Bacteroidetes to Firmicutes dynamically reflects the overall weight condition in mice model. A follow-up study to extend these observations to humans showed that the relative proportion of Bacteroidetes is reduced in obese participants compared to the lean controls (Ley et al., 2006). However, it is very unlikely that the compositional shifts in Firmicutes to Bacteroidetes ratio can be solely responsible for development of obese phenotype. Other researchers, for example (Schwiertz et al., 2010), have failed to find similar trends when comparing the Firmicutes to Bacteroidetes ratio between the microbiota of obese and lean individuals, suggesting that the relationship between MetS, obesity, nutrition, and the microbiota is complex and multifactorial. Bacteroidetes genomes have a significantly higher level of functional diversity and, as demonstrated by the oft-cited study conducted by Turnbaugh et al. (2009a), the gut microbiome of human subjects from obese phenotype is functionally less diverse compared to the Bacteroidetes-rich gut microbiota of lean subjects. In the latter study, the majority of genes that were discriminating between the two phenotypes were annotated to carbohydrate, lipid and amino-acid metabolism pathways, with 75% of the obesity-enriched genes belonging to Actinobacteria and 42% of the lean-enriched genes to Bacteroidetes. These shifts in the functional properties of microbiota can change the energy harvesting capacities of the host via different mechanisms, for example, by altering the amount and profile of microbiota-derived short chain fatty acids (SCFAs), and/or, by modulating the expression of host genes that affect energy deposition in adipocytes (Bäckhed et al., 2007; Murphy et al., 2010). In the present study, subjects’ recruitment was based on the presence of at least one of the MetS risk factors including waist circumference, blood pressure, TG, HDL-C, and glucose levels. As a result, we didn’t have a balanced number of subjects at different BMI groups, and thus, we only observed a tendency of higher ratio of Bacteroidetes to Firmicutes in normal and overweight people compared to that of obese subjects. Nevertheless, when the obese states were classified as subgroups, the correlation coefficients were significantly improved indicating higher predictability of our models for detection of microbiota pattern in response to MUFA- or PUFA-rich diets.
Comparing the MUFA- and PUFA-rich diets, one interesting finding was the shift in g. Parabacteroides and Isobaculum populations suggesting that additional dietary OA can result in an increase in the population of Parabacteroides and a decrease in Isobaculum. When we chose the overweight and obese groups for subgroup analysis, contrasting MUFA vs. PUFA, a positive correlation between Parabacteroides and MUFA-rich diets was found to exist only in obese but not in overweight subjects. Further analysis comparing across CanolaOleic and CornSaff or FlaxSaff oil phases also showed consistent outcomes that this potential correlation with Parabacteroides was significant but was only observed in obese but not overweight volunteers. Thus, our study indicated that the MUFA-rich diets (high OA contents) can increase the population of Parabacteroides in human gut, and people with obese status may express a stronger response to such shift. Unfortunately, there is not much known about the function of Parabacteroides and Isobaculum and further investigation on the potential roles of these genera in development of obesity is needed.
Our results showed that g. Faecalibacterium population was higher in CanolaOleic treatment compared to CanolaDHA. Among different species of g. Faecalibacterium, F. prausnitzii is considered as an indicator of intestinal health due to its anti-inflammatory effects (Carlsson et al., 2013). Based on our findings, we can, therefore, hypothesize that additional intake of OA may increase the anti-inflammatory properties in human intestines by increasing the population of F. prausnitzii although further investigation is required.
The gut microbiome profile in humans is host-specific. Fecal communities in the human gut can change dramatically in response to new diets. Recent studies have shown that gut microbiome can rapidly and reproducibly respond to short-term dietary changes (David et al., 2014). These studies suggest that the human gut microbiome can rapidly switch between herbivorous and carnivorous functional profiles, which might be reflective of past selective pressures during human evolution (David et al., 2014). Turnbaugh et al. (2009b) using a humanized mice model showed that switching from a low-fat, plant polysaccharide-rich diet to a high-fat/high-sugar “Western” diet shifted the structure of the microbiota within a single day, changed the representation of metabolic pathways in the microbiome, and altered microbiome gene expression. The authors also indicated that the gut microbiota composition was stabilized 7 days following the dietary switch. As such, the 4-week dietary interventions used in this study should have been sufficient to stabilize diet-induced changes in the microbiota. It is clear now that each individual is relatively unique in terms of bacterial species of microbiota in their gut, no matter what obese state the individual belongs to, or what dietary oil treatments are provided. This may partly explain the variations in the microbiota profiles among the participants of our study and the variations in their responses to the dietary treatments.
Evidence showed dietary DHA oil intends to increase LDL-C levels in humans (Theobald et al., 2004; Jacobson et al., 2012), although the reason remains unclear. We observed similar elevated LDL-C levels after CanolaDHA treatment in the present study, which was significantly correlated with the increase in Firmicutes and the decrease in Bacteroidetes within this treatment. These observations are supported by a mouse study (Mujico et al., 2013), which reported that a diet supplemented with n-3 fatty acids (Eicosapentaenoic acid + DHA) significantly increased the abundance of Firmicutes and reduced the percentage of Bacteroidetes compared to a diet supplemented with OA. It can be speculated that the changes in metabolic parameters after DHA oil intake could be the result of interactions between gut microbiota and DHA metabolites potentially through the enterohepatic circulation of bile salts (Yokota et al., 2012; Sayin et al., 2013). However, our study was not designed to test this hypothesis and further investigation on the bile excretion analysis is required to proof or disproof this hypothesis.
All being said, our study also had several limitations. First, given an intervention trial period of more than 9 months of duration during which multiple samples were collected at different time points, some participants forgot to provide a fecal sample on the specific sampling day. Second, the small proportion of fecal samples might not be able to represent the total fecal materials. Even after the homogenization, the specimens could still produce large variations in the bacterial composition and, in addition, were limited to reflect the dynamic condition in the intestine. Third, the inclusion criteria of our study were based on the risks of MetS, but not the obese state. Given the lower numbers of participants with normal weight and overweight, the statistical analyses were not optimal to compare the direct variations among all participants.
In summary, the current study investigated the impact of five different healthy unsaturated fatty acid oil treatments on volunteers with risks of MetS. Results indicate that the human gut microbiota profiles strongly relate to the body obesity states, but also exhibit host-specificity. MUFA-rich and PUFA-rich diets failed to shift the composition of gut microbiota at the phylum level in a 30-day treatment period, but the populations of specific genera can slightly alter in response to the feedings due to unclear mechanisms. Further studies may clarify whether there is a cause-and-effect relationship between dietary fatty types perhaps through altered cholesterol and primary and secondary bile acid concentration, and shifts in human gut microbiota.
SP, PJ, and EK conceived and designed the experiment. SP and HK performed lab analyses. EK developed the bioinformatics and statistical models. SP and HK analyzed the data. SP, PJ, and EK wrote the manuscript.
The research was funded by Agri-Food Research and Development Initiative (ARDI) program 10-1057. The clinical trial was supported by Canola Council of Canada, Flax Council of Canada, Agriculture and Agri Food Canada, and Dow Agrosciences through the Growing Forward program of Agriculture and Agri Food Canada.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Authors are indebted to the late Dr. Denis Krause for the special contribution to the design of the study. We also would like to thank Wenhua Tang for technical assistance on human fecal sample analysis. Authors also acknowledge Shucong Li, Navjot Nandha, and Hooman Derakhshani for their helpful discussions and interpretation on statistics and bioinformatics. Special thanks to the study participants, clinical coordinators, and kitchen staff at RCFFN for their contributions in the clinical study.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.01612
FIGURE S1| Principal coordinates analysis (PCoA) of Bray-Curtis distance between the MUFA and PUFA oil treatments.
FIGURE S2| Principal coordinates analysis (PCoA) of the Bray-Curtis distance between the normal weight and overweight subjects.
FIGURE S3| Principal coordinates analysis (PCoA) of the Bray-Curtis distance between the normal weight and obese subjects.
FIGURE S4| Principal coordinates analysis (PCoA) of the Bray-Curtis distance between the overweight and obese subjects.
FIGURE S5| Principal coordinates analysis (PCoA) of the Bray-Curtis distance between the MUFA and PUFA oil treatments within overweight group.
FIGURE S6| Principal coordinates analysis (PCoA) of the Bray-Curtis distance between the MUFA and PUFA oil treatments within obese group.
FIGURE S7| Serum lipid profiles at the endpoint of five experimental diets. Significant differences between treatments are highlighted within each parameter.
TABLE S1 | List of fecal samples collected from the COMIT study.
TABLE S2 | Example of a 7-day rotation isocaloric menu during the intervention.
TABLE S3 | Summary of sequenced OTUs at each taxonomic level.
TABLE S4 | Phylogenetic composition of putative bacterial genera (above 0.1% of community) determined using 16S rRNA gene pyrosequencing.
TABLE S5 | List of bacterial genera in low abundance (below 0.1% of community).
TABLE S6 | Correlation between bacterial phyla and metabolic parameters.
Andersen, A. D., Molbak, L., Thymann, T., Michaelsen, K. F., and Lauritzen, L. (2011). Dietary long-chain n-3 PUFA, gut microbiota and fat mass in early postnatal piglet development–exploring a potential interplay. Prostaglandins Leukot. Essent. Fatty Acids 85, 345–351. doi: 10.1016/j.plefa.2011.08.004
Bäckhed, F., Manchester, J. K., Semenkovich, C. F., and Gordon, J. I. (2007). Mechanisms underlying the resistance to diet-induced obesity in germ-free mice. Proc. Natl. Acad. Sci. U.S.A. 104, 979–984. doi: 10.1073/pnas.0605374104
Bergouignan, A., Momken, I., Schoeller, D. A., Simon, C., and Blanc, S. (2009). Metabolic fate of saturated and monounsaturated dietary fats: the Mediterranean diet revisited from epidemiological evidence to cellular mechanisms. Prog. Lipid Res. 48, 128–147. doi: 10.1016/j.plipres.2009.02.004
Carlsson, A. H., Yakymenko, O., Olivier, I., Hakansson, F., Postma, E., Keita, A. V., et al. (2013). Faecalibacterium prausnitzii supernatant improves intestinal barrier function in mice DSS colitis. Scand. J. Gastroenterol. 48, 1136–1144. doi: 10.3109/00365521.2013.828773
Clarke, K. R., Somerfiled, P. J., and Chapman, M. G. (2006). On resemblance measures for ecological studies, including taxonomic dissimilarities and a zero-adjusted Bray-Curtis coefficient for denuded assemblages. J. Exp. Mar. Biol. Ecol. 330, 55–80. doi: 10.1016/j.jembe.2005.12.017
Conterno, L., Fava, F., Viola, R., and Tuohy, K. M. (2011). Obesity and the gut microbiota: does up-regulating colonic fermentation protect against obesity and metabolic disease? Genes Nutr. 6, 241–260. doi: 10.1007/s12263-011-0230-1
David, L. A., Maurice, C. F., Carmody, R. N., Gootenberg, D. B., Button, J. E., Wolfe, B. E., et al. (2014). Diet rapidly and reproducibly alters the human gut microbiome. Nature 505, 559–563. doi: 10.1038/nature12820
De Filippo, C., Cavalieri, D., Di Paola, M., Ramazzotti, M., Poullet, J. B., Massart, S., et al. (2010). Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc. Natl. Acad. Sci. U.S.A. 107, 14691–14696. doi: 10.1073/pnas.1005963107
de Wit, N., Derrien, M., Bosch-Vermeulen, H., Oosterink, E., Keshtkar, S., Duval, C., et al. (2012). Saturated fat stimulates obesity and hepatic steatosis and affects gut microbiota composition by an enhanced overflow of dietary fat to the distal intestine. Am. J. Physiol. Gastrointest. Liver Physiol. 303, G589–G599. doi: 10.1152/ajpgi.00488.2011
Dowd, S. E., Callaway, T. R., Wolcott, R. D., Sun, Y., Mckeehan, T., Hagevoort, R. G., et al. (2008). Evaluation of the bacterial diversity in the feces of cattle using 16S rDNA bacterial tag-encoded FLX amplicon pyrosequencing (bTEFAP). BMC Microbiol. 8:125. doi: 10.1186/1471-2180-8-125
Grundy, S. M., Abate, N., and Chandalia, M. (2002). Diet composition and the metabolic syndrome: what is the optimal fat intake? Am. J. Med. 113(Suppl. 9B), 25S–29S. doi: 10.1016/S0002-9343(01)00988-3
Jacobson, T. A., Glickstein, S. B., Rowe, J. D., and Soni, P. N. (2012). Effects of eicosapentaenoic acid and docosahexaenoic acid on low-density lipoprotein cholesterol and other lipids: a review. J. Clin. Lipidol. 6, 5–18. doi: 10.1016/j.jacl.2011.10.018
Jones, P. J., Senanayake, V. K., Pu, S., Jenkins, D. J., Connelly, P. W., Lamarche, B., et al. (2014). DHA-enriched high-oleic acid canola oil improves lipid profile and lowers predicted cardiovascular disease risk in the canola oil multicenter randomized controlled trial. Am. J. Clin. Nutr. 100, 88–97. doi: 10.3945/ajcn.113.081133
Khafipour, E., Li, S., Plaizier, J. C., and Krause, D. O. (2009). Rumen microbiome composition determined using two nutritional models of subacute ruminal acidosis. Appl. Environ. Microbiol. 75, 7115–7124. doi: 10.1128/AEM.00739-09
Kim, J., Li, Y., and Watkins, B. A. (2013). Fat to treat fat: emerging relationship between dietary PUFA, endocannabinoids, and obesity. Prostaglandins Other Lipid Mediat. 10, 32–41. doi: 10.1016/j.prostaglandins.2012.11.005
Ley, R. E., Backhed, F., Turnbaugh, P., Lozupone, C. A., Knight, R. D., and Gordon, J. I. (2005). Obesity alters gut microbial ecology. Proc. Natl. Acad. Sci. U.S.A. 102, 11070–11075. doi: 10.1073/pnas.0504978102
Li, R., Khafipour, E., Krause, D. O., Entz, M. H., De Kievit, T. R., and Fernando, W. G. (2012). Pyrosequencing reveals the influence of organic and conventional farming systems on bacterial communities. PLoS ONE 7:e51897. doi: 10.1371/journal.pone.0051897
Liu, T., Hougen, H., Vollmer, A. C., and Hiebert, S. M. (2012). Gut bacteria profiles of Mus musculus at the phylum and family levels are influenced by saturation of dietary fatty acids. Anaerobe 18, 331–337. doi: 10.1016/j.anaerobe.2012.02.004
Madsen, L., Petersen, R. K., and Kristiansen, K. (2005). Regulation of adipocyte differentiation and function by polyunsaturated fatty acids. Biochim. Biophys. Acta 1740, 266–286. doi: 10.1016/j.bbadis.2005.03.001
Mente, A., De Koning, L., Shannon, H. S., and Anand, S. S. (2009). A systematic review of the evidence supporting a causal link between dietary factors and coronary heart disease. Arch. Intern. Med. 169, 659–669. doi: 10.1001/archinternmed.2009.38
Mujico, J. R., Baccan, G. C., Gheorghe, A., Diaz, L. E., and Marcos, A. (2013). Changes in gut microbiota due to supplemented fatty acids in diet-induced obese mice. Br. J. Nutr. 110, 711–720. doi: 10.1017/S0007114512005612
Murphy, E., Cotter, P., Healy, S., Marques, T., O’sullivan, O., Fouhy, F., et al. (2010). Composition and energy harvesting capacity of the gut microbiota: relationship to diet, obesity and time in mouse models. Gut 59, 1635–1642. doi: 10.1136/gut.2010.215665
Myles, I. A., Pincus, N. B., Fontecilla, N. M., and Datta, S. K. (2014). Effects of parental omega-3 Fatty Acid intake on offspring microbiome and immunity. PLoS ONE 9:e87181. doi: 10.1371/journal.pone.0087181
Nielsen, S., Nielsen, D. S., Lauritzen, L., Jakobsen, M., and Michaelsen, K. F. (2007). Impact of diet on the intestinal microbiota in 10-month-old infants. J. Pediatr. Gastroenterol. Nutr. 44, 613–618. doi: 10.1097/MPG.0b013e3180406a11
Pruesse, E., Quast, C., Knittel, K., Fuchs, B. M., Ludwig, W., Peplies, J., et al. (2007). SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 35, 7188–7196. doi: 10.1093/nar/gkm864
Sayin, S. I., Wahlstrom, A., Felin, J., Jantti, S., Marschall, H. U., Bamberg, K., et al. (2013). Gut microbiota regulates bile acid metabolism by reducing the levels of tauro-beta-muricholic acid, a naturally occurring FXR antagonist. Cell Metab. 17, 225–235. doi: 10.1016/j.cmet.2013.01.003
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Senanayake, V. K., Pu, S., Jenkins, D. A., Lamarche, B., Kris-Etherton, P. M., West, S. G., et al. (2014). Plasma fatty acid changes following consumption of dietary oils containing n-3, n-6, and n-9 fatty acids at different proportions: preliminary findings of the Canola Oil Multicenter Intervention Trial (COMIT). Trials 15:136. doi: 10.1186/1745-6215-15-136
Sepehri, S., Kotlowski, R., Bernstein, C. N., and Krause, D. O. (2007). Microbial diversity of inflamed and noninflamed gut biopsy tissues in inflammatory bowel disease. Inflamm. Bowel Dis. 13, 675–683. doi: 10.1002/ibd.20101
Theobald, H. E., Chowienczyk, P. J., Whittall, R., Humphries, S. E., and Sanders, T. A. (2004). LDL cholesterol-raising effect of low-dose docosahexaenoic acid in middle-aged men and women. Am. J. Clin. Nutr. 79, 558–563.
Turnbaugh, P. J., Backhed, F., Fulton, L., and Gordon, J. I. (2008). Diet-induced obesity is linked to marked but reversible alterations in the mouse distal gut microbiome. Cell Host Microbe 3, 213–223. doi: 10.1016/j.chom.2008.02.015
Turnbaugh, P. J., Ridaura, V. K., Faith, J. J., Rey, F. E., Knight, R., and Gordon, J. I. (2009b). The effect of diet on the human gut microbiome: a metagenomic analysis in humanized gnotobiotic mice. Sci. Transl. Med. 1:6ra14. doi: 10.1126/scitranslmed.3000322
Keywords: monounsaturated fatty acid (MUFA), polyunsaturated fatty acid (PUFA), blended oil, gut microbiota, BMI
Citation: Pu S, Khazanehei H, Jones PJ and Khafipour E (2016) Interactions between Obesity Status and Dietary Intake of Monounsaturated and Polyunsaturated Oils on Human Gut Microbiome Profiles in the Canola Oil Multicenter Intervention Trial (COMIT). Front. Microbiol. 7:1612. doi: 10.3389/fmicb.2016.01612
Received: 01 April 2016; Accepted: 26 September 2016;
Published: 10 October 2016.
Edited by:George Tsiamis, University of Patras, Greece
Reviewed by:Spyridon Ntougias, Democritus University of Thrace, Greece
Dimitris Tsaltas, Cyprus University of Technology, Cyprus
Copyright © 2016 Pu, Khazanehei, Jones and Khafipour. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.