- 1Department of Animal Behaviour, Bielefeld University, Bielefeld, Germany
- 2Department of Evolutionary Population Genetics, Bielefeld University, Bielefeld, Germany
- 3Joint Institute for Individualisation in a Changing Environment, Bielefeld University and University of Münster, Bielefeld, Germany
- 4British Antarctic Survey, Cambridge, United Kingdom
- 5Center for Biotechnology (CeBiTec), Faculty of Biology, Bielefeld University, Bielefeld, Germany
- 6Department of Wildlife Diseases, Leibniz Institute for Zoo and Wildlife Research, Berlin, Germany
- 7Department of Behavioural Ecology, Bielefeld University, Bielefeld, Germany
- 8Department of Behavioural Biology, University of Münster, Münster, Germany
The gut microbiota-immune-brain axis is a feedback network which influences diverse physiological processes and plays a pivotal role in overall health and wellbeing. Although research in humans and laboratory mice has shed light into the associations and mechanisms governing this communication network, evidence of such interactions in wild, especially in young animals, is lacking. We therefore investigated these interactions during early development in a population of common buzzards (Buteo buteo) and their effects on individual condition. In a longitudinal study, we used a multi-marker approach to establish potential links between the bacterial and eukaryotic gut microbiota, a panel of immune assays and feather corticosterone measurements as a proxy for long-term stress. Using Bayesian structural equation modeling, we found no support for feedback between gut microbial diversity and immune or stress parameters. However, we did find strong relationships in the feedback network. Immunity was negatively correlated with corticosterone levels, and microbial diversity was positively associated with nestling body condition. Furthermore, corticosterone levels and eukaryotic microbiota diversity decreased with age while immune activity increased. The absence of conclusive support for the microbiota-immune-brain axis in common buzzard nestlings, coupled with the evidence for stress mediated immunosuppression, suggests a dominating role of stress-dominated maturation of the immune system during early development. Confounding factors inherent to wild systems and developing animals might override associations known from adult laboratory model subjects. The positive association between microbial diversity and body condition indicates the potential health benefits of possessing a diverse and stable microbiota.
1 Introduction
A growing body of evidence highlights a crucial relationship among gut microbes, the environment and host immunity, which also molds the maturation and function of the central nervous system (CNS) (Fung, 2020; Chin et al., 2020). This intricate feedback system encompasses the brain, the autonomic and enteric nervous systems, the endocrine and immune systems, and the gut microbiome (Martin et al., 2018). Gut-brain communication occurs through various mechanisms, including neuroanatomical pathways and the neuroendocrine axis (hypothalamic-pituitary-adrenal or HPA axis). However, increasing importance has been given to the interactions between the microbiota and the immune system as a fundamental pathway regulating microbiota-gut-brain communication (Fachi et al., 2019). For example, the products resulting from the metabolic activities of the gut microbiota include bioactive peptides such as neurotransmitters, short-chain fatty acids, intestinal hormones, and branched-chain amino acids that play pivotal roles in regulating communication between the gut, brain and the immune system (Sperringer et al., 2017; Cryan et al., 2019; Ding et al., 2020; Krishnamurthy et al., 2023). These metabolites are able to enter the circulatory system to communicate with the brain, eliciting stimulation of the HPA axis (Lynch et al., 2024; Sudo et al., 2004). Additionally, they exert a direct influence on the mucosal immune system, where specific metabolites can function as immune signaling molecules, activating or enhancing systemic immune responses (Ding et al., 2020; Long-Smith et al., 2020). In turn, the brain regulates the gastrointestinal tract and overall organismal homeostasis (Mindus et al., 2021).
The gut microbiota contributes to the development of the host immune system during early life stages, regulates it to maintain gut homeostasis and protects the host from the colonization of potential pathogens (Lo et al., 2021; Berman et al., 2023). As such increasing importance has been given to the neonatal microbiome and its developmental trajectory during early life which affects host immunity and neural activity maturation (Jašarević and Bale, 2019; Ratsika et al., 2023). For example, in humans, the mode of delivery significantly influences the normal colonization of the gut microbiota (Ratsika et al., 2023). C-section deliveries have been linked to immune disturbances (allergies and asthma) (Roduit et al., 2009; Bisgaard et al., 2011) and disruptions in the structure and function of the CNS, with indications of increased neuronal cell death (Castillo-Ruiz et al., 2018) and deficits in early communication skills and social behavior later in life (Morais et al., 2020). In murine models, penicillin treatment during late pregnancy and early postnatal life has been shown to induce alterations in the gut microbiota, leading to heightened cytokine expression, altered integrity of the blood–brain barrier, and significant differences in behavior (Ratsika et al., 2023). Interactions among these functional systems, particularly during sensitive developmental periods, carry consequences later in life, influencing behavior, disease susceptibility, general health and survival (Strange et al., 2016; Francella et al., 2022; Ratsika et al., 2023).
Research also shows the importance of glucocorticoids (GC) in gut-brain communication, with increased GC levels affecting microbial diversity and composition (Petrullo et al., 2022). Various factors, including, predation (Mohring et al., 2023), heat stress (Lin et al., 2006), high energetic demands during reproduction (Fletcher et al., 2015), and food availability (Romero and Wikelski, 2001) can affect host homeostasis and are known to induce elevated stress levels in animals. Corticosterone, a key glucocorticoid in various vertebrate species, is recognized as a classic endocrine response to stress but also for its role in energy regulation (Almasi et al., 2009). GCs mediate ongoing stress responses, either via maintaining basal levels, allowing other aspects of the stress response to act efficiently, or by actively triggering the stress response (Sapolsky et al., 2000). An alternative view suggests that GCs may suppress the stress response, preventing detrimental over-activation (Sapolsky et al., 2000).
Substantial evidence emphasizes the importance of the microbiota-immune-brain axis in humans and laboratory animals (Fung, 2020). Yet, the understanding of these links in wild animals remains limited (Hird, 2017; Davidson et al., 2020). Studies in birds associated exploratory behavior with microbiota diversity, while learning and memory performance have been correlated with compositional differences (Florkowski and Yorzinski, 2023; Slevin et al., 2020). The gut microbiota's links to stress have also been explored: in common toad tadpoles (Bufo bufo), elevated baseline corticosterone associates with higher microbial diversity (Gabor et al., 2022), while American red squirrels (Tamiasciurus hudsonicus) show lower alpha diversity and fewer gastrointestinal pathogens in response to elevated glucocorticoids (Petrullo et al., 2022). Studies of the microbiota-immune axis in wild barn swallows (Hirundo rustica) and Egyptian fruit bats (Rousettus aegyptiacus) demonstrate that antigen challenges (phytohaemagglutinin and lipopolysaccharides, respectively) can induce changes in gut microbiota composition, which in turn predict the strength of the immune response (Kreisinger et al., 2018; Berman et al., 2023). However, since the majority of these studies have involved adult individuals, it is uncertain whether similar patterns would be observed in nestlings and juveniles. In a rare exception, Stoffel et al. (2020) found that a small proportion of the variation in beta diversity among northern elephant seal juveniles was explained by health status (assessed by counting various white blood cell populations) yet a clear pattern emerged where healthier individuals exhibited higher microbiota diversity.
Here, we investigated the gut microbiota-immune-brain axis in a wild vertebrate population, drawing on predictions based on human and murine models. Using a wild population of common buzzards (Buteo buteo), we sampled nestlings and collected information on bacterial and eukaryotic microbiota diversity (microbiota component). Long-term stress was measured from feather corticosterone (brain component). In addition we performed a series of immune assays (immune component) and estimated the body condition of each individual. Furthermore, we assessed each component at two distinct time points throughout nestling development to incorporate developmental trajectories and dynamic changes. Our approach diverges from traditional methodologies, which often rely on single-point estimates. All components were subsequently integrated into a structural equation modeling framework (Grace, 2006; Lefcheck, 2016) incorporating the following assumptions: 1. All components of the axis exert an influence on the body condition of the nestlings; 2. The immune system is influenced by both stress levels and the microbiota; 3. Corticosterone levels affect microbiota diversity. In common buzzards we have previously shown that as individuals mature, microbiota diversity declines suggesting a shift toward a stable community following uncontrolled colonization after hatching (Pereira et al., 2024). Building on that knowledge we expect that, as individuals age their gut microbiota diversity will decline, their immune system will mature (increased immune capacity), and stress levels will decrease, leading to a state of homeostasis and improved body condition.
2 Methods
2.1 Study system
The common buzzard (Buteo buteo) is a long-lived raptor (up to 25 years) widely distributed across Europe (Walls and Kenward, 2020). This study examines a population that has been closely monitored for over 25 years within a 300 km2 area (Krüger, 2004; Jonker et al., 2014). Predominantly a resident species, common buzzards breed in tall trees (over 10 meters) from March to July. Adults typically form long-term pairs, build stick nests, and rear on average two nestlings per brood, with reproductive success closely tied to vole availability, their primary food (Kostrzewa and Kostrzewa, 1990; Panek, 2021). During their ≈5-week nestling period (≈35 days), nestlings are entirely dependent on parental care and feeding. Even after fledging (≈47 days), they stay close to the nest and continue receiving parental support (Fachi et al., 2019). Parental care (mostly by the female) decreases during this period: initially, females provide active brooding beyond just feeding times (0–8 days), but as the nestlings grow, female presence decreases (9–30 days), eventually mirroring the care level of males (Hubert et al., 1995). While diet seems to remain consistent (as inferred from observational data), there is a shift in feeding practices: parents transition from directly provisioning nestlings to depositing prey in the nest, allowing the nestlings to gradually consume whole prey items.
We used existing microbiome data from Pereira et al. (2024) and complemented these with new data on feather corticosterone and a set of immune markers. We successfully obtained all sample types from a set of 43 common buzzard individuals, each sampled at two different time points during their nestling phase, totalling 86 samples. Individuals were sampled from 23 nests with an average brood size per nest of two nestlings (Table 1; see Supplementary Tables S1, S2 for n° of ASVs per marker). These were sampled across two different habitats: North of the Teutoburg Forest and south of the Teutoburg Forest (8° 25′E and 52° 6′N; Eastern Westphalia, Germany), as described in Krüger (2004) and Jonker et al. (2014).
2.2 Sample collection
In brief, body weight and wing length were recorded at each sampling event, with a nine-day interval on average. The first sampling occurred on average at 19 days of age (mean ± s.d. = 19.3 ± 5.29 days) and the second at 28 days (mean ± s.d. = 27.8 ± 5.16 days). Due to difficulties in precise age estimation prior to sampling (the hatching date was estimated by the number of droppings on the ground below the nest; age was determined following the initial visit to the nest and subsequent measurement of wing length), it was impossible to sample individuals at exactly the same age at sampling points 1 and 2. Nestling age was calculated (post-first sampling) using a sex-specific polynomial regression model on wing length (Bijlmsa, 1999). Body condition index (BCI) was determined by extracting the residuals of a logarithmic regression of weight on wing length, adjusting for sex. Blood samples (500 μl) were collected from the ulnar vein and stored in heparinized tubes. A small blood drop was used for smears, a portion was stored in ethanol for sex determination, and the remaining volume was centrifuged in the field until separation of plasma was visible. Separated plasma was transferred to a new tube, stored in dry ice, and subsequently stored at −80°C until further analysis. The remaining red blood cells were resuspended in PBS solution and stored at −20°C. Cloacal swabs for gut microbiota analysis were obtained and stored in RNAlater, first in dry ice and then long term at −80°C. In order to assess corticosterone levels, one interscapular feather was pulled from each bird and these were individually stored in paper envelopes.
2.3 Microbiome profiling: DNA isolation, sequencing, and data processing
For detailed procedures on microbiome sequencing data processing see Pereira et al. (2024).
2.3.1 DNA isolation and sequencing
Cloacal swabs underwent DNA extraction using a modified phenol-chloroform protocol. For gene library preparation, the “Illumina 16S Metagenomic Library Preparation Guide” (15044223 Rev.B) was followed. A multimarker approach was used, targeting the V4 region of the 16S ribosomal RNA (rRNA) gene with the primers 515F (Parada) (Parada et al., 2016) and 806R (Apprill) (Apprill et al., 2015). To be able to capture not only bacteria but also eukaryotes, the D8-D9 region of the 28S rRNA gene was targeted with the primers GA20F (Van der Auwera et al., 1994)/RM9Rb (Machida and Knowlton, 2012). PCRs were conducted in 25 μl reaction volumes containing 5 μl DNA, 12.5 μl KAPA HiFi HotStart ReadyMix, 1 μl of each primer (1 μM), and 6 μl of PCR-grade water. Index-PCRs utilized Illumina Nextera XT V2 index kits. Libraries were equimolarly pooled, and sequenced on an Illumina MiSeq platform (0.4% MiSeq run) with a 2 × 300 bp paired-end reads protocol.
2.3.2 16S rRNA gene sequence data processing
Sequence data were imported into QIIME2 (Quantitative Insights Into Microbial Ecology 2, version 2022.11) (Bolyen et al., 2019). Quality assessment was done by visualizing quality plots and Amplicon Sequencing Variants (ASVs) were inferred using the Divisive Amplicon Denoising Algorithm pipeline (DADA2) (Callahan et al., 2016). Taxonomy was assigned using a SILVA 138.1-trained naive Bayes taxonomic classifier (Quast et al., 2012). Contaminants were identified and removed with the decontam package version 1.18 (Davis et al., 2018). ASVs assigned to Mitochondria, Chloroplast, Vertebrata, Eukaryota, and unassigned taxa were filtered out using QIIME2. Singletons and samples with a minimum frequency below 500 reads were removed. The remaining ASVs were aligned with MAFFT (Katoh et al., 2002) and used to construct a phylogeny with FastTree (Price et al., 2010), both implemented in QIIME2.
2.3.3 28S rRNA gene sequence data processing
Demultiplexed Illumina sequence data were imported into R version 4.2.2 (R Core Team, 2022). Locus-specific primers were removed using Cutadapt (version 4.4) (Martin et al., 2018). In QIIME2, quality assessment was performed through visualization of quality plots. ASV inference was conducted using DADA2 following the methodology outlined by Callahan et al. (2016). Sequences were trimmed to eliminate low quality regions and paired-end reads were concatenated. In QIIME2, taxonomy was assigned using the naive Bayes taxonomic classifier trained on the SILVA 138.1 database. The decontam pipeline in R was applied to remove contaminants and to perform 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 in QIIME2. The resulting ASVs were aligned using MAFFT, and a phylogeny was constructed using FastTree.
2.4 Assessment of immunity
Due to the complexity of the immune system, we measured four innate immune parameters [bacterial killing activity against Escherichia coli (BKA), lysozyme and haptoglobin concentrations, natural antibodies (HA) and complement (HL) titers] and one component of the acquired immune system [total immunoglobulin Y (IgY) concentration]. All these assays are regularly used in wild bird species, including raptor nestlings, both in comparative and within-species studies (see references below).
2.4.1 Bacterial killing activity
The bacterial killing activity (BKA) against Escherichia coli (ATCC No 8739) was used to characterize the functional activity of a bird's constitutive innate immune system (Irene Tieleman et al., 2005) using a spectrophotometric version of the assay (Nebel et al., 2021; Brust et al., 2022; Vincze et al., 2022). This assay measures the plasma's capacity to kill microbes ex vivo, determining the organism's ability to eliminate bacterial pathogens encountered, providing an environmentally relevant immune response (Millet et al., 2007; Tieleman, 2018). The assay evaluates the synergic function of several immune components (humoral ones in case of working with plasma), including antibacterial enzymes, complement components, and natural antibodies (French and Neuman-Lee, 2012). Briefly, 12 μl of 1:7 PBS-diluted sample was pipetted in duplicate into 96-well-plate and mixed with 4 μl of ≈1.5 × 105 colony-forming units (CFU)/ml. A positive (not containing any plasma) and a negative control (not containing any E. coli or plasma) was run on each plate. After incubation for 30 min at 37°C, 83μl of tryptic soy broth (#22092, Fluka) was added to each well. Absorbance at 300 nm was measured with a spectrophotometer (Biotek; μQuant Microplate Spectrophotometer) to determine background absorbance and again after the plates had been incubated for 12 h at 37°C. The BKA was quantified as the bacteria growth in plasma after 12 h (in %) subtracted by the background absorption in relation to the positive control (Brust et al., 2022).
2.4.2 Lysozyme
Lysozyme is an antibacterial enzyme that causes rapid cell lysis, especially in Gram-positive bacteria. It is part of the constitutive innate immune system and is often measured to assess inflammation-induced levels in plasma (Millet et al., 2007). To measure its concentration in plasma, we used the lysoplate assay (Prüter et al., 2020; Brust et al., 2022): 10 μl of sample was inoculated in the test holes of a 1% Noble agar gel (A5431, Sigma) containing 50 mg/100 ml lyophilized Micrococcus lysodeikticus (M3770, Sigma), a bacterium which is particularly sensitive to lysozyme concentration. Crystalline hen egg white lysozyme (L6876, Sigma) (concentration: 0.5, 0.8, 1, 2, 4, 8, 10, 20, and 40 μg/ml) was used to prepare a standard curve for each plate. After 20h incubation at 37°C, a clear zone developed in the area of the gel surrounding the sample inoculation site corresponding to the bacterial lysis. The diameters of the cleared zones are proportional to the log of the lysozyme concentration. This area was measured three times digitally using the software ImageJ (version 1.48, ImageJ) and the mean was converted to a semi-logarithmic plot into hen egg lysozyme equivalents (HEL equivalents, expressed in μg/mL) according to the standard curve (Prüter et al., 2020; Brust et al., 2022).
2.4.3 Haptoglobin
Haptoglobin, is an acute-phase protein in birds and is part of the induced innate immune system. Acute phase proteins are key indicators of immunological function. Their concentrations fluctuate over time, reflecting changes in health and physiological condition (Hõrak et al., 2002, 2003). Levels can rise quickly in response to infection, inflammation, or trauma (Millet et al., 2007). Elevated plasma haptoglobin often signifies the onset of a non-specific immune response (Matson et al., 2012). We measured haptoglobin concentrations with a commercial kit (TP801, Tri-Delta Diagnostics, Inc.) following the instructions of the manufacturer. Haptoglobin concentrations (mg/ml) in undiluted plasma samples were calculated according to the standard curve on each plate (Prüter et al., 2020; Nebel et al., 2021).
2.4.4 Haemolysis–haemagglutination assay
The levels of the natural antibodies and complement were assessed by using a haemolysis-haemagglutination assay as described by Matson et al. (2005), Nebel et al. (2021), and Brust et al. (2022). Natural antibodies, quantified as Haemagglutination (HA) titers, bind non-specifically to various antigens and play crucial roles in opsonization (Ochsenbein et al., 1999; Forthal, 2014). Haemolysis (HL), facilitated by the complement system, is a component of the innate immune system, its activation leads to cell lysis, particularly of cells that have been opsonized (Nauta et al., 2004). After pipetting 10 μl of plasma into the first two columns of a U-shaped 96-well microtiter plate, 10 μl sterile PBS was added to columns 2–12. The content of the second column wells was serially diluted (1:2) until the 11th column, resulting in a dilution series for each sample from 1/1 to 1/,1024. The last column of the plate was used as negative control, containing PBS only. Ten microliter of 1% rabbit red blood cells (supplied by Innovative Research) suspension was added to all wells and incubated at 37°C for 90 min. After incubation, in order to increase the visualization of agglutination, the plates were tilted at a 45° angle at room temperature. Agglutination and lysis, which reflect the activity of the natural antibodies and the interaction between these antibodies and complement (Matson et al., 2005; Prüter et al., 2020), was recorded after 20 and 90 min, respectively. Haemagglutination is characterized by the appearance of clumped red blood cells as a result of antibodies binding multiple antigens, while during haemolysis, the red blood cells are destroyed. Titers of the natural antibodies and complement were given as the log2 of the reciprocal of the highest dilution of serum showing positive haemagglutination or lysis, respectively (Matson et al., 2005; Prüter et al., 2020; Brust et al., 2022).
2.4.5 Total immunoglobulin Y concentration
Total IgY, the avian equivalent to mammalian immunoglobulin G, is the primary humoral effector of the adaptive immune system, playing a critical role in neutralizing pathogens (Warr et al., 1995). IgY is essential for long-term immunity, providing protection against recurrent infections. Total IgY was measured using an ELISA with commercial anti-chicken antibodies (Bourgeon et al., 2010; Hanssen et al., 2013). Ninety-six well high-binding ELISA plates (82.1581.200, Sarstedt) were coated with 100 μl of diluted plasma sample (1:4,000 diluted in carbonate–bicarbonate buffer, in duplicates) and incubated first for 1 h at 37°C and then overnight at 4°C. After incubation, the plates were washed with a 200 μl solution of PBS and PBS–Tween, before 100 μl of a solution of 1% gelatine in PBS–Tween was added. Plates were then incubated at 37°C for 1 h, washed with PBS–Tween and 100 μl of polyclonal rabbit anti-chicken IgY conjugated with peroxidase (A-9046, Sigma) at 1:250 (v/v) was added. Following 2 h incubation at 37°C, the plates were washed again with PBS–Tween three times. After washing, 100 μl of revealing solution [peroxide diluted 1:1,000 in ABTS (2,20-azino-bis-(3-ethylbenzthiazoline-6-sulphonic acid))] was added, and the plates were incubated for 1 h at 37°C. The final absorbance was measured at 405 nm using a photometric microplate reader (μQuant Microplate Spectrophotometer, Biotek) and subsequently defined as total serum IgY levels (Brust et al., 2022).
2.5 Feather-corticosterone determination
The entire feather was weighed and placed in a tube. For every 10 mg of feather, 1 ml of 100% methanol (p.a.) was added, followed by crushing with scissors. Subsequently, the samples were subjected to an ultrasonic bath incubation at 30°C for 30 min, followed by overnight incubation in an overhead shaker (25 rpm). On the following day, centrifugation (2 min, 19,800 × g) was performed, and the supernatant transferred to a new tube. The feather pellet underwent two washes, each with half of the extraction volume (2 min each at 19,800 × g); the resulting supernatants were then combined and one further centrifugation was performed for 10 min at 19,800 × g. The supernatant was then filtered through a PTFE membrane (Merck Millipore: Ultrafree-CL), and the filtrate was temporarily stored at −20°C. A defined amount of methanol (500 μl) from the temporarily stored filtrate was evaporated in a vacuum concentrator. The samples were re-suspended in 250 μl ELISA buffer (Cayman Chemicals Inc. [# 501320]) and stored at −20°C for at least overnight before conducting the ELISA. Corticosterone concentrations were determined in triplicate following the manufacturer's instructions. f-CORT assay validation is presented in Supplementary Data Sheet 6.
2.6 Statistical analysis
In order to evaluate sequencing depth and sample coverage, rarefaction curves were constructed in QIIME2. Rarefaction was then applied to the 16S rRNA dataset at 4,000 reads and to the 28S rRNA dataset at 2,000 reads. Utilizing the q2-diversity alpha plugin, three alpha diversity metrics were calculated: Shannon diversity index (Shannon, 1948); Faith's Phylogenetic Diversity (Faith PD) (Faith and Baker, 2006) and number of observed ASVs. Prior to constructing the Structural Equation Model (SEM), we explored associations of the variables of interest (BCI, f-CORT, immune assays) with host and environmental factors examined by Pereira et al. (2024). These included sex, habitat and rank (dominance hierarchy within the brood). We have previously established that none of these variables affects microbiota alpha diversity (Pereira et al., 2024). This preliminary analysis aimed to assess the relevance of these variables for the construction of the SEM. For each variable of interest, individual linear mixed models (LMMs) with a Gaussian distribution were constructed using the lmer function from the lme4 package in R (Bates et al., 2015). The significance of factors was assessed through analysis of variance (ANOVA). To account for repeated samples and nest sharing, Nest ID and Individual ID were incorporated as nested random effects (Individual ID nested within Nest ID). Additionally, the Benjamini-Hochberg method was applied to correct the p-values for multiple hypothesis testing (Benjamini and Hochberg, 1995). Scripts and results are presented in Supplementary Data Sheet 2.
2.6.1 Structural equation modeling approach
In order to tackle the complex interactions within the gut microbiota-immune-brain axis and their impact on individual condition, we employed a SEM strategy. SEM is a statistical methodology that allows for the simultaneous testing of complex relationships among multiple variables. Integrating factor analysis and multivariate regression analysis, SEM provides enhanced flexibility, enabling variables to both depend on and influence other variables. Furthermore, SEM allows for the incorporation of mediation effects, facilitating the quantification of direct, indirect, and total effects (Grace, 2006; Lefcheck, 2016). Of particular relevance to this study is SEM's capability to construct and model latent variables: ones that cannot be directly measured but are hypothesized to exist (Bollen, 2002; Lefcheck, 2016).
2.6.2 Exploratory factor analysis
We started by constructing a latent variable representing Immunity; for this, an exploratory factor analysis (EFA) was performed using the “lavaan” (version 0.6-16) package (Rosseel, 2012) in R. The latent variable model incorporated all immune assays as indicators, with “cluster= Individual ID” set to address repeated sampling and Full Information Maximum Likelihood (FIML) utilized for handling missing data. Estimation was performed using Maximum Likelihood with robust standard errors (MLR), and standardization of values for the latent variable was implemented (“std.lv =T”). Model fit was evaluated through the examination of fit indices, including the chi-squared p-value, Comparative Fit Index (CFI), Root Mean Square Error of Approximation (RMSEA), and Standardized Root Mean Square Residual (SRMR). The conventional “rule of thumb”: CFI > 0.9; RMSEA < 0.08; SRMR < 0.08 was used for model fit evaluation (Hu and Bentler, 1999; MacCallum et al., 1996; Sharma et al., 2005). Predicted values for the latent construct were extracted using the “lav.predict” function from “lavaan” with default settings. During the initial assessment, it became evident that haptoglobin concentration did not significantly load on the latent variable Immunity (Figure 1), leading to lower model fit (Table 2). As a result, haptoglobin concentration was subsequently removed and treated individually.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1413976/fmicb-15-1413976-HTML/image_m/fmicb-15-1413976-g001.jpg)
Figure 1. Latent variable model for “Immunity”. Arrows point to different immune assays, which serve as components of the latent variable. Values are indicating the strength of the relationship between the latent variable and its indicators. Red values indicate variables that loaded significantly in the model while black values denote non-significant loadings. ***p < 0.01; *p < 0.05.
2.6.3 Bayesian SEM
A Bayesian structural equation model (Kaplan and Depaoli, 2012) was fitted using the brms R package (Bürkner, 2017; Bürkner, 2018), which incorporates the probabilistic coding language Stan (Carpenter et al., 2017). Four distinct models were inputted into the multivariate analysis in order to test our main hypotheses:
Path 1:
Path 2:
Path 3:
Path 4:
Path 1 allowed us to infer the contributions of the different components of the axis for nestling body condition, while the remaining paths were designed to capture the dynamics of the gut microbiota-immune-brain. Statistical constraints impose that we assume directionality a priori, although we are aware that these are bidirectional relationships. To account for repeated samples and nest sharing, a nested random effect (individual ID nested within nest ID) was incorporated into the model. Given our sampling design, where individuals of different ages were sampled at each time point (age as a continuous variable), the only way to account for age effects was to include age as a fixed effect. The model was run with four chains, each run with 100k iterations, a warm-up phase of 50k iterations and default priors. Model fit was assessed by examining the convergence of the runs, mixing of chains and performing posterior predictive checks (comparing predicted vs. observed posterior distribution). Marginal and conditional Bayes R2 (Gelman et al., 2019) were calculated using the “bayes R2” function in brms. Scripts and intermediate results can be found in Supplementary Data Sheets 3, 4.
2.7 Differential abundance analysis
Analysis of Compositions of Microbiomes with Bias Correction 2 (ANCOM-BC2) with default parameters was implemented in the R package ANCOMBC version 2.0.2 (Lin and Peddada, 2020; Lin et al., 2022). Body condition index, f-CORT, Immunity and age were specified as fixed effects. A nested random effect to account for Individual ID within Nest ID was fitted with the option “rand formula”. The Holm-Bonferroni method with a significance cutoff of padj < 0.05 was used to correct P-values for multiple testing (Holm, 1979). Detailed scripts and results are presented in Supplementary Data Sheet 5.
3 Results
We examined the potential contribution of sex, habitat and brood rank to the variation of BCI, f-CORT, and the various immune assays. The analysis revealed no substantial evidence linking body condition, corticosterone levels, or immune capacity to these variables. However, slight variations in lysozyme and haptoglobin levels were observed among different habitats. These patterns appear to be largely influenced by imbalances in the sample design and the presence of outliers (Supplementary Figure S9). Thus, we assume that the tested variables should not play a significant role in the SEM construct. For a comprehensive overview of the analysis pipeline and detailed results, see Supplementary Data Sheet 2.
3.1 Exploratory factor analysis
Exploratory factor analysis (EFA) derived a latent variable representing the underlying structure of immune parameter values. We synthesized the results of the immune assays into a composite variable denoted as Immunity, which included all of the immune parameters except haptoglobin concentration (see above). Components of the innate immune system (Hemolysis and Hemagglutination) displayed the strongest factor loadings (λHA = 0.94, λHL = 0.90), followed by IgY (λIgY = 0.50) representing the adaptive immune system (Figure 1).
3.2 Bayesian SEM
The Bayesian SEM revealed no evidence for a relationship between gut microbiota diversity and either Immunity/Haptoglobin or f-CORT (Figures 2, 3; Supplementary Figures S1–S3). However, a negative correlation emerged between Immunity and f-CORT, showing a decrease in immune capacity with rising levels of f-CORT (Figures 2–4A; Supplementary Figures S1, S2). Conversely, no evidence was found for a link between Haptoglobin and f-CORT (Figures 2, 3; Supplementary Figure S3). A negative association was evident between f-CORT and BCI (Figures 2–4B; Supplementary Figures S1–S3). Shannon bacterial diversity was associated with BCI, indicating that higher diversity levels corresponded to elevated BCI (μ = 0.18; CI [0.01, 0.35]; Figures 2, 4C). Conversely, no evidence was found for the effects of Immunity/Haptoglobin on body condition (Figures 2, 3; Supplementary Figures S1–S3). Eukaryotic microbiota diversity exhibited no connection with BCI (Figures 2, 3; Supplementary Figures S1–S3). We found no evidence for age-related effects on bacterial microbiota diversity (Figures 2, 3; Supplementary Figures S1, S2). However, there was credible support for a decrease of eukaryotic Faith PD and the number of observed ASVs with age [Faith PD: μ = −0.28; CI (−0.55, −0.01); Number of ASVs: μ = −0.28 CI (0.56, −0.01); Supplementary Figure S4]. Additionally, f-CORT levels decreased as individuals matured, while Immunity increased with age (Figures 2, 3; Supplementary Figures S1–S4). As the results for the number of observed ASVs were similar to those for Faith PD, they are presented in Supplementary Figure S1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1413976/fmicb-15-1413976-HTML/image_m/fmicb-15-1413976-g002.jpg)
Figure 2. Diagrams representing the Bayesian structural equation models for the different bacterial microbiota diversity metrics. Illustrated are the results of the models incorporating the latent variable, with results for the Haptoglobin assay superimposed. Values represent the estimated mean effect of a predictor (μ) on the outcome variable. Credible intervals (95% CI) are provided for predictors that explain the outcome variable.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1413976/fmicb-15-1413976-HTML/image_m/fmicb-15-1413976-g003.jpg)
Figure 3. Diagrams representing the Bayesian structural equation models for the different eukaryotic microbiota diversity metrics. Illustrated are the results of the models incorporating the latent variable, with results for the Haptoglobin assay superimposed. Values represent the estimated mean effect of a predictor (μ) on the outcome variable. Credible intervals (95% CI) are provided for predictors that explain the outcome variable.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1413976/fmicb-15-1413976-HTML/image_m/fmicb-15-1413976-g004.jpg)
Figure 4. Regression plots illustrating the main results of the SEM model. (A) Negative relationship between immunity and f-CORT; (B) Negative correlation between body condition and f-CORT; (C) Increase in bacterial microbiota Shannon diversity with body condition.
3.3 Deferentially abundant taxa
In the ANCOM-BC2 models incorporating the latent variable Immunity, two eukaryotic ASVs exhibited differential abundance. Specifically, one ASV was linked to BCI (Genus: Dothideales), while another was associated with age (Phylum: Phragmoplastophyta) (see Supplementary Figure S5). However, the latter failed to retain statistical significance following sensitivity analysis (see Supplementary Table S21). The genus Dothideales demonstrated a log-fold decrease per unit of BCI (Supplementary Figure S5). Notably, no bacterial taxa displayed differential abundance correlated with the studied variables (Supplementary Figure S5). In contrast, when incorporating the Haptoglobin immune assay into the models, no deferentially abundant taxa were identified (Supplementary Figure S6; Supplementary Table S22).
4 Discussion
The gut microbiota-immune-brain axis shapes a variety of physiological responses through multidirectional communication among the gut microbiota, immune system, and the central nervous system (CNS) (Sylvia and Demas, 2018). This feedback system is known to influence immune function, neuroendocrine pathways, and behavior (Martin et al., 2018; Cryan et al., 2019). Once established, disruptions to this axis may lead to the development of disorders (e.g., irritable bowel syndrome; ulcerative colitis; Alzheimer's and Parkinson's disease in humans) significantly impacting individual health (Rhee et al., 2009; Ghaisas et al., 2016). Here, we used gut microbial diversity measurements, various innate and adaptive immune markers and feather-corticosterone levels (a proxy for stress) to explore relationships among the three components of the axis and the repercussions to body condition in a wild population of common buzzard nestlings.
4.1 Gut microbial diversity, not associated with Immunity and stress
In contrast to our assumptions, our results revealed no association between microbial diversity, Immunity and stress (f-CORT) (Figures 2, 3). However, several parts of the feedback cascade appeared to be functional (see below). Despite the well-established influence of the gut microbiota on immune system modulation and its impact on the CNS, these tripartite connections remain elusive in wild populations (Hird, 2017; Davidson et al., 2020; Madden et al., 2022; Pereira et al., 2023), particularly during the early stages of development, in contrast to the studies in humans and mice (Martin et al., 2018; Francella et al., 2022; Lynch et al., 2024).
Contrary to our findings, recent research by Francella et al. (2022) demonstrated clear links between gut microbiota, the immune system and stress in the early-life stages of laboratory-reared mice. Their study showed that immunocompromized mice had increased stress levels, decreased microbial diversity, and alterations in gut microbiome composition post-weaning; additionally, stress impacted the abundance of specific taxa that, in turn, were associated with specific behavioral traits. Furthermore, these behavioral changes were observed after just 1 week of age, demonstrating that stress can interact with host immunity during early development. Discrepancies between the results of the mice study and ours can be explained by various non-mutually exclusive factors. It is important to consider that their research was conducted under laboratory conditions, which lack ecologically relevant environmental factors, whereas our study was performed under ecologically realistic, natural circumstances. Indeed captivity is known to have an effect on each components of the axis (Slevin et al., 2020; Florkowski and Yorzinski, 2023), the studies also focused on different aspects of the immunity (cellular vs. humoral) and made use of immunocompromised/gene knockout mice. Furthermore species have different diets and life-history strategies in terms of fast–slow life-history continuum (Réale et al., 2010; Bing et al., 2022).
Despite the scarcity of such distinct gut microbiota-immune-brain studies in other vertebrate species with lower degrees of experimental manipulation, studies on captive organisms have demonstrated a connection between the brain and the gut microbiota. For instance, higher gut microbiota alpha diversity in zebra finches (Taeniopygia guttata) has been correlated with elevated exploratory behavior, and in house sparrows (Passer domesticus), beta diversity was associated with enhanced cognitive performance (Slevin et al., 2020; Florkowski and Yorzinski, 2023). It has been recognized that the gut microbiome modulates stress responses, particularly through the hypothalamic-pituitary-adrenal (HPA) axis (Foster and Neufeld, 2013); in general, higher stress levels tend to be associated with reduced microbiota diversity, specifically a reduction in the amount of rare and pathogenic taxa (Petrullo et al., 2022). While the notion of diverse microbiota being crucial for robust immunity is widely accepted (Hooper et al., 2012; Lozupone et al., 2012), extending these conclusions beyond a small number of experimental study systems is not straightforward, as several studies have found no association between microbiota diversity and immune indices. Instead, they observed that compositional differences or changes in specific taxa were more closely linked with variations in immune markers (Kreisinger et al., 2018; van Veelen et al., 2020; Fleischer et al., 2024). Our focus on microbial diversity, without considering composition (beta diversity), limits the depth of our analysis (Shade, 2017; Reese and Dunn, 2018). However, incorporating the complexities of multivariate metrics into structural equation modeling remains challenging. Nevertheless, we hope that ongoing advancements in methodologies for analyzing compositional data (Sweeny et al., 2023; Fountain-Jones et al., 2024) will facilitate this approach in future studies. It is also important to consider timing in the context of stress responses, as short-term and prolonged stress can induce contrasting physiological reactions (Martin, 2009). Evaluating different measures of HPA axis activity would offer a more comprehensive understanding of stress response dynamics (Stothart et al., 2019). Short-term acute stress triggers an organism's immune response (Martin, 2009) and enhances intestinal mucus secretion (Castagliuolo et al., 1996). Conversely, prolonged stress typically suppresses immunity and reduces mucus production, impacting the microbiota differently (Estienne et al., 2010; Shigeshiro et al., 2012).
Several studies of wild populations have delved into some components of the gut microbiota-immune-brain axis (Noguera et al., 2018; Petrullo et al., 2022; Berman et al., 2023), yet few have comprehensively addressed the entire system. An exception is a recent study on eastern newts (Notophthalmus viridescens), which explored the effects of experimentally elevated CORT levels on various immune indices and skin microbiome but did not find any clear evidence for a relationship among the three components (Pereira et al., 2023).
The examples provided demonstrate the challenges encountered in establishing associations among the components of the gut microbiome-immune-brain axis, either due to the complexities of wild settings or the specific time window investigated, or indeed, because there are no discernible associations in the corresponding systems. These challenges increase in less controlled study systems (further aggravated by the fact that studies of completely wild populations are almost non-existent), particularly when dealing with systems where certain components of the axis lack experimental challenge, alteration, or complete knockout. Given this, and considering the absence of an evident connection with the microbiome axis, we propose the following possible, not mutually exclusive mechanisms to be at play: 1. The gut microbiota may not exert a significant influence in the initial stages of life, especially when compared to the impacts of stress and the immune system. This is supported by the robust link between f-CORT and Immunity, which we find despite the correlative nature of our study system (Figures 2–4A). Consequently, we suggest that stress and immune regulation may hold greater importance for maintaining homeostasis during early development. 2. The complexity of wild environments introduces many background effects, potentially including diverse diets, exposure to various pathogens, and the unpredictable nature of ecological interactions. These factors might overshadow the subtle and context-dependent relationships observed in more controlled settings. In reality, the connections found in these controlled environments may not be as crucial or representative of the complete picture encountered in the more realistic complexities of a wild setting. 3. Direct comparisons between species and generalizations across species might be challenging due to ecological differences between the systems.
4.2 Stress, Immunity and body condition
The relationship between stress, immunity, and body condition is a strong feedback mechanism in animal physiology (Vagasi et al., 2018). Baseline levels of corticosterone play important roles in metabolism, development, reproduction, behavior, and immunity (Sapolsky et al., 2000). While beneficial in the short term for resolving inflammation and preventing an overshoot of the immune responses (Dhabhar, 2018), prolonged stress-induced elevation of corticosterone compromises the immune system over time (Dhabhar and Mcewen, 1997). A study of house sparrows revealed that prolonged activation of the stress response inhibits components of the innate immune system, such as complement-mediated lysis, bacteria-killing ability, and agglutination. Indeed, numerous studies have demonstrated the detrimental effects of high stress levels, for prolonged periods of time, on overall body condition, survival, and reproductive success (Angelier et al., 2010; Mikkelsen et al., 2023; Quirici et al., 2021).
We used feather CORT as an indicator of long-term stress. In birds, circulating CORT accumulates in developing feathers, serving as a cumulative gauge of hormone concentrations during feather growth (Jenni-Eiermann et al., 2015; Romero and Fairhurst, 2016). Our results show a negative association between Immunity and f-CORT, supporting the established concept that immunosuppression is expected in the face of allostatic overload (chronic stress) (Romero et al., 2009). Additionally, prolonged stress can incur costs for developing individuals, may divert resources away from essential physiological processes, and impair an individual's ability to fend off infections and maintain overall wellbeing, as evidenced here by a decline in body condition (Figures 2–4A, B).
4.3 Bacterial microbiota diversity and body condition
Our results show that nestlings with higher Shannon diversity have better body condition (Figures 2, 4C). This positive association can be attributed to the enhanced resistance of more diverse gut communities against pathogenic invasion, increased stability and resilience to disturbance (Buffie and Pamer, 2013). Additionally, a more diverse microbiota can potentially offer a broader range of functions performed by various bacterial taxa, leading to benefits for the host (Heiman and Greenway, 2016). Conversely, lower microbiota diversity is typically viewed as being detrimental to hosts (Le Chatelier et al., 2013), as it implies a loss of essential functions that could result in reduced nutrient assimilation or immunodeficiency (Round and Mazmanian, 2009; Hanning and Diaz-Sanchez, 2015). However, it has also been shown that high diversity might be associated with a state of dysbiosis, so a reduction in diversity could signal a return to homeostasis (Johnson and Burnet, 2016; Kohl et al., 2018; Coyte and Rakoff-Nahoum, 2019). Our results align with the Anna Karenina principle, which states that changes in microbiota due to disturbances lead to shifts from stable to unstable community states (Zaneveld et al., 2017). We propose that the increase in Shannon diversity in buzzard nestlings indicates an increase in stable and abundant taxa, which may offer greater benefits to the host (Hanning and Diaz-Sanchez, 2015). Metrics like Faith's PD and the number of ASVs treat all taxa equally (Faith and Baker, 2006), potentially explaining why they do not show associations with body condition (Figure 3; Supplementary Figure S1).
The absence of clear links between eukaryotic microbiota diversity and body condition may indicate a slower rate of change in the eukaryotic microbiota. Bacterial and eukaryotic taxa are likely to have distinct roles within the gut ecosystem (Oever and Netea, 2014; Vemuri et al., 2020). Fast-changing bacterial taxa may experience stronger competition, leading to positive selection on the core functionally-relevant taxa (Abt and Pamer, 2014; Coyte and Rakoff-Nahoum, 2019), whereas changes in the more complex and less abundant eukaryotic taxa (Chin et al., 2020) may occur at a slower pace and therefore not be observed during early stages of host development.
4.4 Age effects on stress, Immunity and eukaryotic microbiota diversity
The positive relationship between age and Immunity, and the negative association of age and stress levels (f-CORT; Figures 2, 3), most likely reflect nestling development and the maturation of the immune system. Initially, nestlings rely on innate and maternal transferred immunity (Klasing and Leshchinsky, 1999; Palacios et al., 2009). As they mature, exposure to pathogens and foreign microbes through the environment and diet, contact with nest material and siblings, challenges and stimulates the immune system, thus aiding its maturation (Lochmiller and Deerenberg, 2000; Morais et al., 2020; Oldereid et al., 2023). During the ontogenetic development of birds, exposure to corticosterone occurs in both the prenatal (embryonic) and postnatal (nestling and fledgling) periods (Henriksen et al., 2011; Strange et al., 2016). In the prenatal stage, corticosterone is transferred from the mother to the embryo through the egg yolk, and is influenced by the maternal environment, e.g., exposure to predators, competitors and other stress-inducing cues during egg production, and differences in environmental quality (Hayward and Wingfield, 2004; Saino et al., 2005; Love et al., 2008). In the postnatal stage, an initial surge in corticosterone levels may serve as an adaptive response to stressors associated with hatching, exposure to a new environment, and nutritional demands, being beneficial in the short term (Chin et al., 2009; Spencer et al., 2009; Strange et al., 2016). As nestlings mature, the HPA axis undergoes maturation, leading to improved stress response regulation. Initially, nestlings invest significant energy in growth and development. As they progress through early life stages, maturation enables more efficient allocation and prioritization of resources (Smulders, 2021; Spencer et al., 2009). As mentioned earlier, the maturation of both physiological systems is not independent; instead, a bidirectional interaction regulates both the immune system and the HPA axis (Francella et al., 2022).
We find a decline in eukaryotic microbiota phylogenetic diversity and the number of ASVs with age. Dominant taxa show higher resilience to disturbance, securing their positions via selective filtering and the occupation of core niches (Costello et al., 2012; Abt and Pamer, 2014; Coyte and Rakoff-Nahoum, 2019). Rare taxa are primarily acquired through stochastic processes but still significantly contribute to certain community alpha diversity measures (like Faith PD and n° of observed ASVs). Over time, these rare taxa are progressively replaced by the more dominant taxa (Shade et al., 2014). It's important to emphasize that microbiota colonization can be either deterministic or stochastic. Studies in gnotobiotic animals highlight this variability, demonstrating that some bacterial strains colonize in a deterministic manner, while others do so stochastically (Vega and Gore, 2017; Jones et al., 2022; Hayashi et al., 2024). Understanding how these processes contribute to community assembly during development in wild populations remains challenging. The dominant taxa also show closer phylogenetic relationships to one another, suggesting potential specialization or competitive advantages driving their prevalence (Janiak et al., 2021; West et al., 2022; Davies et al., 2022). This hints at an initial diverse gut microbiota in newborns due to rapid post-hatching colonization (Trevelline et al., 2018; Pereira et al., 2024). Additionally the decline in microbiota diversity could be linked to decreased parental care (Hubert et al., 1995), resulting in less microbial transmission from parents to nestlings, and to dietary changes as nestlings shift from being fed individual prey pieces by their parents to consuming whole prey items left in the nest (Hoffmann et al., 2013; David et al., 2014).
5 Conclusion
We investigated the interactions of the gut microbiota-immune-brain axis in raptor nestlings. As far as we are aware, this represents one of few studies exploring this axis in a wild vertebrate population and incorporating different time points during the nestling phase in a longitudinal study design. While there was no conclusive evidence for the microbiota-immune-brain axis in our study system, we did find evidence for the hypothesized relationships among stress, Immunity, and body condition. Elevated f-CORT levels were linked to immunosuppression and adverse effects on overall body condition, suggesting that immune and stress regulation play a dominant role in early nestling development. Shannon bacterial microbiota diversity was positively correlated with nestling body condition, suggesting a potential benefit of a diverse and stable gut microbiota. Age plays a crucial role in influencing immune development, stress responses, and eukaryotic microbial diversity. The decline in eukaryotic microbial diversity with age implies an early uncontrolled gut colonization, followed by selective removal of non-relevant taxa. Our study thereby contributes to a growing body of knowledge on the dynamics of the gut microbiota-immune-brain axis in wild populations.
Data availability statement
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-immune-brain-axis-common-buzzard-nestlings.
Ethics statement
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 and 84-02-04.2017.A147. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
HP: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. JH: Conceptualization, Funding acquisition, Investigation, Resources, Supervision, Validation, Writing – review & editing. OK: Writing – review & editing, Funding acquisition, Project administration, Validation, Investigation, Conceptualization, Supervision. GC: Conceptualization, Investigation, Methodology, Resources, Supervision, Writing – review & editing. TR: Data curation, Validation, Writing – review & editing. MO: Data curation, Validation, Writing – review & editing. K-PG: Data curation, Validation, Writing – review & editing. BC: Validation, Writing – review & editing, Methodology. ÖM: Validation, Writing – review & editing, Conceptualization, Methodology. SK: Methodology, Validation, Writing – review & editing, Resources. NC: Conceptualization, Funding acquisition, Investigation, Methodology, Resources, Supervision, Validation, Writing – review & editing, Project administration.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the German Research Foundation (Deutsche Forschungsgemeinschaft—DFG), as part of the SFB TRR 212 (NC3)—project numbers 316099922 and 396780709 and DFG project numbers 398434413, 433069365, and 233740704. The funding body had no involvement in the study design, data collection, analysis, interpretation, or manuscript writing.
Acknowledgments
We express our gratitude to Thomas Grünkorn and Justine Léauté for their dedicated data collection efforts during the field season. Special thanks to Rebecca Chen for her valuable guidance and assistance in statistical analysis. We acknowledge the contributions of Sabine Kruse and Edith Ossendorf in the corticosterone determination process and Katja Pohle during the immunological analyses. We acknowledge the financial support of the German Research Foundation (DFG).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1413976/full#supplementary-material
References
Abt, M. C., and Pamer, E. G. (2014). Commensal bacteria mediated defenses against pathogens. Curr. Opin. Immunol. 29, 16–22. doi: 10.1016/j.coi.2014.03.003
Almasi, B., Roulin, A., Jenni-Eiermann, S., Breuner, C. W., and Jenni, L. (2009). Regulation of free corticosterone and CBG capacity under different environmental conditions in altricial nestlings. Gen. Comp. Endocrinol. 164, 117–124. doi: 10.1016/j.ygcen.2009.05.011
Angelier, F., Wingfield, J. C., Weimerskirch, H., and Chastel, O. (2010). Hormonal correlates of individual quality in a long-lived bird: a test of the “corticosterone-fitness hypothesis”. Biol. Lett. 6, 846–849. doi: 10.1098/rsbl.2010.0376
Apprill, A., McNally, S., Parsons, R., and Weber, L. (2015). Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquatic Microbial Ecology 75, 129–137. doi: 10.3354/ame01753
Bates, D., Mächler, M., Bolker, B., and Walker, S. (2015). Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48. doi: 10.18637/jss.v067.i01
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
Berman, T. S., Weinberg, M., Moreno, K. R., Czirjál, G. Á., and Yovel, Y. (2023). In sickness and in health: the dynamics of the fruit bat gut microbiota under a bacterial antigen challenge and its association with the immune response. Front. Immunol. 14:1152107. doi: 10.3389/fimmu.2023.1152107
Bing, Z., Niu, C., Yang, C., Qi, Y., Zhao, Y., Cao, S., et al. (2022). Beneficial role of gut microbes in maintenance of pace-of-life traits in phrynocephalus vlangalii. Front. Microb. 1:962761. doi: 10.3389/frmbi.2022.962761
Bisgaard, H., Li, N., Bonnelykke, K., Chawes, B. L. K., Skov, T., Paludan-Müller, G., et al. (2011). Reduced diversity of the intestinal microbiota during infancy is associated with increased risk of allergic disease at school age. J. Aller. Clin. Immunol. 128, 646–652. doi: 10.1016/j.jaci.2011.04.060
Bollen, K. A. (2002). Latent variables in psychology and the social sciences. Annu. Rev. Psychol. 53, 605–634. doi: 10.1146/annurev.psych.53.100901.135239
Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37, 852–857. doi: 10.1038/s41587-019-0209-9
Bourgeon, S., Kauffmann, M., Geiger, S., Raclot, T., and Robin, J.-P. (2010). Relationships between metabolic status, corticosterone secretion and maintenance of innate and adaptive humoral immunities in fasted re-fed mallards. J. Exp. Biol. 213, 3810–3818. doi: 10.1242/jeb.045484
Brust, V., Eikenaar, C., Packmor, F., Schmaljohann, H., Hüppop, O., Czirják, G., et al. (2022). Do departure and flight route decisions correlate with immune parameters in migratory songbirds? Funct. Ecol. 36, 3007–3021. doi: 10.1111/1365-2435.14187
Buffie, C. G., and Pamer, E. G. (2013). Microbiota-mediated colonization resistance against intestinal pathogens. Nat. Rev. Immunol. 13, 790–801. doi: 10.1038/nri3535
Bükner, P.-C. (2018). Advanced Bayesian multilevel modeling with the R package brms. R J. 10, 395–411. doi: 10.32614/RJ-2018-017
Bürkner, P.-C. (2017). brms: an R package for bayesian multilevel models using stan. J. Stat. Softw. 80, 1–28. doi: 10.18637/jss.v080.i01
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: High-resolution sample inference from illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Carpenter, B., Gelman, A., Hoffman, M. D., Lee, D., Goodrich, B., Betancourt, M., et al. (2017). Stan: a probabilistic programming language. J. Stat. Softw. 76, 1–32. doi: 10.18637/jss.v076.i01
Castagliuolo, I., Lamont, J. T., Qiu, B., Fleming, S. M., Bhaskar, K., Nikulasson, S., et al. (1996). Acute stress causes mucin release from rat colon: role of corticotropin releasing factor and mast cells. Am. J. Physiol. Gastrointest. Liver Physiol. 271, G884–G892. doi: 10.1152/ajpgi.1996.271.5.G884
Castillo-Ruiz, A., Mosley, M., Jacobs, A. J., Hoffiz, Y. C., and Forger, N. G. (2018). Birth delivery mode alters perinatal cell death in the mouse brain. Proc. Nat. Acad. Sci. U. S. A. 115, 11826–11831. doi: 10.1073/pnas.1811962115
Chin, E. H., Love, O. P., Verspoor, J. J., Williams, T. D., Rowley, K., and Burness, G. (2009). Juveniles exposed to embryonic corticosterone have enhanced flight performance. Proc. R. Soc. B Biol. Sci. 276, 499–505. doi: 10.1098/rspb.2008.1294
Chin, V. K., Yong, V. C., Chong, P. P., Amin Nordin, S., Basir, R., Abdullah, M., et al. (2020). Mycobiome in the gut: a multiperspective review. Mediat. Inflamm. 2020:9560684. doi: 10.1155/2020/9560684
Costello, E. K., Stagaman, K., Dethlefsen, L., Bohannan, B. J., and Relman, D. A. (2012). The application of ecological theory toward an understanding of the human microbiome. Science 336, 1255–1262. doi: 10.1126/science.1224203
Coyte, K. Z., and Rakoff-Nahoum, S. (2019). Understanding competition and cooperation within the mammalian gut microbiome. Curr. Biol. 29, R538–R544. doi: 10.1016/j.cub.2019.04.017
Cryan, J. F., O'Riordan, K. J., Cowan, C. S., Sandhu, K. V., Bastiaanssen, T. F., Boehme, M., et al. (2019). The microbiota-gut-brain axis. Physiol. Rev. 99, 1877–2013. doi: 10.1152/physrev.00018.2018
David, L. A., Maurice, C. F., Carmody, R. N., Gootenberg, D. B., Button, J. E., Wolfe, B. E., et al. (2014). Diet rapidly and reproducibly alters the human gut microbiome. Nature 505, 559–563. doi: 10.1038/nature12820
Davidson, G. L., Raulo, A., and Knowles, S. C. (2020). Identifying microbiome-mediated behaviour in wild vertebrates. Trends Ecol. Evol. 35, 972–980. doi: 10.1016/j.tree.2020.06.014
Davies, C. S., Worsley, S. F., Maher, K. H., Komdeur, J., Burke, T., Dugdale, H. L., et al. (2022). Immunogenetic variation shapes the gut microbiome in a natural vertebrate population. Microbiome 10:41. doi: 10.1186/s40168-022-01233-y
Davis, N. M., Proctor, D. M., Holmes, S. P., Relman, D. A., and Callahan, B. J. (2018). Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data. Microbiome 6, 1–14. doi: 10.1186/s40168-018-0605-2
Dhabhar, F. S. (2018). The short-term stress response-mother nature's mechanism for enhancing protection and performance under conditions of threat, challenge, and opportunity. Front. Neuroendocrinol. 49, 175–192. doi: 10.1016/j.yfrne.2018.03.004
Dhabhar, F. S., and Mcewen, B. S. (1997). Acute stress enhances while chronic stress suppresses cell-mediated immunity in vivo: a potential role for leukocyte trafficking. Brain Behav. Immun. 11, 286–306. doi: 10.1006/brbi.1997.0508
Ding, J.-H., Jin, Z., Yang, X.-X., Lou, J., Shan, W.-X., Hu, Y.-X., et al. (2020). Role of gut microbiota via the gut-liver-brain axis in digestive diseases. World J. Gastroenterol. 26:6141. doi: 10.3748/wjg.v26.i40.6141
Estienne, M., Claustre, J., Clain-Gardechaux, G., Paquet, A., Taché, Y., Fioramonti, J., et al. (2010). Maternal deprivation alters epithelial secretory cell lineages in rat duodenum: role of crf-related peptides. Gut 59, 744–751. doi: 10.1136/gut.2009.190728
Fachi, J. L., de Souza Felipe, J., Pral, L. P., da Silva, B. K., Corrêa, R. O., de Andrade, M. C. P., et al. (2019). Butyrate protects mice from clostridium difficile-induced colitis through an HIF-1-dependent mechanism. Cell Rep. 27, 750–761. doi: 10.1016/j.celrep.2019.03.054
Faith, D. P., and Baker, A. M. (2006). Phylogenetic diversity (PD) and biodiversity conservation: some bioinformatics challenges. Evol. Bioinf. 2:117693430600200007. doi: 10.1177/117693430600200007
Fleischer, R., Jones, C., Ledezma-Campos, P., Czirják, G. Á, Sommer, S., et al. (2024). Gut microbial shifts in vampire bats linked to immunity due to changed diet in human disturbed landscapes. Sci. Total Environ. 907:167815. doi: 10.1016/j.scitotenv.2023.167815
Fletcher, Q. E., Dantzer, B., and Boonstra, R. (2015). The impact of reproduction on the stress axis of free-living male northern red backed voles (myodes rutilus). Gen. Comp. Endocrinol. 224, 136–147. doi: 10.1016/j.ygcen.2015.07.004
Florkowski, M. R., and Yorzinski, J. L. (2023). Gut microbiome diversity and composition is associated with exploratory behavior in a wild-caught songbird. Anim. Microb. 5:8. doi: 10.1186/s42523-023-00227-x
Forthal, D. N. (2014). Functions of antibodies. Microbiol. Spect. 2, 10–1128. doi: 10.1128/microbiolspec.AID-0019-2014
Foster, J. A., and Neufeld, K.-A. M. (2013). Gut-brain axis: how the microbiome influences anxiety and depression. Trends Neurosci. 36, 305–312. doi: 10.1016/j.tins.2013.01.005
Fountain-Jones, N. M., Giraud, T., Zinger, L., Bik, H., Creer, S., and Videvall, E. (2024). Molecular ecology of microbiomes in the wild: common pitfalls, methodological advances and future directions. Mol. Ecol. 33:e17223. doi: 10.1111/mec.17223
Francella, C., Green, M., Caspani, G., Lai, J. K., Rilett, K. C., and Foster, J. A. (2022). Microbe-immune-stress interactions impact behaviour during postnatal development. Int. J. Mol. Sci. 23:15064. doi: 10.3390/ijms232315064
French, S. S., and Neuman-Lee, L. A. (2012). Improved ex vivo method for microbiocidal activity across vertebrate species. Biol. Open 1, 482–487. doi: 10.1242/bio.2012919
Fung, T. C. (2020). The microbiota-immune axis as a central mediator of gut-brain communication. Neurobiol. Dis. 136:104714. doi: 10.1016/j.nbd.2019.104714
Gabor, C. R., Villatoro-Castañeda, M., Carlos-Shanley, C., Ujhegyi, N., and Bókony, V. (2022). Gut bacterial communities vary across habitats and their diversity increases with increasing glucocorticoids in toad tadpoles. Diversity 15:23. doi: 10.3390/d15010023
Gelman, A., Goodrich, B., Gabry, J., and Vehtari, A. (2019). R-squared for bayesian regression models. Am. Stat. 73, 307–309. doi: 10.1080/00031305.2018.1549100
Ghaisas, S., Maher, J., and Kanthasamy, A. (2016). Gut microbiome in health and disease: linking the microbiome-gut-brain axis and environmental factors in the pathogenesis of systemic and neurodegenerative diseases. Pharmacol. Ther. 158, 52–62. doi: 10.1016/j.pharmthera.2015.11.012
Grace, J. B. (2006). Structural Equation Modeling and Natural Systems. Cambridge: Cambridge University Press.
Hanning, I., and Diaz-Sanchez, S. (2015). The functionality of the gastrointestinal microbiome in non-human animals. Microbiome 3, 1–11. doi: 10.1186/s40168-015-0113-6
Hanssen, S. A., Bustnes, J. O., Schnug, L., Bourgeon, S., Johnsen, T. V., Ballesteros, M., et al. (2013). Antiparasite treatments reduce humoral immunity and impact oxidative status in raptor nestlings. Ecol. Evol. 3, 5157–5166. doi: 10.1002/ece3.891
Hayashi, I., Fujita, H., and Toju, H. (2024). Deterministic and stochastic processes generating alternative states of microbiomes. ISME Commun. 4:ycae007. doi: 10.1093/ismeco/ycae007
Hayward, L. S., and Wingfield, J. C. (2004). Maternal corticosterone is transferred to avian yolk and may alter offspring growth and adult phenotype. Gen. Comp. Endocrinol. 135, 365–371. doi: 10.1016/j.ygcen.2003.11.002
Heiman, M. L., and Greenway, F. L. (2016). A healthy gastrointestinal microbiome is dependent on dietary diversity. Mol. Metab. 5, 317–320. doi: 10.1016/j.molmet.2016.02.005
Henriksen, R., Rettenbacher, S., and Groothuis, T. G. (2011). Prenatal stress in birds: pathways, effects, function and perspectives. Neurosci. Biobehav. Rev. 35, 1484–1501. doi: 10.1016/j.neubiorev.2011.04.010
Hird, S. M. (2017). Evolutionary biology needs wild microbiomes. Front. Microbiol. 8:725. doi: 10.3389/fmicb.2017.00725
Hoffmann, C., Dollive, S., Grunberg, S., Chen, J., Li, H., Wu, G. D., et al. (2013). Archaea and fungi of the human gut microbiome: correlations with diet and bacterial residents. PLoS ONE 8:e66019. doi: 10.1371/journal.pone.0066019
Hooper, L. V., Littman, D. R., and Macpherson, A. J. (2012). Interactions between the microbiota and the immune system. Science 336, 1268–1273. doi: 10.1126/science.1223490
Hõrak, P., Saks, L., Ots, I., and Kollist, H. (2002). Repeatability of condition indices in captive greenfinches (carduelis chloris). Can. J. Zool. 80, 636–643. doi: 10.1139/z02-038
Hõrak, P., Saks, L., Ots, I., Kullissaar, T., Kollist, H., and Zilmer, M. (2003). Physiological effects of immune challenge in captive greenfinches (carduelis chloris). Can. J. Zool. 81, 371–379. doi: 10.1139/z03-020
Hu, L.-T., and Bentler, P. M. (1999). Cutoff criteria for fit indexes in covariance structure analysis: conventional criteria versus new alternatives. Struct. Eq. Model. 6, 1–55. doi: 10.1080/10705519909540118
Irene Tieleman, B., Williams, J. B., Ricklefs, R. E., and Klasing, K. C. (2005). Constitutive innate immunity is a component of the pace-of-life syndrome in tropical birds. Proc. R. Soc. B Biol. Sci. 272, 1715–1720. doi: 10.1098/rspb.2005.3155
Janiak, M. C., Montague, M. J., Villamil, C. I., Stock, M. K., Trujillo, A. E., DePasquale, A. N., et al. (2021). Age and sex-associated variation in the multi-site microbiome of an entire social group of free-ranging rhesus macaques. Microbiome 9:1–17. doi: 10.1186/s40168-021-01009-w
Jašarevi ć, E., and Bale, T. L. (2019). Prenatal and postnatal contributions of the maternal microbiome on offspring programming. Front. Neuroendocrinol. 55:100797. doi: 10.1016/j.yfrne.2019.100797
Jenni-Eiermann, S., Helfenstein, F., Vallat, A., Glauser, G., and Jenni, L. (2015). Corticosterone: effects on feather quality and deposition into feathers. Methods Ecol. Evol. 6, 237–246. doi: 10.1111/2041-210X.12314
Johnson, K. V.-A., and Burnet, P. W. (2016). Microbiome: should we diversify from diversity, Gut Microb. 7, 455–458. doi: 10.1080/19490976.2016.1241933
Jones, E. W., Carlson, J. M., Sivak, D. A., and Ludington, W. B. (2022). Stochastic microbiome assembly depends on context. Proc. Nat. Acad. Sci. U. S. A. 119:e2115877119. doi: 10.1073/pnas.2115877119
Jonker, R. M., Chakarov, N., and Krüger, O. (2014). Climate change and habitat heterogeneity drive a population increase in common b uzzards b uteo buteo through effects on survival. Ibis 156, 97–106. doi: 10.1111/ibi.12124
Kaplan, D., and Depaoli, S. (2012). “Bayesian structural equation modeling,” in Handbook of Structural Equation Modeling, Chapter 18, ed. R. H. Hoyle (New York, NY: The Guilford Press), 433–454.
Katoh, K., Misawa, K., Kuma, K.-I., and Miyata, T. (2002). MAFFT: a novel method for rapid multiple sequence alignment based on fast fourier transform. Nucleic Acids Res. 30, 3059–3066. doi: 10.1093/nar/gkf436
Klasing, K., and Leshchinsky, T. (1999). Functions, costs, and benefits of the immune system during development and growth. Ostrich 69:32.
Kohl, K. D., Brun, A., Bordenstein, S. R., Caviedes-Vidal, E., and Karasov, W. H. (2018). Gut microbes limit growth in house sparrow nestlings (passer domesticus) but not through limitations in digestive capacity. Integr. Zool. 13, 139–151. doi: 10.1111/1749-4877.12289
Kostrzewa, A., and Kostrzewa, R. (1990). The relationship of spring and summer weather with density and breeding performance of the buzzard buteo buteo, goshawk accipiter gentilis and kestrel falco tinnunculus. Ibis 132, 550–559. doi: 10.1111/j.1474-919X.1990.tb00278.x
Kreisinger, J., Schmiedová, L., Petržzelková, A., Tomášek, O., Adámková, M., Michálková, R., et al. (2018). Fecal microbiota associated with phytohaemagglutinin-induced immune response in nestlings of a passerine bird. Ecol. Evol. 8, 9793–9802. doi: 10.1002/ece3.4454
Krishnamurthy, H. K., Pereira, M., Bosco, J., and Wang, T. (2023). Gut commensals and their metabolites in health and disease. Front. Microbiol. 14:1244293. doi: 10.3389/fmicb.2023.1244293
Krüger, O. (2004). The importance of competition, food, habitat, weather and phenotype for the reproduction of buzzard buteo buteo. Bird Study 51, 125–132. doi: 10.1080/00063650409461344
Le Chatelier, E., Nielsen, T., Qin, J., Prifti, E., Hildebrand, F., Falony, G., et al. (2013). Richness of human gut microbiome correlates with metabolic markers. Nature 500, 541–546. doi: 10.1038/nature12506
Lefcheck, J. S. (2016). piecewisesem: piecewise structural equation modelling in r for ecology, evolution, and systematics. Methods Ecol. Evol. 7, 573–579. doi: 10.1111/2041-210X.12512
Lin, H., Decuypere, E., and Buyse, J. (2006). Acute heat stress induces oxidative stress in broiler chickens. Comp. Biochem. Physiol. Part A 144, 11–17. doi: 10.1016/j.cbpa.2006.01.032
Lin, H., Eggesbo, M., and Peddada, S. D. (2022). Linear and nonlinear correlation estimators unveil undescribed taxa interactions in microbiome data. Nat. Commun. 13, 1–16. doi: 10.1038/s41467-022-32243-x
Lin, H., and Peddada, S. D. (2020). Analysis of compositions of microbiomes with bias correction. Nat. Commun. 11, 1–11. doi: 10.1038/s41467-020-17041-7
Lo, B. C., Chen, G. Y., Núñez, G., and Caruso, R. (2021). Gut microbiota and systemic immunity in health and disease. Int. Immunol. 33, 197–209. doi: 10.1093/intimm/dxaa079
Lochmiller, R. L., and Deerenberg, C. (2000). Trade-offs in evolutionary immunology: just what is the cost of immunity. Oikos 88, 87–98. doi: 10.1034/j.1600-0706.2000.880110.x
Long-Smith, C., O'Riordan, K. J., Clarke, G., Stanton, C., Dinan, T. G., and Cryan, J. F. (2020). Microbiota-gut-brain axis: new therapeutic opportunities. Annu. Rev. Pharmacol. Toxicol. 60, 477–502. doi: 10.1146/annurev-pharmtox-010919-023628
Love, O., Wynne-Edwards, K., Bond, L., and Williams, T. (2008). Determinants of within-and among-clutch variation in yolk corticosterone in the european starling. Horm. Behav. 53, 104–111. doi: 10.1016/j.yhbeh.2007.09.007
Lozupone, C. A., Stombaugh, J. I., Gordon, J. I., Jansson, J. K., and Knight, R. (2012). Diversity, stability and resilience of the human gut microbiota. Nature 489, 220–230. doi: 10.1038/nature11550
Lynch, C. M. K., O'Riordan, K. J., Clarke, G., and Cryan, J. F. (2024). “Gut Microbes: The gut brain connection,” in Clinical Understanding of the Human Gut Microbiome, eds. M. Pimentel, R. Mathur, and G. M. Barlow (Cham: Springer). doi: 10.1007/978-3-031-46712-7_4
MacCallum, R. C., Browne, M. W., and Sugawara, H. M. (1996). Power analysis and determination of sample size for covariance structure modeling. Psychol. Methods 1:130. doi: 10.1037/1082-989X.1.2.130
Machida, R. J., and Knowlton, N. (2012). PCR primers for metazoan nuclear 18S and 28S ribosomal dna sequences. PLoS ONE 7:e46180. doi: 10.1371/journal.pone.0046180
Madden, A. A., Oliverio, A. M., Kearns, P. J., Henley, J. B., Fierer, N., Starks, P. T., et al. (2022). Chronic stress and captivity alter the cloacal microbiome of a wild songbird. J. Exp. Biol. 225:jeb243176. doi: 10.1242/jeb.243176
Martin, C. R., Osadchiy, V., Kalani, A., and Mayer, E. A. (2018). The brain-gut-microbiome axis. Cell. Mol. Gastroenterol. Hepatol. 6, 133–148. doi: 10.1016/j.jcmgh.2018.04.003
Martin, L. B. (2009). Stress and immunity in wild vertebrates: timing is everything. Gen. Comp. Endocrinol. 163, 70–76. doi: 10.1016/j.ygcen.2009.03.008
Matson, K. D., Horrocks, N. P., Versteegh, M. A., and Tieleman, B. I. (2012). Baseline haptoglobin concentrations are repeatable and predictive of certain aspects of a subsequent experimentally-induced inflammatory response. Comp. Biochem. Physiol. 162, 7–15. doi: 10.1016/j.cbpa.2012.01.010
Matson, K. D., Ricklefs, R. E., and Klasing, K. C. (2005). A hemolysis-hemagglutination assay for characterizing constitutive innate humoral immunity in wild and domestic birds. Dev. Comp. Immunol. 29, 275–286. doi: 10.1016/j.dci.2004.07.006
Mikkelsen, A. J., Lesmeister, D. B., O'Reilly, K. M., and Dugger, K. M. (2023). Juvenile northern spotted owls with higher mass and intermediate levels of corticosterone have greater long-term survival. Ornithol. Appl. 125:duad015. doi: 10.1093/ornithapp/duad015
Millet, S., Bennett, J., Lee, K. A., Hau, M., and Klasing, K. C. (2007). Quantifying and comparing constitutive immunity across avian species. Dev. Comp. Immunol. 31, 188–201. doi: 10.1016/j.dci.2006.05.013
Mindus, C., Ellis, J., Van Staaveren, N., and Harlander-Matauschek, A. (2021). Lactobacillus-based probiotics reduce the adverse effects of stress in rodents: a meta-analysis. Front. Behav. Neurosci. 15:642757. doi: 10.3389/fnbeh.2021.642757
Mohring, B., Öst, M., Jaatinen, K., Parenteau, C., Grimaud, E., and Angelier, F. (2023). Breeding under pressure: corticosterone is associated with reproductive investment under fluctuating predation risk in a long-lived sea duck. Funct. Ecol. 37, 2868–2882. doi: 10.1111/1365-2435.14435
Morais, L. H., Golubeva, A. V., Moloney, G. M., Moya-Pérez, A., Ventura-Silva, A. P., Arboleya, S., et al. (2020). Enduring behavioral effects induced by birth by caesarean section in the mouse. Curr. Biol. 30, 3761–3774. doi: 10.1016/j.cub.2020.07.044
Nauta, A. J., Roos, A., and Daha, M. R. (2004). A regulatory role for complement in innate immunity and autoimmunity. Int. Arch. Allergy Immunol. 134, 310–323. doi: 10.1159/000079261
Nebel, C., Amar, A., Hegemann, A., Isaksson, C., and Sumasgutner, P. (2021). Parental morph combination does not influence innate immune function in nestlings of a colour-polymorphic african raptor. Sci. Rep. 11:11053. doi: 10.1038/s41598-021-90291-7
Noguera, J. C., Aira, M., Pérez-Losada, M., Domínguez, J., and Velando, A. (2018). Glucocorticoids modulate gastrointestinal microbiome in a wild bird. R. Soc. Open Sci. 5:171743. doi: 10.1098/rsos.171743
Ochsenbein, A. F., Fehr, T., Lutz, C., Suter, M., Brombacher, F., Hengartner, H., et al. (1999). Control of early viral and bacterial distribution and disease by natural antibodies. Science 286, 2156–2159. doi: 10.1126/science.286.5447.2156
Oever, J. T., and Netea, M. G. (2014). The bacteriome-mycobiome interaction and antifungal host defense. Eur. J. Immunol. 44, 3182–3191. doi: 10.1002/eji.201344405
Oldereid, T. S., Jiang, X., Nordhus, K. S., Ponzetta, A., Bjørnholt, J. V., Björkström, N. K., et al. (2023). Role of bacteria and microbial metabolites in immune modulation during early life. Scand. J. Immunol. 99:e13336. doi: 10.1111/sji.13336
Palacios, M. G., Cunnick, J. E., Vleck, D., and Vleck, C. M. (2009). Ontogeny of innate and adaptive immune defense components in free-living tree swallows, tachycineta bicolor. Dev. Comp. Immunol. 33, 456–463. doi: 10.1016/j.dci.2008.09.006
Panek, M. (2021). Does habitat diversity modify the dietary and reproductive response to prey fluctuations in a generalist raptor predator, the eurasian buzzard buteo buteo? Birds 2, 114–126. doi: 10.3390/birds2010008
Parada, A. E., Needham, D. M., and Fuhrman, J. A. (2016). Every base matters: assessing small subunit rrna primers for marine microbiomes with mock communities, time series and global field samples. Environ. Microbiol. 18, 1403–1414. doi: 10.1111/1462-2920.13023
Pereira, H., Chakarov, N., Hoffman, J. I., Rinaud, T., Ottensmann, M., Gladow, K.-P., et al. (2024). Early-life factors shaping the gut microbiota of common buzzard nestlings. Anim. Microb. 6:27. doi: 10.1186/s42523-024-00313-8
Pereira, K. E., Bletz, M. C., McCartney, J. A., Woodhams, D. C., and Woodley, S. K. (2023). Effects of exogenous elevation of corticosterone on immunity and the skin microbiome of eastern newts (notophthalmus viridescens). Philos. Transact. R. Soc. B 378:20220120. doi: 10.1098/rstb.2022.0120
Petrullo, L., Ren, T., Wu, M., Boonstra, R., Palme, R., Boutin, S., et al. (2022). Glucocorticoids coordinate changes in gut microbiome composition in wild north american red squirrels. Sci. Rep. 12:2605. doi: 10.1038/s41598-022-06359-5
Price, M. N., Dehal, P. S., and Arkin, A. P. (2010). Fasttree 2-approximately maximum-likelihood trees for large alignments. PLoS ONE 5:e9490. doi: 10.1371/journal.pone.0009490
Prüter, H., Franz, M., Twietmeyer, S., Böhm, N., Middendorff, G., Portas, R., et al. (2020). Increased immune marker variance in a population of invasive birds. Sci. Rep. 10:21764. doi: 10.1038/s41598-020-78427-7
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2012). The silva ribosomal rna gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Quirici, V., Botero-Delgadillo, E., González-Gómez, P. L., Espíndola-Hernández, P., Zambrano, B., Cuevas, E., et al. (2021). On the relationship between baseline corticosterone levels and annual survival of the thorn-tailed rayadito. Gen. Comp. Endocrinol. 300:113635. doi: 10.1016/j.ygcen.2020.113635
R Core Team (2022). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Ratsika, A., Pereira, J. S. C., Lynch, C. M., Clarke, G., and Cryan, J. F. (2023). Microbiota-immune-brain interactions: a lifespan perspective. Curr. Opin. Neurobiol. 78:102652. doi: 10.1016/j.conb.2022.102652
Réale, D., Garant, D., Humphries, M. M., Bergeron, P., Careau, V., and Montiglio, P.-O. (2010). Personality and the emergence of the pace-of-life syndrome concept at the population level. Philos. Transact. R. Soc. B 365, 4051–4063. doi: 10.1098/rstb.2010.0208
Reese, A. T., and Dunn, R. R. (2018). Drivers of microbiome biodiversity: a review of general rules, feces, and ignorance. MBio 9, 10–1128. doi: 10.1128/mBio.01294-18
Rhee, S. H., Pothoulakis, C., and Mayer, E. A. (2009). Principles and clinical implications of the brain-gut-enteric microbiota axis. Nat. Rev. Gastroenterol. Hepatol. 6, 306–314. doi: 10.1038/nrgastro.2009.35
Roduit, C., Scholtens, S., de Jongste, J. C., Wijga, A. H., Gerritsen, J., Postma, D. S., et al. (2009). Asthma at 8 years of age in children born by caesarean section. Thorax 64, 107–113. doi: 10.1136/thx.2008.100875
Romero, L. M., Dickens, M. J., and Cyr, N. E. (2009). The reactive scope model—A new model integrating homeostasis, allostasis, and stress. Horm. Behav. 55, 375–389. doi: 10.1016/j.yhbeh.2008.12.009
Romero, L. M., and Fairhurst, G. D. (2016). Measuring corticosterone in feathers: strengths, limitations, and suggestions for the future. Comp. Biochem. Physiol. Part A 202, 112–122. doi: 10.1016/j.cbpa.2016.05.002
Romero, L. M., and Wikelski, M. (2001). Corticosterone levels predict survival probabilities of galapagos marine iguanas during el nino events. Proc. Nat. Acad. Sci. U. S. A. 98, 7366–7370. doi: 10.1073/pnas.131091498
Rosseel, Y. (2012). lavaan: An R package for structural equation modeling. J. Stat. Softw. 48, 1–36. doi: 10.18637/jss.v048.i02
Round, J. L., and Mazmanian, S. K. (2009). The gut microbiota shapes intestinal immune responses during health and disease. Nat. Rev. Immunol. 9, 313–323. doi: 10.1038/nri2515
Saino, N., Romano, M., Ferrari, R. P., Martinelli, R., and Möller, A. P. (2005). Stressed mothers lay eggs with high corticosterone levels which produce low-quality offspring. J. Exp. Zool. Part A 303, 998–1006. doi: 10.1002/jez.a.224
Sapolsky, R. M., Romero, L. M., and Munck, A. U. (2000). How do glucocorticoids influence stress responses? Integrating permissive, suppressive, stimulatory, and preparative actions. Endocr. Rev. 21, 55–89. doi: 10.1210/edrv.21.1.0389
Shade, A. (2017). Diversity is the question, not the answer. ISME J. 11, 1–6. doi: 10.1038/ismej.2016.118
Shade, A., Jones, S. E., Caporaso, J. G., Handelsman, J., Knight, R., Fierer, N., et al. (2014). Conditionally rare taxa disproportionately contribute to temporal changes in microbial diversity. MBio 5, 10–1128. doi: 10.1128/mBio.01371-14
Shannon, C. E. (1948). A mathematical theory of communication. Bell Syst. Tech. J. 27, 379–423. doi: 10.1002/j.1538-7305.1948.tb01338.x
Sharma, S., Mukherjee, S., Kumar, A., and Dillon, W. R. (2005). A simulation study to investigate the use of cutoff values for assessing model fit in covariance structure models. J. Bus. Res. 58, 935–943. doi: 10.1016/j.jbusres.2003.10.007
Shigeshiro, M., Tanabe, S., and Suzuki, T. (2012). Repeated exposure to water immersion stress reduces the MUC2 gene level in the rat colon via two distinct mechanisms. Brain Behav. Immun. 26, 1061–1065. doi: 10.1016/j.bbi.2012.05.016
Slevin, M. C., Houtz, J. L., Bradshaw, D. J., and Anderson, R. C. (2020). Evidence supporting the microbiota-gut-brain axis in a songbird. Biol. Lett. 16:20200430. doi: 10.1098/rsbl.2020.0430
Smulders, T. V. (2021). Telencephalic regulation of the hpa axis in birds. Neurobiol. Stress 15:100351. doi: 10.1016/j.ynstr.2021.100351
Spencer, K. A., Evans, N. P., and Monaghan, P. (2009). Postnatal stress in birds: a novel model of glucocorticoid programming of the hypothalamic-pituitary-adrenal axis. Endocrinology 150, 1931–1934. doi: 10.1210/en.2008-1471
Sperringer, J. E., Addington, A., and Hutson, S. M. (2017). Branched-chain amino acids and brain metabolism. Neurochem. Res. 42:1697–1709. doi: 10.1007/s11064-017-2261-5
Stoffel, M. A., Acevedo-Whitehouse, K., Morales-Durán, N., Grosser, S., Chakarov, N., Krüger, O., et al. (2020). Early sexual dimorphism in the developing gut microbiome of northern elephant seals. Mol. Ecol. 29, 2109–2122. doi: 10.1111/mec.15385
Stothart, M. R., Palme, R., and Newman, A. E. (2019). It's what's on the inside that counts: stress physiology and the bacterial microbiome of a wild urban mammal. Proc. R. Soc. B 286:20192111. doi: 10.1098/rspb.2019.2111
Strange, M. S., Bowden, R. M., Thompson, C. F., and Sakaluk, S. K. (2016). Pre-and postnatal effects of corticosterone on fitness-related traits and the timing of endogenous corticosterone production in a songbird. J. Exp. Zool. Part A 325, 347–359. doi: 10.1002/jez.2022
Sudo, N., Chida, Y., Aiba, Y., Sonoda, J., Oyama, N., Yu, X.-N., et al. (2004). Postnatal microbial colonization programs the hypothalamic-pituitary-adrenal system for stress response in mice. J. Physiol. 558, 263–275. doi: 10.1113/jphysiol.2004.063388
Sweeny, A. R., Lemon, H., Ibrahim, A., Watt, K. A., Wilson, K., Childs, D. Z., et al. (2023). A mixed-model approach for estimating drivers of microbiota community composition and differential taxonomic abundance. Msystems 8, e00040–e00023. doi: 10.1128/msystems.00040-23
Sylvia, K. E., and Demas, G. E. (2018). A gut feeling: microbiome-brain-immune interactions modulate social and affective behaviors. Horm. Behav. 99, 41–49. doi: 10.1016/j.yhbeh.2018.02.001
Tieleman, B. I. (2018). Understanding immune function as a pace of life trait requires environmental context. Behav. Ecol. Sociobiol. 72:55. doi: 10.1007/s00265-018-2464-z
Trevelline, B. K., MacLeod, K. J., Knutie, S. A., Langkilde, T., and Kohl, K. D. (2018). In ovo microbial communities: a potential mechanism for the initial acquisition of gut microbiota among oviparous birds and lizards. Biol. Lett. 14:20180225. doi: 10.1098/rsbl.2018.0225
Vagasi, C. I., Pătras, L., Pap, P. L., Vincze, O., Muresan, C., Nemeth, J., et al. (2018). Experimental increase in baseline corticosterone level reduces oxidative damage and enhances innate immune response. PLoS ONE 13:e0192701. doi: 10.1371/journal.pone.0192701
Van der Auwera, G., Chapelle, S., and De Wächter, R. (1994). Structure of the large ribosomal subunit rna of phytophthora megasperma, and phylogeny of the oomycetes. FEBS Lett. 338, 133–136. doi: 10.1016/0014-5793(94)80350-1
van Veelen, H. P. J., Falcão Salles, J., Matson, K. D., van der Velde, M., and Tieleman, B. I. (2020). Microbial environment shapes immune function and cloacal microbiota dynamics in zebra finches taeniopygia guttata. Anim. Microb. 2, 1–17. doi: 10.1186/s42523-020-00039-3
Vega, N. M., and Gore, J. (2017). Stochastic assembly produces heterogeneous communities in the caenorhabditis elegans intestine. PLoS Biol. 15:e2000633. doi: 10.1371/journal.pbio.2000633
Vemuri, R., Shankar, E. M., Chieppa, M., Eri, R., and Kavanagh, K. (2020). Beyond just bacteria: functional biomes in the gut ecosystem including virome, mycobiome, archaeome and helminths. Microorganisms 8:483. doi: 10.3390/microorganisms8040483
Vincze, O., Vágási, C. I., Pénzes, J., Szabó, K., Magonyi, N. M., Czirják, G., et al. (2022). Sexual dimorphism in immune function and oxidative physiology across birds: the role of sexual selection. Ecol. Lett. 25, 958–970. doi: 10.1111/ele.13973
Warr, G. W., Magor, K. E., and Higgins, D. A. (1995). IGY: clues to the origins of modern antibodies. Immunol. Today 16, 392–398. doi: 10.1016/0167-5699(95)80008-5
West, A. G., Digby, A., Lear, G., Team, K. R., and Taylor, M. W. (2022). 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. Anim. Microb. 4:55. doi: 10.1186/s42523-022-00204-w
Keywords: eukaryotic microbiota, bacterial microbiota, stress, immune system, body condition, longitudinal study, Buteo buteo
Citation: Pereira H, Hoffman JI, Krüger O, Czirják GÁ, Rinaud T, Ottensmann M, Gladow K-P, Caspers BA, Maraci Ö, Kaiser S and Chakarov N (2024) The gut microbiota-immune-brain axis in a wild vertebrate: dynamic interactions and health impacts. Front. Microbiol. 15:1413976. doi: 10.3389/fmicb.2024.1413976
Received: 08 April 2024; Accepted: 20 August 2024;
Published: 10 September 2024.
Edited by:
Lifeng Zhu, Nanjing University of Chinese Medicine, ChinaReviewed by:
Mason Stothart, University of Oxford, United KingdomLeigh Combrink, University of Arizona, United States
Copyright © 2024 Pereira, Hoffman, Krüger, Czirják, Rinaud, Ottensmann, Gladow, Caspers, Maraci, Kaiser and Chakarov. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hugo Pereira, hugo.eira@uni-bielefeld.de
†These authors have contributed equally to this work