Next Article in Journal
Oxidative Stress and Placental Pathogenesis: A Contemporary Overview of Potential Biomarkers and Emerging Therapeutics
Previous Article in Journal
Association Between Zinc Status and Insulin Resistance/Sensitivity Check Indexes in Gestational Diabetes Mellitus
Previous Article in Special Issue
Molecular Chaperonin HSP60: Current Understanding and Future Prospects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prenatal Acoustic Signals Influence Nestling Heat Shock Protein Response to Heat and Heterophil-to-Lymphocyte Ratio in a Desert Bird

1
School of Life & Environmental Sciences, Deakin University, 75 Pigdons Road, Geelong, VIC 3216, Australia
2
Max Planck Institute for Biological Intelligence, Eberhard-Gwinner-Straße, 82319 Seewiesen, Germany
3
Centre d’Écologie et des Sciences de la Conservation (CESCO), Muséum National d’Histoire Naturelle, Centre National de la Recherche Scientifique, Sorbonne Université, CP 135, 57 rue Cuvier, 75005 Paris, France
4
College of Science and Engineering, Flinders University, Adelaide, SA 5042, Australia
5
Doñana Biological Station EBC-CSIC, Calle Américo-Vespucio 26, 41092 Sevilla, Spain
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(22), 12194; https://doi.org/10.3390/ijms252212194
Submission received: 14 September 2024 / Revised: 9 November 2024 / Accepted: 11 November 2024 / Published: 13 November 2024
(This article belongs to the Special Issue Advances in Heat-Shock Response and Heat-Shock Proteins)

Abstract

:
Heat shock proteins (HSPs) are essential to cellular protection against heat stress. However, the causes of inter-individual variation in HSP regulation remain unclear. This study aimed to test the impact of early-life conditions on the HSP response to heat in zebra finches. In this arid-adapted bird, incubating parents emit “heat-calls” at high temperatures, which adaptively alter offspring’s phenotypes. Embryos were exposed to heat-calls or control-calls, and at 13 days post-hatch nestlings were separated into two different experiments to test responses to either chronic nest temperature (“in-nest” experiment) or an acute “heat-challenge”. Blood samples were collected to measure levels of heat shock cognate 70, heat shock protein 90α, corticosterone and the heterophil-to-lymphocyte (H/L) ratio. In the in-nest experiment, both HSPs were upregulated in response to increasing nest temperatures only in control-calls nestlings (HSC70: p = 0.010, HSP90α: p = 0.050), which also had a marginally higher H/L ratio overall than heat-call birds (p = 0.066). These results point to a higher heat sensitivity in control-call nestlings. Furthermore, comparing across experiments, only the H/L ratio differed, being higher in heat-challenged than in in-nest nestlings (p = 0.009). Overall, this study shows for the first time that a prenatal acoustic signal of heat affects the nestling HSP response to postnatal temperature.

1. Introduction

Animals respond to environmental fluctuations by adjusting their physiology to current conditions. In particular, environmental challenges such as changes in weather, predator presence or food scarcity trigger a complex and multi-level “stress response”, which allows animals to maintain homeostasis in the face of such stressors [1,2]. Understanding organism responses to stressors is essential, because these responses define the fitness consequences of environmental disturbances, and thus ultimately the impact of disturbances on population persistence [3,4]. However, it remains unclear how individuals vary in their stress responses, from the cellular to organismal levels, and what factors may contribute to such inter-individual variation [2,5].
Within cells, heat shock proteins (HSPs) are essential molecular chaperones that play important functions to maintain proteostasis (i.e., protein homeostasis) and, therefore, cell survival under stress [6,7]. These highly conserved molecules respond to multiple stressors, including heat stress [8,9,10]. In particular, HSP70 and HSP90 families are well studied for their role in thermotolerance [11,12,13]. The stress-inducible HSPs (e.g., HSP70, HSP90α) are upregulated under acute stress [12,13], while constitutive HSPs (e.g., heat shock cognate 70, HSC70) are always expressed but modulated by long-term environmental conditions, such as thermal acclimatization (reviews: [4,10]; studies: [14,15]). Overall, however, compared to the wealth of studies on the molecular mechanisms of HSP regulation in vitro or under controlled laboratory conditions, much less is known about HSP regulation in an ecological context and its role in adaptation to environmental stressors [2,16].
At the whole-organism level, acute environmental stressors trigger the so-called “stress response”, which is mediated by the activation of the hypothalamic–pituitary–adrenal (HPA) axis and the subsequent release of glucocorticoid (GC) hormones into the peripheral blood within minutes [17,18]. These hormones bind primarily to glucocorticoid receptors (GR) within cells, whose maturation and functioning involves the participation of HSP70 and HSP90, respectively [19,20]. This suggests HSPs may also modulate the impact of GC release on the organism, at least in mammals [19,20]. Furthermore, across taxa, when exposure to a stressor or elevated GC levels becomes chronic, the relative proportion of circulating white blood cells in peripheral blood is altered [21,22]. Namely, heterophil numbers (in birds, or neutrophil in mammals) increase while lymphocyte numbers decrease, resulting in an elevation of the heterophil-to-lymphocyte (H/L) ratio [23,24]. The H/L ratio is a robust biomarker of medium-to-long-term stress because its time course is slower, steadier and longer lasting than the acute GC response [25,26]. Considering these organismal responses may help with interpreting inter-individual differences in HSP regulation.
Variability in HSP responses is observed between species [4,13], and between populations that differ in environmental harshness or predictability (e.g., [14,27]). However, less is known about how much and why HSP responses vary between individuals. Heritability of the HSP response is moderate (e.g., turtle: [28]; drosophila: [29]; human: [30]) and therefore only explains part of inter-individual variation. Early-life conditions, based on their well-known impact (especially for maternal GCs) on the organismal stress response and HPA axis regulation across taxa (including humans) [31,32], may also contribute to HSP variation. Indeed, at least in ectotherms, maternal HSP levels have been shown to increase offspring HSP levels and/or heat tolerance, (e.g., drosophila: [33]; zebrafish: [34]). In addition, incubation temperature in oviparous species may shape the HSP response (already functional during embryogenesis [28,35]). Specifically in chicken, controlled heat exposure during incubation can increase individual HSP levels during embryogenesis, but later decrease HSP levels and oxidative stress post-hatching during a heat-challenge [36]. However, exposure to high temperatures at the juvenile stage in zebra finches (Taeniopygia guttata castanotis) had a limited effect on HSP responses in adulthood [37]. Overall, therefore, while there is evidence that it can happen, the impact of early life experience on HSP responses remains little understood.
Among all organisms, birds are particularly exposed to high thermal extremes from being diurnal and not using underground thermal refuges [38,39]. Although the great majority of avian research focuses on poultry, there is evidence for the impact of high temperatures on HSP regulation (poultry: [11,12]; other birds: [13,40]), CORT levels [24,41] and the H/L ratio [42,43]. Nonetheless, high heat exposure in birds may also have favored specific adaptations to withstand heat, including at the cellular level [2]. This is the case in the zebra finch, a desert-adapted species, in which a prenatal acoustic signal of heat programs offspring development for high temperatures [44]. Zebra finch parents produce specific “heat-calls” during incubation at high temperatures, which adaptively shapes offspring growth to hot conditions [44], as well as nestling cellular metabolism and its response to temperature [45]. These changes led in adulthood to higher reproductive success [44] and generally lower chronic stress levels (measured as the H/L ratio) [46]. Prenatal heat-call exposure also led to higher heat tolerance in adulthood, although not through improved evaporative cooling efficiency [47]. Whether changes in HSP regulation could underlie some of these effects remains to be determined.
Therefore, this study tested the hypothesis that early life conditions in a wild bird species influence HSP regulation later in life, in a way that may reflect changes towards heat adaption. Specifically, it tested whether prenatal exposure to heat-calls in the zebra finch affects individual HSP responses to chronic and acute thermal conditions, in parallel to other components of nestling stress physiology (CORT and the H/L ratio). To this aim, embryos were exposed to a prenatal playback of either heat-calls (treatment) or control-calls (control), and after hatching postnatal nest temperatures were manipulated. At 13 days post-hatch, the impact of early acoustic experience on responses to either chronic or acute thermal conditions were investigated in two separate experiments. Specifically, (i) in the “in-nest” experiment, nestlings were subjected to manipulated chronic nest temperatures (average daytime temperature from hatching to day 12, 12D-Tnest, ranging from 23 to 32 °C), or, (ii) in the acute heat-challenge experiment, to a 2.5 h gradual temperature increase up to 44 °C for 20 min in a temperature-controlled chamber (Figure 1). Blood samples were then collected within three minutes of nest disturbance or at the end of the heat-challenge, to measure protein levels of the constitutive heat shock cognate 70 (HSC70) and stress-inducible HSP90α (in red blood cells), baseline CORT (in plasma) and the H/L ratio (from blood smears). We predicted that if prenatal heat-calls affect subsequent HSP regulation and prepare individuals for high temperatures in postnatal life, then in the in-nest experiment (i), a higher 12D-Tnest will be chronically challenging for control nestlings, but not for heat-call nestlings. Thus, an increase in both HSPs and both stress biomarkers with 12D-Tnest is expected in control-call birds only. In addition, in the heat-challenge experiment (ii), we predicted that control-call nestlings will be more affected by the heat-challenge than heat-call nestlings, resulting in higher stress-inducible HSP90α levels and H/L ratios (but not higher HSC70 [constitutively expressed] nor CORT [back to baseline levels after acute response: [25]]). Then, (iii) comparing across experiments, we expected that heat-challenged nestlings may have higher HSP90α levels and H/L ratios than in-nest nestlings, but similar levels of HSC70 and CORT. Lastly, (iv) given that HSP regulation may affect growth [9], we assessed how nestling mass relates to the stress biomarkers.

