Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Physiological Condition of Juvenile Wading Birds in Relation to Multiple Landscape Stressors in the Florida Everglades: Effects of Hydrology, Prey Availability, and Mercury Bioaccumulation

  • Garth Herring ,

    gherring@usgs.gov

    Affiliation United States Geological Survey, Forest and Rangeland Ecosystem Science Center, Corvallis, Oregon, United States of America

  • Collin A. Eagles-Smith,

    Affiliation United States Geological Survey, Forest and Rangeland Ecosystem Science Center, Corvallis, Oregon, United States of America

  • Dale E. Gawlik,

    Affiliation Department of Biological Sciences, Florida Atlantic University, Boca Raton, Florida, United States of America

  • James M. Beerens,

    Affiliation Department of Biological Sciences, Florida Atlantic University, Boca Raton, Florida, United States of America

  • Joshua T. Ackerman

    Affiliation United States Geological Survey, Western Ecological Research Center, Dixon Field Station, Dixon, California, United States of America

Abstract

The physiological condition of juvenile birds can be influenced by multiple ecological stressors, and few studies have concurrently considered the effects of environmental contaminants in combination with ecological attributes that can influence foraging conditions and prey availability. Using three temporally distinct indices of physiological condition, we compared the physiological response of nestling great egrets (Ardea alba) and white ibises (Eudocimus albus) to changing prey availability, hydrology (water depth, recession rate), and mercury exposure in the Florida Everglades. We found that the physiological response of chicks varied between species and among environmental variables. Chick body condition (short-term index) and fecal corticosterone levels (medium-term) were influenced by wetland water depth, prey availability, region, and age, but not by mercury contamination. However, mercury exposure did influence heat shock protein 70 (HSP70) in egret chicks, indicating a longer-term physiological response to contamination. Our results indicate that the physiological condition of egret and ibis chicks were influenced by several environmental stressors, and the time frame of the effect may depend on the specialized foraging behavior of the adults provisioning the chicks.

Introduction

Environmental stressors, such as contaminants, prey availability, and adverse weather, can elicit strong physiological responses from wild birds in order to help them overcome short-term challenges. These responses can be particularly harmful during the early stages of life if the response exceeds physiological norms [1][3]. Physiological [4] and behavioral [5] responses can manifest throughout the annual cycle by many avian species in order to cope with seasonal environmental stochasticity. Understanding avian responses to multiple environmental stressors, when one of those stressors is an environmental contaminant, is poorly understood, yet this situation is prevalent in many anthropogenically-impacted environments [6], [7].

Mercury (Hg) contamination and subsequent bioaccumulation in waterbirds is problematic in wetlands throughout the world. Inorganic Hg is converted into the toxic and bioaccumulative form methylmercury (MeHg) under biogeochemical conditions that are common in wetland ecosystems [6], [7]. Once methylated, MeHg biomagnifies as it is transferred through the food chain [8][10]. Environmentally relevant Hg concentrations in birds have been associated with a suite of impaired physiological and reproductive responses, including altered reproductive hormone levels [11], [12], compromised hypothalamic-pituitary-adrenal axis [13][15], and reduced reproductive performance [16].

The decline of several species of Everglades' wading bird between the 1930s and 2001 has been related to changes in prey availability and alterations to the ecosystem's hydrology [17][19] and possibly due to Hg exposure [20], [21]. Great egrets (Ardea alba; hereafter egret) and white ibis (Eudocimus albus; hereafter ibis) are two of the most common breeding waterbirds in the Everglades, and they exhibit different foraging strategies that influence their exposure to ecological stressors. Ibis are more limited in their use of habitats and available prey than egrets [22]. Specifically, Ibis select high-quality foraging patches that they then abandon relatively quickly once prey availability drops, whereas egrets tend to remain at foraging sites even when prey densities are low [22]. Although both species select different foraging habitats [23] nest survival for both species is similarly influenced by hydrological conditions and prey availability [24]. Furthermore, Hg concentrations differ between egrets and ibises in the Florida Everglades [10], [25][28], which is likely a function of their foraging ecology and prey selection [28].

Herein, we assessed the physiological response of juvenile egrets and ibis to changing prey availability, hydrology (water depth, recession rate), and mercury exposure. We measured physiological biomarkers and body condition for egret and ibis nestlings in two consecutive years that differed greatly in hydrologic conditions, prey availability, and Hg exposure. Additionally, we measured physiological biomarkers of environmental stressors that manifest across a range of different time frames. In particular, we used measures of physiological condition that we expected would represent a temporal continuum of responses relative to the age of wading bird chicks in our study, including a body-condition index (short-time frame; 1–2 days; [29]), fecal corticosterone metabolites [FCORT] (medium-time frame; 2–7 days; [30], [31]), and heat shock proteins 60 and 70 [HSP60, HSP70] (long-time frame: weeks - >1 month; [28], [32][34]). Corticosterone (CORT) is a hormone that serves as a signal to modify both behavior and metabolism during a period of acute stress. Corticosterone is released by the hypothalamo-pituitary-adrenal axis into the blood stream when birds experience acute stress, facilitating a rapid response to overcome the stressor [2], [35]. Variation in CORT levels in waterbirds have been correlated with a variety of factors including prey availability [31], [36], Hg exposure [12], [15], and changing hydrological conditions [37]. In contrast, heat shock proteins are highly conserved molecular chaperones that function to maintain optimal cell function and homeostasis [38], [39] by being amplified through up-regulation during periods of stress to minimize cell protein damage [40], [41]. Elevated heat shock proteins are commonly associated with responses to damage associated with heavy metals including Hg [42], [43], decreased food availability [31], and rapid changes in hydrological conditions at foraging sites for birds [37]. In order to better understand the potential conservation implications of exposure to multiple environmental stressors we simultaneously tested those effects on these physiological biomarkers and chick body condition.

Methods

Ethics Statement

We made every attempt to reduce disturbance, stress, and other impacts to target and non-target species during the course of this research. Florida Atlantic University Institutional Animal Care and Use Committee (Protocol A0534) approved this study, and we conducted the research under U.S. Geological Survey banding permit 23354 and Florida Fish and Wildlife Conservation Commission Scientific Research Permit WX04487.

Study Area

