- Research article
- Open Access
Development of the duodenal, ileal, jejunal and caecal microbiota in chickens
Animal Microbiome volume 1, Article number: 17 (2019)
The chicken intestinal microbiota plays a large role in chicken health and productivity and a greater understanding of its development may lead to interventions to improve chicken nutrition, disease resistance and welfare.
In this study we examine the duodenal, jejunal, ileal and caecal microbiota of chickens from day of hatch to 5 weeks of age (day 1, 3, 7, 14 and week 5). DNA was extracted from intestinal content samples and the V4 region of the 16S rRNA gene was amplified and sequenced. We identified significant differences in microbial community composition, diversity and richness between samples taken from different locations within the chicken intestinal tract. We also characterised the development of the microbiota at each intestinal site over time.
Our study builds upon existing literature to further characterise the development of the chicken intestinal microbiota.
Improvements in sequencing technologies have led to a better understanding of the microbiota of many livestock species [1,2,3,4], leading some to suggest that we could optimise the composition of microbiota in these economically important animals to improve production and sustainability . However, to do so would require a good understanding of the types of microbes which naturally occur in these animals and the role they play in nutrition and health.
Many studies have used 16S rRNA gene data to characterise the microbial communities which colonise the gastrointestinal tracts of chickens and to characterise the development of these communities over time. The vast majority of these studies have focussed on the chicken caeca as this is where the largest concentration of microbes can be found. The caecal microbiota has been suggested to play an important role in nutrition via the production of short chain fatty acids, nitrogen recycling and amino acid production [6,7,8]. In early life it is generally observed that the caeca contain high abundances of Enterobacteriales [9, 10] and over the first few weeks of life these decline and members of the Clostridiales come to predominate [9,10,11,12,13,14,15,16], with some studies also showing a large increase in Bacteroidetes [12, 17, 18]. However, the results from some studies do not entirely follow this pattern [19, 20] and variability in microbiota composition between flocks can be high [21, 22]. Several studies have also examined samples from the small intestine which are less rich and diverse than caecal samples and contain a high abundance of Lactobacilli [12, 14, 17, 18, 20, 23,24,25].
Several studies have directly compared samples taken from the small intestine with those from the caeca at specific timepoints [17, 18], and at various life stages. Lu et al. compared ileal and caecal samples from commercial Ross-hybrids, taking samples at day 1, 3, 7 14, 21, 28, and 49 days of age . A separate study examined ileal and caecal samples from Hy-Line W-36 commercial layers at 9 timepoints, starting at 1 week of age . Ileal and caecal samples were also characterised at various timepoints starting at 1 week of age by Johnson et al. in Cobb 500 birds .
In this study we compared the bacterial microbiota of duodenal, jejunal, ileal and caecal samples of Ross 308 broilers at 5 timepoints: 1 day, 3 days, 7 days, 14 days and 5 weeks of age. Similar to the studies above, we found significant differences in the microbiota compositions of specific sample types with age and significant differences between sample types within timepoints.
The V4 region of the 16S rRNA gene was amplified and sequenced, producing a total of 11,115,696 paired-end reads from 164 samples. 59.9% of sequences were removed during quality control. For each sample the average number of reads after quality control was 27,184 ± 44,159 (mean ± standard deviation). A total of 6015 operational taxonomic units (OTUs) were identified. Additional file 1 contains the complete OTU table and Additional file 2 contains the taxonomic assignment of OTUs. The most abundant bacteria in our negative controls were Clostridium_sensu_stricto_1 (0.23 ± 0.20), Lactobacillus (0.17 ± 0.12), unclassified members of the Lachnospiraceae (0.14 ± 0.19) and Enterococcus (0.10 ± 0.09). The composition of the mock community control can be found in Additional file 3: Tables S1 and S2. Prior to analysis, samples were subsampled to 10,000 reads; 35 samples were discarded as they had less than 10,000 reads, including all reagent only controls. This also resulted in only one duodenum sample from the day 1 timepoint remaining. The lowest Good’s coverage value for any of the remaining samples was 0.985, meaning that at least 98.5% of the bacteria in these samples were identified. Samples were compared statistically to see if there were differences in the microbiota between sample types at specific timepoints and if there were changes in specific sample types at different timepoints (Figs. 1 and 2). The number of OTUs shared between sample types within timepoints and between timepoints within sample types can be found in Additional file 3: Figures S1 and S2. Richness and diversity measurements for each subsampled sample can be found in Additional file 3: Table S3.
It is clear that time is an important factor affecting the intestinal microbiota composition, with the day 1 timepoint clustering particularly separately from other timepoints (Fig. 2). It is also clear that position in the digestive tract is important, with caecal samples clearly separating from other all other parts of the digestive tract throughout development, except at day 14 (Fig. 1).
At day 1 in jejunal, ileal and caecal samples the most abundant genus on average was Clostridium_sensu_stricto_1 (jejunum: 0.76 ± 0.32, ileum: 0.86 ± 0.22, caeca 0.95 ± 0.13) (Fig. 3). Sample types did cluster significantly separately from one another by their microbiota compositions (Permutational multivariate analysis of variance (PERMANOVA): P = 0.00478) (Fig. 1). At day 1, using the Pairwise Wilcoxon Rank Sum Test with Bonferroni corrections, microbial communities from jejunal and caecal samples were found to be significantly differently rich (P = 0.0018) and diverse (P = 0.033), with jejunum samples having greater diversity and richness (Figs. 4 and 5). Several OTUs were also found to be more abundant in the small intestinal samples than in caecal samples (Additional file 4).
Between day 1 and day 3, a significant difference was observed in the microbiota composition of all sample types (jejunum: P = 0.01, ileum: P = 0.03, caeca: P = 0.01) except duodenal, due to the lack of duodenal samples at the day 1 timepoint. In caecal samples this correlates with a decrease in the abundance of several OTUs including Otu0001 (Clostridium_sensu_stricto_1) and an increase in several OTUs belonging to the order Clostridiales (Additional file 5). This significant decrease in Otu0001 was also observed in duodenal, jejunal and ileal samples (Additional files 6, 7 and 8). Significant increases were found in caecal microbiota richness and ileal and caecal microbiota diversity between these timepoints (Richness - caeca: P = 0.0011; Diversity - ileum: P = 0.023 caeca: P = 0.0011) but not for duodenal or jejunal samples.
At day 3, caecal samples clustered significantly separately by PERMANOVA from all other sample types (duodenum: P = 0.018, jejunum: P = 0.006, ileum: P = 0.024) (Fig. 1). Duodenal, jejunal and ileal samples did not cluster significantly separately; however, a small number of OTUs were found to be differently abundant (Additional file 9). The most common genera found in caecal samples were Enterobacteriaceae_unclassified (0.24 ± 0.010), Lachnospiraceae_unclassified (0.19 ± 0.13), Enterococcus (0.11 ± 0.058), Escherichia-Shigella (0.080 ± 0.087) and Lactobacillus (0.064 ± 0.033). All three types of small intestinal samples were dominated by Lactobacilli and Enterococci. Caecal samples at day 3 were significantly more diverse than other sample types (duodenum: P = 0.013, jejunum: P = 0.013, ileum: P = 0.013) but were only significantly more rich than ileal (P = 0.026) and duodenal (P = 0.013) samples. There was a significant difference in richness but not diversity between duodenal samples and jejunal (richness: P = 0.013) and ileal samples (richness: P = 0.013). Between day 3 and day 7, no significant differences could be observed in the microbiota composition (PERMANOVA), diversity or richness of any sample type except for an increase in richness in duodenum samples (P = 0.043) and an increase in diversity in caecal samples (P = 0.0433). However, significant changes in the abundance of specific OTUs were observed (Additional files 5, 6, 7 and 8).
At day 7, caecal samples still clustered significantly separately by PERMANOVA from all other sample types (duodenum: P = 0.024, jejunum: P = 0.018, ileum: P = 0.024) (Fig. 1, Additional file 10). Duodenal, jejunal and ileal samples still did not cluster significantly separately from each other by PERMANOVA. The most common genera found in caecal samples were Lachnospiraceae_unclassified (0.34 ± 0.031), Lactobacillus (0.117 ± 0.073), Escherichia-Shigella (0.10 ± 0.052), Lachnoclostridium (0.071 ± 0.031) and Erysipelatoclostridium (0.063 ± 0.035). Duodenal, jejunal and ileal samples were all dominated by Lactobacillus (duodenum: 0.57 ± 0.35, jejunum: 0.5 ± 0.33, ileum: 0.46 ± 0.26) and Enterococcus (duodenum: 0.36 ± 0.35, jejunum: 0.28 ± 0.33, ileum: 0.41 ± 0.32). Caecal samples remained significantly more diverse than duodenal (P = 0.013) and ileal (P = 0.013) samples, but not jejunal samples. There remained a significant difference in richness but not diversity between duodenal samples and jejunal (richness: P = 0.013) samples but not ileal samples. A significant difference in the richness but not diversity of samples was observed between the jejunum and ileum (richness: P = 0.026).
Between day 7 and day 14, significant changes in community composition were observed for caecal (P = 0.03), ileal (P = 0.05) and duodenal (P = 0.03) samples; a difference in diversity was observed in ileal samples (P = 0.022); and differences in richness were observed in duodenal (P = 0.022) and jejunal (P = 0.022) samples. Several OTUs were also found to have changed in abundance within sample types (Additional files 5, 6, 7 and 8). These changes appear to have led to increased homogeneity of microbial communities across all sample types (Fig. 1). No significant differences in community composition, richness or diversity were observed between the sample types at this timepoint and very few specific OTUs were found to be different between sample types (Additional file 11). This seems to have been driven by caecal samples attaining a “small-intestine like” microbiota dominated by Lactobacillus (0.39 ± 0.24) with smaller numbers of Lachnospiraceae_unclassified (0.10 ± 0.13), Candidatus_Arthromitus (0.1 ± 0.087), Escherichia-Shigella (0.098 ± 0.093) and Romboutsia (0.073 ± 0.018).
Between day 14 and 5 weeks, a significant difference can be observed in the microbiota composition of the ileum (P = 0.04, Additional file 8) and caeca (P = 0.03, Additional file 5) but not for jejunal and duodenal samples, despite significant changes in the abundance of specific OTUs (Additional files 6 and 7). Significant increases in richness and diversity were observed for caecal samples (P = 0.0216). The most common genera in caecal samples were Lachnospiraceae_unclassified (0.23 ± 0.050), Ruminococcaceae_unclassified (0.11 ± 0.016), Faecalibacterium (0.11 ± 0.022), Bifidobacterium (0.080 ± 0.041), Lactobacillus (0.066 ± 0.034), Clostridiales_vadinBB60_group_ge (0.055 ± 0.027). Ileal, duodenal and jejunal samples were predominated by Lactobacillus (duodenum: 0.81 ± 0.19, jejunum: 0.77 ± 0.11, ileum: 0.65 ± 0.20) with the next most common genus being Romboutsia (0.090 ± 0.15) in ileal samples and Staphylococcus in duodenal (0.060 ± 0.090) and jejunal (0.097 ± 0.090) samples. Large intestine samples included genera found in both small intestinal and caecal samples: Lactobacillus (0.21 ± 0.074), Lachnospiraceae_unclassified (0.15 ± 0.066), Escherichia-Shigella (0.12 ± 0.14), Romboutsia (0.092 ± 0.071), Faecalibacterium (0.067 ± 0.043), Ruminococcaceae_unclassified (0.061 ± 0.036) and Bifidobacterium (0.060 ± 0.036).
A significant difference was found by PERMANOVA between caecal samples and duodenal (P = 0.04) and ileal (P = 0.02) samples but not large intestine or jejunal samples (Fig. 1). The bacterial community composition of large intestine samples was found to be significantly different to all other sample types except the caeca (duodenum: P = 0.03, jejunum: P = 0.04, ileum: P = 0.04). This is reflected in our Deseq analysis where no OTUs were found to be significantly differently abundant between large intestine and caecal samples but many OTUs were differently abundant between the large intestine and other sample types (Additional file 12). No significant differences by PERMANOVA were observed between duodenal, jejunal and ileal samples, although several OTUs were found to be significantly differently abundant (Additional file 12).
In this study we compared samples taken from the duodenum, jejunum, ileum and caeca of broiler chickens at five timepoints from day of hatch to 5 weeks of age. Gaining an insight into how the microbiota changes over time at specific sites in the small intestine may lead to a better understanding of the microbial ecology of the chicken gut. We found changes in community composition, diversity and richness for all sample types over time. More specifically, we observed an increase in the richness of microbial communities in all gut sections, and a general increase in diversity except at the day 14 timepoint. While the succession of communities was different in each section of the gut, we can broadly say that the intestinal microbiota of our chickens was initially formed by a low diversity community of predominantly Clostridium_sensu_stricto_1 which diversified over time to contain a far greater variety of Clostridiales, with smaller numbers of other taxonomies. Samples taken from different locations within the small intestine were found to be similar to one another, whereas in the majority of timepoints caecal samples clustered significantly separately by composition from small intestinal samples. However, at several timepoints specific OTUs were found to be significantly more abundant in certain small intestinal locations in comparison to other sites, suggesting that there are differences between the microbial communities across the small intestine and therefore a sample taken from one site should not be taken as representative of the small intestine as a whole.
At the day 1 timepoint samples from the caeca and the small intestine were dominated by Clostridium_sensu_stricto_1. The lack of diversity observed at this timepoint was expected as previous studies have also shown that microbial communities in the chicken intestine in the first days of life have little diversity and are usually dominated by either members of the Enterobacteriaceae or the Clostridiaceae [9,10,11, 14, 15, 26]. It is possible that members of these communities act as founding species for chicken gut microbial communities, due to their ubiquitousness in the environment. It has been demonstrated that these early life communities can have a significant impact on later microbiota communities and bird phenotypes [10, 23] and it is therefore important to understand the origin of these founder microbes and the impact they may have.
Between day 1 and day 3 the proportion of Clostridium_sensu_stricto_1 greatly decreased in all of our sample types and while it remained present at later timepoints it was never highly abundant. At day 3 all small intestinal sample types were dominated by two genera: Lactobacillus and Enterococcus. These remained the dominant genera in small intestinal samples at all future timepoints. Lactobacillus has previously been noted as being highly abundant in the small intestine [12, 14, 17, 18, 20, 23] as has Enterococcus . In our study at this timepoint the caeca were dominated by Enterococcus, Escherichia-Shigella, Lactobacillus and unclassified members of the Enterobacteriaceae and the Lachnospiraceae.
No significant change was observed in the community composition of any sample type between day 3 and 7; however, at timepoint 14 days changes were observed in the caecal microbial communities which led them to no longer cluster separately from small intestinal samples by their community compositions. This homogenisation correlated with an increase in the relative abundance of members of Lactobacillus and Candidatus Arthromitus which have previously been associated with the small intestine [12, 23]. We are unsure of the cause of this seeming “homogenisation” of our samples at this timepoint and this phenomenon does not previously appear in the chicken literature. No changes were made to the diet or environment of the birds around this timepoint, and sampling, DNA extraction and PCR were conducted using the same methodologies as for other timepoints, and were carried out by the same individual. It is possible that these results represent a true change in the microbiota; however, we are unable to exclude the possibility that these results occurred due to some unknown error.
By the 5 week timepoint further changes in the microbial community composition led to the caecal microbiota again clustering separately from all other sample types, except to samples from the large intestine. It is likely that the composition of the large intestinal microbiota is highly dependent on when the caeca were last emptied, and the similarity between caecal and large intestinal samples in our study may have been caused by the chickens having recently ejected their caecal contents into the large intestine. Bias can easily be introduced into 16S studies at various methodological stages . However, good experimental design, including the use of appropriate controls, can greatly reduce its impact. We ran several reagent only controls alongside our samples. The most abundant genera in these controls were also found in our samples: Clostridium_sensu_stricto_1, Lactobacillus, Enterococcus and unclassified members of the Lachnospiraceae. It seems unlikely that the presence of these bacteria is due to reagent contamination as they are common intestinal colonisers, unlike the skin and environmental bacteria which tend to occur as contaminants in reagents . Two possible alternative explanations for the presence of these bacteria are index switching or cross-well contamination [29, 30]. Due to the low numbers of reads seen in the negative controls, contamination is unlikely to have driven any of the observed differences, though as always with microbiota studies, our results should be interpreted with care. It is also important to note that care must be taken when generalising our findings, due to the variability of the microbiota between chicken flocks [21, 22]. Our study was also conducted within a research facility with very high biosecurity standards which may not be representative of the conditions on a commercial farm.
While there is high variability present between the microbiota of chickens in different trials  it does seem as though there are taxonomies which are consistently present in the chicken gut . Identifying these shared taxa and elucidating their function within birds could lead to a greater understanding of the microbial ecology of the chicken gut and the development of technologies which could be widely applied to improve chicken health or productivity.
Ross 308 (Aviagen) chickens were hatched and housed at the National Avian Research Facility in Edinburgh (UK). Chickens were housed in groups in floor pens with wood shaving bedding and received food (wheat-based) and water ad libitum. Chickens received the Marek’s-Rispens vaccine and feed included coccidiostats. Samples were collected from chickens at the following ages: day 1 (n = 12), day 3 (n = 6), day 7 (n = 6), day 14 (n = 6) and day 35 (5 weeks: n = 6)). Collected samples included intestinal contents from the duodenum, jejunum, ileum and caeca. Large intestinal contents were also collected at 5 weeks. Animals were killed by cervical dislocation and their gut sections were immediately collected. The contents from each section were transferred into separate tubes within 5 min of collection, and then kept on ice.
Our experiments followed the general principles set out in Pollock et al. . Samples were stored at 4 °C for a maximum of 24 h until DNA extraction, except for those from DNA extraction batch 11 which were frozen at − 20 °C for 9 days prior to DNA extraction. The extraction batches to which each sample belonged can be found in Additional file 13. DNA extraction was performed as described previously using the DNeasy PowerLyzer PowerSoil Kit (Qiagen) . Reagent only controls were included for each batch of DNA extractions. DNA was also extracted from a mock community control (20 Strain Even Mix Whole Cell Material (ATCC® MSA-2002™)). The V4 region of the 16S rRNA gene was amplified from the extracted DNA, following a previously described method . The sequencing reaction consisted of 12.5 μl of Q5 High Fidelity DNA Polymerase (New England Biolabs), 1.25 μl each of custom 10 nM forward (5′–TATGGTAATTGTGTGCCAGCMGCCGCGGTAA–3′) and reverse primers (5′–AGTCAGTCAGCCGGACTACHVGGGTWTCTAAT–3′) with Illumina adaptor sequences, spacers and barcodes, 9 μl of nuclease free water (Qiagen) and 1 μl of DNA template. Cycling conditions were: 95 °C for 2 min followed by 30 cycles of 95 °C for 20 s, 55 °C for 15 s, 72 °C for 5 min followed by 72 °C for 10 min. Amplicons were purified using the AMPure XP System (Beckman Coulter), according to the manufacturer’s instructions, except that a 1:1 ratio of AMPure beads to sample was used. The amplicon concentrations were determined using the Qubit 3.0 Fluorometer (Life Technologies) with the Qubit dsDNA HS Assay Kit (Life Technologies), according to the manufacturer’s instructions. Samples were then pooled into an equimolar library, which was sequenced on an Illumina MiSeq v.2 (Illumina Inc.) producing 250 bp reads (1st read primer: 5′–TATGGTAATTGTGTGCCAGCMGCCGCGGTAA–3′, 2nd read primer: 5′– AGTCAGTCAGCCGGACTACHVGGGTWTCTAAT–3′, index primer: 5′– ATTAGAWACCCBDGTAGTCCGGCTGACTGACT–3′.
Mothur  was used for quality control of sequences, alignment, taxonomic assignment and OTU clustering, following a modified version of the MiSeq pipeline supplied on the mothur website . Sequences were removed if they were > 275 bp in length, contained ambiguous bases, had homopolymers of > 9 bp in length, did not align to the V4 region of the 16S rRNA gene or they did not originate from bacteria. Chimeras were identified using UCHIME  and were removed. The SILVA database , trimmed to the V4 16S rRNA gene, was used for alignment and taxonomic identification of sequences. Where the taxonomy of an OTU is labelled as ###_unclassified, this indicates that the OTU was unable to be identified to that level of taxonomy and the ### indicates the lowest taxonomy this OTU could be assigned to. OTUs were clustered by similarity using the dist.seqs and cluster commands from within mothur, using the default parameters.
Richness and diversity were calculated within mothur. Community richness was measured using the Chao 1 index. Community diversity was measured using the Inverse Simpson’s Index. The following analyses were performed in R (version 3.5.1.)  with the seed 8765. To compare the diversity and richness of microbial communities between groups, the Pairwise Wilcoxon Rank Sum Test was used with Bonferroni corrections to correct for multiple comparisons. Non-metric Multidimensional Scaling (NMDS) graphs were constructed using the Vegan  package and ggplot2 , using the Bray–Curtis dissimilarity. UpSet graphs were constructed using the UpSetR package . Beanplots were constructed using the beanplot package . PERMANOVA analyses were performed using the adonis function from the Vegan package; pairwise comparisons were performed using an adjusted version of the adonis function which output pairwise comparisons with Bonferroni corrections to correct for multiple comparisons. The package DESeq2  was used to calculate differences in abundances between groups for individual OTUs. Prior to statistical analysis, samples were subsampled to 10,000 reads, except for DEseq2 analysis.
Availability of data and materials
The paired-read fastq data that support the findings of this study have been deposited in in the European Nucleotide Archive with the accession code PRJEB33615.
- 16S rRNA:
16 Svedberg ribosomal ribonucleic acid
Non-metric Multidimensional Scaling
Operational taxonomic units
Permutational multivariate analysis of variance
Glendinning L, Stewart RD, Pallen MJ, Watson KA, Watson M. Assembly of hundreds of novel bacterial genomes from the chicken caecum. Preprint bioRxiv. 2019:699843.
Huws SA, Creevey CJ, Oyama LB, Mizrahi I, Denman SE, Popova M, Munoz-Tamayo R, Forano E, Waters SM, Hess M, et al. Addressing global ruminant agricultural challenges through understanding the rumen microbiome: past, present, and future. Front Microbiol. 2018;9:2161.
Aluthge ND, Van Sambeek DM, Carney-Hinkle EE, Li YS, Fernando SC, Burkey TE. The pig microbiota and the potential for harnessing the power of the microbiome to improve growth and health. J Animal Sci. 2019:skz208.
Shang Y, Kumar S, Oakley B, Kim WK. Chicken gut microbiota: importance and detection technology. Front Vet Sci. 2018;5:254.
Ikeda-Ohtsubo W, Brugman S, Warden CH, Rebel JMJ, Folkerts G, Pieterse CMJ. How can we define "optimal microbiota?": a comparative review of structure and functions of microbiota of animals, fish, and plants in agriculture. Front Nutr. 2018;5:90.
Jozefiak D, Rutkowski A, Martin SA. Carbohydrate fermentation in the avian ceca: a review. Anim Feed Sci Technol. 2004;113:1–15.
Karasawa Y. Significant role of the nitrogen recycling system through the ceca occurs in protein-depleted chickens. J Exp Zool. 1999;283:418–25.
Parsons CM, Potter LM, Brown RD. Effects of dietary carbohydrate and of intestinal microflora on excretion of endogenous amino-acids by poultry. Poult Sci. 1983;62:483–9.
Ocejo M, Oporto B, Hurtado A. 16S rRNA amplicon sequencing characterization of caecal microbiome composition of broilers and free-range slow-growing chickens throughout their productive lifespan. Sci Rep. 2019;9:2506.
Ballou AL, Ali RA, Mendoza MA, Ellis JC, Hassan HM, Croom WJ, Koci MD. Development of the chick microbiome: how early exposure influences future microbial diversity. Front Vet Sci. 2016;3:2.
Oakley BB, Buhr RJ, Ritz CW, Kiepper BH, Berrang ME, Seal BS, Cox NA. Successional changes in the chicken cecal microbiome during 42 days of growth are independent of organic acid feed additives. BMC Vet Res. 2014;10:282.
Kumar S, Chen CX, Indugu N, Werlang GO, Singh M, Kim WK, Thippareddi H. Effect of antibiotic withdrawal in feed on chicken gut microbial dynamics, immunity, growth performance and prevalence of foodborne pathogens. PLoS One. 2018;13:e0192450.
Danzeisen JL, Kim HB, Isaacson RE, Tu ZJ, Johnson TJ. Modulations of the chicken cecal microbiome and metagenome in response to anticoccidial and growth promoter treatment. PLoS One. 2011;6:e27949.
Lu JR, Idris U, Harmon B, Hofacre C, Maurer JJ, Lee MD. Diversity and succession of the intestinal bacterial community of the maturing broiler chicken. Appl Environ Microbiol. 2003;69:6816–24.
Richards P, Fothergill J, Bernardeau M, Wigley P. Development of the caecal microbiota in three broiler breeds. Front Vet Sci. 2019;6:201.
Oakley BB, Kogut MH. Spatial and temporal changes in the broiler chicken cecal and fecal microbiomes and correlations of bacterial taxa with cytokine gene expression. Front Vet Sci. 2016;3:11.
Choi JH, Lee K, Kim DW, Kil DY, Kim GB, Cha CJ. Influence of dietary avilamycin on ileal and cecal microbiota in broiler chickens. Poult Sci. 2018;97:970–9.
Fu SJ, Guo SJ, Wang JJ, Wang YM, Zhang ZM, Shen ZQ. Microbial community diversity of Jinghong laying hens at peak production based on 16S rRNA sequencing. J Appl Anim Res. 2018;46:1430–6.
Xi Y, Shuling N, Kunyuan T, Qiuyang Z, Hewen D, ChenCheng G, Tianhe Y, Liancheng L, Xin F. Characteristics of the intestinal flora of specific pathogen free chickens with age. Microb Pathog. 2019;132:325–34.
Ngunjiri JM, Taylor KJM, Abundo MC, Jang H, Elaish M, Mahesh KC, Ghorbani A, Wijeratne S, Weber BP, Johnson TJ, Lee CW. Farm stage, bird age, and body site dominantly affect the quantity, taxonomic composition, and dynamics of respiratory and gut microbiota of commercial layer chickens. Appl Environ Microbiol. 2019;85:e03137–18.
Stanley D, Geier MS, Hughes RJ, Denman SE, Moore RJ. Highly variable microbiota development in the chicken gastrointestinal tract. PLoS One. 2013;8:e84290.
Pandit RJ, Hinsu AT, Patel NV, Koringa PG, Jakhesara SJ, Thakkar JR, Shah TM, Limon G, Psifidi A, Guitian J, et al. Microbial diversity and community composition of caecal microbiota in commercial and indigenous Indian chickens determined using 16s rDNA amplicon sequencing. Microbiome. 2018;6:115.
Johnson TJ, Youmans BP, Noll S, Cardona C, Evans NP, Karnezos TP, Ngunjiri JM, Abundo MC, Lee CW. A consistent and predictable commercial broiler chicken bacterial microbiota in antibiotic-free production displays strong correlations with performance. Appl Environ Microbiol. 2018;84:e00362–18.
Bjerrum L, Engberg RM, Leser TD, Jensen BB, Finster K, Pedersen K. Microbial community composition of the ileum and cecum of broiler chickens as revealed by molecular and culture-based techniques. Poult Sci. 2006;85:1151–64.
Konsak BM, Stanley D, Haring VR, Geier MS, Hughes RJ, Howarth GS, Crowley TM, Moore RJ. Identification of differential duodenal gene expression levels and microbiota abundance correlated with differences in energy utilisation in chickens. Anim Prod Sci. 2013;53:1269–75.
Schokker D, Veninga G, Vastenhouw SA, Bossers A, de Bree FM, Kaal-Lansbergen L, Rebel JMJ, Smits MA. Early life microbial colonization of the gut and intestinal development differ between genetically divergent broiler lines. BMC Genomics. 2015;16:418.
Pollock J, Glendinning L, Wisedchanwet T, Watson M. The madness of microbiome: attempting to find consensus "best practice" for 16S microbiome studies. Appl Environ Microbiol. 2018;84:e02627–17.
Salter SJ, Cox MJ, Turek EM, Calus ST, Cookson WO, Moffatt MF, Turner P, Parkhill J, Loman NJ, Walker AW. Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biol. 2014;12:87.
Minich JJ, Sanders JG, Amir A, Humphrey G, Gilbert JA, Knight R. Quantifying and understanding well-to-well contamination in microbiome research. mSystems. 2019;4:e00186–19.
Larsson AJM, Stanley G, Sinha R, Weissman IL, Sandberg R. Computational correction of index switching in multiplexed sequencing libraries. Nat Methods. 2018;15:305–7.
Glendinning L, Wright S, Pollock J, Tennant P, Collie D, McLachlan G. Variability of the sheep lung microbiota. Appl Environ Microbiol. 2016;82:3225–38.
Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq illumina sequencing platform. Appl Environ Microbiol. 2013;79:5112–20.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glockner FO. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.
Team R Core: R: a language and environment for statistical computing. 2013.
Oksanen J, Blanchet F, Friendly M, Kindt R, Legendre P, McGlinn D, Minchin P, O'Hara R, Simpson G, Solymos P, et al: vegan: Community ecology package. R package version 25–2 edition.
Wickham H: ggplot2: Elegant graphics for data analysis. 2016.
Lex A, Gehlenborg N, Strobelt H, Vuillemot R, Pfister H. UpSet: visualization of intersecting sets. IEEE Trans Vis Comput Graph. 2014;20:1983–92.
Kampstra P. Beanplot: a boxplot alternative for visual comparison of distributions. J Stat Softw. 2008;28:9.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:38.
We would like to thank the staff at the Greenwood Building, Roslin Institute for the care of our animals. Sequencing was carried out by Edinburgh Genomics, The University of Edinburgh.
The Roslin Institute forms part of the Royal (Dick) School of Veterinary Studies, University of Edinburgh. This project was supported by the Biotechnology and Biological Sciences Research Council, including institute strategic programme and national capability awards to The Roslin Institute (BBSRC: BB/P013759/1, BB/P013732/1, BB/J004235/1, BB/J004243/1). The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Animals were housed in premises licensed under a UK Home Office Establishment License within the terms of the UK Home Office Animals (Scientific Procedures) Act 1986 (Study number PUN056: approved by the Roslin Institute named veterinary surgeon). Housing and husbandry complied with the Code of Practice for Housing and Care of Animals Bred, Supplied or Used for Scientific Purposes and were overseen by the Roslin Institute Animal Welfare and Ethical Review Board. Animals were culled by schedule one methods authorized by the Animals (Scientific Procedures) Act 1986. Birds were euthanized by cervical dislocation.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Abundance of OTUs in chicken intestinal samples (not sub-sampled).
Taxonomy of bacterial OTUs.
Contains Table 1 (Mock community composition by genus), Table 2 (Mock community composition by order), Table 3 (Diversity (Inverse Simpson) and richness (Chao 1 index) of subsampled samples), Figure 1 (UpSet graph – timepoints) and Figure 1 (UpSet graph – timepoints).
OTUs which were identified as being significantly more abundant by DESeq2 between locations within the chicken intestinal tract at day 1 of age.
OTUs which were identified as being significantly more abundant by DESeq2 between chicken caecal samples at different timepoints.
OTUs which were identified as being significantly more abundant by DESeq2 between chicken duodenal samples at different timepoints.
OTUs which were identified as being significantly more abundant by DESeq2 between chicken jejunal samples at different timepoints.
OTUs which were identified as being significantly more abundant by DESeq2 between chicken ileal samples at different timepoints.
OTUs which were identified as being significantly more abundant by DESeq2 between locations within the chicken intestinal tract at day 3 of age.
OTUs which were identified as being significantly more abundant by DESeq2 between locations within the chicken intestinal tract at day 7 of age.
OTUs which were identified as being significantly more abundant by DESeq2 between locations within the chicken intestinal tract at day 14 of age.
OTUs which were identified as being significantly more abundant by DESeq2 between locations within the chicken intestinal tract at 5 weeks of age.
About this article
Cite this article
Glendinning, L., Watson, K.A. & Watson, M. Development of the duodenal, ileal, jejunal and caecal microbiota in chickens. anim microbiome 1, 17 (2019). https://doi.org/10.1186/s42523-019-0017-z