2. Results

2.1. Nestling Response to Nest Temperature: In-Nest Experiment

Under “in-nest” conditions, there was a significant interaction between the prenatal playback and the acclimatization or chronic temperature experienced by nestlings in the nest (i.e., 12D-Tnest, averaged from hatching to day 12) on heat shock protein levels (HSC70: p = 0.010, HSP90α: p = 0.050; Table 1, Figure 2). Specifically, levels of both the constitutive HSC70 and stress-inducible HSP90α increased with 12D-Tnest in control-call nestlings, whereas in heat-call nestlings HSC70 levels decreased slightly with 12D-Tnest, and HSP90α levels did not vary (Figure 2). This conforms to our prediction that hotter nests are more challenging for control-call birds and trigger an HSP upregulation.
The H/L ratio in in-nest nestlings was also affected by the prenatal playback, but not by nest temperature (Table 1). Indeed, the H/L ratio was marginally higher in control-call than heat-call nestlings (p = 0.066, Table 1, Figure 3C) across nest temperatures (no interaction playback by 12D-Tnest; p = 0.283, Table S1). As for HSPs, this may indicate that the range of summer nest temperatures was overall more challenging for control nestlings.
By contrast, baseline CORT levels (taken within 3 min of disturbance) did not vary with prenatal acoustic experience (p = 0.860, Table 1, Figure 3D), nor nest temperature (i.e., no effect of 12D-Tnest, p = 0.860, Table 1), including for control-call birds (no significant interaction playback by 12D-Tnest; p = 0.997, Table S1). This may be in part because thermal conditions in mid-morning at the time of sampling were never extreme.
Lastly, nestling mass was not explained by any of the four biomarkers under in-nest conditions (all p > 0.151, Table 2).

2.2. Nestling Stress Response to a Controlled Acute Heat-Challenge

After a 2.5 h heat-challenge in a temperature-controlled chamber, contrary to our prediction, HSP90α levels were overall not higher than under in-nest conditions (p = 0.344, Table 3, Figure 3B), possibly due to the relatively mild, but ecologically realistic, heating procedure we used (progressive temperature increase to a maximum of 44 °C for 20 min). HSC70 and CORT levels were also similar to those in in-nest nestlings, as predicted (p > 0.596, Table 3, Figure 3A,D). Nonetheless, the H/L ratio was significantly higher after the heat-challenge, as expected (p = 0.009, Table 3, Figure 3C).
Within the heat-challenge experiment, the prenatal playback and nest temperature that nestlings had previously experienced (12D-Tnest) did not affect any of the stress biomarkers at the end of the heat-challenge (for all, p > 0.557, Table 1, Figure 3). This suggests that developmental conditions did not explain inter-individual variations in response to an acute heat-challenge.
Instead, biomarker levels under acute heat exposure were weakly related to nestling mass at day 12 (a proxy for individual “quality” or condition). Namely, a lower H/L ratio after a heat-challenge tended to occur in heavier nestlings (p = 0.072, Table 2, Figure 4), which points to a possibly lower sensitivity to acute heat in individuals in better body condition. Nonetheless, after the heat-challenge, HSC70 levels tended to be higher in heavier nestlings (p = 0.031, Table 2, Figure 4), although this effect was potentially only driven by one individual with high HSC70 levels (linear model result without the outlier: t = 2.03, p = 0.051; Figure 4B).
Lastly, as expected based on the different time courses of the various responses [25,48], we found that the stress biomarkers were not correlated to each other, considering each experiment separately or together (Table S2).

3. Discussion

This study is the first to demonstrate that prenatal acoustic experience affects the heat shock protein (HSP) response to thermal environment. In addition to maternal ovarian HSPs or incubation temperature [34,36], this study therefore reveals a novel parental mechanism that leads to alterations in HSP levels in their offspring. It shows that a prenatal acoustic signal of heat affects physiological responses to postnatal chronic temperatures in zebra finch nestlings. While hot nest conditions were usually not extreme for this species (mean daytime temperature < 32 °C; maximum temperature reached daily ranging from 18.0 to 43.4 °C), they were more challenging for nestlings (i.e., the control-call group), which had not received prenatal anticipatory cues about upcoming hot conditions. Specifically, control-call individuals had higher levels than heat-call individuals of both the constitutive heat shock cognate 70 (HSC70) and stress-inducible HSP90α in hotter nests, and a marginally higher heterophil-to-lymphocyte (H/L) ratio overall. This study thereby demonstrates that prenatal acoustic experience contributes to biologically relevant inter-individual variation in HSP regulation. By contrast, corticosterone (CORT) levels, even though they varied between individuals, were not related to prenatal acoustic experience or current thermal conditions in either experiment. This highlights the value of considering multiple stress biomarkers in parallel, to gain better insights into the physiological impact of environmental fluctuations. Overall, our study suggests that a prenatal acoustic signal of heat may improve nestling resilience to chronically high postnatal temperatures.

3.1. Acoustic Programming of HSPs and Stress Physiology

This study found that prenatal sound was capable of altering nestling stress physiology two weeks later, in response to chronic nest temperature. This aligns with previous research on the zebra finch showing that prenatal heat-call exposure leads to inter-individual variation later in life, specifically in traits related to thermoregulation and heat responses (e.g., [44,45,47,49]), but not in other traits (e.g., foraging skills [50]). This study therefore adds to the emerging evidence for the role of prenatal sound in developmental programming, with both positive and negative fitness consequences [44,51,52]. Even though the fitness impact was not measured here, differences in HSP levels are likely to reflect an advantage for heat-call birds under chronic heat. Indeed, in hot nests protein levels for both HSPs were higher in control-call birds than in heat-call birds, which is potentially indicative of higher protein unfolding triggering the HSP upregulation. By contrast, heat-call individuals in hot nests presumably invested less energy into HSP synthesis, leaving more energy for other functions (e.g., the immune system [53]), or lowering overall energy requirements in the heat when the parental provisioning rate is reduced [54,55]. While the molecular mechanisms leading from prenatal sound exposure to differential HSP responsiveness in nestlings are unknown, it is likely that this effect occurred through an indirect pathway, via changes in individual thermal physiology. One previous study nonetheless reported that loud prenatal music increased HSP levels in chicken embryos and hatchlings [56], which also points to an acoustic sensitivity of HSP regulation and possibly other cellular processes [57].
In agreement with HSP level differences reflecting different sensitivities to heat, we found that heat-calls birds had a marginally lower H/L ratio than control-call nestlings in summer. This is consistent with results in adult zebra finches (from a different cohort), where H/L ratios were also lower in heat-call birds in summer, but also, surprisingly, in winter [46]. Both of these results, in nestlings and adults, therefore may indicate a higher resilience of heat-call individuals to thermal fluctuations, rather than solely to hot conditions. This may be related to the ecological conditions in the Australian desert, where heatwaves can trigger storms associated with a drop in temperature. Lastly, unlike in yellow-legged gulls where alarm calls signaling predator presence tends to lead to higher baseline CORT levels at hatching [58,59], we found no effect of prenatal acoustic experience on baseline CORT levels, possibly because nest temperatures were always mild at the time of sampling (in the late morning).

3.2. HSP Response to Acute and Chronic Heat

In individuals that had not received the prenatal acoustic signal of heat (i.e., control-call birds), a high chronic temperature in the nest led to increased levels of both HSP90α and HSC70. Even though HSC70 is generally classified as a constitutive HSP, there is evidence that gene expression for HSC70 (HSPA8) can be upregulated in the heat, although less readily so than HSP90-B1 [13,35,60]. Against our expectations, however, the heat-challenge in our second experiment did not trigger an increase in either HSP compared to nestlings in the first experiment. This may be because, rather than a sudden increase to an extreme temperature lasting several hours, we used a staggered temperature increase (starting at 27 °C), with the most extreme temperature stage (at 44 °C) only lasting 20 min and ending 18 min before blood sampling. This protocol, which represents a more natural heat exposure, may not have triggered significant protein unfolding (which would trigger the HSP upregulation [11]) or it may not have allowed sufficient time to detect HSP upregulation at the protein levels (which takes longer than HSP gene upregulation). The time course of the HSP response, especially for protein translation, is still unclear in birds [11,40], partly because results on cell cultures [61,62], which lack organismal regulatory pathways, may not directly extrapolate to the whole-organism kinetic. Nonetheless, our results may also be indicative of a lower sensitivity of zebra finches to acute heat than other species. For example, in three Australian desert birds, HSP gene expression was upregulated in all species, but the extent of the response across HSP genes and organs was much lower in the zebra finch than in the diamond dove (Geopelia cuneate) or budgerigar (Melopsittacus undulates) [13].

