Skip to main content

No evidence for associations between brood size, gut microbiome diversity and survival in great tit (Parus major) nestlings

Abstract

Background

The gut microbiome forms at an early stage, yet data on the environmental factors influencing the development of wild avian microbiomes is limited. As the gut microbiome is a vital part of organismal health, it is important to understand how it may connect to host performance. The early studies with wild gut microbiome have shown that the rearing environment may be of importance in gut microbiome formation, yet the results vary across taxa, and the effects of specific environmental factors have not been characterized. Here, wild great tit (Parus major) broods were manipulated to either reduce or enlarge the original brood soon after hatching. We investigated if brood size was associated with nestling bacterial gut microbiome, and whether gut microbiome diversity predicted survival. Fecal samples were collected at mid-nestling stage and sequenced with the 16S rRNA gene amplicon sequencing, and nestling growth and survival were measured.

Results

Gut microbiome diversity showed high variation between individuals, but this variation was not significantly explained by brood size or body mass. Additionally, we did not find a significant effect of brood size on body mass or gut microbiome composition. We also demonstrated that early handling had no impact on nestling performance or gut microbiome. Furthermore, we found no significant association between gut microbiome diversity and short-term (survival to fledging) or mid-term (apparent juvenile) survival.

Conclusions

We found no clear association between early-life environment, offspring condition and gut microbiome. This suggests that brood size is not a significantly contributing factor to great tit nestling condition, and that other environmental and genetic factors may be more strongly linked to offspring condition and gut microbiome. Future studies should expand into other early-life environmental factors e.g., diet composition and quality, and parental influences.

Introduction

The digestive tract hosts a large community of different microorganisms (i.e., gut microbiome) and is known to be a fundamental part of organismal health and a powerful proximate mechanism affecting host performance [1, 2]. The gut microbiome has been studied across a wide range of animal taxa e.g., humans [3,4,5], fish [6], and economically important species such as poultry [7], and data from wild populations is slowly increasing [8]. Generally, a more diverse gut microbiome is considered beneficial for individual health [9], but there are also community structure effects that define the functionality [10]. For example, laboratory-bred mice with a less diverse gut microbiome have a substantially lower chance of surviving an influenza infection compared to their wild counterparts unless receiving a gut microbiota transplant from their wild counterparts [11, 12]. Moreover, gut microbiome had been linked to host fitness and survival in the Seychelles warbler (Acrocephalus sechellensis). Individuals that harbored opportunistic pathogens, i.e., microbes that usually do not cause disease in healthy individuals, but may become harmful in individuals that are immunocompromised, in their gut microbiome showed higher mortality [13, 14]. Therefore, understanding how gut microbiome affects fitness within and between individuals is necessary for not only understanding species survival but also evolution [15,16,17].

Gut microbiome forms at a young age and remains somewhat stable in adulthood as found for example in laboratory bred mice [18,19,20]. Disruption in the gut microbiome that leads to a microbiome imbalance at a young age could result in both short-term and long-term changes in the gut microbiome [21, 22]. Of the environmental effects, diet [23], including e.g., macronutrient balance (carbohydrates, fats, amino acids) [3, 24] has been concluded to be major determinants of rat and mouse gut microbiome, and this effect has recently been seen in avian models as well [25,26,27,28]. Moreover, macronutrient balance has been linked to intestinal microbiome composition [3, 24] and the functioning of individual immune response [29, 30]. However, as a large part of the prior research has focused strictly on humans or species living in controlled environments in which environmental effects on both the microbiome and host are sidelined [31, 32], many species, including most birds [8], have only started to attract attention [33].

The mechanisms of bacterial colonization of the bird gut are somewhat unique as avian life-histories differ significantly from those of e.g., mammals [34]. In mammals, the offspring are exposed to bacterial colonization during vaginal birth [35] and lactation [36, 37], whereas bird hatchlings are exposed to bacteria first upon hatching [20, 38]. Few studies have investigated the possibility of bacterial colonization in ovo, but results are still lacking [39]. Genetics [40,41,42] as well as the post-hatch environment [20, 43,44,45,46] have a significant effect on the formation of the avian gut microbiome. Once hatched, most altricial birds feed their young, which exposes the hatchlings to various bacteria that originate from the parents i.e., via vertical transmission [47]. It has also been shown that environmental factors are major contributors in the formation of gut microbiome [48,49,50,51], one of these being the rearing environment in the nest [44].

As early-life environment is connected to the establishment of gut microbiome, brood size may affect gut microbiome [52]. Brood size is often associated with parents’ performance and ability to feed their young [53], and the trade-off between offspring quality and quantity has been studied widely [54, 55]. Food quantity per nestling can decrease in enlarged broods, as parents may not be able to fully compensate for the additional amount of food an enlarged brood requires [56, 57]. For example, in great tits (Parus major) it has been shown that nestlings from reduced broods may have a higher body mass [58] and tend to survive better [59]. Importantly, great tit nestling body mass has been connected to gut microbiome diversity and composition: body mass positively correlates with gut microbiome richness [52]. This could imply that good physiological condition and high food availability would allow the host to have a diverse gut microbiome that promotes a healthy gut.

Alterations in early-life gut microbiome could have long-term consequences on individual performance [60], yet such effects have rarely been studied in wild organisms. In wild birds, some bacterial taxa have been linked to better survival. For example, a high abundance of bacteria in the order Lactobacillales of the phylum Firmicutes is related to higher individual fitness in Seychelles warblers [14] and great tits [61]. These bacteria are also known for the benefits for bird health in economically important species such as poultry, in which Lactobacilli are used as probiotics to boost immune functioning [62]. Besides Lactobacillales, gut bacteria belonging to other genera such as Clostridium and Streptococcus are important for the degradation of non-starch polysaccharides and for the synthesis of essential molecules such as the short-chain fatty acids [63, 64]. Short-chain fatty acids are important in host energy metabolism [65] and therefore crucial for performance. Changes in nestling’s early-life gut microbiome could affect such key physiological processes that could influence for example nestling body mass, which is tightly linked to survival to fledging [58, 59]. Because the gut microbiome establishes at a young age and is less plastic later in life [18,19,20], gut microbiome and changes to its richness can have long-term effects on juvenile and adult survival [21, 22]. For example, antibiotic treatment at infancy can affect the expression of genes involved in immune system functioning and lead to long-term effects on host metabolism [20]. Moreover, changes in the rearing environment can affect individual physiology and these effects can carry over to later stages of an individual's life such as survival to fledging and lifetime reproductive success [66].

Here, we use an experimental approach to investigate whether brood size manipulation influenced wild great tit nestlings’ bacterial gut microbiome diversity on day 7 post-hatch. We also investigated whether brood size influenced nestling body mass on day 7 or on day 14 post-hatch, and if the gut microbiome predicts short-term (i.e., survival to fledging) and mid-term (i.e., apparent juvenile) survival. The great tit is a well-studied species in the fields of ecology and evolution, and it is easy to monitor in the wild due to its habit of breeding in nest boxes. Great tit nestlings’ gut microbiome undergoes profound shifts during early life [52], and it has been linked to nestling natal body mass and body size [52, 61], yet studies focusing on gut microbiome associations with survival are still scarce. Here, we manipulated wild great tit broods by reducing or enlarging the original brood size in order to analyze if this affected the gut microbiome. In large broods, nestlings need to compete for their food more [67, 68], and the lower food availability could result in a lower gut microbiome diversity. This might impair nestling body mass and fitness prospects [13, 52]. We used a partial cross-fostering design that enabled us to disentangle the relative contributions of genetic background, early maternal effects, and rearing environment such as parents, nest and nestmates on gut microbiome. Furthermore, we used an unmanipulated control group in which no nestling was cross-fostered to control for the possible effects of moving the nestlings between nests. For example, early human handling such as marking and weighing at day 2 post-hatch could influence gut microbiome later on. We hypothesized that (1) in reduced broods nestlings would have a higher body mass, (2) in reduced broods nestling gut microbiome would be more diverse than in enlarged broods, and (3) higher gut microbiome diversity on day 7 post-hatch would increase survival to fledging and potentially reflect apparent juvenile survival. Such knowledge could provide new information about gut microbiome in wild passerine bird population and how the early-life environment may associate with nestling gut microbiome, body mass, and short-term and mid-term survival.

Methods

Study area and species

The great tit is a small passerine bird, which breeds in secondary holes and artificial nest-boxes, making it a suitable model species. Great tits breed throughout Europe and inhabit parts of Northern Africa and Asia as well, and the breeding areas differ in environment and diet [69]. In Finland the great tit is a common species with an estimate of 1.5 to 2 million breeding pairs. They lay 6 to 12 eggs between April and May and the female incubates the eggs for 12–15 days. The nestlings fledge approximately 16 to 21 days after hatching. The study was conducted during the breeding season (May–July 2020) on Ruissalo island (60°25′59.99″ N 22°09′60.00″ E). Ruissalo island habitat is a mostly temperate deciduous forest and meadows, and some areas have small patches of coniferous trees.

Brood size manipulation experiment

