Discrete patterns of microbiome variability across timescales in a wild rodent population
BMC Microbiology volume 23, Article number: 87 (2023)
Mammalian gastrointestinal microbiomes are highly variable, both within individuals and across populations, with changes linked to time and ageing being widely reported. Discerning patterns of change in wild mammal populations can therefore prove challenging. We used high-throughput community sequencing methods to characterise the microbiome of wild field voles (Microtus agrestis) from faecal samples collected across 12 live-trapping field sessions, and then at cull. Changes in α- and β-diversity were modelled over three timescales. Short-term differences (following 1–2 days captivity) were analysed between capture and cull, to ascertain the degree to which the microbiome can change following a rapid change in environment. Medium-term changes were measured between successive trapping sessions (12–16 days apart), and long-term changes between the first and final capture of an individual (from 24 to 129 days). The short period between capture and cull was characterised by a marked loss of species richness, while over medium and long-term in the field, richness slightly increased. Changes across both short and long timescales indicated shifts from a Firmicutes-dominant to a Bacteroidetes-dominant microbiome. Dramatic changes following captivity indicate that changes in microbiome diversity can be rapid, following a change of environment (food sources, temperature, lighting etc.). Medium- and long-term patterns of change indicate an accrual of gut bacteria associated with ageing, with these new bacteria being predominately represented by Bacteroidetes. While the patterns of change observed are unlikely to be universal to wild mammal populations, the potential for analogous shifts across timescales should be considered whenever studying wild animal microbiomes. This is especially true if studies involve animal captivity, as there are potential ramifications both for animal health, and the validity of the data itself as a reflection of a ‘natural’ state of an animal.
Gastrointestinal microbiome composition is complex and understanding the causes and consequences of changes in microbiome structure can be challenging. The microbiome can be influenced by a range of environmental factors, including infection status, nutrition, and life history, all of which can be drivers of microbiome structure and diversity . Changes in the microbiome can in turn have impacts on a variety of host phenotypes  including infection responses, food metabolism, pathogenicity of the bacterial taxa themselves, and ultimately host fitness. While some microbial populations may be intrinsically more dynamic in abundance than others, some may also show specific shifts clearly associated with factors such as ageing [3,4,5], helminth infection [6, 7], and diet [8,9,10]. In contrast, a significant proportion of the bacterial community is comparatively stable, comprising a ‘core-microbiome’ of established taxa. The term ‘core microbiome’ is used to refer to communities that are stable either within, or between individuals (or both), and so high-resolution sampling both within and between individuals is required to expose any shifts in this core . Although the important role of archaea and fungal species in microbiome communities is becoming increasingly well understood, this study focusses solely on bacterial communities, and as such ‘microbiome’ will hereafter refer to bacterial taxa.
While temporal variation around this core has been extensively examined in human populations [12,13,14] and laboratory model species [15, 16], the prominence and nature of consistent patterns of change across varying timescales in wild mammal populations is currently underexplored and invites further investigation. These shifts can affect both α-diversity (the local diversity of species in a community, and β-diversity (the difference in taxa between communities.) Human populations are analogous to wild animal populations, being genetically diverse, and experiencing a diversity of lifestyles, diets and immune challenges. The fact that modern humans are typically far more long-lived means ageing effects may not be as prevalent in typical wild animal populations. Despite this, age-related shifts to Bacteroidetes-dominant microbiomes have been recorded in laboratory studies of mice ZHANG, and non-model species including marmosets , and pigs . Ageing effects may, therefore, scale somewhat to the ages of the host species, as is observed in metabolism , physical function , and mutation rates [21, 22].
Age-associated changes reported in human populations vary in nature, for instance depending on existing infections , but are typically characterised by increased proportions of Bacteroidetes species [24,25,26]. Such shifts may not be a linear, as the Bacteroidetes/Firmicutes ratio has been shown to initially decrease from infancy to adulthood, before increasing again in old age . These changes are associated with dysbiosis, characterised by increased abundance of pathogenic taxa, and have been implicated or associated with poorer health , typically coincident with increased levels of frailty and inflammation , and more specific pathologies, such as kidney disease, which is associated with reduced Firmicutes and increased Fusobacteria and Proteobacteria . α-diversity of gut microbiomes increases with age in humans , lab mice , and non-model species [18, 33]. Despite this, studies on some species have shown the opposite effect [17, 34]. These discrepancies may be species-specific, but may also be due to conflation of ‘chronological age’ (a simple measure of time), which is associated with increased richness and ‘biological age’ (maturation of host physiology and increasing physical frailty) which is associated with reductions in richness and associated pathology—while biological and chronological age are strongly associated, they will not always necessarily increase together in a linear fashion [3, 35].
Short-term changes, over the scale of days and hours, have been demonstrated in lab mice in response to stress , changes in diet [37,38,39], as a result of host diurnal rhythmicity , and in human patients following severe injury . Rapid changes in faecal community structure can also occur external to the host, depending on exposure to varying environmental conditions . Microbiome studies of wild populations commonly involve some element of captivity  (Morgan and Tromborg 2007)to allow for ease of monitoring or sample/tissue collection, and so it is important to know whether such changes are occurring following capture. Marked differences between the microbiome of laboratory and wild populations have been recorded in a range of species, with domestication or captive breeding commonly associated with reduction in microbial diversity or loss of taxa, relative to wild populations, and changes in translocated individuals occurring relatively rapidly [44,45,46,47,48]. Again, these shifts are not universal, with some species showing increased microbiome richness under captivity , and more robust microbiomes, showing reduced captivity-associated differences over time [50, 51].
Studies exploring these patterns in detail across timescales in wild populations are lacking, or are often limited. For example, a longitudinal analysis of mouse lemur microbiomes (Microcebus rufus) showed reduced α-diversity in older individuals, but was restricted to 15 individuals . Studies of populations thatare more reliably trappable at high numbers allow us to build upon work on human, lab and wild populations, and more reliably understand the prevalence and significance of characteristic patterns of change in microbiome communities across timescales.
Microbiome studies are typically cross-sectional, allowing for invasive sampling of the gut post-cull, and thus providing what is hopefully the most accurate snapshot of the live animal’s GI microbiome community. While longitudinal microbiome studies must typically rely on faecal samples rather than direct sampling of the gut, diversity metrics have been shown to be highly correlated between faeces and caecum samples, making faeces a suitable representation of microbiome communities in the live animal . Dynamic changes in the microbiome, and their environmental and host-intrinsic causes, are becoming increasingly well characterised, and so the importance of emphasising longitudinal experimental design in wild animal microbiome studies is becoming more apparent .
We investigated the gastrointestinal microbiome of a wild population of the field vole, Microtus agrestis. Using both longitudinal faecal samples from mark-recapture trapping, and faecal samples taken after capture and dissection, we examined the level of between- and within-individual variation associated with different bacterial phyla. We explored how levels of species richness, and the balance between representation of Bacteroidetes and Firmicutes as dominant phyl constitute patterns of change in microbiome structure over different timescales, as well as what ecological factors such as host condition and infection status might be associated with this variation. Patterns of changes over different timescales in wild animals are likely to be complex and somewhat population-specific, but by characterising those patterns in a well-studied population, we aim to provide an example of how such changes might manifest, and how they might be studied.
Trapping & dissection
Field methods are based on previous studies of this population [55, 56], and are also described elsewhere . Voles were live-trapped using a grid composed of 197 Ugglan traps, over approximately one hectare of a clear-felled area in Kielder Forest, Northumbria, UK. Trapping was conducted over twelve 3-day sessions between March and August 2017, with traps checked twice each day, in the morning and evening. Traps were baited with mixed seed and chopped carrot. Newly-trapped animals were tagged with Passive Integrated Transponders (PIT Tags) on first capture to allow for re-identification. Upon first capture, animals were visually inspected, and sex and reproductive status were recorded. At each capture, snout-vent body length was measured, and total body mass was recorded. These values were used to calculate the scaled mass index (SMI), a measure of body condition .
For our longitudinal sampling, 428 faecal samples were collected from 206 voles. The mean number of samples per individual was 2.28, and 3.26 when excluding individuals with only one capture. In most cases these animals were then released, but at the end of each trapping session a small number of the animals were sent to the University of Nottingham for dissection, to perform gastrointestinal helminth surveys, collect caecum samples, and extract eye lenses as a proxy measure of animal age . A total of 60 voles were sampled in this way, forming the cross-sectional dataset. They were housed for either one or two nights prior to dissection, fed bird seed mix and chopped carrot, and provided with water ad libitum. Animals were killed by increasing CO2 concentration in a sealed chamber, with death confirmed by exsanguination. Procedures were performed with approval from the University of Liverpool Animal Welfare Committee, under a UK Home Office license (PPL 40/3235 to MB. Field-to-lab workflow is shown in Fig S1.)
Morphometric measurements taken at cull include mass, snout-vent length and tail length. Eyes were removed and stored in formalin. Later, eye lenses were removed, desiccated at 60 °C for 48 h, and weighed with an electronic balance for use as an age proxy [59, 60]. Gastrointestinal tracts were removed and stored in 80% ethanol. These animals had further faecal samples taken at cull, and for 44 of those animals, caecum tissue was also taken. All faeces and caecum samples were stored at -80 °C.
Gastrointesintal helminth surveys
Gastrointestinal tracts were dissected under dissection microscopes, and any helminth macroparasites identified morphologically and counted. Descriptions, keys or other literature on which the identifications were based are given alongside the respective helminth taxa below. Of the animals included in this study, two types of macroparasite were commonly observed – the pinworm Syphacia nigeriana (64.3% prevalence) (Chromadorea: Oxyuridae) [61, 62] & tapeworms (50% prevalence) (Cestoda: Hymenolepidae and Anoplocephalidae) . Other species recorded in this population include Heligmosomoides laevis (Chromadorea: Heligmosomidae) and Trichuris arvicolae (Enoplea: Trichuridae) , but were not present in the animals recorded in this study.
DNA extraction & microbiome sequencing
DNA was extracted from faeces and caecum tissue using the DNeasy Powersoil extraction kit (Qiagen Cat. 47,016) and sent for 16S community sequencing at University of Liverpool Centre for Genomic Research. Alongside these samples, positive and negative controls were included, provided in-house at the Centre for Genomic Research, Liverpool. Primers described by Caporaso et al., 2011  were used to amplify and barcode the V4 region of 16 s (detail of primers for every stage are provided in Table S1.) A total of 658 samples were submitted and 5 μl of each DNA sample at 1 ng/μl was entered into the first-round PCR with total reaction volume of 20ul, and the following conditions: 98 °C for 2 min, 20 s at 95 °C, 15 s at 65 °C, 30 s at 70 °C for 10 cycles then a 5 min extension at 72 °C. Samples were then purified with AMPure SPRI beads in a 1:1 volume ratio (Beckman Coulter, Indiana, USA), and a secondary, nested PCR was then performed to incorporate i5 & i7 Illumina adapter sequences, using the same conditions for a further 15 cycles. Samples were again purified with SPRI beads in a 1:1 volume ratio and quantified by Qubit dsDNA HS Assay (Thermo Fisher Scientific, Massachusetts, US) using the Agilent Fragment Analyzer (Agilent Technologies, Santa Clara, US). Samples which failed to amplify were not sequenced.
Final amplified libraries and controls were pooled in equimolar amounts into 8 pools according to the Qubit data, followed by size selection on a 1.5% Pippin prep gel (Sage Science Inc., Massachusetts) using a size range of 300-600 bp. Quantity of the size selected pools of amplicon libraries was completed using a Qubit dsDNA HS Assay Kit (Thermo Fisher Scientific, Massachusetts, US), while the quality and average size was assessed using the High Sensitivity DNA Kit on the Agilent Bioanalyzer (Agilent Technologies, California, US).
Subsequently, a quantitative real-time PCR (qPCR) assay, designed to specifically detect adapter sequences flanking the Illumina libraries, was performed using an Illumina KAPA Library Quantification Kit (Kapa Biosystems, Wilmington, USA). A 20 μl PCR reaction (performed in triplicate for each pooled library) was prepared on ice with 12 μl SYBR Green Master Mix and 4 μl diluted pooled DNA (1:1000 to 1:100,000 depending on the initial concentration determined by the Qubit dsDNA HS Assay Kit). PCR thermal cycling conditions consisted of initial denaturation at 95 °C for 5 min, 35 cycles of 95 °C for 30 s (denaturation) and 60 °C for 45 s (annealing and extension), melt curve analysis to 95 °C (continuous) and cooling at 37 °C (LightCycler LC48011, Roche Diagnostics Ltd, Burgess Hill, UK). Template DNA was denatured for 5 min at room temperature using freshly diluted 0.1 N-sodium hydroxide (NaOH) and the reaction was subsequently terminated by the addition of hybridization buffer. Following calculation of the molarity using qPCR data, template DNA was diluted to a loading concentration of 8 pM using the hybridization buffer. The amplicon libraries were sequenced on an Illumina HiSeq 2500 platform (Illumina Inc., California, US) with version 2 chemistry using sequencing by synthesis (SBS) technology to generate 2 × 300 bp paired-end reads. Fragmented PhiX phage genome was added to the sequence library to increase the sequence complexity.
Base-calling and de-multiplexing of indexed reads was performed by CASAVA version 1.8.2 (Illumina) to produce 658 samples across the two runs, in FASTQ format. The raw FASTQ files were trimmed to remove Illumina adapter sequences using Cutadapt version 1.2.1 . Any reads which matched the adapter sequence over at least 3 bp were trimmed off. The reads were further trimmed to remove low quality bases, using Sickle version 1.200  with a minimum window quality score of 20. After trimming, reads shorter than 20 bp were removed. If both reads from a pair passed this filter, each was included in the R1 (forward reads) or R2 (reverse reads) file. If only one of a read pair passed this filter, it was included in the R0 (unpaired reads) file. To improve base quality in both read pairs, sequencing errors were corrected in both forward and reverse reads using the error-correct module within SPAdes assembler, version 3.1.0  using options '–careful' and '–only-error-correction'. The average number of paired-end reads per sample was 512,850 (SD = 161,805, IQR = 145,881).
Read pairs were merged to produce a single sequence for each pair that would entirely span the amplicon using PEAR, version 0.9.10 . Additionally, sequences with uncalled bases (Ns) were removed. To remove sequences originating from potential PCR primer dimers or from any spurious amplification events, a size selection was applied to select sequences between 200 and 600 bp. To remove PhiX sequences associated with indices, each sample was compared with the complete PhiX sequence (GenBank gi9626372) using BLASTN . Sequences matching PhiX (E-value < 10–5) were filtered out of the dataset. An average of 99.67% of reads per-sample, were successfully aligned, which was reduced to 98.57% post PhiX-filtering.
Sequences passing the filters for each sample were concatenated into a single file, which was used subject to a custom analysis pipeline based on QIIME 1.9.1 . The RDP classifier was used against the GreenGenes database (version 13.8).throughout the analysis. To identify sequence variability in each sample, amplicon sequences and assigned to clusters according to sequence similarity, using SWARM version 2.2.1, . To calculate the abundance of each cluster, sequences were aligned on the identified centroid clusters sequences, using a minimum similarity threshold of 97% for the entire length of the sequence.
Taxonomic assignment of each cluster (now referred to as operational taxonomic unit, OTU) was carried out using the QIIME script ‘assign_taxonomy.py’, using the RDP classifier  to match the centroid sequence of each cluster obtained by swarm, to a sequence from the database. The abundance table was post-processed to remove any OTU below 0.005% of the total sequence count of sequences .
All statistical analyses were carried out in R 3.6.2 , and read counts were centered log-ratio (CLR) transformed using the ‘SleuthALR’ package . The package ‘phyloseq’ was used to calculate measures of α-diversity. Three metrics were chosen to assess different aspects of α-diversity, with Shannon index values emphasising taxa evenness, phylogenetic diversity (hereafter, ‘PD’) being correlated with evenness  and Chao1 being suitable for datasets skewed to low-abundance taxa  Bray–Curtis and weighted UniFrac (wUniFrac) distances were calculated with phyloseqand used in Non-Metric Multidimensional Scaling (NMDS) to provide individual (site) scores (wUniFrac K = 5, stress = 0.0012, Bray–Curtis K = 3, stress = 0.164. In addition, robust principal component analysis (RPCA) was performed using the ‘rospca’ package . UniFrac distance incorporates OTU relatedness data from a provided phylogenetic tree, and wUniFrac adjusts this distance to reduce the influence of rare OTUs and alleviate any oversized influence of rare taxa by taking abundances into account. Bray–Curtis is an abundance-based metric, whereby distance values give a measure of between-sample dissimilarity, but which are sensitive to presence of rare taxa. RPCA is also abundance-based, but can better deal with sparse, highly-dimensional datasets. Three OTUs of a total 1321 were excluded due to severe overrepresentation in specific individuals, and causing issues with MDS and RPCA ordination. For both α- and β-diversity measures, multiple indices were calculated to provide a more comprehensive picture of microbial diversity, and capture any variation more weighted to specific diversity metrics.
Using CLR abundance data, the coefficient of variance (CoV) was calculated for each OTU, based on within-individual variation from longitudinal faecal samples, and cross-sectionally between individuals at cull. These values were normalised using Box-Cox transformation, a power-transformation for positive non-normal data (Longitudinal λ = 0.51, Cross-Sectional Lambda λ = -0.42) in R  using the ‘car’ and ‘MASS’ packages [80, 81]. Transformed data was then used to measure correlation between variation in OTUs from cross-sectional faecal samples and longitudinal faecal samples, to ascertain bywhether OTUs which were variable between individuals were the same as those which were variable within the same individuals over time. Average OTU variance per phylum was compared by linear model, with post-hoc Tukey analysis to determine differences in variability between major phyla (phyla with fewer than 10 OTUs present, and unclassified OTUs, were excluded).
Both cross-sectional faecal and longitudinal faecal community datasets were subject to variance partitioning using the ‘VariancePartition’ R package to ascertain the ecological factors associated with variation in CLR abundance . These factors include age category (designated as ‘mature’ or ‘juvenile’ from physical inspection of body size), the within-year Julian date (i.e. day number 1–365), sex and scaled mass index (SMI)  as a measure of condition. Julian date was chosen to account for seasonal affects as it offers a more precise continuous measure of time than simply including the trapping session, and as all samples in this analysis were collected in 2017, there was no need to account for inter-annual patterns of variation. For cross-sectional cull samples, the days kept in captivity and prevalence of gastrointestinal helminth infections were included, and for longitudinal samples, the individual ID.
Subsequent analyses were carried out to determine changes in α- and β-diversity over three time-scales; short-term changes from field to lab following capture, medium-term changes between sequential trapping sessions, and long-term changes between the first and last longitudinal samples of individual voles. All changes were assessed in R through Gaussian linear mixed models (LMER) using the lme4 package  and tested for significance with ‘lmerTest’ package , with individual ID as a random effect, and incorporating the Julian date from the start of the year as a fixed effect to account for any seasonal changes. These models were performed on α-diversity metrics, including total unique OTU count as a measure of richness, and Bray–Curtis, wUniFracand RPCA scores to measure changes in β-diversity.
Short-term models compared samples taken at cull with the most recent longitudinal sample from the same animal (1–2 days prior). Medium-term models incorporated all consecutive longitudinal faecal samples with no missed trapping sessions between sessions (12–16 days between each sample), and assessed the size of changes observed from the first capture in the sequence as explained by the number of sessions passed. Long-term models compared paired longitudinal samples from the first and last captures, excluding all cases where that time difference was only one trapping session, with time between first and last samples also included as a factor. Trapping date was included as a factor in all models, in an effort to show changes in diversity, regardless of the time of year in which the trapping sequence occurred. However, seasonal effects can still not be completely discounted, as all traps were recorded within 2017, and as such, a general seasonal change in microbiome diversity affecting all animals from early Spring to Autumn would be difficult to distinguish from any effects due purely to ageing.
In order to assess associations between diversity metrics and other ecological and environmental factors, linear models were performed on measures of cull faecal α-diversity and Bray–Curtis and wUniFrac site scores to measure changes in β-diversity, which incorporated Julian date, scaled mass index (SMI) as a measure of condition, sex and prevalence of two gastrointestinal helminths which are common in the population – the pinworm Syphacia obvelata and tapeworms (Class: Cestoda). In order to incorporate infection data from gut dissections, these analyses were performed on a reduced cross-sectional dataset, using only cull faecal data.
Taxonomic structure & variance
Across all sample samples, Firmicutes and Bacteroidetes constituted the majority of OTUs sequenced, (summarised in Figs. 4 and 6). Bray–Curtis MDS1, wUniFrac MDS2 and RPC1 allprovided an axis which clearly distinguished relative prevalence of these two phyla, with higher scores of both Bray–Curtis MDS1 and wUniFrac MDS2 and RPC1 representing a significant shift towards a more Bacteroidetes-dominant microbiome (Fig. 1). (OTUs most strongly represented in loadings are in Tables S3, S4 & S5).
Within-individual OTU variance, measured from longitudinal faecal samples, was consistent with between-individual variance of OTUs from cull faecal samples, indicating that bacterial taxa which are more variable between individuals are also more variable within an individual over time (Linear model, p = < 2.6 × 10–16, Adjusted R2 = 0.78; see Fig S2). OTUs of phylum Firmicutes were significantly more variable than those of Bacteroidetes (GLM with Tukey post-hoc: within-individual p < 0.0001, faeces p < 0.0001) (Cross-sectional differences shown in Fig. 2).
Variance partitioning was performed to show the percentage of variation in theCLR abundance of each OTU that is explained by each variable used in statistical models (Fig. 3). For the longitudinal data, these variables include individual ID, the Julian date, body condition (SMI), while for the cross-sectional data they include Julian date, the days spent in captivity prior to cull, SMI, eye lens mass as a proxy for age, sex, and prevalence of pinworm (Syphacia obvelata) and tapeworm infections. In the longitudinal data individual ID typically explained 8–25% of variation in most OTUs, indicating a substantial level of within-individual consistency in microbiome structure. A small number of OTUs were strongly associated with sex – of 46 OTUs which had over 50% of their variation explained by sex in the longitudinal dataset, over half  were of order Clostridiales. Variance partitioning of cross-sectional data found that ecological factors explained a small amount of variation for most OTUs, with 3 OTUs having over 50% of their variation explained (2 of order Bacteroidales, 1 of phylum Tenericutes). There is, however, still a great deal of unexplained residual variation across both datasets for the majority of OTUs.
Short-term changes (Field to Lab)
(Summary of model outputs is provided in Table S2). Shifts in microbial α-diversity were recorded between faecal samples from the field and matched samples from the lab 1–2 days later. Faecal samples taken at cull showed significant decreases in α-diversity compared to paired longitudinal samples (LMER Chao1 p = 1.60 × 10–8, PD p = 2.35 × 10–12, Shannon p = 6.29 × 10–15). Significant differences inBray-Curtis (GLMER p < 2 × 10–16),wUniFrac (GLMER p = 2.33 × 10–7) and RPCA (GLMER p = 6.12 × 10–8) indicate shifts in β-diversity corresponding to a more uniform, Bacteroidetes-dominant microbiome upon arrival in captivity. Observed short-term changes are summarised in Fig. 4.
Medium-term changes (between successive trapping sessions)
Microbiomes showed increases in measures of α-diversity between successive trapping sessions, though only in Chao1(LMER, Chao1 p = 0.0177, PD = 0.279, Shannon p = 0.97). The inclusion of individual ID as a random effect means this is independent of any population level effects such as survival bias, and less likely to be due to seasonal effects, as any observed shifts are relative to each animal’s starting state, regardless of how it is being influenced by environmental variables at that time. While all three measures of β-diversity did show slight shifts towards Bacteroidetes dominance between sessions, these changes were not statistically significant. Observed medium-term changes are summarised in Fig. 5.
α-diversity was significantly higher at the final longitudinal capture relative to the first, indicating an accrual of bacterial species over long periods (LMER, p = 0.0035, PD p = 7.97 × 10–4, Shannon p = 0.0248). A significant difference is observed in Bray–Curtis distance between first and last captures (GLMER, p = 0.043), indicating a shift towards a more Bacteroidetes-dominant microbiome. No significant difference was observed in wUniFrac RPCA distances (GLMER: wUniFrac p = 0.058, RPCA p = 0.149). Observed long-term changes are summarised in Fig. 6.
Cull faeces α-diversity was positively associated with body condition, measured as SMI, though not in Shannon diversity (GLM, Chao1 p = 0.00323, PD p = 0.00861, Shannon p = 0.51884), indicating that richer microbiome communities are associated with better host health. SMI was also significantly associated with the Bray–Curtis distance, with Bacteroidetes dominance associated with lower condition (GLM, p = 0.0126), but not wUniFrac (GLM p = 0.162) or RPCA (GLM p = 0.798) distance. Tapeworm infection was positively associated with Bacteroidetes dominance as explained by Bray–Curtis (GLM, p = 0.0313), and RPCA (GLM, p = 0.00931) distances, but not wUniFrac (GLM, p = 0.162) (These associations are summarised in Fig S4).
Differences between bacterial taxa in how they are affected by temporal patterns are a crucial element of wild animal ecology, as taxon-specific shifts in microbiome structure have been implicated in many aspects of host health . Here we have shown clear and defined patterns of change in gastrointestinal microbial diversity across different time scales, including long-term ageing-associated increases in α-diversity and shifts in β-diversity towards Bacteroidetes-dominance, and rapid reductions in α-diversity associated with capture and captivity. Firmicutes taxa showed significantly higher levels of variation than those of Bacteroidetes, and so while increases in α-diversity suggest an accrual of bacteria over time, it is likely that shifts in β-diversity are at least in part explained by loss of Firmicutes species. These phylum-level differences in stability have been reported in human populations, presumably due to differences in susceptibility to changes in environmental conditions . These observed changes increase our understanding of natural variation in microbiome structures in wild populations, show clear parallels to analogous health-related shifts observed in ageing humans, and provide important context to consider when performing any analysis of wild animals involving captivity or rehousing. Increases in α-diversity with age, over both medium-term and long-term timescales show a pattern of constant species accrual throughout an animal’s lifetime, while positive associations with species richness and body condition suggest an important role in this process for maintaining host health.
The underlying taxonomic composition was broadly what would be expected of mammalian gastrointestinal microbiomes [87, 88], and has been recorded in other wild rodent populations  with Firmicutes and Bacteroidetes being the dominant phyla. While Bacteroidetes was generally more prevalent than Firmicutes in this population, in domesticated populations of Microtus ochrogaster Firmicutes was the more dominant phylum . Most variation in OTU abundance was either associated with the specific animal, or was unexplained by any of the accompanying data. This unexplained variation may arise from purely stochastic processes, and/or there may be factors not captured in this study which play a role in shaping microbiome communities, such as dietary variation . While field voles are predominantly herbivorous, variation within the plant species consumed, or occasional consumption of other foods like insect larvae may influence microbiome composition [91,92,93], and so alteration in diet after capture may have an impact on microbiome communities. Recent research in wild rodent population has also highlighted a potentially significant role for social network structure in determining population-level microbiome diversity . While host-genetic influences not accounted for in this study, the impact of the genome on microbiome structure has shown to be secondary to environmental and temporal factors in both humans  and rodents .
While we have characterised distinct patterns of variability across different timescales, it should be noted that the observed short-term changes are likely to have limited relevance to ecological processes that would occur without human intervention, as the conditions and circumstances of animal capture and rehousing are not directly analogous to any naturally occurring process. Despite this, they do still illustrate that dramatic shifts in microbiome structure can occur in short timescales dependent on surrounding conditions, and that any microbiome study involving captivity and capture may be significantly confounded, even when efforts are made to minimise interference. With these caveats, we have been able to identify specific and directed changes in the microbiome across multiple time scales.
Short-term changes and effects of capture
α-diversity was reduced between paired final live-capture and cull samples, across multiple richness indices, indicating rapid diversity loss associated with capture and 1–2 days of captivity. Reduced α-diversity associated with captivity has also been reported in a number of mammals , though this effect is not universal, and in some instance the opposite may be observed [97, 98]. The changes observed in this study both confirm that captivity-associated changes in microbiome are common across species and populations, and show how rapidly these changes can arise. There are many potential direct causes of these changes; changes in diet have been suggested as a key contributor [38, 96], but other factors such as stress [99, 100], disruption to diurnal sleep cycles , or combinations of these factors  may also play a role. These changes confirm that the microbiome can be highly volatile, and captivity-associated changes are essential to consider, both in terms of health of captive animals, and for studies of wild populations. This is particularly relevant for studies involving live trapping and rehousing of animals prior to sample collection, as changes from ‘natural’ conditions could significantly impact the microbiome of the animal in question within a short window of time. It is unlikely that these changes will be universal in nature or magnitude across study species, as captivity has been shown to have differing impacts on the gut microbiome of different host species, even within the same genus, potentially having a reduced impact on generalists compared to specialists [103, 104].
A significant taxonomic shift in microbiome communities was also observed following capture. While such shifts are in the same direction as what would be expected from ageing (increased Bacteroidetes-dominance), the magnitude and speed of the change suggests that capture and captivity itself is impacting the microbiome. Alongside the decreases in species, this indicates a loss of Firmicutes, resulting in increased prominence of Bacteroidetes. Shifts towards Bacteroidetes-dominant microbiome following captivity have been observed in deer mice (Peromyscus maniculatus) examined pre- and post-captivity , as well as in comparisons between captive mammals and counterpart wild populations [96, 97, 105, 106].
Medium- & long-term changes in microbiome structure
Both the medium-term (between successive trapping sessions) and long-term (between first and last live-traps) models showed that microbiome α-diversity increases with age, regardless of the time of year, with no obvious drop-off observed in the oldest voles. They also showed that β-diversity shifted toward a more Bacteroidetes-dominant microbiome throughout an animal’s life. Due to the diversity metrics used in this study, some of these patterns may be driven in part by rare taxa having an outsized impact. One Bacteroidetes OTU in particular has a very strong representation in Bray–Curtis MDS1, but was not classified below phylum level (Table S 3). OTUs which strongly contributed to Firmicutes dominance in Bray–Curtis and RPCA distances were of the genus Ruminococcus, a taxon which is important in cellulolytic metabolism, and so loss of these bacteria over time may have an impact on host digestion and health  Age-associated increases in α-diversity are recorded in humans, from infancy to adulthood . Age-related shifts in taxonomic structure of the microbiome have been recorded in laboratory mice  and captive mammals [17, 109], with multiple human studies showing specific age-related shifts to Bacteroidetes [26, 27]. The combination of changes in both α- and β-diversity observed in this study suggest that as the voles age, they accrue new bacterial OTUs, primarily of Phylum Bacteroidetes.
Other correlates of microbiome diversity and fitness implications
While most ecological factors were not associated with microbiome diversity, body condition, measured as SMI, was found to be positively associated with both α-diversity, and a more Firmicutes-dominant microbiome. The association with α-diversity may suggest that while animals in different life stages and under different constraints may harbour qualitatively different microbiome communities, it is the richness of those communities which is most significant for host condition. The positive impacts of age and infection on α-diversity highlight how the state of the gut can be associated with multiple ecological factors which may be having indirect effects on host condition via the microbiome. On the other hand, rapid reductions in α-diversity indices associated with capture and captivity are essential to consider as context for microbiome analyses of wild and captive populations, and when considering questions of captive animal health and welfare. The negative association between Bacteroidetes and scaled mass index mirrors what has been observed in human populations, although in humans, high-BMI, particularly in obese individual, is associated with Firmicutes-dominance, and thus higher mass indices broadly corresponded to pathology, rather than healthy condition [110, 111].
The Bacteroidetes/Firmicutes ratio is commonly implicated as a key factor affecting gastrointestinal, and general, health, and so this change may be relevant to health and survival over time in wild populations. Many studies have looked at the Bacteroidetes/Firmicutes ratio as an important correlate of gut health, with alterations in this ratio linked to obesity and other pathologies [112, 113]. How these associations relate to ageing-related shifts in microbiome, in which both richness and Bacteroidetes-dominance increase is unclear, and this element of the analysis is somewhat limited by the reduced cross-sectional dataset. Further investigation could elucidate whether microbiome-associated variation in condition is primarily the result of inter-individual differences, or could be relevant to within-individual changes throughout an animal’s life.
The gastrointestinal microbiome is complex and dynamic, particularly in wild, heterogeneous populations, and the factors underlying temporal changes in its composition are often difficult to determine. Using both longitudinal and cross-sectional data from a well-characterised wild rodent population, we have established robust temporal patterns of changes in microbiome structure, both in short periods following change in environmental conditions, and over the course of an animal’s life in the wild.
While microbiome structure was found to be highly individual, broad patterns of change over different timescales can be observed in both α- and β-diversity. Ageing is associated with slight shifts towards Bacteroidetes-dominance and increase in species richness, while captivity was associated with marked and larger short-term increases in Bacteroidetes-dominance and drops in species richness. Firmicutes OTUs, being significantly more variable than those of Bacteroidetes, are likely responsible for many of the temporal patterns of change in microbiome structure which were observed.
These results provide a useful framework with which to understand time-sensitive, consistent changes in microbiome structure within wild animal populations. The specific characteristics of changes observed in this population invite direct comparisons with analogous findings from laboratory rodent models and human populations, and suggest a potentially significant role for temporal microbiome dynamics on host health and fitness.
Availability of data and materials
Full read count data metadata, taxonomy data and phylogenetic tree are publicly available on the Dryad repository: https://doi.org/10.5061/dryad.08kprr559.
Rothschild D, Weissbrod O, Barkan E, Kurilshikov A, Korem T, Zeevi D, et al. Environment dominates over host genetics in shaping human gut microbiota. Nature 2018 555:7695. 2018 28;555(7695):210–5.
Suzuki TA. Links between natural variation in the microbiome and host fitness in wild mammals. Integr Comp Biol. 2017;57(4):756–69.
Kim S, Jazwinski SM. The Gut Microbiota and Healthy Aging: A Mini-Review. Gerontology. 2018;64(6):513–20.
García-Peña C, Álvarez-Cisneros T, Quiroz-Baez R, Friedland RP. Microbiota and Aging. A Review and Commentary. Arch Med Res. 2017 48(8):681–9.
Bana B, Cabreiro F. The Microbiome and Aging. 101146/annurev-genet-112618-043650. 2019 53:239–61.
Zaiss MM, Harris NL. Interactions between the intestinal microbiome and helminth parasites. Parasite Immunol. 2016;38(1):5–11.
Rapin A, Harris NL. Helminth-Bacterial Interactions: Cause and Consequence. Trends Immunol. 2018;39(9):724–33.
Voreades N, Kozil A, Weir TL. Diet and the development of the human intestinal microbiome. Front Microbiol. 2014;5(SEP):494.
Singh RK, Chang HW, Yan D, Lee KM, Ucmak D, Wong K, et al. Influence of diet on the gut microbiome and implications for human health. J Transl Med. 2017;15:73.
David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2014;505(7484):559.
Risely A. Applying the core microbiome to understand host-microbe systems. J Anim Ecol. 2020;89:1549–58.
David LA, Materna AC, Friedman J, Campos-Baptista MI, Blackburn MC, Perrotta A, et al. Host lifestyle affects human microbiota on daily timescales. Genome Biol. 2014;15(R89):1–15.
Caporaso JG, Lauber CL, Costello EK, Berg-Lyons D, Gonzalez A, Stombaugh J, et al. Moving pictures of the human microbiome. Genome Biol. 2011;12(5):1–8.
Costello EK, Lauber CL, Hamady M, Fierer N, Gordon JI, Knight R. Bacterial community variation in human body habitats across space and time. Science. 2009;326(5960):1694–7.
Wang J, Lang T, Shen J, Dai J, Tian L, Wang X. Core gut bacteria analysis of healthy mice. Front Microbiol. 2019;10(APR):887.
Schloss PD, Schubert AM, Zackular JP, Iverson KD, Young VB, Petrosino JF. Stabilization of the murine gut microbiome following weaning. Gut Microbes. 2012;3(4):383–93.
Reveles KR, Patel S, Forney L, Ross CN. Age-related changes in the marmoset gut microbiome. Am J Primatol. 2019;81:e22960.
Lim MY, Song EJ, Kang KS, Do Nam Y. Age-related compositional and functional changes in micro-pig gut microbiome. Geroscience. 2019;41(6):935–44.
Ma S, Yim SH, Lee SG, Kim EB, Lee SR, Chang KT, et al. Organization of the Mammalian Metabolome according to Organ Function, Lineage Specialization, and Longevity. Cell Metab. 2015;22(2).
Justice JN, Cesari M, Seals DR, Shively CA, Carter CS. Comparative Approaches to Understanding the Relation Between Aging and Physical Function. J Gerontol A Biol Sci Med Sci. 2016 Oct 1 [cited 2022 Nov 28];71(10):1243. Available from: /pmc/articles/PMC5018556/
Cagan A, Baez-Ortega A, Brzozowska N, Abascal F, CoorensSanders THHMA, et al. Somatic mutation rates scale with lifespan across mammals. Nature. 2022;604(7906):517–24.
MacRae SL, Croken MMK, Calder RB, Aliper A, Milholland B, White RR, et al. DNA repair in species with extreme lifespan differences. Aging. 2015;7(12).
Liu J, Johnson R, Dillon S, Kroehl M, Frank DN, Tuncil YE, et al. Among older adults, age-related changes in the stool microbiome differ by HIV-1 serostatus. EBioMedicine. 2019;40:583–94.
Claesson MJ, Cusack S, O’Sullivan O, Greene-Diniz R, De Weerd H, Flannery E, et al. Composition, variability, and temporal stability of the intestinal microbiota of the elderly. Proc Natl Acad Sci U S A. 2011;108(Supplement 1):4586–91.
La-ongkham O, Nakphaichit M, Nakayama J, Keawsompong S, Nitisinprasert S. Age-related changes in the gut microbiota and the core gut microbiome of healthy Thai humans. 3 Biotech. 2020;10(276):276.
Odamaki T, Kato K, Sugahara H, Hashikura N, Takahashi S, Xiao JZ, et al. Age-related changes in gut microbiota composition from newborn to centenarian: A cross-sectional study. BMC Microbiol. 2016;16(90):1–12.
Mariat D, Firmesse O, Levenez F, Guimarǎes VD, Sokol H, Doré J, et al. The firmicutes/bacteroidetes ratio of the human microbiota changes with age. BMC Microbiol. 2009;9(123):1–6.
O’Toole PW, Jeffery IB. Gut microbiota and aging. Science. 2015;350(6265):1214–5.
Claesson MJ, Jeffery IB, Conde S, Power SE, O’connor EM, Cusack S, et al. Gut microbiota composition correlates with diet and health in the elderly. 2012;
Zhang J, Luo D, Lin Z, Zhou W, Rao J, Li Y, et al. Dysbiosis of gut microbiota in adult idiopathic membranous nephropathy with nephrotic syndrome. Microb Pathog. 2020;147:104359.
De La Cuesta-Zuluaga J, Kelley ST, Chen Y, Escobar JS, Mueller NT, Ley RE, et al. Age- and Sex-Dependent Patterns of Gut Microbial Diversity in Human Adults. mSystems. 2019;4(4):e00261-e319.
Scott KA, Ida M, Peterson VL, Prenderville JA, Moloney GM, Izumo T, et al. Revisiting Metchnikoff: Age-related alterations in microbiota-gut-brain axis in the mouse. Brain Behav Immun. 2017;65:20–32.
Ke S, Fang S, He M, Huang X, Yang H, Yang B, et al. Age-based dynamic changes of phylogenetic composition and interaction networks of health pig gut microbiome feeding in a uniformed condition. BMC Vet Res. 2019;15(1):1–13.
Smith P, Willemsen D, Popkes M, Metge F, Gandiwa E, Reichard M, et al. Regulation of life span by the gut microbiota in the short-lived african turquoise killifish. Elife. 2017;22:6.
Jackson MA, Jeffery IB, Beaumont M, Bell JT, Clark AG, Ley RE, et al. Signatures of early frailty in the gut microbiota. Genome Med. 2016;8(8):1–11.
Bastiaanssen TFS, Gururajan A, van de Wouw M, Moloney GM, Ritz NL, Long-Smith CM, et al. Volatility as a Concept to Understand the Impact of Stress on the Microbiome. Psychoneuroendocrinology. 2021;124:105047.
Carmody RN, Gerber GK, Luevano JM, Gatti DM, Somes L, Svenson KL, et al. Diet dominates host genotype in shaping the murine gut microbiota. Cell Host Microbe. 2015;17(1):72–84.
Zarrinpar A, Chaix A, Yooseph S, Panda S. Diet and feeding pattern affect the diurnal dynamics of the gut microbiome. Cell Metab. 2014;20(6):1006–17.
Shang Y, Khafipour E, Derakhshani H, Sarna LK, Woo CW, Siow YL, et al. Short Term High Fat Diet Induces Obesity-Enhancing Changes in Mouse Gut Microbiota That are Partially Reversed by Cessation of the High Fat Diet. Lipids. 2017;52(6):499–511.
Liang X, Bushman FD, FitzGerald GA. Rhythmicity of the intestinal microbiota is regulated by gender and the host circadian clock. Proc Natl Acad Sci U S A. 2015;112(33):10479–84.
Howard BM, Kornblith LZ, Christie SA, Conroy AS, Nelson MF, Campion EM, et al. Characterizing the gut microbiome in trauma: significant changes in microbial diversity occur early after severe injury. Trauma Surg Acute Care Open. 2017;2(e000108):1–6.
Wong K, Shaw TI, Oladeinde A, Glenn TC, Oakley B, Molina M. Rapid Microbiome Changes in Freshly Deposited Cow Feces under Field Conditions. Front Microbiol. 2016;7(500):1–12.
Morgan KN, Tromborg CT. Sources of stress in captivity. Appl Anim Behav Sci. 2007;102(3–4):262–302.
Uenishi G, Fujita S, Ohashi G, Kato A, Yamauchi S, Matsuzawa T, et al. Molecular analyses of the intestinal microbiota of chimpanzees in the wild and in captivity. Am J Primatol. 2007;69(4):367–76.
Wienemann T, Schmitt-Wagner D, Meuser K, Segelbacher G, Schink B, Brune A, et al. The bacterial microbiota in the ceca of Capercaillie (Tetrao urogallus) differs between wild and captive birds. Syst Appl Microbiol. 2011;34(7):542–51.
Xenoulis PG, Gray PL, Brightsmith D, Palculict B, Hoppes S, Steiner JM, et al. Molecular characterization of the cloacal microbiota of wild and captive parrots. Vet Microbiol. 2010;146(3–4):320–5.
Villers LM, Jang SS, Lent CL, Lewin-Koh SC, Norosoarinaivo JA. Survey and comparison of major intestinal flora in captive and wild ring-tailed lemur (Lemur catta) populations. Am J Primatol. 2008;70(2):175–84.
Schmidt E, Mykytczuk N, Schulte-Hostedde AI. Effects of the captive and wild environment on diversity of the gut microbiome of deer mice (Peromyscus maniculatus). ISME J. 2019;13(5):1293–305.
Edenborough KM, Mu A, Mühldorfer K, Lechner J, Lander A, Bokelmann M, et al. Microbiomes in the insectivorous bat species Mops condylurus rapidly converge in captivity. PLoS ONE. 2020;15(3): e0223629.
Kohl KD, Dearing MD. Wild-caught rodents retain a majority of their natural gut microbiota upon entrance into captivity. Environ Microbiol Rep. 2014;6(2):191–5.
Dhanasiri AKS, Brunvold L, Brinchmann MF, Korsnes K, Bergh Ø, Kiron V. Changes in the Intestinal Microbiota of Wild Atlantic cod Gadus morhua L. Upon Captive Rearing Microb Ecol. 2011;61:20–30.
Aivelo T, Laakkonen J, Jernvall J. Population-and individual-level dynamics of the intestinal microbiota of a small primate. Appl Environ Microbiol. 2016;82(12):3537–45.
Čížková D, Ďureje Ľ, Piálek J, Kreisinger J. Experimental validation of small mammal gut microbiota sampling from faeces and from the caecum after death. Heredity (Edinb). 2021;127(2):141–50.
Björk JR, Dasari M, Grieneisen L, Archie EA. Primate microbiomes over time: Longitudinal answers to standing questions in microbiome research. Am J Primatol. 2019;81(10–11):e22970.
Jackson JA, Hall AJ, Friberg IM, Ralli C, Lowe A, Zawadzka M, et al. An Immunological Marker of Tolerance to Infection in Wild Rodents. PLoS Biol. 2014;12(7):1–13.
Jackson JA, Begon M, Birtles R, Paterson S, Friberg IM, Hall A, et al. The analysis of immunological profiles in wild animals: A case study on immunodynamics in the field vole. Microtus agrestis Mol Ecol. 2011;20(5):893–909.
Wanelik KM, Begon M, Bradley JE, Friberg IM, Jackson JA, Taylor CH, et al. Effects of an IgE receptor polymorphism acting on immunity, susceptibility to infection, and reproduction in a wild rodent. Elife. 2023;16:12.
Peig J, Green AJ. New perspectives for estimating body condition from mass/length data: The scaled mass index as an alternative method. Oikos. 2009;118(12):1883–91.
Rowe FP, Bradfield A, Quy RJ, Swinney T. Relationship Between Eye Lens Weight and Age in the Wild House Mouse (Mus musculus). J Appl Ecol. 1985;22(1):55–61.
Augusteyn RC. Growth of the eye lens: I. Weight accumulation in multiple species. Mol Vis. 2014;20(October 2013):410–26.
Behnke JM, Stewart A, Smales L, Cooper G, Lowe A, Kinsella JM, et al. Parasitic nematodes of the genus Syphacia Seurat, 1916 infecting Cricetidae in the British Isles: the enigmatic status of Syphacia nigeriana. Parasitology. 2022 Jan 5 [cited 2023 Jan 26];149(1):76–94. Available from: https://pubmed.ncbi.nlm.nih.gov/34608855/
Baylis HA. On a Collection of Nematodes from Nigerian Mammals (chiefly Rodents). Parasitology [Internet]. 1928 [cited 2023 Jan 26];20(3):280–304. Available from: https://www.cambridge.org/core/journals/parasitology/article/abs/on-a-collection-of-nematodes-from-nigerian-mammals-chiefly-rodents/1F1E5D14DFF6766E9171D14B7C06B238
Schmidt GD. CRC Handbook of Tapeworm Identification. [Internet]. Vol. 61, https://doi.org/10.1086/415201. Stony Brook Foundation, Inc.; 1986 [cited 2023 Jan 26]. 556–556 p. Available from: https://www.journals.uchicago.edu/doi/10.1086/415201
Mészáros F. Parasitic Nematodes of Microtus arvalis (Rodentia) in Hungary. Parasitologica hungarica. 1977;10:67–83.
Feliu C, et al. Genetic and morphological heterogeneity in small rodent whipworms in southwestern Europe: characterization of Trichuris muris and description of Trichuris arvicolae n. sp. (Nematoda: Trichuridae), Journal of Parasitology. American Society of Parasitologists. 2000;86(3):442–9.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci. 2011;108(Supplement 1):4516–22.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17(1):10–2.
Joshi NA, Fass JN. Sickle: A sliding-window, adaptive, quality-based trimming tool for FastQ files. 2011.
Prjibelski A, Antipov D, Meleshko D, Lapidus A, Korobeynikov A. Using SPAdes De Novo Assembler. Curr Protoc Bioinformatics. 2020;70(1):e102.
Zhang J, Kobert K, Flouri T, Stamatakis A. PEAR: a fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics. 2014;5:614–20.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Mahé F, Rognes T, Quince C, de Vargas C, Dunthorn M. Swarm: robust and fast clustering method for amplicon-based studies. PeerJ. 2014;2(1):e593.
Wang Q, Garrity G, Tiedje J, Cole J. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–7.
Bokulich N, Subramanian S, Faith J, Gevers D, Gordon J, Knight R, et al. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat Methods. 2013;10(1):57–9.
Team RC. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.; 2018.
Mcgee WA, Pimentel H, Pachter L, Wu JY. Compositional Data Analysis is necessary for simulating and analyzing RNA-Seq data. [cited 2023 Jan 12]; Available from: https://doi.org/10.1101/564955
Faith DP. Phylogenetic diversity and conservation evaluation: Perspectives on multiple values, indices, and scales of application. In: Phylogenetic Diversity: Applications and Challenges in Biodiversity Science. Springer International Publishing; 2018. p. 1–26.
Kim BR, Shin J, Guevarra RB, Lee JH, Kim DW, Seol KH, et al. Deciphering Diversity Indices for a Better Understanding of Microbial Communities. J Microbiol Biotechnol. 2017;27(12):2089–93.
Hubert M, Reynkens T, Schmitt E, Verdonck T. Sparse PCA for High-Dimensional Data With Outliers. Technometrics. 2016;58(4):424–34.
Fox J, Weisberg S. An R Companion to Applied Regression. 3rd ed. Thousand Oaks, CA: Sage; 2019.
Venables WN, Ripley BD. Modern Applied Statistics with S. 4th ed. New York: Springer; 2002.
Hoffman GE, Schadt EE. variancePartition: Interpreting drivers of variation in complex gene expression studies. BMC Bioinformatics. 2016;17(1):1–13.
Bates D, Mächler M, Bolker BM, Walker SC. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015 Oct 1;67(1).
Kuznetsova A, Brockhoff PB, Christensen RHB. lmerTest Package: Tests in Linear Mixed Effects Models. J Stat Softw. 2017;82(13):1–26.
Hanning I, Diaz-Sanchez S. The functionality of the gastrointestinal microbiome in non-human animals. Microbiome. 2015;3(1):51.
Faith JJ, Guruge JL, Charbonneau M, Subramanian S, Seedorf H, Goodman AL, et al. The long-term stability of the human gut microbiota. Science. 2013;341(6141):1237439.
Ley RE, Hamady M, Lozupone C, Turnbaugh P, Ramey RR, Bircher JS, et al. Evolution of mammals and their gut microbes. Science. 2008;320(5883):1647–51.
Ley RE, Bäckhed F, Turnbaugh P, Lozupone CA, Knight RD, Gordon JI. Obesity alters gut microbial ecology. Proc Natl Acad Sci U S A. 2005;102(31):11070–5.
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:2423–34.
Curtis JT, Assefa S, Francis A, Kö Hler GA. Fecal microbiota in the female prairie vole (Microtus ochrogaster). 2018 [cited 2022 Nov 29]; Available from: https://doi.org/10.1371/journal.pone.0190648
IUCN Red List - Field Vole [Internet]. [cited 2022 Nov 29]. Available from: https://www.iucnredlist.org/species/13426/115112050
The Wildlife Trusts - Mammals: Field Vole [Internet]. [cited 2022 Nov 29]. Available from: https://www.wildlifetrusts.org/wildlife-explorer/mammals/field-vole
Mammal Society - Discover Mammals: Field Vole [Internet]. [cited 2022 Nov 29]. Available from: https://www.mammal.org.uk/species-hub/full-species-hub/discover-mammals/species-field-vole/
Raulo A, Allen BE, Troitsky T, Husby A, Firth JA, Coulson T, et al. Social networks strongly predict the gut microbiota of wild mice. ISME J. 2021;15(9):2601–13.
Lees H, Swann J, Poucher SM, Nicholson JK, Holmes E, Wilson ID, et al. Age and Microenvironment Outweigh Genetic Influence on the Zucker Rat Microbiome. PLoS ONE. 2014;9(e100916):1–11.
McKenzie VJ, Song SJ, Delsuc F, Prest TL, Oliverio AM, Korpita TM, et al. The effects of captivity on the mammalian gut microbiome. Integr Comp Biol. 2017;57(4):690–704.
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:212.
Li Y, Zhang K, Liu Y, Li K, Hu D, Wronski T. Community Composition and Diversity of Intestinal Microbiota in Captive and Reintroduced Przewalski’s Horse (Equus ferus przewalskii). Front Microbiol. 1821;2019:10.
Bharwani A, Firoz Mian M, Foster JA, Surette MG, Bienenstock J, Forsythe P. Structural & functional consequences of chronic psychosocial stress on the microbiome & host. Psychoneuroendocrinology. 2016;63:217–27.
Bailey MT, Dowd SE, Parry NMA, Galley JD, Schauer DB, Lyte M. Stressor Exposure Disrupts Commensal Microbial Populations in the Intestines and Leads to Increased Colonization by Citrobacter rodentium. Infect Immun. 2010;78(4):1509–19.
Maki KA, Burke LA, Calik MW, Watanabe-Chailland M, Sweeney D, Romick-Rosendale LE, et al. Sleep fragmentation increases blood pressure and is associated with alterations in the gut microbiome and fecal metabolome in rats. Physiol Genomics. 2020;52:280–92.
Voigt RM, Forsyth CB, Green SJ, Mutlu E, Engen P. Circadian Disorganization Alters Intestinal Microbiota. Vol. 9, PLoS ONE. 2014.
Weinstein SB, Mart Inez-Mota R, Stapleton TE, Klure DM, Greenhalgh R, Orr TJ, et al. Microbiome stability and structure is governed by host phylogeny over diet and geography in woodrats (Neotoma spp.). The Proceedings of the National Academy of Sciences. 2021;118(47):e2108787118:1–9.
Kohl KD, Skopec MM, Dearing MD. Captivity results in disparate loss of gut microbial diversity in closely related hosts. Conserv Physiol [Internet]. 2014 [cited 2022 Jun 21];2(1):1–11. Available from: /pmc/articles/PMC4806740/
Gibson KM, Nguyen BN, Neumann LM, Miller M, Buss P, Daniels S, et al. Gut microbiome differences between wild and captive black rhinoceros – implications for rhino health. Sci Rep. 2019;9:7570.
Clayton JB, Vangay P, Huang H, Ward T, Hillmann BM, Al-Ghalith GA, et al. Captivity humanizes the primate microbiome. PNAS. 2016;113(37):10376–81.
Helaszek CT, White BA. Cellobiose uptake and metabolism by Ruminococcus flavefaciens. Appl Environ Microbiol [Internet]. 1991 [cited 2022 Dec 15];57(1):64–8. Available from: https://journals.asm.org/doi/10.1128/aem.57.1.64-68.1991
Conley MN, Wong CP, Duyck KM, Hord N, Ho E, Sharpton TJ. Aging and serum MCP-1 are associated with gut microbiome composition in a murine model. PeerJ. 2016;4:e1854.
Adriansjach J, Baum ST, Lefkowitz EJ, van Der Pol WJ, Buford TW, Colman RJ. Age-related differences in the gut microbiome of rhesus macaques. J Gerontology A Biol Sci Med Sci. 2020;75(7):1293–8.
Koliada A, Syzenko G, Moseiko V, Budovska L, Puchkov K, Perederiy V, et al. Association between body mass index and Firmicutes/Bacteroidetes ratio in an adult Ukrainian population. BMC Microbiol. 2017;17(120):1–6.
Indiani CMDSP, Rizzardi KF, Castelo PM, Ferraz LFC, DarrieuxParisotto MTM. Childhood Obesity and Firmicutes/Bacteroidetes Ratio in the Gut Microbiota: A Systematic Review. Child Obes. 2018;14(8):501–9.
Stojanov S, Berlec A, Štrukelj B. The influence of probiotics on the firmicutes/bacteroidetes ratio in the treatment of obesity and inflammatory bowel disease. Microorganisms. 2020;8:1715.
Indiani CMDSP, Rizzardi KF, Castelo PM, Ferraz LFC, DarrieuxParisotto MTM. Childhood Obesity and Firmicutes/Bacteroidetes Ratio in the Gut Microbiota: A Systematic Review. Childhood Obesity. 2018;14(8):501–9.
The authors wish to thank Rebecca Turner, Lukasz Lukomski, Benoit Poulin, Stephen Price, Sarah Gore, Ed Parker, William Foster, Ann Lowe, Anna Tomason, Maria Capstick, Noelia Dominguez Alvarez and Susan Withenshaw for their work collecting and processing field samples, as well as Lucille Rainbow, Anita Lucaci, John Kenny, Luca Lenzi, Sam Haldenby for microbiome sequencing and processing. They also wish to thank the Forestry Commission for access to the study sites, and the Natural Environment Research Council for funding.
The work in this manuscript was funded by the Natural Environment Research Council, award number NE/L013452/1.
Ethics approval and consent to participate
All procedures were performed with approval from the University of Liverpool Animal Welfare Committee, under a UK Home Office license (PPL 40/3235 to MB), and according to all related welfare guidelines and regulations. Where relevant, all procedures and reporting are in line with ARRIVE guidelines for the reporting of animal experiments.
Consent to publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Fenn, J., Taylor, C., Goertz, S. et al. Discrete patterns of microbiome variability across timescales in a wild rodent population. BMC Microbiol 23, 87 (2023). https://doi.org/10.1186/s12866-023-02824-x
- Wild Populations
- Alpha Diversity
- Beta Diversity