3.3. Stress Response to Acute Heat

Exposure to the acute heat-challenge nonetheless increased the H/L ratio, as also found in adults in two other desert species [24] and in poultry [42,43,63]. In addition, nestlings with a lower body mass tended to have a higher H/L ratio, as also reported in adults, in the zebra finch [46] and other avian species (e.g., [64]), although this was only seen in the heat-challenge, and not in in-nest conditions. This may indicate a higher sensitivity of smaller nestlings to acute heat but not to long-term high temperatures. Furthermore, CORT levels from heat-challenged nestlings did not differ from those sampled in in-nest conditions. Even if heat itself did not increase CORT levels, handling just before the heat-challenge would have triggered a CORT response, as observed in other studies on zebra finch nestlings exposed to a standard restraint protocol [65,66]. The similar CORT levels across the two experiments therefore do not indicate that nestlings in the heat treatment did not experience stress, but rather that their CORT levels had returned to baseline by the end of the 2.5 h heat-challenge, which is consistent with the expected time course of CORT in response to an acute stressor [48]. This short response may explain why previous studies also failed to find a CORT increase after an acute heat-challenge [24], or during heatwaves in the wild [67].

3.4. Stress Biomarker Covariation

Overall, this study suggests that the response profiles of the organismal and cellular stress biomarkers provide complementary information on an individual’s physiological state. This is also indicated by the absence of correlation between the stress biomarkers. The lack of correlation between plasma CORT levels and the H/L ratio is common [24,25,68]. By contrast, to our knowledge, relationships between HSPs and CORT or the H/L ratio have been tested only in one previous study, which also reported an absence of correlation in blackbirds (Turdus merula) [69]. That, overall, we found differences in the H/L ratio and HSPs but not in CORT, supports the usefulness of measuring other stress biomarkers in addition to CORT, in order to better interpret the effect of a stressor or chronic conditions on individuals [70,71].

3.5. Conclusions

In conclusion, this study shows that prenatal exposure to “heat-calls”, an acoustic signal of heat, affects the stress physiology response of zebra finch nestlings to postnatal chronic temperatures. It showed that heat-call-exposed individuals did not upregulate heat shock proteins in hot nests, which, combined with lower H/L ratios, is indicative of a lesser sensitivity to high temperatures than control birds. This study therefore adds to the increasing evidence in this system showing that heat-call exposure shifts offspring traits towards a heat-adapted phenotype. It also adds to the small number of studies to date on the sources of inter-individual variation in HSP regulation in wildlife, and the role of developmental programming in generating such variation.

4. Materials and Methods

4.1. Prenatal Acoustic Playback

In this study, nestlings were obtained from the same breeding experiment described in [45]. The breeding experiment (see Figure 1 for a summary diagram), was carried out at Deakin University, Geelong, Australia, from January to April 2019. In total, 59 male and 52 female wild-derived zebra finches were allowed to pair and breed freely in an outdoor aviary (12 × 6 × 3 m). Eggs were collected on laying day and replaced with dummy eggs. Freshly collected eggs were incubated in a fan-less artificial incubator (Bellsouth 100 electronic incubator) at 37.5 °C and 60% humidity. On the 10th incubation day, eggs were randomly transferred into one of two experimental incubators until hatching, 4–5 days later, and were exposed daily to either heat-calls (treatment) or control-calls (control). Heat-calls are fast, high-pitch calls given by adult zebra finches under heat stress, whereas as control-calls, “tet-calls” were used, which are contact calls commonly uttered by parents at the nest, including during incubation, but which are not related to temperature [44,72]. Both groups were also exposed to whine calls, characterized by a complex acoustic structure, to ensure normal stimulation of the auditory system. Prenatal playbacks were broadcast from 10:00 to 18:00 at 65 dB from two speakers (Sennheiser HD439, Video Guys, Melbourne, Australia) inside the incubators and externally connected to an amplifier (Digitech 18 W, Jaycar, Geelong, Australia) and an audio player (Zoom H4nSP and Marantz PMD670, Video Guys, Melbourne, Australia). Egg trays and audio players were swapped daily between the two playback incubators to prevent incubator-specific effects. Upon hatching, nestlings were returned to foster parents and broods of three-to-four nestlings were composed from both prenatal playback groups.

4.2. Postnatal Nest Temperature

The nest temperature varied naturally with air temperature, and was, in addition, experimentally manipulated during daytime (10:00–18:00) throughout nestling development (from hatching to 13 days post-hatch) to obtain a gradient of cool-to-hot nests. Cool nests were provided shade using a shading cloth on the roof of the aviary and cool pads (for air temperature (Ta) between 25 °C and 30 °C) or frozen pads (for Ta > 30 °C) placed under the nest-box roof. Hotter nests were provided heat pads: one placed under the nest-box roof (for Ta < 30 °C, Medi Heat Pack®, Melbourne, Australia), and a second one at the back of the nest, inside the nest-box (for Ta < 25 °C, Hotteeze®, Melbourne, Australia). All pads were placed at 10:00, replaced at 13:30 and 15:30, and removed at 18:00. The nest temperatures were continuously recorded with thermometers (Minnow-1.0TH, Senonics, Instrument Choice, Adelaide, Australia) placed in the nest boxes. The average nest temperature experienced by each nestling from hatching to 12 days post-hatch during daytime (11:00–18:00) was calculated, referred to as “12D-Tnest”. 12D-Tnest per nestling varied from 23.1 to 33.8 °C. Within the daytime temperature variation, the highest nest temperature reached per day varied from 18.0 to 43.4 °C (mean = 30.0 °C).

4.3. In-Nest and Heat-Challenge Experiments

All individuals were weighed at 12 days post-hatch between 15:00 and 15:30 to determine their body mass. At day 13, nestlings were allocated to either the “in-nest” or “heat-challenge” experiment. In the in-nest experiment, nestlings were picked up from their nests, between 11.30 and 12.45, immediately prior blood sampling. In the heat-challenge experiment, nestlings were exposed to an acute heat-challenge. Specifically, they were picked up from their nests at 09:30 and food-deprived (for metabolic measurement) until blood sampling. They were acclimated at room temperature (25 °C) for 30 min, then at 29 °C in a temperature-controlled chamber. Then, they were exposed to 20–30 min stages at 35, 40, 42 and 44 °C in the chamber, before the temperature was dropped to 35 °C for 15 min. The heat-challenge lasted 2.5 h in total, except for two individuals showing distress signs that required the experiment to stop at 42 °C.
Within 3 min of collection from the nest or the end of the heat-challenge, nestlings were euthanized using isoflurane and decapitation to collect tissues as part of a separate study. Up to 250 µL of blood from the jugular vein was collected with heparinized capillary tubes, a blood smear was prepared (see below) and the remaining blood was expelled in microcentrifuge tubes. Blood samples were kept on ice for up to 30 min until centrifugation, at 2000× g for 10 min at 4 °C. Plasma and red blood cell (RBC) samples were stored at −80 °C until assayed for corticosterone (CORT) and heat shock proteins (HSPs). Up to two samples were collected per day, and, in total, 75 blood samples and smears were collected.
All procedures described in this study were approved by the Animal Ethics Committee of Deakin University (G23-2018) and were carried out in accordance with the Australian guidelines and regulations for the use of animals in research.

4.4. Corticosterone Assays

Corticosterone levels were quantified in 53 plasma samples using Enzyme Immunoassay (EIA) kits (ADI-901-097, Enzo Life Sciences, Farmingdale NY, USA). All samples were spiked with tritiated CORT (NET399250UC; Perkin Elmer, Sydney, Australia) to determine the recovery percentage. CORT was extracted twice from plasma using dichloromethane. After evaporation under nitrogen gas, samples were reconstituted in buffer solution (1:25 dilution). We ran the samples following an adjusted protocol using half the volume of the reagents provided in the EIA kit [72]. Every plate had a six-point standard curve ranging from 20,000 pg/mL to 6.4 pg/mL and an external CORT standard of 500 pg/mL to calculate the inter-plate variation. All samples and standards were run in triplicates. The plates were read using a Varioskan Lux microplate reader (Thermo Fisher Scientific, Melbourne, Australia) at 405 nm corrected at 580 nm. The CORT concentrations were corrected by the sample recoveries (the mean recovery was 69%). The intra- and inter-plate coefficients of variation were 4.15% and 9%, respectively.

