- Research article
- Open Access
High-throughput sequencing technology to reveal the composition and function of cecal microbiota in Dagu chicken
BMC Microbiology volume 16, Article number: 259 (2016)
The chicken gut microbiota is an important and complicated ecosystem for the host. They play an important role in converting food into nutrient and energy. The coding capacity of microbiome vastly surpasses that of the host’s genome, encoding biochemical pathways that the host has not developed. An optimal gut microbiota can increase agricultural productivity. This study aims to explore the composition and function of cecal microbiota in Dagu chicken under two feeding modes, free-range (outdoor, OD) and cage (indoor, ID) raising.
Cecal samples were collected from 24 chickens across 4 groups (12-w OD, 12-w ID, 18-w OD, and 18-w ID). We performed high-throughput sequencing of the 16S rRNA genes V4 hypervariable regions to characterize the cecal microbiota of Dagu chicken and compare the difference of cecal microbiota between free-range and cage raising chickens. It was found that 34 special operational taxonomic units (OTUs) in OD groups and 4 special OTUs in ID groups. 24 phyla were shared by the 24 samples. Bacteroidetes was the most abundant phylum with the largest proportion, followed by Firmicutes and Proteobacteria. The OD groups showed a higher proportion of Bacteroidetes (>50 %) in cecum, but a lower Firmicutes/Bacteroidetes ratio in both 12-w old (0.42, 0.62) and 18-w old groups (0.37, 0.49) compared with the ID groups. Cecal microbiota in the OD groups have higher abundance of functions involved in amino acids and glycan metabolic pathway.
The composition and function of cecal microbiota in Dagu chicken under two feeding modes, free-range and cage raising are different. The cage raising mode showed a lower proportion of Bacteroidetes in cecum, but a higher Firmicutes/Bacteroidetes ratio compared with free-range mode. Cecal microbiota in free-range mode have higher abundance of functions involved in amino acids and glycan metabolic pathway.
Chickens have proportionally smaller intestines and shorter transit digestion times than mammals, but do not appear to any less efficient at digestion than their mammalian counterparts [1, 2]. Their digestive system is adapted to extract energy from difficult to digest food sources. This may be explained, in part, by the fact that the chicken gastrointestinal tract is home to a complex microbial community, the chicken gut microbiota, which underpins the links between diet and health [3, 4]. The host is unable to digest and utilize the complicated polysaccharide substance from the feedstuff in the absence of microbial fermentation . Particularly relevant to the intensive farming of chickens is the cecum’s role in digestion of non-starch polysaccaharides NSPs , which are found in the grains used in commercial chicken feed. The gut microbiota has one of the highest cell densities for any ecosystem and ranges from 107 to 1011 bacteria per g of gut content in poultry . The most densely populated microbial community within the chicken gut is found in the ceca, a pair of blind-ended sacs that open off the large intestine . This microbiota is also home to a rich collection of genes, the chicken gut microbiome, likely to include many sequences of scientific interest and biotechnological potential . The coding capacity of microbiome vastly surpasses that of the host’s genome, encoding biochemical pathways that the host has not developed.
An optimal gut microbiota can increase agricultural productivity, as evidenced by the ability of antibiotics to promote growth in chickens . Studies on rumen microbe in ruminants have revealed that Ruminococcus and Fibrobacter species are important members of the rumen microbial community that enable the host to degrade and utilize fibrous plant materials efficiently as nutrients [10–12]. As a result, animal productivity has been improved through refining the animals’ ability to degrade fiber by these microorganisms. Energy and nutrient extraction from feed requires interplay between the biochemical functions provided by the chicken and the microbiota present within the gastrointestinal tract (GIT). Highly productive chickens have been developed by selection for elite genetic traits; it is possible that in the future, gains in productivity and health outcomes could be influenced by selection of elite GIT microbiota . Therefore, studies on the composition and functions of gut microbiota in animals raised in different feeding modes is significant for the improvement of feedstuff efficiency and animal productivity. At present, our ability to culture intestinal bacteria is limited, and hence, there is a need to profile and investigate this community using culture-independent techniques. Culture-independent analysis of the chicken cecal microbiota estimated 900 species of bacteria in 100 genera existing in the cecum of chickens, with most of them belonging to uncategorized genera [7, 14]. Previous studies have shown that the caeca microbial communities were more diverse in comparison to ilea . Left and right ceca of chickens are harbouring similar bacterial communities . But, the composition and function of cecal microbiota under different feeding mode are unknown.
Consumer interest in free-range and organic poultry is growing. The meat of the outdoor chickens had more protein than the indoor chickens . Dagu chicken is a well-known local breed in China. Dagu chicken is native to Zhuanghe City, Liaoning and is free-range. This chicken has been called Cao Chicken (cao means grass in Chinese) because of the favorable living environment and fine feed resources of water and grass. Whether living habits influence the formation of gut microflora in Dagu chicken are unknown.
This study aims to explore the composition and function of cecal microbiota in Dagu chicken under two feeding modes, free-range and cage raising. Thus, providing base informations for designing high efficiency feed formula, developing applicable probiotics and regulating chicken meat quality.
Chicken farm and sampling
Zhuanghe City is located in the south end of the Liaodong peninsula. Its location, with coordinates N39.32′–40.5′, E122′–124.5′, indicates a typical mountainous hilly terrain. Dagu chicken in free-range farming is a traditional feeding in Zhuanghe and relies on abundant rivers and flourish pasture. Thus, this study was carried out in Zhuanghe Dagu chicken breeding center.
A total of 1000 1-day-old male Dagu chickens were selected. The chickens were raised in plastic mesh floors (80 cm above ground) for 6 weeks. The chickens were provided access to feed and water ad libitum. The house temperature was maintained at 35 °C during the first week, and it was reduced 2 °C per week until reaching the temperature of 23 °C. Six weeks later, 300 chickens with similar weights were randomly selected. Among them, 150 chickens were raised outside, which are in the outdoor group (OD group), while the other 150 chickens were raised inside their respective cages (50 cm × 50 cm × 50 cm,80 cm above ground) which are in the indoor group (ID group). The house temperature was maintained at 23 °C. The chickens were provided access to feed and water ad libitum. The difference of body weight between the two groups was not significant (P >0.05). Two groups were given the same compound feed (Additional file 1) as well as other environmental factors. The difference is for the OD group, each chicken in the OD group was let out every 5 am for self-help feeding in >30 m2 area, where abundant water and grasses are found. The chickens were given supplementary feed at 1 pm, and kept indoors at 7 pm. When the chickens were 12 weeks old and 18 weeks old, weighed one by one, six of them with an average weight were randomly selected in each group,and then slaughtered. The cecum contents removed, preserved in liquid nitrogen, used for DNA extraction and PCR amplification. These samples were divided into four groups, namely, 12-w OD group, 12-w ID group, 18-w OD group, and 18-w ID group.
Gut microbes 16S rRNA sequencing
Microbial genomic DNA was extracted from cecal content samples by using the TIANGEN DNA stool mini kit (TIANGEN, cat#DP328) according to the producer’s instructions (http://www.tiangen.com/asset/imsupload/up0921879001368428871.pdf). Variable region of 16S rRNA V4 was amplified using its universal primer sequence 520 F: AYTGGGYDTAAAGNG; 802R: TACNVGGGTATCTAATCC . The PCR conditions were as follows: initial denaturation at 98 °C for 5 min; 98 °C denaturation for 30 s, 50 °C annealing for 30 s, and 72 °C extension for 30 s, which is repeated for 28 cycles; and a final extension at 72 °C for 5 min. PCR production was purified using QIAGEN Quick Gel Extraction Kit (QIAGEN, cat# 28706). PCR production from each sample was applied to construct a sequencing library by using Illumina TruSeq DNA Sample Preparation Kit (library was constructed using TruSeq Library Construction Kit). For each sample, barcoded V4 PCR amplicons were sequenced by the Illumina MiSeq PE250 platform.
Sequence reads were removed if sequence length was shorter than 150 bp, if average phred score was lower than 20, if containing ambiguous bases, if homopolymer run exceeded 6, or if there were mismatches in primers. Afterward, the sequences passed the quality filter that were assembled by Flash (http://www.genomics.jhu.edu), which required that the overlap of read 1 and read 2 ≥ 10 bp, and without any mismatches. The reads which could not be assembled were discarded. Chimera sequences were removed using UCHIME in mothur (version 1.31.2, http://www.mothur.org/). Amplification and sequencing of 16S rRNA v4 variable region was completed by Personal Biotechnology Co., Ltd. (Shanghai, China).
OTU clustering and statistical analysis
Sequences clustering was performed by uclust algorithm in QIIME (http://qiime.org/scripts/pick_otus.html), and clustered into operational taxonomic units (OTUs). The longest sequence in each cluster was selected as the representative. Taxonomy of each OTU was assigned by blasting the representative sequence against Greengenes reference database (Release 13.8, http://greengenes.secondgenome.com/). Unknown archaeal or eukaryotic sequences were filtered and removed. Ace, Chao, Simpson index were calculated using summary.single command in MOTHUR. A Venn diagram of between-group OTU was generated through R. The relative abundance of OTUs or taxa was compared between samples.
Diversity index data were analyzed statistically using analysis of variance (ANOVA) and significant differences between group means were determined using the least significant difference (LSD) test. Data of body weight and abundance at the phylum level between groups were analyzed statistically using T test. All values for diversity index and body weight are expressed as means ± standard errors (SE). Non-metric multidimensional scaling (NMDS) plots of sequence read abundance were generated with Vegan in R. All statistical analyses were performed using the SPSS 16.0 software.
Microbial function prediction
Functional genes were predicted through PICRUSt according to the abundance of OTU level . The OTUs were mapped in gg13.5 database at 97 % similarity by QIIME’s command “pick_closed_otus”. The abundance of the OTUs was normalized automatically by using 16S rRNA gene copy numbers from known bacterial genomes in the Integrated Microbial Genomes (IMG). The predicted genes and their function were aligned to the Kyoto Encyclopedia of Genes and Genomes (KEGG) database and the differences among groups were compared using STAMP (http://kiwi.cs.dal.ca/Software/STAMP) . Two-side Welch’s t-test and Benjamini–Hochberg FDR correction were employed in the two-group analysis. The relative abundance of KEGG metabolic pathways is referred to as a metabolic profile.
OTU clustering and annotation
The trimmed and assembled sequences were clustered at 97 % similarity by calling uclust from Qiime. 1217 OTUs were obtained through database alignment by blast in Qiime. The total of OTUs obtained in each group were as follows: 1188 in the 12-w OD group, 1089 in the 12-w ID group, 1186 in the 18-w OD group, and 1158 in the 18-w ID group (Fig. 1). Figure 1 shows 34 special OTUs in OD groups (including 12-w OD and 18-w OD) and 4 special OTUs in ID groups (including 12-w OD and 18-w OD). The number of OTUs in each group slightly changed in the OD groups, whereas that increased in the ID groups within days. The diversity of cecal microbiota in OD groups can be established earlier. The Chao and ACE in the 12-w OD group were significantly higher (P < 0.05) than those in the three other groups, but the Simpson in the OD groups was significantly lower (P < 0.05) than that in the ID groups. These results revealed that the richness of cecum microorganism in the 12-w OD group was higher than those in the three other groups, the evenness of cecum microorganism in the ID groups was higher than those in the OD groups (Table 1).
Differences of body weight and cecal microbiota in chickens raised in different feeding modes
In this study, chicken body weight in different feeding modes has obvious differences. Chicken body weight in the ID group was significantly higher than that in the OD group both 12-w or 18-w stage (Table 2).
A total of 24 phyla were shared by the 24 samples. Bacteroidetes was the most abundant phylum with the largest proportion, followed by Firmicutes and Proteobacteria (Fig. 2). Three significant differences (P < 0.05) in the 12-w groups and five significant (P < 0.05) differences in the 18-w groups were found (Table 3).
Spirochaetes had dynamic changes in the ID groups; its proportion was 5.73 % in the 12-w ID group, but it reduced to 1.4 % in the 18-w ID group. However, the proportion had a slight change in the OD groups (3.6 %, 3.05 %).
In the 12-w groups, Bacteroidetes, Firmicutes and Proteobacteria accounted for 83 % and 53 %, and Bacteroidetes for 52 % and 26.7 % in the OD and ID groups, respectively. In the 18-w groups, the three phyla accounted for 84.9 % and Bacteroidetes accounted for 53.66 % in the OD group, which exhibited a slight difference from that in the 12-w OD group. The proportion was 60.5 %, and Bacteroidetes had the largest share of 35 % in the ID group.
SAR406 mainly existed in the ID groups, accounting for 22.1 % in the 12-w group and 15.54 % in the 18-w group. SAR406 accounted for 0.17 % in the 12-w group and 0.37 % in the 18-w group in the OD groups, respectively.
At the genus level we detected 60 genera. 10 genera were significantly different (P < 0.01) between the 12-w OD and 12-w ID groups (Additional file 2), 6 genera were significantly different (P < 0.01) between the 18-w OD and 18-w ID groups (Additional file 3).
NMDS results showed the difference in microorganism distributions in the four groups. The distribution was evidently different in the OD groups compared with that in the ID groups (Fig. 3). The microorganisms in the OD groups concentrated on one group whereas those in the ID groups concentrated on another. Numerical values in correlation analysis revealed that the cecal microbiota in the 12-w OD groups were quite different from those in the 12-w ID group (0.5729). However, the cecal microbiota in the 18-w OD group were remarkably similar to those in the 18-w ID group (0.9626) (Table 4). The results show that the richness and evenness of cecal microbiota in chickens raised in cages were noticeably different from those in chickens from free-range farming, especially at 12 weeks.
Microbial function analysis through PICTUSt was conducted to determine the differences in the functions of microbiota between the OD and ID groups. Numerous functions are involved in metabolic pathways. At KEGG level 2, cecal microbiota in the OD groups have higher abundance of functions involved in amino acids metabolic pathway (Fig. 4). At KEGG level 3, cecal microbiota in the 12-w OD group have higher abundance of functions involved in metabolic pathway such as metabolism of arginine, praline, histidine, glycine, serine, threonine, alanine, aspartate and glutamate, starch and sucrose, galactose, amino sugar and nucleotide sugar, and transcription machinery, DNA replication proteins than those in the 12-w ID group. Cecal microbiota in the 18-w OD group have higher abundance of functions involved in metabolic pathway such as metabolism of glycine, serine, threonine, arginine, praline, tryptophan, phenylalanine, tyrosine, and valine, leucine and isoleucine biosynthesis, amino acid related enzymes than those in the 18-w ID group. In the OD groups, cecum contained more microbiota associated with glycosaminoglycan degradation and other glycan degradation (Additional file 4).
Digestion and nutrient absorption are the basic function of the intestine, where gut microbiota play an important role. These microbiota have a significant influence on intestinal tract movement, growth and development, physiological functions, and non-specific immunity [20–25]. The diversity of gut microorganism is the foundation for animals’ digestion and nutrient uptake, maintenance of biochemical functions and the intestine’s physiological functions, and promotion of the immune system’s development. Medical researches discovered that obesity is related to the changes of gut microbiota, diversity of gut microbiota apparently decreases in obese patients . The results of this study show that the body weight of caged chicken was significantly higher than that of free-range groups (Table 2). This is consistent with the results of other studies [27, 28]. Figure 1 shows 34 special OTUs in OD groups and 4 special OTUs in ID groups. The diversity of cecal microbiota in the OD groups was remarkably higher than that in the ID groups (Table 1). Bailey et al. discovered that long-term stress could reduce the diversity of gut microbiota in mice . Chickens raised in OD and ID groups were exposed to distinct stresses and microbiota. Chickens raised in ID groups were exposed to more stresses, such as feeding density and space . Chickens raised in OD groups may be due to the earlier contact to the natural environment; thus, the diversity can be established earlier.
Host and environmental factors influence the gut microbiota. The environmental factor is more important than the host factor [31–34]. The phyla Firmicutes and Bacteroidetes dominate the intestine of mammals, followed by Fusobacteria, Proteobacteria, and Actinobacteria . Bacteroidetes and Firmicutes have attracted considerable attention and are bounded to the host’s metabolism. Numerous studies investigated the probiotic effect of Bacteroidetes; they found that Bacteroidetes help the host in polysaccharide decomposition to improve nutrient utilization , promote immune system development, improve host’s immunity [37, 38], and maintain intestinal microecological balance [39, 40]. Results in this paper show that more (>50 %) Bacteroidetes existed in chickens in the OD groups (Fig. 2). and that the Firmicutes/Bacteroidetes ratio was smaller in the OD groups, with 0.42 and 0.62 in the 12-w groups and 0.37 and 0.49 in the 18-w groups. Research has shown that adding more dietary fiber can increase the amount of Bacteroidetes and lower the Firmicutes/Bacteroidetes ratio [41, 42]. The results show that compositions of cecal microbiota in chickens raised in two feeding modes were apparently different (Table 3, Fig. 3), especially at 12 weeks (Table 4). The difference may be attributed to the access of chickens from free-range farming to abundant microbiota in the outdoor environment; these chickens have abundance of food source and are able to intake more feedstuff containing fiber, which directly affects the composition of gut microbiota, increasing the Bacteroidetes content and lowering the Firmicutes/Bacteroidetes ratio.
Obesity is related to the distribution of gut bacteria. High ratio of Firmicutes/Bacteroidetes causes obesity because more energy has been absorbed . The small intestine is mainly involved in digestion and uptake of food, while a large amount of microorganisms related to microbial fermentation exists in the large intestine, especially the cecum . Food rapidly passes the front of the intestinal tract but stays for several hours in the tail end of the tract . Fat deposits mainly in the large intestine , which is closely related to the composition of microorganisms. In chicken production, bacteria related to productivity mainly include the phylum Firmicutes, along with Bacteroidetes and Proteobacteria . Researchers suggested that fat pigs have more Firmicutes but fewer Bacteroidetes, especially fewer Bacteroides that are crucial in carbohydrate degradation [48, 49]. A study revealed that free-range farming can evidently reduce the growth performance and abdominal fat of chickens . However, the efficiency of converting feedstuff to energy together with the chickens’ productivity attracts increasing attention in the chickens production. In this paper, body weight of caged chickens was significantly higher than that of free-rage chickens (Table 2). We speculate that this may be due to that more Firmicutes and higher ratio of Firmicutes/Bacteroidetes in cecal microbiota improve the utilization efficiency of feed energy, of course, this needs further study.
Gut microbiota contains about 600,000 genes that are 25 times more compared with the genes in host’s genome. Therefore, gut microbiota is usually regarded as one organ of the host and creates a gut microecosystem with the host’s eucells [50, 51]. This microecosystem can execute numerous metabolic functions that alter with the change of microbiota’ composition. In this paper, numerous functions are involved in metabolic pathways, such as metabolism of amino acid, carbohydrates, energy, lipid, replication and repair, nucleotides, and cofactors and vitamins. At KEGG level 2, there are 5 significant differences (P < 0.05) in abundance of functional categories between OD and ID group at 12-w, whereas 7 significant differences (P < 0.05) were found in between at 18-w (Fig. 2). At KEGG level 3, there are 42 significant differences (P < 0.05) in abundance of functional categories between OD and ID group at 12-w, among them 34 in OD group was significantly higher than that in ID group (P < 0.05). There are 72 significant differences (P < 0.05) in abundance of functional categories between the OD and ID group at 18-w, among them 44 in OD group was significantly higher than that in ID group (P < 0.05) (Additional file 4).
Cecal microbiota of OD group at 12-w and 18-w both has higher abundance of functions involved in metabolic pathway for certain amino acids, sugar compounds. Significant difference in amino sugar and nucleotide sugar metabolism pathways were observed in 12-w groups. Utilization of amino sugar and nucleotide sugar is important in chicken metabolism and growth. Amino sugar metabolism specifically is responsible for breaking down protein present in feed to amino acids or di- or tri-peptides . These were then transported from intestinal lumen to epithelial cell for energy. Nucleotide sugar metabolism on the other hand is important for purine and pyrimidine synthesis which is vital substrate for deoxyribonucleic acids derivatives. In addition, these components are also needed for producing high-energy nucleotides needed for cellular metabolism . In this study, we observed that the genes responsible for amino sugar and nucleotide sugar metabolism were up-regulated in 12-w OD group compared to 12-w ID group (P < 0.05) (Additional file 4). This may be the outdoor chickens needs more energy due to the large amount of movement. And movement promotes muscle development, and therefore the synthesis of more body protein. In contrast, the genes related to the metabolism of amino acids, amino sugars and nucleotide sugar were up-regulated in the cecum. Previous studies show that outdoor activities could make an improvement on the meat quality. The meat of chickens with outdoor access is darker, it has more protein contents and a better water-holding capacity [54, 55]. In addition, studies have revealed that feeding chickens with probiotics can improve meat quality and increase the output of breast and leg muscles . All of these are likely to be related to the changes in compositions of gut microbiota. But, more scientific research is needed to confirm this.
Based on the research above, the many metabolic functions are involved in chickens’ gut microbiota and these functions may vary because of the different compositions of gut microbiota. The compositions of chickens’ cecal microbiota varied because the chickens were raised in different feeding modes. In-depth studies on the functions of dominant gut microbiota, such as Bacteroidetes and Firmicutes and their interaction, can help us develop a special probiotics and guide us to use the special probiotics to achieve the anticipated breed goals.
The composition and function of cecal microbiota in Dagu chicken under two feeding modes, free-range and cage raising are different. The cage raising mode showed a lower proportion of Bacteroidetes in cecum, but a higher Firmicutes/Bacteroidetes ratio compared with free-range mode. Cecal microbiota in free-range mode have higher abundance of functions involved in amino acids and glycan metabolic pathway. The results in this paper can provide relevant information for making strategies in raising Dagu chickens. This also provided valuable information for the study on microbiota in chicken gut.
McWhorter TJ, Caviedes-Vidal E, Karasov WH. The integration of digestion and osmoregulation in the avian gut. Biol Rev Camb Philos Soc. 2009;84(4):533–65.
Stanley D, Geier MS, Chen H, Hughes RJ, Moore RJ. Comparison of fecal and cecal microbiotas reveals qualitative similarities but quantitative differences. BMC Microbiol. 2015;15:51.
Björkstén B. The gut microbiota: a complex ecosystem. Clin Exp Allergy. 2006;36(10):1215–17.
Sergeant MJ, Constantinidou C, Cogan TA, Bedford MR, Penn CW, Pallen MJ. Extensive Microbial and Functional Diversity within the Chicken Cecal Microbiome. PLoS One. 2014;9(3):e91941. doi:10.1371/journal.pone.0091941.
Vrieze A, Holleman F, Zoetendal EG, de Vos WM, Hoekstra JB, Nieuwdorp M. The environment within: how gut microbiota may influence metabolism and body composition. Diabetologia. 2010;53(4):606–13. doi:10.1007/s00125-010-1662-7.
Jozefiak D, Rutkowski A, Martin SA. Carbohydrate fermentation in the avian ceca: a review. Anim Feed Sci Technol. 2004;113:1–15.
Apajalahti J, Kettunen A, Graham H. Characteristics of the gastrointestinal microbial communities, with special reference to the chicken. Worlds Poult Sci J. 2004;60:223–32. doi:10.1079/WPS200415.
Clench MH, Mathias JR. The avian cecum: a review. Wilson Bull. 1995;107(l):93–121.
Huyghebaert G, Ducatelle R, Van Immerseel F. An update on alternatives to antimicrobial growth promoters for broilers. Vet J. 2011;187:182–8.
Mackie RI, White BA. Recent advances in rumen microbial ecology and metabolism: potential impact on nutrient output. J Dairy Sci. 1990;73(10):2971–95. PMID:2178174.
Perumbakkam S, Mitchell EA, Craig AM. Changes to the rumen bacterial population of sheep with the addition of 2,4,6-trinitrotoluene to their diet. Antonie Van Leeuwenhoek. 2011;99(2):231–40. doi:10.1007/s10482-010-9481-x. PMID:20607404.
Weimer PJ. Redundancy, resilience, and host specificity of the ruminal microbiota: implications for engineering improved ruminal fermentations. Front Microbiol. 2015;6:296. doi:10.3389/fmicb.2015.00296. PMID:25914693.
Stanley D, Hughes RJ, Moore RJ. Microbiota of the chicken gastrointestinal tract: influence on health, productivity and disease. Appl Microbiol Biotechnol. 2014;98(10):4301–10.
Wei S, Morrison M, Yu Z. Bacterial census of poultry intestinal microbiome. Poult Sci. 2013;92:671–83.
Mohd Shaufi MA, Sieo CC, Chong CW, Gan HM, Ho YW. Deciphering chicken gut microbial dynamics based on high-throughput 16S rRNA metagenomics analyses. Gut Pathog. 2015;7:4. doi:10.1186/s13099-015-0051-7.
Fanatico AC, Pillai PB, Emmert JL, Owens CM. Meat quality of slow- and fast-growing chicken genotypes fed low-nutrient or standard diets and raised indoors or with outdoor access. Poult Sci. 2007;86(10):2245–55.
Blanton LV, Charbonneau MR, Salih T, Barratt MJ, Venkatesh S, Ilkaveya O, et al. Gut bacteria that prevent growth impairments transmitted by microbiota from malnourished children. Science. 2016; doi:10.1126/science.aad3311.
Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31(9):814–21. doi:10.1038/nbt.2676.
Parks DH, Tyson GW, Hugenholtz P, Beiko RG. STAMP: Statistical analysis of taxonomic and functional profiles. Bioinformatics. 2014;30(21):3123–4. doi:10.1093/bioinformatics/btu494.
De Angelis M, Siragusa S, Berloco M, Caputo L, Settanni L, Alfonsi G, et al. Selection of potential probiotic lactobacilli from pig feces to be used as additives in pelleted feeding. Res Microbiol. 2006;157(8):792–801.
Forder RE, Howarth GS, Tivey DR, Hughes RJ. Bacterial modulation of small intestinal goblet cells and mucin compositionduring early posthatch development of poultry. Poult Sci. 2007;86(11):2396–403.
Klasing KC. Nutrition and the immune system. Br Poult Sci. 2007;48(5):525–37.
Niba AT, Beal JD, Kudi AC, Brooks PH. Bacterial fermentation in the gastrointestinal tract of non-ruminants: influence of fermented feeds and fermentable carbohydrates. Trop Anim Health Prod. 2009;41(7):1393–407. doi:10.1007/s11250-009-9327-6.
Shira EB, Sklan D, Friedman A. Impaired immune responses in broiler hatchling hindgut following delayed access to feed. Vet Immunol Immunopathol. 2005;105(1-2):33–45.
Yin Y, Lei F, Zhu L, Li S, Wu Z, Zhang R, et al. Exposure of different bacterial inocula to newborn chicken affects gut microbiota development and ileum gene expression. ISME J. 2010;4(3):367–76. doi:10.1038/ismej.2009.128.
Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457(7228):480–4.
Wang KH, Shi SR, Dou TC, Sun HJ. Effect of a free-range raising system on growth performance, carcass yield, and meat quality of slow-growing chicken. Poult Sci. 2009;88(10):2219–23. doi:10.3382/ps.2008-00423.
Dou TC, Shi SR, Sun HJ, Wang KH. Growth rate, carcass traits and meat quality of slow-growing chicken grown according to three raising systems. Anim Sci Pap Rep. 2009;27:361–9.
Bailey MT, Dowd SE, Parry NM, Galley JD, Schauer DB, Lyte M. Stressor exposure disrupts commensal microbial populations in the intestines and leads to increased colonization by Citrobacter rodentium. Infect Immun. 2010;78(4):1509–19. doi:10.1128/IAI.00862-09.
Jones M, Millis AD. Divergent selection for social reinstatement and behaviors in Japanese quail: effects on sociality and social discrimination. Poult Avian Biol Rev. 1999;10(4):213–23.
Carmody RN, Gerber GK, Luevano Jr JM, Gatti DM, Somes L, Svenson KL. Diet dominates host genotype in shaping the murine gut microbiota. Cell Host Microbe. 2015;17(1):72–84. doi:10.1016/j.chom.2014.11.010.
Turnbaugh PJ, Ridaura VK, Faith JJ, Rey FE, Knight R, Gordon JI. The effect of diet on the human gut microbiome: a metagenomic analysis in humanized gnotobiotic mice. Sci Transl Med. 2009;1(6):6ra14. doi:10.1126/scitranslmed.3000322.
Tachon S, Zhou J, Keenan M, Martin R, Marco ML. The intestinal microbiota in aged mice is modulated by dietary resistant starch and correlated with improvements in host responses. FEMS Microbiol Ecol. 2013;83(2):299–309. doi:10.1111/j.1574-6941.2012.01475.x.
Fava F, Gitau R, Griffin BA, Gibson GR, Tuohy KM, Lovegrove JA. The type and quantity of dietary fat and carbohydrate alter faecal microbiome and short-chain fatty acid excretion in a metabolic syndrome ‘at-risk’ population. Int J Obes (Lond). 2013;37(2):216–23. doi:10.1038/ijo.2012.33.
Ley RE, Hamady M, Lozupone C, Turnbaugh PJ, Ramey RR, Bircher JS, et al. Evolution of mammals and their gut microbes. Science. 2008;320(5883):1647–51. doi:10.1126/science.1155725.
Bäckhed F, Ding H, Wang T, Hooper LV, Koh GY, Nagy A, et al. The gut microbiota as an environmental factor that regulates fat storage. Proc Natl Acad Sci. 2004;101(44):15718–23.
Stappenbeck TS, Hooper LV, Gordon JI. Developmental regulation of intestinal angiogenesis by indigenous microbes via Paneth cells. Proc Natl Acad Sci. 2002;99(24):15451–5.
Hooper LV. Bacterial contributions to mammalian gut development. Trends Microbiol. 2004;12(3):129–34.
Sears CL. A dynamic partnership: celebrating our gut flora. Anaerobe. 2005;11(5):247–51.
Hooper LV, Wong MH, Thelin A, Hansson L, Falk PG, Gordon JI. Molecular analysis of commensal host-microbial relationships in the intestine. Science. 2001;291(5505):881–4. PMID:11157169.
Trompette A, Gollwitzer ES, Yadava K, Sichelstiel AK, Sprenger N, Ngom-Bru C, et al. Gut microbiota metabolism of dietary fiber influences allergic airway disease and hematopoiesis. Nat Med. 2014;20(2):159–66.
Parnell JA, Reimer RA. Prebiotic fiber modulation of the gut microbiota improves risk factors for obesity and the metabolic syndrome. Gut Microbes. 2012;3(1):29–34.
Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V, Mardis ER, Gordon JI. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. 2006;444(7122):1027–31.
DiBaise JK, Zhang H, Crowell MD, Krajmalnik-Brown R, Decker GA, Rittmann BE. Gut microbiota and its possible relationship with obesity. Mayo Clin Proc. 2008;83(4):460–9. doi:10.4065/83.4.460.
Kohl KD, Miller AW, Marvin JE, Mackie R, Dearing MD. Herbivorous rodents (Neotoma spp.) harbour abundant and active foregut microbiota. Environ Microbiol. 2014;16(9):2869–78. doi:10.1111/1462-2920.12376.
Zhao W, Wang Y, Liu S, Huang J, Zhai Z, He C, et al. The Dynamic distribution of porcine microbiota across different ages and gastrointestinal tract segments. PLoS One. 2015;10(2):e0117441. doi:10.1371/journal.pone.0117441.
Torok VA, Hughes RJ, Mikkelsen LL, Perez-Maldonado R, Balding K, MacAlpine R, et al. Identification and characterization of potential performance-related gut microbiotas in broiler chickens across various feeding trials. Appl Environ Microbiol. 2011;77(17):5868–78. doi:10.1128/AEM.00165-11.
Ley RE, Turnbaugh PJ, Klein S, Gordon JI. Microbial ecology: human gut microbes associated with obesity. Nature. 2006;444(7122):1022–3.
Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, et al. Enterotypes of the human gut microbiome. Nature. 2011;473(7346):174–80. doi:10.1038/nature09944.
Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, et al. A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010;464(7285):59–65. doi:10.1038/nature08821.
Lederberg J. Infectious history. Science. 2000;288(5464):287–93.
Miska KB, Fetterer RH, Wong EA. The mRNA expression of amino acid transporters, aminopeptidase N, and the di- and tri-peptide transporter PepT1 in the embryo of the domesticated chicken (Gallus gallus) shows developmental regulation. Poult Sci. 2014;93:2262–70.
Rengaraj D, Lee BR, Jang HJ, Kim YM, Han JY. Comparative metabolic pathway analysis with special reference to nucleotide metabolism-related genes in chicken primordial germ cells. Theriogenology. 2013;79:28–39.
Mikulski D, Celej J, Jankowski J, Majewska T, Mikulska M. Growth performance, carcass traits and meat quality of slower-growing and fast-growing chickens raised with and without outdoor access. Asian Aust J Anim Sci. 2011;24(10):1407–16. doi:10.5713/ajas.2011.11038.
Chen X, Jiang W, Tan HZ, Xu GF, Zhang XB, Wei S, et al. Effects of outdoor access on growth performance, carcass composition, and meat characteristics of broiler chickens. Poult Sci. 2013;92(2):435–43. doi:10.3382/ps.2012-02360. PMID:23300311.
Kabir SML, Rahman MM, Rahman MB, Rahman MM, Ahmed SU. The dynamics of probiotics on growth performance and immune response in broilers. Int J Poult Sci. 2004;3(5):361–4.
We thank T.-Y. Wang for help with sample collection. We also thank Y. Liu for technical assistance. All authors participated in the writing of the manuscript and agreed with the final format.
This study was funded by the National Natural Science Foundation of China (Grant no. 31172242; http://www.nsfc.gov.cn/) and Science and Technology Innovation Team Project of Liaoning Medical University (Grant no. 30420130104; http://www.lnmu.edu.cn/).
Availability of data and material
Sequences of this project have been deposited in the NCBI sequence read archive (http://www.ncbi.nlm.nih.gov/Traces/sra/) under SRA accession number: SRP078216.
YT, YS and YX carried out the whole experiment and wrote the paper. LZ, HY, HX and DS participated in the sample collecting and analysis. All authors read and approved the final manuscript.
The authors declare that they have no competing interest.
Consent for publication
Ethics approval and consent to participate
Protocols in this experiment were in accordance with the Guidelines for the Care and Use of Laboratory Animals established by the Beijing Association for Laboratory Animal Science and approved by the Animal Ethics Committee of the Institute of Zoology, Chinese Academy of Sciences. Client-owner consent was obtained to collect the samples from the free-range and cage raising chickens.
Diet ingredients. (DOC 38 kb)
Comparisons for abundance between the 12-w OD and 12-w ID groups at the genus level. (XLS 19 kb)
Comparisons for abundance between the 18-w OD and 18-w ID groups at the genus level. (XLS 19 kb)
Mean proportion and their differences in predicted functional metagenomes of the cecal microbiota at KEGG level 3. Comparison of functional pathway between microbes of 12-w OD group and 12-w ID group (a), 18-w OD group and 18-w ID group (b). (TIF 3907 kb)