ORIGINAL RESEARCH article
Sec. Microbial Symbioses
Volume 11 - 2020 | https://doi.org/10.3389/fmicb.2020.564160
High-Fat Diets Led to OTU-Level Shifts in Fecal Samples of Healthy Adult Dogs
- 1Department of Animal Science, Iowa State University, Ames, IA, United States
- 2Interdepartmental Microbiology Graduate Program, Iowa State University, Ames, IA, United States
- 3Department of Veterinary Microbiology and Preventive Medicine, Iowa State University, Ames, IA, United States
High fat diets have been reported to negatively affect the microbiota in both mice and humans. However, there is a lack of studies in canine models. The variation among the gastrointestinal (GI) tract anatomy/physiology and typical diet compositions of these animal species may lead to vastly different results. Due to the large inclusion rate of dietary fat in pet food, it is critical to understand its effects in a canine model. Therefore, the study objective was to report the effects of high fat, low carbohydrate diets on the fecal microbiota in healthy adult dogs. Eight adult beagles were randomly assigned to one of four dietary treatments within each 15-day period of a replicated 4x4 Latin Square design. Diets contained 32% (T1), 37% (T2), 42% (T3), and 47% (T4) fat. T2, T3, and T4 were created by adding increasing levels of canola oil to T1, a commercially manufactured canned canine diet, which served as the control diet. Fresh fecal samples were collected during the last 5 days of each period for microbial analysis. DNA was extracted from fecal samples and paired-end 16S rRNA gene amplicon sequencing was performed using the Illumina MiSeq platform. When comparing whole microbial communities using PERMANOVA, no significant differences were observed among treatments (P = 0.735). Individual OTUs were analyzed using the GLIMMIX procedure of SAS with fixed effects of diet and room, and the random effects of period and animal. Out of the 100 most abundant individual OTUs, 36 showed significant differences in abundance based on treatment (q < 0.05). Overall, OTUs assigned to genera related to fat digestion increased while OTUs assigned to genera involved in carbohydrate digestion decreased. In conclusion, the microbial community adapted to dietary intervention without jeopardizing the health of the animals, evaluated by body condition score, fecal characteristics, and blood parameters.
The increased interest in the gastrointestinal (GI) tract microbiota in humans has extended to companion animals. This may be due to the idea that a balanced relationship between the GI microbes and the host animal is critical for host health (Mackie et al., 1999; Hooper et al., 2001). The GI tract microbiota is comprised of thousands of interdependent and/or competing microbial species (Eckburg et al., 2005; Ley et al., 2008; Spor et al., 2011), many of them are still not fully characterized (Hand et al., 2013). The GI tract microbiota can benefit the host in many ways; it can enhance metabolic capabilities, protect against pathogens, develop the immune system, and modulate gastrointestinal development (Backhed et al., 2005; McKenney and Pamer, 2015; Rooks and Garrett, 2016; Wernimont et al., 2020). In addition, GI tract microorganisms contain enzymes that digest fiber and carbohydrates that cannot be digested by the host, producing for example short chain fatty acids (SCFA) (Sunvold et al., 1995), which can be used as an additional energy source for the host. SCFA account for ~10% of human caloric requirement and ~80% of maintenance energy for ruminants (Bergman, 1990). Unlike other species, dogs do not rely heavily on microbial fermentation to meet daily energy requirements, even when fed high fiber diets (Swanson et al., 2010; Hooda et al., 2012; Deng and Swanson, 2015). Even though dogs do not rely on this energy source, a balanced microbiota is nevertheless critical for GI health (Swanson et al., 2010). A disruption, or dysbiosis, of the GI tract microbiota has been associated with disease in both humans and dogs including chronic diarrhea (Bell et al., 2008; Jia et al., 2010) and inflammatory bowel disease (IBD) (Nobaek et al., 2000; Janeczko et al., 2008; Xenoulis et al., 2008; Suchodolski et al., 2010). With disease, there is usually specific shifts in microbial population or a decrease in overall diversity, making disturbances in the GI tract microbiota a possible early warning sign for disease (Deng and Swanson, 2015).
To date, most of the research investigating the dog microbiota has analyzed fecal samples from healthy (free of disease with ideal body condition score) Beagle dogs in controlled laboratory settings (Vanhoutte et al., 2005; Middelbos et al., 2010; Handl et al., 2011; Hang et al., 2012; Beloshapka et al., 2013; Deng and Swanson, 2015; Panasevich et al., 2015; Herstad et al., 2017). These studies have shown that bacteria dominate the canine gut microbiota accounting for ~99% of total sequences with archaea accounting for the remaining 1% (Middelbos et al., 2010; Swanson et al., 2010; Garcia-Mazcorro et al., 2011; Handl et al., 2011; Hand et al., 2013). The predominant phyla found in the GI tract of healthy dogs are Firmicutes, Bacteroidetes, Proteobacteria, Fusobacteria, and Actinobacteria (Suchodolski et al., 2008; Middelbos et al., 2010; Swanson et al., 2010; Hooda et al., 2012; Deng and Swanson, 2015; Herstad et al., 2017; Li et al., 2017; Schauf et al., 2018). However, bacterial species typically indicated as pathogens such as Clostridium difficile, Clostridium perfringens, Enterococcus spp., E. coli, and Helicobacter are often considered part of a dog's healthy microbiota (Jia et al., 2010; Handl et al., 2011; Goldstein et al., 2012). The fecal microbiome of the dog reflects the high concentrations of protein and fat in their diets (Moon et al., 2018). For reference, a typical canned dog food contains 20–32% fat, 28–50% protein, and 18–57% carbohydrate on a dry matter basis (Case et al., 2011).
Compared to other nutrients, dietary fat and its effect on the microbiota have been underestimated due to the argument that little dietary fat reaches the colon where the highest density of bacteria reside (Cândido et al., 2018). However, Gabert et al. (2011) showed that free fatty acids were being excreted in healthy people. Free fatty acids are known to have potent antimicrobial effects even at small doses (Huang et al., 2011; Cândido et al., 2018). Therefore, the small amounts of fat reaching the colon could interact with the resident microbiota. In addition, a higher fat content will require an increased amount of bile acids for digestion, which are also known to have an antimicrobial effect (Stacey and Webb, 1947). Specific bacteria are known to be involved with the digestion and absorption of dietary fat. Lactobacillus, Bifidobacterium, Enterobacter, Bacteroides, and Clostridium are involved in bile acid metabolism and affect the absorption of dietary fats and lipid-soluble vitamins (Ridlon et al., 2006, Swann et al., 2011). Faecalibacterium prausnitzii and Bifidobacterium are associated with choline metabolism to modulate lipid metabolism and glucose homeostasis (Martin et al., 2010; Wang et al., 2011).
In recent studies, high fat diets are typically associated with a decrease in overall microbial abundance and diversity with a shift from Bacteroidetes to Firmicutes (Hildebrandt et al., 2009; Zhang et al., 2012; Murphy et al., 2015). This shift may lead to increased gut permeability, inflammation, and disease (Murphy et al., 2015). These aforementioned studies have been conducted using murine models with a lack of evidence in canine or other large animal models. Due to the high inclusion rate of dietary fat in pet foods and the demand for less processed diets (decreased carbohydrates), it is important to understand the role of dietary fat on the canine microbiota. Additionally, a recent study by Coelho et al. (2018) using metagenome shotgun sequencing has determined the genetic potential of canine gut microbial communities (consisting of 1,247,405 non-redundant genes) to be more similar to human gut microbial communities, compared to swine and murine microbiota. The authors suggest that canine models may be more accurate in estimating the impact of dietary intervention on human microbial communities.
The study objective within this manuscript was to evaluate the effects of feeding adult dogs increasing levels of fat in low carbohydrate diets on the fecal microbiome. The hypothesis of this study was that microbial shifts would occur based on microbial adaptation to dietary intervention, but dogs would maintain health status due to their ability to efficiently digest fat. Results concerning diet digestibility, fecal characteristics, and blood parameters of the dogs used in this study have been previously published (Kilburn et al., 2020).
Materials and Methods
The protocol for this experiment was reviewed and approved by the Iowa State University Institutional Animal Care and Use Committee (IACUC).
Animals and Housing
Eight female beagles, 1 year of age with an average baseline body weight of 8.57 ± 0.93 kg (mean ± SD) were enrolled in this study. Dogs were spayed prior to the study to prevent any confounding hormonal effects. To ensure that all dogs were healthy before conducting the trial, complete blood count and chemistry panels were performed. In addition, fecal samples were analyzed for parasite presence. All dogs were housed in pairs at the College of Veterinary Medicine at Iowa State University (Ames, IA, USA) in temperature-controlled rooms (20°C) on a 12:12 h light: dark schedule. During feeding and collection periods, dogs were separated by gate closure.
Diets and Feeding
Dietary compositions are presented in Table 1. A commercially manufactured canned canine diet (Supplementary Table 1) was used as a control. Canola oil was then added to the control at 2%, 4%, and 6%, as fed, to create the additional treatments. Treatment diets contained 32% (T1), 37% (T2), 42% (T3), and 47% (T4) total dietary fat on a dry matter basis. The control diet (T1) was selected as it was already higher in fat and lower in carbohydrates compared to other commercially manufactured diets. As canola oil was added, the estimated carbohydrate content (nitrogen-free extract) of the diets decreased. Diets contained an estimated 6.61% (T1), 7.52% (T2), 5.84% (T3), and 3.71% (T4) nitrogen-free extract. Of note, the calculated nitrogen-free extract of T2 was greater than expected.
Dogs were fed twice daily (0800 and 1700 h) to meet their daily energy requirements. Total daily energy requirements were calculated per treatment for each individual dog based on body weight at the beginning of each period. In other words, as dietary fat (energy) of the diet increased less diet was offered. Therefore, as dietary fat increased total carbohydrate consumed decreased. Weight and body condition score (BCS) were recorded weekly. If needed, feed intake was adjusted during the adaption phase to maintain ideal BCS. Water was provided ad libitum throughout the study.
Experimental Design and Sample Collection
Dogs were randomly assigned to one of four dietary treatments in a replicated 4 × 4 Latin Square design. This design allowed each dog to serve as its own control. Each period included a 10-day diet adaption phase followed by a 5-day total collection phase.
During the collection phase, 2 g of fresh feces (defecated within 15 min) were placed into a cryovial tube and immediately stored in −80°C for microbiota analysis for each dog per treatment.
Fecal DNA Extraction
Fecal samples were thawed, and DNA was extracted from ~0.25 g of feces using the Qiagen DNeasy Powerlyzer Powersoil kit (Germantown, MD) following the manufacturer's instructions. Mechanical cell lysis was performed using a Fischer Scientific Beadmill 24. DNA concentrations were determined using a spectrophotometer (ND-100; NanoDrop Technologies, Inc., Rockland, DE) prior to sequencing. DNA concentrations for all samples are shown in Supplementary Table 2.
After extraction, DNA was sent to the Iowa State University DNA facility for paired-end, 16S rRNA gene amplicon sequencing (V4) using the Illumina MiSeq platform. Briefly, the genomic DNA from each sample was amplified using the Platinum™ Taq DNA Polymerase (Thermo Fisher Scientific, Waltham, MA) with one replicate per sample using universal 16S rRNA gene bacterial primers [515F (5′ GTGYCAGCMGCCGCGGTAA-3′; Parada et al., 2016) and 806R (5′-GGACTACNVGGGTWTCTAAT-3′; Apprill et al., 2015)] amplifying the variable region V4. All samples underwent PCR with an initial denaturation step at 94°C for 3 min, followed by 45 s of denaturing at 94°C, 20 s of annealing at 50°C, and 90 s of extension at 72°C. This was repeated for 35 total PCR cycles and finished with a 10 min extension at 72°C. DNA was then purified of primers, nucleotides, enzymes, mineral oil, salts, agarose, ethidium bromide, and other impurities using the QIAquick 96 PCR Purification Kit (Qiagen Sciences Inc, Germantown, MD). PCR bar-coded amplicons were mixed at equal molar ratios and used for Illumina MiSeq paired-end sequencing with 250 bp read length and cluster generation with 10% PhiX control DNA on an Illumina MiSeq platform (Illumina Inc., San Diego, CA).
Sequence analysis was done with Mothur V1.40.5 following the Mothur MiSeq SOP (Kozich et al., 2013). Barcode sequences, primer and low-quality sequences were trimmed using a minimum average quality score of 35, with a sliding window size of 50 bp. Chimeric sequences were removed with the “Chimera.uchime” command. For alignment and taxonomic classification of sequences, the SILVA SSU NR reference database v132 provided by the Mothur website was used. The sequences were then clustered into operational taxonomic units (OTUs) based on 99% 16S rRNA gene similarity (=0.01 distance).
Entire microbial communities of each sample were either rarefied to the lowest sequencing depth per sample (20,900), or non-rarefied prior to assigning Bray-Curtis dissimilarity coefficients to perform statistical comparisons between treatment groups. After dissimilarity coefficients were assigned to each sample, treatment groups were compared using the Adonis (PERMANOVA) command from the vegan package in R (Oksanen et al., 2019).
Both rarefied and non-rarefied microbial communities were visualized by plotting (ggplot2 v2_3.1.1 graphing package in R 3.6.0; Wickham, 2009; R Core Team, 2019) principle coordinate analysis (PCoA) generated with the Phyloseq (v1.28.0, McMurdie and Holmes, 2013) and vegan (v2.5-5) packages using the shared and taxonomy file generated in Mothur. Bray-Curtis dissimilarity measures were used to generate distances between samples for the PCoA plot, then each sample value was plotted.
Canonical analysis of principle coordinates (CAP) (Anderson and Willis, 2003) was conducted to detect any differences in whole microbial communities based on treatment in relation to animal health parameters collected in Kilburn et al. (2020). Once again, Bray-Curtis dissimilarity measures were used to generate distances between samples, which were then constrained based on model effects (treatment, room, and period) and certain animal health measurements (fat digestibility, fecal dry matter, feed intake (as fed), and red blood cell distribution width) collected in the previous study. Animal health measurements used to constrain the data were selected based on (1) significance detected in Kilburn et al. (2020) and (2) correlations detected between these animal parameters. If variables were correlated, a single variable was selected to represent all correlated variables based on study relevance (i.e., fat digestibility representing gross energy digestibility and dry matter digestibility).
To compare alpha diversity between experimental groups, reads were either rarefied to accommodate the sample with the lowest number of reads (20,900 sequences), or non-rarefied similar to entire microbial community comparisons. Measurements of Chao species richness, Shannon diversity, and Simpson evenness were generated within Phyloseq to compare community characteristics between experimental groups. The means of the treatment group alpha diversity measures were compared with ANOVA assuming equal variance.
To create phylum, class, and genus level comparison bar graphs between treatment groups, all sequences agglomerated (tax_glom command in Phyloseq) based on their taxonomic classification assigned via the classify.seqs command using the Silva reference database in Mothur. These sequences were then adjusted to relative abundance values and plotted.
Phyla and individual OTUs were analyzed using the GLIMMIX procedure of SAS (Version 9.4, SAS Inst., Cary, NC) with fixed effects of diet and room, and the random effects of period and animal. A negative binomial was used to determine the distribution with an offset of log library size. P-values were transformed to q-values using false discovery rate (FDR) correction (Storey, 2002). Q-values were used to determine significance (q < 0.05). Orthogonal contrasts were performed on significant phyla and OTUs to determine linear, quadratic, and/or cubic relationships among treatments.
Statistical Analysis of Body Weight and Body Condition Score
Body weight and BCS were analyzed using the MIXED procedure of SAS (Version 9.4, SAS Inst., Cary, NC) with fixed effects of diet and room, and the random effects of period and animal. Initial body weight or initial body condition score were used as a covariate for their respective analysis. Differences between diets were determined using least squared means. A probability of P < 0.05 was considered statistically significant and standard error of the means (SEM) were determined.
The 16S rRNA gene sequences have been submitted to the NCBI Sequence Read Archive SRA and are available under the BioProject ID PRJNA630443.
Body Weight and Body Condition Score
Body weight and BCS are presented in Table 2. Mean body weight (P = 0.199) and BCS (P = 0.907) of dogs were maintained throughout treatments.
Fecal Microbial Communities
Overall, 2,438 OTUs were generated after quality control and removal of OTUs representing <10 sequences. The average number of sequences per samples was 59,783 with a standard deviation of 25,370. 99.9% of the reads were bacterial while only 0.1% were archaeal. From the 2,438 OTUs, 25 phyla were identified with Firmicutes (40%), Bacteroidetes (34%), Fusobacteria (17%), Proteobacteria (7%), and Actinobacteria (1%) being the most abundant. The most abundant phyla, classes, and genera per treatment are presented in Figure 1. The classes Bacteroidia, Clostridia, Fusobacteriia, and Erysipelotrichia accounted for 33%, 30%, 17%, and 6% of total reads, respectively. Additionally, the Fusobacterium genus accounted for 18% of all reads. OTU 1 was classified into the Peptoclostridium genus which accounted for 14% of total reads. Several OTUs were classified within the genera Bacteroides and Alloprevotella, with each accounting for 10% of all reads. In addition, the genus Allobaculum accounted for 3% of the total reads. The assigned classifications of the 50 most abundant OTUs are presented in Supplementary Table 3.
Figure 1. Relative abundance of the 10 most abundant (A) phyla, (B) classes, and (C) genera in dog fecal samples per dietary treatment.
When comparing entire bacterial communities of treatment groups using PERMANOVA, no significant differences were observed in either the rarefied or non-rarefied data (P = 0.735 and P = 0.834, respectively, Supplementary Table 4). This result was supported by the lack of apparent clustering of samples based on treatment types seen in the unconstrained PCoA (Figure 2). Additionally, although we selected animal health measurements that relate to GI microbiota, no clear trends were detected across treatment when constrained by these parameters (Figure 3). Finally, no significant treatment differences were detected across treatment for alpha diversity estimators (Table 3) either.
Figure 2. Beta diversity of dog fecal microbial communities revealed by a PCoA plot based on Bray-Curtis dissimilarities of the overall composition of microbial communities among dietary treatment groups. (A) Rarefied to 20,900 sequences and (B) non-rarefied. Each point corresponds to a community from a single dog. Colors represent each treatment. A summary of PERMANOVA results (P-value and R2) are displayed with the boxes. Additional PERMANOVA information can be found in Supplementary Table 4.
Figure 3. Beta diversity of dog fecal microbial communities revealed by a CAP plot based on Bray-Curtis dissimilarities of the overall composition of microbial communities among dietary treatment groups adjusted by animal health parameters collected by Kilburn et al. (2020). (A) Rarefied to 20,900 sequences and (B) non-rarefied. Each point corresponds to a community from a single dog. Colors represent each treatment.
However, significant differences (q < 0.05) in the relative abundance of certain phyla were detected between treatment groups (Table 4). Treatment differences were reported for phyla Tenericutes (q = 0.0052), Spirochaetes (q < 0.001), Euryarchaeota (q = 0.0007), Fibrobacteres (q < 0.001), Kiritimatiellaeota (q < 0.001), Deferribacteres (q < 0.001), and Planctomycetes (q < 0.001). Tenericutes, Spirochaetes, Fibrobacteres, and Planctomycetes increased linearly with the increase in dietary fat (P < 0.05).
Similarly, when comparing the abundance of the 100 most abundant individual OTUs between treatment levels, 36 OTUs showed significant (q < 0.05) differences (Table 5). For example, treatment differences were reported among the genera Catenibacterium (q < 0.037), Paeniclostridium (q = 0.010), Romboutsia (q = 0.024), Blautia (q = 0.024), and Lactobacillus (q < 0.001). A visual representation of significant shifts in OTU abundance per treatment among the 50 most abundant OTUs is presented in Figure 4.
Table 5. Relative abundance of significanta OTUs among treatments out of 100 most abundant OTUs of dog fecal samples.
Figure 4. Bar charts presenting a visual representation of the shifts in relative abundance of significant (q < 0.05) OTUs per dietary treatment among the 50 most abundant OTUs of dog fecal samples. Error bars represent the standard error among samples (see Table 4 for details).
OTUs 12, 15, 31, 32, 35, 36, 46, 48, 57, 62, 72, 90, and 99 significantly (P < 0.05) increased in abundance from T1 to T4 (Table 5). OTUs 12 and 32 (P < 0.001) were assigned to the genus Allobaculum, of the Firmicutes phylum. Within the Proteobacteria phylum, OTU 15 (P < 0.001), OTU 31 (P < 0.001), OTU 48 (P = 0.001), and OTU 62 (P < 0.001) were assigned to Sphingomonadaceae unclassified, Parasutterela, Histophilus, and Succinivibrionaceae_UCG-001, respectively. OTU 35 (P = 0.005) was assigned to the Paeniclostridium genus of the Firmicutes phylum.
Other OTUs showed a significant (P < 0.05) linear decrease in abundance with increasing fat content in the diet, these included OTUs 2, 3, 9, 17, 21, 29, 37, 41, 42, 43, 44, 45, 58, 66, 86, and 96 (Table 5). OTU 2 was assigned to Prevotellaceae Ga6A1 group (P = 0.029). OTU 3 was assigned to the genus Fusobacterium (P < 0.001) of the Fusobacteria phylum. OTU 9 was assigned to Catenibacterium (P = 0.010). OTUs 17, 21, and 96 were assigned to the genus Alloprevotella (P = 0.017, P < 0.001, P = 0.005, respectively). OTU 29 (P = 0.001) and OTU 86 (P < 0.001) were assigned to the Lachnospiraceae family. OTU 37 and OTU 45 were both assigned to the genus Bacteroides (P = 0.002, P = 0.006). Most of the genera which significantly decreased in abundance belong to the Bacteroidetes or Firmicutes phylum. OTUs 2, 17, 21, 37, 45, and 96 were all classified to genera within the Bacteroidetes phylum while OTUs 9, 29, 43, 44, 58, 66, and 86 all belong to the Firmicutes phylum.
Body Weight and Body Condition Score
It is important that dogs maintained ideal body weight and BCS due to known changes in the microbiota with obesity which may confound diet effects (de La Serre et al., 2010; Cândido et al., 2018). In addition, high fat diets and obesity have been shown to have similar effects on the microbiota making it difficult to determine which caused the microbiota to change if modeled together (Ley et al., 2005; Turnbaugh et al., 2006, 2008; Cândido et al., 2018). The maintenance of ideal BCS allowed this study to measure high fat diets independently of obesity.
Fecal Microbial Communities
We chose to study fecal microbiota as a proxy for gut microbiota as fecal samples make non-invasive periodic measurements possible. One major limitation of using fecal samples is the fact that fecal samples are more representative of the digesta in the lower gut and might not adequately represent important host-microbial interactions specific to different regions of the GI tract (Leite et al., 2020). In addition, 16S rRNA gene amplicon datasets offer only limited information about functional roles of members of these microbial communities. In this manuscript, any connections between the microbial community and a host phenotype were made using previously published work regarding assigned taxonomy.
Overall, on a whole community level, the high-fat diet did not result in significant shifts in the fecal microbial communities regarding alpha and beta diversity. Shannon diversity indices in this study were greater than those previously reported by Handl et al. (2013) but lower than those reported by Schauf et al. (2018) in dogs. Similar to the study reported here, Schauf et al. (2018) found no difference in bacterial species richness or Shannon diversity comparing a low-fat, high-starch diet to a high-fat, low-starch diet in dogs. Additionally, Coelho et al. (2018) observed a more stable microbiota in non-obese dogs similar to the animals used in this experiment. Other studies using mice have reported a decrease in GI tract microbial abundance when fed high fat diets (Hildebrandt et al., 2009; Zhang et al., 2012). Martinez et al. (2009) reported a decrease in species richness in hamsters fed grain sorghum lipid extract. The decrease in species richness may be caused by the antimicrobial effect of fatty acids and/or bile acids (Stacey and Webb, 1947; Huang et al., 2011; Cândido et al., 2018). To summarize, studies that have fed high fat diets to mice have indicated a decrease in microbial richness while those in dogs have reported no difference with increased dietary fat levels. The varying results in species richness may be due to the differences in physiology and typical diet composition of the animal species used in the studies.
Previous studies have shown a range of abundance for the dominant phyla in healthy dogs with 14–48% Firmicutes, 12–38% Bacteroidetes, 7–44% Fusobacteria, 5–23% Proteobacteria, and 0.8–1.4% Actinobacteria (Suchodolski et al., 2008; Middelbos et al., 2010; Swanson et al., 2010; Herstad et al., 2017; Li et al., 2017; Coelho et al., 2018; Salas-Mani et al., 2018). The total abundance of each phylum found in this study fall within those ranges indicating normal values. A common metric used in obesity studies from humans, mice, and canines is to report significant differences detected in Bacteroidetes and Firmicutes. Coelho et al. (2018) reported an increase in the Firmicutes:Bacteroides ratio with the switch to a diet with increased protein levels and reduced carbohydrate levels. In addition, high fat diets have been reported to decrease Bacteroidetes and increase Firmicutes in mice (Hildebrandt et al., 2009; Zhang et al., 2012; Murphy et al., 2015) possibly leading to dysbiosis. No differences in Firmicutes abundance or Bacteroidetes abundance were seen with the shift to the higher fat, lower carbohydrate diet (T4) in this study (Table 4).
The following section of the manuscript will discuss potential contributions of specific OTUs that demonstrated a linear trend (either increasing or decreasing) with treatment. The data reported in this manuscript is gathered from additional peer-review publications, and a summary of the points discussed can be found in Supplementary Table 5.
Bacterial Populations That Increased With Increasing Levels of Dietary Fat may Utilize Fat or Bile as Key Metabolic Substrate
OTUs classified as Allobaculum showed differing trends with treatment. The abundance of OTUs 12 and 32 increased with increasing fat levels, whereas the abundance of OTU 58 decreased. Using a 99% OTU clustering threshold, OTUs 12, 32, and 58 could be considered different species which would explain the differences in behavior. In other words, even though these OTUs are classified within the Allobaculum genus, they may contain different genetic potential, and thus different metabolic capabilities. Allobaculum has been suggested to have beneficial effects and contribute to mucus formation (Everard et al., 2014). The presence of Allobaculum and its effect on the host differ among studies. Martinez et al. (2009) reported an increase in Allobaculum in hamsters fed grain sorghum lipid extract. In contrast, Everard et al. (2014) and Ravussin et al. (2012) found an increase in Allobaculum in lower fat diets compared to high fat diets. Jakobsson et al. (2015) questioned the beneficial role of Allobaculum by showing that mice had increased mucus permeability with increased abundance. Different physiological properties of closely related phylotypes have been revealed in previous studies (Berry et al., 2012) and highlight the importance of OTU-level analyses of microbiome data.
Four significant OTUs (15, 31, 48, and 62) were classified into genera within the Proteobacteria phylum. Previous studies have reported an increase in Proteobacteria with increased consumption of a more natural diet (higher in fat) compared to kibble diets in dogs (Herstad et al., 2017; Sandri et al., 2017). In addition, Parasutterella (OTU 31) has a potential role in bile acid maintenance and cholesterol metabolism (Ju et al., 2019). The increase in Parasutterella may be explained by the higher fat content which will require an increased amount of bile acid secretion for digestion (Di Ciaula et al., 2017).
OTUs 46 and 57 were classified as uncultured or unclassified Erysipelotrichaceae. Another study reported a decreased abundance of unclassified members of the Erysipelotrichaceae family in hamsters consuming diets containing grain sorghum lipid extract (Martinez et al., 2009). The differing results may be due to the difference in animal species consuming such diets as the GI tract of a hamster is unlike that of a dog.
Bifidobacterium (OTU 90) increased in abundance with increased dietary fat level and has been described as beneficial due to its ability to reduce intestinal endotoxins and improve mucosal barrier function (Griffiths et al., 2004; Wang et al., 2006). In addition, the genus Bifidobacterium has bile salt hydrolase (BSH) activity (Ridlon et al., 2006), which might be important when growing on high-fat diets. Martinez et al. (2009) reported an increase in Bifidobacterium in hamsters fed grain sorghum lipid extract. However, a study feeding mice a high fat diet found a decrease in bifidobacteria (Cani et al., 2007).
The genus Lactobacillus (OTU 99) is known to metabolize carbohydrates resulting in lactic acid production (Walter, 2008). The increase in Lactobacillus with increase in dietary fat is interesting due its known role in carbohydrate metabolism (Walter, 2008). Related to the study reported in this manuscript, Salas-Mani et al. (2018) and Coelho et al. (2018) reported a decrease in canine GI tract Lactobacillus abundance when fat levels were decreased in canine diets. Lactobacillus has been reported to have BSH activity (Ridlon et al., 2006). Therefore, the increase in this genus may be due to the increased need for BSH activity for fat digestion and not its proposed role in carbohydrate metabolism. In addition, an increase in Lactobacillus is considered beneficial as some Lactobacillus species are commonly used as probiotics due to their health-promoting properties. For example, this genus is thought to modulate the immune system and protect the epithelial barrier in the gut (Lebeer et al., 2008).
Bacterial Populations That Decreased With Increasing Levels of Dietary Fat may Prefer Other Nutrients as Key Metabolic Substrate
Fusobacterium (OTU 3) is reported to utilize amino acids and produce butyrate (Barcenilla et al., 2000; Butowski et al., 2019). In previous studies, Fusobacterium has been linked to high protein diets in dogs (Beloshapka et al., 2013; Bermingham et al., 2017) as well as a general high level of presence within canine fecal microbial communities (Coelho et al., 2018; Wernimont et al., 2020).
Catenibacterium (OTU 9), Anaerobiospirillum (OTU 41), Romboutsia (OTU 43), and Blautia (OTU 44) can utilize many types of carbohydrates to produce succinic acid, acetic acid, lactic acid, butyric acid, iso-butyric acid, and ethanol depending on the genera (Davis et al., 1976; Kageyama and Benno, 2000; Liu et al., 2008; Gerritsen et al., 2014). In addition, Alloprevotella (OTU 17, OTU 21, and OTU 96) is reported to be a saccharolytic bacteria (Qu et al., 2017). A trend was discovered by Salas-Mani et al. (2018), who reported a decrease in Blautia species with the switch from a high fat and high carbohydrate to a low fat and high protein diet in canines, but this may have been driven by the decrease in carbohydrates as opposed to fat levels. Pilla and Suchodolski (2020) point out that a decrease in Blautia species may be associated with both acute diarrhea and canine IBD. Neither of these symptoms were detected in this trial, however. Yan et al. (2013) suggested an increase in Catenibacterium led to increased SCFA production. In this study, the decreased Catenibacterium did not affect short chain fatty acid production among treatments (Kilburn et al., 2020).
The genus Bacteroides (OTU 37 and OTU 45) is involved in the fermentation of indigestible carbohydrates (Handl et al., 2013), which were minimal in dietary treatments in this study. Bacteroides has been reported to increase in humans consuming a Western diet, which is high in fat and sugar (De Filippo et al., 2010; David et al., 2014). Additionally, the Bacteroides genus has been reported by Wu et al. (2011), Coelho et al. (2018), and Salas-Mani et al. (2018) to be positively associated with higher levels of dietary protein. The decrease in OTUs 37 and 45 in response to the increase in dietary fat within this study may indicate that Bacteroides do indeed rely mainly on protein-rich substrates.
Lachnoclostridium has been reported to produce secondary bile acids through bile acid dihydroxylation activity (Ridlon et al., 2015). In the current study, the decrease in this genus is interesting due to the thought that an increase in dietary fat would require greater production of bile acids for digestion. Consequently, an increase in this genus would rather be expected.
Potential Use of Canine Models for Human Disease and Nutrition
As noted in the introduction, studying the dog microbiota and its shifts due to dietary intervention or disease is becoming increasingly relevant for understanding both pet and human health. It is only recently that studies have reported more so on canine microbial communities and their genetic potential rather than that of mice and humans. Additionally, Coelho et al. (2018) reported a high similarity between the genetic catalog of canine microbial communities to that of human microbial communities (63% genes mapped from the dog microbial communities to human gene catalog, as opposed to 32.9% for swine and 19.9% for murine microbial communities). More generally, the high similarity of the dog and human GI tract microbiota suggests possible benefits from using canines as models for human disease or dietary intervention.
In conclusion, the increase in dietary fat and subsequent decrease in carbohydrate levels did not impact the overall microbial diversity in dogs fed dietary treatment. However, the microbiota did shift based on available diet substrate. Even with this microbial shift, dogs remained healthy during the time of the study. This may indicate that the dog microbiota can adapt to high fat diets without creating a dysbiosis. Further research is needed to analyze the functional characteristics of these changes in microbial communities from dogs fed similar dietary treatments.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, PRJNA630443.
The animal study was reviewed and approved by Iowa State University IACUC.
MR and LKi conceived and designed the study. LKi conducted the animal trial and sampling. LKo prepared the samples and analyzed the sequence data. LKo and SS-E assisted in interpretation of results. NS contributed to statistical analysis. All authors have reviewed and approved the final version of the manuscript.
Conflict of Interest
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.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.564160/full#supplementary-material
Anderson, M. J., and Willis, T. J. (2003). Canonical analysis of principal coordinates: a useful method of constrained ordination for ecology. Ecology 84, 511–525. doi: 10.1890/0012-9658(2003)084[0511:CAOPCA]2.0.CO;2
Apprill, A., Mcnally, S., Parsons, R., and Weber, L. (2015). Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat. Microbiol. Ecol. 75, 129–137. doi: 10.3354/ame01753
Backhed, F., Ley, R. E., Sonnenburg, J. L., Peterson, D. A., and Gordon, J. I. (2005). Host-bacterial mutualism in the human intestine. Science 307, 1915–1919. doi: 10.1126/science.1104816
Barcenilla, A., Pryde, S. E., Martin, J. C., Duncan, S. H., Stewart, C. S., Henderson, C., et al. (2000). Phylogenetic relationships of butyrate-producing bacteria from the human gut. Appl. Environ. Microbiol. 66, 1654–1661. doi: 10.1128/AEM.66.4.1654-1661.2000
Bell, J. A., Kopper, J. J., Turnbull, J. A., Barbu, N. I., Murphy, A. J., and Mansfield, L. S. (2008). Ecological characterization of the colonic microbiota of normal and diarrheic dogs. Interdiscip. Perspect. Infect. Dis. 2008, 1–17. doi: 10.1155/2008/149694
Beloshapka, A. N., Dowd, S. E., Suchodolski, J. S., Steiner, J. M., Duclos, L., and Swanson, K. S. (2013). Fecal microbial communities of healthy adult dogs fed raw meat-based diets with or without inulin or yeast cell wall extracts as assessed by 454 Pyrosequencing. FEMS Microbiol. Ecol. 84, 532–541. doi: 10.1111/1574-6941.12081
Bergman, E. N. (1990). Energy contributions of volatile fatty acids from the gastrointestinal tract in various species. Physiol. Rev. 70, 567–590. doi: 10.1152/physrev.19126.96.36.1997
Bermingham, E. N., Maclean, P., Thomas, D. G., Cave, N. J., and Young, W. (2017). Key bacterial families (Clostridiaceae, Erysipelotrichaceae, and Bacteroidaceae) are related to the digestion of protein and energy in dogs. PeerJ. 5:e3019. doi: 10.7717/peerj.3019
Berry, D., Schwab, C., Milinovich, G., Reichert, J., Ben Mahfoudh, K., Decker, T., et al. (2012). Phylotype-level 16S rRNA analysis reveals new bacterial indicators of health state in acute murine colitis. ISME J. 6, 2091–2106. doi: 10.1038/ismej.2012.39
Butowski, C. F., Thomas, D. G., Young, W., Cave, N. J., McKenzie, C. M., Rosendale, D. I., et al. (2019). Addition of plant dietary fibre to a raw red meat high protein, high fat diet, alters the faecal bacteriome, and organic acid profiles of the domestic cat (Felis catus). PLoS ONE 14:e216072. doi: 10.1371/journal.pone.0216072
Cândido, F. G., Valente, F. X., Grześkowiak, Ł. M., Moreira, A. P. B., Rocha, D. M. U. P., Alfenas, R., et al. (2018). Impact of dietary fat on gut microbiota and low-grade systemic inflammation: mechanisms and clinical implications on obesity. Int. J. Food Sci. Nutr. 69, 125–143. doi: 10.1080/09637486.2017.1343286
Cani, P. D., Amar, J., Iglesias, M. A., Poggi, M., Knauf, C., Bastelica, D., et al. (2007). Metabolic endotoxemia initiates obesity and insulin resistance. Diabetes 56, 1761–1772. doi: 10.2337/db06-1491
Case, L. P., Daristotle, L., Hayek, M. G., and Raasch, M. F. (2011). Canine and Feline Nutrition: A Resource for Companion Animal Professionals, 3rd Edn. Maryland Heights, MO: Mosby Elsevier.
Coelho, L. P., Kultima, J. R., Costea, P. I., Fournier, C., Pan, Y., Czarnecki-Maulden, G., et al. (2018). Similarity of the dog and human gut microbiomes in gene content and response to diet. Microbiome 6, 1–11. doi: 10.1186/s40168-018-0450-3
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
Davis, C. P., Cleven, D., Brown, J., and Balish, E. (1976). Anaerobiospirillum, a new genus of spiral-shaped bacteria. Int. J. Syst. Bacteriol. 4, 498–504. doi: 10.1099/00207713-26-4-498
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 La Serre, C. B., Ellis, C. L., Lee, J., Hartman, A. L., Rutledge, J. C., and Raybould, H. E. (2010). Propensity to high-fat diet-induced obesity in rats is associated with changes in the gut microbiota and gut inflammation. Am. J. Physiol. Liver Physiol. 299, 440–448. doi: 10.1152/ajpgi.00098.2010
Deng, P., and Swanson, K. S. (2015). Gut microbiota of humans, dogs, and cats: current knowledge and future opportunities and challenges. Br. J. Nutr. 113, S6–S17. doi: 10.1017/S0007114514002943
Di Ciaula, A., Garruti, G., Baccetto, R. L., Molina-Molina, E., Bonfrate, L., Wang, D. Q. H., et al. (2017). Bile acid physiology. Ann. Hepatol. 16, s4–s14. doi: 10.5604/01.3001.0010.5493
Eckburg, P. B., Bik, E. M., Bernstein, C. N., Purdom, E., Dethlefsen, L., Sargent, M., et al. (2005). Diversity of the human intestinal microbial flora. Science 308, 1635–1638. doi: 10.1126/science.1110591
Everard, A., Lazarevic, V., Gaia, N., Johansson, M., Stahlman, M., Backhed, F., et al. (2014). Microbiome of prebiotic-treated mice reveals novel targets involved in host response during obesity. ISME 8, 2116–2130. doi: 10.1038/ismej.2014.45
Gabert, L., Vors, C., Louche-Pélissier, C., Sauvinet, V., Lambert-Porcheron, S., Drai, J., et al. (2011). 13C tracer recovery in human stools after digestion of a fat-rich meal labelled with [1, 1, 1-13 C3] tripalmitin and [1, 1, 1-13 C3] triolein. Rapid Commun. Mass Spectrom. 25, 2697–2703. doi: 10.1002/rcm.5067
Garcia-Mazcorro, J. F., Lanerie, D. J., Dowd, S. E., Paddock, C. G., Grutzner, N., Steiner, J. M., et al. (2011). Effect of a multi-species synbiotic formulation on fecal bacterial microbiota of healthy cats and dogs as evaluated by pyrosequencing. FEMS Microbiol. Ecol. 78, 542–554. doi: 10.1111/j.1574-6941.2011.01185.x
Gerritsen, J., Fuentes, S., Grievink, W., van Niftrik, L., Tindall, B. J., Timmerman, H. M., et al. (2014). Characterization of Romboutsia ilealis gen. nov., sp. nov., isolated from the gastro-intestinal tract of a rat, and proposal for the reclassification of five closely related members of the genus Clostridium into the genera Romboutsia gen. nov., Intestinibacter gen. nov., Terrisporobacter gen. nov., and Asaccharospora gen. nov. Int. J. Syst. Evol. Microbiol. 64, 1600–1616. doi: 10.1099/ijs.0.059543-0
Goldstein, M. R., Kruth, S. A., Bersenas, A. M. E., Holowaychuk, M. K., and Weese, J. S. (2012). Detection and characterization of Clostridium perfringens in the feces of healthy and diarrheic dogs. Can. J. Vet. Res. 76, 161–165.
Griffiths, E. A., Duffy, L. C., Schanbacher, F. L., Qiao, H., Dryja, D., Leavens, A., et al. (2004). In vivo effects of bifidobacteria and lactoferrin on gut endotoxin concentration and mucosal immunity in Balb/c mice. Dig. Dis. Sci. 49, 579–589. doi: 10.1023/B:DDAS.0000026302.92898.ae
Hand, D., Wallis, C., Colyer, A., and Penn, C. W. (2013). Pyrosequencing the canine faecal microbiota: breadth and depth of biodiversity. PLoS ONE 8:53115. doi: 10.1371/journal.pone.0053115
Handl, S., Dowd, S. E., Garcia-Mazcorro, J. F., Steiner, J. M., and Suchodolski, J. S. (2011). Massive parallel 16S rRNA gene pyrosequencing reveals highly diverse fecal bacterial and fungal communities in healthy dogs and cats. FEMS Microbiol. Ecol. 76, 301–310. doi: 10.1111/j.1574-6941.2011.01058.x
Handl, S., German, A. J., Holden, S. L., Dowd, S. E., Steiner, J. M., Heilmann, R. M., et al. (2013). Faecal microbiota in lean and obese dogs. FEMS Microbiol. Ecol. 84, 332–343. doi: 10.1111/1574-6941.12067
Hang, I., Rinttila, T., Zentek, J., Kettunen, A., Alaja, S., Apajalahti, J., et al. (2012). Effect of high contents of dietary animal-derived protein or carbohydrates on canine faecal microbiota. BMC Vet. Res. 8:90. doi: 10.1186/1746-6148-8-90
Herstad, K. M. V., Gajardo, K., Bakke, A. M., Moe, L., Ludvigsen, J., Rudi, K., et al. (2017). A diet change from dry food to beef induces reversible changes on the faecal microbiota in healthy, adult client-owned dogs. BMC Vet. Res. 13, 1–13. doi: 10.1186/s12917-017-1073-9
Hildebrandt, M. A., Hoffmann, C., Sherrill-Mix, S. A., Keilbaugh, S. A., Hamady, M., Chen, Y. Y. U., et al. (2009). High-fat diet determines the composition of the murine gut microbiome independently of obesity. Gastroenterology 137, 1716–1724.e2. doi: 10.1053/j.gastro.2009.08.042
Hooda, S., Minamoto, Y., Suchodolski, J. S., and Swanson, K. S. (2012). Current state of knowledge: the canine gastrointestinal microbiome. Anim. Heal. Res. Rev. 13, 78–88. doi: 10.1017/S1466252312000059
Hooper, L. V., Wong, M. H., Thelin, A., Hansson, L., Falk, P. G., and Gordon, J. I. (2001). Molecular analysis of commensal host-microbial relationships in the intestine. Science 291, 881–884. doi: 10.1126/science.291.5505.881
Huang, C. B., George, B., and Ebersole, J. L. (2011). Antimicrobial activity of n-6, n-7, and n-9 fatty acids and their esters for oral microorganisms. Arch Oral Biol. 55, 555–560. doi: 10.1016/j.archoralbio.2010.05.009
Jakobsson, H. E., Rodríguez-Piñeiro, A. M., Schütte, A., Ermund, A., Boysen, P., Bemark, M., et al. (2015). The composition of the gut microbiota shapes the colon mucus barrier. EMBO Rep. 16, 164–177. doi: 10.15252/embr.201439263
Janeczko, S., Atwater, D., Bogel, E., Greiter-wilke, A., Gerold, A., Baumgart, M., et al. (2008). The relationship of mucosal bacteria to duodenal histopathology, cytokine mRNA, and clinical disease activity in cats with inflammatory bowel disease. Vet. Microbiol. 128, 178–193. doi: 10.1016/j.vetmic.2007.10.014
Jia, J., Frantz, N., Khoo, C., Gibson, G. R., Rastall, R. A., and Mccartney, A. L. (2010). Investigation of the faecal microbiota associated with canine chronic diarrhoea. FEMS Microbiol. Ecol. 71, 304–312. doi: 10.1111/j.1574-6941.2009.00812.x
Ju, T., Kong, J. Y., Stothard, P., and Willing, B. P. (2019). Defining the role of Parasutterella, a previously uncharacterized member of the core microbiota. ISME 13, 1520–1534. doi: 10.1038/s41396-019-0364-5
Kageyama, A., and Benno, Y. (2000). Catenibacterium mitsuokai gen. nov., sp. nov., a gram-positive anaerobic bacterium isolated from human faeces. Int. J. Syst. Evol. Microbiol. 50, 1595–1599. doi: 10.1099/00207713-50-4-1595
Kilburn, L. R., Allenspach, K., Jergens, A. E., Bourgois-Mochel, A., Mochel, J. P., and Serao, M. C. R. (2020). Apparent total tract digestibility, fecal characteristics, and blood parameters of healthy adult dogs fed high-fat diets. J. Anim. Sci. 98, 1–8. doi: 10.1093/jas/skaa043
Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K., and Schloss, P. D. (2013). Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the miseq illumina sequencing platform. Appl. Environ. Microbiol. 79, 5112–5120. doi: 10.1128/AEM.01043-13
Lebeer, S., Vanderleyden, J., and De Keersmaecker, S. C. J. (2008). Genes and molecules of lactobacilli supporting probiotic action. Microbiol. Mol. Biol. Rev. 72, 728–764. doi: 10.1128/MMBR.00017-08
Leite, G. G. S., Weitsman, S., Parodi, G., Celly, S., Sedighi, R., Sanchez, M., et al. (2020). Mapping the segmental microbiomes in the human small bowel in comparison with stool: a REIMAGINE study. Dig. Dis. Sci. 65, 2595–2604. doi: 10.1007/s10620-020-06173-x
Ley, R. E., Backhed, F., Turnbaugh, P., Lozupone, C. A., Knight, R. D., and Gordon, J. I. (2005). Obesity alters gut microbial ecology. PNAS 102, 11070–11075. doi: 10.1073/pnas.0504978102
Ley, R. E., Hamady, M., Lozupone, C., Turnbaugh, P., Ramey, R. R., Bircher, J. S., et al. (2008). Evolution of mammals and their gut microbes. Science 320, 1647–1651. doi: 10.1126/science.1155725
Li, Q., Lauber, C. L., Czarnecki-Maulden, G., Pan, Y., and Hannah, S. S. (2017). Effects of the dietary protein and carbohydrate ratio on gut microbiomes in dogs of different body conditions. MBio 8, 1–14. doi: 10.1128/mBio.01703-16
Liu, C., Finegold, S. M., Song, Y., and Lawson, P. A. (2008). Reclassification of Clostridium coccoides, Ruminococcus hansenii, Ruminococcus hydrogenotrophicus, Ruminococcus luti, Ruminococcus productus, and Ruminococcus schinkii as Blautia coccoides gen. nov., comb. nov., Blautia hansenii comb. nov., Blautia hydroge. Int. J. Syst. Evol. Microbiol. 58, 1896–1902. doi: 10.1099/ijs.0.65208-0
Mackie, R. I., Sghir, A., and Gaskins, H. R. (1999). Developmental microbial ecology of the neonatal gastrointestinal. Am. J. Clin. Nutr. 69, 1035S−1045S. doi: 10.1093/ajcn/69.5.1035s
Martin, F.-P. J., Sprenger, N., Montoliu, I., Rezzi, S., Kochhar, S., and Nicholson, J. K. (2010). Dietary modulation of gut functional ecology studied by fecal metabonomics. J. Proteome Res. 9, 5284–5295. doi: 10.1021/pr100554m
Martinez, I., Wallace, G., Zhang, C., Legge, R., Benson, A. K., Carr, T. P., et al. (2009). Diet-Induced metabolic improvements in a hamster model of hypercholesterolemia are strongly linked to alterations of the gut microbiota. Appl. Environ. Microbiol. 75, 4175–4184. doi: 10.1128/AEM.00380-09
McKenney, P. T., and Pamer, E. G. (2015). From hype to hope: the gut microbiota in enteric infectious disease. Cell 163, 1326–1332. doi: 10.1016/j.cell.2015.11.032
McMurdie, P. J., and Holmes, S. (2013). Phyloseq: an r package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 8:e61217. doi: 10.1371/journal.pone.0061217
Middelbos, I. S., Vester Boler, B. M., Qu, A., White, B. A., Swanson, K. S., and Fahey, G. C. J. (2010). Phylogenetic characterization of fecal microbial communities of dogs fed diets with or without supplemental dietary fiber using 454 pyrosequencing. PLoS ONE 5:e09768. doi: 10.1371/journal.pone.0009768
Moon, C. D., Young, W., Maclean, P. H., Cookson, A. L., and Bermingham, E. N. (2018). Metagenomic insights into the roles of Proteobacteria in the gastrointestinal microbiomes of healthy dogs and cats. Microbiol. Open 7, 1–20. doi: 10.1002/mbo3.677
Murphy, A. E., Velazquez, K. T., and Herbert, K. M. (2015). Influence of high-fat-diet on gut microbiota: a driving force for chronic disease risk. Curr. Opin. Clin. Nutr. Metab. Care 18, 515–520. doi: 10.1097/MCO.0000000000000209
Nobaek, S., Johansson, M., Molin, G., Ahrne, S., and Jeppsson, B. (2000). Alteration of intestinal microflora is associated with reduction in abdominal bloating and pain in patients with irritable bowel syndrome. Am. J. Gastroenterol. 95, 1231–1238. doi: 10.1111/j.1572-0241.2000.02015.x
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2019). Vegan: Community Ecology Package. R Package Version 2.5–5. Available online at: https://CRAN.R-project.org/package=vegan
Panasevich, M. R., Kerr, K. R., Dilger, R. N., Fahey, G. C. J., Guerin-Deremaux, L., Lynch, G. L., et al. (2015). Modulation of the faecal microbiome of healthy adult dogs by inclusion of potato fibre in the diet. Br. J. Nutr. 113, 125–133. doi: 10.1017/S0007114514003274
Parada, A. E., Needham, D. M., and Fuhrman, J. A. (2016). Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series, and global field samples. Environ. Microbiol. 18, 1403–1414. doi: 10.1111/1462-2920.13023
Pilla, R., and Suchodolski, J. S. (2020). The role of the canine gut microbiome and metabolome in health and gastrointestinal disease. Front. Vet. Sci. 6:498. doi: 10.3389/fvets.2019.00498
Qu, W., Yuan, X., Zhao, J., Zhang, Y., Hu, J., Wang, J., et al. (2017). Dietary advanced glycation end products modify gut microbial composition and partially increase colon permeability in rats. Mol. Nutr. Food Res. 61, 1–14. doi: 10.1002/mnfr.201700118
R Core Team (2019). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Ravussin, Y., Koren, O., Spor, A., LeDuc, C., Gutman, R., Stombaugh, J., et al. (2012). Responses of gut microbiota to diet composition and weight loss in lean and obese mice. Obesity 20, 738–747. doi: 10.1038/oby.2011.111
Ridlon, J. M., Kang, D.-J., and Hylemon, P. B. (2006). Bile salt biotransformations by human intestinal bacteria. J. Lipid Res. 47, 241–259. doi: 10.1194/jlr.R500013-JLR200
Ridlon, J. M., Kang, D. J., Hylemon, P. B., and Bajaj, J. S. (2015). Gut microbiota, cirrhosis, and alcohol regulate bile acid metabolism in the gut. Dig. Dis. 33, 338–345. doi: 10.1159/000371678
Rooks, M. G., and Garrett, W. S. (2016). Gut microbiota, metabolites, and host immunity. Nat. Rev. Immunol. 16, 341–352. doi: 10.1038/nri.2016.42
Salas-Mani, A., Jeusette, I., Castillo, I., Manuelian, C. L., Lionnet, C., Iraculis, N., et al. (2018). Fecal microbiota composition changes after a BW loss diet in Beagle dogs. J. Anim. Sci. 96, 3102–3111. doi: 10.1093/jas/sky193
Sandri, M., Dal Monego, S., Conte, G., Sgorlon, S., and Stefanon, B. (2017). Raw meat based diet influences faecal microbiome and end products of fermentation in healthy dogs. BMC Vet. Res. 13, 1–11. doi: 10.1186/s12917-017-0981-z
Schauf, S., de la Fuente, G., Newbold, C. J., Salas-Mani, A., Torre, C., Abecia, L., et al. (2018). Effect of dietary fat to starch content on fecal microbiota composition and activity. J. Anim. Sci. 96, 3684–3698. doi: 10.1093/jas/sky264
Spor, A., Koren, O., and Ley, R. (2011). Unravelling the effects of the environment and host genotype on the gut microbiome. Nat. Rev. 9, 279–290. doi: 10.1038/nrmicro2540
Stacey, M., and Webb, M. (1947). Studies on the antibacterial properties of the bile acids and some compounds derived from cholanic acid. Proc. R. Soc. Med. 134, 523–537. doi: 10.1098/rspb.1947.0029
Storey, J. D. (2002). A direct approach to false discovery rates. J. R. Statist. Soc. B. 64, 479–498. doi: 10.1111/1467-9868.00346
Suchodolski, J. S., Camacho, J., and Steiner, J. M. (2008). Analysis of bacterial diversity in the canine duodenum, jejunum, ileum, and colon by comparative 16S rRNA gene analysis. FEMS Microbiol. Ecol. 66, 567–578. doi: 10.1111/j.1574-6941.2008.00521.x
Suchodolski, J. S., Xenoulis, P. G., Paddock, G. C., Steiner, J. M., and Jergens, A. E. (2010). Molecular analysis of the bacterial microbiota in duodenal biopsies from dogs with idiopathic inflammatory bowel disease. Vet. Microbiol. 142, 394–400. doi: 10.1016/j.vetmic.2009.11.002
Sunvold, G. D., Fahey, J. G. C., Merchen, N. R., Titgemeyer, E. C., Bourquin, L. D., Bauer, L. L., et al. (1995). Dietary fiber for dogs : IV . In Vitro fermentation of selected fiber sources by dog fecal inoculum and in vivo digestion and metabolism of fiber-supplemented diets. J. Anim. Sci. 73, 1099–1109. doi: 10.2527/1995.7341099x
Swann, J. R., Want, E. J., Geier, F. M., Spagou, K., Wilson, I. D., Sidaway, J. E., et al. (2011). Systemic gut microbial modulation of bile acid metabolism in host tissue compartments. PNAS 108, 4523–4530. doi: 10.1073/pnas.1006734107
Swanson, K. S., Dowd, S. E., Suchodolski, J. S., Middelbos, I. S., Vester, B. M., Barry, K. A., et al. (2010). Phylogenetic and gene-centric metagenomics of the canine intestinal microbiome reveals similarities with humans and mice. ISME J. 5, 639–649. doi: 10.1038/ismej,.2010.162
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., Ley, R. E., Mahowald, M. A., Magrini, V., Mardis, E. R., and Gordon, J. I. (2006). An obesity-associated gut microbiome with increased capacity for energy harvest. Nature 444, 1027–1031. doi: 10.1038/nature05414
Vanhoutte, T., Huys, G., Brandt, E., De F.ahey, G. C. J., and Swings, J. (2005). Molecular monitoring and characterization of the faecal microbiota of healthy dogs during fructan supplementation. FEMS Microbiol. Lett. 249, 65–71. doi: 10.1016/j.femsle.2005.06.003
Walter, J. (2008). Ecological role of lactobacilli in the gastrointestinal tract: implications for fundamental and biomedical research. Appl. Environ. Microbiol. 74, 4985–4996. doi: 10.1128/AEM.00753-08
Wang, Z., Klipfell, E., Bennett, B. J., Koeth, R., Levison, B. S., DuGar, B., et al. (2011). Gut flora metabolism of phosphatidylcholine promotes cardiovascular disease. Nature 472, 57–63. doi: 10.1038/nature09922
Wang, Z., Xiao, G., Yao, Y., Guo, S., Lu, K., and Sheng, Z. (2006). The role of bifidobacteria in gut barrier function after thermal injury in rats. J. Trauma. 61, 650–657. doi: 10.1097/01.ta.0000196574.70614.27
Wernimont, S. M., Radosevich, J., Jackson, M. I., Ephraim, E., Badri, D. V., MacLeay, J. M., et al. (2020). The effects of nutrition on the gastrointestinal microbiome of cats and dogs: impact on health and disease. Front. Microbiol. 11:1266. doi: 10.3389/fmicb.2020.01266
Wickham, H. (2009). ggplot2: Elegant Graphics for Data Analysis. New York, NY: Springer Publishing Company, Incorporated.
Wu, G. D., Chen, J., Hoffmann, C., Bittinger, K., Chen, Y., Sue, A., et al. (2011). Linking long-term dietary patterns with gut microbial enterotypes gary. Science 334, 105–108. doi: 10.1126/science.1208344
Xenoulis, P. G., Palculict, B., Allenspach, K., Steiner, J. M., Van House, A. M., and Suchodolski, J. S. (2008). Molecular-phylogenetic characterization of microbial communities imbalances in the small intestine of dogs with inflammatory bowel disease. FEMS Microbiol. Ecol. 66, 579–589. doi: 10.1111/j.1574-6941.2008.00556.x
Yan, H., Potu, R., Lu, H., Vezzoni De Almeida, V., Stewart, T., Ragland, D., et al. (2013). Dietary fat content and fiber type modulate hind gut microbial community and metabolic markers in the pig. PLoS ONE 8:e059581. doi: 10.1371/journal.pone.0059581
Zhang, C., Zhang, M., Pang, X., Zhao, Y., Wang, L., and Zhao, L. (2012). Structural resilience of the gut microbiota in adult mice under high-fat dietary perturbations. ISME J. 6, 1848–1857. doi: 10.1038/ismej.2012.27
Keywords: dogs, high fat diet, low carbohydrate, health, microbiota
Citation: Kilburn LR, Koester LR, Schmitz-Esser S, Serão NVL and Rossoni Serão MC (2020) High-Fat Diets Led to OTU-Level Shifts in Fecal Samples of Healthy Adult Dogs. Front. Microbiol. 11:564160. doi: 10.3389/fmicb.2020.564160
Received: 20 May 2020; Accepted: 13 November 2020;
Published: 08 December 2020.
Edited by:Luisa I. Falcon, National Autonomous University of Mexico, Mexico
Reviewed by:Andres Sanchez Quinto, National Autonomous University of Mexico, Mexico
M. Pilar Francino, Fundación para el Fomento de la Investigación Sanitaria y Biomédica de la Comunitat Valenciana (FISABIO), Spain
Copyright © 2020 Kilburn, Koester, Schmitz-Esser, Serão and Rossoni Serão. 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) and the copyright owner(s) 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.
*Correspondence: Mariana C. Rossoni Serão, firstname.lastname@example.org