Nest boxes were first monitored weekly and later daily when clutches were close to the estimated hatching date. Brood size manipulation took place on day 2 after hatching. Increases in great tit brood size can lead to lowered weight in both the nestlings and adults [70,71,72,73,74,75], and our decision on the number (i.e., + 2 or − 2) of manipulated nestlings (i.e, + 2 or − 2) followed the cited studies. We had four treatment groups (see Fig. 1): in the ‘enlarged group (henceforward called E)’, we increased the brood size by two individuals that were taken from a ‘reduced brood’. Correspondingly, in the ‘reduced group (henceforward called R)', we decreased the brood size by two individuals, that were added to the enlarged broods. In the ‘control group (henceforward called C)', we swapped nestlings between nests but did not change the brood size. And lastly, in the ‘unmanipulated control group (henceforward called COU)’, we only weighed and collected fecal samples on day 7 but did not move the nestlings between nests. We also moved nestlings between the reduced nests to ensure that all nests except for COU had both original and fostered nestlings. Control nests were used to control for potential cross-fostering effects unrelated to brood size. Additionally, in the unmanipulated control group nestlings were not moved or weighed on day 2 in order to control for any handling effects per se. This study design enabled us to test the potential impacts of handling nestlings and swapping the nest early after hatching. We aimed to move approximately half of the chicks in the manipulated nests, so that the number of original and the fostered nestlings would be the same in each nest after manipulation.

Fig. 1
figure 1

Brood size manipulation experiment schematic diagram. 2-day-old nestlings were moved between nestboxes to enlarge or to reduce original brood size (an example with brood size of seven is given). Some nests were kept as control nests (nestlings were moved but brood size remained the same) and some were kept as unmanipulated control nests (nestlings were not moved at all to test whether early-life handling affects gut microbiome). The original brood size varied between nests

Before they were moved, nestlings were weighed using a digital scale with a precision of 0.1 g and identified by clipping selected toenails. We aimed to add/remove nestlings that were of similar weight to avoid changing the sibling hierarchy in the brood. The moving procedure was performed as quickly as possible to minimize the risk of stress and the nestlings were kept in a warmed box during transportation. For each pair of nests in the brood size manipulation experiment, we selected nests that had a similar hatching date. In case of uneven number of nests hatching within a day, one or three nest(s) was/were allocated to the COU group. To avoid potential bias from hatching date, we allocated nests in any given day evenly to each treatment. We also checked that the treatments had an equal brood size on average i.e., we did not want to only reduce the larger clutches and enlarge the smaller clutches. These is also a significant bias towards COU nests being later in the season on average (Table 1).

Table 1 (A) Brood size before and after manipulation, (B) hatching date across treatments

Fecal sample collection

To study the effects that brood size may have on the nestling gut microbiome and its links to individual nestling body mass, survival to fledging and apparent juvenile survival, we used a subset of data from a larger experiment (Cossin-Sevrin et al., unpublished data). In this subset, we use individuals from which fecal samples were collected on day 7 after hatching and analyzed for microbiome diversity and composition (C = 23 nestlings/15 nests, COU = 22/13, E = 23/15, R = 24/16) We aimed to collect two samples (one from original and one from foster nestlings) per nest. Fecal samples from the nestlings were collected gently by stimulating the cloaca with the collection tube. Samples were collected straight into a sterile 1.5 ml Eppendorf tube to avoid possible contamination of the sample. At time of sampling, each nestling was weighed (0.1 g), and the nestlings were ringed for individual identification using aluminum bands. The samples were stored in cool bags onsite and afterwards moved to a -80 °C freezer for storage until DNA extraction.

Apparent juvenile survival

We monitored all study nests until fledging to measure short-term survival. On day 14 post-hatch, the sampled nestlings were weighed, and wing-length was measured to detect if the manipulation had any effects on nestling growth. Nests were subsequently monitored for fledging success. Additionally, we monitored our study population for apparent juvenile survival (i.e., mid-term survival) after the breeding season (i.e., approximately 3 months after fledging) to assess the association between gut microbiome and post-fledging survival. We captured juvenile great tits by mist netting during the autumn–winter 2020 at six different feeding stations that had a continuous supply of sunflower seeds and suet blocks. Feeding stations were located within the previously mentioned nest box population areas. For each site mist netting with playback was conducted on three separate days during October–November 2020 for three hours at a time, leading to a total of 69 h of mist netting. A total of 88 individuals from the brood size manipulation experiment were caught, and the caught juvenile great tits were weighed, and wing length was measured. Our catching method provides an estimate of post-fledging survival yet, it could be slightly biased based by dispersal. In a previous study in our population [76], none of the birds ringed as nestlings were recaptured outside the study area, suggesting that dispersal is likely limited.

DNA extraction and sequencing

We chose two samples per nest for DNA extraction, yet in such a way that both fledged and not-fledged nestlings would be included in the dataset. DNA was extracted from nestling fecal samples using the Qiagen QIAamp PowerFecal Pro DNA Kit (Qiagen; Germany) following the manufacturer’s protocols. Additionally, we included negative (RNAse and DNAse free ddH2O) controls to control for contamination during DNA extraction and additional controls to confirm successful amplification during PCR. A short fragment of hypervariable V4 region in the 16S rRNA gene was amplified using the purified DNA samples as template with the following primers: 515F_Parada (5’-GTGYCAGCMGCCGCGGTAA-3’) and 806R_Apprill (5’-GGACTACNVGGGTWTCTAAT-3’) [77, 78]. PCRs were performed in a total volume of 12 µL using MyTaq RedMix DNA polymerase (Meridian Bioscience; Cincinnati, OH, USA). The PCR cycling conditions were as follows: first, an initial denaturation at 95 °C for 3 min followed by 30 cycles of 95 °C for 45 s, 55 °C for 60 s, and 72 °C for 90 s, and finished with a 10-min extension at 72 °C. After the first round of PCR, a second round was conducted to apply barcodes for sample identification [79]. For this, PCR cycling conditions were as follows: first, an initial denaturation at 95 °C for 4 min followed by 18 cycles of 98 °C for 20 s, 60 °C for 15 s, and 72 °C for 30 s, and finished with a 3-min extension at 72 °C. We performed replicate PCR reactions to control for errors during the amplification. Further on, the PCR products were measured for DNA concentration with Quant-IT PicoGreen dsDNA Assay Kit (ThermoFischer Scientific; Waltham, MA, USA) and for quality with TapeStation 4200 (Agilent; Santa Clara, CA, USA). The samples from each of the PCR replicates were pooled equimolarly creating two separate pools and purified using NucleoMag NGS Clean-up and Size Select beads (Macherey–Nagel; Düren, Germany). Finally, pooled samples were sequenced (2 × 300 bp) on the Illumina MiSeq platform (San Diego, CA, USA) at the Finnish Functional Genomic Center at the University of Turku (Turku, Finland).

Sequence processing

All statistical analyses were performed with R (v. 4.11.0; R Development Core Team 2021) unless otherwise stated. The demultiplexed Illumina sequence data was first processed with Cutadapt version 2.7 [80] to remove locus-specific primers from both R1 and R2 reads. Then, the DADA2 pipeline (v. 1.24.0; [81]) was used to filter the reads based on quality, merge the paired-end (R1 and R2) reads, to define the single DNA sequences i.e., Amplicon Sequence Variants (henceforward ASV), and to construct a ‘seqtab’. Seqtab is a matrix also known as otutable or readtable: ASVs in columns, samples in rows, number of reads in each cell, using default parameter settings. In total, our seqtab consisted of 6,929,537 high-quality reads. Reads were assigned to taxa against the SILVA v132 reference database [82] resulting in 8658 ASVs. To control for contamination, negative DNA extraction and PCR controls were used to identify contaminants (60 ASVs) using the decontam package (v. 1.12; [83]) and all were removed from the dataset. Sequencing runs (replicate PCR’s) were merged using the phyloseq package (v. 1.32.0) and non-bacterial sequences (mainly Chlorophyta) were removed from the data as they were not of interest in this study resulting in a total of 6566 ASVs (a total of 4,107,028 high-quality reads in all samples; mean per sample: 15,155.085; mean range per sample: 0–97,264). Singleton reads were removed from the dataset by the DADA2 pipeline. Data was further analyzed with the phyloseq package (v. 1.32.0; [84]), and the microbiome package (v. 1.18.0; [85]) and visualized with the ggplot2 package (v. 3.3.6; [86]).

The final dataset contained 92 samples from great tit nestlings resulting in a total of 3,161,696 reads (mean per sample: 34,366.261; mean range per sample 108 – 189,300 reads), which belonged to 6,505 ASVs. The dataset was then rarefied for alpha diversity analyses at a depth of 5000, as this was where the rarefaction curves plateaued (see Additional file 2). The rarefied dataset contained 4,791 ASVs in 88 samples. For beta diversity, the unrarefied dataset was used after confirming that the beta diversity statistics were quantitatively similar for the rarefied and unrarefied datasets. Bacterial relative abundances were summarized at the phylum and genus level and plotted based on relative abundance for all phyla and genera. A Newick format phylogenetic tree with the UPGMA algorithm to cluster treatment groups together was used to visualize sample relatedness (see Additional file 3) and was constructed using the DECIPHER (v. 2.24.0; [87]), phangorn (v. 2.8.1; [88]), and visualized with ape (v. 5.6-2; [89]), and ggtree (v. 3.4.0; [90]) packages.

Statistical analyses

Nestling body mass

First, to analyze whether brood size manipulation affected nestling body mass in the C, E, and R treatment groups, we ran two linear mixed-effects models with the lme4 package (v. 1.1-29; [91]). In these models we used either body mass on day 7 or 14 as the dependent variable and brood size manipulation treatment, hatching date, body mass on day 2 post-hatch and original brood size as predicting variables. Hatching date is used as a predicting variable because it is known to affect nestling body mass during the breeding season [92] and there were differences in hatching date between the COU and other treatment groups (see Table 1). We included the interaction between original brood size and brood size manipulation treatment in both models as the effect of manipulation may depend on the original brood size. For example, there could be stronger effect of enlargement in already large broods. Nest of origin and nest of rearing were used as random intercepts to control for the non-independence of nestlings sharing the same original or foster nests. Here, we did not include the COU group in the analysis because we wanted to measure the effects of treatment on nestling body mass, and only enlarged, reduced or control broods’ nestlings were moved between nests.

