Original Research ARTICLE
The effect of diet and host genotype on ceca microbiota of Japanese quail fed a cholesterol enriched diet
- 1The State Key Laboratory of Animal Nutrition, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, China
- 2Faculty of Land and Food Systems, Avian Research Centre, The University of British Columbia, Vancouver, BC, Canada
- 3School of Biological Sciences, The University of Hong Kong, Hong Kong, Hong Kong
Two Japanese quail strains, respectively atherosclerosis-susceptible (SUS) and –resistant (RES), have been shown to be good models to study cholesterol metabolism and transportation associated with atherosclerosis. Our objective was to examine possible difference in cecal microbiota between these strains when fed a control diet and a cholesterol enriched diet, to determine how host genotype and diet could affect the cecal microbiome that may play a part in cholesterol metabolism. A factorial study with both strains and two diets (control, cholesterol) was carried out. Cecal content was collected from 12 week old quail that have been on their respective diets for 6 weeks. DNA was extracted from the samples and the variable region 3–5 of the bacterial 16S rRNA gene was amplified. The amplicon libraries were subjected to pyrosequencing. Principal Component Analysis (PCA) of β-diversity showed four distinct microbiota communities that can be assigned to the 4 treatment groups (RES/control, RES/cholesterol, SUS/control, SUS/cholesterol). At the Phylum level, the 4 treatment groups has distinct Firmicutes community characteristics but no significant difference in Bacteroidetes. Eubacterium dolichum was rare in RES/control but became overabundant in RES/cholesterol. An unclassified species of Lactobacillaceae was found in abundance in SUS/control but the same species was rare in RES/cholesterol. On the other hand, two Lactobacillus species were only found in RES/control and an unclassified Lachnospiraceae species was abundant in RES/cholesterol but rare in SUS/control. The abundance of four species of Lachnospiraceae, three species of Ruminococcaceae and one species of Coprobacillaceae was positively correlated with plasma Total Cholesterol, plasma LDL, and LDL/HDL ratio. Our study of cecal microbiota in these quail has demonstrated that selection for susceptibility/resistance to diet induced atherosclerosis has also affected the quail's cecal environment to host distinctly different cecal microbiome.
Despite many measures are available for the management of cardiovascular disorders, this common disease is still associated with high incidence of mortality and morbidity (Go et al., 2014). This is a reflection that atherosclerosis is a complex pathological process, affected by both genetic and environmental factors that we still do not fully understand. Over the past several years, researchers have turned their attention to the effects of gastro-intestinal microbiota on the development of metabolic diseases (Cani and Delzenne, 2009; Caesar et al., 2010; Sekirov et al., 2010; Karlsson et al., 2012; Strowig et al., 2012). In particular, Wang et al. (2011) found that in mice, intestinal microbiota metabolism of dietary choline produced trimethylamine (TMA). TMA is further converted to trimethylamine N-oxide (TMAO) by liver enzymes Flavin monooxygenases, and TMAO has been found to promote endogenous cholesterol-laden macrophage foam cell formation, one of the earliest hallmarks of the atherosclerotic process. Koeth et al. (2013) found that L-carnitine, a TMA abundantly found in red meat, could also be converted to TMAO by intestinal microbiota. They also found that TMAO significantly reduced reverse cholesterol transport in mice. Never-the-less, most of the studies carried out so far lacked defined host genotypes and defined diets at the same time to examine their interaction effect on gut flora diversity.
The Japanese quail was introduced as a laboratory animal in the 1960s (Padgett and Ivey, 1959) and proved to be useful in many areas of biomedical studies (Minvielle et al., 2007; Cheng et al., 2010). At the University of British Columbia (UBC) Avian Research Centre, we have been maintaining two related strains of Japanese quail for research in atherosclerosis. The two strains were developed by divergent selection from a common foundation population (Shih et al., 1983). One strain is susceptible to diet-induced atherosclerosis (SUS), and the other is resistant (RES). When challenged with a high cholesterol diet, about 80% of the SUS males will develop atherosclerosis whereas only about 4% of the RES males will. Li et al. (2012) examined differential mRNA expression of seven genes involved in cholesterol metabolism and transport in the liver of the SUS and RES and concluded that these quail are good models for studying cholesterol metabolism and transport in relationship to atherosclerosis development. We therefore think that this quail model would be useful for studying the interaction of host genotype and diet in affecting the gut flora diversity in association with the development of atherosclerosis.
The objectives of our study were (1) to characterize the phylogenetic diversity of the cecal microbiota of the SUS and RES males fed a regular (control) diet, using 454 pyrosequencing after amplification for V3–V5 region of bacterial 16S rRNA gene, (2) to characterize the cecal microbiota diversity of the SUS and RES males fed the control diet dosed with cholesterol (0.5% w/w) using the same pyrosequencing procedure, and (3) to identify gut bacteria that are key to each of the four treatment groups (SUS/control diet, RES/control diet, SUS/cholesterol diet, RES/cholesterol diet) and to examine the association of these bacteria with the development of atherosclerosis.
Materials and Methods
The two strains of Japanese quail, SUS and RES, have been acquired by the UBC Quail Genetic Resource Centre from North Carolina State University in 1989. The history of their selective breeding has been described by Shih et al. (1983). Since their transfer to UBC, they have undergone further divergent selection for susceptibility and resistance to atherosclerotic plaque formation induced by dietary cholesterol (0.5%w/w) (Cheng et al., 1997).
After hatching, both SUS (N = 80) and RES (N = 80) males were fed a semi-synthetic diet (Li et al., 2012) (Table 1) prepared by the feed mill at the Agriculture and Agri-Food Canada Poultry Research Station at Agassiz, British Columbia, according to the NRC nutrient requirements standards recommended for Japanese quail (http://www.nap.edu/catalog/2114.html). At 6 weeks of age, they were divided into two dietary treatment groups and fed either a regular synthetic diet (control) or a synthetic diet with added cholesterol (0.5%w/w) for another 6 weeks (Li et al., 2012). Individually marked birds (both RES and SUS) fed the same diet were kept in the same pen. Birds on the alternative diet were kept in a neighboring pen. At 12 weeks of age, 6 birds from each of the treatment groups with body weight closest to the mean of the population were euthanized by decapitation and trunk blood was collected into Vacutainer tubes (Becton–Dickinson, Mississauga, ON, Canada), containing lithium heparin, and centrifuged at 4°C for 10 min at 3000 × g. Plasma was stored at −20°C until it was later used for lipid analysis. Sections of ceca, including gut content were collected from each bird. All samples were quick frozen on dry ice immediately after collection and stored at −70°C until processed for DNA extraction. The aortic tree (the brachycephalic arteries to their bifurcations and the aorta to the iliac branching) of each bird was dissected out, opened longitudinally and examined under a 10–30X dissecting microscope for a semi-quantitative scoring of the seriousness of the atherosclerotic lesions on the interior wall. The scoring system was adopted from Godin et al. (1995). A score of 0 (normal) to 4 (presence of severe atherosclerotic lesions) was assigned by two independent scorers who were blind to the genetic and diet status of the bird. Four quail from each treatment group were selected, based on their atherosclerotic lesion scores, for examination of their cecal microbiota. This research was carried out with the approval of the UBC Animal Care Committee (Certificate # A12-0087).
Plasma Lipid Analysis
Plasma samples were sent to the Department of Pathology and Laboratory Medicine at St. Paul's Hospital (Vancouver, BC) and assayed for total cholesterol, HDL, and triglycerides using enzymatic methods on an ADVIA 1650 Chemistry System. Any sample with 3+ lipemia or greater (as measured by the analyzer), was cleared by Airfuge® Air-Driven ultracentrifugation (Beckman Coulter). Lipemia at that level can cause interference with the HDL assay. HDL was assessed by the direct method without precipitation of apolipoprotein B (Warnick and Albers, 1978; Warnick et al., 2001; Gootjes et al., 2009). LDL values were calculated by Friedewald's formula, using measured values for total cholesterol, HDL and triglycerides (Friedewald et al., 1972; Okada et al., 1998).
DNA Extraction and Pyrosequencing
The intestinal segments were thawed and the contents were gently scraped from the intestinal wall. The surgical tools and vials were autoclaved and the bench area was wiped clean with 70% ethanol to minimize contamination. Genomic DNA was isolated using the PowerMax Soil DNA Isolation Kit (Mo. Bio laboratories. Inc., Carlsbad, CA) according to the instructions of the manufacturer with 200 mg as starting material. PCR amplifications were performed using the FastStart high fidelity PCR system (Roche Molecular Diagnostics, Branchburg, NJ, USA). The variable region 3–5 (V3–V5) of the bacterial 16S rRNA gene was amplified with a primer set of 341F (5′- ACTCCTACGG GAGGCAGCAG-3′) and 926R (5′- CCGTCAATTCMTTTGAGTTT-3′) with the sample specific forward primer bearing a multiplex identifier (MID) sequences. All 341F and 926R primers modified with adaptor A and B sequences respectively for pyrotag sequencing. The amplification program consisted of an initial denaturation step at 94°C for 2 min; 32 cycles of denaturation at 94°C for 30 s, annealing at 60°C for 30 s, and elongation at 72°C for 30 s; and a final extension step at 72°C for 7 min. The size of the PCR products was confirmed by gel electrophoresis. The PCR products was then purified using Gel extraction kit (Invitrogen) and were quantified using the NanoDrop 2000 (Thermo Scientific, Wilmington, DE, USA). The Amplicon libraries were subjected to pyrotag sequencing using a bench-top 454 GS Junior (454 Life Sciences-a Roche Company, Branford, CT, USA) with the GS Junior Titanium Sequencing Kit (https://lifescience.roche.com/shop/en/us/products/gs-junior-titanium-sequencing-kit).
Sequences obtained from pyrosequencing were processed using the QIIME (quantitative insights into microbial ecology) software package (Caporaso et al., 2010b). Quality trimming of dataset removed sequences if a mean quality score was ≤ 25; lengths were < 150 or >900 bp; sequences were without primer, uncorrectable, or contained ambiguous characters; or homopolymer run exceeding 8 nt. De-noising of dataset was performed using DENOISER v. 0.9.1 (Quince et al., 2011) as implemented in QIIME platform. Chimeric sequences were removed using Chimera Slayer. The sequences were assigned to groups basing on their respective barcode sequences. Similar sequences were assigned into operational taxonomic units (OTUs) at a pairwise identify of 97% using UCLUST (http://www.drive5.com/usearch/). Representative sequence was the most abundant sequence in each OTU. Representative sequences (at 97% similarity) were then classified taxonomically using Ribosomal Database Project (RDP) classifier 2.0.1 (Cole et al., 2009). The OTUs were aligned using PyNAST with a minimum alignment length of 150 bp and a minimum percent identity of 75% (Caporaso et al., 2010a). After alignment, PH LANE mask (http://greengenes.lbl.gov/) was conducted to screen out the hypervariable regions.
Richness and Diversity Indices
Rarefaction plots were constructed and diversity indices (Chao1 richness, Simpson's Diversity) were estimated as implemented in QIIME (Caporaso et al., 2010b). For the comparison of β-diversity among microbial communities, we used principal component analysis (PCA) to visualize all OTUs and OTUs in phylum level (Firmicutes) differences. Results of the PCA were then statistically tested by permutational multivariate analysis (PERMANOVA) of variance (Anderson et al., 2008) for ceca microbiota compositions differences among four treatment groups. Mahalanobis distance (de Maesschalck et al., 2000) was calculated to confirm the difference between every two cluster.
Comparison of Microbial Communities
Bacterial abundance difference on phylum-, family-, and species- levels were examined using multivariate analysis, and further using Tukey's HSD for mean separation (SPSS 13.0; SPSS Institute, 2001) and expressed as means ± SE.
Least squares analysis of variance was performed to compare the diversity and richness parameters using JMP 8.0 (SAS Institute, North Carolina, 2008). The statistical model as follow:
Where Yijk represents the measure for the kth individual of the ith strain from jth diet; Si = RES or SUS bird; Dj = control or a cholesterol diet; (SDij) = the two-way interaction term and Eijk = the error term. The results were reported as the least square mean ± standard error of means (SEM). Tukey's HSD was used for mean separation and statistical significance was defined at P < 0.05.
Key Microbiota in Each Treatment Group
We used Venn Diagram (Oliveros, 2007) and a “Nearest-shrunken Centroid” (NSC) classification approach (Tibshirani et al., 2002; Koren et al., 2011) to detect core microbiota community which best characterize each group. The amount of shrinkage was determined by cross-validation and test error was minimized. We selected the OTUs from all 6 Phyla that are common in at least 3 of 4 samples examined per treatment, and with a minimum sequence count of 100 sequences per OTU to generate a filtered OTU table for the Venn diagram and NSC analyses. NSC analysis was performed on normalized Z-score profiles of OTUs. The misclassification error was 0.25 and threshold was 1.13.
Correlation of Abundance of OTUs with Blood Lipid Parameters
A correlation heatmap was generated to examine the correlation between the abundance of particular OTUs with the levels of blood lipid parameters—plasma Total Cholesterol, plasma LDL, plasma HDL, plasma Triglycerides, and LDL/HDL ratio. Because the levels of blood lipid parameters were not independent of dietary cholesterol, only RE and SE were included in the analysis. Because of small sample size, correlation plots of significant correlations were examined for data point distribution to eliminate correlations due to outliers.
Atherosclerotic Lesions on the Intimal Surface of the Aortae
All SUS and RES fed the control diet scored 0. All four SUS on cholesterol diet scored 4, while two RES on the same diet score 0 and two scored 1.
Richness of Ceca Microbiota
After trimming, assembly and quality filtering, a total of 257,860 sequence reads were obtained with a mean reading length of 545 bp and 16,116 ± 4269 reads/sample. There was no significant difference among the four treatment groups in the number of sequence reads.
Rarefaction curves (Figure 1A) of the quality filtered sequences showed that our sequencing was deep enough to recover almost all of the OTUs in the sampled population [RES/control diet (RC) 150.0 ± 25.1 OTUs/sample, SUS/control diet (SC) 125.8 ± 14.4 OTUs /sample, RES/cholesterol diet (RE) 94.0 ± 8.5 OTUs /sample, SUS/cholesterol diet (SE) 122.2 ± 6.9 OTUs /sample]. Rarefaction curves (Figure 1B) also indicated that four samples per treatment group was a big enough sample size to detect most of the core OTUs in the population (Hughes and Hellmann, 2005).
Figure 1. Rarefaction analysis, calculated at 97% dissimilarity, for the assessment of operational taxonomic unit (OTU) coverage within the16S rRNA gene-based cecal bacterial communities in the RES and SUS quail fed the control (RC, blue open square; SC, green open square) or cholesterol (RE, orange closed square; SE, red closed square) diets. (A) The number of OTUs as a function of the number of sequence reads. (B) The number of OTUs as a function of the number of individual quail sampled.
The sequences were classifiable into 366 species-level operational taxonomic units (OTUs) (123 ± 24.7 OTUs/sample) belonging to 6 bacterial phyla. The vast majority (98%) of the sequences belonged to two bacterial phyla: Firmicutes (77%) and Bacteroidetes (21%).
The remaining sequences were identified as Spirochaetes, Tenericutes, Proteobacteria, and Actinobacteria. Because of very low number of sequences, we did not compare bacterial abundance differences in these phyla. The sequences have been submitted to Sequence Read Archive (SRE) with accession number SRR2537231.
A Comparison of Ceca Microbiota Diversity
Community Level Variations
Chao1 estimator indicated that there was a significant (P = 0.002) diet × genotype interaction affecting richness. RES on the control diet (RC) had significantly higher OTU richness (164 ± 11.7) than RES on the cholesterol diet (RE: 104 ± 2.5), whereas there was no significant difference in richness between the two SUS dietary groups (SC: 138 ± 7.2, SE; 133 ± 2.1). There was no significant difference (P > 0.05) in diversity among the 4 treatment groups according to the Simpson estimate of diversity.
Principal components Analysis (PCA) of β-diversity showed four distinct microbiota communities that can be assigned to the four treatment groups (Figures 2, 3). There was a significant (PERMANOVA; P = 0.02) diet × genotype interaction indicating these four treatment groups had distinct microbial community characteristics.
Figure 2. Three-dimensional projection of PCA of whole cecal microbial community. The variance explained by the PCs is indicated in parentheses on the axes. Each symbol represents a single sample.
Figure 3. Clustering of gut microbiota based on distances between different groups calculated with multivariate analysis of variance test of the first six PCs of OTUs data. The Mahalanobis distances between group means are shown. **P < 0.01.
Phylum Level Variations
PCA was also conducted for OTUs within a Phylum. Figure 4 showed the PCA of Firmicutes OTUs. There was a significant (PERMANOVA; P = 0.02) diet × genotype interaction indicating these four treatment groups had distinct Firmicutes community characteristics. Figure 5 showed the Mahalanobis distance between the groups to confirm their clustering. All Bacteroidetes belonged to the family Rikenellaceae. There was no significant difference among the treatment groups when OTUs within Bacteroidetes were analyzed.
Figure 4. Joint Plot of PCA of Phylum Firmicutes in cecal content. The variance explained by the PCs is indicated in parentheses on the axes. Each symbol represents a single sample. The vector projections showed the relationship between family-level variables in Firmicutes and four treatment groups. The angle and length of the arrows indicated the direction and strength of the relationship. f1, Bacilli Lactobacillales Lactobacillaceae; f2, Bacilli Lactobacillales Streptococcaceae; f3, Unclassified Bacilli; f4, Clostridia Clostridiales Clostridiaceae; f5, Clostridia Clostridiales Lachnospiraceae; f6, Clostridia Clostridiales; f7, Unclassified Clostridia Clostridiales; f8, Unclassified Clostridia; f9, Erysipelotrichi Erysipelotrichales Coprobacillaceae; f10, Erysipelotrichi Erysipelotrichales Erysipelotrichaceae.
Figure 5. Clustering of Firmicutes based on distances between different groups calculated with multivariate analysis of variance test of the first six PCs of OTUs data. The Mahalanobis distances between group means are shown. **P < 0.01.
Family Level Variations
We used Joint Plot of PCA (PC-ORD) to examine the correlation of microbiota families (within Firmicutes) with the four treatment groups (indicated by the arrows in Figure 4). In general, birds on control diet were characterized by the abundance of Lactobacillaceae and Streptococcaceae. Birds on cholesterol diet were characterized by the abundance of Erysipelotrichaceae and Clostridiaceae. Specifically, cecal microbiota in RC was characterized by the relative abundance of unclassified Clostridiales, Bacilli and Clostridia, and SC was characterized by Coprobacillaceae. SE was characterized by Lachnospiraceae, while RE by Erysipelotrichaceae.
Using multivariate analysis, we found a significant (P < 0.012) diet × genotype interaction in affecting the abundance of Ruminococcaceae (Table 2). RC had significantly more than RE, but there was no significant difference for SUS on the two different diets. The abundance of Erysipelotrichaceae, unclassified family of Bacilli (UB), and unclassified family of Clostridia (UC) was significantly affected by diet. Compared to birds (SUS and RES) on control diet, birds on cholesterol diet had significantly more Erysipelotrichaceae, significantly less UB (P < 0.04) and UC (P < 0.03).
Genus Level Variations
Using multivariate analysis, we found a significant (P < 0.0001) diet effect in the abundance of Ruminococcus (Table 3) (Control diet: 125.63 ± 17.88, Cholesterol diet: 16.13 ± 15.36). We have also found a significant (P < 0.002) diet effect in the abundance of Cc_115 (belonging to Family Erysipelotrichaceae). There was also a significant (P < 0.05) diet × genotype interaction in Unclassified Ruminococcaceae (Table 3).
Venn diagram showed the distribution of 66 filtered OTUs (Figure 6). With much overlap, birds on the control diet harbored 54 OTUs while birds on the cholesterol diet harbored 44 OTUs. Twelve OTUs were unique to birds on the cholesterol diet, while 22 OTUs were unique to birds on the control diet (Table 4). Twelve OTUs were unique to RES birds and 13 were unique to SUS birds. Ten OTUs (7.36%) were unique to RC, 6 (1.57%) unique to SC, 2 (0.71%) to RE, and 7 (2.36%) to SE.
NSC shown rarity of an unclassified species of Ruminococcus (ID 182245) and overabundance of two unclassified species of Rikenellaceae (ID 4336943, 157573, respectively) in RC; Overabundance of ID 182245 in RE; Rarity of ID 182245 and overabundance of a different unclassified species of Ruminococcus (ID 185972) in SC; Overabundance of ID 182245 and a third species of unclassified Ruminococcus (ID 548503) in SE.
Table 5 summarized the results of the Venn Diagram/NSC analyses. Eubacterium dolichum was rare in RC but became overabundant when RES was fed the cholesterol diet (RE). An unclassified species of Lactobacillaceae was found in abundance in SC but the same species was rare in RE. On the other hand, two Lactobacillus species (also in the Lactobacillaceae Family) were only found in RC, and an unidentified Lachnospiraceae species was abundant in RE but rare in SC.
Association of Key Bacteria Species with Plasma Lipid Parameters
There was a significant diet × genotype interaction in plasma Total Cholesterol (TC) (P < 0.006; Table 6), LDL (P < 0.004; Table 7) levels, and LDL/HDL ratio (P < 0.004 Table 8). SE was significantly higher in these parameters than the other three treatment groups.
Table 6. Significant (P < 0.006) Diet × Genotype interaction in plasma total cholesterol level (mmol/L).
The abundance of four species of ceca bacteria in Lachnospiraceae, three species in Ruminococcaceae, and one species in Coprobacillaceae were positively correlated with the plasma TC, plasma LDL, and LDL/HDL ratio (Figure 7 and Table 9). The abundance of one species of Lachnospiraceae was positively correlated with plasma HDL while the abundance of two species of Lachnospiraceae and Eubacterium dolichum was negatively correlated with plasma HDL level. The abundance of four species (Unclassified Rikenellaceae (157573), Unclassified Oscillospira (366143), Unclassified Ruminococcaceae (295861), and Unclassified Cloacamonaceae (NCUR 160)) were found to be significantly and positively correlated with plasma triglycerides level due to a single outlier point and were eliminated.
Figure 7. Correlations heat map demonstrating the association between the abundances of different cecal microbial species and plasma lipid parameters. Correlation coefficients are represented by color ranging from blue, negative correlation, to red, positive correlation. Significant correlations are noted by *p < 0.05 and **p < 0.01.
Table 9. Significant Pearson's Correlations§ between the abundance of cecal bacteria species and plasma lipid parameters.
We examined the cecal microbiota of 12 week old quail that had been fed their respective diets for at least 6 weeks. The cecal microbiota should be mature and stable by that time (Lu et al., 2003). Taxonomic analysis showed that composition of the quail's cecal microbiota at various levels is similar to that of human, mice, hamsters, chickens, emu, and Bobwhite quail (Eckburg et al., 2005; Ley et al., 2005; Turnbaugh et al., 2008; Karlsson et al., 2012; Bennett et al., 2013; Su et al., 2014; Videnska et al., 2014). At the Phylum level, Bacteroidetes and Firmicutes dominated the cecal microbiota, representing over 98% of all the sequences detected. An unusual feature of the quail cecal microbiota was that all (>99.9%) Bacteroidetes were of the Family Rikenellaceae. At the Family level, the predominant taxons were Rikenellaceae, Lactobacillaceae, Streptococcaceae, Lachnospiraceae, Coprobacillaceae, and Erysipelotrichaceae. At the Genus level, Ruminococcus, Blautia, Coprococcus, Eubacterium were abundant. While in most cases, it was not possible to identify the OTUs down to the species level because of high similarity of the 16S sequences among different bacterial species, we were able to identify a few OTUs down to the species level (e.g., see Table 4) because of their unique 16S sequences. We have BLAST searched these OTU sequences and found 100% sequence similarity with the species identified. Since even the closest species members of those bacteria have at least 5% sequence dissimilarity, we are confident that we have correct taxonomic identification for these OTUs.
The Effect of Host Genotype on Cecal Microbiota When Birds were on Control Diet
The SUS and RES quail strains are a result of divergent selection from a common foundation population. The selection criteria being the highest and lowest atherosclerotic plaque scores, respectively, when the birds were fed a diet containing 1% cholesterol (w/w) (Shih et al., 1983). As a correlated response, we found that the selection not only modified the quail cecal environment to host a different cecal microbiome when the birds were fed a regular diet, but they also reacted differently when fed a high cholesterol diet. In birds, embryos develop in the egg outside the mother's body and parental influence should be minimum. In our study, the eggs of the four treatment groups were artificially incubated at the same time in close proximity in the same incubator. The birds of the two strains fed the same diet were raised in the same pen. The housing density was low enough to minimize any possibility of one strain dominating the other and affected their feed intake. This was confirmed by the fact that body weight and the mortality rate were not different between the two strains. The difference in cecal microbiota between RC and SC can therefore be attributed to host genetic differences in affecting the cecal environment.
While there was no significant difference in CHAO1 richness and Simpson diversity between RC and SC cecal microbiome, PCA analysis of cecal microbiota at the Phylum level detected two distinct Firmicutes communities between RC and SC. At the family level, RC was characterized by significant abundance of Ruminococcaceae. At the OTU level, 10 OTUs were unique to RC. Combining the results from the Venn diagram and NSC analyses, we concluded that RC hosted two unique species, Lactobacillus reuteri (ID 137043) and an unclassified Lactobacillus (ID 1375808), in significant abundance. In comparison with other treatment groups, RC has less abundance of Eubacterium dolichum (ID 229069) but more abundance in an unclassified species of Clostridia (ID 178304). Ruminococcaceae is one of the two most abundant families from the order Clostridiales found in the mammalian gut environment, and have been associated with the maintenance of gut health . Ruminococcaceae is also the most common family of microbes in chicken cecum (Apajalahti and Kettunen, 2006; Torok et al., 2011). L. reuteri has been commonly used as a probiotic to suppress GI tract inflammation in human (Shornikova et al., 1997). Resting cells of this species convert glycerol into a potent, broad-spectrum antimicrobial substance termed reuterin (Axelsson et al., 1989). Dietary supplementation with L. reuteri ATCC 55730 resulted in significant colonization in the stomach, duodenum, and ileum of healthy humans, and this is associated with significant improvements of the immune response in the gastrointestinal mucosa (Valeur et al., 2004). Selection for resistance to diet induced atherosclerosis may have improved the general gut (cecum) health of RES quail.
At the family level, SC was characterized by abundance of Coprobacillaceae. Out of the six OTUs identified by Venn Diagram Analysis as unique in SC, only one was deemed significantly abundant by NSC analysis: an unclassified Clostridia species (ID 188057). Comparing with other treatment groups, SC has more abundance of four Lachnospiraceae species, including Blautia producta (ID 158211) and unclassified Coprococcus (ID 186319). SC also has more abundance of two unclassified Coprobacillaceae species (ID 136526 and 592616), an unclassified Lactobacillaceae species (ID 292057), and an unclassified Rikenellaceae (Bacteroidetes) species (ID 4476780). In the less abundant category, SC has an unclassified Clostridia (ID 178304) which was significantly more abundant in RC, and an unclassified Lachnospiraceae species (ID 211212). Although Coprobacillaceae is found in most of the microbiome data set of different hosts, it is not a well-studied Family and not much information is available about activities of members of this Family (Verbarg et al., 2014). Lachnospiraceae species are cellulose-degrading bacteria prevalent in bovine gut samples. They produce butyric acid to degrade plant fiber such as xylans. In human, butyrate arising from such microbial fermentation is important for the energy metabolism and normal development of colonic epithelial cells and has a mainly protective role (activating immune/inflammatory responses) in relation to colonic disease (Pryde et al., 2002; Maslowski et al., 2009; Vinolo et al., 2011). Mice precolonized with a murine Lachnospiraceae isolate had significantly decreased Clostrium difficile (a pathogen) colonization, lower intestinal cytotoxin levels and exhibited less severe clinical signs and colonic histopathology (Reeves et al., 2012). Many Clostridia species are toxigenic (Hatheway, 1990). The question remains whether the unclassified Clostridia species (ID 188057) that was only found in SC is pathogenic. Videnska et al. (2014) examined the succession and replacement of bacterial population in the cecum of laying hens at various ages and found that cecal microbiota in young chicks were dominated by Firmicutes, but as the hen became sexually mature and started egg production, a gradual succession of the representatives of Firmicutes and also their replacement with the representatives of Bacteroidetes was observed (Videnska et al., 2014). Rikenellaceae only started to colonize the cecum when the hen was in full egg production. Atherosclerosis is also accelerated by aging (Weingand et al., 1986; Clarkson et al., 1987; Collins et al., 2009). The hypothesis that selection for susceptibility to diet-induced atherosclerosis has also pre-maturely aged the cecal environment is worth testing.
Divergent selection for resistance and susceptibility to dietary cholesterol induced atherosclerosis has shifted the cecal microbiome of Japanese quail in different ways.
The Effect of Dietary Cholesterol on Cecal Microbiota
The diet we used for the study was a synthetic diet where all the dietary ingredients were known and standardized. The composition of the control and experimental diets was identical except for the added cholesterol and a small amount (0.02% w/w) of Cholic acid to aid the digestion and absorption of cholesterol. The difference in cecal microbiota between birds fed the control diet and the experimental diet can therefore be attributed to dietary cholesterol.
Dietary cholesterol has major effects on cecal microbiota in RES and SUS alike. At the Family level, birds on cholesterol diet had significantly more Erysipelotrichaceae, but significantly less unclassified Bacilli and unclassified Clostridia than birds on control diet. At the genus level, birds fed the cholesterol diet had significantly less abundance of unclassified Ruminococcus. All members of the Erysipelotrichaceae Family were associated with one or several hosts including mammals, birds, fish and marine invertebrates, and most members were found as opportunistic pathogens affecting various parts of the body (Verbarg et al., 2014). Dietary cholesterol has reduced the abundance of Ruminococcus and facilitated the abundance of opportunistic pathogens in the quail ceca and may have increased the risk of assaults by these opportunistic pathogens.
The Effect of Genotype × Diet Interaction on Cecal Microbiota
When SUS and RES were put on a high cholesterol diet, they also reacted differently, in terms of the microbiota that they were hosting, to the dietary cholesterol. CHAO1 richness of cecal microbiome in RE was significantly reduced when compared with RC, but CHAO1 richness of cecal microbiome in SUS was not affected by diet.
At the family level, cecal microbiota in RC was characterized by the relative abundance of Ruminococcaceae, unclassified Clostridiales, Bacilli and Clostridia, and RE by abundance of Erysipelotrichaceae. RC had significantly more Ruminococcaceae than RE. SC was characterized by abundance of Coprobacillaceae, but SE was characterized by abundance of Lachnospiraceae. The abundance of Ruminococcaceae was not affected by diet in SUS. At the OTU level, Eubacterium dolichum was rare in RC but became overabundant in RE. Both SC and SE showed rarity of E. dolichum.
When mice were fed a “Western diet” which was high in fat and cholesterol, the overall diversity of their gut microbiota dropped significantly due to a bloom of a class of Firmicutes called Mollicutes, a member of which is E. dolichum (Turnbaugh et al., 2008). E. dolichum has a number of genomic features that could promote their own fitness in competition with other microbes in the cecal nutrient metabolic milieu created by the host's consumption of the Western diet (Turnbaugh et al., 2008). Their abundance is associated with obesity in mice. A similar situation may have occurred in RE in their reaction to dietary cholesterol. SE has abundance of Lachnospiraceae. At the same time, their abundance of Ruminococcaceae was not compromised by dietary cholesterol. Lachnospiraceae and Ruminococcaceae have been associated with the maintenance of gut health (Place et al., 2005; Huws et al., 2011; Vinolo et al., 2011; Reeves et al., 2012; Biddle et al., 2013; Greer et al., 2013). These two families are specialists for degrading cellulose and hemicellulose components of plant materials which are fermented and converted into short chain fatty acids (SCFAs) be absorbed and used by the host (Biddle et al., 2013). SCFAs have an important roles in maintaining intestinal homeostasis (Pryde et al., 2002; Place et al., 2005; Cotta and Forster, 2006; Wong et al., 2006; Greer et al., 2013). The results from our study seem to indicate that the divergent selection for susceptibility/resistant to diet induced atherosclerosis has adversely affected the cecal health of RE but not SE, via their cecal microbiome. Whether this change in the cecal environment has effects on the metabolism and absorption of dietary cholesterol remains to be studied.
Cecal Microbiota and Atherosclerosis
Recently there has been a flourish of studies on the relationship between gut microbiota and cardiovascular diseases in human and in animal models. In human, about 50% of dietary cholesterol is absorbed in the duodenum. All cholesterol arriving in the large intestine can be metabolized by Eubacterium bacteria to coprostanol and minor amounts of coprostanone (Macdonald et al., 1983). Coprostanol, unlike cholesterol, is poorly absorbed by the human intestine and hence, conversion of cholesterol to coprostanol might be a way to lower serum cholesterol in human and rodents (Sekimoto et al., 1983; Li, 1995; Stepankova et al., 2010). However, feeding Eubacterium coprostanoligenes to laying hens failed to lower plasma cholesterol (Li et al., 1996). In our study, Eubacterium dolichum was found in abundance in the RE cecum but not in SE. However, the ability of E. dolichum to convert cholesterol to coprostanol has not been demonstrated although we have found a significant but negative correlation of E. dolichum abundance with plasma HDL level. On the other hand, the primary cholesterol absorption sites are in the small intestine and it will be worthwhile to examine the microbiota in duodenum and ileum (S Liu and KM Cheng, study in progress).
In human and mice, intestinal microbes (Clostridium, Peptostreptococcaceae, Tenerites, and Clostridiaceae) can catabolize choline and L-carnitine to gaseous trimethylamine (TMA) (Al-Waiz et al., 1992; Koeth et al., 2013) which can be efficiently absorbed and metabolized by hepatic enzymes, Flavin monooxygenase 3 (FMO3), to form TMAO, an oxidized product of TMA (Cashman et al., 2003). TMAO promotes atherosclerosis by up-regulation of macrophage scavenger receptors (Wang et al., 2011) and by down-regulating genes involved in reverse cholesterol transportation (Koeth et al., 2013). The dietary sources of choline are foods such as eggs, milk, red meat, liver, shell fish, and fish, which are rich in lecithin (Wang et al., 2011). The dietary source of L-carnitine is red meat (Koeth et al., 2013). None of these food items are in the natural diet of Japanese quail. The synthetic diets that we have prepared for the study are also very low in choline (0.04% w/w of choline chloride) and L-carnitine. The atherosclerosis we found in the SE may not be promoted through this pathway but Shih et al. (1983) and Godin et al. (2003) reported that after cholesterol feeding, plasma cholesterol levels remained high for a significantly longer time in SE than in RE (Shih et al., 1983). Li et al. (2012) also found down regulation of some of the cholesterol transport genes in the SE liver (Li et al., 2012). It will be worthwhile to examine the expression of hepatic FMO3 and RCT associated genes which can be regulated by TMAO, in the SE compared with RE (JE Kim and KM Cheng, study in progress).
It has been proposed that bacterial lipopolysaccharide (LPS), a constituent of Gram negative bacteria present in the gut microbiota, can be transported from the intestine to target tissue and combine with CD14 and the toll-like receptor 4 (TLR4) at the surface of innate immune cells such as macrophages. Such “metabolic endotoxemia” can trigger the secretion of proinflammatory cytokines. Efflux of cholesterol from vessel wall macrophages is believed to be a critical first step by which RCT protects against atherosclerosis. TLR4 inhibits RCT and thus may modulate cholesterol metabolism. These findings suggest that gut microbiota may be important for RCT but it is not clear whether gut microbiota contributes to atherosclerosis through this pathway (Caesar et al., 2010). In our study, RE is the only group that has abundant gram negative bacteria (E. dolichum) in their cecal microbiota and yet they were resistant to diet induced atherosclerosis. On the other hand, the abundance of several gram positive species (Lachnospiraceae and Ruminococcaceae) was positively correlated with plasma TG and plasma LDL levels. The hypothesis that selection for susceptibility to diet induced atherosclerosis in Japanese quail has also shifted their gut microbiota to enhance the metabolism and absorption of cholesterol remains to be tested.
Several animal models have been developed for studying atherosclerosis, but each has advantages and limitations (Jokinen et al., 1985; Getz and Reardon, 2012; Kapourchali et al., 2014). The Japanese quail model may have advantages over others because quail are naturally deficient in apolipoprotein E. When fed a high cholesterol diet, males of the SUS strain developed lesions exhibiting structural features (e.g., focal hemorrhage, calcification and fibrosis) that closely resemble those in the human disorder (Shih et al., 1983). Atherogenic diets increase the LDL and VLDL fraction of cholesterol with minimal effects on HDL, thus facilitating the study of metabolism and transport of cholesterol in relationship to atherosclerosis (Li et al., 2012). Our study of cecal microbiota in these quail has demonstrated that selection for susceptibility/resistance to diet induced atherosclerosis has also affected the quail's cecal environment to host distinctly different cecal microbiome. Moreover, the SUS and RES quail also reacted differently, in terms of the cecal microbiota that they are hosting, to dietary cholesterol. Our study allowed us to raise new questions about the relationship between gut microbiota and cholesterol metabolism.
This manuscript is an extension of the thesis research carried out by SL. DB provided expertise in quail nutrition and gut microbiota. HT provided expertise in Bioinformatics. JK provided expertise in DNA extraction. FL provided expertise and laboratory facilities for pyrosequencing. HZ was SL's thesis supervisor at CAAS. KC was the thesis research supervisor and provided expertise in genetics, cholesterol metabolism and atherosclerosis, and experimental design.
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.
This paper was submitted by SL as part of a thesis for the partial fulfillment of the requirements for a Ph.D. degree, Chinese Academy of Agricultural Sciences/China Agricultural University. The study was carried out while SL was a visiting student at the UBC Avian Research Centre. The funding for this research was provided by the BC Ministry of Agriculture (funds administered by the UBC Specialty Birds Research Committee) to KC, and by the Strategic Initiative (Biomedical Engineering) Fund to FL. SL was supported by a scholarship from the China Scholarship Council. The feeding trial was carried out at the Agriculture and Agri-Food Canada Research Station at Agassiz, under the supervision of Dr. Frederick Silversides. We would like to thank Dr. Raymond Kin-hi Hui, Dr. Tae-Jin Park, Ms. Charis Chan, and Ms Angel Ma at The University of Hong Kong for technical assistance; Dr. Xiaonan Lu (Food, Nutrition, and Health Program, UBC) for a critical review of an earlier draft of the manuscript.
FMO3, Flavin monooxygenase 3; LPS, lipopolysaccharide; OUTs, operational taxonomic units; PCA, Principal Components Analysis; RC, atherosclerosis-resistant quail on control diet; RCT, reverse cholesterol transport; RE, atherosclerosis-resistant quail on cholesterol diet; RES, atherosclerosis-resistant quail; SC, atherosclerosis-susceptible quail on control diet; SE, atherosclerosis-susceptible quail on cholesterol diet; SUS, atherosclerosis-susceptible quail; TG, triglycerides; TLR4, toll-like receptor 4; TMA, trimethylamine; TMAO, trimethylamine N-oxide.
Bennett, D. C., Tun, H. M., Kim, J. E., Leung, F. C., and Cheng, K. M. (2013). Characterization of cecal microbiota of the emu (Dromaius novaehollandiae). Vet. Microbiol. 166, 304–310. doi: 10.1016/j.vetmic.2013.05.018
Biddle, A., Stewart, L., Blanchard, J., and Leschine, S. (2013). Untangling the genetic basis of fibrolytic specialization by Lachnospiraceae and Ruminococcaceae in diverse gut communities. Diversity 5, 627–640. doi: 10.3390/d5030627
Caesar, R., Fåk, F., and Bäckhed, F. (2010). Effects of gut microbiota on obesity and atherosclerosis via modulation of inflammation and lipid metabolism. J. Intern. Med. 268, 320–328. doi: 10.1111/j.1365-2796.2010.02270.x
Caporaso, G., Bittinger, K., Bushman, F. D., DeSantis, T. Z., Andersen, G. L., and Knight, R. (2010a). PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics 26, 266–267. doi: 10.1093/bioinformatics/btp636
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., et al. (2010b). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Cashman, J. R., Camp, K., Fakharzadeh, S. S., Fennessey, P. V., Hines, R. N., Mamer, O. A., et al. (2003). Biochemical and clinical aspects of the human flavin-containing monooxygenase form 3 (FMO3) related to trimethylaminuria. Curr. Drug Metab. 4, 151–170. doi: 10.2174/1389200033489505
Cheng, K. M., Aggrey, S., Nichols, C. R., Garnett, M. E., and Godin, D. V. (1997). Antioxidant enzymes and atherosclerosis in Japanese quail: heritability and genetic correlation estimates. Can. J. Cardiol. 13, 669–676.
Cheng, K. M., Bennett, D. C., and Mills, A. D. (2010). “Chapter 42: The Japanese quail,” in UFAW Handbook on the Care and Management of Laboratory Animals, 8th Edn., eds R. Hurbrecht and J. Kirkwood (London: Blackwell), 655–673.
Clarkson, T. B., Adams, M. R., Weingand, K. W., Miller, L. C., and Heydrick, S. (1987). “Effect of age on atherosclerosis progression in nonhuman primates,” in Atherogenesis and Aging, eds S. R. Bates and E. C. Gangloff (New York, NY: Springer), 57–71.
Cole, J. R., Wang, Q., Cardenas, E., Fish, J., Chai, B., Farris, R. J., et al. (2009). The ribosomal database project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 37, D141–D145. doi: 10.1093/nar/gkn879
Collins, A. R., Lyon, C. J., Xia, X., Liu, J. Z., Tangirala, R. K., Yin, F., et al. (2009). Age-accelerated atherosclerosis correlates with failure to upregulate antioxidant genes. Circ. Res. 104, e42–e54. doi: 10.1161/CIRCRESAHA.108.188771
Cotta, M., and Forster, R. (2006). “Chapter 1.3.5: The family Lachnospiraceae, including the genera Butyrivibrio, Lachnospira and Roseburia,” in The Prokaryotes, eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackbrandt, and F. Thompson (New York, NY: Springer), 1002–1021.
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
Friedewald, W. T., Levy, R. I., and Fredrickson, D. S. (1972). Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin. Chem. 18, 499–502.
Go, A. S., Mozaffarian, D., Roger, V. L., Benjamin, E. J., Berry, J. D., Blaha, M. J., et al. (2014). Heart disease and stroke statistics—2014 update: a report from the American Heart Association. Circulation 129, e28. doi: 10.1161/01.cir.0000441139.02102.80
Godin, D. V., Garnett, M. E., Cheng, K. M., and Nichols, C. R. (1995). Sex-related alterations in antioxidant status and susceptibility to atherosclerosis in Japanese quail. Can. J. Cardiol. 11, 945–951.
Godin, D. V., Nichols, C. R., Hoekstra, K. A., Garnett, M. E., and Cheng, K. M. (2003). Alterations in aortic antioxidant components in an experimental model of atherosclerosis: a time-course study. Mol. Cell. Biochem. 252, 193–203. doi: 10.1023/A:1025548111491
Gootjes, J., Tel, R. M., Bergkamp, F. J., and Gorgels, J. P. (2009). Laboratory evaluation of a novel capillary blood sampling device for measuring eight clinical chemistry parameters and HbA1c. Clin. Chim. Acta 401, 152–157. doi: 10.1016/j.cca.2008.12.016
Huws, S. A., Kim, E. J., Lee, M. R. F., Scott, M. B., Tweed, J. K. S., Pinloche, E., et al. (2011). As yet uncultured bacteria phylogenetically classified as Prevotella, Lachnospiraceae incertae sedis and unclassified Bacteroidales, Clostridiales and Ruminococcaceae may play a predominant role in ruminal biohydrogenation. Environ. Microbiol. 13, 1500–1512. doi: 10.1111/j.1462-2920.2011.02452.x
Karlsson, F. H., Fåk, F., Nookaew, I., Tremaroli, V., Fagerberg, B., Petranovic, D., et al. (2012). Symptomatic atherosclerosis is associated with an altered gut metagenome. Nat. Commun. 3:1245. doi: 10.1038/ncomms2266
Koeth, R. A., Wang, Z., Levison, B. S., Buffa, J. A., Org, E., Sheehy, B. T., et al. (2013). Intestinal microbiota metabolism of L-carnitine, a nutrient in red meat, promotes atherosclerosis. Nat. Med. 19, 576–585. doi: 10.1038/nm.3145
Koren, O., Spor, A., Felin, J., Fåk, F., Stombaugh, J., Tremaroli, V., et al. (2011). Human oral, gut, and plaque microbiota in patients with atherosclerosis. Proc. Natl. Acad. Sci. U.S.A. 108, 4592–4598. doi: 10.1073/pnas.1011383107
Ley, R. E., Bäckhed, 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, L. (1995). Characterization and Application of a Novel Cholesterol-reducing Anaerobe, Eubacterium Coprostanoligenes ATCC 51222. Retrospective theses and dissertations, Paper 10957, Digital Repository @ Iowa State University. Available online at: http://lib.dr.iastate.edu?utm_source=lib.dr.iastate.edu%2frtd%2f6460&utm_medium=pdf&utm_campaign=pdfcoverpages/
Li, L., Baumann, C. A., Meling, D. D., Sell, J. L., and Beitz, D. C. (1996). Effect of orally administered Eubacterium coprostanoligenes ATCC 51222 on plasma cholesterol concentration in laying hens. Poult. Sci. 75, 743–745. doi: 10.3382/ps.0750743
Li, X., Schulte, P., Godin, D. V., and Cheng, K. M. (2012). Differential mRNA expression of seven genes involved in cholesterol metabolism and transport in the liver of atherosclerosis-susceptible and-resistant Japanese quail strains. Gen. Sel. Evol. 44:20. doi: 10.1186/1297-9686-44-20
Lu, J., Idris, U., Harmon, B., Hofacre, C., Maurer, J., and Lee, M. D. (2003). Diversity and succession of the intestinal bacterial community of the maturing broiler chicken. Appl. Environ. Microbiol. 69, 6816–6824. doi: 10.1128/AEM.69.11.6816-6824.2003
Maslowski, K. M., Vieira, A. T., Ng, A., Kranich, J., Sierro, F., Yu, D., et al. (2009). Regulation of inflammatory responses by gut microbiota and chemoattractant receptor GPR43. Nature 461, 1282–1286. doi: 10.1038/nature08530
Minvielle, F., Gourichon, D., Ito, S., Inoue-Murayama, M., and Rivière, S. (2007). Effects of the dominant lethal yellow mutation on reproduction, growth, feed consumption, body temperature, and body composition of the Japanese quail. Poult. Sci. 86, 1646–1650. doi: 10.1093/ps/86.8.1646
Oliveros, J. C. (2007). Venny. An Interactive Tool for Comparing Lists with Venn Diagrams. Available online at: http://bioinfogp.cnb.csic.es/tools/venny/index.html
Place, R. F., Noonan, E. J., and Giardina, C. (2005). HDAC inhibition prevents NF-κB activation by suppressing proteasome activity: down-regulation of proteasome subunit expression stabilizes IκBα. Biochem. Pharm. 70, 394–406. doi: 10.1016/j.bcp.2005.04.030
Pryde, S. E., Duncan, S. H., Hold, G. L., Stewart, C. S., and Flint, H. J. (2002). The microbiology of butyrate formation in the human colon. FEMS Microbiol. Lett. 217, 133–139. doi: 10.1111/j.1574-6968.2002.tb11467.x
Reeves, A. E., Koenigsknecht, M. J., Bergin, I. L., and Young, V. B. (2012). Suppression of Clostridium difficile in the gastrointestinal tracts of germfree mice inoculated with a murine isolate from the family Lachnospiraceae. Infect. Immunol. 80, 3786–3794. doi: 10.1128/IAI.00647-12
Sekimoto, H., Goto, Y., Goto, Y., Naito, C., Yasugi, T., Okido, M., et al. (1983). Changes of serum total cholesterol and triglyceride levels in normal subjects in Japan in the past twenty years. Jpn. Circ. J. 47, 1351–1358. doi: 10.1253/jcj.47.1351
Shih, J. C., Pullman, E., and Kao, K. (1983). Genetic selection, general characterization, and histology of atherosclerosis-susceptible and-resistant Japanese quail. Atherosclerosis 49, 41–53. doi: 10.1016/0021-9150(83)90006-0
Shornikova, A.-V., Casas, I. A., Mykkänen, H., Salo, E., and Vesikari, T. (1997). Bacteriotherapy with Lactobacillus reuteri in rotavirus gastroenteritis. Pediatr. Infect. Dis. J. 16, 1103–1107. doi: 10.1097/00006454-199712000-00002
Stepankova, R., Tonar, Z., Bartova, J., Nedorost, L., Rossman, P., Poledne, R., et al. (2010). Absence of microbiota (germ-free conditions) accelerates the atherosclerosis in ApoE-deficient mice fed standard low cholesterol diet. J. Atheroscler. Thromb. 17, 796–804. doi: 10.5551/jat.3285
Su, H., McKelvey, J., Rollins, D., Zhang, M., Brightsmith, D. J., Derr, J., et al. (2014). Cultivable bacterial microbiota of northern bobwhite (Colinus virginianus): a new reservoir of antimicrobial resistance? PLoS ONE 9:e99826. doi: 10.1371/journal.pone.0099826
Tibshirani, R., Hastie, T., Narasimhan, B., and Chu, G. (2002). Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc. Natl. Acad. Sci. U.S.A. 99, 6567–6572. doi: 10.1073/pnas.082099299
Torok, V. A., Hughes, R. J., Mikkelsen, L. L., Perez-Maldonado, R., Balding, K., MacAlpine, R., et al. (2011). Identification and characterization of potential performance-related gut microbiotas in broiler chickens across various feeding trials. Appl. Environ. Microbiol. 77, 5868–5878. doi: 10.1128/AEM.00165-11
Turnbaugh, P. J., Bäckhed, 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
Valeur, N., Engel, P., Carbajal, N., Connolly, E., and Ladefoged, K. (2004). Colonization and immunomodulation by Lactobacillus reuteri ATCC 55730 in the human gastrointestinal tract. Appl. Environ. Microbiol. 70, 1176–1181. doi: 10.1128/AEM.70.2.1176-1181.2004
Verbarg, S., Göker, M., Scheuner, C., Schumann, P., and Stackebrandt, E. (2014). “Chapter 7: The families Erysipelotrichaceae emend., Coprobacillaceae fam. nov., and Turicibacteraceae fam. nov.,” in The Prokaryotes, eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackbrandt, and F. Thompson (New York, NY: Springer), 79–105.
Videnska, P., Rahman, M. M., Faldynova, M., Babak, V., Matulova, M. E., Prukner-Radovcic, E., et al. (2014). Characterization of egg laying hen and broiler fecal microbiota in poultry farms in Croatia, Czech Republic, Hungary and Slovenia. PLoS ONE 9:e110076. doi: 10.1371/journal.pone.0110076
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
Weingand, K. W., Clarkson, T. B., Adams, M. R., and Bostrom, A. D. (1986). Effects of age and/or puberty on coronary artery atherosclerosis in cynomolgus monkeys. Atherosclerosis 62, 137–144. doi: 10.1016/0021-9150(86)90059-6
Keywords: quail model, cecal microbiome, cholesterol, atherosclerosis, genotype by diet interaction
Citation: Liu S, Bennett DC, Tun HM, Kim J-E, Cheng KM, Zhang H and Leung FC (2015) The effect of diet and host genotype on ceca microbiota of Japanese quail fed a cholesterol enriched diet. Front. Microbiol. 6:1092. doi: 10.3389/fmicb.2015.01092
Received: 01 April 2015; Accepted: 22 September 2015;
Published: 07 October 2015.
Edited by:Biswarup Mukhopadhyay, Virginia Tech, USA
Reviewed by:Sherry Layton, Vetanco, Argentina
Guillermo Tellez, University of Arkansas, USA
Maria Luisa Callegari, Università Cattolica del Sacro Cuore, Italy
Copyright © 2015 Liu, Bennett, Tun, Kim, Cheng, Zhang and Leung. 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.
*Correspondence: Kimberly M. Cheng, Faculty of Land and Food Systems, Avian Research Centre, University of British Columbia, 187-MacMillan Building, 2357 Main Mall, Vancouver, BC V6T 1Z4, Canada, email@example.com;
Frederick C. Leung, Faculty of Science, School of Biological Sciences, The University of Hong Kong, Kadoorie Biological Sciences Building, Pok Fu Lam Road, Hong Kong SAR, Hong Kong, firstname.lastname@example.org
†Present Address: Darin C. Bennett, Animal Science Department, California Polytechnic State University, San Luis Obispo, CA, USA