BACKGROUND
the avian ceca generally have a more important role in digestion than the cecum in most mammals.
avian ceca are usually finger-shaped blind pouches, presenting as lateral extensions at the junction of the small and large intestine and are commonly present in pairs.
in birds the ceca vary considerably in size and morphology; ranging from very long, such as in most domestic poultry, to completely absent, such as in pigeons and parrots  <cit> .
in the chicken they reach 16–18 cm long in adult birds.
it is a multi-purpose organ vital to the bird’s physiology; a complex system inhabited by a very dense microbial community that converts the cecal pouches into fermentation powerhouses.
members of the cecal microbiota have the ability to digest cellulose, starch and other resistant polysaccharides  <cit> .
ceca are not only a major site of water absorption  <cit>  but are also a site of nutrient transport and absorption  <cit> .
recent advances in culture-free technologies for microbiota characterization have facilitated an increasing number of studies investigating changes in chicken microbiota within the gastrointestinal tract  following manipulation of feed, environment, or health.
it has been documented that chicken microbiota responds to changes in feed  <cit> , litter composition  <cit> , antibiotics  <cit>  and probiotic addition to feed  <cit> , disease  <cit>  and stress  <cit> .
recently the use of next generation sequencing techniques to study microbiota composition has been extended beyond phylogenetic analysis to also include functional analysis using metagenomics of whole cecal microbiotas  <cit> .
the important biological role of the cecum has resulted in it being a major focus of research into the influences of chicken microbiota on bird health and productivity  <cit> .
the ceca sample microbes from both descending and ascending microbiota via normal peristalsis and retrograde gut movements  <cit> .
the cecal content is emptied several times per day into the gut lumen and then regrows to fill the ceca.
this cyclic emptying of the ceca means that elements of the fecal microbiota must be directly derived from the ceca.
we were interested to examine the relationship between the population structure of cecal and fecal microbiotas to determine if fecal sampling via cloacal swabbing is an effective proxy for cecal sampling.
cloacal sampling, which is effectively a reliable and fast method of collecting fresh fecal samples from a specific bird, has the great advantage that it is easy to obtain a series of samples from the one bird over time, unlike the case with cecal samples where it is usual to sacrifice the bird to recover a sample and so only a single time-point snapshot can be obtained for any one bird.
this has meant that temporal studies of cecal microbiota have had to rely on population sampling with different birds studied at each time point.
the analysis of such data is complicated by the large inter-bird variation in microbiota structure that has recently been documented  <cit> .
if the relationship between fecal and cecal samples was understood and consistent then fecal samples could be used to determine at least some elements of the microbiota present in the ceca.
the fecal and cecal microbiota profiles of mice have been compared using a gel based analysis method and no obvious correlations were noted  <cit>  however the biology of the chicken ceca is very different to that of the mouse and so it is likely that this study gives little indication of the relationship of the two sites in chickens.
sekelja et al.
<cit>  used conventional sequencing of cloned 16s fragments to investigate the variability of chicken fecal microbiota and how this was related to emptying of different regions of the git.
here we report a large chicken study that used next generation sequencing technology to characterize git microbiota.
the study builds on previously reported work that analyzed the cecal microbiota of birds in a series of production efficiency trials  <cit> .
here we analyse further samples and extend the microbiota analysis to generate a statistically powerful data set to compare fecal and cecal sample pairs from  <dig> birds to determine the level and nature of similarity in microbiota structure between the two sample sites.
the relationship between the microbiota in the two cecal compartments within each bird was also investigated.
RESULTS
cecal microbiota is richer in otus and has fewer dominant otus compared to fecal microbiota
a total of  <dig> , <dig> sequences were obtained from the  <dig> samples included in this analysis .
after quality trimming and chimera checking  <dig> , <dig> sequences remained with a minimum length of  <dig> nucleotides, a maximum length of  <dig>  and an average read length of  <dig> nucleotides.
there were on average over  <dig>  sequences per sample.
a number of alpha diversity measures were inspected to compare within sample diversity of cecal and fecal samples.
all of the diversity indicators showed statistically significant differences, all with p < 1e− <dig>  the lowest possible p-value based on the 1e <dig> monte carlo permutations used .
both non-phylogenetic estimator chao <dig>  and phylogenetic diversity  indicator whole tree analysis showed higher diversity in cecal samples relative to fecal samples.
cecal richness, based on number of observed species, was significantly higher than richness in fecal communities.
the singles and doubles estimator suggested more rare otus in cecal samples.
simpson’s evenness  and strong’s dominance index confirmed that fecal samples have more dominant otus.
shannon entropy  and the equitability indices  showed cecal phylotypes as more evenly distributed within samples.figure  <dig>
alpha diversity is significantly different between fecal and cecal microbiota.
alpha diversity indicators of cecal  and fecal  samples all showed significant  difference based on 1e <dig> monte carlo permutations.
the indices plotted in the four panels are; a, chao 1; b, simpson; c, shannon; and d, equitability.
a range of beta diversity measures were inspected  and, similar to alpha diversity, they indicated significant differences between the fecal and cecal samples.
unweighted and weighted unifrac  both showed differences in microbial communities with adonis statistics p-values lower than p < 1e− <dig>  the lowest possible p-value based on the 1e <dig> permutations.
separation in both weighted and unweighted unifrac indicated that the microbiota profiles group on origin  based on presence/absence as well as abundance.figure  <dig>
unifrac distance shows microbiota differences between fecal and cecal origins.
a tree representing cecal  and fecal  samples.
the sample clustering is based on jacknifed weighted  and unweighted  unifrac.
tree files were generated in qiime and visualised in figtree.
the two communities were significantly  different based on both weighted and unweighted unifrac and adonis statistics using 1e <dig> permutations.
members of cecal and fecal microbiota
most of the otus were present in both cecal and fecal samples .
the shared otus, represented by  <dig> % of all otus, accounted for  <dig> % of all sequences.
the  <dig> % of otus that were exclusive to cecal samples comprised only  <dig> % of all sequences, indicating that they are mostly rare, low abundance otus.
fewer otus were exclusive to fecal samples , however they made up a greater proportion of all sequences  than the cecum exclusive otus.
otus present only in cecum were all, with the exception of one otu, unknown and uncultured relatives of bacteroides fragilis, faecalibacterium prausnitzii, and a number of unknown lactobacillus and clostridia.
on the other hand, 27% of otus exclusive to fecal origin shared sequence similarity higher than 97% to known type strains of bacteria not commonly found in chicken cecum such as streptococcus minor, vagococcus fluvialis, streptococcus henryi, staphylococcus gallinarum and staphylococcus aureus, all with 100% sequence identity to type strains.figure  <dig>
venn diagrams showing percentage of shared otus.
diagrams are show the percentage of shared otus  and the percentage of sequence reads that they represented  of cecal and fecal origin at 3% divergence.
<dig> % of all otus are shared between the sections and those otus represented  <dig> % of all sequences.
tables  <dig> and  <dig> show the  <dig> most abundant otus, based on sequence counts, from cecal and cloacal samples respectively, classified using eztaxon to their nearest culturable isolate.
the most abundant otu in chicken cecum, across the  <dig> trials, was  <dig> % identical to a type strain of bacteroides fragilis and represented  <dig> % of all cecal sequences based on the 100-times rarefied data .
this otu was the most abundant only in trial  <dig> with  <dig> %, third most abundant in trial  <dig>  and second most abundant in trial  <dig> with  <dig> % of total sequences.
the next  <dig> most abundant otus were confidently  classified as lactobacillus crispatus, lactobacillus johnsonii, lactobacillus salivarius and lactobacillus reuteri.
there were  <dig> more otus in the cecal top  <dig> classified with sequence similarity >97% to lactobacillus helveticus, lactobacillus vaginalis and parabacteroides distasonis.
the remaining otus, except for lactobacillus helveticus  were of unknown species with sequence similarity to the closest known isolates in the range of  <dig> - <dig> %.
cecal microbiota contained more unknown otus  than the fecal microbiota; those otus comprised 47% of cecal microbiota as opposed to  <dig> % in fecal samples, based on total cecal and total fecal sequences.table  <dig>
the  <dig> most abundant otus in chicken cecum
closest culturable isolate
isolate accession number
% similarity
% in cecal
% in fecal
otu id
bacteroides fragilis
lactobacillus crispatus
lactobacillus johnsonii
lactobacillus salivarius
lactobacillus reuteri
acholeplasma palmae
lactobacillus helveticus
butyricicoccus pullicaecorum
faecalibacterium prausnitzii
lactobacillus vaginalis
lactobacillus helveticus
pontibacillus litoralis
exiguobacterium acetylicum
ruminococcus albus
ruminococcus albus
parabacteroides distasonis
ruminococcus flavefaciens
clostridium cellobioparum
clostridium termitidis
clostridium symbiosum
the representative sequences for all abundant otus are publically available in embl database with otu id as identifier with accession numbers hg810851- hg <dig>
the  <dig> most abundant otus in chicken feces
closest culturable isolate
isolate accession number
% similarity
% in cecal
% in fecal
otu id
lactobacillus crispatus
lactobacillus salivarius
lactobacillus johnsonii
lactobacillus helveticus
lactobacillus reuteri jc
clostridium perfringens
lactobacillus crispatus
clostridium lituseburense
lactobacillus salivarius
lactobacillus helveticus
candidatus arthromitus sp.
lactobacillus salivarius
lactobacillus crispatus
lactobacillus johnsonii
bacteroides fragilis
lactobacillus pontis
lactobacillus helveticus
lactobacillus vaginalis
lactobacillus crispatus
lactobacillus gallinarum
the representative sequences for all abundant otus are publically available in embl database with otu id as identifier with accession numbers hg810851- hg <dig>
fecal microbiota was dominated by lactobacillus with the  <dig> most abundant otus making up  <dig> % of all fecal sequences.
the most abundant classified with >97% sequence similarity to type strains were l. crispatus, l. salivarius, l. johnsonii, l. helveticus, l. reuteri and l. vaginalis.
moreover, only  <dig> out of  <dig> most abundant otus were not lactobacilli: bacteroides fragilis, candidatus arthromitus, clostridium perfringens and unknown clostridium similar to c. lituseburense.
the remaining lactobacilli otus, which showed sequence similarity to the closest culturable isolates of between  <dig>  and  <dig> %, most likely comprise a number of novel members of this genus.
there were  <dig> otus out of  <dig> total fecal and cecal, classified as most closely related to lactobacillus strains,  <dig> of those with sequence similarity <95% to known culturable isolates.
fecal samples contained clostridium perfringens at an average level of  <dig> %, the 6th most abundant otu classified, with  <dig> % similarity to the type strain.
however, cecal samples contained on average only  <dig> % of this known chicken pathogen.
the carriage of c. perfringens varied markedly across the  <dig> trials; in trial  <dig> the fecal samples carried  <dig> %, trial  <dig> had  <dig> % and none was detected in trial  <dig>  although c. perfringens carriage was lower in the cecal samples the same relative trend between trials was noted with  <dig>   <dig>  and  <dig> % detected in trials 1– <dig> respectively.
low abundance fecal otus more closely correlate with cecal levels than high abundance otus
the main question we aimed to answer in this study was whether fecal samples could provide a reliable snapshot of cecal community structure.
the high number of sequences and samples across the three trials provides sufficient statistical power to attempt to make these predictions.
the three trials were inspected separately and as a whole set for correlations between fecal and cecal abundance of all shared otus present in both cecal and fecal samples across all of the birds.
all  <dig> trials showed positive correlation between cecal and fecal abundances .
to inspect if the level of correlation is influenced by abundance, i.e., if more abundant otus show better correlation, we inspected correlations for subsets of data with different cecal and separately fecal minimal abundance for all trials.
although we expected that more abundant otus would have higher cecal-fecal correlation, we found the opposite trend for both cecal and fecal abundances .
rare otus seem to be of similar low abundance in both cecal and fecal while higher abundance taxa tended to differ more.
this is especially notable in the fecal dominant otus.
otus comprising more than 50% of cecal sequences in one bird  had a wide range of fecal abundances.
the number of birds with cecal microbiota dominated by more than 50% was much lower than in fecal samples.
figure  <dig> simplifies the message from detailed figure 5a and b; shared otus present in lower abundance correlated better, while correlation for more abundant otus was very low.figure  <dig>
correlations between cecal and fecal abundance.
line graph representing correlations between cecal and fecal abundance in trial  <dig> , trial  <dig> , trial  <dig>  and complete set .
there is a reproducible positive correlation across the  <dig> trials.
all individual otus present in both cecal and fecal sample for all  <dig> birds were used in calculation, however, due to a very high number they are not displayed in the plot.
correlation plots for the different minimal cecal  and fecal  abundances.
note that each shared otu was plotted for all  <dig> birds and may be highly abundant in many of the birds, thus total percentages are not expected to add up to  <dig>  for example, there were  <dig> birds that had otus with cecal abundance higher than 50% vs  <dig> birds that contained fecal otus with abundance over 50%.
relationship between fecal/cecal correlation and levels of otu abundance.
correlations were plotted for both cecal  and fecal  samples across the range of otu abundances.
there was linear relationship between correlation and minimum abundance  showing higher correlation for low abundance otus and a very low correlation for high abundance otus.
the microbiota of cecal pairs are similar
the interpretation of cecal microbiota data could be influenced by sampling procedure if the pair of cecal pouches within a bird had radically different microbiota populations.
to address this issue, in an independent analysis, each of the pair of ceca within  <dig> birds were sampled and analysed according to the same methods detailed for the main experiment.
there was no statistically significant difference in alpha diversity between cecal pairs with shannon, simpson or observed species indices .
beta diversity was also unaffected by choice of ceca; unweighted unifrac differences using adonis statistics showed a p-value of  <dig>  and unweighted unifrac p-value of  <dig> .
in figure  <dig> it can be seen that the pairs of cecal samples from each bird tend to be closely related to each other and hence map close to each other in both the unweighted  and weighted  unifrac plots.figure  <dig>
differences between left and right ceca samples from the same bird represented using unifrac distance.
the unweighted  and weighted  unifrac pcoa plots demonstrate that left and right ceca are harbouring similar bacterial communities.
left and right ceca from each bird are represented with the same unique symbol and colour combination.
discussion
chicken cecal microbiota has been widely investigated due to the significant role of the cecum in heath, performance, and disease  <cit> .
the first insights into cecal microbiota that suggested the complexity of the microbial populations came from culture-based studies.
it has been suggested that the microbial diversity of chicken cecum has been altered with the common use of antibiotics in feed worldwide and that it has lost much of its natural chicken specialised microbiota and with it microbial potential and metabolic capabilities  <cit> .
ceca have been found to harbour microbiota capable of degrading cellulose and other indigestible carbohydrates and producing high amounts of beneficial metabolites such as short chain fatty acids   <cit> .
along the git the cecal community takes the most time to develop and mature  <cit>  and it contains the highest microbial diversity within the git  <cit> .
although it was reported that differences between trials and flocks in chickens can be extensive, overall cecal composition identified in the present study agrees with previously published data.
lu et al.
<cit>  found clostridia to dominate chicken cecum in their trial.
in the present study  <dig> of the  <dig> most abundant otus, based on sequence counts and not adjusting for rrna operon copy number, are most closely related to the clostridiales order.
however all of these had <95% sequence similarity to known culturable isolates, indicating that there is a high degree of cecal microbiota diversity that is yet to be revealed.
in general the clostridia have a reputation as bad, undesirable members of the microbiota, based on pathogenic potential of strains of c. perfringens, c. difficile, c. tetani and a few others.
however, this order also encompasses many beneficial bacteria such as cellulose and starch degraders; clostridium clusters iv and xiv have been linked to prevention of inflammatory bowel disease  and maintenance of mucosal homeostasis  <cit>  and clostridia protect from allergy and autoimmune disorders  <cit>  in mice.
atarashi et.
al.
<cit>  reported the colonisation of germ free mice with a mix of clostridium strains and also inspected mice enriched in clostridium abundance.
clostridium strains promoted treg cell accumulation.
oral administration to conventionally grown mice, to increase clostridia abundance, was beneficial for mice health, for example in inducing resistance to colitis by the means of suppressing weight loss, bleeding, colon shortening, edema and other symptoms of colitis allowing clostridia supplemented mice to show mild, if any, symptoms.
chicken cecum is a major source of uncultured clostridia that may represent enormous microbial potential.
the abundance of bacteroides, lactobacillus and clostridiales in the cecum in the present study is in broad agreement with previous studies  <cit> .
the benefits of lactobacillus are well known, cecal samples in this study have shown high diversity in this genus and indicated potential for further probiotic research by estimating  <dig> potentially novel lactobacillus-related species with similarity to the closest lactobacillus database matches of <95%.
the sequence count data indicate that the most abundant otu in the cecal microbiota was b. fragilis due to extreme values in one of the trials.
bacteroides are known to have beneficial effects on the host, being effective degraders of indigestible carbohydrates, especially cellulose and starch  <cit> .
b. fragilis produces scfa  <cit>  and when colonising germ free animals aids treg differentiation and il <dig> production  <cit> .
similarly another abundant phylotype related to faecalibacterium prausnitzii belongs to a butyrate producing cluster  <cit> .
absence of faecalibacterium prausnitzii is linked to crohn’s disease  <cit>  while cecal richness in ruminococcus is in line with the previously reported  <cit>  cellulose degrading potential  <cit>  of the ceca.
fecal microbiota was found to be dominated by lactobacillus.
this is consistent with previous studies  <cit> .
the high abundance of c. perfringens was attributed to one trial with extreme values and b. fragilis was present but at much lower abundance than in the cecum.
surprisingly we identified an abundant fecal otu, comprising  <dig> % of fecal sequences, as 100% identical to candidatus arthromitus sp.
lsfo <dig> ,lsfo <dig>  .
gong et al.
<cit>  reported high abundance of candidatus division arthromitus in chicken jejunum and ileum mucosal samples, representing 34% and 28% of sequences in these regions respectively.
they are found in gut microbiota of humans, chicken, rodents and fish, where they anchor to the intestinal epithelial cells in the ileum to act as immune system modulators  <cit> .
they are major activators of t-cells  <cit> , epithelial lymphocytes,  <cit>  and iga  <cit>  and could thus also be important in chicken health and performance.
until recently candidatus arthromitus was the term used for two physiologically similar groups of segmented filamentous bacteria , one commonly found in the gut of arthropods and others commonly found in git of vertebrates  <cit> .
although there were striking similarities between the two groups, thompson et al.
<cit>  demonstrated them to be distinct and unrelated.
arthropod inhabiting filamentous bacteria have yet unknown function while the ones inhabiting git of vertebrates, including chicken, play absolutely critical roles in immune function of the host  <cit> .
thompson et al.
suggested that git originating sfb form a monophyletic group in the clostridiaceae based on 16s sequence analysis, are now renamed “candidatus savagella”  <cit> .
however, arthropods specific filamentous bacteria were identified as members of a lachnospiraceae arthromitus cluster and should keep the name “candidatus arthromitus”  <cit> .
the abundance of candidatus savagella in chicken fecal microbiota is not reproducible across the studies.
in our own opinion its identification is influenced by the taxonomic databases and algorithms used and in many cases they fall under unknown and uncultured bacteria.
it is possible that their role in chicken health, especially immunity is significant therefore more attention should be given to deeper taxonomic identification of significant otus in chicken studies.
sekelja et al.
<cit>  inspected the influence of other git sections on chicken fecal microbiota over 16 days.
they detected massive temporal variations in fecal microbiota and different profiles to match different sections of git.
they proposed that the major reason for the temporal variation was periodic emptying of different git sections.
the emptying of different git areas may influence fecal profile but fecal analysis still remains a useful and powerful approach for microbiota studies in animals with proven success record in human and mammalian studies.
it is possible that the timing of sampling and cecal emptying is the reason for the existence of a number of outlier birds with higher similarity between fecal and cecal communities.
in an auxiliary study we found that there were strong similarities between the microbiota compositions of pairs of ceca from birds, indicating that the choice of ceca is unlikely to be a concern in sampling procedures.
perhaps a major lesson to be learnt is that it may be misleading to draw conclusions from just a few samples; a large number of samples will represent a range of different emptying events from different regions of the git and may allow a more complete and representative overall picture to emerge of microbiota composition.
in the present study we compared fecal and cecal microbiota across many birds and three independent trials in order to investigate the relationship between the two most commonly used and reported sample types for microbiota analysis in birds.
we were particularly interested in determining the extent to which the fecal microbiota may reflect the content of the cecal microbiota.
from this extensive sample set we can conclude that for the majority of otus , comprising  <dig> % of sequences, presence in fecal implies presence in cecal community.
our results show the highest correlation between fecal and cecal samples within the rare biosphere species, the correlation analysis however, included only shared otus, but all otus that were not shared  were also from the low abundance group.
therefore, high correlation among rare otus cannot be used to exactly predict cecal from fecal abundance since some rare otus are also likely not to be detected in the other community at all.
the data shows that otus exclusive to cecum, which would be missed if fecal origin was sampled instead, were all of low abundance.
based on the positive correlation for nearly all shared otus it is possible that the low abundance otus that appear to be unique to the cecal or to fecal samples would be identified in the opposite group if the depth of sequencing was higher as they are near the detection limit and the limits of the exclusion criteria used in quality filtering parameters to generate the data set.
although the alpha diversity plots indicated satisfactory sequence coverage an increase in sequencing depth may identify more shared otus but is also likely to detect additional rare otus.
there is a positive overall correlation between cecal and fecal abundance as shown in figure 4; this correlation is however negligible in the more dominant otus comprising over 10% of total sequences.
this trend was found to be reproducible over the three trials that harbour very different microbial communities.
CONCLUSIONS
based on the findings of the present study, we anticipate that with enough sequencing depth, fecal samples can be used to reflect the presence and absence of the vast majority of the members of the cecal community.
however, the two communities are generally very distinct based on alpha and beta diversity.
positive correlation cannot be used to accurately predict otu counts in the low abundance part of the microbiota due to it also carrying the highest number of unshared otus, nor in high abundance otus where correlation is minimal.
however, regardless of community structure differences, i.e., species showing different abundance, dominance and equitability, both cecal and fecal microbiota analyses are likely to accurately report if a treatment or condition has induced changes in microbiota.
this conclusion is based on the high number of shared species that represent  <dig> % of all community members, which would be part of the community response to treatment.
although control/treatment differences could be detected using either cecal or fecal samples, conclusions drawn from the separate analyses are likely to differ.
hence choice of sampling site remains critical in experimental design as fecal microbiota do not provide a complete indication of cecal community structure.

 the avian ceca generally have a more important role in digestion than the cecum in most mammals. avian ceca are usually finger-shaped blind pouches, presenting as lateral extensions at the junction of the small and large intestine and are commonly present in pairs. in birds the ceca vary considerably in size and morphology; ranging from very long, such as in most domestic poultry, to completely absent, such as in pigeons and parrots. in the chicken they reach 16–18 cm long in adult birds. it is a multi-purpose organ vital to the bird’s physiology; a complex system inhabited by a very dense microbial community that converts the cecal pouches into fermentation powerhouses. members of the cecal microbiota have the ability to digest cellulose, starch and other resistant polysaccharides. ceca are not only a major site of water absorption but are also a site of nutrient transport and absorption. recent advances in culture-free technologies for microbiota characterization have facilitated an increasing number of studies investigating changes in chicken microbiota within the gastrointestinal tract following manipulation of feed, environment, or health. it has been documented that chicken microbiota responds to changes in feed, litter composition, antibiotics and probiotic addition to feed, disease and stress. recently the use of next generation sequencing techniques to study microbiota composition has been extended beyond phylogenetic analysis to also include functional analysis using metagenomics of whole cecal microbiotas. the important biological role of the cecum has resulted in it being a major focus of research into the influences of chicken microbiota on bird health and productivity. the ceca sample microbes from both descending and ascending microbiota via normal peristalsis and retrograde gut movements. the cecal content is emptied several times per day into the gut lumen and then regrows to fill the ceca. this cyclic emptying of the ceca means that elements of the fecal microbiota must be directly derived from the ceca. we were interested to examine the relationship between the population structure of cecal and fecal microbiotas to determine if fecal sampling via cloacal swabbing is an effective proxy for cecal sampling. cloacal sampling, which is effectively a reliable and fast method of collecting fresh fecal samples from a specific bird, has the great advantage that it is easy to obtain a series of samples from the one bird over time, unlike the case with cecal samples where it is usual to sacrifice the bird to recover a sample and so only a single time-point snapshot can be obtained for any one bird. this has meant that temporal studies of cecal microbiota have had to rely on population sampling with different birds studied at each time point. the analysis of such data is complicated by the large inter-bird variation in microbiota structure that has recently been documented. if the relationship between fecal and cecal samples was understood and consistent then fecal samples could be used to determine at least some elements of the microbiota present in the ceca. the fecal and cecal microbiota profiles of mice have been compared using a gel based analysis method and no obvious correlations were noted however the biology of the chicken ceca is very different to that of the mouse and so it is likely that this study gives little indication of the relationship of the two sites in chickens. sekelja et al. used conventional sequencing of cloned 16s fragments to investigate the variability of chicken fecal microbiota and how this was related to emptying of different regions of the git. here we report a large chicken study that used next generation sequencing technology to characterize git microbiota. the study builds on previously reported work that analyzed the cecal microbiota of birds in a series of production efficiency trials. here we analyse further samples and extend the microbiota analysis to generate a statistically powerful data set to compare fecal and cecal sample pairs from birds to determine the level and nature of similarity in microbiota structure between the two sample sites. the relationship between the microbiota in the two cecal compartments within each bird was also investigated.  cecal microbiota is richer in otus and has fewer dominant otus compared to fecal microbiota a total of, sequences were obtained from the samples included in this analysis. after quality trimming and chimera checking, sequences remained with a minimum length of nucleotides, a maximum length of and an average read length of nucleotides. there were on average over sequences per sample. a number of alpha diversity measures were inspected to compare within sample diversity of cecal and fecal samples. all of the diversity indicators showed statistically significant differences, all with p < 1e− the lowest possible p-value based on the 1e monte carlo permutations used. both non-phylogenetic estimator chao and phylogenetic diversity indicator whole tree analysis showed higher diversity in cecal samples relative to fecal samples. cecal richness, based on number of observed species, was significantly higher than richness in fecal communities. the singles and doubles estimator suggested more rare otus in cecal samples. simpson’s evenness and strong’s dominance index confirmed that fecal samples have more dominant otus. shannon entropy and the equitability indices showed cecal phylotypes as more evenly distributed within samples.figure alpha diversity is significantly different between fecal and cecal microbiota. alpha diversity indicators of cecal and fecal samples all showed significant difference based on 1e monte carlo permutations. the indices plotted in the four panels are; a, chao 1; b, simpson; c, shannon; and d, equitability. a range of beta diversity measures were inspected and, similar to alpha diversity, they indicated significant differences between the fecal and cecal samples. unweighted and weighted unifrac both showed differences in microbial communities with adonis statistics p-values lower than p < 1e− the lowest possible p-value based on the 1e permutations. separation in both weighted and unweighted unifrac indicated that the microbiota profiles group on origin based on presence/absence as well as abundance.figure unifrac distance shows microbiota differences between fecal and cecal origins. a tree representing cecal and fecal samples. the sample clustering is based on jacknifed weighted and unweighted unifrac. tree files were generated in qiime and visualised in figtree. the two communities were significantly different based on both weighted and unweighted unifrac and adonis statistics using 1e permutations. members of cecal and fecal microbiota most of the otus were present in both cecal and fecal samples. the shared otus, represented by % of all otus, accounted for % of all sequences. the % of otus that were exclusive to cecal samples comprised only % of all sequences, indicating that they are mostly rare, low abundance otus. fewer otus were exclusive to fecal samples, however they made up a greater proportion of all sequences than the cecum exclusive otus. otus present only in cecum were all, with the exception of one otu, unknown and uncultured relatives of bacteroides fragilis, faecalibacterium prausnitzii, and a number of unknown lactobacillus and clostridia. on the other hand, 27% of otus exclusive to fecal origin shared sequence similarity higher than 97% to known type strains of bacteria not commonly found in chicken cecum such as streptococcus minor, vagococcus fluvialis, streptococcus henryi, staphylococcus gallinarum and staphylococcus aureus, all with 100% sequence identity to type strains.figure venn diagrams showing percentage of shared otus. diagrams are show the percentage of shared otus and the percentage of sequence reads that they represented of cecal and fecal origin at 3% divergence. % of all otus are shared between the sections and those otus represented % of all sequences. tables  and show the most abundant otus, based on sequence counts, from cecal and cloacal samples respectively, classified using eztaxon to their nearest culturable isolate. the most abundant otu in chicken cecum, across the trials, was % identical to a type strain of bacteroides fragilis and represented % of all cecal sequences based on the 100-times rarefied data. this otu was the most abundant only in trial with %, third most abundant in trial and second most abundant in trial with % of total sequences. the next most abundant otus were confidently classified as lactobacillus crispatus, lactobacillus johnsonii, lactobacillus salivarius and lactobacillus reuteri. there were more otus in the cecal top classified with sequence similarity >97% to lactobacillus helveticus, lactobacillus vaginalis and parabacteroides distasonis. the remaining otus, except for lactobacillus helveticus were of unknown species with sequence similarity to the closest known isolates in the range of - %. cecal microbiota contained more unknown otus than the fecal microbiota; those otus comprised 47% of cecal microbiota as opposed to % in fecal samples, based on total cecal and total fecal sequences.table the most abundant otus in chicken cecum closest culturable isolate isolate accession number % similarity % in cecal % in fecal otu id bacteroides fragilis lactobacillus crispatus lactobacillus johnsonii lactobacillus salivarius lactobacillus reuteri acholeplasma palmae lactobacillus helveticus butyricicoccus pullicaecorum faecalibacterium prausnitzii lactobacillus vaginalis lactobacillus helveticus pontibacillus litoralis exiguobacterium acetylicum ruminococcus albus ruminococcus albus parabacteroides distasonis ruminococcus flavefaciens clostridium cellobioparum clostridium termitidis clostridium symbiosum the representative sequences for all abundant otus are publically available in embl database with otu id as identifier with accession numbers hg810851- hg the most abundant otus in chicken feces closest culturable isolate isolate accession number % similarity % in cecal % in fecal otu id lactobacillus crispatus lactobacillus salivarius lactobacillus johnsonii lactobacillus helveticus lactobacillus reuteri jc clostridium perfringens lactobacillus crispatus clostridium lituseburense lactobacillus salivarius lactobacillus helveticus candidatus arthromitus sp. lactobacillus salivarius lactobacillus crispatus lactobacillus johnsonii bacteroides fragilis lactobacillus pontis lactobacillus helveticus lactobacillus vaginalis lactobacillus crispatus lactobacillus gallinarum the representative sequences for all abundant otus are publically available in embl database with otu id as identifier with accession numbers hg810851- hg fecal microbiota was dominated by lactobacillus with the most abundant otus making up % of all fecal sequences. the most abundant classified with >97% sequence similarity to type strains were l. crispatus, l. salivarius, l. johnsonii, l. helveticus, l. reuteri and l. vaginalis. moreover, only out of most abundant otus were not lactobacilli: bacteroides fragilis, candidatus arthromitus, clostridium perfringens and unknown clostridium similar to c. lituseburense. the remaining lactobacilli otus, which showed sequence similarity to the closest culturable isolates of between and %, most likely comprise a number of novel members of this genus. there were otus out of total fecal and cecal, classified as most closely related to lactobacillus strains, of those with sequence similarity <95% to known culturable isolates. fecal samples contained clostridium perfringens at an average level of %, the 6th most abundant otu classified, with % similarity to the type strain. however, cecal samples contained on average only % of this known chicken pathogen. the carriage of c. perfringens varied markedly across the trials; in trial the fecal samples carried %, trial had % and none was detected in trial although c. perfringens carriage was lower in the cecal samples the same relative trend between trials was noted with and % detected in trials 1– respectively. low abundance fecal otus more closely correlate with cecal levels than high abundance otus the main question we aimed to answer in this study was whether fecal samples could provide a reliable snapshot of cecal community structure. the high number of sequences and samples across the three trials provides sufficient statistical power to attempt to make these predictions. the three trials were inspected separately and as a whole set for correlations between fecal and cecal abundance of all shared otus present in both cecal and fecal samples across all of the birds. all trials showed positive correlation between cecal and fecal abundances. to inspect if the level of correlation is influenced by abundance, i.e., if more abundant otus show better correlation, we inspected correlations for subsets of data with different cecal and separately fecal minimal abundance for all trials. although we expected that more abundant otus would have higher cecal-fecal correlation, we found the opposite trend for both cecal and fecal abundances. rare otus seem to be of similar low abundance in both cecal and fecal while higher abundance taxa tended to differ more. this is especially notable in the fecal dominant otus. otus comprising more than 50% of cecal sequences in one bird had a wide range of fecal abundances. the number of birds with cecal microbiota dominated by more than 50% was much lower than in fecal samples. figure  simplifies the message from detailed figure 5a and b; shared otus present in lower abundance correlated better, while correlation for more abundant otus was very low.figure correlations between cecal and fecal abundance. line graph representing correlations between cecal and fecal abundance in trial, trial, trial and complete set. there is a reproducible positive correlation across the trials. all individual otus present in both cecal and fecal sample for all birds were used in calculation, however, due to a very high number they are not displayed in the plot. correlation plots for the different minimal cecal and fecal abundances. note that each shared otu was plotted for all birds and may be highly abundant in many of the birds, thus total percentages are not expected to add up to for example, there were birds that had otus with cecal abundance higher than 50% vs birds that contained fecal otus with abundance over 50%. relationship between fecal/cecal correlation and levels of otu abundance. correlations were plotted for both cecal and fecal samples across the range of otu abundances. there was linear relationship between correlation and minimum abundance showing higher correlation for low abundance otus and a very low correlation for high abundance otus. the microbiota of cecal pairs are similar the interpretation of cecal microbiota data could be influenced by sampling procedure if the pair of cecal pouches within a bird had radically different microbiota populations. to address this issue, in an independent analysis, each of the pair of ceca within birds were sampled and analysed according to the same methods detailed for the main experiment. there was no statistically significant difference in alpha diversity between cecal pairs with shannon, simpson or observed species indices. beta diversity was also unaffected by choice of ceca; unweighted unifrac differences using adonis statistics showed a p-value of and unweighted unifrac p-value of. in figure  it can be seen that the pairs of cecal samples from each bird tend to be closely related to each other and hence map close to each other in both the unweighted and weighted unifrac plots.figure differences between left and right ceca samples from the same bird represented using unifrac distance. the unweighted and weighted unifrac pcoa plots demonstrate that left and right ceca are harbouring similar bacterial communities. left and right ceca from each bird are represented with the same unique symbol and colour combination. discussion chicken cecal microbiota has been widely investigated due to the significant role of the cecum in heath, performance, and disease. the first insights into cecal microbiota that suggested the complexity of the microbial populations came from culture-based studies. it has been suggested that the microbial diversity of chicken cecum has been altered with the common use of antibiotics in feed worldwide and that it has lost much of its natural chicken specialised microbiota and with it microbial potential and metabolic capabilities. ceca have been found to harbour microbiota capable of degrading cellulose and other indigestible carbohydrates and producing high amounts of beneficial metabolites such as short chain fatty acids. along the git the cecal community takes the most time to develop and mature and it contains the highest microbial diversity within the git. although it was reported that differences between trials and flocks in chickens can be extensive, overall cecal composition identified in the present study agrees with previously published data. lu et al. found clostridia to dominate chicken cecum in their trial. in the present study of the most abundant otus, based on sequence counts and not adjusting for rrna operon copy number, are most closely related to the clostridiales order. however all of these had <95% sequence similarity to known culturable isolates, indicating that there is a high degree of cecal microbiota diversity that is yet to be revealed. in general the clostridia have a reputation as bad, undesirable members of the microbiota, based on pathogenic potential of strains of c. perfringens, c. difficile, c. tetani and a few others. however, this order also encompasses many beneficial bacteria such as cellulose and starch degraders; clostridium clusters iv and xiv have been linked to prevention of inflammatory bowel disease and maintenance of mucosal homeostasis and clostridia protect from allergy and autoimmune disorders in mice. atarashi et. al. reported the colonisation of germ free mice with a mix of clostridium strains and also inspected mice enriched in clostridium abundance. clostridium strains promoted treg cell accumulation. oral administration to conventionally grown mice, to increase clostridia abundance, was beneficial for mice health, for example in inducing resistance to colitis by the means of suppressing weight loss, bleeding, colon shortening, edema and other symptoms of colitis allowing clostridia supplemented mice to show mild, if any, symptoms. chicken cecum is a major source of uncultured clostridia that may represent enormous microbial potential. the abundance of bacteroides, lactobacillus and clostridiales in the cecum in the present study is in broad agreement with previous studies. the benefits of lactobacillus are well known, cecal samples in this study have shown high diversity in this genus and indicated potential for further probiotic research by estimating potentially novel lactobacillus-related species with similarity to the closest lactobacillus database matches of <95%. the sequence count data indicate that the most abundant otu in the cecal microbiota was b. fragilis due to extreme values in one of the trials. bacteroides are known to have beneficial effects on the host, being effective degraders of indigestible carbohydrates, especially cellulose and starch. b. fragilis produces scfa and when colonising germ free animals aids treg differentiation and il production. similarly another abundant phylotype related to faecalibacterium prausnitzii belongs to a butyrate producing cluster. absence of faecalibacterium prausnitzii is linked to crohn’s disease while cecal richness in ruminococcus is in line with the previously reported cellulose degrading potential of the ceca. fecal microbiota was found to be dominated by lactobacillus. this is consistent with previous studies. the high abundance of c. perfringens was attributed to one trial with extreme values and b. fragilis was present but at much lower abundance than in the cecum. surprisingly we identified an abundant fecal otu, comprising % of fecal sequences, as 100% identical to candidatus arthromitus sp. lsfo,lsfo. gong et al. reported high abundance of candidatus division arthromitus in chicken jejunum and ileum mucosal samples, representing 34% and 28% of sequences in these regions respectively. they are found in gut microbiota of humans, chicken, rodents and fish, where they anchor to the intestinal epithelial cells in the ileum to act as immune system modulators. they are major activators of t-cells, epithelial lymphocytes, and iga and could thus also be important in chicken health and performance. until recently candidatus arthromitus was the term used for two physiologically similar groups of segmented filamentous bacteria, one commonly found in the gut of arthropods and others commonly found in git of vertebrates. although there were striking similarities between the two groups, thompson et al. demonstrated them to be distinct and unrelated. arthropod inhabiting filamentous bacteria have yet unknown function while the ones inhabiting git of vertebrates, including chicken, play absolutely critical roles in immune function of the host. thompson et al. suggested that git originating sfb form a monophyletic group in the clostridiaceae based on 16s sequence analysis, are now renamed “candidatus savagella”. however, arthropods specific filamentous bacteria were identified as members of a lachnospiraceae arthromitus cluster and should keep the name “candidatus arthromitus”. the abundance of candidatus savagella in chicken fecal microbiota is not reproducible across the studies. in our own opinion its identification is influenced by the taxonomic databases and algorithms used and in many cases they fall under unknown and uncultured bacteria. it is possible that their role in chicken health, especially immunity is significant therefore more attention should be given to deeper taxonomic identification of significant otus in chicken studies. sekelja et al. inspected the influence of other git sections on chicken fecal microbiota over 16 days. they detected massive temporal variations in fecal microbiota and different profiles to match different sections of git. they proposed that the major reason for the temporal variation was periodic emptying of different git sections. the emptying of different git areas may influence fecal profile but fecal analysis still remains a useful and powerful approach for microbiota studies in animals with proven success record in human and mammalian studies. it is possible that the timing of sampling and cecal emptying is the reason for the existence of a number of outlier birds with higher similarity between fecal and cecal communities. in an auxiliary study we found that there were strong similarities between the microbiota compositions of pairs of ceca from birds, indicating that the choice of ceca is unlikely to be a concern in sampling procedures. perhaps a major lesson to be learnt is that it may be misleading to draw conclusions from just a few samples; a large number of samples will represent a range of different emptying events from different regions of the git and may allow a more complete and representative overall picture to emerge of microbiota composition. in the present study we compared fecal and cecal microbiota across many birds and three independent trials in order to investigate the relationship between the two most commonly used and reported sample types for microbiota analysis in birds. we were particularly interested in determining the extent to which the fecal microbiota may reflect the content of the cecal microbiota. from this extensive sample set we can conclude that for the majority of otus, comprising % of sequences, presence in fecal implies presence in cecal community. our results show the highest correlation between fecal and cecal samples within the rare biosphere species, the correlation analysis however, included only shared otus, but all otus that were not shared were also from the low abundance group. therefore, high correlation among rare otus cannot be used to exactly predict cecal from fecal abundance since some rare otus are also likely not to be detected in the other community at all. the data shows that otus exclusive to cecum, which would be missed if fecal origin was sampled instead, were all of low abundance. based on the positive correlation for nearly all shared otus it is possible that the low abundance otus that appear to be unique to the cecal or to fecal samples would be identified in the opposite group if the depth of sequencing was higher as they are near the detection limit and the limits of the exclusion criteria used in quality filtering parameters to generate the data set. although the alpha diversity plots indicated satisfactory sequence coverage an increase in sequencing depth may identify more shared otus but is also likely to detect additional rare otus. there is a positive overall correlation between cecal and fecal abundance as shown in figure 4; this correlation is however negligible in the more dominant otus comprising over 10% of total sequences. this trend was found to be reproducible over the three trials that harbour very different microbial communities.  based on the findings of the present study, we anticipate that with enough sequencing depth, fecal samples can be used to reflect the presence and absence of the vast majority of the members of the cecal community. however, the two communities are generally very distinct based on alpha and beta diversity. positive correlation cannot be used to accurately predict otu counts in the low abundance part of the microbiota due to it also carrying the highest number of unshared otus, nor in high abundance otus where correlation is minimal. however, regardless of community structure differences, i.e., species showing different abundance, dominance and equitability, both cecal and fecal microbiota analyses are likely to accurately report if a treatment or condition has induced changes in microbiota. this conclusion is based on the high number of shared species that represent % of all community members, which would be part of the community response to treatment. although control/treatment differences could be detected using either cecal or fecal samples, conclusions drawn from the separate analyses are likely to differ. hence choice of sampling site remains critical in experimental design as fecal microbiota do not provide a complete indication of cecal community structure. 