In 2006 and 2007, we monitored egret and ibis nestlings in Florida Everglades wading bird colonies at the Arthur R. Marshall Loxahatchee National Wildlife Refuge (Lox; n = 8 colonies) in Palm Beach County, Water Conservation Area (WCA) 2A (n = 2 colonies) in Broward County, and WCA 3A (n = 5 colonies) in Broward and Miami-Dade Counties, Florida, USA (Figure S1). Colonies were located by using a combination of ground and aerial surveys [24], by tracking radio-tagged adults [23], and by visiting known colony sites (e.g., South Florida Wading Bird Reports, South Florida Water Management District). Although we attempted to sample colonies in which both species nested, only one species nested in five of the twelve colonies.

Chick sampling

We surveyed nests within egret and ibis colonies that were associated with our radio-tagged adults ([23], [24], Figure S1). We revisited nests and resampled chicks (n = 61 egret and 65 ibis nests) approximately every 7 days, and determined the hatching order and post-hatch age of chicks.

We used FCORT as an index of stress because this approach minimizes handling-induced bias that results during bird sampling. We collected fecal material as described in [37], and fecal samples were stored in labeled cryovial tubes and placed on ice. We recorded mass for both species to the nearest 5 g using a spring scale and recorded all other morphometric measurements to the nearest 1 mm. We sampled all nestlings during the morning to minimize heat stress and time of day effects on physiological metrics.

We sampled blood from chicks to measure both HSPs and Hg concentrations. From each chick, we collected up to 2 ml of blood from the brachial vein using a 27-gauge needle. Blood samples were stored on ice in heparinized cryovials until transport to the laboratory. Following The Ornithological Council guidelines [44], we limited the total blood volume of blood collection to <1% of bird mass and thus had limited blood volume available for determining physiological parameters and Hg concentrations. Therefore, we separated plasma and RBC fractions in the laboratory via centrifugation (15 min, 5,000 rpm) and determined HSPs and THg concentrations in RBCs. Each fraction was stored frozen at −20°C until analyses could be determined.

Laboratory Methods

Fecal samples were homogenized in the laboratory, divided into 2 equal 1-ml wet portions, and dried using a Labconco CentriVap Concentro (Labconco, Kansas City, MO). Dried samples (∼0.25 g) were then mixed with 5 ml of 95% ethanol and vortexed for 30 min. After centrifugation (15 min, 2500 × g) the supernatant was transferred to a new vial, and evaporated under a stream of nitrogen gas. Corticosterone metabolites were then resuspended in diluted extraction buffer and measured using the Correlate-EIA Corticosterone Enzyme Immunoassay Kit (EIA [45];) following the manufacturer's instructions (Assay Design, Inc., Ann Arbor, MI). Inter- and intra-assay coefficients of variation for FCORT internal standards were 7% and 11%, and 8% and 9% respectively for egrets and ibises. We validated that FCORT levels did not change after freezing as in the case of mammals [46] by freezing and measuring FCORT levels monthly for 6 months [45].

A subset of RBCs was washed three times using phosphate buffered saline, centrifuged and the pellet was removed after the final wash. The RBC pellet was then mixed with 1× extraction reagent and a protease inhibitor cocktail (Sigma), vortexed for 5 min and then sonicated for 1 min. Samples were again centrifuged (15 min, 2500 g) and the supernatant removed. We measured HSP60 (HSPD1) and HSP70 (HSP72/HSPA1A) in the supernatant using EIA kits specific to just those stress proteins and not all other HSP60 and HSP70 family members EIA kits (Assay Designs, Inc., Ann Arbor, MI). Inter- and intra-assay coefficients of variation for HSP60 and HSP70 internal standards were 5% and 7% and 6% and 7% respectively for egrets and ibis. All samples were run in duplicate, and means of duplicates were used in all analyses. All EIA kits were validated using serial dilutions and spike tests to determine percent recovery [45].

Mercury Determination

We analyzed the remaining portion of RBCs for total mercury (THg) concentrations because the majority of Hg in blood is in the RBC fraction [47], [48] and RBC Hg concentrations are correlated with whole blood Hg concentrations [49], [50]. Further, almost all (>95%) of the mercury in avian blood is in the methylmercury form and THg concentrations are highly correlated with methylmercury levels [51], [52]. To minimize any potential variability in THg concentrations associated with processing and differential moisture loss during long-term storage of samples, RBC fractions were dried (48 hours at 50°C) to a constant weight before THg analysis. We measured up to 0.05 g of dried blood into nickel sample vessels (weighed to the nearest 0.00001 g; Mettler Toledo model XS105, Mettler Toledo, Columbus, Ohio, USA). Following U.S. Environmental Protection Agency method 7473 [53], we determined THg concentrations in each RBC sample on a Milestone DMA-80 direct mercury analyzer (Milestone, Monroe, Connecticut, USA). Analytical equipment was calibrated using certified standard solutions prior to analysis, and accuracy and precision were evaluated within each analytical batch through the inclusion of certified reference materials (either dogfish muscle tissue [DORM-3] or dogfish liver [DOLT-4] by the National Research Council of Canada, Ottawa, Canada), calibration verifications (liquid standards), duplicates, and blanks. Recoveries averaged 102.4±0.06% (n = 17) and 97.7±0.04% (n = 33) for certified reference materials and calibration checks, respectively. Absolute relative percent difference for all duplicates averaged 1.7%±0.08% (n = 16).

Hydrological Conditions

We used the Everglades Depth Estimation Network (EDEN [54];) to estimate water depths and water level recession rates within foraging ranges of adults associated with breeding colonies where chicks were sampled. The EDEN uses a system of water level gauges to produce a water surface model coupled with a ground elevation model to estimate water depth for the entire freshwater portion of the Greater Everglades [55]. The EDEN calculates water level depths in 400 m×400 m grid cells at daily time steps accounting for evapotranspiration, rainfall, and sheet flow. The estimated water depths have been found to be accurate to within 5 cm [55].

To estimate recession rates and water depths, we used protocols developed during previous studies [23], [24], [28]. Briefly we used fixed radius buffers (hereafter foraging range) associated with breeding colonies based on mean distances flown plus one standard error by radio tagged adult egrets and ibises in a concurrent study. We defined foraging range as the average distance that radio tagged adults flew from breeding colonies to a foraging site throughout the breeding season. We extracted the water depth from each EDEN grid cell within the foraging ranges for each specific colony, year, and species combination. For each day that we sampled chicks, we calculated water recession rate and water depth as the mean value for the preceding 2-week period [24], [28]. Positive recession rates indicate decreasing water levels. Surface water could still be present when estimated cell depth <0 cm; however, the 5-cm error in EDEN depths suggests that a majority of the ground surface would be exposed in a given 400 m×400 m grid cell.

