Spatiotemporal variation in the fecal microbiota of mule deer is associated with proximate and future measures of host health
BMC Veterinary Research volume 17, Article number: 258 (2021)
Mule deer rely on fat and protein stored prior to the winter season as an energy source during the winter months when other food sources are sparse. Since associated microorganisms (‘microbiota’) play a significant role in nutrient metabolism of their hosts, we predicted that variation in the microbiota might be associated with nutrient storage and overwintering in mule deer populations. To test this hypothesis we performed a 16S rRNA marker gene survey of fecal samples from two deer populations in the western United States before and after onset of winter.
PERMANOVA analysis revealed the deer microbiota varied interactively with geography and season. Further, using metadata collected at the time of sampling, we were able to identify different fecal bacterial taxa that could potentially act as bioindicators of mule deer health outcomes. First, we identified the abundance of Collinsella (family: Coriobacteriaceae) reads as a possible predictor of poor overwintering outcomes for deer herds in multiple locations. Second, we showed that reads assigned to the Bacteroides and Mollicutes Order RF39 were both positively correlated with deer protein levels, leading to the idea that these sequences might be useful in predicting mule deer protein storage.
These analyses confirm that variation in the microbiota is associated with season-dependent health outcomes in mule deer, which may have useful implications for herd management strategies.
Mule Deer (Odocoileus hemionus) are indigenous to the western United States and Canada. Populations have experienced significant fluctuations over the past one-half century and there is evidence of an overall decline in numbers . Many factors contribute to fluctuations in herd size including maternal investment in reproduction, predation, disease, accidents (i.e.., road kill), winter severity, and others [2,3,4,5,6]. However, it is likely that the primary driver of population dynamics is herd health or the availability of sufficient food for growth, maintenance, and reproduction . In fact, condition (e.g., ingesta free body fat (IFBF) of female mule deer is linked to life-history characteristics including over-winter survival, birthweight of offspring, and survival of offspring [8,9,10,11]. Understanding the effects of management decisions on health or condition of mule deer are major initiatives of governmental agencies and non-governmental organizations throughout the range of this species.
Fat acquisition and storage is a crucial aspect of individual survival for mule deer and other ungulates. With lower access to energy and nutrients during the winter months, deer first utilize fat stores and then protein mass for energy during winter. For example, fat levels may increase from between 9 to 25% from June to December, with a then steady decrease from January to April . This pattern of fat deposition in deer is well-established and tied to the scarcity of forage during winter months. In addition to food scarcity, seasonal changes in mule deer fat deposition may also be related to the microorganisms associated with the deer, especially in the gastrointestinal tract (the ‘microbiota’). The gut microbiota of ruminants, including deer, can directly impact the health of an individual deer and variation in the identity and abundance of specific microoorganisms can be affected by differences in host genotype, geography, age, gender, and diet ([13,14,15,16,17,18]). Ruminants depend on their gut microbiota for up to 70% of the energy extracted from their diet  and 95% of the rumen microbiome is made up of bacteria, mostly anaerobes or facultative anaerobes . Two of the most common bacterial phyla in ruminant guts are involved in digestion: the Firmicutes break down fibers such as cellulose and the Bacteroidetes digest proteins and polysaccharides . Ruminants lack several enzymes necessary for digestion and it is believed that gut symbionts perform these functions. Thus, the microbiota is essential for deer health and metabolism. Currently there are no known microbial markers that directly link the gut microbiota and mule deer overwintering success.
Our purpose was to test if the deer gastrointestinal microbiota varies with season (winter vs. spring) and to identify specific microorganisms, if any, that are related to individual deer survival. We hypothesized that the microbiota would vary significantly with season and expected that specific microorganisms that varied in abundance, if detected, could be potential bioindicators of mule deer health. As part of the approach we asked two specific questions about the microbiota of two Utah mule deer herds: 1) Does the microbiota of individual deer vary with geography and season? and 2) Can variation in the microbiota predict proximate or future measures of deer health?
The microbiota of two Utah mule deer populations varies with geography and season
To better understand how location and season contribute to differences in the deer microbiota, we performed a 16S rRNA gene survey of 108 deer samples collected from herds located in Monroe (n = 48) and Cache (n = 60) valley (location) during December and March (season) (Table 1). A total of 3,404,323 bacterial reads and 4266 amplicon sequence variants (ASVs) were obtained on a partial Illumina HiSeq lane, with a median of 24,621 reads per sample. The most obvious differences in the microbiota were the increased abundance of Coriobacteriales in the Cache herd versus the Monroe herd (Fig. 1, Supplementary Tables 1–2). To identify factors associated with differences in the overall microbiota composition the microbiota of the different populations was compared using principal coordinate analysis (PCoA) of Unifrac distances in 13,000-read rarefied samples (Fig. 2A-B), which was near-saturating sampling (Figure S1). Microbiota composition varied significantly with both geography and season by identity of the associated microbes, and only by geography when microbial abundance was considered (Fig. 2, Supplementary Table 3). The significant interaction term in the unweighted Unifrac distance analysis suggests that the microbiota responses to seasonal variation were geography specific (p < 0.05). Together, these data showed a significant effect of both geography and season on the microbiota of mule deer.
We defined the core microbiota in different groups of sampled animals. Some bacterial phyla were present in every sample, including Bacteroidetes, Firmicutes, Actinobacteria, Verrucomicrobia, Cyanobacteria, and Proteobacteria. No ASVs were detected in every individual, but 4–12 core ASVs were detected in 90% of the animals in the different pools of animals grouped by season and geography (Table 1). In particular, samples from Cache Valley in both December and March showed multiple Coriobacteriaceae ASVs that were not present in the core microbiota of the Monroe Valley deer. As a whole, a search for a core microbiota revealed uniformity at the phylum level and ASV-level distinctions based on geographic location and season.
Coriobacteriaceae are associated with post-winter fat levels
To test if certain microbes might predict mule deer health, we calculated the relationship between the relative abundances of ASVs detected in deer in December and three health and nutritional measures in the same deer approximately 3 months later: rump fat, loin thickness, and body condition score (BCS) (Supplementary Tables 4–8). No bacterial ASVs or higher-level designations were significantly associated with variation in any of these three traits. However, when immature deer (< 2 years old) were removed from the analysis, read counts assigned to an unspecified species in the genus Collinsella, family Coriobacteriaceae, (Fig. 3) in December deer were positively associated with rump fat levels of the same individuals 3 months later (Fig. 3, FDR-corrected p < 0.05). This finding suggests that this ASV from the family Coriobacteriaceae may be a bioindicator of fat storage and consequently linked to outcomes such as survival and reproduction that are associated with increased levels of deer rump fat.
Individual ASVs are associated with mule deer protein storage
In addition to future measure of mule deer condition, we also tested if variation in the microbiota was associated with proximate measures of mule deer nutritional health. Spearman rank tests were used to test for correlations between microbial ASVs and deer metadata collected at the same time as the fecal samples. The only correlations that were significant after false-discovery rate correction were for mule deer loin thickness (Supplementary Tables 9–13). The lone correlated ASV was of relatively low abundance (88 reads in the full dataset) and could not be assigned past the family level as a member of the Ruminococcaceae (Fig. 4A). This correlation would not be useful to predict loin thickness from ASV content since the maximum value of the trendline was based on fewer than four of 13,000 reads sampled per deer. Therefore, the analysis was extended to ASVs grouped at higher taxonomic levels, based on the expectation that functional redundancy among conphyletic taxa might mask the detection of significant correlations in ASVs. At higher levels, 243 ASVs belonging to the order RF39 (phylum Tenericutes, class Mollicutes) drove significant positive correlations between loin thickness (a surrogate protein measure) at the order and phylum levels (Fig. 4B, see also Supplementary Tables 9–13). Taken together, these findings suggest that total RF39 ASV abundance in freshly collected fecal pellets is correlated with and can be used to predict proximate measures of protein storage in mule deer across multiple geographies and seasons.
Our goals were to survey the gastrointestinal microbiota of wild mule deer populations in Utah and test for possible relationships between the detected microbial communities and measures of mule deer health. A survey of two different populations of deer in Cache County and Monroe Valley, UT, USA revealed differences in the microbiota between these two herds that varied with the sampling time. Key biological differences between the two herds both before and after winter included fat levels, protein levels, and BCS, suggesting a relationship between herd location and health, accompanied by geographic variation in the microbiota. The identitification of microbial taxa whose abundance was correlated with individual post-winter rump fat levels provides an additional link between microbiota composition and herd health. However, as our data are insufficient to establish which factors are causal, we cannot infer if the associated microbes are responsible for or respond to the differences in BCS; or how deer behaviors (feeding) or environment (diet) may contribute to the observation variation in the microbiota.
Sequencing reads assigned to Collinsella, a Coriobacteriaceae genus, predicted future rump fat levels in the deer in this study, suggesting this taxon may be associated with condition of wintering animals. Coriobacteriaceae are known lipid metabolizers , and one explanation for the association we detected could be that these bacteria grow better in deer that begin the winter with high fat stores. An alternative view is that if Coriobacteriaceae metabolize lipids then their abundance should be negatively, not positively associated with host fat levels. Instead, the correlation may be related to bacterial utilization of other (e.g. dietary) fat resources, and their greater abundance in some animals may reflect that animals’ consumption of a diet elevated in fat. In either scenario, the data suggest that pre-winter sampling for Coriobactericeae abudances from deer droppings might provide useful insights to predicted deer performance after the winter season and aid in herd management plans by predicting populations likely to require overwintering interventions. One limitation of using sampled feces to predict future rump fat levels is that the correlation was not detected when deer of young ages were included in the analysis, and fecal pellets from 1 to 2 year old deer are not readily distinguished from adult fecal pellets. Overall, further research is necessary to fully test the idea that fecal pellet sampling can be used to predict deer overwintering and to define the relationship between Coriobacteriaceae and deer fat-content.
As loin thickness reflects protein content in individual deer , its correlation with RF39 and Ruminoccocaeceae abundance suggests bacteria belonging to these taxa may be related to protein storage or utilization. In other ruminants, RF39 abundance has been positively associated with residual feed intake (RFI), which measures the efficiency of an animal’s feed consumption, and gross feed intake (GFE), which reports an animal’s milk output normalized to dry matter intake [24, 25]. Larger RFI values are associated with excessive food intake and high GFE values suggest increased milk output per unit of food. If these bacteria have the same effects in mule deer, then the positive correlation between RF39 and loin thickness could be associated with increased accumulation or decreased utilization of protein storage, such as through excess feed consumption or increased efficiency. Alternatively, these microbes might causally shape the change in mule deer protein content. Regardless of the basis for the effect, our intitial analysis suggests fecal RF39 counts are a readily accessible bioindicator of mule deer loin thickness. Comprehensive support for this idea awaits further replication, where deer that were not used to calculate the regression are tested simultaneously for loin thickness and microbiota content, and matched for their fit to predictions.
An ideal goal is to economize time and finances by incorporating bioindicators of mule deer health or overwintering outcomes into established herd management strategies. Historical strategies for monitoring deer health included checking zyphoid fat levels at deer check stations or from roadkill animals, which provide relatively poor estimates of condition . More recently, a reliable method of measuring body fat has been developed , but the new method requires researchers to capture live animals and estimate condition based on body palpation and multiple measures acquired via ultrasonography. Because of the expense to capture life animals and the expertise required to correctly palpate animals and to make ultrasonography measures, it is difficult to implement such methods at a large scale. In this study, we show that sequencing the microbiota in fresh feces from wild-caught deer identifies taxa that predict proximate and future outcomes of deer. A key remaining gap is to test if older, wild-collected droppings can predict deer health outcomes with similar patterns to fresh samples. If fecal samples that are several days old retain predictive signatures (either those identified here, or possibly distinct ones), it points to the success of this approach.
In summary, this study surveyed the gut microbiota of mule deer in two locations before and after the winter season. The significant variation in mule deer microbiota composition with geography suggests environmental or dietary characters may contribute to differences in the mule deer microbiota. In addition, when the abundances of specific microbes were different between groups or significantly correlated with mule deer health outcomes, these findings identify ASVs that could putatively function as bioindicators of mule deer health. Further studies that clarify the causal basis for the correlative relationship between these taxa and deer health have the potential to provide further insight into the relationship between season, mule deer fat and protein storage, and the microbiota. Such work is important to better understand narrow issues addressed by this work, such as if mule deer feces can predict ot report mule deer health, as well as broader areas including the relationship between the microbiota and their animal hosts.
Capture and handling procedures were reviewed, evaluated, and approved by by the Brigham Young University’s (BYU) Animal Care and Use Committee (protocol 150,110). All capture and handling procedures complied with the approved methods. Samples were obtained at two different times from each of the two herds in Cache and Monroe valley, once in December and once in March, for 108 total samples. Individual deer sampled from the herds were randomly selected, so that some of the deer were captured twice (see the ‘Animal_ID’ column of Supplementary Table 2; i.e. 75 unique animals were captured of which 32 were captured twice). Measurements and samples were only collected for captured female deer. During handling, deer were weighed and age was estimated via tooth wear and eruption pattern [27, 28]. Preliminary results suggest > 80% accuracy within 2 years when aging mule deer in Utah from tooth wear compared to cementum analysis (Hinton, unpublished data). For condition of mule deer, ingesta free body fat (IFBF) was estimated using the BCS method and equations developed by Cook et al. . This method estimates IFBF from measurements of rump fat via ultrasonography (E.I. Medical Imaging portable ultrasound), estimates of age, body mass, and a BCS from 1 to 6 . In addition, loin thickness was estimated as suggested by Cook et al. . Fecal samples were collected directly from each animal with a gloved hand, and samples from each animal were stored in individual plastic bags and a cooler prior to freezing at − 20 °C. Captures took place in remote locations throughout the State of Utah. Nearly all animals were released near the site of capture. The only animals that were euthanized (~ 2%) were injured during capture such that they could not be released back into the wild. Those individuals were euthanized by authorized State of Utah personnel under the supervision of the State of Utah Wildlife Veterinarian. The procedure typically consisted of shooting the animal with a 0.22 long-rifle mushroom bullet such that the bullet would pass through the brain toward the spine.
Molecular biology and sequencing
For each of the 108 fecal samples, microbial DNA was isolated from 60 to 80 mg of a homogenized fecal pellet using the ZR-96 Fecal DNA Kit™ (Zymo Research, Irvine, CA; the kit includes a homogenization step) according to manufacturer instructions. 16S rRNA marker gene libraries (V4 region) were prepared exactly as described previously . Briefly, the V4 region was amplified using the AccuPrime Pfx SuperMix (Invitrogen, Carlsbad, CA) in combination with dual-barcoded primers (barcodes in Supplementary Table 2; custom (not from Illumina) sequencing primers were used, as described in ). PCR amplicons were normalized using the SequalPrep Normalization kit (ThermoFisher Scientific, Waltham, MA) according to manufacturer instructions. Sequencing was performed on a partial lane of a HiSeq 2500 at the BYU DNA Sequencing center. Samples plus 10% PhiX control DNA were sequenced using 2x250bp v2 Illumina sequencing kits. Sequences were deposited to the National Center for Biotechnological Information’s Short Read Archive under study number PRJNA743316 at https://www.ncbi.nlm.nih.gov/bioproject/PRJNA743316.
Demultiplexed samples were obtained from a partial lane of a paired-end 2 × 250 Illumina HiSeq 2500 and analyzed using QIIME2 . Reads were imported into a paired-end QIIME2 artefact, cleaned via denoising and chimera removal with the DADA2 plug-in , and clustered into ASVs. Taxonomy was assigned using pre-fitted, sklearn-based, Naïve Bayes classifier  trained on the GreenGenes reference base [34, 35] and a phylogenetic tree was built using the QIIME2 Fast Tree phylogeny pipeline . ASV tables were filtered to exclude ASVs assigned to Archaea. Samples were rarefied to 13,000 reads (Fig. 1) and unweighted and weighted Unifrac distances were used to test for statistically significant differences between treatments. The complete ASV table and taxonomic assignments are shown in Supplementary Table 1.
Statistical analyses were performed using QIIME2 and R . Beta diversity was calculated using unweighted and weighted Unifrac distance  and differences between samples were confirmed by PERMANOVA using the R package vegan . Differences in ASV abundance between samples were performed using ANCOM . Spearman rank correlations and Kruskal-wallis test for categorical variables were calculated in R between the values of traits estimating deer health (Rump fat, loin thickness, BCS) and the absolute abundance of bacterial ASVs in a 13,000-read rarefied ASV table. The cutoff for p-value significance was 0.05 or, where applicable, a false discovery rate (FDR)-corrected p-value of 0.05.
Availability of data and materials
Raw data are included as supplementary material or are available under accession numbers listed in the methods section.
ribosomal ribonucleic acid
Ingensta free body fat
Principal coordinate analysis
Amplicon sequence variant
Residual feed intake
Gross feed intake
Brigham Young University
Body condition score
False discovery rate
Degrees of freedom
Sum of squares
Mean sum of squares
- RS :
deVos JC, Conover MR, Headrick NE. Mule deer conservation: Issues and management strategies. Logan: Jack H. Berryman Institute Press; 2003.
Ballard WB, Lutz D, Keegan TW, Carpenter LH, deVos JC Jr. Deer-predator relationships: a review of recent north American studies with emphasis on mule and black-tailed deer. Wildl Soc Bull. 2001;29:99–115.
Bishop CJ, Unsworth JW, Garton EO. Mule deer survival among adjacent populations in Southwest Idaho. J Wildl Manag. 2005;69(1):311–21. https://doi.org/10.2193/0022-541X(2005)069<0311:MDSAAP>2.0.CO;2.
Unsworth JA, Pac DF, White GC, Bartmann RM. Mule deer survival in Colorado, Idaho, and Montana. J Wildl Manag. 1999;63(1):315–26. https://doi.org/10.2307/3802515.
Lukacs PM, White GC, Watkins BE, Kahn RH, Banulis BA, Finley DJ, et al. Separating components of variation in survival of mule deer in Colorado. J Wildl Manag. 2009;73(6):817–26. https://doi.org/10.2193/2008-480.
Hurley MA, Unsworth JW, Zager P, Hebblewhite M, Garton EO, Montgomery DM, et al. Demographic response of mule deer to experimental reduction of coyotes and mountain lions in southeastern Idaho. Wildl Monogr. 2011;178(1):1–33. https://doi.org/10.1002/wmon.4.
Bishop CJ, White GC, Freddy DJ, Watkins BE, Stephenson TR. Effect of enhanced nutrition on mule deer population rate of change. Wildl Monogr. 2009;172:1–28. https://doi.org/10.2193/2008-107.
Bender LC, Lomas LA, Browning J. Condition, survival, and cause-specific mortality of adult female mule deer in north-Central New Mexico. J Wildl Manag. 2007;71(4):1118–24. https://doi.org/10.2193/2006-226.
Tollefson TN, Shipley LA, Myers WL, Dasgupta N. Forage quality’s influence on mule deer fawns. J Wildl Manag. 2011;75(4):919–28. https://doi.org/10.1002/jwmg.113.
Tollefson TN, Shipley LA, Myers WL, Keisler DH, Dasgupta N. Influence of summer and autumn nutrition on body condition and reproduction in lactating mule deer. J Wildl Manag. 2010;74(5):974–86. https://doi.org/10.2193/2008-529.
Monteith KL, Bleich VC, Stephenson TR, Pierce BM, Conner MM, Kie JG, et al. Life-history characteristics of mule deer: effects of nutrition in a variable environment. Wildl Monogr. 2014;186(1):1–62. https://doi.org/10.1002/wmon.1011.
Worden K, Pekins P. Seasonal change in feed intake, body composition, and metabolic rate of white-tailed deer. Can J Zool. 2011;73:452–7.
Guan Y, Yang H, Han S, Feng L, Wang T, Ge J. Comparison of the gut microbiota composition between wild and captive sika deer (Cervus nippon hortulorum) from feces by high-throughput sequencing. AMB Express. 2017;7(1):212. https://doi.org/10.1186/s13568-017-0517-8.
Li Z, Zhang Z, Xu C, Zhao J, Liu H, Fan Z, et al. Bacteria and methanogens differ along the gastrointestinal tract of Chinese roe deer (Capreolus pygargus). PLoS One. 2014;9(12):e114513. https://doi.org/10.1371/journal.pone.0114513.
Li Y, Hu X, Yang S, Zhou J, Zhang T, Qi L, et al. Comparative analysis of the gut microbiota composition between captive and wild Forest musk deer. Front Microbiol. 2017;8:1705. https://doi.org/10.3389/fmicb.2017.01705.
Zhao G, Ma T, Tang W, Li D, Mishra SK, Xu Z, et al. Gut microbiome of Chinese Forest musk deer examined across gender and age. Biomed Res Int. 2019;2019:9291216.
Wang L, Ding J, Yang Z, Chen H, Yao R, Dai Q, et al. Pere David's deer gut microbiome changes across captive and translocated populations: implications for conservation. Evol Appl. 2019;12(3):622–35. https://doi.org/10.1111/eva.12743.
Hu X, Liu G, Li Y, Wei Y, Lin S, Liu S, et al. High-throughput analysis reveals seasonal variation of the gut microbiota composition within Forest musk deer (Moschus berezovskii). Front Microbiol. 2018;9:1674. https://doi.org/10.3389/fmicb.2018.01674.
Flint HJ, Bayer EA. Plant cell wall breakdown by anaerobic microorganisms from the mammalian digestive tract. Ann N Y Acad Sci. 2008;1125(1):280–8. https://doi.org/10.1196/annals.1419.022.
Puniya AK, Singh R, Kamra DN. Rumen microbiology: from evolution to revolution. New Delhi: Springer (India) Private Limited; 2015.
Flint HJ, Scott KP, Duncan SH, Louis P, Forano E. Microbial degradation of complex carbohydrates in the gut. Gut Microbes. 2012;3(4):289–306. https://doi.org/10.4161/gmic.19897.
Clavel T, Lepage P, Charrier C. The Family Coriobacteriaceae. In: Rosenberg E, EF DL, Lory S, Stackebrandt E, Thompson F, editors. The Prokaryotes: Actinobacteria. Berlin: Springer Berlin Heidelberg; 2014. p. 201–38.
Bender L, Cook J, Cook R, Hall P. Relations between nutritional condition and survival of north American elk Cervus elaphus. Wildl Biol. 2009;14:70–80.
Jami E, White BA, Mizrahi I. Potential role of the bovine rumen microbiome in modulating milk composition and feed efficiency. PLoS One. 2014;9(1):e85423. https://doi.org/10.1371/journal.pone.0085423.
Omoniyi LA, Jewell KA, Isah OA, Neumann AP, Onwuka CF, Onagbesan OM, et al. An analysis of the ruminal bacterial microbiota in west African dwarf sheep fed grass- and tree-based diets. J Appl Microbiol. 2014;116(5):1094–105. https://doi.org/10.1111/jam.12450.
Cook RC, Cook JG, Stephenson TR, Myers WL, Mccorquodale SM, Vales DJ, et al. Revisions of rump fat and body scoring indices for deer, elk, and moose. J Wildl Manag. 2010;74(4):880–96, 817. https://doi.org/10.2193/2009-031.
Robinette WL, Jones DA, Rogers G, Gashwiler JS. Notes on tooth development and wear for Rocky Mountain mule deer. J Wildl Manag. 1957;21(2):134–53. https://doi.org/10.2307/3797579.
Severinghaus CW. Tooth development and wear as criteria of age in white-tailed deer. J Wildlife Manag. 1949;13(2):195.
Cook RC, Stephenson TR, Myers WL, Cook JG, Shipley LA. Validating predictive models of nutritional condition for mule deer. J Wildl Manag. 2007;71(6):1934–43. https://doi.org/10.2193/2006-262.
Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79(17):5112–20. https://doi.org/10.1128/AEM.01043-13.
Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7. https://doi.org/10.1038/s41587-019-0209-9.
Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13(7):581–3. https://doi.org/10.1038/nmeth.3869.
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(16):5261–7. https://doi.org/10.1128/AEM.00062-07.
McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012;6(3):610–8. https://doi.org/10.1038/ismej.2011.139.
Werner JJ, Koren O, Hugenholtz P, DeSantis TZ, Walters WA, Caporaso JG, et al. Impact of training sets on classification of high-throughput bacterial 16s rRNA gene surveys. ISME J. 2012;6(1):94–103. https://doi.org/10.1038/ismej.2011.82.
Price MN, Dehal PS, Arkin AP. FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5(3):e9490. https://doi.org/10.1371/journal.pone.0009490.
Core Team R. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2018.
Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. J Appl Environ Microbiol. 2005;71(12):8228–35. https://doi.org/10.1128/AEM.71.12.8228-8235.2005.
Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, Minchin PR, O'Hara RB, Simpson GL, Solymos P et al: vegan: Community Ecology Package. 2017.
Mandal S, Van Treuren W, White RA, Eggesbo M, Knight R, Peddada SD. Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb Ecol Health Dis. 2015;26:27663.
We recognize the contributions of many biologists from state and federal agencies who helped capture and process mule deer. We also thank two anonymous reviewers for helpful comments that guided improvement of this manuscript.
Funding was provided by the Utah Division of Wildlife Resources and partners including Sportsmen for Fish and Wildlife and the Mule Deer Foundation. The funding body had no role in the design, analysis, interpretation, or writing of the data in the manuscript in any way.
Ethics approval and consent to participate
As stated in the methods, “Capture and handling procedures were reviewed, and evaluated, and approved by by the Brigham Young University’s (BYU) Animal Care and Use Committee (protocol 150110). All capture and handling procedures complied with the approved methods” Also, all animals were captured and handled by Utah Division of Wildlife Resources, the legal regulatory and management authority for mule deer in the State of Utah, which approved and granted permission for all activities.
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.
. Rarefied (13,000) OTU’s with counts per sample. Table S2. Primer Barcodes and Relevant Metadata for All Samples. Table S3. Weighted and unweighted Adonis tables for MaxFat, LoinThickness, Age, and RumpFat. Table S4. Phylum level correlations for taxa (December samples only) with March health outcomes in repeat deer. Table S5. Order level correlations for taxa (December samples only) with March health outcomes in repeat deer. Table S6. Family level correlations for taxa (December samples only) with March health outcomes in repeat deer. Table S7. Genus level correlations for taxa (December samples only) with March health outcomes in repeat deer. Table S8. ASV level correlations for taxa (December samples only) with March health outcomes in repeat deer. Table S9. Phylum level correlations for taxa (all samples) with the health measures of deer at time of sample collection. Table S10. Order level correlations for taxa (all samples) with the health measures of deer at time of sample collection. Table S11. Family level correlations for taxa (all samples) with the health measures of deer at time of sample collection. Table S12. Genus level correlations for taxa (all samples) with the health measures of deer at time of sample collection. Table S13. ASV level correlations for taxa (all samples) with the health measures of deer at time of sample collection.
. Alpha Rarefaction Curve. Rarefaction curve suggests optimal sequencing depth for rareified samples to be approximately 13,000 reads. 104 samples out of 108 are retained.
About this article
Cite this article
Eddington, H.S., Carroll, C., Larsen, R.T. et al. Spatiotemporal variation in the fecal microbiota of mule deer is associated with proximate and future measures of host health. BMC Vet Res 17, 258 (2021). https://doi.org/10.1186/s12917-021-02972-0