Skip to main content

Early-life factors shaping the gut microbiota of Common buzzard nestlings

Abstract

Background

Exploring the dynamics of gut microbiome colonisation during early-life stages is important for understanding the potential impact of microbes on host development and fitness. Evidence from model organisms suggests a crucial early-life phase when shifts in gut microbiota can lead to immune dysregulation and reduced host condition. However, our understanding of gut microbiota colonisation in long-lived vertebrates, especially during early development, remains limited. We therefore used a wild population of common buzzard nestlings (Buteo buteo) to investigate connections between the early-life gut microbiota colonisation, environmental and host factors.

Results

We targeted both bacterial and eukaryotic microbiota using the 16S and 28S rRNA genes. We sampled the individuals during early developmental stages in a longitudinal design. Our data revealed that age significantly affected microbial diversity and composition. Nest environment was a notable predictor of microbiota composition, with particularly eukaryotic communities differing between habitats occupied by the hosts. Nestling condition and infection with the blood parasite Leucocytozoon predicted microbial community composition.

Conclusion

Our findings emphasise the importance of studying microbiome dynamics to capture changes occurring during ontogeny. They highlight the role of microbial communities in reflecting host health and the importance of the nest environment for the developing nestling microbiome. Overall, this study contributes to understanding the complex interplay between microbial communities, host factors, and environmental variables, and sheds light on the ecological processes governing gut microbial colonisation during early-life stages.

Introduction

The gut microbiome constitutes a diverse assemblage of microorganisms encompassing bacteria, archaea, viruses, and microbial eukaryotes, which collectively exert significant influences on various host processes [59, 94]. Numerous studies have demonstrated the importance of gut microbial communities in regulating physiological functions such as digestion, absorption, metabolism, and immune response, fostering host health across diverse animal species [105].

A fundamental aspect of studying animal microbiota communities involves understanding how they are acquired and how these communities evolve over time. A comprehensive study of humans conducted by Martino et al [74] demonstrated a sequence of transformations in gut microbial diversity and composition throughout different life stages. The process commences with initial colonisation by pioneer species at birth, followed by the emergence of body site-specific microbial communities (”primary succession”). Over time, these communities grow in complexity until they reach a stable structure (”secondary succession”). Ultimately, as the host organism senesces, the community undergoes a final transformation marked by a decline in diversity (”third succession”).

In avian species, the chick gastrointestinal tract is initially colonised by various transient species, with bacterial communities gradually transitioning into a more stable adult state [37]. Avian embryos develop within a closed and essentially sterile environment (the egg); gut microbial colonisation commences shortly after birth and is influenced by interactions among early-life environment, parental microbial transmission, and diet [26, 42, 54]. This early establishment phase, aptly captured by the “nidobiome” concept, encompasses the collective impacts of parents, the nest, and nestlings on the initial microbiome assembly [24, 26, 36]. This concept considers their interactions over time, as well as the inherent feedback mechanisms between hosts, their microbiome, and the surrounding environment [24]. Disruptions of these early-life microbiome assembly processes can adversely affect the host [55, 99] by influencing immune system development [41, 54], and diminishing resistance to parasitic infections [7, 78]. Understanding the factors influencing avian gut microbial composition in early life is therefore of great importance due to the potential long-term impacts of the gut microbiota on host development and body condition [55, 77].

The relationship between the gut microbiome and host condition is complex, exhibiting variation across different host taxa [42, 52, 55]. While many studies have shown that the gut microbiome undergoes dynamic changes throughout development, important questions persist regarding its predictive role in host health [105, 106]. Specifically, little is known about the developmental stages at which the microbiome may influence host condition [33, 114]. In general, higher gut microbial diversity has been associated with improved host health [95, 101] whereas low diversity has been linked to disease states [111]. A diverse gut microbiota exhibits a broader spectrum of metabolic capabilities, enhances immune signalling and competition against harmful pathogens [17, 29, 79]. It is also possible that high diversity might be related to a state of dysbiosis and thus a reduction in diversity might be indicative of a return to a state of homeostasis [30, 49]. On the other hand, emerging evidence suggests that the composition of the gut microbiota, rather than diversity alone, may have an even greater impact on host health [118].

In addition to host and environmental factors, parasite infection can trigger changes in the host microbiome [67, 70, 112]. For example, Madlala et al [67] showed that Eimeria spp. infection disrupts the gut environment increasing susceptibility to diseases that significantly threaten the health and productivity of chickens. On the other hand, Plasmodium infection appears not to markedly alter microbial diversity and composition in canaries (Serinus canaria domestica), Hawaiian honeycreepers (subfamily Hemithraupinae) and Eurasian tree sparrows (Passer montanus), although correlations with specific bacterial taxa were identified [7, 78, 91].

Despite this body of knowledge, the majority of gut microbiome research focuses exclusively on bacteria, the dominant taxonomic group in vertebrate gut ecosystems [47]. Relatively few studies have investigated the structure and diversity of other microbial kingdoms within the gut microbiome, such as viruses, archaea and microbial eukaryotes, and the underlying drivers of such variation [12, 104, 119]. Nonetheless, mounting evidence suggests that host-associated eukaryotic communities, particularly fungal communities, contribute to host health by participating in vital host processes. For example, within the gut, eukaryotic species play pivotal roles in digestion [120] and the development and modulation of the host immune system [122]. While most of the evidence supporting the importance of gut eukaryotic communities on host health comes from captive and model organisms [122], scarce research has explored patterns of eukaryotic variation in wild populations [104, 116, 119].

Here, we aim to elucidate the factors that shape the early-life gut microbiota of common buzzard (Buteo buteo) nestlings and the effects of the gut microbiota on body condition, focusing on both bacterial and eukaryotic components of the microbiome. On average, breeding pairs of common buzzards typically have brood sizes of approximately two nestlings [57] and there is a dominance hierarchy within the brood, with older nestlings usually being stronger at competing for food resources [39]. During their approximately five-week nestling period (\(\approx\) 35 days), common buzzard nestlings rely entirely on parental care and provisioning, followed by a post-nestling phase where fledglings remain in the nest vicinity, receiving parental care and being provisioned even after they are fully fledged (\(\approx\) 47 days) [110]. We collected data at two distinct time points during the nestling period of each nestling. Our dataset includes comprehensive information on environmental factors, such as nest sharing, habitat type, and sampling year, as well as key host factors including age, sex, rank (dominance hierarchy within the brood), blood parasite infection status, weight, and wing length. Additionally, we collected microbiome samples via cloacal swabs to thoroughly assess the microbial diversity and composition. We adopt a multimarker approach utilising both the conventional 16S rRNA gene and the less explored 28S rRNA gene. This strategy aims not only to identify bacteria but also to capture diverse eukaryotic communities beyond fungal populations [56], thereby contributing to a deeper understanding of the broader eukaryotic microbiome.

We hypothesise that factors such as nest sharing, age, and body condition, play pivotal roles in shaping the diversity and composition of the gut microbiome in common buzzard nestlings. Concerning microbiome diversity, two plausible scenarios emerge: 1. An increase in diversity with increasing age, suggesting that maturing individuals acquire an increasingly diverse microbiome, enabling them to take advantage of diverse functional capabilities provided by a broader array of microorganisms [17]; 2. As individuals mature, diversity declines, indicating a gradual transition to a more stable community, after a rapid and uncontrolled colonisation which occurred shortly after hatching [30]. By unravelling the complex relationships among microbial communities, host and environmental factors, we aim to shed light on the mechanisms driving gut microbial colonisation in the first stages of life.

Methods

Study area and sample collection

This study was conducted within a designated 300 km\(^{2}\) study area located in North Rhine-Westphalia, Germany (8\(^\circ\) 25′ E and 52\(^\circ\) 06′ N). For a detailed description of the study system see Krüger [57] and Jonker et al [50]. A total of 117 common buzzard nestlings were sampled, with 107 individuals being sampled in 2020 and 10 individuals being sampled in 2021. Out of these, 5 individuals were sampled once, while the majority were sampled twice, and one individual was sampled three times. Individuals were sampled from 54 nests with a brood size per nest of approximately 2.17 ± 0.853 nestlings. The sampling was distributed across three different habitats: North, South, and Teuto (Teutoburg Forest), as described in Krüger [57] and Jonker et al [50] (Appendix A Table S1).

Morphometric measurements were recorded during each sampling event, including body weight (to the nearest 5 g) and wing length (to the nearest mm). Individuals were re-sampled on average nine days after the initial sampling (mean ± s.d. = 8.7 ± 5.23 days). The average age at the first sampling point was 19 days (mean ± s.d. = 19.3 ± 5.29 days), while the average age at the second sampling point was 28 days (mean ± s.d. = 27.8 ± 5.16 days). Due to the inherent difficulty of precisely determining nestling ages without the daily disturbance of nests, it was not possible to sample individuals at the exact same ages.

Nestling age estimation (after first sampling) was done by employing a sex-specific polynomial regression model on wing length. This model was constructed based on growth curve data published by Bijlmsa [13]. Body condition index (BCI) was calculated by obtaining the residuals from a logarithmic regression of weight on wing length, while accounting for sex.