Statistical procedures

We used linear mixed-effects models in JMP using maximum likelihood estimation (JMP 2001 [56];) to test the effects of environmental stressors (water depth, recession rate, prey availability, and Hg) on the physiological condition (chick body condition index, FCORT, HSP70, and HSP60) of egret and ibis chicks. We included age of chicks in models to account for variation in prey selection behavior as chicks aged [57][59]. The hatch order of chicks within the nest was included to account for competition for food within nests [24], [28]. Lastly, we included region because our previous research found differences in adult physiological condition among Everglades' regions [37]. We did not include a year effect because it was confounded with large differences in prey biomass observed between years. For egret chick models, we statistically nested the chick identification code within the nest identification code and added it as a random effect to the model to account for the non-independence of sampling chicks multiple times and different numbers of chicks within the same nest. For ibis chick models, we included the nest identification code as a random effect to account for the non-independence of sampling multiple chicks within a nest. We used the maximum likelihood estimation method in our linear models. Because we intentionally studied two species with contrasting foraging strategies, we ran separate models for each species. Natural log transformations where used to improve normality of residuals for prey biomass, THg, FCORT, HSP60, and HSP70 data. We assessed the level of multicollinearity using the variance inflation factor (VIF; 60); across all models the VIFs ranged from 1.04–2.41 suggesting no significant multicollinearity. We calculated a marginal and conditional coefficient of determination (R2) for each model to examine goodness-of-fit in R [GLMM R2; 61]. Briefly, for mixed effects models, the marginal R2 (R2GLMM(m)) reports the goodness-of-fit for just the fixed effects, while the conditional R2 (R2GLMM(c)) incorporates both the random and fixed effects as a measure of explained variation [61].

In the case of chick mass, we did not have repeated measurements of all chicks across time and subsequently could not calculate growth rates. Instead, we calculated a chick body condition index that allowed us to examine the potential effects of multiple landscape stressors on a relative measure of the condition of the chick that was independent of the mass of a chick [29], [62]. To do so, we first we calculated a structural body size index using principal components analysis on the morphometric measurements (e.g. tarsus, bill length, wing length) for each species separately. We then regressed chick mass on tarsus length (the variable our principal components analysis identified as the best predictor of structural size for both egrets and ibis), including the first principal component as a covariate. We then used the residuals expressed as a percentage of the predicted mass, to obtain our chick body condition index [29], [62]. Benson et al. [29] found that measurements taken during a single visit to a colony were useful to detect variation in chick body condition across time.

Results

During the 2006 and 2007 nesting seasons we collected blood samples from 254 chicks (169 egret; 2006 = 27, 2007 = 142, and 85 ibis; 2006 = 30, 2007 = 55). The geometric mean across all dates, regions, and hatch order was 6.99±0.63 (SE) in egrets and 4.57±0.83 in ibises for FCORT concentration (µg/g), 13.73±0.96 in egrets and 14.73±2.06 in ibis for HSP60 concentration (ng/ml), 4.44±0.44 in egrets and 2.53±0.35 in ibis for HSP70 concentration (ng/ml), and 4.13±0.16 in egrets and 0.58±0.07 in ibis for THg concentration (µg/g dw).

Chick body condition index

Egret chick body condition was influenced by water depth (F1, 128.8 = 6.60, P = 0.01), but not region (F2, 114.7 = 2.30, P = 0.10), water recession rate (F1, 119.2 = 2.28, P = 0.14), prey biomass (F1, 142.7 = 0.38, P = 0.53), age (F1, 52.4 = 0.43, P = 0.51), hatch order (F1, 79.9 = 1.95, P = 0.17), nor Hg concentration (F1, 150.8 = 0.02, P = 0.88). The low marginal R2 value of our egret chick body condition model (R2GLMM(m) = 0.10, R2GLMM(c) = 0.21) suggests that the model was a relatively poor predictor of body condition. However, egret chick body condition increased with increasing water depth (β = 0.20±0.13–95% confidence interval) on average by 23% across the range of water depths observed (range = −46 to 17 cm; Figure 1A).

thumbnail
Figure 1. Relationship between partial residuals of chick body condition index and water depth (A; cm/2) for great egret chicks and age (B; days) for white ibis chicks in the Florida Everglades, after accounting for other model variables.

Dashed lines indicated 95% confidence.

https://doi.org/10.1371/journal.pone.0106447.g001

Conversely, ibis chick body condition was influenced by both region (F1, 37.6 = 9.09, P = 0.004) and age (F1, 75 = 5.26, P = 0.02), but not by water depth (F1, 68.8 = 0.0001, P = 0.99), water recession rate (F1, 72.5 = 0.005, P = 0.98), prey biomass (F1, 55.6 = 0.06, P = 0.80), hatch order (F1, 38.6 = 33, P = 0.56), nor Hg concentration (F1, 75 = 0.10, P = 0.75). The ibis chick body condition model marginal R2 value (R2GLMM(m) = 0.23), suggested that the model was a modest predictor of ibis chick body condition. Similarity between the marginal R2 value and the conditional R2 (R2GLMM(c) = 0.23) suggested that little additional variance was explained by the random effect associated with sampling multiple chicks within nests. Ibis chick body condition was 1.9 times higher in Lox (mean = 5.36±2.68 SE) than WCA3A Lox (−10.02±4.31 SE; β = 7.69±4.99) and decreased by 26% across the range of ages (β = −0.64±0.54, range = 3–30 days; Figure 1B).

Fecal corticosterone metabolites

Egret FCORT metabolite concentrations were influenced by age (F1, 13.5 = 13.50, P<0.001), water depth (F1, 7.5 = 7.50, P = 0.007), and region (F1, 108.6 = 3.06, P = 0.007), but not by water recession rate (F1, 106.8 = 0.01, P = 0.97), prey biomass (F1, 124.9 = 2.82, P = 0.09), hatch order (F1, 80.5 = 0.78, P = 0.38), nor Hg concentration (F1, 133.6 = 1.65, P = 0.20). The similarity between the marginal and conditional R2 values (R2GLMM(m) = 0.34, R2GLMM(c) = 0.37) suggests that the fixed effects alone were considerably more influential than the random effects associated with sampling multiple chicks within the same nest on more than one occasion. Egret FCORT metabolite concentrations decreased with increasing age (β = −0.04±0.02) and water depth (β = −0.02±0.01) on average by 79% across the range of ages (range = 3–42; Figure 2A) and 72% across the range of water depths (range = −46 to 17 cm; Figure 2B). Egret FCORT metabolite concentrations were higher in WCA3A (11.74±2.45 SE) than either WCA2A (t118.3 = 2.09, P = 0.04; mean = 4.71±1.60 SE) or Lox (t102.3 = 2.25, P = 0.03; mean  = 6.34±0.76 SE) respectively, but similar between WCA2A and Lox (t114.5 = −0.83, P = 0.40).