Second, to analyze whether the actual brood size affected nestling body mass, we ran two models where we used it as a continuous dependent variable to explain body mass either on day 7 or on day 14 post-hatch. Hatching date and body mass on day 2 post-hatch were used as predicting variables and nest of origin and nest of rearing as random intercepts to control for the non-independency of samples. We included the interaction between manipulated brood size and hatching date in the models because the effect of brood size may depend on the hatching date. For example, hatching date can reflect environmental conditions and large broods may perform poorly late in the season due to poorer food availability. The COU group was initially excluded from this model to see which of the two random effects, nest of origin or nest of rearing, explained a larger portion of variation in the treatment groups. In the COU group, nest of origin and nest or rearing were the same, which meant we could not include both random effects in models where all treatment groups were present due to the model failing to converge. Nest of origin explained more of the variation in the first model (see Additional file 4) and therefore, we used it in the full models with all treatment groups: C, COU, E and R. In these models, nestling body mass either on day 7 and or on day 14 post-hatch was used as a dependent variable and manipulated brood size as the explanatory variable. Hatching date and body mass on day 2 post-hatch were set as predicting variables. Nest of rearing was used as a random intercept to control for the non-independence of nestlings sharing the same foster nests. The significance of factors included in the models were tested using the F-test ratios in analysis of variance (Type III ANOVA).

Alpha diversity

For alpha diversity analyses, which measures within-sample species diversity, we ran two linear mixed-effects models with the lme4 package (v. 1.1-29; [91]) to measure if either brood size manipulation or manipulated brood size as a continuous variable were associated with gut microbiome diversity. We used two alpha diversity metrics: the Shannon Diversity Index, which measures the number of bacterial ASVs and their abundance evenness within a sample, and Chao1 Richness, which is an estimation of the number of different bacterial ASVs in a sample. Both metrics were used to check if alpha diversity results were consistent across different metrics. Each diversity index was used as the dependent variable at a time and either brood size manipulation treatment or manipulated brood size as a predicting variable. In both models we included original brood size, weight on day 7 post-hatch and hatching date as covariates. We included interaction between brood size manipulation treatment and original brood size as there could be a stronger effect of enlargement in initially large broods. We also included interaction between manipulated brood size and weight on day 7 post-hatch because effect of brood size on microbiome may depend on nestling weight. We also tested whether alpha diversity predicted weight on day 7 post-hatch, as weight and gut microbiome diversity have been connected in previous studies. In this analysis we used weight on day 7 post-hatch as the dependent variable and alpha diversity (Shannon Diversity Index and Chao1 Richness), treatment and hatching date as predicting variables and nest of rearing as the random effect. In these sets of models, we first excluded the COU group to see which of the two random effects, nest of origin or nest of rearing, explained a larger proportion of variation in the treatment groups. Nest of rearing explained more of the variation in this model (see Additional file 4) and therefore, we used it in the full model with all treatment groups: C, COU, E and R. The significance of factors included in the models were tested using the F-test ratios in analysis of variance (ANOVA).

Short-term survival

To explore whether alpha diversity associated with survival to fledging (i.e., short-term survival) and with apparent juvenile survival in Autumn 2020 (i.e., mid-term survival), we used generalized linear models with binomial model (v. 1.1-29; lme4 package, [91]), and then tested the significance of factors with type 2 ANOVA from the car package (v. 3.0-13; [93]). Type 2 ANOVA was used because the model did not contain interaction between predicting and there was no order between covariates, as they could not be ranked. Survival to fledging and recapture in Autumn 2020 were used as the binomial response variable (yes–no) in each model. Alpha diversity (Shannon Diversity Index and Chao1 Richness) was the main predicting variable, and weight on day 7 post-hatch (same time as sampling the fecal gut microbiome), hatching date and manipulated brood size were included as covariates in the model. We did not include brood size manipulation treatment in the survival models as not enough birds from each treatment group were recorded for fledging and juvenile survival. Moreover, we excluded random effects from this model as the model failed to converge. 65 nestlings fledged successfully, while 8 nestlings were found dead in nest boxes. For 15 nestlings we had no fledging record, so these were excluded from the survival to fledging analysis. In apparent juvenile survival, 19 birds out of 92 (with data on microbiome diversity) were recaptured as juveniles. For all analyses, the R package car (v. 3.0-13; [93]) was used to test Variance Inflation Factors (VIFs) and the package DHARMa (v. 0.4.5; [94]) to test model diagnostics for linear mixed-effects and generalized linear models.

Beta diversity

For visualizing beta diversity, i.e., the similarity or dissimilarity between the treatment group gut microbiomes, non-metric multidimensional scaling (NMDS) was used with three distance matrices: Bray–Curtis [95], weighted UniFrac, and unweighted UniFrac [96]. Permutational multivariate analysis of variance (PERMANOVA) using the Euclidean distance matrix and 9999 permutations was tested with the R package vegan (adonis2 function; v. 2.6-2; [97]) to investigate if any variables affected to the variation in gut microbiome composition. Nest of rearing was set as a blocking factor in the PERMANOVA to control for the non-desirable effects of the repeated sampling of foster siblings. The test for homogeneity of multivariate dispersions was used to measure the homogeneity of group dispersion values. We used the phyloseq package (v. 1.32.0; [84]) to run a differential abundance analysis with a significance cut-off p < 0.01 to test the differential abundance of ASVs between the treatment groups.

Results

The effects of brood size manipulation on nestling body mass

Brood size manipulation did not significantly affect nestling body mass on day 7 post-hatch (ANOVA: F2, 25.832 = 0.441, p = 0.648; see Additional file 5). Moreover, there was no significant interaction between brood size manipulation and original brood size (ANOVA: F2, 24.610 = 0.678, p = 0.517; see Additional file 5). On day 14 post-hatch, brood size manipulation did not significantly affect nestling body mass (ANOVA: F2, 24.335 = 0.831, p = 0.448; see Additional file 5). However, body mass increased with increasing hatching date (ANOVA: F1 24.070 = 13.367, p = 0.001; see Additional file 5). Next, we did not find any significant associations between manipulated brood size and nestling body mass (ANOVA for weight on day 7: F1, 35.149 = 1.777, p = 0.191; ANOVA for weight on day 14: F1, 29.491 = 2.156, p = 0.153; see Additional file 6). Nest of origin explained a larger proportion of the variation in weight than the nest of rearing on both day 7 (nest of origin 41.1% and nest of rearing 24.4%) and day 14 (nest of origin 65.5% and nest of rearing 21.9%) post-hatch, but this result was not statistically significant (Pr > χ2 = 1) (see Additional file 4).

Alpha diversity

As 7-day-old nestlings, most bacterial taxa belonged to the phyla Proteobacteria, Firmicutes, and Actinobacteria (Fig. 2).

Fig. 2
figure 2

Bacterial relative abundances on Phylum level across the four treatment groups. Each bar represents an individual sample. Treatment groups are control (C), unmanipulated control (COU), enlarged (E), and reduced (R). N = 88 samples divided into treatment groups as follows: C = 23, COU = 21, E = 20, R = 24. Phyla with less that 10% in relative abundance is collapsed into the category “< 10% abundance”

Brood size manipulation did not significantly influence alpha diversity (Shannon Diversity Index) (ANOVA: F3, 47.488 = 1.026, p = 0.390, Fig. 3; see Additional file 7). Moreover, original brood size (ANOVA: F1, 50.269 = 0.388, p = 0.536; see Additional file 7), weight on day 7 post-hatch (ANOVA: F1, 80.551 = 0.003, p = 0.959; see Additional file 7), and hatching date (ANOVA: F1, 50.276 = 1.073, p = 0.305; see Additional file 7) did not significantly associate with alpha diversity. There was no significant interaction between brood size manipulation and original brood size (ANOVA: F3, 48.053 = 0.126, p = 0.944; see Additional file 7). Results for Chao1 Richness were quantitatively similar: brood size manipulation did not affect alpha diversity (ANOVA: F3, 45.936 = 0.358 p = 0.784, Fig. 3; see Additional file 7). Nest of rearing explained a larger proportion of the observed variance in alpha diversity (27.7%) than nest of origin (10.8%), but the result was not statistically significant (Pr > χ2 = 1) (see Additional file 4).

Fig. 3
figure 3

The gut microbiome alpha diversity of 7-day-old great tit nestlings across the four treatment groups visualized with two diversity metrics: A Shannon Diversity Index and B Chao1 Richness. The black dots represent each observation within a treatment group. The whiskers represent 95% confidence intervals. Treatment groups are control (C), unmanipulated control (COU), enlarged (E), and reduced (R). N = 88 samples divided into treatment groups as follows: C = 23, COU = 21, E = 20, R = 24

Next, we tested whether the manipulated brood size as a continuous variable was associated with alpha diversity (Shannon Diversity Index), but found no significant association (ANOVA: F1, 63.001 < 0.001, p = 0.984; see Additional file 8) in this analysis either. Weight on day 7 post-hatch (ANOVA: F1, 82.840 = 0.015, p = 0.903; see Additional file 8) and hatching date (ANOVA: F1, 59.734 = 0.137, p = 0.713; see Additional file 8) did not correlate with alpha diversity in this model either. There was no significant interaction between manipulated brood size and weight on day 7 post-hatch (ANOVA: F1, 82.702 < 0.000, p = 0.998; see Additional file 8). Results for Chao1 Richness were quantitatively similar (ANOVA: F1, 65.064 = 0.246, p = 0.622; see Additional file 8): manipulated brood size did not affect alpha diversity, and neither did weight on day 7 post-hatch (ANOVA: F1 83.513 = 0.690, p = 0.409; see Additional file 8) nor hatching date (ANOVA: F1 57.110 = 1.133, p = 0.292; see Additional file 8).