cecal samples and fecal samples via cloacal swab were collected from birds across replicate trials. the fecal and cecal samples were investigated to determine to what extent the microbiota found in fecal samples represented the microbiota of the ceca. fecal samples can be effectively used to detect some shifts and responses of cecal microbiota. fecal sample analysis captures a large percentage of the microbial diversity present in the ceca. the fecal microbiota is qualitatively similar to cecal microbiota but quantitatively different. the microbial populations present within the paired ceca of individual birds were also compared and shown to be similar. there was a positive correlation between cecal and fecal abundance in the shared sequences, however the two communities differed significantly in community structure, represented as either alpha or beta diversity. it was found that % of all operational taxonomic units, containing % of all sequences, were shared between the two sample types, with otus unique for each sample type found to be very rare. this is a more meaningful and preferred way of sampling as the same animals can be monitored and tracked for temporal studies.
 the majority of chicken microbiota studies have used the ceca as a sampling site due to the specific role of ceca in chicken productivity, health and wellbeing. however, sampling from ceca and other gastrointestinal tract sections requires the bird to be sacrificed. in contrast, fecal sampling does not require sacrifice and thus allows the same bird to be sampled repeatedly over time. this is a more meaningful and preferred way of sampling as the same animals can be monitored and tracked for temporal studies. the commonly used practice of selecting a subset of birds at each time-point for sacrifice and sampling introduces added variability due to the known animal to animal variation in microbiota.  cecal samples and fecal samples via cloacal swab were collected from birds across replicate trials. dna was extracted and 16s rrna gene sequences amplified and pyrosequenced to determine and compare the phylogenetic profile of the microbiota within each sample. the fecal and cecal samples were investigated to determine to what extent the microbiota found in fecal samples represented the microbiota of the ceca. it was found that % of all operational taxonomic units, containing % of all sequences, were shared between the two sample types, with otus unique for each sample type found to be very rare. there was a positive correlation between cecal and fecal abundance in the shared sequences, however the two communities differed significantly in community structure, represented as either alpha or beta diversity. the microbial populations present within the paired ceca of individual birds were also compared and shown to be similar.  fecal sample analysis captures a large percentage of the microbial diversity present in the ceca. however, the qualitative similarities in otu presence are not a good representation of the proportions of otus within the microbiota from each sampling site. the fecal microbiota is qualitatively similar to cecal microbiota but quantitatively different. fecal samples can be effectively used to detect some shifts and responses of cecal microbiota. keywords intestinalmicrobiotachickenfecalcecalissue-copyright-statement© the author 2015 