thumbnail
Figure 2. Relationship between partial residuals of fecal corticosterone metabolite (µg/g dw) and chick age (A; days) and water depth (B; cm) for great egret chicks and prey biomass (C; g/m2) for white ibis in the Florida Everglades, after accounting for other model variables.

Dashed lines indicated 95% confidence interval.

https://doi.org/10.1371/journal.pone.0106447.g002

Similar to egrets, region influenced ibis FCORT metabolite concentrations (F1, 54.8 = 4.40, P = 0.04), however unlike egrets, prey biomass also influenced ibis FCORT metabolites (F1, 58.4 = 19.41, P<0.001). Ibis FCORT metabolites were not influenced by water depth (F1, 59.6 = 0.65, P = 0.42), water recession rate (F1, 60.4 = 0.36, P = 0.55), age (F1, 60.8 = 0.25, P = 0.62), hatch order (F1, 46.1 = 1.43, P = 0.23), nor Hg concentration (F1, 59.1 = 0.001, P = 0.96). The ibis FCORT model had a moderate marginal R2 value (R2GLMM(m) = 0.36), but the conditional R2 value (R2GLMM(c) = 0.71) suggested that a similar amount of the variability in ibis FCORT concentrations were associated with just the random effect of sampling multiple chicks within the same nest. FCORT metabolite concentrations were 3.3 times higher in Lox (mean = 5.26±1.68 SE) than WCA3A (1.60±0.70 SE; Β = 0.59±0.54), and decreased by 98% (Β = −1.40±0.62), across the range of prey biomass (range = 1.8–26.7 g/m2; Figure 2C).

Heat shock protein 70

Egret HSP70 concentrations were influenced by Hg concentration (F1, 109.9 = 4.09, P = 0.04), but not region (F2, 100.3 = 1.56, P = 0.21), water depth (F1, 101.2 = 0.16, P = 0.68), water recession rate (F1, 91.0 = 0.37, P = 0.54), prey biomass (F1, 99.36 = 0.23, P = 0.62), age (F1, 110.8 = 0.02, P = 0.87), nor hatch order (F1, 76.6 = 1.63, P = 0.20). The overall amount of variation in egret HSP70 explained by the fixed effects was relatively low (R2GLMM(m) = 0.13), however very little of that variation was explained by our random effects (R2GLMM(c) = 0.16) associated sampling multiple chicks from the same nest or sampling those chicks on several occasions. The beta coefficient estimate for Hg (β = −0.48±0.22) indicated a decrease in HSP70 concentrations of 83% across the range Hg values (range = 0.26–11.01 µg/g dw; Figure 3).

thumbnail
Figure 3. Partial residuals of HSP70 concentrations (ng/ml) and red blood cell THg (µg/g dw) for great egret chicks in the Florida Everglades, after accounting for other model variables.

Dashed lines indicated 95% confidence interval.

https://doi.org/10.1371/journal.pone.0106447.g003

In the case of ibis, neither water depth (F1, 30.9 = 0.12, P = 0.72), region (F1, 28.6 = 0.11, P = 0.73), water recession rate (F1, 32.9 = 0.0009, P = 0.97), prey biomass (F1, 32.2 = 0.12, P = 0.72), age (F1, 31.4 = 0.49, P = 0.48), hatch order (F1, 29.5 = 0.02, P = 0.87), nor Hg concentration (F1, 13.9 = 0.0001, P = 0.99) influenced HSP70 concentrations. Correspondingly, the fixed effects (R2GLMM(m) = 0.06) explained very little of the variation in HSP70, however the conditional effects (R2GLMM(c) = 0.74) associated with sampling chicks from the same nest played a considerably larger role in ibis HSP70 concentrations.

Heat shock protein 60

Egret chick HSP60 concentrations were not influenced by water depth (F1, 118.9 = 0.0002, P = 0.98), region (F2, 111.4 = 0.01, P = 0.98), water recession rate (F1, 119.7 = 0.27, P = 0.60), prey biomass (F1, 86.6 = 0.91, P = 0.34), age (F1, 134.5 = 1.81, P = 0.18), hatch order (F1, 81.7 = 0.05, P = 0.81), nor Hg concentration (F1, 120.2 = 0.86, P = 0.34). The low R2 values associated with both the fixed and random effects suggested that variation in egret HSP60 results from factors other than those measured in this study (R2GLMM(m) = 0.04, R2GLMM(c) = 0.04).

In the case of Ibis HSP60 concentrations were influenced by age (F1, 43.4 = 4.32, P = 0.04), but not prey biomass (F1, 45.5 = 0.89, P = 0.35), water depth (F1, 42.0 = 0.11, P = 0.74), water recession rate (F1, 46.0 = 0.04, P = 0.83), hatch order (F1, 31.6 = 0.02, P = 0.87), nor Hg concentration (F1, 46.6 = 1.61, P = 0.21). The overall amount of variation in ibis HSP60 explained by the fixed effects was relatively low (R2GLMM(m) = 0.15), but the random effect (R2GLMM(c) = 0.89) associated with sampling multiple chicks from the same nest explained considerably more of the variation in ibis HSP60 concentrations. The beta coefficient estimates for age (β = −0.05±0.04) indicated HSP concentrations declined by 74% across the range of ages (range = 3–30 days; Figure 4).

thumbnail
Figure 4. Relationship between partial residuals of HSP60 (ng/ml) and chick age (days) for white ibis chicks in the Florida Everglades, after accounting for other model variables.

Dashed lines indicated 95% confidence interval.

https://doi.org/10.1371/journal.pone.0106447.g004

Discussion