Alpha diversity and short/mid-term survival

Next, we explored whether alpha diversity (Shannon Diversity Index and Chao1 Richness) contributed to predicting short/mid-term survival (survival to fledging and apparent juvenile survival). Survival to fledging was not predicted by alpha diversity (Shannon Diversity Index: χ2 = 0.010, df = 1, p = 0.923; see Additional files 9 and 10), manipulated brood size (χ2 = 0.090, df = 1, p = 0.764; see Additional file 9), weight on day 7 post-hatch (χ2 = 0.388, df = 1, p = 0.533; see Additional file 9) or hatching date (χ2 = 0.438, df = 1, p = 0.508; see Additional file 9).

Apparent juvenile survival was not significantly associated with alpha diversity (Shannon Diversity Index: χ2 = 1.916, df = 1, p = 0.166; see Additional file 9 and Additional file 10). Moreover, there was no significant interaction between alpha diversity and manipulated brood size (χ2 = 1.268, df = 1, p = 0.260; see Additional file 9). However, apparent juvenile survival was negatively associated with hatching date (χ2 = 4.654, df = 1, p = 0.031; see Additional file 9). Additional analyses to check for the consistency of results were tested the following way: survival to fledging with nestlings from the COU group removed and apparent juvenile survival without the nestlings with no recorded survival for fledging (see methods). These results were quantitatively similar as in the whole dataset for both Shannon Diversity Index (survival to fledging: χ2 = 2.285, df = 1, p = 0.131; apparent juvenile survival: χ2 = 1.515, df = 1, p = 0.218; see Additional file 11) and Chao1 Richness (survival to fledging: χ2 = 0.665, df = 1, p = 0.415; apparent juvenile survival: χ2 = 2.654, df = 1, p = 0.103; see Additional file 11).

Beta diversity

Non-metric multidimensional scaling (NMDS) using weighted and unweighted UniFrac and Bray–Curtis dissimilarity did not show clear clustering of samples based on brood size manipulation treatment (see Additional file 3). The test for homogeneity of multivariate dispersions supported the visual assessment of the NMDS (Betadispersion9999 permutations: F3, 0.069 = 0.650, p < 0.001; see Additional file 12). Pairwise PERMANOVA further indicated that the treatment (PERMANOVA: R2 = 0.061, F = 1.951, p = 0.278; see Additional file 12), weight on day 7 post-hatch (PERMANOVA: R2 = 0.015, F = 1.387, p = 0.091) or hatching date (PERMANOVA: R2 = 0.0232, F = 2.214, p = 0.993) did not significantly contribute to the variation in gut microbiome composition between the treatment groups. Differential analysis of ASV abundance between the treatment groups showed that there is variation in taxa abundance. E group showed higher taxa abundance when compared to COU and C groups and was slightly higher than the R group. C and COU groups were generally lower in taxa abundance than R and E groups, and COU group showed lower abundance than the other groups in each comparison (see Additional file 13).

Discussion

In this study, we investigated the associations between great tit nestling gut microbiome, brood size, and nestling body mass by experimentally manipulating wild great tit broods to either reduce or enlarge the original brood size. The results show that even though there was individual variation in the nestling gut microbiome (Fig. 2), brood size did not significantly contribute to gut microbiome diversity. Neither did gut microbiome diversity explain short-term (survival to nestling) nor mid-term (apparent juvenile) survival. Body mass was also not significantly affected by brood size manipulation. The COU group that functioned as a control for moving and handling effects, did not differ in this respect from the other groups. This suggests that human contact or handling nestlings 2 days post-hatch did not influence nestling gut microbiome or body mass. The partial cross-fostering design enabled us to disentangle the relative contributions of rearing environment (i.e., parents, nest and nestmates) from genetic, prenatal such as maternal allocation to egg, and early post-natal effects such as feeding up to day 2. Nest of rearing seemed to explain more of the variation in nestling gut microbiome diversity than the nest of origin (although not statistically significant), which follows previous studies. Contrastingly, nest of origin seemed to be a stronger contributor than the nest of rearing on nestling body mass on day 7 and day 14 post-hatch. This result was also not statistically significant.

Brood size manipulation and nestling body mass

First, we explored whether brood size was associated with nestling body mass, as such changes may explain the underlying patterns in gut microbiome [52]. Against our hypothesis, we found no significant association between nestling body mass and brood size: neither reduction nor enlargement of the broods resulted in significant body mass differences in the nestlings on day 7 and day 14 post-hatch. While the result is supported by some studies in which associations between nestling body mass and brood size have been tested [61, 98], the majority of the literature shows that brood size negatively correlates with nestling body mass: in larger broods nestlings are generally of lower mass [52, 53, 57, 67, 99,100,101,102,103,104].

There are a few possible explanations why brood size manipulation did not affect nestling body mass. Firstly, if environmental conditions were good, parents may have been able to provide enough food even for the enlarged nests and thus, variance in brood size may not result in differences in nestling body mass between reduced and enlarged nests. In that case the number of nestlings transferred between enlarged and reduced nests should probably have been larger to create differences in nestling body mass between the two treatments. Still, we think that the decision to transfer + 2/− 2 was reasonable since it was based on extensive evidence from previous studies [103]. Secondly, it could be that the enlarged brood size negatively influences some other physiological traits while body mass was retained at the expense of these other traits e.g., immune system functioning [105, 106]. Moreover, our analysis showed that hatching date had a significant effect on nestling body mass: nestlings that hatched later in the season were of lower weight. This could be a result of changes in the food items that great tits use, changes in temperature conditions or in parental investment during the breeding season. As the season progresses, the abundance of insect taxa varies, and this can result in changes in nutrient rich food [103, 107]. For example, great tits can select certain lepidopteran larvae that vary in their abundance during the great tit breeding season [108]. Thirdly, it could be that the change in brood size was influencing the parents’ condition instead of the nestlings [109, 110]. In enlarged broods, parents are required to forage more which can lead to higher energy expenditure and increased stress levels in parents [72, 73, 109].

Brood size manipulation and gut microbiome

We found large inter-individual differences in gut microbiome diversity, yet this variation was not explained by brood size or nestling body mass. It is possible that brood size did not result in differences in food intake. For example, parents were likely able to provide an equivalent amount of food, given that body mass was not significantly affected by the brood size manipulation. Therefore, brood size manipulation did not affect gut microbiome diversity through differences in nutrient uptake. Alternatively, in this study, fecal sampling took place 5 days after the initial brood size manipulation (day 2 post-hatch). It could be that sampling on a later date or at multiple timepoints [61, 111] would have led to different results. Firstly, the time interval may not have been long enough to detect effects of the brood size manipulation. Secondly, it has been shown in previous studies that the nestling gut microbiome undergoes profound shifts at the nestling stage: overall gut microbiome diversity decreases but relative abundance in some taxa increases [52]. We suggest that fecal samples could be collected on multiple days post-hatch to understand the potential day to day changes in the nestling gut microbiome.

Our results suggest that the variance in gut microbiome is a result of other factors than those linked to brood size. Firstly, one of these factors could be diet (i.e., food quality) which has gained attention in gut microbiome studies during the past years [25, 27, 112,113,114,115]. The overall diversity in gut microbiome could be explained by adaptive phenotypic plasticity because it is sensitive to changes in the environment e.g., changes in diet [116, 117]. The food provided by the parents can vary between broods in different environments [118], and this variation in diet can lead to differences in gut microbiome diversity [114,115,116,117,118,119]. For example, abundance in certain dietary items such as insects or larvae can result in lower gut microbiome diversity than other dietary items [113,114,115,116]. As great tits have been reported to adapt their diet along the breeding season due to changes in insect taxa frequency [103, 107] this could affect the between-nestling and between-nest gut microbiome diversity. However, using wild bird populations in gut microbiome studies limits the ability to control the consumed dietary items because parents may use variable food resources and there can be variance in dietary between sexes and even individuals. Visual assessment of dietary items [116] and metabarcoding could be of use here as they enable the identification of food items on genus and even species level from e.g., fecal samples [119].

Secondly, breeding habitat may lead to differences in gut microbiome diversity [120]: adult birds living in deciduous forests have shown to harbor different gut microbiome diversity than their counterparts living in open forested hay meadows. Here, we used a cross-fostering design to study if the rearing environment contributed to the variation in gut microbiome diversity: Our study indicated that the nest of rearing seemed to explain more of the gut microbiome variation than the nest of origin (although not significant), which follows some previous results [43, 44, 52]. For example, a study with great and blue tit (Cyanistes caerulaeus) nestlings showed that the nest of rearing contributed more to the gut microbiome than the nest of origin [43], and another study with the brown-headed cowbird (Molothrus ater) concluded that the sampling locality had a significant contribution to the gut microbiome [44]. Teyssier et al. [52] conducted cross-fostering at day 8 post-hatch in great tits and found that the nest of rearing influenced the gut microbiome more than the nest of origin. Additionally, parents can pass down their bill and feather microbiome through vertical transmission, which could influence nestling gut microbiome [20].

Results from beta diversity analysis were similar to that of alpha diversity: brood size manipulation did not contribute to the variation in gut microbiome composition. Overall, variation in gut microbiome composition could be a result of different genetic and environmental contributors. Firstly, great tit nestling gut microbiome composition could be explained by underlying genetic effects that we did not measure in this study. Phylosymbiosis i.e., the matching of gut microbiome composition to host genetic structure, could be explained by underlying genetics that may translate into physiological differences that affect the gut microbiome e.g., founder effects or genetic drift [121]. Davies et al. [14] found that MHC genes correlate with gut microbiome composition: the expression of specific alleles in the MHC genes was connected to the abundance of specific bacterial taxa such as Lactobacillales and Bacteroidales that influenced host health. In a study by Benskin et al. [41] captive zebra finches (Taeniopygia guttata) showed significant variation in gut microbiome composition between individuals even though their diet and housing conditions were standardized. The study suggested that individual homeostatic mechanisms linking to naturally occurring differences in individual gut microbiome could be why gut microbiome composition varied even with standardized housing conditions [41]. Secondly, gut microbiome composition could have been affected by the same environmental effects that may have linked to the variation in gut microbiome diversity: diet and feeding behavior [115, 116].

Differential analysis of ASV abundance showed variation in differential abundance of taxa between the treatment groups. However, several ASVs were not taxonomically assigned beyond family level making it difficult to draw conclusions about the significance of these results. All treatment groups had taxa belonging to the order Firmicutes, Proteobacteria and Actinobacteria, which was to be expected because they are usually the most core phyla in passerine gut microbiomes [33]. Nestlings belonging to E, R or C group showed higher taxa abundance than the COU group in each comparison. This result could be a result of the COU nestlings generally hatching later in the season and potentially having a less diverse diet [103, 107]. Of the E, R and C groups, C group was less abundant than E and R groups. Both E and R group showed high taxa abundance, which is interesting because we hypothesized that nestlings belonging to the E group would potentially experience less parental investment per nestling and have lower gut microbiome diversity and therefore, be less abundant [56, 57, 67, 68]. We did not observe differential abundance in e.g., the order Lactobacillales which would have been of interest, because the order hosts taxa that are beneficial for gut microbiome health [14, 62]. The genus Staphylococcus was differentially abundant in the E group, but not in the other groups. Staphylococcus is a gram-positive genus of bacteria and known to cause infections in its host species [122]. Curiously, the COU group was differentially abundant in the genus Dietzia, which is a human pathogen [123].

Gut microbiome and short-term and mid-term survival

Our results showed that gut microbiome diversity and brood size were not significantly associated with short-term (survival to fledging) or mid-term (apparent juvenile) survival. However, while a more diverse gut microbiome is considered a possible indicator of a healthy gut microbiome, the effects of the gut microbiome on the host health may often be more complex and related to specific taxa [9, 10]. For example, Worsley et al. [13] did not find a correlation between body condition and gut microbiome diversity, yet they found that specific taxa in the gut microbiome linked with individual body condition and survival. Not only environment, but also genetic background of the individual may contribute to gut microbiome and survival. In a study by Davies et al. [14], Ase-ua4 allele of the MHC genes was linked to lower gut microbiome diversity and it was suspected that the variation in the MHC genes could affect the sensitivity to pathogens that could lead to variation in gut microbiome diversity and eventually, host survival.

To gain a better understanding of gut microbiome diversity and the contribution of different taxa to host survival, functional analyses of the gut microbiome should be included in gut microbiome studies. Different bacterial taxa can have similar functions in the gut microbiome [5, 124] and therefore, the absence of some taxa may be covered by other functionally similar taxa, resulting in a gut microbiome that is functionally more stable [125]. Similarity in functions may also contribute to host’s local adaptation e.g., to the changes in the host’s early-life environment [124]: changes in brood size or dietary items could result in variation in the gut microbiome diversity, yet there may be no effects on host body condition.

The lack of association between brood size, nestling size and survival contrasts with previous studies, but it should be noted that the majority of previous studies have been done with adult birds and not nestlings. Because nestling gut microbiome is still quite flexible compared to that of the adults [20], it is possible that our experiment did not result in a strong enough effect on the gut microbiome. In future studies, it would be important to study the parents as well as it could be more likely to find an association between adult microbiome and fitness than with nestling gut microbiome and survival. Also, our sample size in the survival analyses was small, and it is hard to determine if the result was affected by the sample size. Firstly, nestling survival is often found to correlate with brood size and more specifically, with fledging mass and in particular, the ability to forage for food [61, 126]. Intra-brood competition may explain survival to fledging, as competition between nestlings can limit food availability and thus, leading to lower nestling body condition [68, 127]. A study with blackbirds (Turdus merula) showed that nestling body mass explained juvenile survival [128], and similar results have been shown with great tits and collared flycatchers (Ficedula albicollis; [31]). Contrastingly, Ringsby et al. [129] observed that in house sparrows (Passer domesticus) juvenile survival was independent of nestling mass and brood size. Moreover, natal body mass is often positively correlated with survival to fledging and juvenile survival as heavier nestlings are more likely to be recruited [92, 130, 131], yet we failed to demonstrate this in our study. Hatching date is also often positively correlated with fledging success [132] yet we did not find this association in our study, but instead found a significant association between hatching date and apparent juvenile survival.

Conclusions

Offspring condition can be affected by the early-life environment and early-life gut microbiome, thus highlighting the importance of understanding how changes in the rearing environment affect individual body mass and survival. Even though our results showed between-individual variation in nestling gut microbiome diversity, we did not find a significant link between brood size and nestling gut microbiome. Moreover, we did not find a significant association between nestling gut microbiome diversity and short-term or mid-term survival. This suggests that other environmental factors (e.g., diet quality) may contribute more to variation in nestling gut microbiome. Further research is needed to uncover the environmental factors that contribute to nestling gut microbiome in wild bird populations, and how gut microbiome may be linked to nestling survival. Gut microbiome can adapt faster to environmental changes than the host, which makes it important to understand the causes of inter-individual variation in microbiome, and how variation in microbiome possibly mediate adaptation to environmental changes.

Availability of data and materials

Sequence files are available at NCBI database (BioProject PRJNA877058). Metadata and R scripts used to run the analyses are available in the Github Repository, https://github.com/marannli/bsm_analyses.