4.5. HSPs Assays

Levels of HSC70 and HSP90α were quantified in 73 RBC samples using EIA kits (SKT-106, SKT-107, StressMarq Biosciences Inc., Victoria, BC, Canada) following the manufacturer’s instructions and a published protocol [40] that was optimized in the laboratory to our samples. For protein extraction, 2 µL of RBCs were diluted with 1× extraction reagent solution (1:500 dilution) containing protease inhibitors (Cat. No. 04693159001, Sigma Aldrich, Melbourne, Australia). The samples were incubated on ice for 30 min and vortexed for 10 s at T0, T10, T20 and T30 min. Then, the samples were centrifuged at 4500 RPM (4520× g) for 10 min at 4 °C. The supernatant was collected and stored at −80 °C until further processing.
Bradford assays [73] were then performed to quantify the total protein (TP) content of the samples. Briefly, 200 µL of Bradford reagent (B6916, Sigma-Aldrich) was added to 5 µL of the sample into the wells of a plate (NUN475094, Thermo Fisher Scientific). The standard curves were generated using a known concentration of bovine serum albumin (P0834, Sigma-Aldrich). Plates were incubated for 5 min at room temperature and read at 595 nm with a Varioskan Lux microplate reader.
To assay HSC70 and HSP90α, samples were run undiluted and diluted (1:2) with the sample diluent supplied in the kit, respectively. Every plate had a seven-point standard curve ranging from 150 ng/mL to 2.34 ng/mL for HSC70, and 28 ng/mL to 0.44 ng/mL for HSP90α, as well as an inter-assay standard to calculate inter-plate variation. All samples and standards were run in duplicate. The plates were read at 450 nm using a Varioskan Lux microplate reader. The HSP concentrations were corrected by TP content. The intra- and inter-plates variations were 7.7% and 12.9% for HSC70, and 8.7% and 25% for HSP90α.

4.6. Heterophil-to-Lymphocyte Ratio

Blood smears were prepared by smearing one drop of fresh blood on a glass slide to obtain a thin blood film. All smears were immediately air-dried and stained using Quick Dip staining (Fronine, Sydney, Australia). Specifically, the smears were fixed with the fixative solution (methanol) and dipped seven times successively into the two stain solutions (eosin Y and methylene azure B). Slides were rinsed (water), air-dried and stored in the dark at room temperature until analysis by a single observer (EU).
Leucocyte counts were performed using a compound microscope (Olympus Australia, Melbourne, Australia) at 1000× magnification using oil immersion. The first 100 leucocytes were counted, differentiating heterophils, eosinophils, basophils, lymphocytes and monocytes according to standard avian guidelines [74]. Then, the H/L ratio was calculated by dividing the number of heterophils by the number of lymphocytes. In total, leucocyte counts from 56 individuals were obtained.
To assess the repeatability of the measurements, a subset of 18 slides was examined a second time. Eight additional slides that showed an elevated H/L ratio (>3) were also counted a second time to verify if the first count was correct.

4.7. Statistical Analyses

All statistical analyses were conducted using R version 4.4.1. Physiological traits were not all measured in all individuals, and 13 blood smears were too small and/or thin to count 100 leucocytes, resulting in varying sample sizes: in the in-nest experiment, N = 25, 35, 35 and 28 birds for CORT, HSC70, HSP90α and the H/L ratio, respectively; in the heat-challenge experiment, sample sizes were 28, 38, 38 and 28 birds, respectively.
All models were fitted using the packages lme4 (v 1.1-35-5) and lmerTest (v 3.1-3). Continuous predictors were scaled (i.e., mean = 0, std = 1), and the normality and homoscedasticity of residuals were visually inspected. Furthermore, collinearity among predictors was tested in the full models by calculating the variance inflation factors (VIFs), which were always lower than three, indicating no significant collinearity [75].

4.7.1. H/L Ratio Distribution

While CORT and HSP levels were centered around the mean, the H/L ratio distribution was bimodal, with a great majority of nestlings with a ratio < 2 (i.e., 82% in the in-nest experiment [range: 0.18–1.15] and 75% in the heat-challenge experiment [0.10–1.78]) and a few individuals with high values > 2 (in-nest: 5 nestlings [2.22–6.91], heat-challenge: 7 nestlings [2.17–13.17]; Figure S1). These high H/L ratios were not the result of miscounts, as the H/L ratio, including for nestlings with high H/L ratios, was repeatable among the blood smears counted twice [intra-class correlations ICCs: N = 26, R = 0.73, CI = (0.49, 0.87), p = 0.001]. It is not known whether methodological issues other than cell identification and counting could lead to extreme H/L ratios. Although similarly extreme values in a few individuals have been reported in nestlings of another avian species [76], values > 2 were higher than any observed in zebra finch adults [46]. Given that the number of circulating heterophils and lymphocytes vary with other factors than stress [23,77,78], this suggests that nestlings with H/L ratios > 2 were experiencing an unknown challenge or condition (such as, for example, an ongoing bacterial infection), unrelated to the experimental manipulations. Accordingly, even after the heat-challenge, the H/L ratio distribution was still bimodal, which confirms that extreme H/L ratio values were unrelated to the experiment, and to the heat stress response. Given the uncertainty around the interpretation of these values, and since extreme values may mask or bias experimental effects and make data distribution less suitable for mixed models, H/L ratios > 2 were excluded from the analyses using the H/L ratio as a response variable. These individuals had, however, average levels of CORT and of both HSPs (Figure S2), which again points to a different driver than stress leading to extreme H/L ratios.

4.7.2. Developmental Effects Within Each Experiment

First, for the experiment under in-nest conditions (chronic exposure to manipulated nest temperature), the effects of the developmental conditions (i.e., prenatal playback and postnatal nest temperature) were assessed on each stress biomarker after log-transformation (CORT, HSC70, HSP90α and H/L ratio). Linear mixed models (LMMs) were used, including the prenatal playback and nest temperature (12D-Tnest) as predictors, and their interaction. Brood identity (rearing nest; “brood-ID”) was also included as a random factor, unless its variance was null (see Table S1). Brood size was not added as a covariate because brood size was standardized to four nestlings, and was consistently non-significant in preliminary analyses. Considering clutch identity (to account for genetic similarities) instead of brood identity as a random effect led to null or extremely low variance.
Second, for the acute heat-challenge experiment, developmental effects in heat-challenged individuals were tested, using the same models as for the in-nest experiment above. Expressing the nest temperature as the chamber temperature deviation from the nest temperature (i.e., chamber temperature above 12D-Tnest: Δ12D-Tnest, as in [45]) instead of the nest temperature did not change any of the results.
For all models, the interaction (playback × 12D-Tnest) was removed when it was non-significant. Models with the interaction are reported in the Supplementary Materials (Table S1).

4.7.3. Variation Across Experiments (In-Nest vs. Heat-Challenge)

To test whether stress biomarkers were affected by the heat-challenge, LMMs were run using each stress biomarker as the response variable, and the experiment (in-nest or heat-challenge) as a predictor, with brood identity (rearing nest; “brood-ID”) as a random factor, unless its variance was null (see Table 3).

4.7.4. Link Between Nestling Mass and Stress Physiology

In each experiment, to test how the stress biomarkers may impact nestling mass, two separate LMMs were run, including mass as the response variable, brood-ID as a random factor and either CORT and the H/L ratio or both HSPs as predictors (log-transformed). Not all four stress biomarkers were included in a single analysis as this would considerably reduce the sample size because of missing values, although it led to the same results. As reported elsewhere [24,25,68,69], because of different response time courses, that the stress biomarkers were not correlated to each other was verified, considering each experiment separately or together (Table S2).

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms252212194/s1.

Author Contributions

Conceptualization, M.M.M. and E.U.; methodology, O.L.C., B.A. and K.L.B.; validation, O.L.C., B.A. and E.U.; formal analysis, E.U., M.M.M. and O.L.C.; investigation, E.U., A.P. and M.M.M.; writing—original draft preparation, M.M.M. and E.U.; writing—review and editing, all authors; visualization, E.U.; supervision, M.M.M., O.L.C., B.A. and K.L.B.; project administration, M.M.M. and K.L.B.; funding acquisition, M.M.M. and K.L.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Australian Research Council (grant numbers DE170100824 (M.M.M.), DP180101207 (K.L.B., M.M.M.) and DP210101238 (K.L.B., M.M.M.)) and the Spanish Ministry of Science, Innovation and Universities (grant RYC2019-028066-I (M.M.M.)). The APC was covered by the Max Plank Institute for Biological Intelligence.

Institutional Review Board Statement

The animal study protocol was approved by the Animal Ethics Committee of Deakin University (G23-2018, approved on 22 January 2019).

Data Availability Statement