Different environmental stressors were associated with different time frames and biomarkers of physiological condition in egret and ibis chicks. The short and moderate time frame markers, chick body condition and FCORT, respectively, were correlated with landscape variables such as prey biomass and water depth, whereas there was no measurable effect of Hg. In contrast, the long-term physiological markers were not influenced by any of the measured landscape stressors, but one (HSP70 in egrets) was influenced by Hg concentration. One potential confounding aspect in our study is the fact that the same hydrology variables used to define the landscape stressors that we tested to determine their influence on physiological condition are also know to influence concentrations of Hg in egret and ibis adults and chicks [28]. Previously, Herring et al. [28] found that water depths and recession rates, while not the primary variables that influenced Hg concentrations, did play a secondary role in influencing Hg concentrations in egrets and ibis. Disentangling those effects is beyond the scope of this study, but illustrates the complexity of examining different environmental stressors on the physiological condition of birds.

Egret HSP70 decreased by 83% across the range (low-high) of Hg concentrations, suggesting that Hg may have down-regulated these heat shock proteins. Previous research on heat shock protein regulation in relation to both inorganic Hg and MeHg have found that heat shock protein regulation can be either inhibited (down-regulated) or induced (up-regulated) depending on the chemical form of Hg [63]. Up-regulation of heat shock proteins are associated with protective stress responses to minimize cell protein damage [40], [41], [64]. Methylmercury exposure in birds is associated with oxidative stress in liver, kidney, and brain tissue [65], and up-regulation of heat shock proteins has recently been identified as a potential biological defense mechanism against MeHg toxicity [66], [67]. However, it is unclear what the consequences are for down-regulation of heat shock proteins as we observed with egrets. Presumably, a reduction of heat shock proteins could result in decreased cell function and homeostatic imbalance, given that this is the primary role heat shock proteins when an individual is not in a stressful situation [38], [39]. Opposite response patterns of other stress biomarkers in relation to Hg exposure have been documented in other birds (e.g., [15]). Herring et al. [15] found that Hg exposure in Forster's tern (Sterna forsteri) chicks depressed FCORT concentrations by up to 81% via the downregulation of the HPA axis. These results suggest that the physiological stress response systems of birds are dynamic and that Hg exposure can negatively influence critical physiological processes.

Current Hg exposure levels of egrets and ibises throughout much of the Florida Everglades are low relative to historical reports [28], [68]. However, the fact that for egrets, one of their HSPs was influenced by Hg, demonstrates that the current Hg levels could still influence the physiology of this species. Moreover, Hg hotspots still exist within the Florida Everglades (e.g., interior portions of Everglades National Park not sampled in this study) where Hg levels in fish far exceed levels measured anywhere else within the ecosystem, and still pose a risk to wading birds [10], [21]. Ibis heat shock proteins were not influenced by Hg exposure in this study, which was likely a function of Hg exposure levels being approximately 7 times lower than those of egrets.

We found that ibis but not egret FCORT levels were influenced by prey biomass levels, and that they declined by 98% across low-to-high range of prey biomass at foraging sites. This is consistent with responses in pre-breeding adult birds [37], suggesting that this relationship may be consistent across life stages. Similarly, ibis nest success was also more sensitive to changes in prey biomass than that of egrets [24], demonstrating the importance of prey biomass to ibis. Egret FCORT metabolite concentrations also declined by 79% across the age range (young to old) of chicks. Other studies of CORT stress levels in chicks have found that concentrations often increase as chicks get older, in preparation for fledging [3], [69], [70]. However, there is some evidence that elevated CORT levels can result in decreased protein deposition in growing feathers [71], which may lead to poorer feather quality [72]. Increased baseline CORT levels can also impair growth rates and immune function [73]. Alternatively egret CORT levels may be higher in younger chicks to facilitate improved begging to adults [74], [75] and promote growth at a critical stage in the chicks' development since CORT allows the organism to mobilize needed energy reserves [76].

We found support for the influence of hydrological variables on egret physiology, in particular egret chick body condition increased with increasing water depth. Critical to understanding these responses is the fact that the range of water depths associated with our study did not exceed the observed foraging range and optimal depths for egrets [23] and generally would be considered to be shallow for egret foraging [22]. As such, the increase in chick body condition is likely driven by the fact that egrets prefer larger prey items [57] found in deeper depths relative to those very shallow sites where large prey die rapidly. We expect that body condition would have dropped had depths increased beyond the normal foraging range of egrets.

Neither water depths nor recession rates influenced ibis chick physiological condition in the current study. Although recession rates were found to influence the pre-breeding physiological condition of ibis previously [37], recession rates during the nestling stage may play a smaller role in chick physiological condition because water depths have largely reached optimal foraging depths for adults around breeding colonies [22], [23]. Thus, small variation in water depth and recession rate likely have little influence on accessibility of foraging sites or the availability of prey for ibis. However, if rainfall events occur, rapid negative recession rates that result in increased water depths may result in the loss of foraging site access and reduced prey availability. Although these conditions did not occur in the present study, they can have a substantial effect on nestlings [77] and would be expected to decrease chick physiological condition. Similarly, water depths were largely within the foraging range and were near optimal within the proximity of breeding colonies during this study. In fact, water depths associated with ibises were less variable than those observed in egrets (coefficient of variation  = 20% vs 56% respectively). Under conditions with deeper or more variable water depths, we expect that ibises may have an increased physiological response due to their dependence on high quality (shallow water with high prey densities) foraging sites.

We also found that chick age was an influential variable in explaining physiological responses, including ibis body condition and HSP60, and egret FCORT. The body condition of ibis chicks declined 26% from the youngest to the oldest chicks. Although it is unclear why chick body condition declined as they aged, there was not a concomitant increase in stress biomarkers such as FCORT or HSPs. In fact, HSP60 levels actually declined by 74% across the range of ibis chick ages. One previous study found that chickens (Gallus gallus domesticus) that were experimentally fed lower energy diets, resulting in lower body masses, also had lower levels of HSPs relative to those with higher energy diets and higher masses [78]. Ibis body condition may have declined as the chicks prepared to fledge, a period where many bird species reduce their body mass/condition [79], [80]. It is believed that this reduction in overall mass allows birds to change from plump immobile chicks to slender and proficient flying juveniles [80][82].

Conservation and Management Implications

There were two primary findings of this study that have the potential to benefit conservation of wading bird species. First, we found support for chick physiology to be influenced by hydrology variables such as water depths. Water managers have the potential to regulate the hydrology throughout much of the Everglades through a system of canals, levees, and pumps, theoretically allowing them to adjust water depths for wading birds. This hydrology variable has also been previously linked to pre-breeding adult physiological condition [37], nesting success [24], and foraging site selection of egrets and ibises [23]. Optimal water depths for foraging wading birds differ by species [22], [23], foraging strategies (searchers and exploiters), and years depending on prey biomass conditions; however maintaining water depths between −22 and 21 cm has been suggested to provide suitable foraging conditions for both species and foraging groups [23].