References

  1. Kinross JM, Darzi AW, Nicholson JK. Gut microbiome-host interactions in health and disease. Genome Med. 2011;3(3):1–12.

    Article  Google Scholar 

  2. van Dongen WF, White J, Brandl HB, Moodley Y, Merkling T, Leclaire S, Wagner RH. Age-related differences in the cloacal microbiota of a wild bird species. BMC Ecol. 2013;2013(13):1–12.

    Google Scholar 

  3. David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, Turnbaugh PJ, et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2014;505(7484):559–63.

    Article  CAS  PubMed  Google Scholar 

  4. Sánchez B, Delgado S, Blanco-Míguez A, Lourenço A, Gueimonde M, Margolles A. Probiotics, gut microbiota, and their influence on host health and disease. Mol Nutr Food Res. 2017;61(1):1600240.

    Article  Google Scholar 

  5. Heintz-Buschart A, Wilmes P. Human gut microbiome: function matters. Trends Microbiol. 2018;26(7):563–74.

    Article  CAS  PubMed  Google Scholar 

  6. Egerton S, Culloty S, Whooley J, Stanton C, Ross RP. The gut microbiota of marine fish. Front Microbiol. 2018;9:873.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Pan D, Yu Z. Intestinal microbiome of poultry and its interaction with host and diet. Gut microbes. 2014;5(1):108–19.

    Article  PubMed  Google Scholar 

  8. Hird SM. Evolutionary biology needs wild microbiomes. Front Microbiol. 2017;8:725.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Lozupone CA, Stombaugh JI, Gordon JI, Jansson JK, Knight R. Diversity, stability and resilience of the human gut microbiota. Nature. 2012;489(7415):220–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Zaneveld JR, McMinds R, Vega TR. Stress and stability: applying the Anna Karenina principle to animal microbiomes. Nat Microbiol. 2017;2(9):1–8.

    Article  Google Scholar 

  11. Rosshart SP, Vassallo BG, Angeletti D, Hutchinson DS, Morgan AP, Takeda K, Rehermann B. Wild mouse gut microbiota promotes host fitness and improves disease resistance. Cell. 2017;171(5):1015–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Rosshart SP, Vassallo BG, Angeletti D, Hutchinson DS, Morgan AP, Takeda K, Rehermann B. Wild mouse gut microbiome protects laboratory mice against lethal influenza virus infection and colorectal cancer. J Immunol. 2017;198:68–5.

    Article  Google Scholar 

  13. Worsley SF, Davies CS, Mannarelli ME, Hutchings MI, Komdeur J, Burke T, Richardson DS. Gut microbiome composition, not alpha diversity, is associated with survival in a natural vertebrate population. Anim microbiome. 2021;3(1):1–18.

    Article  Google Scholar 

  14. Davies CS, Worsley SF, Maher KH, Komdeur J, Burke T, Dugdale HL, Richardson DS. Immunogenetic variation shapes the gut microbiome in a natural vertebrate population. Microbiome. 2022;10(1):1–22.

    Article  Google Scholar 

  15. Zilber-Rosenberg I, Rosenberg E. Role of microorganisms in the evolution of animals and plants: the hologenome theory of evolution. FEMS Microbiol Rev. 2008;32:723–35.

    Article  CAS  PubMed  Google Scholar 

  16. Round JL, Mazmanian SK. The gut microbiota shapes intestinal immune responses during health and disease. Nat Rev Immunol. 2009;9:313–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Kopac SM, Klassen JL. Can they make it on their own? Hosts, microbiomes, and the holobiont niche. Front Microbiol. 2016;7:6.

    Article  Google Scholar 

  18. Hansen CHF, Nielsen DS, Kverka M, Zakostelska Z, Klimesova K, Hudcovic T, Hansen AK. Patterns of early gut colonization shape future immune responses of the host. PLoS ONE. 2012;7(3): e34043.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Cox LM, Yamanishi S, Sohn J, Alekseyenko AV, Leung JM, Cho I, Blaser MJ, et al. Altering the intestinal microbiota during a critical developmental window has lasting metabolic consequences. Cell. 2014;158(4):705–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Grond K, Lanctot RB, Jumpponen A, Sandercock BK. Recruitment and establishment of the gut microbiome in arctic shorebirds. FEMS Microbiol Ecol. 2017;93:9.

    Article  Google Scholar 

  21. Knutie SA, Shea LA, Kupselaitis M, Wilkinson CL, Kohl KD, Rohr JR. Early-life diet affects host microbiota and later-life defenses against parasites in frogs. Integr Comp Biol. 2017;57(4):732–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Videvall E, Song SJ, Bensch HM, Strandh M, Engelbrecht A, Serfontein N, Cornwallis CK, et al. Early-life gut dysbiosis linked to juvenile mortality in ostriches. Microbiome. 2020;8(1):1–13.

    Article  Google Scholar 

  23. Carmody RN, Gerber GK, Luevano JM Jr, Gatti DM, Somes L, Svenson KL, Turnbaugh PJ. Diet dominates host genotype in shaping the murine gut microbiota. Cell Host Microbe. 2015;17(1):72–84.

    Article  CAS  PubMed  Google Scholar 

  24. Hooper LV, Littman DR, Macpherson AJ. Interactions between the microbiota and the immune system. Science. 2012;336(6086):1268–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Teyssier A, Matthysen E, Hudin NS, De Neve L, White J, Lens L. Diet contributes to urban-induced alterations in gut microbiota: experimental evidence from a wild passerine. Proc R Soc B. 1920;2020(287):20192182.

    Google Scholar 

  26. Bodawatta KH, Freiberga I, Puzejova K, Sam K, Poulsen M, Jønsson KA. Flexibility and resilience of great tit (Parus major) gut microbiomes to changing diets. Animal microbiome. 2021;3(1):1–14.

    Article  Google Scholar 

  27. Góngora E, Elliott KH, Whyte L. Gut microbiome is affected by inter-sexual and inter-seasonal variation in diet for thick-billed murres (Uria lomvia). Sci Rep. 2021;11(1):1–12.

    Article  Google Scholar 

  28. Schmiedová L, Tomášek O, Pinkasová H, Albrecht T, Kreisinger J. Variation in diet composition and its relation to gut microbiota in a passerine bird. Sci Rep. 2022;12(1):1–13.

    Article  Google Scholar 

  29. Datta FU, Nolana JV, Rowe JB, Gray GD. Protein supplementation improves the performance of parasitised sheep fed a straw-based diet. Int J Parasitol. 1998;28(8):1269–78.

    Article  CAS  PubMed  Google Scholar 

  30. Cotter SC, Simpson SJ, Raubenheimer D, Wilson K. Macronutrient balance mediates trade-offs between immune function and life history traits. Funct Ecol. 2011;25(1):186–98.

    Article  Google Scholar 

  31. Linden M, Gustafsson L, Part T. Selection on fledging mass in the collared flycatcher and the great tit. Ecology. 1992;73(1):336–43.

    Article  Google Scholar 

  32. Spor A, Koren O, Ley R. Unravelling the effects of the environment and host genotype on the gut microbiome. Nat Rev Microbiol. 2011;9(4):279–90.

    Article  CAS  PubMed  Google Scholar 

  33. Bodawatta KH, Hird SM, Grond K, Poulsen M, Jønsson KA. Avian gut microbiomes taking flight. Trends Microbiol. 2022;30:268–80.

    Article  CAS  PubMed  Google Scholar 

  34. Waite DW, Taylor M. Exploring the avian gut microbiota: current trends and future directions. Front Microbiol. 2015;6:673.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO. Development of the human infant intestinal microbiota. PLoS biol. 2007;5(7): e177.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Chen X, Xu J, Ren E, Su Y, Zhu W. Co-occurrence of early gut colonization in neonatal piglets with microbiota in the maternal and surrounding delivery environments. Anaerobe. 2018;49:30–40.

    Article  CAS  PubMed  Google Scholar 

  37. Kalbermatter C, Fernandez Trigo N, Christensen S, Ganal-Vonarburg SC. Maternal microbiota, early life colonization and breast milk drive immune development in the newborn. Front Immunol. 2021;12:1768.

    Article  Google Scholar 

  38. Kohl KD. Diversity and function of the avian gut microbiota. J Comp Physiol B. 2012;182(5):591–602.

    Article  PubMed  Google Scholar 

  39. Trevelline BK, MacLeod KJ, Knutie SA, Langkilde T, Kohl KD. In ovo microbial communities: a potential mechanism for the initial acquisition of gut microbiota among oviparous birds and lizards. Biol Let. 2018;14(7):4.

    Article  Google Scholar 

  40. Ruiz-Rodriguez M, Lucas FS, Heeb P, Soler JJ. Differences in intestinal microbiota between avian brood parasites and their hosts. Biol J Lin Soc. 2009;96(2):406–14.

    Article  Google Scholar 

  41. Benskin CMH, Rhodes G, Pickup RW, Wilson K, Hartley IR. Diversity and temporal stability of bacterial communities in a model passerine bird, the zebra finch. Mol Ecol. 2010;19(24):5531–44.

    Article  PubMed  Google Scholar 

  42. Lee CY, Peralta-Sánchez JM, Martínez-Bueno M, Møller AP, Rabelo-Ruiz M, Zamora-Muñoz C, Soler JJ. The gut microbiota of brood parasite and host nestlings reared within the same environment: disentangling genetic and environmental effects. ISME J. 2020;14(11):2691–702.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Lucas FS, Heeb P. Environmental factors shape cloacal bacterial assemblages in great tit Parus major and blue tit P. caeruleus nestlings. J Avian Biol. 2005;36:510–6.

    Article  Google Scholar 

  44. Hird SM, Carstens BC, Cardiff S, Diffmann DL, Brumfield RT. Sampling locality is more detectable than taxonomy or ecology in the gut microbiota of the brood-parasitic Brown-headed Cowbird (Molothrus ater). PeerJ. 2014;2: e321.

    Article  PubMed  PubMed Central  Google Scholar 

  45. van Veelen HPJ, Salles JF, Tieleman BI. Multi-level comparisons of cloacal, skin, feather and nest-associated microbiota suggest considerable influence of horizontal acquisition on the microbiota assembly of sympatric woodlarks and skylarks. Microbiome. 2017;5(1):1–17.

    Google Scholar 

  46. Devaynes A, Antunes A, Bedford A, Ashton P. Bacterial species richness at three stages of the breeding season in Cyanistes caeruleus (blue tit). Acta Oecologica. 2018;92:123–30.

    Article  Google Scholar 

  47. Ding J, Dai R, Yang L, He C, Xu K, Liu S, et al. Inheritance and establishment of gut microbiota in chickens. Front Microbiol. 2017;8:1967. https://doi.org/10.3389/fmicb.2017.01967.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Preest MR, Folk DG, Beuchat CA. Decomposition of nitrogenous compounds by intestinal bacteria in hummingbirds. Auk. 2003;120(4):1091–101.

    Article  Google Scholar 

  49. Godoy-Vitorino F, Leal SJ, Díaz WA, Rosales J, Goldfarb KC, García-Amado MA, Domínguez-Bello MG, et al. Differences in crop bacterial community structure between hoatzins from different geographical locations. Res Microbiol. 2012;163(3):211–20.

    Article  PubMed  Google Scholar 

  50. Mendoza MLZ, Roggenbuck M, Vargas KM, Hansen LH, Brunak S, Gilbert MTP, Sicheritz-Pontén T. Protective role of the vulture facial skin and gut microbiomes aid adaptation to scavenging. Acta Vet Scand. 2018;60(1):1–19.

    Google Scholar 

  51. Fuirst M, Veit RR, Hahn M, Dheilly N, Thorne LH. Effects of urbanization on the foraging ecology and microbiota of the generalist seabird Larus argentatus. PLoS ONE. 2018;13(12): e0209200.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Teyssier A, Lens L, Matthysen E, White J. Dynamics of gut microbiota diversity during the early development of an avian host: evidence from a cross-foster experiment. Front Microbiol. 2018;9:1524.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Pettifor RA, Perrins CM, McCleery RH. Individual optimization of clutch size in great tits. Nature. 1988;336(6195):160–2.

    Article  Google Scholar 

  54. Wright J, Both C, Cotton PA, Bryant D. Quality vs. quantity: energetic and nutritional trade-offs in parental provisioning strategies. J Anim Ecol 1998;620–634.

  55. Badyaev AV, Ghalambor CK. Evolution of life histories along elevational gradients: trade-off between parental care and fecundity. Ecology. 2001;82(10):2948–60.

    Article  Google Scholar 

  56. Naef-Daenzer B, Keller LF. The foraging performance of great and blue tits (Parus major and P. caeruleus) in relation to caterpillar development, and its consequences for nestling growth and fledging weight. J Anim Ecol. 1999;68:708–18.

    Article  Google Scholar 

  57. de Kogel CH. Long-term effects of brood size manipulation on morphological development and sex-specific mortality of offspring. J Anim Ecol. 1997;66:167–78.

    Article  Google Scholar 

  58. Pettifor RA, Perrins CM, McCleery RH. The individual optimization of fitness: variation in reproductive output, including clutch size, mean nestling mass and offspring recruitment, in manipulated broods of great tits Parus major. J Anim Ecol. 2001;70(1):62–79.

    Google Scholar 

  59. Smith HG, Källander H, Nilsson JA. The trade-off between offspring number and quality in the great tit Parus major. J Anim Ecol. 1989;58:383–401.

    Article  Google Scholar 

  60. Simon K, Verwoolde MB, Zhang J, Smidt H, de Vries RG, Kemp B, Lammers A. Long-term effects of early life microbiota disturbance on adaptive immunity in laying hens. Poult Sci. 2016;95(7):1543–54.

    Article  CAS  PubMed  Google Scholar 

  61. Davidson GL, Somers SE, Wiley N, Johnson CN, Reichert MS, Ross RP, Quinn JL, et al. A time-lagged association between the gut microbiome, nestling weight and nestling survival in wild great tits. J Anim Ecol. 2021;90(4):989–1003.

    Article  PubMed  Google Scholar 

  62. Kabir SM. The role of probiotics in the poultry industry. Int J Mol Sci. 2009;10(8):3531–46.

    Article  CAS  Google Scholar 

  63. Maki JJ, Klima CL, Sylte MJ, Looft T. The microbial pecking order: utilization of intestinal microbiota for poultry health. Microorganisms. 2019;7(10):376.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Du X, Xiang Y, Lou F, Tu P, Zhang X, Hu X, Xiao Y, et al. Microbial community and short-chain fatty acid mapping in the intestinal tract of quail. Animals. 2020;10(6):1006.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Den Besten G, Van Eunen K, Groen AK, Venema K, Reijngoud DJ, Bakker BM. The role of short-chain fatty acids in the interplay between diet, gut microbiota, and host energy metabolism. J Lipid Res. 2013;54(9):2325–40.

    Article  Google Scholar 

  66. Harrison XA, Blount JD, Inger R, Norris DR, Bearhop S. Carry-over effects as drivers of fitness differences in animals. J Anim Ecol. 2011;80(1):4–18.

    Article  PubMed  Google Scholar 

  67. Leonard ML, Horn AG, Gozna A, Ramen S. Brood size and begging intensity in nestling birds. Behav Ecol. 2000;11(2):196–201.

    Article  Google Scholar 

  68. Nicolaus M, Michler SP, Ubels R, Van Der Velde M, Komdeur J, Both C, Tinbergen JM. Sex-specific effects of altered competition on nestling growth and survival: an experimental manipulation of brood size and sex ratio. J Anim Ecol. 2009;78(2):414–26.

    Article  PubMed  Google Scholar 

  69. Perrins CM. Population fluctuations and clutch size in the Great tit, Parus major. J Anim Ecol. 1965;34:601–47.

    Article  Google Scholar 

  70. Smith HG, Källander H, Nilsson JÅ. Effect of experimentally altered brood size on frequency and timing of second clutches in the great tit. Auk. 1987;104(4):700–6.

    Article  Google Scholar 

  71. Hõrak P, Tegelmann L, Ots I, Møller AP. Immune function and survival of great tit nestlings in relation to growth conditions. Oecologia. 1999;121(3):316–22.

    Article  PubMed  Google Scholar 

  72. Sanz JJ, Tinbergen JM. Energy expenditure, nestling age, and brood size: an experimental study of parental behavior in the great tit Parus major. Behav Ecol. 1999;10(5):598–606.

    Article  Google Scholar 

  73. Tinbergen JM, Verhulst S. A fixed energetic ceiling to parental effort in the great tit? J Anim Ecol. 2000;69(2):323–34.

    Article  Google Scholar 

  74. Neuenschwander S, Brinkhof MW, Kölliker M, Richner H. Brood size, sibling competition, and the cost of begging in great tits (Parus major). Behav Ecol. 2003;14(4):457–62.

    Article  Google Scholar 

  75. Tinbergen JM, Sanz JJ. Strong evidence for selection for larger brood size in a great tit population. Behav Ecol. 2004;15(4):525–33.

    Article  Google Scholar 

  76. Cossin-Sevrin N, Hsu BY, Marciau C, Viblanc VA, Ruuskanen S, Stier A. Effect of prenatal glucocorticoids and thyroid hormones on developmental plasticity of mitochondrial aerobic metabolism, growth and survival: an experimental test in wild great tits. J Exp Biol. 2022;225(9): jeb243414.

    Article  PubMed  Google Scholar 

  77. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, Knight R, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6(8):1621–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Parada AE, Needham DM, Fuhrman JA. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ Microbiol. 2016;18(5):1403–14.

    Article  CAS  PubMed  Google Scholar 

  79. de Muinck EJ, Trosvik P, Gilfillan GD, Hov JR, Sundaram AY. A novel ultra high-throughput 16S rRNA gene amplicon sequencing library preparation method for the Illumina HiSeq platform. Microbiome. 2017;5(1):1–15.

    Article  Google Scholar 

  80. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet Journal. 2011;17(1):10–2.

    Article  Google Scholar 

  81. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Glöckner FO, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2012;41(D1):D590–6.

    Article  PubMed  PubMed Central  Google Scholar 

  83. Davis NM, Proctor DM, Holmes SP, Relman DA, Callahan BJ. Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data. Microbiome. 2018;2018(6):1–14.

    Google Scholar 

  84. McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8(4): e61217.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Lahti L & Sudarshan Shetty et al. Tools for microbiome analysis in R. 2017; http://microbiome.github.com/microbiome

  86. Wickham H. Data analysis. In ggplot2 (pp. 189–201). 2016; Springer, Cham.

  87. Wright ES. Using DECIPHER v2. 0 to analyze big biological sequence data in R. R Journal. 2016; 8:1.

  88. Schliep KP. phangorn: phylogenetic analysis in R. Bioinformatics. 2011;27(4):592–3.

    Article  CAS  PubMed  Google Scholar 

  89. Paradis E, Blomberg S, Bolker B, Brown J, Claude J, Cuong HS, Didier G et al. Package ‘ape’. Analyses of phylogenetics and evolution, version. 2019; 2(4).

  90. Yu G, Smith DK, Zhu H, Guan Y, Lam TTY. ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol. 2017;8(1):28–36.

    Article  Google Scholar 

  91. Bates D, Maechler M, Bolker B, Walker S. Fitting Linear Mixed-Effects Models Using lme4. J Stat Softw. 2015;67(1):1–48. https://doi.org/10.18637/jss.v067.i01.

    Article  Google Scholar 

  92. Monrós JS, Belda EJ, Barba E. Post-fledging survival of individual great tits: the effect of hatching date and fledging mass. Oikos. 2002;99(3):481–8.

    Article  Google Scholar 

  93. Fox J, Weisberg S, Adler D, Bates D, Baud-Bovy G, Ellison S, Monette G, et al. Package ‘car.’ Vienna: R Foundation for Statistical Computing; 2012. p. 16.

    Google Scholar 

  94. Hartig F & Hartig MF. Package ‘DHARMa’. R package. 2017.

  95. Bray JR, Curtis JT. An ordination of the upland forest communities of southern Wisconsin. Ecol Monogr. 1957;27(4):326–49.

    Article  Google Scholar 

  96. Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71(12):8228–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  97. Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’hara RB, Oksanen MJ etl al. Package ‘vegan’. Community ecology package, version. 2013; 2(9):1–295.

  98. Moreno J, Cowie RJ, Sanz JJ, Williams RS. Differential response by males and females to brood manipulations in the pied flycatcher: energy expenditure and nestling diet. J Anim Ecol. 1995;64:721–32.

    Article  Google Scholar 

  99. Askenmo C. Effects of addition and removal of nestlings on nestling weight, nestling survival, and female weight loss in the Pied Flycatcher Ficedula hypoleuca (Pallas). Ornis Scand. 1977;8:1–8.

    Article  Google Scholar 

  100. Nur N. The consequences of brood size for breeding blue tits II. Nestling weight, offspring survival and optimal brood size. J Anim Ecol. 1984;53:497–517.

    Article  Google Scholar 

  101. Sanz JJ. Clutch size manipulation in the pied flycatcher: effects on nestling growth, parental care and moult. J Avian Biol. 1998;28:157–62.

    Article  Google Scholar 

  102. Burness GP, McClelland GB, Wardrop SL, Hochachka PW. Effect of brood size manipulation on offspring physiology: an experiment with passerine birds. J Exp Biol. 2000;203(22):3513–20.

    Article  CAS  PubMed  Google Scholar 

  103. Naef-Daenzer L, Naef-Daenzer B, Nager RG. Prey selection and foraging performance of breeding Great Tits Parus major in relation to food availability. J Avian Biol. 2000;31(2):206–14.

    Article  Google Scholar 

  104. Rytkönen S, Orell M. Great tits, Parus major, lay too many eggs:experimental evidence in mid-boreal habitats. Oikos. 2001;93:439–50.

    Article  Google Scholar 

  105. van der Most PJ, de Jong B, Parmentier HK, Verhulst S. Trade-off between growth and immune function: a meta-analysis of selection experiments. Funct Ecol. 2011;25(1):74–80.

    Article  Google Scholar 

  106. Rauw WM. Immune response from a resource allocation perspective. Front Genet. 2012;3:267.

    Article  PubMed  PubMed Central  Google Scholar 

  107. Naef-Daenzer B, Keller LF. The foraging performance of great and blue tits (Parus major and P. caeruleus) in relation to caterpillar development, and its consequences for nestling growth and fledging weight. J Anim Ecol. 1999;6(84):708–18.

    Article  Google Scholar 

  108. Barba E, Gil-Delgado JA. Seasonal variation in nestling diet of the great tit Parus major in orange groves in eastern Spain. Ornis Scand. 1990;21:296–8.

    Article  Google Scholar 

  109. Tinbergen JM, Dietz MW. Parental energy expenditure during brood rearing in the great tit (Parus major) in relation to body mass, temperature, food availability and clutch size. Funct Ecol. 1994;8:563–72.

    Article  Google Scholar 

  110. Ilmonen P, Hasselquist D, Langefors Å, Wiehn J. Stress, immunocompetence and leukocyte profiles of pied flycatchers in relation to brood size manipulation. Oecologia. 2003;136(1):148–54.

    Article  PubMed  Google Scholar 

  111. Kohl KD, Brun A, Caviedes-Vidal E, Karasov WH. Age-related changes in the gut microbiota of wild House Sparrow nestlings. Ibis. 2019;161(1):184–91.

    Article  Google Scholar 

  112. Waite DW, Taylor MW. Characterizing the avian gut microbiota: membership, driving influences, and potential function. Front Microbiol. 2014;5:223.

    Article  PubMed  PubMed Central  Google Scholar 

  113. Davidson GL, Wiley N, Cooke AC, Johnson CN, Fouhy F, Reichert MS, Quinn JL, et al. Diet induces parallel changes to the gut microbiota and problemsolving performance in a wild bird. Sci Rep. 2020;10(1):1–13.

    Article  Google Scholar 

  114. Murray MH, Lankau EW, Kidd AD, Welch CN, Ellison T, Adams HC, Hernandez SM, et al. Gut microbiome shifts with urbanization and potentially facilitates a zoonotic pathogen in a wading bird. PLoS ONE. 2020;15(3): e0220926.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  115. Bodawatta KH, Koane B, Maiah G, Sam K, Poulsen M, Jønsson KA. Species-specific but not phylosymbiotic gut microbiomes of New Guinean passerine birds are shaped by diet and flight-associated gut modifications. Proc R Soc B. 2021;288(1949):20210446.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  116. Bodawatta KH, Klečková I, Klečka J, Pužejová K, Koane B, Poulsen M, Sam K, et al. Specific gut bacterial responses to natural diets of tropical birds. Sci Rep. 2022;12(1):1–15.

    Article  Google Scholar 

  117. Alberdi A, Aizpurua O, Bohmann K, Zepeda-Mendoza ML, Gilbert MTP. Do vertebrate gut metagenomes confer rapid ecological adaptation? Trends Ecol Evol. 2016;31(9):689–99.

    Article  PubMed  Google Scholar 

  118. Wilkin TA, King LE, Sheldon BC. Habitat quality, nestling diet, and provisioning behaviour in great tits Parus major. J Avian Biol. 2009;40(2):135–45.

    Article  Google Scholar 

  119. Rytkönen S, Vesterinen EJ, Westerduin C, Leviäkangas T, Vatka E, Mutanen M, Orell M, et al. From feces to data: a metabarcoding method for analyzing consumed and available prey in a bird-insect food web. Ecol Evol. 2019;9(1):631–9.

    Article  PubMed  Google Scholar 

  120. Drobniak SM, Cichoń M, Janas K, Barczyk J, Gustafsson L, Zagalska-Neubauer M. Habitat shapes diversity of gut microbiomes in a wild population of blue tits Cyanistes caeruleus. J Avian Biol. 2022;2: e02829.

    Google Scholar 

  121. Fleischer R, Risely A, Hoeck PE, Keller LF, Sommer S. Mechanisms governing avian phylosymbiosis: genetic dissimilarity based on neutral and MHC regions exhibits little relationship with gut microbiome distributions of Galápagos mockingbirds. Ecol Evol. 2020;10(23):13345–54.

    Article  PubMed  PubMed Central  Google Scholar 

  122. Kloos WE. Natural populations of the genus Staphylococcus. Annu Rev Microbiol. 1980;34(1):559–92.

    Article  CAS  PubMed  Google Scholar 

  123. Gharibzahedi SMT, Razavi SH, Mousavi SM. Characterization of bacteria of the genus Dietzia: an updated review. Ann Microbiol. 2014;64(1):1–11.

    Article  CAS  Google Scholar 

  124. Hanning I, Diaz-Sanchez S. The functionality of the gastrointestinal microbiome in non-human animals. Microbiome. 2015;3:5. https://doi.org/10.1186/s40168-015-0113-6.

    Article  Google Scholar 

  125. Doak DF, Bigger D, Harding EK, Marvier MA, O’malley RE, Thomson D. The statistical inevitability of stability-diversity relationships in community ecology. Am Nat. 1998;151(3):264–76.

    Article  CAS  PubMed  Google Scholar 

  126. Naef-Daenzer B, Widmer F, Nuber M. Differential post-fledging survival of great and coal tits in relation to their condition and fledging date. J Anim Ecol. 2001;70(5):730–8.

    Article  Google Scholar 

  127. Nilsson JA, Svensson M. Sibling competition affects nestling growth strategies in marsh tits. J Anim Ecol. 1996;65:825–36.

    Article  Google Scholar 

  128. Magrath RD. Nestling weight and juvenile survival in the blackbird, Turdus merula. J Anim Ecol. 1991;60:335–51.

    Article  Google Scholar 

  129. Ringsby TH, Sæther BE, Solberg EJ. Factors affecting juvenile survival in house sparrow Passer domesticus. J Avian Biol. 1998;29:241–7.

    Article  Google Scholar 

  130. Both C, Visser ME, Verboven N. Density–dependent recruitment rates in great tits: the importance of being heavier. Proc R Soc Lond Seri B Biol Sci. 1999;266(1418):465–9.

    Article  Google Scholar 

  131. Cleasby IR, Nakagawa S, Gillespie DO, Burke T. The influence of sex and body size on nestling survival and recruitment in the house sparrow. Biol J Lin Soc. 2010;101(3):680–8.

    Article  Google Scholar 

  132. Tinbergen JM, Boerlijst MC. Nestling weight and survival in individual great tits (Parus major). J Anim Ecol. 1990;59:1113–27.

    Article  Google Scholar 