Blood samples (500–1000 \(\upmu\)l) were collected from the ulnar vein. A small drop of blood was used to prepare blood smears and the remaining volume was placed into 96% ethanol tubes and stored at −20 \(^\circ\)C. Sex determination followed the standard protocol described by Fridolfsson and Ellegren [40], while infection status with the haemosporidian blood parasite Leucocytozoon was determined using the procedure outlined in Chakarov et al [25].

For gut microbiota analysis, adult birds were swabbed in the field by inserting the entire head of the swab into the cloaca and swabbing the inside for approximately 5 s in circular motions. The swab was then transferred to a tube containing RNAlater (Sigma- Aldrich, R0901), which was stored in dry ice and subsequently transferred to long-term storage at −80 \(^\circ\)C. Environmental controls were obtained by swabbing various surfaces around the working station. Additionally, blank swabs were collected as negative controls to account for potential contamination during the sampling process.

DNA extraction from cloacal swabs

DNA extraction from cloacal swabs was performed using a modified phenol-chloroform protocol. Initially, the swabs were centrifuged (10 min at 13000 rpm) and the RNAlater was removed. Subsequently, the swabs were re-suspended in an extraction buffer and subjected to mechanical lysis using 3 mm stainless steel beads (15 min, 50 Hz; TissueLyser LT; Qiagen, Hilden, Germany). Following centrifugation, proteinase K was added, and the samples were incubated at 56 \(^\circ\)C for 2 h. We then used a modified phenol-chloroform procedure to purify the DNA. Adjustments in reagent amounts and centrifugation times were made to accommodate the low biomass of the samples (see Appendix B).

Library preparation

For gene library preparation the ”Illumina 16S Metagenomic Library Preparation Guide” (15044223 Rev.B) was followed. A multi-marker approach was employed, targeting the V4 region of the 16S ribosomal RNA (rRNA) using primers 515F (Parada) [81] and 806R (Apprill) [3]. Additionally, to enhance taxonomic coverage for eukaryotes and to prevent bacterial read contamination [56], the D8-D9 region of the 28S rRNA was targeted using the primers GA20F[6]/RM9Rb[66]. This approach allows for a comprehensive assessment of the eukaryotic microbiota, encompassing not only fungal communities, as with the classic ITS2 (Internal Transcribed Spacer 2), but also other microbial eukaryotes typically targeted by 18S rRNA [56, 75, 120].

Polymerase chain reactions (PCR) were performed were performed separately for each marker in a 25 \(\upmu\)l reaction volume, containing 5 \(\upmu\)l DNA, 12.5 \(\upmu\)l KAPA HiFi HotStart ReadyMix (KAPA Biosystems, MA, United States), 1 \(\upmu\)l of each primer (1 \(\upmu\)M), and 6 \(\upmu\)l of PCR-grade water. The PCR amplification conditions comprised an initial denaturation step at 95 \(^\circ\)C for 3 min, followed by 30 cycles of denaturation at 95 \(^\circ\)C for 30 s, annealing at 55 \(^\circ\)C for 16 S rRNA and 60 \(^\circ\)C for 28 S rRNA for 25 s, extension at 72 \(^\circ\)C for 40 s, and a final extension step of 5 min at 72 \(^\circ\)C. The presence and size of amplicons were verified using a 2% agarose gel. Subsequently, PCR products were purified using the Agencourt AMpure XP PCR purification system (Beckman Coulter, Brea, CA, United States), following the manufacturer’s instructions. To increase amplicon concentration, a second PCR was performed using the purified PCR products from the first PCR (PCR-1). The same PCR conditions were applied, except for reducing the annealing time to 20 s. The PCR-1 and PCR-2 products were pooled and subjected to another round of purification using the Agencourt AMpure XP PCR purification system.

The purified PCR products were sent to the CeBiTec sequencing centre at Bielefeld University for subsequent processing. Index-PCR was performed using Illumina Nextera XT V2 index kits. Quantification of sequencing libraries was conducted using the Fragment Analyzer (Agilent). Following quantification, the libraries were diluted and equimolarly pooled. The pooled amplicon libraries were then sequenced on a Illumina MiSeq platform (0.4% MiSeq run), applying the protocol for 2\(\times\)300 bp paired-end reads. In addition to the 235 biological samples, the final library pool comprised of two replicates of the \(ZymoBIOMICS^{TM}\) Microbial Community Standard (D6300, Irvine, CA, USA), two \(ZymoBIOMICS^{TM}\) Microbial Community DNA Standard (D6305), four environmental controls, five extraction blanks, and four PCR negatives. Negative controls were included to monitor potential contamination throughout the entire procedure, while positive controls served for posterior workflow validation.

16S rRNA sequence data processing

