- Research article
- Open Access
Distinct microbiotas of anatomical gut regions display idiosyncratic seasonal variation in an avian folivore
Animal Microbiome volume 1, Article number: 2 (2019)
Current knowledge about seasonal variation in the gut microbiota of vertebrates is limited to a few studies based on mammalian fecal samples. Seasonal changes in the microbiotas of functionally distinct gut regions remain unexplored. We investigated seasonal variation (summer versus winter) and regionalization of the microbiotas of the crop, ventriculus, duodenum, cecum, and colon of the greater sage-grouse (Centrocercus urophasianus), an avian folivore specialized on the toxic foliage of sagebrush (Artemesia spp.) in western North America.
We sequenced the V4 region of the 16S rRNA gene on an Illumina MiSeq and obtained 6,639,051 sequences with a median of 50,232 per sample. These sequences were assigned to 457 bacterial and 4 archaeal OTUs. Firmicutes (53.0%), Bacteroidetes (15.2%), Actinobacteria (10.7%), and Proteobacteria (10.1%)were the most abundant and diverse phyla. Microbial composition and richness showed significant differences among gut regions and between summer and winter. Gut region explained almost an order of magnitude more variance in our dataset than did season or the gut region × season interaction. The effect of season was uneven among gut regions. Microbiotas of the crop and cecum showed the greatest seasonal differences.
Our data suggest that seasonal differences in gut microbiota reflect seasonal variation in the microbial communities associated with food and water. Strong differentiation and uneven seasonal changes in the composition and richness of the microbiota among functionally distinct gut regions demonstrate the necessity of wider anatomical sampling for studies of composition and dynamics of the gut microbiota.
The vertebrate gut is colonized by diverse microbial communities [1, 2]. The rapidly growing appreciation of the microbial contribution to digestion, immune function, development, and reproduction in vertebrates is primarily based on studies of humans and a relatively small number of domestic and captive species of mammals and birds [2,3,4]. Consequently, our knowledge of the factors that affect the assembly and dynamics of the gut micriobiota of wild species is limited [3, 5]. Vertebrate diets often exhibit marked seasonal variation in the composition, abundance, and nutritional quality of food resources. Concomitant seasonal changes in fecal microbiotas have been demonstrated in indigenous humans , giant pandas , squirrels , wood mice , macaques , and ground squirrels . However, the generality of this pattern has not been investigated in non-mammalian vertebrates.
Seasonal variation of gut microbiota in vertebrates with specialized diets rich in secondary plant compounds is of particular interest because of the hypothesized functional role of prokaryotes in the degradation and metabolism of dietary toxins . The greater sage-grouse (Centrocercus urophasianus; Fig. 1) exhibits one of the most specialized diets among the 10,135 extant species of birds [13, 14]. It feeds predominately on the foliage of chemically-defended sagebrush (Artemesia spp.) and its geographic range coincides with the distribution of sagebrush-dominated habitat (Fig. 1) in western North America . Sagebrush foliage is rich in toxic monoterpenes [16,17,18,19,20,21], phenolics [20, 22], and sesquiterpene lactones  that inhibit browsing by ungulates [24,25,26]. The greater sage-grouse feeds exclusively on evergreen sagebrush foliage in the winter and browses frequently on sagebrush during snow-free months . Kohl et al.  recently found that relative to domestic chickens, the cecal microbiota of the greater sage-grouse was enriched in genes of Bacteroides, Eggerthella and Clostridium associated with the metabolism of plant secondary compounds, especially phenolics . This suggests that sage-grouse rely heavily on specialized microbiota to cope with their toxic sagebrush diet. However, Kohl et al.  study was limited to cecal samples from three specimens collected during November–December, which precludes any conclusions on seasonal variation in microbiota across multiple gut regions. To date, the seasonal variation of the gut microbiota has yet to be studied in sage-grouse or other wild birds.
In this study we investigated the seasonal variation in the bacterial and archaeal communities (hereafter microbiota) of the gut of the greater sage-grouse. Using next-generation amplicon sequencing of the 16S rRNA gene, we sought to (i) characterize the regionalization of microbiota of the crop, ventriculus, duodenum, cecum, and colon within individual sage-grouse; (ii) evaluate the seasonal differences in microbiota associated with summer and winter diets in population samples; and (iii) determine whether the gut microbiota of females and males differed in composition and richness. Finally, we address the implications of seasonality and anatomical regionalization of gut microbiota on current and future microbiome studies.
Diversity and abundance of microbial taxa
Our initial dataset included 8,755,549 sequences obtained from 145 samples (Additional file 1). Five samples of the original 150 were excluded due to poor DNA or PCR yield. Filtering reduced the total number of sequences to 6,639,051. The number of filtered sequences obtained from single samples ranged from 1140 to 102,201 (Additional file 2) with a median of 50,232.
Filtered sequences from the pooled samples of crop, ventriculus, duodenum, cecum, and colon were assigned to 461 OTUs (Additional file 2). Four OTUs were classified as Archaea (Euryarchaeota; Thermoplasmata; Thermoplasmatales; [Thermoplasmatales]) and 457 as Bacteria. All OTUs were taxonomically classified in SILVA v128 16S rRNA gene reference database  to order, 93.7% were classified to family, 76.4% to genus, and 3.9% to species.
Ten bacterial phyla were present: Firmicutes was the richest (229 OTUs) and the most abundant (53.6% of the total CSS + Log2 OTU abundance) phylum. Other common phyla included Bacteroidetes (73 OTUs, 15.2% CSS + Log2 abundance), Actinobacteria (42 OTUs, 10.7%), and Proteobacteria, (65 OTUs, 10.1%, Additional file 3). Of 64 bacterial families, Ruminococcaceae (Firmicutes) was the richest (80 OTUs) and the most abundant (21.0% total CSS + Log2 abundance), followed by Lachnospiraceae (Firmicutes; 42 OTUs, 9.8%), Bacteroidaceae (Bacteroidetes; 31 OTUs, 8.1%), and Lactobacillaceae (Firmicutes; 28 OTUs, 6.4%; Additional file 3). Among the 142 bacterial genera, Ruminococcaceae UCG-014 was the richest and the most abundant (34 OTUs, 9.4%), followed by Bacteroides (31 OTUs, 8.1%) and Lactobacillus (28 OTUs, 6.4%, Additional file 3).
Regionalization of gut microbiota
No single OTU constituted more than 1.1% of all CSS + Log2 sequences in the pooled dataset (all gut regions of all birds; Additional file 3). However, most OTUs occurred in multiple gut regions: 308 of 461 OTUs (66.8%) were found in all five gut regions, 111 OTUs (24.1%) were found in four, 28 OTUs (6.1%) in three, 10 OTUs (2.2%) were observed in two gut regions, and only 4 OTUs (0.9%) were detected in a single gut region (Additional file 3).
The connectivity among different gut regions was further demonstrated by the presence of a positive abundance-occupancy relationship between the overall OTU CSS + Log2 abundance and the number of samples in which they were detected (N samples containing an OTU = 22.58 + 0.12 × [CSS + Log2 OTU count]; adjusted r2 = 0.816, df = 459, P < 2.2 × 1016). Although the majority of OTUs occurred in all gut regions, 387 OTUs showed significantly different CSS + Log2 abundances in different gut regions (Additional file 3).
Our PERMANOVA results indicated that gut region, season of collection, their interaction, latitude and longitude of collection localities had significant effects on microbiota composition (Table 1), whereas sex, body mass and their interactions with other explanatory variables did not. The PERMANOVA r2 value for the effect of gut region (r2 = 0.483) was eight times larger than the values for season (r2 = 0.055), or their interaction (r2 = 0.060), which in turn were six times larger than the values for longitude (r2 = 0.011) and latitude (r2 = 0.010). Gut region, season of collection, and their interaction remained the only significant variables or those with largest effects in PERMANOVAs for individual pairs of gut regions, with the exception of the interaction term in the ventriculus-cecum pair (Additional file 4). Effects of other variables and interactions were much smaller and only significant in some of the pairwise comparisons between gut regions. Season was the only significant variable or had the largest effect in PERMANOVAs for individual gut regions. The effect of the gut region in PERMANOVAs within season was more than 21 times greater than the effect of longitude, which was the only other variable that had a significant effect in both summer and winter samples (Additional file 4).
The clusters of data points for crop and cecum samples did not overlap in the bivariate PCoA plot and displayed limited overlap with those of the ventriculus, duodenum, and colon (Fig. 2). Pairwise PERMANOVA r2 values (Additional file 4) corresponded to the degree of overlap among gut region clusters in the PCoA plot (Fig. 2). Furthermore, LEfSe indicated that crop and cecum microbiotas had larger numbers of distinguishing taxa than other gut regions (Fig. 3, Additional file 5). Significantly higher abundances of the genera Lactobacillus, Mycoplasma, and unclassified genera of Pasteurellaceae and Leptotrichiaceae were likely responsible for the observed differences between microbiotas of the crop and other gut regions. The cecum microbiota was distinguished by overrepresentation of two Clostridia families: Ruminococcaceae and Lachnospiraceae, order Bacteroidales including Alistipes, Synergistaceae, Oxalobacter, and unclassified genera of Coriobacteriaceae and Flavobacteriaceae. LEfSe failed to identify any distinguishing taxa from the ventriculus microbiota. The duodenum microbiota was distinguished by an overrepresentation of Staphylococcus and unclassified genera of Erysipelotrichaceae and Veillonellaceae. The colon microbiota was distinguished by higher abundances of the phylum Actinobacteria, an unclassified genus of Lachnospiraceae, and Helicobacter.
Gut regions also differed in microbial richness (Fig. 4). The cecum had the richest microbiota (x̅ = 247.4 ± 12.7 OTUs per sample), followed by ventriculus (x̅ = 209.1 ± 52.3 OTUs per sample), colon (x̅ = 198.9 ± 69.4 OTUs per sample), duodenum (x̅ = 130.7 ± 54.0 OTUs per sample), and crop (x̅ = 88.0 ± 34.0 OTUs per sample). Linear mixed model regression accounting for the matched design by individual revealed that gut region (P < 2.200 × 10− 16) and season (P = 0.039) were significantly associated with microbiota richness (Fig. 4a) in a model also accounting for sex, longitude, latitude, and body mass (Additional file 6). Neither sex (P = 0.936), longitude (P = 0.339), latitude (P = 0.562), nor body mass (P = 0.829) were significantly associated with microbiota richness. Other α-diversity indices revealed similar patterns of variation among gut regions and seasons (Fig. 4b-d; Additional file 7).
Gut region was significantly associated with microbiota richness in all subsets restricted to pairs of gut regions or single seasons except the ventriculus - colon pair (Table 2). The season was significantly associated with microbiota richness in all five subsets that include duodenum (four gut region pairs and duodenum) and body mass in three of five subsets that include duodenum (crop – duodenum, duodenum – cecum, and duodenum). The effect of sex, geographic coordinates, or individual was not significantly associated with microbiota richness in any subsets.
The cecum samples displayed lower variability than samples from other gut regions. The size of the cecum cluster in the PCoA plot (Fig. 2) was much smaller than the clusters formed by samples from other gut regions. The average distance to centroid (Fig. 5) for the cecal samples (summer: 0.029, winter: 0.031) was significantly lower than the distances (summer: 0.051–0.066, winter: 0.055–0.076) for the other four gut regions during both seasons (Additional file 8). The average distances to centroid for the other four regions did not differ significantly from each other.
Seasonal variation of the greater sage-grouse microbiota
Microbial richness of the crop, duodenum, and cecum was significantly higher in the summer than in in the winter (Fig. 4). The significant interaction between the gut region and season (Table 1) indicated that season disparately affected the microbiota composition of different gut regions. The strength of seasonal effects (Table 1) was lower in the ventriculus (r2 = 0.141), duodenum (r2 = 0.186), and colon (r2 = 0.172) samples than for those from the cecum (r2 = 0.262) and crop (r2 = 0.321). PCoA plots showed that clusters of summer and winter samples were non-overlapping for all gut regions and average distances to centroid did not differ between seasonal clusters except for the duodenum (Fig. 6). However, the distance between centroids of the season-specific clusters was greater than the intra-cluster mean distances only for the crop and cecum.
LEfSe identified 109 genera whose abundance differed significantly between season in at least one gut region (Additional files 9 and 10). Twenty-eight bacterial genera, likely acquired from food, water, or grit, were relatively more abundant in summer samples, primarily in the crop (22 of 28 genera). Among host-associated bacteria, 36 genera were overrepresented in the summer and 38 in the winter. Four of 36 genera (Methanogranum, unclassified Clostridiales vadinBB60 group, Turicibacter, and unclassified Synergistaceae) were overrepresented in all five gut regions in the summer, whereas only two (Olsenella and unclassified Pasteurellaceae) of 38 genera were overrepresented in all five gut regions in the winter.
Our study provides robust evidence that the greater sage-grouse microbiota exhibits significant regionalization among gut regions. Gut region explained an order of magnitude more variance in microbiota composition than seasonal effects. A significant interaction effect between gut region and season revealed by our analyses reflects idiosyncratic patterns of seasonal changes in microbial communities of different gut regions. Microbiota richness also differed among gut regions and was significantly higher in the summer than in the winter in three of the five sampled gut regions. At the same time, latitude and longitude had limited effects on microbiota richness. We failed to identify any differences between the gut microbiota of females and males or a correlation between OTU richness and body mass of grouse specimens.
The overall composition of the sage-grouse gut microbiota appears similar to the microbiota of other birds [3, 5, 30, 31]. Firmicutes, Bacteroidetes, Actinobacteria, and Proteobacteria accounted for 90% of all CSS-normalized and Log2-transformed reads and 89% of all OTUs in our dataset. Unfortunately, the lack of standardized protocols in avian microbiota studies  and nearly exclusive use of fecal or cloacal samples  prevent detailed comparisons of our data with most of previously published studies. The few studies that have sampled multiple gut regions involved poultry or captive birds fed artificial diets: domestic chickens [32,33,34,35], captive Attwater’s prairie chickens Tympanuchus cupido attwateri , and farmed ostriches Struthio camelus . A recent study sampled three gut regions from the frozen carcasses of eight species of wild birds (a total of 32 individuals) in Venezuela . Uncertainties about the effects of freezing prior to microbial sampling, sampling protocols, and small sample sizes make it difficult to draw any comparisons with our sage-grouse results.
The best comparative data are provided in a recent study  of wild Canada geese (Branta canadensis) that employed the same sampling and lab protocols, sequencing platform, and raw data filtering used in our sage-grouse analysis. We found striking similarities in patterns of microbiota richness and regionalization among gut regions in sage-grouse and geese. We observed 461 OTUs in sage-grouse compared to 421 in geese. The cecum had the richest but least variable microbiota in both geese and sage-grouse. The considerably longer retention time of digesta in the cecum, relative to other gut regions [40, 41], permits cecal microbial communities to stabilize and is likely the cause for the reduced variability observed among individuals.
Despite the strong positive abundance-occupancy relationship and low proportion of OTUs restricted to one or two gut regions in sage-grouse and geese, gut regions exhibited significant differences in their microbiota composition. In both studies, the cecum and crop/esophagus had the most distinctive microbiotas. Pasteurellaceae and Leptotrichiaceae were overrepresented in the crop/esophagus, whereas Staphylococcaceae were overrepresented in the duodenum of geese and sage-grouse. Ruminococcaceae were overrepresented in the cecum of both species and so were at least some OTUs from Coriobacteriaceae, Bacteroidaceae, Prevotellaceae, and Rikenellaceae. Compositional differences among gut regions are likely related to their functional differences [40,41,42,43] which impose strong selection on microbiota despite the bidirectional flow of digesta among gut regions [40, 42, 44].
Kohl et al.  investigated the cecal microbiota of three individuals of greater sage-grouse with shotgun metagenomic sequencing. The authors provided little information on the taxonomic composition of the cecal microbiotas, but listed Bacteroides (19.2% of reads), Prevotella (9.6%), and Clostridium (9.1%) as the most abundant genera. Our sage-grouse cecal results for these genera (10.5, 0.4, 0.0%, respectively) were significantly different from those reported by Kohl et al.  but were similar to the abundances reported for Canada geese (9.5, 2.9, and 0.3%)  and Japanese rock ptarmigan (8.1, 2.1, 0.0%) . Furthermore, Arthrobacter (Actinobacteria: Micrococcaceae), which Kohl et al.  reported as the main source of genes degrading toxic phenols and catechols in the sage-grouse cecum, was not present in any of our sage-grouse cecal samples, nor were they detected in rock ptarmigan or geese ceca. These discrepancies highlight the need for additional metagenomic analyses of our samples.
Our study appears to be the first to evaluate seasonal microbiota changes in multiple gut regions of a wild vertebrate. Seasonal changes were most pronounced in the crop and cecum microbiotas. The crop is the most anterior of the sampled gut regions and serves as storage organ for consumed foliage before it is released into the ventriculus. We suspect that the pronounced seasonal changes of crop microbiota may be due to marked seasonal variation in the microbial communities of ingested foliage, water, soil, and arthropods as influenced by temperature and snow cover. The paired intestinal ceca are of crucial importance for nutrition in grouse (Aves: Tetraoninae) and are thought to be an adaptation for processing large quantities of poor quality foods . The finely-ground, soluble, and readily digestible forage is shunted to the ceca while large indigestible fragments are rapidly excreted from the gut. The ceca serve as the primary site for the reabsorption of salts and water, and the breakdown of complex carbohydrates (i.e., dietary fiber) and uric acid into volatile fatty acids and ammonia via microbial fermentation [40, 46, 47]. In contrast to other gut regions, ceca size fluctuates substantially in response to seasonal and even short-term changes in diet, increasing concomitantly with the increase of fiber content in forage [40, 46, 47] or when energy demands increase and food quality decreases . Our data suggest that cecal microbiota composition and richness change in response to seasonal diet variation to a greater degree than those of other gut regions distal to the crop.
Our findings add to a small but rapidly growing body of work that has demonstrated spatial structuring of microbial communities in functionally distinct gut regions of vertebrates: catfish , snakes , rodents , and birds [32,33,34, 36,37,38,39, 51]. Collectively, these analyses suggest that fecal samples alone may be poor predictors of total microbial diversity and abundance in the vertebrate gut. Our study contributes an additional caveat by discovering uneven seasonal changes in microbiota composition and richness among functionally distinct gut regions. These findings together with the independence of region-specific microbiota variation among individuals  strongly suggest that future studies of gut microbiota ecology and evolution should sample functionally distinct gut regions in addition to fecal or cloacal samples.
Greater sage-grouse specimens were collected in Wyoming at the end of phenological summer in Sublette County (19–22 September 2016, n = 15) and during phenological winter (16–20 December 2016) in Natrona (n = 5) and Sublette (n = 10) counties (Fig. 1). Detailed data on voucher specimens deposited in the National Museum of Natural History are provided in the Additional file 1.
Greater sage-grouse show marked sexual dimorphism in external measurements and body mass . We measured body mass with a digital scale to the nearest gram before microbiota sampling. Males (x̅ = 2284 ± 298 g; n = 14) included in this study were significantly heavier than the females (x̅ = 1226 ± 151 g; n = 16; Wilcoxon W = 0, P = 3.541 × 10− 6). Although we did not measure the size or weight of gut regions in individuals, we assumed that luminal volume of different gut regions (e.g., cecum) was correlated with body mass .
We sampled five gut regions: crop, ventriculus, duodenum, cecum, and colon . Specimens were put on ice soon after collecting and were kept on ice and processed the same day. We used sterile single-use polyester-tipped applicators with a plastic shaft (Fisher Scientific, Hampton, NH, USA) to sample microbiota from the crop, ventriculus, and colon. Applicators were inserted through a small incision made with a sterile scalpel blade in the crop and ventriculus. For colon sampling, applicators were inserted approximately 30 mm into the colon through the cloaca. Applicators were gently swirled in the lumen and then rubbed on the mucosal lining. Immediately after sampling, applicator ends were snapped off in sterile collection vials, capped, and submerged in liquid nitrogen.
The duodenum and cecum (one side) were sampled by isolating a ~ 25 mm section of the respective gut region with two sterile surgical hemostats. We then injected 1 mL of sterile water (HyPure molecular biology grade water, GE Healthcare Life Sciences, Logan, UT, USA) into the isolated section with a sterile single-use syringe. The water was drawn in and out of the syringe three times to ensure mixing of the intestinal contents. The resulting lavage mixture was immediately transferred into a sterile collection vial, capped, and submerged in liquid nitrogen. All samples were kept in liquid nitrogen until DNA extraction.
Prior to DNA extraction, vials were centrifuged at 2 × 104 G for 20 min and the supernatant was removed with a pipette using sterile filter tips. Total genomic DNA was extracted using the PowerSoil DNA Isolation Kit (Mo Bio Laboratories, Carlsbad, CA, USA) and cleaned with the PowerClean Pro DNA Clean-Up Kit (Mo Bio Laboratories, Carlsbad, CA, USA) to remove PCR inhibitors.
We amplified and sequenced the V4 region of the 16S rRNA gene (252–254 bp) from Bacteria and Archaea using primers 515FB and 806RB , following procedures outlined in the Earth Microbiome Project 16S protocol [54, 55]. Detailed description of PCR conditions, library preparation, and sequencing on an Illumina MiSeq (Illumina, Inc., San Diego, CA, USA) are provided in Drovetski et al. .
Illumina data processing
Raw Illumina data processing included joining of forward and reverse reads (join_paired_ends.py), demultiplexing (split_libraries_fastq.py), and quality filtering (Phred score Q ≥ 30, max barcode error = 0, min length = 200 bp) using the Quantitative Insights Into Microbial Ecology (QIIME) pipeline v1.9.0 .
We used UCHIME  to identify and remove chimeric sequences against the ChimeraSlayer reference database (version ‘microbiomeutil-r20110519’). The remaining sequences were aligned in PyNAST  and taxonomically classified using Bayesian RDP Classifier  trained with the SILVA v128 16S rRNA gene reference database . Aligned sequences were clustered into operational taxonomic units (OTUs) at 3% divergence [29, 60], GenBank accession numbers of sequences that matched OTUs in our dataset were retrieved from the SILVA v128 database. We classified OTUs as host-associated or likely obtained through environmental sources (e.g., water, soil, plants, arthropods, etc.) from ‘isolation source’ metadata in GenBank entries and references for accession numbers and their closest matches (99% identity).
Singleton sequences and OTUs with an overall relative abundance < 0.01% were filtered out to reduce the likelihood of sequence artifacts affecting downstream diversity metrics . We also excluded OTUs classified as Mitochondria, Eukaryota, Chloroplast, and those of unknown origin. We constructed an UPGMA tree of remaining OTU representative sequences in Geneious 11.1.4 (Biomatters Ltd., Auckland, New Zealand). We cumulative sum scaled  and Log2-transformed OTU abundances (CSS + Log2 OTU abundances) to account for variation in sequencing depth among samples  and non-normal distribution of abundances among OTUs in Calypso 8.10 . All statistical analyses were based on CSS + Log2 OTU abundances.
We constructed rarefaction plots (Additional file 12) with R package Vegan 2.4–4  to ensure our sequencing depth was sufficient and OTU accumulation curves in Microsoft Excel v14.7.7 to evaluate sampling effort of individual birds (Additional file 12). We used QIIME for α and β-diversity analyses. We calculated all available α-diversity indexes (alpha_diversity.py). We also calculated abundance-weighted evolutionary distinctiveness (AED)  in the R package BAT (https://cran.r-project.org/web/packages/BAT/). Differences in OTU abundances among gut regions were tested using non-parametric Kruskal-Wallis tests (group_significance.py).
We calculated weighted UniFrac distances  among individual samples and conducted PCoA (beta_diversity_through_plots.py) to compare microbiota composition. Weighted UniFrac distances were also used in the Permutational Multivariate Analysis of Variance (PERMANOVA)  implemented in R package Vegan 2.4–4  to test for the effect of the gut region (crop, ventriculus, duodenum, cecum, and colon), season (summer and winter), sex (male and female), geographic coordinates, body mass, and their interactions on microbiota variation among samples. Effect of the same variables on microbiota richness was evaluated using linear mixed model regression accounting for the matched design by individual grouse specimen. Their significance was determined by log likelihood tests.We used R version 3.3.3 (http://www.R-project.org) to generate Principal Coordinate Analysis (PCoA) and box plots. We modeled the abundance-occupancy relationship [69, 70] by regressing prevalence of OTUs on their total abundance. Paired Wilcoxon Signed Rank tests were used to compare richness among gut regions, body mass of males and females, and F-tests were used to compare variance of the PCoA scores among gut regions. Radar plots were made in Microsoft Excel v14.7.7.
We used the Linear Discriminant Analysis effect size (LEfSe) algorithm implemented in the LEfSe software package  to identify microbial taxa that best explained microbiota differences among gut regions and seasons. LEfSe consists of three consecutive steps: (i) Kruskal-Wallis tests to identify taxa with significantly different abundances among groups, (ii) a pairwise Wilcoxon Signed Rank test to select the subsets of taxa overrepresented in only a single group, and (iii) Linear Discriminant Analysis to estimate the effect size of each taxon.
Abundance-weighted evolutionary distinctiveness
Cumulative sum scaled
Linear discriminant analysis effect size
Operational taxonomic unit
Principal coordinates analysis
Polymerase chain reaction
Permutational multivariate analysis of variance
Ribosomal ribonucleic acid
Fraune S, Bosch TCG. Why bacteria matter in animal development and evolution. BioEssays. 2010;32(7):571–80.
Hanning I, Diaz-Sanchez S. The functionality of the gastrointestinal microbiome in non-human animals. Microbiome. 2015;3(1):51.
Colston TJ, Jackson CR. Microbiome evolution along divergent branches of the vertebrate tree of life: what is known and unknown. Mol Ecol. 2016;25(16):3776–800.
Hird SM. Evolutionary biology needs wild microbiomes. Front Microbiol. 2017;8:725.
Kohl KD. Diversity and function of the avian gut microbiota. J Comp Physiol B. 2012;182(5):591–602.
Smits SA, Leach J, Sonnenburg ED, Gonzalez CG, Lichtman JS, Reid G, et al. Seasonal cycling in the gut microbiome of the Hadza hunter-gatherers of Tanzania. Science. 2017;357(6353):802.
Wu Q, Wang X, Ding Y, Hu Y, Nie Y, Wei W, et al. Seasonal variation in nutrient utilization shapes gut microbiome structure and function in wild giant pandas. Proc R Soc B Biol Sci. 2017;284(1862):20170955.
Ren T, Boutin S, Humphries MM, Dantzer B, Gorrell JC, Coltman DW, et al. Seasonal, spatial, and maternal effects on gut microbiome in wild red squirrels. Microbiome. 2017;5(1):163.
Maurice CF, Cl Knowles S, Ladau J, Pollard KS, Fenton A, Pedersen AB, et al. Marked seasonal variation in the wild mouse gut microbiota. ISME J. 2015;9(11):2423–34.
Sun B, Wang X, Bernstein S, Huffman MA, Xia D-P, Gu Z, et al. Marked variation between winter and spring gut microbiota in free-ranging Tibetan Macaques (Macaca thibetana). Sci Rep. 2016;6:26035.
Carey HV, Walters WA, Knight R. Seasonal restructuring of the ground squirrel gut microbiota over the annual hibernation cycle. Am J Physiol Regul Integr Comp Physiol. 2013;304(1):R33–42.
Hammer TJ, Bowers MD. Gut microbes may facilitate insect herbivory of chemically defended plants. Oecologia. 2015;179(1):1–14.
Dickinson EC, Christidis L, editors. The Howard and Moore complete checklist of the birds of the world fourth edition, volume 2: passerines. Eastbourne: Aves Press; 2014.
Dickinson EC, Remsen JV, editors. The Howard and Moore complete checklist of the birds of the world fourth edition, volume 1: non-passerines. Eastbourne: Aves Press; 2013.
Schroeder MA, Young JR, Braun CE. Greater sage-grouse (Centrocercus urophasianus), version 2.0. In: Rodewald PG, editor. The birds of North America. Cornell Lab of Ornithology: Ithaca, New York, USA; 1999.
Shafizadeh F, Bhadane NR, Kelsey RG. Sesquiterpene lactones of sagebrush: constituents of Artemisia tripartita. Phytochemistry. 1974;13(3):669–70.
Welch BL, McArthur ED. Variation of monoterpenoid content among subspecies and accessions of Artemisia tridentata grown in a uniform garden. J Range Manag. 1981;34(5):380–4.
Kelsey RG. The chemical constituents of sagebrush foliage and their isolation. —. J Range Manag. 1982;35:617.
Wilt FM, Miller GC. Seasonal variation of coumarin and flavonoid concentrations in persistent leaves of Wyoming big sagebrush (Artemisia tridentata ssp. wyomingensis: Asteraceae). Biochem Syst Ecol. 1992;20(1):53–67.
Wilt FM, Geddes JD, Tamma RV, Miller GC, Everett RL. Interspecific variation of phenolic concentrations in persistent leaves among six taxa from subgenus Tridentatae of Artemisia (Asteraceae). Biochem Syst Ecol. 1992;20(1):41–52.
Kohl KD, Pitman E, Robb BC, Connelly JW, Dearing MD, Forbey JS. Monoterpenes as inhibitors of digestive enzymes and counter-adaptations in a specialist avian herbivore. J Comp Physiol B. 2015;185(4):425–34.
Kohl KD, Connelly JW, Dearing MD, Forbey JS. Microbial detoxification in the gut of a specialist avian herbivore, the Greater Sage-Grouse. FEMS Microbiol Lett. 2016;363(14).
Kelsey RG, Morris MS, Shafizadeh F. The use of sesquiterpene lactones as taxonomic markers in the shrubby species of Artemisia (section Tridentatae) in Montana. J Range Manag. 1976;29(6):502–5.
Olsen FW, Hansen RM. Food relations of wild free-roaming horses to livestock and big game, Red Desert, Wyoming. J Range Manag. 1977;30(1):17–20.
Johnson MK. Foods of primary consumers on cold desert shrub-steppe of southcentral Idaho. J Range Manag. 1979;32(5):365–8.
Hanley TA, Kathleen AH. Food resource partitioning by sympatric ungulates on Great Basin rangeland. J Range Manag. 1982;35(2):152–8.
Wallestad R, Eng RL. Foods of adult sage grouse in Central Montana. J Wildl Manag. 1975;39(3):628–30.
Selma MV, Espín JC, Tomás-Barberán FA. Interaction between Phenolics and gut microbiota: role in human health. J Agric Food Chem. 2009;57(15):6485–501.
Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, et al. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007;35(21):7188–96.
Waite D, Taylor M. Exploring the avian gut microbiota: current trends and future directions. Front Microbiol. 2015;6:673.
Waite DW, Taylor MW. Characterizing the avian gut microbiota: membership, driving influences, and potential function. Front Microbiol. 2014;5:223.
Wang L, Lilburn M, Yu Z. Intestinal microbiota of broiler chickens as affected by litter management regimens. Front Microbiol. 2016;7:593.
Choi JH, Kim GB, Cha CJ. Spatial heterogeneity and stability of bacterial community in the gastrointestinal tracts of broiler chickens. Poult Sci. 2014;93(8):1942–50.
Xiao Y, Xiang Y, Zhou W, Chen J, Li K, Yang H. Microbial community mapping in intestinal tract of broiler chicken. Poult Sci. 2017;96(5):1387–93.
Lu J, 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(11):6816–24.
Zhang Y, Simon SE, Johnson JA, Allen MS. Spatial microbial composition along the gastrointestinal tract of captive Attwater’s prairie chicken. Microb Ecol. 2016;73(4):966–977.
Videvall E, Strandh M, Engelbrecht A, Cloete S, Cornwallis CK. Measuring the gut microbiome in birds: comparison of faecal and cloacal sampling. Mol Ecol Resour. 2018;18(3):424–34.
García-Amado MA, Shin H, Sanz V, Lentino M, Martínez LM, Contreras M, et al. Comparison of gizzard and intestinal microbiota of wild neotropical birds. PLoS One. 2018;13(3):e0194857.
Drovetski SV, O’Mahoney M, Ransome EJ, Matterson KO, Lim HC, Chesser RT, et al. Spatial organization of the gastrointestinal microbiota in urban Canada geese. Sci Rep. 2018;8(1):3713.
Svihus B. Function of the digestive system. The J Applied Poultry Res. 2014;23(2):306–14.
Roto SM, Rubinelli PM, Ricke SC. An introduction to the avian gut microbiota and the effects of yeast-based prebiotic-type compounds as potential feed additives. Frontiers in Veterinary Sci. 2015;2:28.
Duke GE. Gastrointestinal physiology and nutrition in wild birds. Proc Nutr Soc. 1997;56(3):1049–56.
Pan D, Yu Z. Intestinal microbiome of poultry and its interaction with host and diet. Gut Microbes. 2014;5(1):108–19.
Remington TE. Why do grouse have ceca? A test of the fiber digestion theory. J Exp Zool. 1989;252(S3):87–94.
Ushida K, Segawa T, Tsuchida S, Murata K. Cecal bacterial communities in wild Japanese rock ptarmigans and captive Svalbard rock ptarmigans. J Vet Med Sci. 2016;78(2):251–7.
Svihus B, Choct M, Classen HL. Function and nutritional roles of the avian caeca: a review. World’s Poultry Sci J. 2013;69(2):249–64.
Clench MH, Mathias JR. The avian cecum: a review. The Wilson Bulletin. 1995;107(1):93–121.
McDonald R, Schreier HJ, Watts JEM. Phylogenetic analysis of microbial communities in different regions of the gastrointestinal tract in Panaque nigrolineatus, a wood-eating fish. PLoS One. 2012;7(10):e48018.
Colston TJ, Noonan BP, Jackson CR. Phylogenetic analysis of bacterial communities in different regions of the gastrointestinal tract of Agkistrodon piscivorus, the cottonmouth snake. PLoS One. 2015;10(6):e0128793.
Kohl KD, Dearing MD, Bordenstein SR. Microbial communities exhibit host species distinguishability and phylosymbiosis along the length of the gastrointestinal tract. Mol Ecol. 2018;27(8):1874–83.
Han GG, Kim EB, Lee J, Lee J-Y, Jin G, Park J, et al. Relationship between the microbiota in different sections of the gastrointestinal tract, and the body weight of broiler chickens. Springerplus. 2016;5(1):911.
Franz R, Hummel J, Kienzle E, Kölle P, Gunga H-C, Clauss M. Allometry of visceral organs in living amniotes and its implications for sauropod dinosaurs. Proc R Soc B Biol Sci. 2009;276(1662):1731–6.
McLelland J. Apparatus digestorius [systema alimentarium]. In: Baumel, JJ, King, AS, Breazile, JE, Evans, HE, and Vanden Berge, JC. (eds). Handbook of avian anatomy: nomina anatomica avium. Publications of the Nuttall Ornithological Club (USA). no 23. Second edition. Cambridge: The Nuttall Ornithological Club; 1993. p. 301–327.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6(8):1621–4.
EMP. Earth Microbiome Project: 16S Illumina amplicon protocol 2017 [cited 2017 Accessed 17 May]. Available from: http://press.igsb.anl.gov/earthmicrobiome/protocols-and-standards/16s/.
Caporaso JG. Qime allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1.
Caporaso JG, Bittinger K, Bushman FD, DeSantis TZ, Andersen GL, Knight R. PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics. 2010;26(2):266–7.
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41(Database issue):D590–D6.
Bokulich NA, Subramanian S, Faith JJ, Gevers D, Gordon JI, Knight R, et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Meth. 2013;10(1):57–9.
Paulson JN, Stine OC, Bravo HC, Pop M. Robust methods for differential abundance analysis in marker gene surveys. Nat Methods. 2013;10(12):1200–2.
McMurdie PJ, Holmes S. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol. 2014;10(4):e1003531.
Zakrzewski M, Proietti C, Ellis JJ, Hasan S, Brion M-J, Berger B, et al. Calypso: a user-friendly web-server for mining and visualizing microbiome–environment interactions. Bioinformatics. 2017;33(5):782–3.
Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: community ecology package. R package version 2.4–3. 2017. Available from: https://CRAN.R-project.org/package=vegan.
Cadotte MW, Jonathan Davies T, Regetz J, Kembel SW, Cleland E, Oakley TH. Phylogenetic diversity metrics for ecological communities: integrating species richness, abundance and evolutionary history. Ecol Lett. 2009;13(1):96–105.
Lozupone CA, Hamady M, Kelley ST, Knight R. Quantitative and qualitative β diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microbiol. 2007;73(5):1576–85.
Anderson MJ. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26(1):32–46.
Gaston KJ, Blackburn TM, Greenwood JJD, Gregory RD, Quinn RM, Lawton JH. Abundance–occupancy relationships. J Appl Ecol. 2000;37:39–59.
Gaston KJ. The multiple forms of the interspecific abundance-distribution relationship. Oikos. 1996;76:211–20.
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60.
We thank the Wyoming Game and Fish Department for the research permit to sample birds in December 2016. We are grateful to Brandon Scurlock, Hilda Sexauer, Darren Rhea, Joe Bohne, and Justin Binfet (all current or former employees of Wyoming Game and Fish Department), Oscar “Boomer” Freesen (Pinedale), Jacob Saucier (Smithsonian Institution), Paul Rhymer, and dogs (Oscar, Buzz, Karoo, Husker) for contributed field expertise and assistance. All laboratory and computer work were conducted in and with the support of the Laboratories of Analytical Biology (LAB) facilities of the National Museum of Natural History, Smithsonian Institution. GRG acknowledges the continuing support of the Smoketree Trust. Two anonymous reviewers provided valuable comments that helped streamline statistical analyses and improve manuscript clarity.
This study was funded by the Shealy family through the Basis Foundation (https://thebasisfoundation.weebly.com). The funding body played no role in the design of the study, collection, analysis, and interpretation of data, or in writing the manuscript.
Availability of data and materials
The raw sequences archive has been deposited to GenBank: BioProject ID PRJNA478431. GenBank biosample accession numbers for individual samples are listed in Additional file 1. All other data generated or analyzed during this study are included in this published article (and its Additional files 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12).
This project was approved by the Institutional Animal Care and Use Committee (permit 2016–06) of the National Museum of Natural History (Smithsonian Institution). All parts of this study were conducted in accordance with guidelines and regulations of the Wyoming Game and Fish Department (research permit 1085 and hunting licenses).
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.
Voucher specimens deposited at the National Museum of Natural History (Smithsonian Institution), sample IDs for different gut regions, sex, age, and locality data for birds used in this study. (XLSX 57 kb)
Distribution of raw OTU counts among gut regions and individual samples. (XLSX 323 kb)
Taxonomy and distribution of CSS-normalized and Log2-transformed OTU counts among gut regions and individual samples. (XLSX 491 kb)
Results of PERMANOVAs for the effect of the season, sex, gut region, longitude, latitude, body mass, and their interactions on the weighted UniFrac distances among samples in different subsets of data. (XLSX 58 kb)
LEfSe scores and FDR-corrected P-values for microbial taxa most likely responsible for differences among gut regions. (XLSX 15 kb)
Results of the log likelihood tests of the mixed linear models for the association between OTU richness and season, sex, gut region, longitude, latitude, body mass, and individual in different subsets of data. (XLSX 21 kb)
Alpha diversity indices, season, sex, gut region, longitude, latitude, and body mass for individual samples. (XLSX 68 kb)
Results of Tukey’s multiple comparison tests for β-diversity dispersion (differences in average distance to multidimensional median) among gut regions in different seasons. (PDF 44 kb)
Radar plots summarizing LEfSe results at the genus level. Each radial line represents a genus significantly overrepresented in summer (top) and winter (bottom) microbiota. The scale and distance from the center of the plot to a data point represents the Log10 (LDA score). Point colors represent different gut regions according to the legend. Blue font of the generic names indicates environmental genera. Gray font shows genera with different seasonal assignments in different gut regions. Bold black font identifies genera with consistent seasonal differences in all five gut regions. Regular black font identifies genera with seasonal abundance differences in at least a single gut region. LDA scores for individual OTUs and nested higher-level taxa are presented in Additional file 10. (PDF 369 kb)
LEfSe scores and FDR corrected P-values for microbial taxa most likely responsible for differences between summer and winter in each gut region. (XLSX 66 kb)
Rarefaction plots for individual samples in each gut region. (TIF 2081 kb)
OTU accumulation plots. Individual birds were added in the order they were sampled. (PDF 124 kb)