Download references

Acknowledgements

We would like to thank Charli Davies, Petri Papponen, Sari Viinikainen and Elina Virtanen at the University of Jyväskylä for their help with bioinformatics (CD) laboratory optimizations (PP, SV and EV), as well as Jorma Nurmi, Robin Cristofari, Natacha Garcin, Toni Laaksonen at the University of Turku and volunteer bird ringers for their help in the field.

Funding

Open Access funding provided by University of Jyväskylä (JYU). This study was funded by the Emil Aaltonen Foundation. NCS was supported by EDUFI Fellowship (Opetushallitus). AS was funded by the Turku Collegium for Science and Medicine, who contributed to fund the field study. AS acknowledges funding from the European Commission Marie Sklodowska-Curie Postdoctoral Fellowship (#894963) at the time of writing.

Author information

Authors and Affiliations

Authors

Contributions

The idea for this study was by AS, SR and NCS. Sample collection was done by NCS, MH, AS and SR. Laboratory analyses were done by ML with assistance from SR, KG and EV. Sequence processing and statistical analyses were done by ML with assistance from KG, EV (bioinformatics) and SR (statistics). Manuscript was written by ML and all authors commented and approved the manuscript.

Corresponding author

Correspondence to Martta Liukkonen.

Ethics declarations

Ethics approval and consent to participate

All animal work was conducted under relevant national and international guidelines and legislation. The animal work was licensed by the environmental and ethical committee of Varsinais-Suomi (environmental permit license number VARELY/890/2020; animal ethics permit number ESAVI/5454/2020). The birds were ringed from their nest-sites by permission from the Finnish Ringing Centre to SR (Finnish Museum of Natural History).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Brood size before and after manipulation: brood sizes between treatment groups were tested with a linear model to see if the differences were statistically significant.

Additional file 2:

Rarefaction curves for the unrarefied dataset. Species (ASVs) plateaued at about 5000 reads which was used as the rarefying depth.

Additional file 3:

Phylogenetic tree using the Newick-format. The tree describes the dissimilarity among the treatment groups. Each tip represents an individual sample, and each tip is colored and shaped based on treatment. Treatment groups are clustered using the UPGMA algorithm.

Additional file 4:

(A) Linear mixed effects model for gut microbiome diversity (Shannon Diversity Index and Chao1 Richness) and brood size manipulation treatment. (B) Linear mixed effects model for GM diversity (Shannon Diversity Index and Chao1 Richness) and manipulated brood size.

Additional file 5:

A linear mixed effects model investigating the effects of brood size manipulation on nestling body mass on day 7 and day 14 post-hatch.

Additional file 6:

A linear mixed effects model investigating the effects of manipulated brood size on nestling body mass on day 7 and day 14 post-hatch.

Additional file 7:

A linear mixed effects model investigating the associations between alpha diversity (Shannon Diversity Index and Chao1 Richness) and brood size manipulation.

Additional file 8:

A linear mixed effects model investigating the association between alpha diversity (Shannon Diversity Index and Chao1 Richness) and manipulated brood size.

Additional file 9:

A generalized linear model exploration into alpha diversity’s (Shannon Diversity Index and Chao1 Richness) association with short-term (survival to fledging) and mid-term (apparent juvenile) survival.

Additional file 10:

The gut microbiome alpha diversity (Shannon Diversity Index and Chao1 Richness) and short-term survival.

Additional file 11:

Generalized linear model to measure the association between alpha diversity (Shannon Diversity Index and Chao1 Richness) survival to fledging and apparent juvenile survival.

Additional file 12:

Ordination of the gut microbial communities.

Additional file 13:

Differential analysis of abundance (DESeq2) to assess the ASV abundance between the treatment groups

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liukkonen, M., Hukkanen, M., Cossin-Sevrin, N. et al. No evidence for associations between brood size, gut microbiome diversity and survival in great tit (Parus major) nestlings. anim microbiome 5, 19 (2023). https://doi.org/10.1186/s42523-023-00241-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s42523-023-00241-z

Keywords