Demultiplexed Illumina sequence data were imported into QIIME2 (Quantitative Insights Into Microbial Ecology 2, version 2022.11 [16]. Quality assessment of the reads was performed by visualising quality plots. To filter out low-quality bases and infer Amplicon Sequencing Variants (ASVs), the Divisive Amplicon Denoising Algorithm pipeline (DADA2) was employed [22]. The forward and reverse sequences were truncated at 253 and 185 base pairs, respectively, with 20 base pairs trimmed from the 5’ end of the reads. Taxonomy was assigned to the ASVs using a naive Bayes taxonomic classifier trained on the SILVA SSU 138.1 database [86]. The classifier was built and trained using the REference Sequence annotation and CuRatIon Pipeline plugin (RESCRIPt) [89].

The data were imported into R version 4.2.2 [87] using the qiime2R package version 0.99.6 [15]. Negative controls (environmental, extraction, and PCR negatives) were included in the analysis. 578 sequence contaminants were identified and removed using the decontam package version 1.18 [35]. The ”prevalence” method, with a probability threshold of 0.1, was applied for contaminant removal. Data was imported back into QIIME2 and ASVs assigned to Mitochondria, Chloroplast, Vertebrata, Eukaryota, and unassigned taxa were filtered out using QIIME2. Only taxa present in more than one sample and samples containing a minimum frequency of 500 reads were retained for further analysis. The remaining ASVs were aligned using MAFFT [51] as implemented in the q2- alignment plugin. The aligned sequences were then used to construct a phylogeny with FastTree [85], as implemented in the q2-phylogeny plugin. To assess the performance of the pipeline, the q2-quality-control plugin was used. This allowed us to evaluate the accuracy with which the expected taxonomic composition, derived from the community standards, was reconstructed (The results of the microbial community standards analysis are presented in Appendix C). To visualise microbial community composition, taxa-bar plots were generated exclusively for the core taxa, defined as taxa present in at least 70% of the samples. (Detailed scripts and intermediate results can be found in Appendix D.)

28S rRNA sequence data processing

Demultiplexed Illumina sequence data were imported into R version 4.2.2. The data underwent initial processing using Cutadapt version 4.4 [72] to identify and remove locus-specific primers from both reads. Quality assessment was performed through visualising quality plots.

ASV inference was conducted using DADA2, based on the pipeline outlined in Callahan et al [23]. After quality assessment, the sequences were trimmed to eliminate low quality regions: forward reads were truncated to 230 base pairs, and reverse reads to 185 base pairs. Short sequences below 100 base pairs were removed. Paired-end sequences were merged. A high percentage of reads failed to merge due to the larger fragment resulting in a low base-pair overlap. In order to maximise data utilisation, the unmerged paired-end sequences were concatenated. Chimeras were subsequently removed, and the data were imported into QIIME2. A primer-region specific classifier was built and trained using RESCRIPt. Taxonomy was assigned using the naive Bayes taxonomic classifier trained on the SILVA LSU 138.1 database. The processed data was then imported back into R. The decontam pipeline was applied, 218 ASV contaminants were removed using the ”prevalence” method with a probability threshold of 0.1.

Similar to the 16 S rRNA analysis, taxonomy-based filtering (host reads, Mitochondria, Chloroplast, Vertebra and unassigned reads were removed), removal of unique features, and filtering of samples with fewer than 500 reads were performed in QIIME2. ASVs were aligned using MAFFT as implemented in the q2-alignment plugin, and phylogeny was constructed using FastTree through the q2-phylogeny plugin.

The performance of the pipeline was evaluated using the q2-quality-control plugin (Appendix C) and taxa-bar plots (core taxa only) were generated to visualise microbial community composition. (Detailed scripts and intermediate results can be found in Appendix D.)

Alpha diversity

In Qiime2, rarefaction curves were generated using the q2-diversity-alpha-rarefaction plugin to assess sequencing depth and sample coverage. After inspection, the 16S rRNA dataset was rarefied at 4000 reads, while the 28S rRNA dataset was rarefied at 2000 reads (Appendix A Fig. S1–S3). Two measures of alpha diversity, the Shannon diversity index [97] and Faith’s Phylogenetic Diversity (Faith PD) [38], were calculated using the q2-diversity-alpha plugin. Correlations among the variables studied were assessed using the Pearson correlation coefficient as part of the correlation package version 0.8.4 in R [69]. To investigate individual differences in the gut microbiome, linear mixed models (LMMs) with a Gaussian distribution and identity link were performed using the lmer function from the lme4 package [9] in R. The significance of factors included in the models was tested using analysis of variance (ANOVA). This analysis aimed to examine the relationship between gut microbiota diversity, age and body condition index (BCI), while accounting for differences in habitat, sex, year, hatching sequence (rank), and Leucocytozoon infection. To incorporate the longitudinal nature of the dataset and the fact that individuals belonged to the same/different nests, Nest ID and Individual ID were included as nested random effects (Individual ID nested within Nest ID) as follows:

$$\begin{aligned}{} & {} \text {Microbiota diversity }\sim Age + BCI + Rank + Habitat + Sex + Year \\{} & {} \quad + Infection + (1|Nest ID/Individual ID) \end{aligned}$$

The significance of the fixed effects was assessed at \(\alpha\) < 0.05. Significance of the random effects was tested using the ‘ranova” function from the lmerTest package [58]. Marginal and conditional R\(^2\) values were calculated using the MuMIn package [8]. Assumptions of normality and homogeneity of variance of the residuals were assessed through visual inspection of plots using the performance package [65] and furthermore analysed with the Shapiro-Wilk test. When necessary, the data were normality transformed. (For complete scripts and intermediate results, refer to Appendix E)

Beta diversity

Microbiota community analysis was performed using R version 4.2.2. The unrarified dataset, excluding samples with fewer than 4000 reads for 16S rRNA and 2000 reads for 28S rRNA, was subjected to Cumulative Sum Scaling (CSS) normalisation [83] using the metagenomeSeq package version 1.30.0 [82] to account for uneven sequence coverage. Bray-Curtis dissimilarities (BC) [19] and weighted UniFrac distances (WU) [64] were then computed using the phyloseq package [76]. Our dataset violated some of the underlying assumptions of PERMANOVA, namely the exchangeability assumption under the null hypothesis (independence of variables) and the challenge of assessing homogeneity of variances with continuous variables [2]. We therefore adopted a Bayesian framework to model pairwise (dyadic) values, described by Raulo et al [88]. The analysis pipeline described at Analysing-dyadic-data-with-brms was implemented. Bayesian regression models were fitted using the brms package [21]. brms allows for the incorporation of random effect structures to account for dependence in dyadic data and repeated sampling [20]. The models included all pairwise sample comparisons except within-sample comparisons. BC and WU were used as response variables; matrices of nest sharing, habitat, sex, and year (all coded as 0/1 for different/same), age differences (\(\Delta\)Age), BCI differences (\(\Delta\)BCI) and infection status combination (II: infected-infected; INi: infected/non-infected; NiNi: non-infected/non-infected) were included as fixed effects. Continuous variables were either naturally scaled or transformed to range from 0 to 1. To control for data dependency resulting from pairwise comparisons and repeat samples per individual, two multi-membership random effects were included in the model: one capturing the individuals in each dyad (ID A + ID B) and another capturing the samples in each dyad (Sample A + Sample B):

$$\begin{aligned}{} & {} \text {Microbiota composition }\sim \Delta Age + \Delta BCI + Habitat(similarity) + Sex(similarity) \\{} & {} \quad + Nest(sharing) + Year(similarity) \\{} & {} \quad + Infection(combination) + (1|mm(sampleA,sampleB)) +(1|mm(IDA,IDB)) \end{aligned}$$

(Complete scripts and intermediate results in Appendix G)

Differential abundance analysis

To identify key ASVs driving the observed patterns of alpha and beta diversity, we performed a differential abundance analysis (DAA) using an Analysis of Compositions of Microbiomes with Bias Correction 2 (ANCOM-BC2) implemented in the R package ANCOMBC version 2.0.2 [61, 62]. ANCOM-BC2 allows for model fitting to the analysis with the ”fix formula=” option: Age, BCI, Rank, Habitat, Sex, Year, and Infection were specified as fixed effects. Unlike most DAA methods, ANCOM-BC2 also enables control for random effects with the option ”rand formula”. Again a nested random effect to account for Nest ID and Individual ID was fitted to the analysis (1|NestID/Individual ID). As part of ANCOM-BC2, the Holm-Bonferroni method was employed to correct P values for multiple testing [44]. A significance cutoff of \(p_{\text {adj}}\) < 0.05 was used. Default parameter settings were applied, and DAA was performed at ASV level. (Detailed scripts and intermediate results can be found in Appendix H)

Results

Gut microbiota composition

After quality control and filtering steps, the 16S rRNA dataset consisted of 117 individuals (66 males; 51 females) with a total of 230 samples (mean = 2 samples per individual) from 54 different nests. On average 22,268 ± 8,764.14 reads per sample were obtained, resulting in the identification of 2,078 ASVs (amplicon sequencing variants) (Appendix A Table S3). Two microbial kingdoms were identified, Bacteria (99.7%, SD = 0.24%) and Archaea (0.03%, SD = 0.24%). A total of 25 bacterial phyla were detected, of which five phyla were identified as core taxa: Firmicutes (39.3%, SD = 10.9%), Actinobacteria (33.1%, SD = 10.1%), Proteobacteria (19.1%, SD = 8.4%), Bacteroidota (3.4%, SD = 6.0%), and Campylobacterota (2.1%, SD = 3.5%). 219 bacterial families were identified, 19 of which were classified as core taxa. Among these families, the most abundant ones were: Corynebacteriaceae (15.9%, SD = 9.3%), Peptostreptococcaceae (12.0%, SD = 8.5%), Actinomycetaceae (9.5%, SD = 9.1%), Enterobacteriaceae (6.5%, SD = 7.3%), and Gemellaceae (6.1%, SD = 10.7%) (Fig. 1A).

Fig. 1
figure 1

Relative abundances (as percentages) of: A core bacterial gut microbiota phyla and families and B core eukaryotic gut microbiota phyla and families in common buzzard nestlings. Each individual is represented at two different sampling ages. Core taxa are defined as microbial taxa present in at least 70% of the samples. Note: In Figure 1B, P refers to Phylum, C - Classe, O - Order and F - Family

The 28S rRNA final dataset included 109 individuals (64 males; 45 females) and a total of 180 samples, representing 54 nests. The average number of reads per sample were 8,212 ± 4,231.45, resulting in the identification of 1770 ASVs (Appendix A Table S5). A total of 22 eukaryotic phyla were identified, with Ascomycota (57.1%, SD = 24.6%), Basidiomycota (25.4%, SD = 20.9%), and Phragmoplastophyta (7.3%, SD = 13.9%) being the core phyla. The majority of the core taxa could not be identified at family level. Among the 94 families, Cladosporiaceae (8.9%, SD = 17.7%) and unassigned families belonging to the class Exobasidiomycetes (23.8%, SD = 22.9%), order Capnodiales (11.4%, SD = 9.1%), phylum Phragmoplastophyta (8.9%, SD = 17.7%), were identified as the core taxa (Fig. 1B).

Factors driving gut microbiota alpha diversity

Age was the only factor that explained gut bacterial diversity. Both Shannon and Faith phylogenetic diversity decreased with age (Shannon: \(\beta\) = \(-\)0.17, \(\chi ^2\) = 9.76, p = 0.002; Faith PD: \(\beta\) = \(-\)0.03, \(\chi ^2\) = 7.82, p = 0.005) (Fig. 2A and B; Appendix A Table S8 & S9). The fixed effects of the model accounted for only a small portion of the diversity variation, with approximately 10% of the variability explained by the random effects (\(LMM_{\text {Shannon}}\): \(R^2_{\text {marginal}}\)=0.075, \(R^2_{\text {conditional}}\)=0.176; \(LMM_{\text {Faith PD}}\): \(R^2_{\text {marginal}}\)=0.092, \(R^2_{\text {conditional}}\)=0.128). However, the analysis of the significance of random effects revealed no evidence for variation driven by the grouping factors (Table 1). Similar patterns were found for eukaryotic microbial diversity. Shannon diversity decreased with age (\(\beta\) = \(-\)0.036, \(\chi ^2\) = 4.08, p = 0.043), while no evidence was found for Faith PD (\(\beta\) = \(-\)0.020, \(\chi ^2\) = 1.24, p = 0.23; Fig. 2C and D; Appendix A Table S10 & S11). Additionally, no significant evidence for variation driven by random effects was detected (Table 2). Evidence of an effect of habitat on Faith PD was found, but after adjusting for multiple comparisons with Benjamini-Hochberg method, habitat did not remain statistically significant (\(p_{\text {adj}}\) > 0.05; Appendix A Table S12). Overall, our results suggest that age is the primary factor influencing gut microbiota diversity in common buzzard nestlings.

Fig. 2
figure 2

Effect of age on microbiome alpha diversity as obtained from a LMM. p-Values are shown for each diversity metric along with the variance explained by the fixed effects (\(R^2_{\text {m}}\)) and the total variance explained by the model (\(R^2_{\text {c}}\)) A Bacterial Shannon diversity. B Bacterial Faith PD. C Eukaryotic Shannon diversity. D Eukaryotic Faith PD

Table 1 Analysis of significance of the random effects (“ranova”) with bacterial alpha diversity as response variable
Table 2 Analysis of significance of the random effects (”ranova”) with eukaryotic alpha diversity as response variable

Gut microbiota beta diversity

We found a substantial influence of age difference and nest sharing (among compared pairs) on the bacterial community structure (beta diversity). Bacterial microbiota composition became more distinct as the age difference between individuals increased (\(\mu\) = 0.38,CI [0.36, 0.41]). Conversely, dissimilarity decreased among individuals from the same nest (\(\mu _{\text {Same Nest}}\) = \(-\)0.45, CI [\(-\)0.48, \(-\)0.43]; Fig. 3A). In addition to age differences and nest sharing, we found evidence for the effects of \(\Delta\)BCI (\(\mu\) = 0.03, CI [0.01, 0.06]), year (\(\mu _{\text {Same year}}\) = \(-\)0.1, CI [\(-\)0.12, \(-\)0.08]) and habitat similarity (\(\mu _{\text {Same Habitat}}\) = \(-\)0.04, CI [\(-\)0.05, \(-\)0.02]) on bacterial composition, albeit less pronounced (Fig. 3A). Notably, individuals with similar BCI, individuals from the same year, and occupying the same habitat tended to exhibit similar bacterial microbiota compositions. Sex similarity (\(\mu _{\text {Same sex}}\) = \(-\)0.01, CI [\(-\)0.02, 0]) did not explain beta diversity (Fig. 3A). Blood parasite-infected individuals exhibited greater dissimilarities among themselves than did non-infected individuals among them (\(\mu _{\text {Ni-Ni}}\) = \(-\)0.1, CI [\(-\)0.29, \(-\)0.03]). Moreover, within infected individuals dissimilarities were higher than dissimilarities between non-infected/infected pairs (\(\mu _{\text {Ni-I}}\) = \(-\)0.1, CI [\(-\)0.13, \(-\)0.001]) (Fig. 3A). When considering phylogenetic relationships among bacterial communities (WU), credible intervals indicated an effect of age difference (\(\mu\) = 0.12, CI [0.10, 0.14]), nest sharing (\(\mu _{\text {Same Nest}}\) = \(-\)0.15, CI [\(-\)0.17, \(-\)0.13]), and year similarity (\(\mu _{\text {Same year}}\) = \(-\)0.05, CI [\(-\)0.07, \(-\)0.03]) (Fig. 3B).

Fig. 3
figure 3

Effect size estimates (points) and corresponding 95% credible intervals (horizontal lines) derived from Bayesian regression (brms) models with pairwise Bray-Curtis dissimilarities and Weighted UniFrac distances as response variables. All predictors correspond to pairwise comparisons: nest sharing, habitat, sex, and year (different/same), \(\Delta\)Age (age differences), \(\Delta\)BCI (body condition index differences) and infection status (infected-infected; infected/non infected; non-infected/non-infected). A variable significantly predicts microbiota dissimilarity/distance when the credible intervals do not overlap zero

Similar but less pronounced trends were observed in the eukaryotic microbiota. In terms of Bray-Curtis dissimilarity, nest sharing (\(\mu _{\text {Same Nest}}\) = \(-\)0.18, CI [\(-\)0.22, \(-\)0.14]) had the largest impact on compositional differences; \(\Delta\)age (\(\mu\) = 0.09, CI [0.05, 0.12]), \(\Delta\)BCI (\(\mu\) = 0.05, CI [0.0006, 0.09]), and habitat similarity exhibited weaker effects (Fig. 3C). When considering phylogenetic distances, nest sharing (\(\mu _{\text {Same Nest}}\) = \(-\)0.11, CI [\(-\)0.14, \(-\)0.09]) and \(\Delta\)age (\(\mu\) = 0.10, CI [0.08, 0.12]) remained as primary contributors to compositional differences, with \(\Delta\)BCI and habitat similarity (\(\mu _{\text {Same Habitat}}\) = \(-\)0.09, CI [\(-\)0.12, \(-\)0.07]) showing weaker effects (Fig. 3D). No effects of year, sex, or blood parasite infection infection were observed (Fig. 3C and D). Complete Bayesian model diagnostics and results are available in the supplementary information, Appendix A Tables S13 to S16.

Age-associated differential abundant taxa

We identified nine differential abundant bacterial ASVs associated with age (\(p_{\text {adj}}\) < 0.05; Appendix A Table S17). Among these, six ASVs displayed a considerable log-fold decrease with age and belonged to the genera Microbacterium, Pseudoxanthomonas, Bifidobacterium, and Tepidomonas, while one ASV was identified only at the family level, Bacillaceae. Conversely, three ASVs demonstrated a log-fold increase with age and were assigned with the genera Varibaculum and Enterococcus (Fig. 4A). No eukaryotic ASVs showed differential variation with age; however, two ASVs, both belonging to the order Dothideales, exhibited differential abundance between northern and southern habitats (Fig. 4B). No differential abundant taxa were found associated with any of the other studied variables.

Fig. 4
figure 4

Differential abundance analysis using ANCOM-BC2. A Log fold changes (LFC) of deferentially abundant bacterial taxa with increasing age. B Habitat pairwise comparisons showing log-fold changes of deferentially abundant eukaryotic taxa

Discussion

Gut microbial communities acquired during early-life stages significantly impact host condition by influencing metabolism [26, 42] and the development of the nervous and immune systems [41, 52]. Investigating the driving factors behind the colonisation of the gut is therefore essential for gaining a comprehensive understanding of the broader implications of the gut microbiota on the health and fitness of the host. In this study, we analysed longitudinal data from a wild population of common buzzards nestlings and explored relationships among the gut microbiome and environmental and host factors. We found that age shapes both bacterial and eukaryotic gut microbiota diversity and composition, with the nesting environment playing a crucial role in shaping microbiota composition. Nestling condition affects gut microbial composition, while habitat influences both bacterial and eukaryotic communities. Additionally, we observed associations between Leucocytozoon infection and bacterial beta diversity, although no such differences were detected in the eukaryotic microbiota.

Gut microbiota diversity is affected by nestling age

Our study revealed strong evidence for an influence of age on bacterial gut microbiota diversity, as evidenced by declines in both Shannon and Faith PD indices with increasing age (Fig. 2A and B). This suggests that, over time, certain bacterial taxa may become more dominant as nestlings age, leading to a reduction in overall diversity (competitive exclusion) [1, 29]. These dominant taxa also exhibit closer phylogenetic relationships, indicating potential specialisation or competitive advantages that drive their prevalence. While certain raptor studies have not found significant age-related effects [45, 107, 123], our findings align with studies on other species [34, 48, 115], suggesting a diverse gut microbiota in newborns due to rapid colonisation after hatching.

Maraci et al [71] reported fluctuating changes in gut microbiota diversity during ontogeny in zebra finches (Taeniopygia guttata) and Bengalese finches (Lonchura striata domestica). Our results from the differential abundance analysis reveal that some specific taxa exhibit a log-fold increase with age, while others decrease in abundance (Fig. 4A), providing additional support for these dynamic shifts in relation to age. These findings emphasise the importance of sampling an adequately extended time-frame to effectively capture variations in the gut microbiota diversity. Similarly, eukaryotic taxa showed a decline in Shannon diversity with age, although Faith PD showed no detectable effects (Fig. 2C and D). Unlike Shannon diversity, Faith PD is not weighted on abundance, which means that increasingly dominant taxa, responsible for reducing Shannon diversity, do not exert the same influence on Faith PD [38, 97]. These findings align with results from West et al [115], who documented a decline in mycobiota diversity with age in kākāpō chicks, mirroring the pattern previously identified in their bacterial microbiota study [115].

As birds typically hatch with little microbial diversity, this rapid initial colonisation most likely comes from the parents, either by direct contact or by food provisioning and passive exposure to the nesting environment [24, 109]. Buzzard nestlings fledge at around 35 days old [110]. During this period, parental care and presence (mainly by females) decrease gradually: from days 0–8 after hatching, active brooding occurs also outside of feeding sessions; from days 9–30, there is a decline in female presence and care, and by around 25–30 days, female care resembles that of males [46]. The decrease in microbiota diversity during this period may be related to reduced parental contact and dietary shifts [32, 43]. Empirical data suggests no major dietary shifts in the common buzzard, there are however shifts in feeding practices from direct provisioning to depositing prey in the nest resulting in nestlings transitioning from consuming individual pieces of prey to whole prey items. Additional work would be needed to test the impact of changing feeding habits on the development of the gut microbiota of buzzard chicks. Additionally, maturation of the immune system and other physiological systems can drive microbial selection leading to the observed patterns of decreasing diversity with age [108, 122].

Age difference and nest similarity driving gut compositional differences

Age differences and nest similarity exhibited the most pronounced effects on beta diversity. We observed that compositional dissimilarities and phylogenetic distances increase with age difference and between distinct nests, a trend consistently evident in both bacterial and eukaryotic diversity (Fig. 3). Studies by Worsley et al [118, 119] found that age primarily impacted bacterial compositional differences. Similarly, in humans and other primates, compositional differences linked to age have been consistently observed [74, 102, 104]. Coupled with the identification of specific taxa showing increased abundance while others decrease with age (Fig. 4A), it is plausible to hypothesize that these changes may correspond to changing functional requirements across the developmental stages of buzzard nestlings [84, 90].

A plethora of studies have reported robust links between the nesting environment and chick gut microbiota [26, 28, 36, 98]. Our findings align with the ”nidobiome” concept, which integrates the collective effects of parents, nest conditions, and nestling interactions during microbiome assembly and their contribution over the developmental stages [24].

Nestling condition shapes gut microbial composition but not diversity

Differences in nestling body condition impacted gut microbial composition, with both bacterial and eukaryotic communities exhibiting increased similarity among nestlings with similar body condition (Fig. 3). However, we did not detect individually significantly abundant taxa, indicating that compositional changes are relatively small and are distributed across the entire microbial community rather than being concentrated on a small number of taxa. Alterations in the microbiota community can potentially result in dysbiosis, which may have adverse effects on the host’s health [63]. This could explain the observed increased dissimilarities among nestlings with poor body condition compared to those in better health. However, an alternative plausible mechanism would be that changes in microbiota composition are a consequence of variation in individual body condition, considering that factors such as anatomy, physiology, and the immune system play pivotal roles in shaping gut microbial communities [31, 42]. The relationship between gut microbial composition and host body condition is intricate and varies across host species [42, 55]. Certain microbial taxa can either enhance or diminish host condition, depending on their interactions with the current diet, host immunity and other members of the microbiota [93]. A study conducted on the steppe buzzard (Buteo buteo vulpinus) showed no correlation between body condition and microbiota community composition, but it revealed a significant increase in the relative abundance of the genus Escherichia-Shigella with decreasing body condition [107]. Studies highlighting connections between eukaryotic gut microbiota and health primarily focus on economically important species [103], with a scarcity of direct connections observed in the wild. Nevertheless, research on the Seychelles warbler revealed an association between fungal microbiota and differential survival [119]. Our results add to a limited body of evidence suggesting that not only bacterial but also eukaryotic gut microbial communities are associated with host body condition.

In contrast with microbiota composition, no discernible connection was found between body condition and gut microbiota diversity. Gut microbiota diversity is influenced by a multitude of variables and these complex interactions may obscure direct links between body condition and microbiota diversity [42, 52, 105]. Furthermore, the potential time lag between changes in body condition and alterations in the gut microbiota might also contribute to the observed lack of association [33].

Habitat and year affects microbiota composition

Habitat similarity also affected gut microbiota composition with individuals from the same habitat sharing similar compositions (Fig. 3). Interestingly, this effect was more pronounced within the eukaryotic communities, resulting in substantial impacts on both dissimilarity and phylogenetic measurements (Fig. 3C and D). The study site encompasses three distinct habitat types: the Teutoburg Forest (Teuto), a low mountain region where Norway Spruce (Picea abies), Beech (Fagus sylvatica), and Oak (Quercus robur and Q. petraea) are dominant; the northern habitat characterised by a prevalence of Beech and Oak; and the southern habitat marked by the dominance of Pine (Pinus silvestris) and sandy soils [50, 57]. It is plausible that variation in tree and soil types contributes to the observed differences in gut microbiota composition, possibly reflecting diverse ecological niches that influence microbial transmission dynamics [42, 105]. Variation in tree and soil types likely contribute to the observed differences in gut microbiota composition. This suggests that eukaryotic microbiota can be acquired both directly from the environment and indirectly through prey items (particularly the common vole, Microtus arvalis). Indeed fungi often have more specialised substrate requirements compared to bacteria [27]. Different fungal species may have specific host plants or environmental substrates they prefer, leading to a stronger association with certain habitats. In contrast, bacteria might be more generalist and thus adaptable to diverse habitats. Notably, core microbiota taxa such as Exobasidiomycetes species have been linked to abnormal plant tissue outgrowths [10], and Cladosporiaceae species are commonly found in soil and plant materials [11]. Supporting the influence of habitat differences on eukaryotic microbiota, we also observed differential abundance of Dothideales between Northern and Southern habitats. Dothideales species have strong affinities with conifer trees, particularly Pinus species [14], which are dominant in the South habitat.

Year similarity also had an influence on bacterial microbiota composition, with no discernible effects observed within the eukaryotic communities. Seasonal variations, presumably (linked to fluctuations in food availability and weather conditions) [42] and ecological drift (random changes in the frequency of different microbial taxa) likely contribute to these differences [92, 96].

Despite the effects on microbiota composition, sampling year and habitat showed no impact on gut microbiota diversity. Generally, annual variation in prey availability and divergent environmental factors can also contribute to fluctuations in gut microbial alpha diversity [45, 53, 104]. It is, nevertheless, important to acknowledge the limitations of our sample set. Only a small percentage of individuals were sampled in 2021 (10%) compared to 2020, and similarly, few individuals were sampled from the South (10%) and Teuto (<5%) habitats compared to the North habitat. Consequently, this reduces statistical power and our capacity to make inferences about detected effects, while also potentially limiting our ability to detect differences.

Gut microbiota and Leucocytozoon infection

Severe malaria is caused by Plasmodium blood parasites, while related Leucocytozoon parasites can induce similar conditions in a wide range of avian hosts [5]. Despite a high prevalence (55%) of Leucocytozoon infection in buzzard nestlings previously found on this study site, research suggests low pathogenicity and minimal impacts on host condition [117]. Nevertheless, signatures of infection were found in haematological and blood chemistry profiles [117].

Our analysis revealed evidence that bacterial beta diversity, rather than alpha diversity, varied with infection status (Fig. 3A). Infected individuals displayed greater compositional dissimilarity among themselves than compared to non-infected nestlings. Furthermore, the compositional differences do not appear to be influenced by individually differential abundant taxa, but rather by overall changes across the community. No correlation was found between eukaryotic gut microbiota and infection status (Fig. 3C and D). In the context of infection, the host’s immune response can affect the gut microbiota composition; conversely, alterations in the microbiota following infection might signal the immune system to trigger a response [41, 60]. Mateos-Hernández et al [75] showed that altering the gut microbiota can enhance protection against avian aspergillosis by modulating anti-\(\alpha\)-Gal immunity. In the context of infection, the host’s immune response can significantly impact the gut microbiota composition, potentially leading to changes in the host microbiota [41].

Differences in beta diversity between bacterial and eukaryotic microbiota may arise from differential interactions with mucosal immunity (e.g. recognition mechanisms; immune cell responses), as well as differential microbiota-driven [41] signalling of the immune response [4, 18, 60]. Additionally, the higher abundance and broader range of functions performed by bacterial taxa might make small changes easier to detect [100].

While recognizing existing connections between gut microbiome and blood parasite infection [113], the mechanisms behind these relationships, especially among bird species, remain largely unknown. Although many studies find no correlations with alpha or beta diversity metrics, subtle links with specific taxa have been detected [7, 78, 91]. Further studies, particularly with controlled experimental designs, will be necessary to untangle the complex effects of parasite infection alongside other variables, while here we aim to shed light on population dynamics in their natural habitat.

Sex and rank does not impact the gut microbiota

Neither bacterial nor eukaryotic alpha and beta diversity differed between sexes. There is no sex dimorphism in nestling common buzzards and sex-related differences in gut microbial diversity may become more apparent in adulthood, as the birds mature and experience hormonal changes associated with sexual development [42, 119].

Rank exhibited no influence on nestling gut microbiota diversity. Similar to sex, this could be attributed to the relatively short time frame of this study, which may not have allowed hierarchical differences among chicks to become evident. The nesting environment appears to have a homogenising effect on microbiota diversity across the brood, surpassing other potential factors like dominance rank in its influence [26, 36].

Gut microbiota profile

There is a scarcity of studies investigating the microbiota composition of raptor species. Nevertheless, our findings align with a handful of existing studies, which show that Firmicutes, Proteobacteria, and Actinobacteriota (Fig. 1A) are among the most abundant bacterial phyla in raptor gut microbiota [73, 80, 107, 121]. Notably, while studies focusing on raptors lack descriptions of the eukaryotic microbiota, our results are in accordance with research conducted in other bird species that also revealed that the two major fungal phyla, Ascomycota and Basidiomycota, are predominant in the gut microbiota [68, 116, 119].

Our findings show that the core eukaryotic microbiota is predominantly dominated by fungal communities as the core microbiota (Fig. 1B). However, we did identify the presence of other eukaryotes, such as Arthropoda (4.8%, SD = 12.8%, most likely derived from diet, although not typically consumed by common buzzards), Apicomplexa (4.8%, SD = 12.8%), and Annelida (0.6%, SD = 3.0%), albeit in smaller abundances (Appendix A Table S6). These results show the utility of the 28S rRNA marker for studying the gut microbiome, while also highlighting the challenges in fully exploring the less abundant eukaryotic communities in the gut.

Conclusion

We present a study into the gut microbiota of wild common buzzard nestlings, encompassing both bacterial and eukaryotic components. Using a longitudinal approach, we established that age is a critical factor shaping both bacterial and eukaryotic gut microbiota diversity and composition. These age-related shifts underscore the importance of capturing an adequate time frame that can help to disentangle the temporal dynamics of gut microbiota development. Moreover, our study highlights the importance of the nesting environment, with the ”nidobiome” concept aptly summarising the combined effects of parents, nest, and nestlings on microbiome assembly and interactions. Nestling condition emerges as a determinant of gut microbial composition, reflecting host health, and habitat plays a role in shaping not only bacterial but also eukaryotic communities. While we also observed a correlation between Leucocytozoon infection and bacterial beta diversity, no such differences were detected in the eukaryotic microbiota. These different outcomes may result from distinct interactions between the mucosal immune system, bacteria and eukaryotes.

Supplementary information

Supplementary tables and figures can be found in Appendix  A

Availability of data and materials

All 16S and 28S rRNA raw reads have been submitted to the European Nucleotide Archive repository, Project ID: PRJEB70791. The scripts and metadata to reproduce all analyses can be accessed via the GitHub repository: https://github.com/hugoeira/Gut-microbiota-of-buzzard-nestlings.

Code available

Supplementary materials and can be accessed via the GitHub repository: https://github.com/hugoeira/gut-microbiota-of-buzzard-nestlings.

References

  1. Abt MC, Pamer EG. Commensal bacteria mediated defenses against pathogens. Curr Opin Immunol. 2014;29:16–22.

    Article  CAS  PubMed  Google Scholar 

  2. Anderson MJ. Permutational multivariate analysis of variance (permanova). Wiley statsref: statistics reference online 2014; pp 1–15.

  3. Apprill A, McNally S, Parsons R, et al. Minor revision to v4 region ssu rrna 806r gene primer greatly increases detection of sar11 bacterioplankton. Aquat Microb Ecol. 2015;75(2):129–37.

    Article  Google Scholar 

  4. Arrieta MC. Gut mycobiome core species causally modulate metabolic health in mice. researchsquare preprint 2024.

  5. Atkinson CT, Thomas NJ, Hunter DB. Parasitic diseases of wild birds. Hoboken: Wiley; 2009.

    Google Scholar 

  6. Van der Auwera G, Chapelle S, De Wächter R. Structure of the large ribosomal subunit rna of phytophthora megasperma, and phylogeny of the oomycetes. FEBS Lett. 1994;338(2):133–6.

    Article  PubMed  Google Scholar 

  7. Aželytė J, Wu-Chuang A, Maitre A, et al. Avian malaria parasites modulate gut microbiome assembly in canaries. Microorganisms. 2023;11(3):563.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Bartoń K. MuMIn: Multi-Model Inference. r package version 1.47.5 2023; https://CRAN.R-project.org/package=MuMIn.

  9. Bates D, Mächler M, Bolker B, et al. 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 

  10. Begerow D, Stoll M, Bauer R. A phylogenetic hypothesis of ustilaginomycotina based on multiple gene analyses and morphological data. Mycologia. 2006;98(6):906–16.

    Article  PubMed  Google Scholar 

  11. Bensch K, Braun U, Groenewald JZ, et al. The genus cladosporium. Stud Mycol. 2012;72:1–401.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Bergner LM, Orton RJ, Benavides JA, et al. Demographic and environmental drivers of metagenomic viral diversity in vampire bats. Mol Ecol. 2020;29(1):26–39.

    Article  PubMed  Google Scholar 

  13. Bijlmsa R. Sex determination of nestling common buzzards buteo buteo. Limosa. 1999;72:1–10.

    Google Scholar 

  14. Bills GF, Collado J, Ruibal C, et al. Hormonema carpetanum sp. nov., a new lineage of dothideaceous black yeasts from spain. Stud Mycol. 2004;50:149–57.

    Google Scholar 

  15. Bisanz JE. qiime2r: Importing qiime2 artifacts and associated data into r sessions, 2018; https://github.com/jbisanz/qiime2R, v0.99.

  16. Bolyen E, Rideout JR, Dillon MR, et al. Reproducible, interactive, scalable and extensible microbiome data science using qiime 2. Nat Biotechnol. 2019;37(8):852–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Borda-Molina D, Seifert J, Camarinha-Silva A. Current perspectives of the chicken gastrointestinal tract and its microbiome. Comput Struct Biotechnol J. 2018;16:131–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Borghi M, Renga G, Puccetti M, et al. Antifungal th immunity: growing up in family. Front Immunol. 2014;5:506.

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  Google Scholar 

  20. Bürkner PC. Advanced bayesian multilevel modeling with the r package brms. arXiv preprint arXiv:1705.11123, 2017a.

  21. Bürkner PC. brms: an r package for Bayesian multilevel models using stan. J Stat Softw. 2017;80:1–28.

    Article  Google Scholar 

  22. Callahan BJ, McMurdie PJ, Rosen MJ, et al. Dada2: high-resolution sample inference from illumina amplicon data. Nat Methods. 2016;13(7):581–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Callahan BJ, Sankaran K, Fukuyama JA, et al. Bioconductor workflow for microbiome data analysis: from raw reads to community analyses. F1000Research. 2016;5:1492.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Campos-Cerda F, Bohannan BJ. The nidobiome: a framework for understanding microbiome assembly in neonates. Trends Ecol Evol. 2020;35(7):573–82.

    Article  PubMed  Google Scholar 

  25. Chakarov N, Pauli M, Krüger O. Immune responses link parasite genetic diversity, prevalence and plumage morphs in common buzzards. Evol Ecol. 2017;31:51–62.

    Article  Google Scholar 

  26. Chen CY, Chen CK, Chen YY, et al. Maternal gut microbes shape the early-life assembly of gut microbiota in passerine chicks via nests. Microbiome. 2020;8:1–13.

    Article  Google Scholar 

  27. Chen Y, Xi J, Xiao M, et al. Soil fungal communities show more specificity than bacteria for plant species composition in a temperate forest in china. BMC Microbiol. 2022;22(1):208.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Costanzo A, Ambrosini R, Franzetti A, et al. The cloacal microbiome of a cavity-nesting raptor, the lesser kestrel (falco naumanni). PeerJ. 2022;10:e13927.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Coyte KZ, Rakoff-Nahoum S. Understanding competition and cooperation within the mammalian gut microbiome. Curr Biol. 2019;29(11):R538–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Coyte KZ, Schluter J, Foster KR. The ecology of the microbiome: networks, competition, and stability. Science. 2015;350(6261):663–6.

    Article  CAS  PubMed  Google Scholar 

  31. Cui L, Morris A, Ghedin E. The human mycobiome in health and disease. Gen Med. 2013;5(7):1–12.

    Google Scholar 

  32. David LA, Maurice CF, Carmody RN, et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2014;505(7484):559–63.

    Article  CAS  PubMed  Google Scholar 

  33. Davidson GL, Somers SE, Wiley N, 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 

  34. Davies CS, Worsley SF, Maher KH, et al. Immunogenetic variation shapes the gut microbiome in a natural vertebrate population. Microbiome. 2022;10(1):41.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Davis NM, Proctor DM, Holmes SP, et al. Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data. Microbiome. 2018;6:1–14.

    Article  Google Scholar 

  36. Diez-Méndez D, Bodawatta KH, Freiberga I, et al. Indirect maternal effects via nest microbiome composition drive gut colonization in altricial chicks. Mol Ecol. 2023;32(13):3657–71.

    Article  PubMed  Google Scholar 

  37. van Dongen WF, White J, Brandl HB, et al. Age-related differences in the cloacal microbiota of a wild bird species. BMC Ecol. 2013;13:1–12.

    Google Scholar 

  38. Faith DP, Baker AM. Phylogenetic diversity (pd) and biodiversity conservation: some bioinformatics challenges. Evol Bioinforma. 2006;2:117693430600200000.

    Article  Google Scholar 

  39. Fehérvári P, Solt S, Erdélyi K, et al. Hatching rank influences nutritional condition in the common buzzard: evidence from ptilochronology. J Raptor Res. 2014;48(3):280–4.

    Article  Google Scholar 

  40. Fridolfsson AK, Ellegren H. A simple and universal method for molecular sexing of non-ratite birds. J Avian Biol. 2023;30(1):116–21.

    Article  Google Scholar 

  41. Gensollen T, Iyer SS, Kasper DL, et al. How colonization by microbiota in early life shapes the immune system. Science. 2016;352(6285):539–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Grond K, Sandercock BK, Jumpponen A, et al. The avian gut microbiota: community, physiology and function in wild birds. J Avian Biol. 2018;49(11):e01788.

    Article  Google Scholar 

  43. Hoffmann C, Dollive S, Grunberg S, et al. Archaea and fungi of the human gut microbiome: correlations with diet and bacterial residents. PLoS ONE. 2013;8(6):e66019.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Holm S. A simple sequentially rejective multiple test procedure. Scand J Stat. 1979;6(2):65–70.

    Google Scholar 

  45. Houtz JL, Melo M, Therrien JF, et al. Disentangling relationships between physiology, morphology, diet, and gut microbial diversity in American kestrel nestlings. J Avian Biol. 2023;7–8(2023):e03019.

    Article  Google Scholar 

  46. Hubert C, Gallo A, Le Pape G. Buzzard (buteo buteo). J Raptor Res. 1995;29(2):103–9.

    Google Scholar 

  47. Huffnagle GB, Noverr MC. The emerging world of the fungal microbiome. Trends Microbiol. 2013;21(7):334–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Janiak MC, Montague MJ, Villamil CI, et al. Age and sex-associated variation in the multi-site microbiome of an entire social group of free-ranging rhesus macaques. Microbiome. 2021;9:1–17.

    Article  Google Scholar 

  49. Johnson KVA, Burnet PW. Microbiome: should we diversify from diversity? Gut microbes. 2016;7(6):455–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Jonker RM, Chakarov N, Krüger O. Climate change and habitat heterogeneity drive a population increase in common buzzards buteo buteo through effects on survival. Ibis. 2014;156(1):97–106.

    Article  Google Scholar 

  51. Katoh K, Misawa K, Ki Kuma, et al. Mafft: a novel method for rapid multiple sequence alignment based on fast fourier transform. Nucl Acids Res. 2002;30(14):3059–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Kirschman LJ, Khadjinova A, Ireland K, et al. Early life disruption of the microbiota affects organ development and cytokine gene expression in threespine stickleback. Integr Compar Biol. 2023;63(1):250–62.

    Article  CAS  Google Scholar 

  53. Knutie SA. Food supplementation affects gut microbiota and immunological resistance to parasites in a wild bird species. J Appl Ecol. 2020;57(3):536–47.

    Article  CAS  Google Scholar 

  54. Kogut MH. Impact of the gut microbiota on the immune system. In: Avian Immunology. Elsevier, 2022; pp. 353–364.

  55. Kohl KD, Brun A, Bordenstein SR, et al. Gut microbes limit growth in house sparrow nestlings (passer domesticus) but not through limitations in digestive capacity. Integr Zool. 2018;13(2):139–51.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Kounosu A, Murase K, Yoshida A, et al. Improved 18s and 28s rdna primer sets for ngs-based parasite detection. Sci Rep. 2019;9(1):15789.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Krüger O. The importance of competition, food, habitat, weather and phenotype for the reproduction of buzzard buteo buteo. Bird Study. 2004;51(2):125–32.

    Article  Google Scholar 

  58. Kuznetsova A, Brockhoff PB, Christensen RHB, et al. Package ‘lmertest’. R package version. 2015;2:734.

    Google Scholar 

  59. Laforest-Lapointe I, Arrieta MC. Microbial eukaryotes: a missing link in gut microbiome studies. MSystems. 2018;3(2):e00201-17.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Lambring CB, Siraj S, Patel K, et al. Impact of the microbiome on the immune system. Crit Rev Immunol. 2019;39(5):313–28.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Lin H, Peddada SD. Analysis of compositions of microbiomes with bias correction. Nat Commun 2020;11(1):1–11. https://www.nature.com/articles/s41467-020-17041-7.

  62. Lin H, Eggesbo M, Peddada SD. Linear and nonlinear correlation estimators unveil undescribed taxa interactions in microbiome data. Nat Commun 2022;13(1):1–16. https://www.nature.com/articles/s41467-022-32243-x.

  63. Logan AC, Jacka FN, Prescott SL. Immune-microbiota interactions: dysbiosis as a global health issue. Curr Allergy Asthma Rep. 2016;16(2):13.

    Article  PubMed  Google Scholar 

  64. Lozupone CA, Hamady M, Kelley ST, et al. Quantitative and qualitative β diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microbiol. 2007;73(5):1576–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Lüdecke D, Ben-Shachar MS, Patil I, et al. performance: an R package for assessment, comparison and testing of statistical models. J Open Source Softw. 2021;6(60):313. https://doi.org/10.21105/joss.03139.

    Article  Google Scholar 

  66. Machida RJ, Knowlton N. PCR primers for metazoan nuclear 18s and 28s ribosomal DNA sequences. PLOS ONE. 2012. https://doi.org/10.1371/journal.pone.0046180.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Madlala T, Okpeku M, Adeleke MA. Understanding the interactions between eimeria infection and gut microbiota, towards the control of chicken coccidiosis: a review. Parasite. 2021;28.

  68. Mahtab N, Zhou L, Zhang F, et al. Seasonal variations in the gut fungal communities of hooded crane (grus monacha) at wintering and stopover sites in china. Animals. 2021;11(4):941.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Makowski D, Ben-Shachar MS, Patil I, et al. Methods and algorithms for correlation analysis in R. J Open Source Softw. 2020;5(51):2306.

    Article  Google Scholar 

  70. Mammeri M, Obregón DA, Chevillot A, et al. Cryptosporidium parvum infection depletes butyrate producer bacteria in goat kid microbiome. Front Microbiol. 2020;11:548737.

    Article  PubMed  PubMed Central  Google Scholar 

  71. Maraci Ö, Antonatou-Papaioannou A, Jünemann S, et al. Timing matters: age-dependent impacts of the social environment and host selection on the avian gut microbiota. Microbiome. 2022;10(1):1–20.

    Article  Google Scholar 

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

    Article  Google Scholar 

  73. Martinez-Hernandez JE, Berrios P, Santibáñez R, et al. First metagenomic analysis of the andean condor (vultur gryphus) gut microbiome reveals microbial diversity and wide resistome. PeerJ. 2023;11:e15235.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Martino C, Dilmore AH, Burcham ZM, et al. Microbiota succession throughout life from the cradle to the grave. Nat Rev Microbiol. 2022;20(12):707–20.

    Article  CAS  PubMed  Google Scholar 

  75. Mateos-Hernández L, Risco-Castillo V, Torres-Maravilla E, et al. Gut microbiota abrogates anti-α-gal iga response in lungs and protects against experimental aspergillus infection in poultry. Vaccines. 2020;8(2):285.

    Article  PubMed  PubMed Central  Google Scholar 

  76. McMurdie PJ, Holmes S. phyloseq: an r package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8(4):e6121. https://doi.org/10.1371/journal.pone.0061217.

    Article  CAS  Google Scholar 

  77. Moreno J, Briones V, Merino S, et al. Beneficial effects of cloacal bacteria on growth and fledging size in nestling pied flycatchers (ficedula hypoleuca) in spain. Auk. 2003;120(3):784–90.

    Google Scholar 

  78. Navine AK, Paxton KL, Paxton EH, et al. Microbiomes associated with avian malaria survival differ between susceptible Hawaiian honeycreepers and sympatric malaria-resistant introduced birds. Mol Ecol. 2022;32(23):6659–70.

    Article  PubMed  Google Scholar 

  79. Oldereid TS, Jiang X, Nordhus KS, et al. Role of bacteria and microbial metabolites in immune modulation during early life. Scand J Immunol. 2024;99(2):e13336.

    Article  CAS  PubMed  Google Scholar 

  80. Oliveira BC, Murray M, Tseng F, et al. The fecal microbiota of wild and captive raptors. Animal Microb. 2020;2:1–9.

    Article  CAS  Google Scholar 

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

  82. Paulson JN, Pop M, Bravo HC. metagenomeseq: statistical analysis for sparse high-throughput sequencing. Biocond Pack. 2013;1:191.

    Google Scholar 

  83. Paulson JN, Stine OC, Bravo HC, et al. Differential abundance analysis for microbial marker-gene surveys. Nat Methods. 2013;10(12):1200–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Pereira AC, Bandeira V, Fonseca C, et al. Egyptian mongoose (herpestes ichneumon) gut microbiota: taxonomical and functional differences across sex and age classes. Microorganisms. 2020;8(3):392.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Price MN, Dehal PS, Arkin AP. Fasttree 2-approximately maximum-likelihood trees for large alignments. PLoS ONE. 2010;5(3):e9490.

    Article  PubMed  PubMed Central  Google Scholar 

  86. Quast C, Pruesse E, Yilmaz P, et al. The silva ribosomal RNA gene database project: improved data processing and web-based tools. Nucl Acids Res. 2012;41(D1):D590–6.

    Article  PubMed  PubMed Central  Google Scholar 

  87. R Core Team. R: A language and environment for statistical computing. R foundation for statistical computing, Vienna, Austria, 2022; https://www.R-project.org/.

  88. Raulo A, Allen BE, Troitsky T, et al. Social networks strongly predict the gut microbiota of wild mice. ISME J. 2021;15(9):2601–13.

    Article  PubMed  PubMed Central  Google Scholar 

  89. Robeson MS, O’Rourke DR, Kaehler BD, et al. Rescript: Reproducible sequence taxonomy reference database management. PLoS Comput Biol. 2021;17(11):e1009581.

    Article  PubMed  PubMed Central  Google Scholar 

  90. Rohrer SD, Jiménez-Uzcátegui G, Parker PG, et al. Composition and function of the galapagos penguin gut microbiome vary with age, location, and a putative bacterial pathogen. Sci Rep. 2023;13(1):5358.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Rohrer SD, Robertson BQ, Chubiz LM, et al. Gut microbiome composition associated with plasmodium infection in the Eurasian tree sparrow. J Avian Biol. 2023;1–2:e03027.

    Article  Google Scholar 

  92. Rose C, Lund MB, Schramm A, et al. Does ecological drift explain variation in microbiome composition among groups in a social host species? J Evol Biol. 2023;36(12):1684–94.

    Article  PubMed  Google Scholar 

  93. Rosshart SP, Vassallo BG, Angeletti D, et al. Wild mouse gut microbiota promotes host fitness and improves disease resistance. Cell. 2017;171(5):1015–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  95. Ruiz-Rodríguez M, Soler JJ, Lucas FS, et al. Bacterial diversity at the cloaca relates to an immune response in magpie pica pica and to body condition of great spotted cuckoo clamator glandarius nestlings. J Avian Biol. 2009;40(1):42–8.

    Article  Google Scholar 

  96. Seki D, Schauberger C, Hausmann B, et al. Individuality of the extremely premature infant gut microbiota is driven by ecological drift. Msystems. 2022;7(3):e00163-22.

    Article  PubMed  PubMed Central  Google Scholar 

  97. Shannon CE. A mathematical theory of communication. Bell Syst Tech J. 1948;27(3):379–423.

    Article  Google Scholar 

  98. Somers SE, Davidson GL, Johnson CN, et al. Individual variation in the avian gut microbiota: the influence of host state and environmental heterogeneity. Mol Ecol. 2023;32(12):3322–39.

    Article  PubMed  Google Scholar 

  99. Stanley D, Hughes RJ, Moore RJ. Microbiota of the chicken gastrointestinal tract: influence on health, productivity and disease. Appl Microbiol Biotechnol. 2014;98:4301–10.

    Article  CAS  PubMed  Google Scholar 

  100. Stefanaki C, Peppa M, Mastorakos G, et al. Examining the gut bacteriome, virome, and mycobiome in glucose metabolism disorders: are we on the right track? Metabolism. 2017;73:52–66.

    Article  CAS  PubMed  Google Scholar 

  101. Stoffel MA, Acevedo-Whitehouse K, Morales-Durán N, et al. Early sexual dimorphism in the developing gut microbiome of northern elephant seals. Mol Ecol. 2020;29(11):2109–22.

    Article  PubMed  Google Scholar 

  102. Strati F, Di Paola M, Stefanini I, et al. Age and gender affect the composition of fungal population of the human gastrointestinal tract. Front Microbiol. 2016;7:1227.

    Article  PubMed  PubMed Central  Google Scholar 

  103. Summers KL, Arfken AM. The gut mycobiome and animal health. In: Gut microbiota, immunity, and health in production animals. Springer, 2022; p 85–125.

  104. Sun B, Gu Z, Wang X, et al. Season, age, and sex affect the fecal mycobiota of free-ranging tibetan macaques (macaca thibetana). Am J Primatol. 2018;80(7):e22880.

    Article  PubMed  Google Scholar 

  105. Sun F, Chen J, Liu K, et al. The avian gut microbiota: Diversity, influencing factors, and future directions. Front Microbiol. 2022;13:934272.

    Article  PubMed  PubMed Central  Google Scholar 

  106. Teyssier A, Lens L, Matthysen E, et al. 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 

  107. Thie N, Corl A, Turjeman S, et al. Linking migration and microbiota at a major stopover site in a long-distance avian migrant. Mov Ecol. 2022;10(1):46.

    Article  PubMed  PubMed Central  Google Scholar 

  108. Thomson CA, Morgan SC, Ohland C, et al. From germ-free to wild: modulating microbiome complexity to understand mucosal immunology. Mucosal Immunol. 2022;15(6):1085–94.

    Article  CAS  PubMed  Google Scholar 

  109. Trevelline BK, MacLeod KJ, Knutie SA, et al. In ovo microbial communities: a potential mechanism for the initial acquisition of gut microbiota among oviparous birds and lizards. Biol Let. 2018;14(7):20180225.

    Article  Google Scholar 

  110. Tyack A, Walls S, Kenward R. Behaviour in the post-nestling dependence period of radio-tagged common buzzards buteo buteo. Ibis. 1998;140(1):58–63.

    Article  Google Scholar 

  111. Videvall E, Song SJ, Bensch HM, et al. Early-life gut dysbiosis linked to juvenile mortality in ostriches. Microbiome. 2020;8:1–13.

    Article  Google Scholar 

  112. Videvall E, Paxton KL, Campana MG, et al. Transcriptome assembly and differential gene expression of the invasive avian malaria parasite plasmodium relictum in hawaii. Ecol Evol. 2021;11(9):4935–44.

    Article  PubMed  PubMed Central  Google Scholar 

  113. Waide ML, Schmidt NW. The gut microbiome, immunity, and plasmodium severity. Curr Opin Microbiol. 2020;58:56–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  Google Scholar 

  115. West AG, Digby A, Lear G, et al. Influence of management practice on the microbiota of a critically endangered species: a longitudinal study of kākāpō chick faeces and associated nest litter. Animal Microbiome. 2022;4(1):55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  116. West AG, Digby A, Team KR, et al. The mycobiota of faeces from the critically endangered kākāpō and associated nest litter. New Zealand Journal of Zoology 2023; https://doi.org/10.1080/03014223.2023.2170428.

  117. Wiegmann A, Springer A, Rinaud T, et al. The prevalence of leucocytozoon spp. in nestlings of three wild raptor species including implications on haematological and blood chemistry values. Int J Parasitol: Parasit Wildl. 2021;16:236–43.

    Google Scholar 

  118. Worsley SF, Davies CS, Mannarelli ME, et al. Gut microbiome composition, not alpha diversity, is associated with survival in a natural vertebrate population. Animal Microb. 2021;3:1–18.

    Article  Google Scholar 

  119. Worsley SF, Davies CS, Mannarelli ME, et al. Assessing the causes and consequences of gut mycobiome variation in a wild population of the seychelles warbler. Microbiome. 2022;10(1):242.

    Article  PubMed  PubMed Central  Google Scholar 

  120. Yang S, Gao X, Meng J, et al. Metagenomic analysis of bacteria, fungi, bacteriophages, and helminths in the gut of giant pandas. Front Microbiol. 2018;9:1717.

    Article  PubMed  PubMed Central  Google Scholar 

  121. Zhang K, Gong X, Wang X, et al. Comparison of the gut microflora of common kestrel (falco tinnunculus) according to differences in sex, age, and physiological condition. Preprint at 2022; https://www.researchsquare.com/article/rs-1857423/v1.

  122. Zhang L, Zhan H, Xu W, et al. The role of gut mycobiome in health and diseases. Ther Adv Gastroenterol. 2021;14:17562848211047130.

    Article  CAS  Google Scholar 

  123. Zhou L, Huo X, Liu B, et al. Comparative analysis of the gut microbial communities of the eurasian kestrel (falco tinnunculus) at different developmental stages. Front Microbiol. 2020;11:592539.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are grateful to Thomas Grünkorn, and Justine Léauté for their efforts in collecting data during the field season. Thanks to Elke Hippauf for the support during laboratory procedures. We extend our thanks to Julia Hassa and Anika Winkler from the NGS team of the Omics CF NGS Unit and CeBiTec for their meticulous work in the sequencing process. We acknowledge Marc Gilles, Rebecca Chen and Kosmas Hench for general advice and assistance with statistical analysis. We acknowledge the financial support of the German Research Foundation (DFG). The funding body played no role in the design of the study and collection, analysis, interpretation of data, and in writing the manuscript.

Funding

Open Access funding enabled and organized by Projekt DEAL. This work was supported by the DFG, as part of the SFB TRR 212 (NC3) - Project Nos. 316099922 and 396780709; DFG Project No. 398434413; DFG Project No. 433069365 and DFG Project No. 233740704.

Author information

Authors and Affiliations

Authors

Contributions

HP, NC, JIH and OK conceptualised the research project. HP conducted field and laboratory work, carried out bioinformatic and statistical analysis and drafted the manuscript with the input from NC, JIH, OK. NC, MO, TR, KPG and OK contributed to field work/data collection. TB conducted Illumina sequencing. BC and ÖM provided scientific advice and contributed to the development of the methodology. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Hugo Pereira.

Ethics declarations

Ethics approval and consent to participate

Sampling of common buzzards nestlings followed the ARRIVE guidelines and were approved by the Animal Ethics Committee at Bielefeld University and permitted by the local authority Kreis Gütersloh, permit number: 4.5.2-723-Bussard and by the ethics committee of the Animal Care and Use Committee of the German North Rhine-Westphalia State Office for Nature, Environment and Consumer Protection (Landesamt für Nature, Umwelt und Verbraucherschutz Nordrhein-Westfalen) under permit numbers: 84-02.04.2014.A091, 84-02-04.2017.A147.

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

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

Pereira, H., Chakarov, N., Hoffman, J.I. et al. Early-life factors shaping the gut microbiota of Common buzzard nestlings. anim microbiome 6, 27 (2024). https://doi.org/10.1186/s42523-024-00313-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s42523-024-00313-8

Keywords