We also found that current levels of Hg exposure in the Everglades may influence the expression of heat shock proteins in egrets. Understanding what the consequences are for these changes in heat shock protein concentrations may provide insight as to how Hg may have impacted Everglades' wading bird populations in the past. This understanding may be relevant to future patterns of wading bird nesting because there are portions of the Everglades that still contain high levels of Hg in prey species of wading birds [10].

Supporting Information

Figure S1.

Location of study area indicating wading bird nesting colonies sampled during 2006 and 2007. The primary Water Management Areas (WCAs), Arthur R. Marshall Loxahatchee National Wildlife Refuge (Lox), and Everglades National Park (ENP) are indicated.

https://doi.org/10.1371/journal.pone.0106447.s001

(TIF)

Table S1.

Data used in analyses of the Physiological condition of juvenile wading birds in relation to multiple landscape stressors in the Florida Everglades: effects of hydrology, prey availability, and mercury bioaccumulation.

https://doi.org/10.1371/journal.pone.0106447.s002

(PDF)

Acknowledgments

Disclaimer: Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government. We appreciate the support and cooperation of the Arthur R. Marshall Loxahatchee National Wildlife Refuge staff. We thank our field research crews and fellow researchers that assisted in collection of field data and laboratory analyses: T. Anderson, N. Baker, T. Beck, B. Botson, K. Donner, H. Herring, N. Hill, A. Horton, B. Imdieke, B. Johnson, S. Lantz, and J. Pierce.

Author Contributions

Conceived and designed the experiments: GH CES DEG JTA. Performed the experiments: GH JEB. Analyzed the data: GH. Contributed reagents/materials/analysis tools: GH DEG CES. Contributed to the writing of the manuscript: GH CES DEG JEB JTA.