The data supporting reported results are publicly available on the Mendeley repository: https://doi.org/10.17632/mrc9yzngwm.1.

Acknowledgments

We thank Nicolas de Almeida for his help during the breeding experiment, as well as Deakin Animal House staff for logistical support. We also thank Haruka Wada for sharing with us her HSP assay protocol, and Marja A. Oscos-Snowball for providing training in leucocyte identification.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Romero, L.M.; Wingfield, J.C. Tempests, Poxes, Predators, and People: Stress in Wild Animals and How They Cope; Oxford Series in Behavioral Neuroendocrinology; Oxford University Press: New York, NY, USA, 2015; ISBN 978-0-19-536669-3. [Google Scholar]
  2. Sørensen, J.G. Application of Heat Shock Protein Expression for Detecting Natural Adaptation and Exposure to Stress in Natural Populations. Curr. Zool. 2010, 56, 703–713. [Google Scholar] [CrossRef]
  3. Angelier, F.; Wingfield, J.C. Importance of the Glucocorticoid Stress Response in a Changing World: Theory, Hypotheses and Perspectives. Gen. Comp. Endocrinol. 2013, 190, 118–128. [Google Scholar] [CrossRef] [PubMed]
  4. Chen, B.; Feder, M.E.; Kang, L. Evolution of Heat-Shock Protein Expression Underlying Adaptive Responses to Environmental Stress. Mol. Ecol. 2018, 27, 3040–3054. [Google Scholar] [CrossRef] [PubMed]
  5. Guindre-Parker, S. Individual Variation in Glucocorticoid Plasticity: Considerations and Future Directions. Integr. Comp. Biol. 2020, 60, 79–88. [Google Scholar] [CrossRef]
  6. Jolly, C.; Morimoto, R.I. Role of the Heat Shock Response and Molecular Chaperones in Oncogenesis and Cell Death. JNCI J. Natl. Cancer Inst. 2000, 92, 1564–1572. [Google Scholar] [CrossRef]
  7. Richter, K.; Haslbeck, M.; Buchner, J. The Heat Shock Response: Life on the Verge of Death. Mol. Cell 2010, 40, 253–266. [Google Scholar] [CrossRef]
  8. Kültz, D. Molecular and evolutionary basis of the cellular stress response. Annu. Rev. Physiol. 2005, 67, 225–257. [Google Scholar] [CrossRef]
  9. Feder, M.E.; Hofmann, G.E. Heat-Shock Proteins, Molecular Chaperones, and the Stress Response: Evolutionary and Ecological Physiology. Annu. Rev. Physiol. 1999, 61, 243–282. [Google Scholar] [CrossRef]
  10. Kregel, K.C. Invited Review: Heat Shock Proteins: Modifying Factors in Physiological Stress Responses and Acquired Thermotolerance. J. Appl. Physiol. 2002, 92, 2177–2186. [Google Scholar] [CrossRef]
  11. Balakrishnan, K.N.; Ramiah, S.K.; Zulkifli, I. Heat Shock Protein Response to Stress in Poultry: A Review. Animals 2023, 13, 317. [Google Scholar] [CrossRef]
  12. Xie, J.; Tang, L.; Lu, L.; Zhang, L.; Xi, L.; Liu, H.C.; Odle, J.; Luo, X. Differential Expression of Heat Shock Transcription Factors and Heat Shock Proteins after Acute and Chronic Heat Stress in Laying Chickens (Gallus gallus). PLoS ONE 2014, 9, e102204. [Google Scholar] [CrossRef] [PubMed]
  13. Xie, S.; Tearle, R.; McWhorter, T.J. Heat Shock Protein Expression Is Upregulated after Acute Heat Exposure in Three Species of Australian Desert Birds. Avian Biol. Res. 2018, 11, 263–273. [Google Scholar] [CrossRef]
  14. Fangue, N.A.; Hofmeister, M.; Schulte, P.M. Intraspecific Variation in Thermal Tolerance and Heat Shock Protein Gene Expression in Common Killifish, Fundulus Heteroclitus. J. Exp. Biol. 2006, 209, 2859–2872. [Google Scholar] [CrossRef] [PubMed]
  15. Roberts, D.A.; Hofmann, G.E.; Somero, G.N. Heat-Shock Protein Expression in Mytilus californianus: Acclimatization (Seasonal and Tidal-Height Comparisons) and Acclimation Effects. Biol. Bull. 1997, 192, 309–320. [Google Scholar] [CrossRef] [PubMed]
  16. Sørensen, J.G.; Kristensen, T.N.; Loeschcke, V. The Evolutionary and Ecological Role of Heat Shock Proteins. Ecol. Lett. 2003, 6, 1025–1037. [Google Scholar] [CrossRef]
  17. de Bruijn, R.; Romero, L.M. The Role of Glucocorticoids in the Vertebrate Response to Weather. Gen. Comp. Endocrinol. 2018, 269, 11–32. [Google Scholar] [CrossRef]
  18. Romero, L.M. Physiological Stress in Ecology: Lessons from Biomedical Research. Trends Ecol. Evol. 2004, 19, 249–255. [Google Scholar] [CrossRef]
  19. Baker, J.D.; Ozsan, I.; Rodriguez Ospina, S.; Gulick, D.; Blair, L.J. Hsp90 Heterocomplexes Regulate Steroid Hormone Receptors: From Stress Response to Psychiatric Disease. Int. J. Mol. Sci. 2018, 20, 79. [Google Scholar] [CrossRef]
  20. Grad, I.; Picard, D. The Glucocorticoid Responses Are Shaped by Molecular Chaperones. Mol. Cell. Endocrinol. 2007, 275, 2–12. [Google Scholar] [CrossRef]
  21. Dhabhar, F.S.; Miller, A.H.; McEwen, B.S.; Spencer, R.L. Effects of Stress on Immune Cell Distribution. Dynamics and Hormonal Mechanisms. J. Immunol. (Baltim. Md. 1950) 1995, 154, 5511–5527. [Google Scholar] [CrossRef]
  22. Minias, P. Evolution of Heterophil/Lymphocyte Ratios in Response to Ecological and Life-History Traits: A Comparative Analysis across the Avian Tree of Life. J. Anim. Ecol. 2019, 88, 554–565. [Google Scholar] [CrossRef] [PubMed]
  23. Davis, A.K.; Maney, D.L.; Maerz, J.C. The Use of Leukocyte Profiles to Measure Stress in Vertebrates: A Review for Ecologists. Funct. Ecol. 2008, 22, 760–772. [Google Scholar] [CrossRef]
  24. Xie, S.; Romero, L.M.; Htut, Z.W.; McWhorter, T.J. Stress Responses to Heat Exposure in Three Species of Australian Desert Birds. Physiol. Biochem. Zool. 2017, 90, 348–358. [Google Scholar] [CrossRef]
  25. Davis, A.K.; Maney, D.L. The Use of Glucocorticoid Hormones or Leucocyte Profiles to Measure Stress in Vertebrates: What’s the Difference? Methods Ecol. Evol. 2018, 9, 1556–1568. [Google Scholar] [CrossRef]
  26. Gross, W.B.; Siegel, H.S. Evaluation of the Heterophil/Lymphocyte Ratio as a Measure of Stress in Chickens. Avian Dis. 1983, 27, 972–979. [Google Scholar] [CrossRef] [PubMed]
  27. Woodruff, M.J.; Zimmer, C.; Ardia, D.R.; Vitousek, M.N.; Rosvall, K.A. Heat Shock Protein Gene Expression Varies among Tissues and Populations in Free-Living Birds. Ornithology 2022, 139, ukac018. [Google Scholar] [CrossRef]
  28. Tedeschi, J.N.; Kennington, W.J.; Tomkins, J.L.; Berry, O.; Whiting, S.; Meekan, M.G.; Mitchell, N.J. HeriTable Variation in Heat Shock Gene Expression: A Potential Mechanism for Adaptation to Thermal Stress in Embryos of Sea Turtles. Proc. R. Soc. B Biol. Sci. 2016, 283, 20152320. [Google Scholar] [CrossRef]
  29. McColl, G.; Hoffmann, A.A.; McKechnie, S.W. Response of Two Heat Shock Genes to Selection for Knockdown Heat Resistance in Drosophila Melanogaster. Genetics 1996, 143, 1615–1627. [Google Scholar] [CrossRef]
  30. Dixon, A.L.; Liang, L.; Moffatt, M.F.; Chen, W.; Heath, S.; Wong, K.C.C.; Taylor, J.; Burnett, E.; Gut, I.; Farrall, M.; et al. A Genome-Wide Association Study of Global Gene Expression. Nat. Genet. 2007, 39, 1202–1207. [Google Scholar] [CrossRef]
  31. Kapoor, A.; Dunn, E.; Kostaki, A.; Andrews, M.H.; Matthews, S.G. Fetal Programming of Hypothalamo-Pituitary-Adrenal Function: Prenatal Stress and Glucocorticoids. J. Physiol. 2006, 572, 31–44. [Google Scholar] [CrossRef]
  32. Groothuis, T.G.G.; Müller, W.; Von Engelhardt, N.; Carere, C.; Eising, C. Maternal Hormones as a Tool to Adjust Offspring Phenotype in Avian Species. Neurosci. Biobehav. Rev. 2005, 29, 329–352. [Google Scholar] [CrossRef] [PubMed]
  33. Lockwood, B.L.; Julick, C.R.; Montooth, K.L. Maternal Loading of a Small Heat Shock Protein Increases Embryo Thermal Tolerance in Drosophila Melanogaster. J. Exp. Biol. 2017, 220, 4492–4501. [Google Scholar] [CrossRef] [PubMed]
  34. Lim, M.Y.-T.; Bernier, N.J. Zebrafish Parental Progeny Investment in Response to Cycling Thermal Stress and Hypoxia: Deposition of Heat Shock Proteins but Not Cortisol. J. Exp. Biol. 2022, 225, jeb244715. [Google Scholar] [CrossRef] [PubMed]
  35. Galletta, L.; Craven, M.J.; Meillère, A.; Crowley, T.M.; Buchanan, K.L.; Mariette, M.M. Acute Exposure to High Temperature Affects Expression of Heat Shock Proteins in Altricial Avian Embryos. J. Therm. Biol. 2022, 110, 103347. [Google Scholar] [CrossRef]
  36. Vinoth, A.; Thirunalasundari, T.; Tharian, J.A.; Shanmugam, M.; Rajkumar, U. Effect of Thermal Manipulation during Embryogenesis on Liver Heat Shock Protein Expression in Chronic Heat Stressed Colored Broiler Chickens. J. Therm. Biol. 2015, 53, 162–171. [Google Scholar] [CrossRef]
  37. Hoffman, A.J.; Finger, J.W.; Kavazis, A.N.; Wada, H. Early Life Thermal Conditioning Alters Heat-Shock Protein Expression in Response to an Adult Thermal Stressor. J. Exp. Zool. A Ecol. Integr. Physiol. 2024, 341, 1030–1040. [Google Scholar] [CrossRef]
  38. Riddell, E.A.; Iknayan, K.J.; Hargrove, L.; Tremor, S.; Patton, J.L.; Ramirez, R.; Wolf, B.O.; Beissinger, S.R. Exposure to Climate Change Drives Stability or Collapse of Desert Mammal and Bird Communities. Science 2021, 371, 633–636. [Google Scholar] [CrossRef]
  39. Conradie, S.R.; Woodborne, S.M.; Wolf, B.O.; Pessato, A.; Mariette, M.M.; McKechnie, A.E. Avian Mortality Risk during Heat Waves Will Increase Greatly in Arid Australia during the 21st Century. Conserv. Physiol. 2020, 8, coaa048. [Google Scholar] [CrossRef]
  40. Finger, J.W.; Hoffman, A.J.; Wada, H. Temporal Variation in Constitutive and Inducible Heat Shock Proteins in the Zebra Finch Taeniopygia Guttata. Ibis 2018, 160, 390–396. [Google Scholar] [CrossRef]
  41. Soleimani, A.F.; Zulkifli, I.; Omar, A.R.; Raha, A.R. Physiological Responses of 3 Chicken Breeds to Acute Heat Stress. Poult. Sci. 2011, 90, 1435–1440. [Google Scholar] [CrossRef]
  42. Altan, Ö.; Pabuçcuoǧlu, A.; Altan, A.; Konyalioǧlu, S.; Bayraktar, H. Effect of Heat Stress on Oxidative Stress, Lipid Peroxidation and Some Stress Parameters in Broilers. Br. Poult. Sci. 2003, 44, 545–550. [Google Scholar] [CrossRef] [PubMed]
  43. Prieto, M.T.; Campo, J.L. Effect of Heat and Several Additives Related to Stress Levels on Fluctuating Asymmetry, Heterophil: Lymphocyte Ratio, and Tonic Immobility Duration in White Leghorn Chicks. Poult. Sci. 2010, 89, 2071–2077. [Google Scholar] [CrossRef] [PubMed]
  44. Mariette, M.M.; Buchanan, K.L. Prenatal Acoustic Communication Programs Offspring for High Posthatching Temperatures in a Songbird. Science 2016, 353, 812–814. [Google Scholar] [CrossRef]
  45. Udino, E.; George, J.M.; McKenzie, M.; Pessato, A.; Crino, O.L.; Buchanan, K.L.; Mariette, M.M. Prenatal Acoustic Programming of Mitochondrial Function for High Temperatures in an Arid-Adapted Bird. Proc. R. Soc. B 2021, 288, 20211893. [Google Scholar] [CrossRef] [PubMed]
  46. Udino, E.; Oscos-Snowball, M.A.; Buchanan, K.L.; Mariette, M.M. A Prenatal Acoustic Signal of Heat Reduces a Biomarker of Chronic Stress at Adulthood across Seasons. Front. Physiol. 2024, 15, 1348993. [Google Scholar] [CrossRef] [PubMed]
  47. Pessato, A.; McKechnie, A.E.; Mariette, M.M. A Prenatal Acoustic Signal of Heat Affects Thermoregulation Capacities at Adulthood in an Arid-Adapted Bird. Sci. Rep. 2022, 12, 5842. [Google Scholar] [CrossRef]
  48. Gormally, B.M.G.; Romero, L.M. What Are You Actually Measuring? A Review of Techniques That Integrate the Stress Response on Distinct Time-scales. Funct. Ecol. 2020, 34, 2030–2044. [Google Scholar] [CrossRef]
  49. Udino, E.; Mariette, M.M. How to Stay Cool: Early Acoustic and Thermal Experience Alters Individual Behavioural Thermoregulation in the Heat. Front. Ecol. Evol. 2022, 10, 818278. [Google Scholar] [CrossRef]
  50. Katsis, A.C.; Buchanan, K.L.; Kleindorfer, S.; Mariette, M.M. Long-Term Effects of Prenatal Sound Experience on Songbird Behavior and Their Relation to Song Learning. Behav. Ecol. Sociobiol. 2021, 75, 18. [Google Scholar] [CrossRef]
  51. Mariette, M.M.; Clayton, D.F.; Buchanan, K.L. Acoustic Developmental Programming: A Mechanistic and Evolutionary Framework. Trends Ecol. Evol. 2021, 36, 722–736. [Google Scholar] [CrossRef]
  52. Meillère, A.; Buchanan, K.L.; Eastwood, J.R.; Mariette, M.M. Pre- and Postnatal Noise Directly Impairs Avian Development, with Fitness Consequences. Science 2024, 384, 475–480. [Google Scholar] [CrossRef] [PubMed]
  53. Morales, J.; Moreno, J.; Lobato, E.; Merino, S.; Tomás, G.; Puente, J.M.D.L.; Martínez, J. Higher Stress Protein Levels Are Associated with Lower Humoral and Cell-Mediated Immune Responses in Pied Flycatcher Females. Funct. Ecol. 2006, 20, 647–655. [Google Scholar] [CrossRef]
  54. Barras, A.G.; Niffenegger, C.A.; Candolfi, I.; Hunziker, Y.A.; Arlettaz, R. Nestling Diet and Parental Food Provisioning in a Declining Mountain Passerine Reveal High Sensitivity to Climate Change. J. Avian Biol. 2021, 52, jav.02649. [Google Scholar] [CrossRef]
  55. Oswald, K.N.; Smit, B.; Lee, A.T.K.; Peng, C.L.; Brock, C.; Cunningham, S.J. Higher Temperatures Are Associated with Reduced Nestling Body Condition in a Range-Restricted Mountain Bird. J. Avian Biol. 2021, 52, e02756. [Google Scholar] [CrossRef]
  56. Hanafi, S.A.; Zulkifli, I.; Ramiah, S.K.; Chung, E.L.T.; Kamil, R.; Awad, E.A. Prenatal Auditory Stimulation Induces Physiological Stress Responses in Developing Embryos and Newly Hatched Chicks. Poult. Sci. 2023, 102, 102390. [Google Scholar] [CrossRef]
  57. Mariette, M.M. Developmental Programming by Prenatal Sounds: Insights into Possible Mechanisms. J. Exp. Biol. 2024, 227, jeb246696. [Google Scholar] [CrossRef]
  58. Noguera, J.C.; Velando, A. Bird Embryos Perceive Vibratory Cues of Predation Risk from Clutch Mates. Nat. Ecol. Evol. 2019, 3, 1225–1232. [Google Scholar] [CrossRef]
  59. Noguera, J.C.; Velando, A. Reduced Telomere Length in Embryos Exposed to Predator Cues. J. Exp. Biol. 2019, 222, jeb.216176. [Google Scholar] [CrossRef]
  60. Fangue, N.A.; Richards, J.G.; Schulte, P.M. Do Mitochondrial Properties Explain Intraspecific Variation in Thermal Tolerance? J. Exp. Biol. 2009, 212, 514–522. [Google Scholar] [CrossRef] [PubMed]
  61. Lindquist, S. The Heat-Shock Response. Annu. Rev. Biochem. 1986, 55, 1151–1191. [Google Scholar] [CrossRef]
  62. Pawar, S.S.; Kurade, N.P.; Bhendarkar, M.P.; Bhosale, S.V.; Nirmale, A.V.; Kochewad, S.A. Modulation of Heat Shock Protein 70 (HSP70) Gene Expression Ex Vivo in Response to Heat Stress in Chicken. Anim. Biotechnol. 2023, 34, 5168–5172. [Google Scholar] [CrossRef] [PubMed]
  63. Olfati, A.; Mojtahedin, A.; Sadeghi, T.; Akbari, M.; Martínez-Pastor, F. Comparison of Growth Performance and Immune Responses of Broiler Chicks Reared under Heat Stress, Cold Stress and Thermoneutral Conditions. Span. J. Agric. Res. 2018, 16, e0505. [Google Scholar] [CrossRef]
  64. Gladbach, A.; Gladbach, D.J.; Quillfeldt, P. Variations in Leucocyte Profiles and Plasma Biochemistry Are Related to Different Aspects of Parental Investment in Male and Female Upland Geese Chloephaga Picta Leucoptera. Comp. Biochem. Physiol. Part A Mol. Integr. Physiol. 2010, 156, 269–277. [Google Scholar] [CrossRef] [PubMed]
  65. Wada, H.; Kriengwatana, B.; Allen, N.; Schmidt, K.L.; Soma, K.K.; MacDougall-Shackleton, S.A. Transient and Permanent Effects of Suboptimal Incubation Temperatures on Growth, Metabolic Rate, Immune Function and Adrenocortical Responses in Zebra Finches. J. Exp. Biol. 2015, 218, 2847–2855. [Google Scholar] [CrossRef]
  66. Crino, O.L.; Driscoll, S.C.; Brandl, H.B.; Buchanan, K.L.; Griffith, S.C. Under the Weather: Corticosterone Levels in Wild Nestlings Are Associated with Ambient Temperature and Wind. Gen. Comp. Endocrinol. 2020, 285, 113247. [Google Scholar] [CrossRef]
  67. Cooper, C.E.; Hurley, L.L.; Deviche, P.; Griffith, S.C. Physiological Responses of Wild Zebra Finches (Taeniopygia guttata) to Heatwaves. J. Exp. Biol. 2020, 223, jeb225524. [Google Scholar] [CrossRef]
  68. Müller, C.; Jenni-Eiermann, S.; Jenni, L. Heterophils/Lymphocytes-Ratio and Circulating Corticosterone Do Not Indicate the Same Stress Imposed on Eurasian Kestrel Nestlings. Funct. Ecol. 2011, 25, 566–576. [Google Scholar] [CrossRef]
  69. Ibáñez-Álamo, J.D.; Jimeno, B.; Gil, D.; Thomson, R.L.; Aguirre, J.I.; Díez-Fernández, A.; Faivre, B.; Tieleman, B.I.; Figuerola, J. Physiological Stress Does Not Increase with Urbanization in European Blackbirds: Evidence from Hormonal, Immunological and Cellular Indicators. Sci. Total Environ. 2020, 721, 137332. [Google Scholar] [CrossRef]
  70. Breuner, C.W.; Delehanty, B.; Boonstra, R. Evaluating Stress in Natural Populations of Vertebrates: Total CORT Is Not Good Enough. Funct. Ecol. 2013, 27, 24–36. [Google Scholar] [CrossRef]
  71. Telemeco, R.S.; Gangloff, E.J. Analyzing Stress as a Multivariate Phenotype. Integr. Comp. Biol. 2020, 60, 70–78. [Google Scholar] [CrossRef]
  72. Crino, O.L.; Buchanan, K.L.; Trompf, L.; Mainwaring, M.C.; Griffith, S.C. Stress Reactivity, Condition, and Foraging Behavior in Zebra Finches: Effects on Boldness, Exploration, and Sociality. Gen. Comp. Endocrinol. 2017, 244, 101–107. [Google Scholar] [CrossRef] [PubMed]
  73. Bradford, M.M. A Rapid and Sensitive Method for the Quantitation of Microgram Quantities of Protein Utilizing the Principle of Protein-Dye Binding. Anal. Biochem. 1976, 72, 248–254. [Google Scholar] [CrossRef] [PubMed]
  74. Campbell, T.W. (Ed.) Exotic Animal Hematology and Cytology, 4th ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2015; ISBN 978-1-118-99370-5. [Google Scholar]
  75. Zuur, A.F.; Ieno, E.N.; Elphick, C.S. A Protocol for Data Exploration to Avoid Common Statistical Problems. Methods Ecol. Evol. 2010, 1, 3–14. [Google Scholar] [CrossRef]
  76. Masello, J.F.; Choconi, R.G.; Helmer, M.; Kremberg, T.; Lubjuhn, T.; Quillfeldt, P. Do Leucocytes Reflect Condition in Nestling Burrowing Parrots Cyanoliseus Patagonus in the Wild? Comp. Biochem. Physiol.—A Mol. Integr. Physiol. 2009, 152, 176–181. [Google Scholar] [CrossRef] [PubMed]
  77. Colominas-Ciuró, R.; Masero, J.A.; Benzal, J.; Bertellotti, M.; Barbosa, A. Oxidative Status and Stress during Highly Energetic Life-History Stages of Chinstrap Penguins: Breeding versus Molting. J. Field Ornithol. 2019, 90, 190–199. [Google Scholar] [CrossRef]
  78. Skwarska, J. Variation of Heterophil-to-Lymphocyte Ratio in the Great Tit Parus Major—A Review. Acta Ornithol. 2019, 53, 103. [Google Scholar] [CrossRef]