References

  1. 1. Wingfield JC (1988) Changes in reproductive function of free-living birds in direct response to environmental perturbations. In: Stetson, J. H. (ed) Processing of environmental information in vertebrates. Berlin; Springer-Verlag: 121–148.
  2. 2. Astheimer LB, Buttemer WA, Wingfield JC (1992) Interactions of corticosterone with feeding, activity and metabolism in passerine birds. Ornis Scandinavia 23: 355–365.
  3. 3. Kitaysky AS, Kitaiskaia EV, Piatt JF, Wingfield JC (2003) Benefits and costs of increased levels of corticosterone in seabird chicks. Hormone Behav 43: 140–149.
  4. 4. Hahn TP (1998) Reproductive seasonality in an opportunistic breeder, the Red Crossbill, Loxia curvirostra. Ecology 79: 2365–2375.
  5. 5. Kacelnik A (1984) Central place foraging in starlings (Sturnus vulgaris) I: Patch residence time. J Anim Ecol 53: 283–299.
  6. 6. Marvin-DiPasquale MC, Agee JL, Bouse RM, Jaffe BE (2003) Microbial cycling of mercury in contaminated pelagic and wetland sediments of San Pablo Bay, California. Environ Geol 43: 260–267.
  7. 7. Hall BD, Aiken GR, Krabbenhoft D, Marvin-DiPasquale M, Swarzenski CM (2004) Wetlands as principal zones of methylmercury production in southern Louisiana and the Gulf of Mexico region. Environ Pollut 154: 124–134.
  8. 8. Sundlof SF, Spalding MG, Wentworth JD, Steible CK (1994) Mercury in livers of wading birds (Ciconiiformes) in southern Florida. Arch Environ Contam Toxicol 27: 299–305.
  9. 9. Frederick PC, Spalding MG, Dusek K (2002) Wading birds as bioindicators of mercury contamination in Florida: annual and geographic variation. Environ Toxicol Chem 21: 163–167.
  10. 10. Rumbold DG, Lange TR, Axelrad DM, Atkeson TD (2008) Ecological risk of methylmercury in Everglades National Park, Florida, USA. Ecotoxicology 17: 632–641.
  11. 11. Adams EM, Frederick PC, Larkin ILV, Guillette LJ (2009) The effects of chronic methylmercury exposure upon fecal testostosterone, estradiol, and corticosterone levels in juvenile White Ibises (Eudocimus albus). Environ Toxicol Chem 28: 982–989.
  12. 12. Jayasena N, Frederick PC, Larkin ILV (2011) Endocrine disruption in white ibises (Eudocimus albus) caused by exposure to environmentally relevant levels of methylmercury. Aquat Toxicol 105: 321–327.
  13. 13. Franceschini MD, Lane OP, Evers DC, Reed JM, Hoskins B, et al. (2009) The corticosterone stress response and mercury contamination in free-living tree swallows, Tachycineta bicolor. Ecotoxicology 18: 514–521.
  14. 14. Wada HD, Cristol A, McNabb FMA, Hopkins WA (2009) Suppressed adenocortical responses and thyroid hormone levels in birds near a mercury-contaminated river. Environ Sci Technol 43: 6031–6038.
  15. 15. Herring G, Ackerman JT, Herzog M (2012) Mercury exposure suppresses baseline stress levels in juvenile birds. Environ Sci Technol 46: 6339–6346.
  16. 16. Jackson AK, Evers DC, Etterson MA, Condon AM, Folsom SB, et al. (2011) Mercury exposure affects the reproductive success of a free-living terrestrial songbird, the Carolina wren (ThryothorusLudovicianus). Auk 128: 759–769.
  17. 17. Kushlan JA (1986) Response of wading birds to seasonally fluctuating water levels: strategies and their limits. Waterbirds 9: 155–162.
  18. 18. Frederick PC, Spalding MG (1994) Factors affecting reproductive success of wading birds (Ciconiiformes) in the Everglades ecosystem. Pages 659–692 in Everglades: The Ecosystem and Its Restoration (Davis SM, Ogden JC, Eds.). St. Lucie Press, Delray Beach, FL.
  19. 19. Ogden JC (1994) A comparison of wading bird nesting colony dynamics (1931–1946 and 1974–1989) as an indication of ecosystem conditions in the southern Everglades. Pages 533–570 in Everglades: The Ecosystem and Its Restoration (Davis SM, Ogden JC, Eds.). St. Lucie Press, Delray Beach, FL.
  20. 20. Heath JA, Frederick PC (2005) Relationships among mercury concentrations, hormones, and nesting effort of White Ibises in the Florida Everglades. Auk 122: 255–267.
  21. 21. Frederick P, Jayasena N (2010) Altered pairing behavior and reproductive success in white ibises exposed to environmentally relevant concentrations of methylmercury. Proc Royal Soc Lond B 278: 1851–1857.
  22. 22. Gawlik DE (2002) The effects of prey availability on the numerical response of wading birds. Ecol Monog 72: 329–346.
  23. 23. Beerens JM, Gawlik DE, Herring G, Cook MI (2011) Dynamic habitat selection by two wading bird species with divergent foraging strategies in a seasonally fluctuating wetland. Auk 128: 651–661.
  24. 24. Herring G, Gawlik DE, Cook MI, Beerens JM (2010) Sensitivity of Great Egret and White Ibis nesting to reduced prey density. Auk 127: 660–670.
  25. 25. Frederick PC, Spalding MG, Sepulveda MS, Williams Jr GE, Nico L, et al. (1998) Exposure of great egret nestlings to mercury through diet in the Everglades of Florida. Environ Toxicol Chem 18: 1940–7.
  26. 26. Rumbold DG, Niemczyk SL, Flink LE, Chandrasekhar T, Harkanson B, et al. (2001) Mercury in eggs and feathers of great egrets (Ardea alba) from the Florida Everglades. Arch Environ Contam Toxicol 41: 501–507.
  27. 27. Herring G, Gawlik DE, Rumbold DA (2009) Feather mercury concentrations and physiological condition of Great Egret and White Ibis nestlings in the Florida Everglades. Science of the Total Environment 407: 2641–2649.
  28. 28. Herring G, Eagles-Smith C, Ackerman J, Gawlik D, Beerens J (2013) Landscape factors and hydrology influence mercury concentrations in wading birds breeding the Florida Everglades, USA. Sci Total Environ 460: 637–646.
  29. 29. Benson J, Suryan RM, Piatt JF (2003) Assessing chick growth from a single colony visit to a seabird colony. Marine Ornithol 31: 181–184.
  30. 30. Wasser SK, Hunt KE, Brown JL, Cooper K, Crockett CM (2000) A generalized fecal glucocorticoid assay for use in a diverse array of nondomestic mammalian and avian species. Gen Comp Endocrinol 120: 260–275.
  31. 31. Herring G, Cook MI, Gawlik DE, Call EM (2011) Food availability is expressed through physiological stress indicators in nestling white ibis: a food supplementation experiment. Func Ecol 25: 682–690.
  32. 32. Sørensen JG, Kristensen TN, Loeschke V (2003) The evolutionary and ecological role of heat shock proteins. Ecol Letters 6: 1025–1037.
  33. 33. Tomás G, Martínez J, Merino S (2004) Collection and analysis of blood samples to detect stress proteins in wild birds. J Field Ornith 75: 281–278.
  34. 34. Herring G, Gawlik DE (2007) The role of stress proteins in the study of allostatic overload in birds, use and applicability to current studies in avian ecology. Sci World J 7: 1596–1602.
  35. 35. Wingfield JC (1994) Modulation of the adrenocortical response to stress in birds. In Davey KG, Peter RE, Tobe SS (eds). Perspectives in Comparative Endocrinology National Research Council of Canada, Ottawa, Ontario, Canada.
  36. 36. Kitaysky AS, Piatt JF, Wingfield JC, Romano M (1999) The adrenocortical stress response of black-legged kittiwake chicks in relation to dietary restrictions. J Comp Phys B 69: 303–310.
  37. 37. Herring G, Gawlik DE (2012) Differential physiological response to landscape prey availability by the great egret and white ibis. J Wild Manage 77: 58–67.
  38. 38. Linquist S (1986) The heat shock response. Ann Rev Biochem 55: 1151–1191.
  39. 39. Bierkens JGEA (2000) Applications and pitfalls of stress-proteins in biomonitoring. Toxicology 153: 61–72.
  40. 40. Locke M, Noble EG (1995) Stress proteins: the exercise response. Can J Appl Phys 20: 155–167.
  41. 41. Råberg L, Grahn M, Hasselquist D, Svensson E (1998) On the adaptive significance of stress-induced immunosupression. Proc Royal Soc Lon B 265: 1637–1641.
  42. 42. Papaconstantinou AD, Brown KM, Noren BT, McAlister T, Fisher BR, et al. (2003) Mercury, cadmium, and arsenite enhance heat shock protein synthesis in chicken embryos prior to embryotoxicity. Birth Defects Res B 68: 456–464.
  43. 43. Liu YR, Zheng YM, Zhang LM, Luan YX, He JZ (2010) Effects of mercury on reproduction, avoidance, and heat shock protein gene expression of the soil springtail Folsomia candida. Environ Toxicol Chem 29: 654–659.
  44. 44. Gaunt AS, Oring LW (1997) Guidelines for the Use of Wild Birds in Research. Washington: Ornithological Council Special Publication.
  45. 45. Herring G, Gawlik DE (2009) Stability of avian fecal corticosterone metabolite levels in frozen avian feces. J Wild Manage 73: 1010–1013.
  46. 46. Khan MZ, Altmann I, Isani SS, Yu J (2002) A matter of time: evaluating the storage of fecal samples for steroid analysis. Gen Comp Endocrinol 128: 57–64.
  47. 47. Myers GJ, Davidson PW (1998) Prenatal methylmercury exposure and children: neurological development, and behavioral research. Environ Health Perspect 106: 841–847.
  48. 48. Day RD, Segars AL, Arendt MD, Lee AM, Peden-Adams MM (2007) Relationship of blood mercury levels to health parameters in the loggerhead sea turtle (Caretta caretta). Environ Health Perspect 115: 1421–1428.
  49. 49. Chen JW, Chang CC, Wang SL (2011) Mercury concentration of whole blood and red blood cell in Taiwan residents. Epidemiology 22: S289.
  50. 50. Lieske CL, Moses SK, Castellini JM, Kleka J, Hueffer K, et al. (2011) Toxicokinetics of mercury in blood compartments and hair of fish-fed sled dogs. Acta Vet Scanda 2011 53: 66.
  51. 51. Thompson DR, Furness RW (1989) Comparison of the levels of total and organic mercury in seabird feathers. Mar Pollut Bull 20: 577–579.
  52. 52. Fournier F, Karasov WH, Kenow KP, Meyer MW, Hines RK (2002) The oral bioavailability and toxicokinetics of methylmercury in common loons (Gavia immer) chicks. Comp Biochem Physiol A 133: 703–714.
  53. 53. U.S. Environmental Protection Agency (2000) Method 7473, Mercury in solids and solutions by thermal decomposition, amalgamation, and atomic absorption spectrophotometry, test methods for evaluating solid waste, physical/chemical methods SW 846, Update IVA. U.S. Government Printing Office, Washington, DC.
  54. 54. U. S. Geological Survey (2006) The Everglades Depth Estimation Network (EDEN) for Support of Ecological and Biological Assessments. U. S. Department of the Interior. Fact Sheet 2006–3087.
  55. 55. Liu Z, Volin JC, Owen VD, Pearlstine LG, Allen JR (2009) Validation and ecosystem applications of the EDEN water-surface model for the Florida Everglades. Ecohydrology 2: 182–194.
  56. 56. JMP (2001) A guide to statistics and data analysis using JMP and JMP IN software. Second Edition SAS Institute, Cary, North Carolina, USA.
  57. 57. Smith JP (1997) Nesting season food habits of 4 species of herons and egrets at Lake Okeechobee, Florida. Waterbirds 20: 198–220.
  58. 58. Loftus WF (2000) Accumulation and fate of mercury in an Everglade aquatic food web; Doctoral Thesis. Florida International University: Miami, FL.
  59. 59. Dorn NJ, Cook MI, Herring G, Boyle R, Nelson J, et al. (2011) Diet variation among White Ibis (Eudocimus albus) chicks: prey composition depends on recent hydrologic conditions. Ibis 153: 323–335.
  60. 60. Neter J, Kutner MH, Nachtsheim CJ, Wasserman W (1996) Applied linear regression models. Third Edition. Times Mirror Higher Education Group, Chicago, Illinois, USA
  61. 61. Nakagawa S, Schielzeth H (2013) A general and simple method for obtaining R2 from generalized linear mixed-effects models. Method Ecol Evol 4: 133–142.
  62. 62. Beauliea M, Thierry AM, Raclot T, Maho YL, Robert-Coudert Y, et al. (2009) Sex-specific paraetal startegies according to the sex of offspring in the Adélie penguin. Behav Ecol 20: 878–883.
  63. 63. Franzellitti S, Fabbri E (2006) Cytoprotective responses in the Mediterranean mussel exposed to Hg2+ and CH3Hg+. Biochem Biophys Res Comm 351: 719–725.
  64. 64. Sørensen JG (2010) Application of heat shock protein expression for detecting natural adaptation and exposure to stress in natural populations. Curr Zool 56: 703–713.
  65. 65. Hoffman DJ, Eagles-Smith CA, Ackerman JT, Adelsbach TL, Stebbins KR (2011) Oxidative stress response of Forster's terns (Sterna forsteri) and Caspian terns (Hydroprogne caspia) to mercury and selenium bioaccumulation in liver, kidney, and brain. Environ Toxicol Chem 30: 920–929.
  66. 66. Brkljačić J, Perišić T, Dundjerski J, Matić G (2007) Interaction of rat renal glucorticoid receptor with Hsp90 and Hsp70 upon stress provoked by mercury. J Appl Toxicol 27: 43–50.
  67. 67. Hwang GW (2012) Role of intracellular defense factors against methylmercury toxicity. Biol Pharm Bull 35: 1881–1884.
  68. 68. Sepulveda MS, Frederick PC, Spalding MG, Williams GE (1998) Mercury contamination in free-ranging great egret nestlings (Ardea albus) from southern Florida. Environ Toxicol Chem 18: 985–992.
  69. 69. Schwabl H (1999) Developmental changes and among-sibling variation of corticosterone levels in an altricial avian species. Gen Comp Endocrinol 116: 403–408.
  70. 70. Kern M, Bacon W, Long D, Cowie RJ (2001) Possible roles for corticosterone and critical size in the fledging of nestling pied flycatchers. Physiol Biochem Zool 74: 651.
  71. 71. Romero LM, Strochlic D, Wingfield JC (2005) Corticosterone inhibits feather growth: potential mechanism explaining seasonal down regulation of corticosterone during molt. Comp Biochem Phys A 142: 65–73.
  72. 72. Desrochers DW, Reed JM, Awerman J, Kluge JA, Wilkison J (2009) Exogenous and endogenous corticosterone alter feather quality. Comp Biochem Phys A 152: 46–52.
  73. 73. Butler MW, Leppert LL, Duffy AL (2010) Effects of small increases in corticosterone levels on morphology, immune function, and feather development. Phys Biogeo Zool 83: 78–86.
  74. 74. Kitaysky AS, Wingfield JC, Piatt JF (1999) Corticosterone facilitates begging and affects resource allocation in the Black-legged Kittiwake. Behav Ecol 12: 619–625.
  75. 75. Loiseau C, Sorci G, Dano S, Chastel O (2008) Effects of experimental increase of corticosterone levels on begging behavior, immunity and parental provisioning rate in house sparrows. Gen Comp Endocrinol 155: 101–108.
  76. 76. Wingfield JC (2005) Modulation of the adrenocortical response to acute stress in breeding birds. In: Dawson A, Sharp PJ (eds.), Functional Avian Endocrinology. Narosa Publishing House, New Delhi, India.
  77. 77. Frederick PC, Collopy MW (1989) Nesting success of Ciconiiform species in relation to water conditions in the Florida Everglades. Auk 106: 625–634.
  78. 78. Gabriel JE, Ferro MIT, Ferro JA, Silva MM, Givisiez PEN, et al. (2000) Influence of dietary energy level on hepatic 70-kDa heat shock protein expression in broiler chicks submitted to acute heat stress. Braz J Poult Scie 2: 259–266.
  79. 79. Gray CM, Hamer KC (2001) Prefledging mass recession in Manz shearwaters: parental desertion or nestling anorexia? Anim Behav 62: 705–709.
  80. 80. Mauck RA, Ricklefs RE (2005) Control of fledging age in Leach's storm petrel Oceanodroma leucorhoa: chick development and pre-fledging mass loss. Func Ecol 19: 73–80.
  81. 81. Martins TLF (1997) Fledging in common swift, Apus apus; weight watching with a difference. Anim Behav 54: 99–108.
  82. 82. Wright J, Markham S, Denney SM (2014) Facultative adjustment of pre-fledging mass loss by nestling swifts preparing for flight. Proc Royal Soc Lon B 273: 1895–1900.