Figure 1. Experimental design: eggs were collected on the day of lay (embryonic day 0, E0) and artificially incubated until hatching (day 0, D0). Prenatal playback exposure (in blue) to either heat-calls or control-calls occurred in the last 5 days before hatching (E10 to E14), and the postnatal nest temperature (in green) was manipulated daily from D0 to D12 post-hatch. Nestlings’ body masses were measured on D12, and on D13 nestlings were separated into two experiments (in yellow), and a blood sample collected, to test the impact of prenatal acoustic experience on either (i) Experiment 1 (Exp1) in-nest: chronic response to nest temperature, or (ii) Experiment 2 (Exp2) heat-challenge: acute response to a heat-challenge in a temperature-controlled chamber.
Figure 1. Experimental design: eggs were collected on the day of lay (embryonic day 0, E0) and artificially incubated until hatching (day 0, D0). Prenatal playback exposure (in blue) to either heat-calls or control-calls occurred in the last 5 days before hatching (E10 to E14), and the postnatal nest temperature (in green) was manipulated daily from D0 to D12 post-hatch. Nestlings’ body masses were measured on D12, and on D13 nestlings were separated into two experiments (in yellow), and a blood sample collected, to test the impact of prenatal acoustic experience on either (i) Experiment 1 (Exp1) in-nest: chronic response to nest temperature, or (ii) Experiment 2 (Exp2) heat-challenge: acute response to a heat-challenge in a temperature-controlled chamber.
Ijms 25 12194 g001
Figure 2. Levels of heat shock proteins across nest temperature in in-nest nestlings. Individuals were exposed to prenatal control-calls (black circles) or heat-calls (red triangles) and, upon hatching until day 12, experienced a mean daytime nest temperature (12D-Tnest). At day 13, the nestlings were taken from their nest to collect blood samples. (A) Levels of constitutive heat shock cognate 70 (HSC70; N = 35) and (B) stress-inducible heat shock protein 90α (HSP90α; N = 35) were measured in the red blood cells. Open circles show individual data and regression lines are represented with 95% CIs.
Figure 2. Levels of heat shock proteins across nest temperature in in-nest nestlings. Individuals were exposed to prenatal control-calls (black circles) or heat-calls (red triangles) and, upon hatching until day 12, experienced a mean daytime nest temperature (12D-Tnest). At day 13, the nestlings were taken from their nest to collect blood samples. (A) Levels of constitutive heat shock cognate 70 (HSC70; N = 35) and (B) stress-inducible heat shock protein 90α (HSP90α; N = 35) were measured in the red blood cells. Open circles show individual data and regression lines are represented with 95% CIs.
Ijms 25 12194 g002
Figure 3. Means (±SE) of the stress biomarkers in in-nest or heat-challenged nestlings. Individuals were exposed to prenatal control-calls (black circles) or heat-calls (red circles) and, at day 13 post-hatch, blood samples were collected from individuals in their nest (in-nest) or following a heat-challenge. The panels show levels of (A) heat shock cognate 70 (HSC70; N = 73), (B) heat shock protein 90α (HSP90α; N = 73), (C) the heterophil-to-lymphocyte (H/L) ratio (N = 44) and (D) corticosterone (CORT; N = 53). Removing the high HSC70 value (A) or the high CORT value (D) did not change the results. Open circles show individual data and solid circles show group means ± SE.
Figure 3. Means (±SE) of the stress biomarkers in in-nest or heat-challenged nestlings. Individuals were exposed to prenatal control-calls (black circles) or heat-calls (red circles) and, at day 13 post-hatch, blood samples were collected from individuals in their nest (in-nest) or following a heat-challenge. The panels show levels of (A) heat shock cognate 70 (HSC70; N = 73), (B) heat shock protein 90α (HSP90α; N = 73), (C) the heterophil-to-lymphocyte (H/L) ratio (N = 44) and (D) corticosterone (CORT; N = 53). Removing the high HSC70 value (A) or the high CORT value (D) did not change the results. Open circles show individual data and solid circles show group means ± SE.
Ijms 25 12194 g003
Figure 4. Body mass of heat-challenged nestlings in relation to their (A) heterophil-to-lymphocyte (H/L) ratio (N = 21) and (B) heat shock cognate 70 (HSC70) levels (N = 38). In (B), the p-value is marginal (p = 0.051) when the high HSC70 value is excluded. Regression lines are represented with 95% CIs.
Figure 4. Body mass of heat-challenged nestlings in relation to their (A) heterophil-to-lymphocyte (H/L) ratio (N = 21) and (B) heat shock cognate 70 (HSC70) levels (N = 38). In (B), the p-value is marginal (p = 0.051) when the high HSC70 value is excluded. Regression lines are represented with 95% CIs.
Ijms 25 12194 g004
Table 1. Linear mixed models 1 fitting levels of heat shock protein 90α (HSP90α), heat shock cognate 70 (HSC70), corticosterone (CORT) or the heterophil-to-lymphocyte (H/L) ratio as a function of the prenatal playback (heat-call or control-call), nest temperature (12D-Tnest 2) and their interaction in zebra finch nestlings.
Table 1. Linear mixed models 1 fitting levels of heat shock protein 90α (HSP90α), heat shock cognate 70 (HSC70), corticosterone (CORT) or the heterophil-to-lymphocyte (H/L) ratio as a function of the prenatal playback (heat-call or control-call), nest temperature (12D-Tnest 2) and their interaction in zebra finch nestlings.
Response VariableFixed EffectEst. 3SE 4t Valuep-Value 5
In-nest condition
HSP90α (N = 35) *intercept1.270.0342.93<0.001
playback (heat-calls)−0.040.04−0.880.388
12D-Tnest0.080.032.410.022
playback × 12D-Tnest−0.090.04−2.040.050
HSC70 (N = 35)intercept1.730.0442.32<0.001
playback (heat-calls)0.070.051.310.202
12D-Tnest0.080.041.780.085
playback × 12D-Tnest−0.150.06−2.760.010
CORT (N = 25) *intercept0.680.097.58<0.001
playback (heat-calls)0.020.120.180.860
12D-Tnest0.010.060.180.860
H/L ratio (N = 23) *intercept−0.230.06−3.930.001
playback (heat-calls)−0.160.08−1.940.066
12D-Tnest0.050.041.090.289
Heat-challenge condition
HSP90α (N = 38)intercept1.260.0433.62<0.001
playback (heat-calls)0.020.0525.140.656
12D-Tnest0.010.0320.750.799
HSC70 (N = 38) *intercept1.780.0444.77<0.001
playback (heat-calls)0.000.060.070.949
12D-Tnest−0.020.03−0.590.557
CORT (N = 28) *intercept0.680.079.72<0.001
playback (heat-calls)−0.030.10−0.270.789
12D-Tnest−0.020.05−0.490.632
H/L ratio (N = 21)intercept−0.100.10−1.0570.305
playback (heat-calls)−0.070.14−0.5120.615
12D-Tnest0.040.080.4850.634
* The random factor brood identity (brood-ID) was removed because its variance was null. 1 full model: response ~ prenatal playback + 12D-Tnest + playback × 12D-Tnest + (1|brood-ID); 2 12D-Tnest: mean daytime nest temperature experienced from hatching to day 12; 3 Est.: estimate; 4 SE: standard error; 5 p-values below 0.05 (alpha level threshold) are significant and indicated in bold.
Table 2. Full linear models of nestling mass as a function of heat shock cognate 70 (HSC70) and heat shock protein 90α (HSP90α), or corticosterone (CORT) and the heterophil-to-lymphocyte (H/L) ratio.
Table 2. Full linear models of nestling mass as a function of heat shock cognate 70 (HSC70) and heat shock protein 90α (HSP90α), or corticosterone (CORT) and the heterophil-to-lymphocyte (H/L) ratio.
Physiological PredictorFixed EffectEst. 1SE 2t Valuep-Value 3Variance
[brood-ID 4] (±SE)
In-nest condition
HSC70 and HSP90α (N = 35) Intercept9.260.2143.82<0.0010.722 (±0.850)
HSC70−0.270.18−1.480.151
HSP90α0.040.170.260.801
CORT and H/L ratio (N =17) Intercept9.520.2834.14<0.0010.674 (±0.822)
CORT0.200.260.750.464
H/L ratio−0.290.23−1.240.252
Heat-challenge condition
HSC70 and HSP90α (N = 35) Intercept9.450.1755.55<0.0010.121 (±0.348)
HSC700.370.172.250.031
HSP90α−0.280.17−1.640.109
CORT and H/L ratio (N = 17) Intercept9.680.2637.80<0.001null (removed)
CORT−0.020.28−0.090.930
H/L ratio−0.540.28−1.950.072
1 Est.: estimate; 2 SE: standard error; 3 p-values below 0.05 (alpha level threshold) are significant and indicated in bold; 4 brood-ID: brood identity.
Table 3. Linear models 1 fitting each stress biomarker of zebra finch nestlings as a function of the experiment (in-nest or heat-challenge) in a pooled dataset.
Table 3. Linear models 1 fitting each stress biomarker of zebra finch nestlings as a function of the experiment (in-nest or heat-challenge) in a pooled dataset.
Response VariableFixed EffectEst. 2SE 3tp-Value 4Variance
[brood-ID 5] (±SE)
CORT 6 (N = 53)intercept0.700.0512.92<0.001null (removed)
experiment (heat-challenge)−0.030.07−0.370.712
HSC70 7 (N = 73)intercept1.760.0358.26<0.0010.003 (±0.05)
experiment (heat-challenge)0.020.040.530.596
HSP90α 8 (N = 73)intercept1.240.0348.75<0.001null (removed)
experiment (heat-challenge)0.030.040.950.344
H/L ratio 9 (N = 44)intercept−0.330.06−5.70<0.0010.03 (±0.19)
experiment (heat-challenge)0.210.082.760.009
1 model: response ~ experiment + (1|brood-ID); 2 Est.: estimate; 3 SE: standard error; 4 p-values below 0.05 (alpha level threshold) are significant and indicated in bold; 5 brood-ID: brood identity; 6 CORT: corticosterone; 7 HSC70: heat shock cognate 70; 8 HSP90α: heat shock protein 90α; 9 H/L ratio: heterophil-to-lymphocyte ratio.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Udino, E.; Pessato, A.; Addison, B.; Crino, O.L.; Buchanan, K.L.; Mariette, M.M. Prenatal Acoustic Signals Influence Nestling Heat Shock Protein Response to Heat and Heterophil-to-Lymphocyte Ratio in a Desert Bird. Int. J. Mol. Sci. 2024, 25, 12194. https://doi.org/10.3390/ijms252212194

AMA Style

Udino E, Pessato A, Addison B, Crino OL, Buchanan KL, Mariette MM. Prenatal Acoustic Signals Influence Nestling Heat Shock Protein Response to Heat and Heterophil-to-Lymphocyte Ratio in a Desert Bird. International Journal of Molecular Sciences. 2024; 25(22):12194. https://doi.org/10.3390/ijms252212194

Chicago/Turabian Style

Udino, Eve, Anaïs Pessato, BriAnne Addison, Ondi L. Crino, Katherine L. Buchanan, and Mylene M. Mariette. 2024. "Prenatal Acoustic Signals Influence Nestling Heat Shock Protein Response to Heat and Heterophil-to-Lymphocyte Ratio in a Desert Bird" International Journal of Molecular Sciences 25, no. 22: 12194. https://doi.org/10.3390/ijms252212194

APA Style

Udino, E., Pessato, A., Addison, B., Crino, O. L., Buchanan, K. L., & Mariette, M. M. (2024). Prenatal Acoustic Signals Influence Nestling Heat Shock Protein Response to Heat and Heterophil-to-Lymphocyte Ratio in a Desert Bird. International Journal of Molecular Sciences, 25(22), 12194. https://doi.org/10.3390/ijms252212194

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop