Amphibian breeding phenology influences offspring size and response to a common wetland contaminant
Frontiers in Zoology volume 18, Article number: 31 (2021)
Increases in temperature variability associated with climate change have critical implications for the phenology of wildlife across the globe. For example, warmer winter temperatures can induce forward shifts in breeding phenology across taxa (“false springs”), which can put organisms at risk of freezing conditions during reproduction or vulnerable early life stages. As human activities continue to encroach on natural ecosystems, it is also important to consider how breeding phenology interacts with other anthropogenic stressors (e.g., pollutants). Using 14 populations of a widespread amphibian (wood frog; Rana sylvatica), we compared 1) growth; 2) tolerance to a common wetland contaminant (NaCl); and 3) the ability of tadpoles to acclimate to lethal NaCl exposure following sublethal exposure earlier in life. We evaluated these metrics across two breeding seasons (2018 and 2019) and across populations of tadpoles whose parents differed in breeding phenology (earlier- versus later-breeding cohorts). In both years, the earlier-breeding cohorts completed breeding activity prior to a winter storm and later-breeding cohorts completed breeding activities after a winter storm. The freezing conditions that later-breeding cohorts were exposed to in 2018 were more severe in both magnitude and duration than those in 2019.
In 2018, offspring of the later-breeding cohort were larger but less tolerant of NaCl compared to offspring of the earlier-breeding cohort. The offspring of the earlier-breeding cohort additionally were able to acclimate to a lethal concentration of NaCl following sublethal exposure earlier in life, while the later-breeding cohort became less tolerant of NaCl following acclimation. Interestingly, in 2019, the warmer of the two breeding seasons, we did not detect the negative effects of later breeding phenology on responses to NaCl.
These results suggest that phenological shifts that expose breeding amphibians to freezing conditions can have cascading consequences on offspring mass and ability to tolerate future stressors but likely depends on the severity of the freeze event.
Identifying how species are affected by anthropogenic climate change is an urgent challenge in ecology. One prominent ecological consequence of climate change is phenological shifts in critical life-history events associated with rising average temperatures and increased temperature variability [14, 53]. In temperate zones, warming temperatures have resulted in progressively earlier starts to spring breeding and migratory activities since the 1960’s ; this pattern is strong for species across taxa (e.g., marine animals – ; amphibians – ; plants – ; birds – ; insects – ; mammals – ). One cost of such temperature-induced shifts in the phenology of organisms are the risks associated with “false springs” [6, 37]. As climate change causes warmer winter temperatures, annual reproduction events may be prematurely cued (i.e., a false spring) exposing breeding individuals and/or their offspring at vulnerable early life stages to season-typical frost or freezing temperatures. The risk of false springs was identified early in our understanding of the effects of climate change on species’ phenology . Since then, there have been numerous documented instances of freeze-induced plant death and damage due to premature budding and shooting (e.g., in the winters of 2007  and 2010 ).
Of all the plant and animal taxa examined in a comprehensive meta-analysis, amphibians were identified as the taxon with the most advanced phenological shifts. Amphibian spring activities, on average, have shifted more than twice as early as that of trees, birds, or butterflies, and eight times as early as non-woody plants; the most extreme amphibian species exhibited early shifts of roughly 30 days/decade over the study period . Notably, in a study on eastern North American amphibians, those that bred earliest (i.e., Rana sylvatica, Pseudacris crucifer) were also those shown to have the strongest climate-change related phenological shifts toward even earlier breeding . In light of the increased variability in temperature and associated false spring events, early spring breeding amphibians are likely to fall victim to direct effects from post-emergence freezing. Though temperate adult amphibians have strategies to cope with freezing temperatures , repeated freeze-thaw cycles carry physiological costs  that can potentially lead to carry over effects that negatively affect offspring. Investigating potential carry over effects associated with phenological shifts have important implications for amphibian population persistence in the face of climate change.
Wildlife populations are also faced with threats from exposure to environmental contaminants . Of growing concern is the salinization of freshwater systems, which has increased over the past several decades due to agricultural irrigation, coastal flooding and the application of road salts [12, 29]. Elevated salinities can result in direct mortality of freshwater species leading to deleterious outcomes for wildlife populations . For wood frogs, salinization can negatively influence population demography , rates of growth and development via altered activity levels , and can reduce survival . Unlike some contaminants, salts are conserved within the aquatic environment  and can accumulate over time. This issue is further compounded by the predicted effects of climate warming on wetland evaporation, where increased evaporation is predicted to reduce water volume, potentially further elevating salt concentrations in freshwater systems . These concurrent issues are particularly problematic for spring-breeding amphibians, who breed following spring rains that can flush road salts applied during the winter months into breeding habitats , exposing amphibians to increased salinities across multiple life stages [2, 38]. Given the potential for increased freshwater salinization under climate change, there is a need to understand how altered temperature regimes and salts interactively effect wildlife health .
While salinization can negatively impact wildlife health, previous research suggests that wildlife populations are able to acclimate to contaminant exposure [31, 32]. For example, Wu et al.  found that pre-acclimation of larval Indian rice frogs (Fejervarya limnocharis) in saltwater concentrations of 7 g L− 1 for 48 h increased their survival to an otherwise lethal dose of saltwater at a concentration of 11 g L− 1. In addition to contaminant acclimation, research suggests that varying temperatures, such as those expected under climate change  also play a role in the ability of organisms to withstand contaminant exposure [43, 46]. For example, Kimberly and Salice  exposed Physa polmin snails to different temperature regimes (22 or 28 °C) and then measured their sensitivity to cadmium. No differences in cadmium toxicity were found when snails were acclimated to either temperature alone. However, when snails initially acclimated to one temperature were switched to the other (i.e., 22 to 28 °C, or 28 to 22 °C) at the time of cadmium exposure, snail mortality greatly increased . Thus, contaminant acclimation and temperature variability can both individually alter the toxicity of contaminants to wildlife. However, as temperature variability under climate change scenarios can additionally influence the breeding phenology of organisms, there is a need to not only understand the influence of temperature on the acclimation of contaminant tolerance, but also how breeding phenology may interact with temperature to influence the ability of organisms to acclimate to contaminant exposure. Indeed, by potentially exposing adults to stressful freezing conditions prior to breeding, shifts in breeding phenology may lead to carry over effects that influence the ability for their offspring to tolerate or acclimate to contaminants.
Using a widespread early spring-breeding amphibian (wood frog; Rana sylvatica) as our model, the goal of this study was to evaluate how adult exposure to different temperature conditions due to phenological variation in breeding times influences offspring growth and performance (baseline tolerance to NaCl and ability to acclimate to NaCl). Specifically, across two breeding seasons, we compared the growth and performance of offspring from earlier-breeding parental cohorts vs. later-breeding parental cohorts. In both years, later-breeding parental cohorts experience an additional freeze-thaw cycle prior to breeding compared to earlier-breeding parental cohorts. Because repeated freeze-thaw cycles carry physiological costs , we predicted that tadpoles of parents that delay breeding in our study populations and encountered a winter storm would experience carryover effects making them smaller and less able to tolerate and less able to acclimate to NaCl exposure.
Methods and materials
To evaluate how differences in adult breeding phenology influences offspring growth and responses to NaCl, we compared wood frog tadpoles from populations that bred prior to freezing conditions (earlier breeding cohorts) vs. tadpoles from populations that bred after freezing conditions (later breeding cohorts) across two breeding seasons (2018 and 2019). In 2018, we collected egg masses from seven populations of wood frogs that bred on 6 April (earlier cohort) and egg masses from seven different populations that bred after a winter storm on 13 April (late cohort; Fig. 1; Fig. 2; Supplementary Materials Table 1). Similarly, in 2019, we collected egg masses from five populations of wood frogs that bred on 30 March (earlier cohort) and egg masses from each of 5 populations that bred after a winter storm on 7 April (late cohort; Fig. 1; Fig. 2; Supplementary Materials Table 1). All egg masses were collected from populations in western Pennsylvania, USA (Fig. 1; Supplementary Materials Table 1). The breeding window of wood frog populations in this region is characterized by breeding beginning in early March, continuing into April . Thus, the seven-day difference in breeding times between the “earlier” and “late” populations represents a relatively significant difference in time of breeding given the short breeding window of wood frog populations in this region. Further, in addition to the seven-day difference in breeding dates, earlier- and- late-breeders are distinct from one another as only the later-breeding cohorts bred immediately after a winter storm. Thus, in both years, later-breeding parental cohorts experienced an additional freeze-thaw cycle prior to breeding compared to earlier-breeding parental cohorts.
To ensure broad sampling of each population, in both 2018 and 2019, we collected 10 partial egg masses from each site. To minimize any effect of environment on the egg masses, we collected egg masses within 48 h of oviposition.
Natal pond conductivity
To account for variation that may arise in tadpole susceptibility to NaCl due to differences in natal pond exposure to NaCl, we measured pond conductivity at the time of egg collection. Water samples were taken from the edge and middle of each wetland, and specific conductance was measured (mean ± SE) using a YSI Multiparameter Sonde (Xylem, OH, USA). Specific conductance is a measure of water’s ability to pass electrical flow and is directly related to total ion concentration within the water . Thus, our conductivity measurements do not distinguish between different salt ions within our water samples (e.g., chloride, sodium, magnesium, calcium, etc.). However, specific conductance has been shown to be strongly correlated with concentrations of sodium and chloride ions within bodies of water , particularly when road salts, such as NaCl, are thought to be a contributing factor to conductance . Road salt applications within the counties that our populations reside in are applied largely with pure NaCl or salt brine (a mixture of NaCl and water ;). Thus, our measurements of specific conductance are a reliable predictor of salinization resulting from inputs of NaCl across our pond sites.
To assess mass, baseline tolerance to NaCl, and ability to acclimate to a lethal concentration of NaCl, we used a two-phase experimental design following Hua et al. . We standardized our methodology across years and between earlier and late cohorts. In Phase 1, we employed a fully factorial design placing 50 embryos (Gosner stage 8; Gosner 1960) from each of the respective populations into 0 (background Cl− concentration of 0.034 g L− 1), 0.50, or 1 g L− 1 NaCl solutions in 2 L plastic tubs, filled with 1.5 L of the respective Phase 1 pre-treatment solutions (Fig. 3). Once tadpoles reached the larval stage (Gosner stage 25), we measured the mass of a subset of tadpoles (see Supplementary Materials Table 2 for sample sizes) from all populations across the three treatment concentrations. We then initiated Phase 2, the time-to-death assay (hereafter TTD). Across both years and cohorts, all animals used in the mass measurements and TTD assay reached Gosner stage 25 within 12 h of one another. To start the TTD assay, we removed individuals from their respective Phase 1 exposure treatment and moved them to either a control (no salts added; N = 10) or a lethal NaCl treatment of 8 g L− 1 NaCl (N = 15). Experimental units were 100 mL cups filled with 80 mL of either control or NaCl treatment water, with a single tadpole/ experimental unit and were randomized across a single shelving unit. All experiments were carried out under identical temperature-controlled conditions (20 °C) between cohorts and years.
TTD assays are designed using concentrations of contaminants that cause mortality, while still allowing for the detection of variation in timing of mortality across treatments and populations . Thus, we chose 8 g L− 1 NaCl, as it has been shown to be lethal to wood frog tadpoles [15, 56]. Additionally, 8 g L− 1 (~ 4.85 g L− 1 Cl−), albeit on the higher end of salt concentrations seen in natural systems, is environmentally relevant, with urban ponds, creeks, and rivers reaching concentrations of up to 13.50 g Cl− L− 1 . For information on the creation of the Phase 1 and Phase 2 NaCl solutions, please see Supplementary Materials. We measured the TTD of all individuals by assessing tadpole mortality every 2 hours until 100% mortality. In accordance with standard toxicology procedure, we did not feed any individuals during the TTD assay . We euthanized all tadpoles that survived (i.e., control animals) using an overdose of unbuffered 5% MS-222 and then preserved all tadpoles used in the experiment in a 10% formalin solution.
Natal pond conductivity
We used a linear mixed effect model to test for fixed effects of year, breeding cohort, and their interaction on pond conductivity. Following significant main effects, we conducted adjusted Bonferroni tests to evaluate pairwise comparisons.
Tadpole mass at Gosner 25
We used a linear mixed effect model to evaluate the mass of tadpoles from earlier vs. late cohorts after being reared in the various NaCl-acclimation treatments (0 g L− 1, 0.5 g L− 1, 1 g L− 1) in 2018 and 2019. We measured all tadpoles at Gosner stage 25 to control for differences in mass associated with development. Our model included cohort, acclimation concentration, year, and their interactions as fixed effects and pond conductivity as a random effect. We included pond conductivity in our model because the wood frog natal ponds varied in their conductivity (Fig. 4). For all significant main effects and interactions, we conducted Bonferroni-pairwise comparisons. Due to lost samples, we do not have mass data for tadpoles from the earlier-breeding cohort in 2018. These analyses were performed using SPSS version 26.
We used Cox Mixed Effect models to compare the survival rates of tadpoles from earlier vs. late cohorts in 2018 and 2019. For baseline tolerance, we only considered tadpoles that were reared in the 0 g L− 1 NaCl acclimation treatment. In the first analysis, we included cohort, year, and their interaction as fixed effects and population and pond conductivity as random effects. Because there was a significant cohort*year interaction, we conducted two additional analyses, one for 2018 and one for 2019. These analyses were performed in R version 3.6. 2.
Effect of acclimation on mortality rate
We used Cox regression analysis to calculate beta coefficient (β), a measure of change in mortality risk for tadpoles exposed to the acclimation treatments relative to a one-unit change in mortality rate of tadpoles exposed to the control. For both 2018 and 2019, we calculated a separate β for each population by comparing the survival curve of tadpoles exposed to each acclimation treatment (0.5 g L− 1 and 1 g L− 1) to survival curves of tadpoles from the same population that were not exposed to NaCl (0 g L− 1). A negative β -value indicates that tadpoles exposed to the acclimation treatment had lower risk of mortality compared to tadpoles not exposed to the acclimation treatment. A positive β-value indicates that tadpoles exposed to the acclimation treatment had higher risk of mortality compared to tadpoles not exposed to the acclimation treatment. After calculating β as our metric for acclimation, we conducted a linear mixed effect model to compare acclimation in tadpoles (β) from earlier vs. late cohorts exposed to the various NaCl-acclimation treatments (0.5 g L− 1 vs. 1 g L− 1) in 2018 and 2019. Our model included cohort, NaCl-acclimation concentration, year, and their interactions as fixed effects and pond conductivity as a random effect. For all significant main effects and interactions, we conducted Bonferroni-pairwise comparisons. These analyses were performed using SPSS version 26.
Natal pond conductivity
Using a linear mixed effect model, we found no effect of year (F1, 7.8 = 1.7; p = 0.23) or breeding cohort*year interaction (F1, 8.5 = 0.001; p = 0.98) on natal pond conductivity. However, there was a significant main effect of breeding cohort (F1, 12.5 = 12.5; p = 0.013) on natal pond conductivity. The conductivity of ponds for populations of tadpoles from later breeding cohorts was higher than the conductivity of ponds for populations of tadpoles from earlier breeding cohorts (Fig. 4).
To better understand whether patterns of natal pond conductivity influenced tadpole tolerance to NaCl, we conducted correlation analyses. We found no relationship between natal pond conductivity and tadpole tolerance to NaCl for tadpoles reared in 0 g L− 1 (r = − 0.22; p = 0.31), 0.5 g L− 1 (r = 0.34; p = 0.12), or 1 g L− 1 (r = − 0.18; p = 0.41). These data suggest that while pond conductivity differs between earlier and later cohorts, this variation does not appear to influence tadpole tolerance to NaCl. However, as a conservative approach, we include natal pond conductivity as a random factor in all analyses.
Tadpole mass at Gosner 25
Using a linear mixed effect model, we found a significant main effect of cohort (F1,21.3 = 5.0; p = 0.04), NaCl-acclimation treatment (F2,599.9 = 14.9; p < 0.001), and year (F1,11.5 = 27.6; p < 0.001) on tadpole mass at Gosner stage 25. Additionally, while there was no significant interaction between Cohort*year (F1,21.8 = 1.3; p = 0.25), there were significant interactions between Cohort*NaCl-acclimation (F2,600.1 = 10.2; p < 0.001), Year* NaCl-acclimation (F2,600.5 = 12.6; p < 0.001), and Year*NaCl-acclimation*Cohort (F1,600.5 = 6.6; p = 0.01). To facilitate the description of the significant interactions and in accordance with our a priori hypotheses, we separately describe the effect of earlier vs. late cohorts on tadpole mass at each year for each of the three NaCl-acclimation treatments. For tadpoles exposed to the 0 g L− 1 rearing condition, we do not have data for 2018 but for 2019, there was no significant difference in mass between earlier vs. late cohorts (p = 0.44; Fig. 5). For tadpoles exposed to the 0.5 g L− 1 rearing condition, in 2018, there was no significant difference in mass between earlier vs. late cohorts (p = 0.44; Fig. 5). However, in 2019, there was a significant difference in mass between earlier vs. late cohorts (p = 0.04; Fig. 5). Tadpoles from earlier breeding cohorts were smaller than tadpoles in later breeding cohorts. Finally, for tadpoles exposed to the 1 g L− 1 rearing condition, there was a significant difference in mass between earlier vs. late cohorts in both 2018 and 2019 (p = 0.004 and p = 0.009, respectively; Fig. 5).
In the overall Cox Mixed Effect model, while there was no effect of year (β = 0.02; p = 0.9), we found a significant effect of cohort (β = 5919.3; p < 0.001) and cohort*year (β = − 2.9; p < 0.001). To understand the interaction, we additionally conducted two separate Cox Mixed Effects analyses for 2018 vs. 2019. In 2018, we found a significant effect of cohort (β = 4.4; p < 0.001). Tadpoles from earlier breeding cohorts were significantly more tolerant to lethal concentrations of NaCl than tadpoles from later breeding cohorts. However, in 2019, we did not find a significant effect of cohort (β = 0.94; p = 0.7; Fig. 6).
Effect of acclimation on mortality rate
Using a linear mixed effect model, we found no significant main effect of cohort (F1,24 = 1.2; p = 0.28), NaCl-acclimation treatment (F1,24 = 1.2; p = 0.93), or year (F1,24 = 0.9; p = 0.35) on tadpole acclimation (β) to NaCl (for details on β (Beta coefficient) calculations and figures, see Supplementary Materials Figures a2, a3, a4, a5, a6, a7, a8 and a9). There was also no significant interaction between Cohort*NaCl-acclimation (F1,24 = 1.2; p = 0.93), Year* NaCl-acclimation (F = 0.09; p = 0.76), or Year*NaCl-acclimation*Cohort (F1,24 = 1.1; p = 0.31). In contrast, we did find a significant interaction between Cohort*year (F1,24 = 4.8; p = 0.038). To better understand this interaction, we conducted Bonferroni-corrected pairwise comparisons. In 2018, we found a significant effect of cohort (F1,24 = 4.4; p < 0.001) on the ability for tadpoles from earlier vs. later-breeding cohorts to acclimate to higher concentrations of NaCl. Tadpoles from earlier-breeding cohorts reared in either 0.5 g L− 1 or 1 g L− 1 NaCl had lower β -values compared to tadpoles from later breeding cohorts (Fig. 7). This indicates that tadpoles from earlier breeding cohorts that were reared in 0.5 g L− 1 or 1 g L− 1 became more tolerant to higher concentrations of NaCl (negative β) whereas tadpoles from later-breeding cohorts became less tolerant to higher concentrations of NaCl (positive β).
Anthropogenic activities have serious consequences for species’ survival, particularly those species directly affected by both climate change and environmental contaminants. Here, we evaluated how breeding phenology (earlier versus later-breeding cohorts) and a common environmental pollutant interact to affect a widespread spring-breeding amphibian. At a single study region over 2 years, several populations of R. sylvatica began breeding activities early in response to spring-like temperature and precipitation cues; other populations did not breed until later, after the year’s final severe winter weather events were completed. In both years, when reared in the highest NaCl acclimation treatment, the earlier-breeding cohort produced smaller offspring than the later-breeding cohort. Interestingly, the effect of breeding phenology on tadpole mass seemed to be exacerbated as NaCl levels increased in the rearing environment. Although we do not have data for the earlier cohort of 2018, comparisons between the earlier and late cohorts of 2019 show no difference in their mass. However, when reared at an intermediate concentration of NaCl (0.5 g L− 1), the earlier cohort of 2019 produced significantly smaller individuals than that of the later-breeding cohort of the same year. Finally, when reared in the highest NaCl acclimation treatment (1 g L− 1 NaCl), the offspring of the later-breeding cohorts of 2018 and 2019 weighed ~ 37% and ~ 40% more on average, respectively, than offspring of the earlier-breeding cohorts within their respective years. Consistent with other pollutant studies, this suggests that the interactive effect of phenology and NaCl on mass is likely dose dependent .
We additionally found that wood frog tadpoles from the later-breeding cohort of 2018 were less tolerant of NaCl than tadpoles from the earlier-breeding cohort of the same year. While direct exposure to cold temperatures has been documented to lower immune responses [16, 54], reduce larval growth and development , and interact synergistically with other stressors , it is important to note that in this study, tadpoles were not directly exposed to the cold temperatures of the parental environment. Thus, the reduced baseline tolerance to NaCl of the late breeding cohort may reflect cross-generational consequences of parental exposure to suboptimal temperatures . For example, Piiroinen et al. , found that exposure of Colorado potato beetles (Leptinotarsa decemlineata) to cold stress resulted in deleterious cross-generational effects, where offspring of cold-stressed parents were less tolerant of pyrethroid insecticides than those whose parents were reared at more optimal temperatures. Under climate change, we expect wildlife to contend with increasingly variable temperatures that can influence breeding phenology as well as facing challenges associated with pollutant exposure . As such, our findings underscore the importance of evaluating the consequences of phenological shifts in concert with other stressors, such as NaCl exposure.
To better cope with contaminants in their environment, some individuals can increase their future tolerance to pollutants following an acclimation period to sublethal concentrations of the pollutant. Consistent with past work , we demonstrate that some populations of wood frogs are also capable of inducing increased tolerance to NaCl following exposure to NaCl early in life as embryos and hatchlings. Though previous research demonstrates that wood frogs can induce increased tolerance to pesticides , to our knowledge, this is the first evidence of NaCl-acclimation in this species. Whether this phenomenon is broadly generalizable across toxicants is still unknown. Though there is some evidence that wood frog populations that can induce tolerance to one pollutant are capable of inducing cross-tolerance to pollutants that share similar modes of action . Interestingly, the ability for wood frogs to acclimate to NaCl differed between earlier versus later-breeding cohorts in 2018. Specifically, the later-breeding cohort of 2018 whose parents experienced a freezing event prior to breeding gave rise to offspring that were less capable of acclimating to NaCl exposure compared to the earlier-breeding cohort of the same year whose parents did not experience such events. We found that this pattern held regardless of the sublethal NaCl concentration at which tadpoles were exposed to prior to lethal exposure (0.5 or 1 g L− 1 NaCl). While the mechanisms underlying the reduced ability of the later-breeding cohort to acclimate to NaCl exposure were beyond the scope of this study, several studies suggest that amphibian acclimation to saline conditions is related to shifts in stress hormones, such as corticosterone (CORT) . Extreme weather events, such as freezing temperatures, can upregulate CORT production in wildlife [10, 66], which can then be passed on to offspring . While the expression of stress hormones contributes to addressing short term challenges, chronically elevated CORT is often associated with negative fitness effects . Thus, future work might evaluate whether the reduced ability of the later-breeding cohort of 2018 to acclimate to NaCl exposure could be due to cross-generational fitness costs of parental expression of stress hormones. Collectively, while there is strong evidence demonstrating that amphibians can acclimate to diverse pollutants, this study suggests that considering parental environmental conditions may be critical to understanding the relative importance of acclimation as a mechanism for population persistence in the face of global change.
Interestingly, we found that the influence of breeding phenology (earlier vs. late) on baseline tolerance and acclimation differed between years. Temperature averages and variability throughout the two winters of our study can help to explain the differences in NaCl baseline tolerance and acclimation in 2018 and 2019. The later-breeding R. sylvatica cohort in 2018 experienced a more severe winter weather event than the later-breeding cohort in 2019 (5 versus 2 days of freezing temperatures in 2018 and 2019, respectively). The later-breeding cohort in 2018 therefore likely experienced greater physiological stress from the freeze than the later-breeding cohort in 2019; there are costly physiological responses to low temperatures [19, 63] and maternal stress is often transferred to developing eggs [18, 62] in a variety of herpetofauna. The higher level of stress experienced by the later-breeding cohort in 2018 may be one potential cause of reduced offspring baseline tolerance to NaCl compared to that in 2019. The difference in offspring NaCl tolerance (following NaCl acclimation) from earlier and later-breeding cohorts likewise reflects the pattern shown in the baseline tolerance: the offspring of later breeders were more susceptible than those of earlier breeders in 2018, whereas the 2019 earlier- and later-breeders’ offspring did not differ in susceptibility. Our findings suggest that small shifts in temperature in the parental environment may impact the viability of offspring in salinized systems. As shifting temperature regimes are likely to continue under global climate change and subsequently increase NaCl pollution across some regions of the globe [5, 41], our findings suggest a need for cross-generational studies to identify whether the patterns between earlier and late cohort survival following NaCl exposure seen here were indeed a product of freezing temperatures in the parental environment.
In this system, breeding timing (earlier vs. late) does not appear to be a locally adapted trait. Instead, it is well established that breeding behavior in wood frogs is induced by environmental cues . More specifically, in the system we tested, in some breeding seasons (i.e., years with high variability in environmental cues), wood frog breeding times across populations are variable often splitting into earlier vs. later breeding populations [32, 34]. In contrast, in other years, these same wood frog populations may all breed on the same day . The differential breeding times amongst and within wood frog populations across years used in our study are in line with those of other species, where differential breeding phenology has been posited to be due to differences in microclimate and/or weather [40, 44]. Importantly, as seen in this study, even within a population, there is evidence of annual variation in breeding time. Indeed, for the three populations (MIN, SKN, and MAL; Fig. 1), wood frogs differed in whether they laid earlier or later between 2018 and 2019. Because of this switch in relative breeding phenology between years, we were able to compare NaCl tolerance in early vs. later-breeding cohorts within a particular population. Consistent with the overall conclusions across populations, we found that tadpoles of earlier-breeders were more tolerant of NaCl than tadpoles of late-breeders even while holding population constant (SI- Figure A10). Collectively, this suggests that the early-late phenotypic differences in NaCl tolerance seen between the 2 years of study are unlikely to be completely due to underlying population-specific traits.
Lastly, we found that the conductivity of natal pond sites utilized by the wood frog populations used in our experiments differed. For some populations of amphibians, natal pond conductivity is associated with tolerance to NaCl exposure [8, 9, 30]. However, our correlation analyses revealed that natal pond conductivity was not a causal factor for the differential tolerance to NaCl seen in our study populations. Further, as eggs were removed from natal ponds ≤48 h following oviposition and were then placed immediately into filtered well water, the direct effects of natal pond environment on wood frog embryos were likely limited. Given the lack of correlation between pond conductivity and NaCl tolerance, along with the limited exposure of embryos to the natal environment, this altogether suggests the main contributor to the responses seen in our study were due to differences in breeding phenology and/or temperature in the parental environment.
Our experimental design was highly conservative. By removing eggs from vernal ponds following breeding and rearing them in a temperature-controlled common garden environment, developing embryos from the early-breeding cohorts were protected from extreme low temperatures during the freeze. Thus, although offspring from the early-breeding cohort may have enhanced tolerance to NaCl compared to the later-breeding cohort, we did not account for the fact that these embryos likely would have experienced reduced growth and increased mortality in freezing ponds. In addition, the average temperature in the winter of 2019 was 1.03 °C warmer than that of 2018, and 2019 populations experienced an earlier phenological shift (i.e., both earlier and late breeders began breeding ~ 1 week earlier in 2019 than their counterparts in 2018). Warmer winter weather and earlier breeding (as documented at our sites in 2019) results in colder developmental temperatures for R. sylvatica larvae and corresponding reduced rates of development . Overall, had offspring been reared in a field setting, it is likely that a) earlier cohorts in 2018 and 2019 would have experienced greater freeze-related mortality, and b) all larvae in 2019 would have experienced slowed development and exposure to stressful (cold) conditions. The consequences of these freezing conditions are evaluated in a companion paper . Taken together, our data suggests that freezes following false springs (exacerbated by warmer average winter temperatures), particularly when combined with chemical stressors, are likely to be a major concern for spring breeding amphibian populations. Despite some debate on climate change’s direct role in amphibian extinctions , amphibians are increasingly identified as being highly vulnerable to climate change . We suggest that the interactive effect of climate-induced phenological shifts and environmental contaminants is an underexplored and direct threat to amphibian populations.
Availability of data and materials
Data associated with this manuscript can be found at Binghamton University’s Open Repository (https://orb.binghamton.edu/bio_students/4/).
Acevedo-Whitehouse K, Duffus ALJ. Effects of environmental change on wildlife health. PhilosTrans R Soc Lond B Biol Sci. 2009;364(1534):3429–38. https://doi.org/10.1098/rstb.2009.0128.
Albecker MA, McCoy MW. Adaptive responses to salinity stress across multiple life stages in anuran amphibians. Front Zool. 2017;14(1):40. https://doi.org/10.1186/s12983-017-0222-0.
Allstadt AJ, Vavrus SJ, Heglund PJ, Pidgeon AM, Thogmartin WE, Radeloff VC. Spring plant phenology and false springs in the conterminous US during the 21st century. Environ Res Lett. 2015;10(10):104008. https://doi.org/10.1088/1748-9326/10/10/104008.
Andrews KM, Nanjappa P, Riley SPD. Roads and ecological infrastructure: concepts and applications for small animals. Baltimore: John's Hopkins University Press; 2015.
Arvidsson A, Blomqvist G, Öberg G. Impact of climate change on use of anti-icing and deicing salt in Sweden. In: Winter maintenance and surface transportation weather: international conference on winter maintenance and surface transportation weather, vol. 30. Coralville: IA; 2012. p. 3–10.
Ault TR, Henebry GM, De Beurs KM, Schwartz MD, Betancourt JL, Moore D. The false spring of 2012, earliest in north American record. Eos. 2013;94(20):181–2. https://doi.org/10.1002/2013EO200001.
Benard MF. Warmer winters reduce frog fecundity and shift breeding phenology, which consequently alters larval development and metamorphic timing. Glob Chang Biol. 2015;21(3):1058–65. https://doi.org/10.1111/gcb.12720.
Brady SP. Road to evolution? Local adaptation to road adjacency in an amphibian (Ambystoma maculatum). Sci Rep. 2012;2(1):235. https://doi.org/10.1038/srep00235.
Brady SP. Microgeographic maladaptive performance and deme depression in response to roads and runoff. PeerJ. 2013;1:e163. https://doi.org/10.7717/peerj.163.
Breuner CW, Wingfield JC, Romero LM. Diel rhythms of basal and stress-induced corticosterone in a wild, seasonal vertebrate, Gambel’s white-crowned sparrow. J Exp Zool. 1999;284(3):334–42. https://doi.org/10.1002/(SICI)1097-010X(19990801)284:3<334::AID-JEZ11>3.0.CO;2-#.
Budischak SA, Belden LK, Hopkins WA. Effects of malathion on embryonic development and latent susceptibility to trematode parasites in ranid tadpoles. Environ Toxicol Chem. 2008;27(12):2496–500. https://doi.org/10.1897/08-018.1.
Cañedo-Argüelles M, Hawkins CP, Kefford BJ, Schäfer RB, Dyack BJ, Brucet S, et al. Saving freshwater from salts. Science. 2016;351(6276):914–6. https://doi.org/10.1126/science.aad3488.
Cannell MGR, Smith RI. Climatic warming, spring budburst, and forest damage on trees. J Appl Ecol. 1986;23(1):177–91. https://doi.org/10.2307/2403090.
Cohen JM, Lajeunesse MJ, Rohr JR. A global synthesis of animal phenological responses to climate change. Nat Clim Chang. 2018;8(3):224–8. https://doi.org/10.1038/s41558-018-0067-3.
Collins SJ, Russell RW. Toxicity of road salt to Nova Scotia amphibians. Environ Pollut. 2009;157(1):320–4. https://doi.org/10.1016/j.envpol.2008.06.032.
Cone RE, Marchalonis JJ. Cellular and humoral aspects of the influence of environmental temperature on the immune response of poikilothermic vertebrates. J Immunol. 1972;108(4):952–7.
Corn PS. Climate change and amphibians. Anim Biodivers Conserv. 2005;28(1):59–67.
Cree A, Tyrrell CL, Preest MR, Thorburn D, Guillette LJ. Protecting embryos from stress: corticosterone effects and the corticosterone response to capture and confinement during pregnancy in a live-bearing lizard (Hoplodactylus maculatus). Gen Comp Endocrinol. 2003;134(3):316–29. https://doi.org/10.1016/S0016-6480(03)00282-X.
Dupoué A, Brischoux F, Lourdais O, Angelier F. Influence of temperature on the corticosterone stress-response: an experiment in the children’s python (antaresia childreni). Gen Comp Endocrinol. 2013;193(2013):178–84. https://doi.org/10.1016/j.ygcen.2013.08.004.
Evans M, Frick C. The effects of road salts on aquatic ecosystems; 2001.
Findlay SEG, Kelly VR. Emerging indirect and long-term road salt effects on ecosystems. Ann N Y Acad Sci. 2011;1223(1):58–68. https://doi.org/10.1111/j.1749-6632.2010.05942.x.
Flemming MH. PennDOT LTAP technical information sheet #167. Pennsylvania: Department of Transportation; 2015.
Foden WB, Butchart SHM, Stuart SN, Vié JC, Akçakaya HR, Angulo A, et al. Identifying the World’s Most climate change vulnerable species: a systematic trait-based assessment of all birds, amphibians and corals. PLoS One. 2013;8(6):e65427. https://doi.org/10.1371/journal.pone.0065427.
Gibbs JP, Breisch AR. Climate warming and calling phenology of frogs near Ithaca, New York, 1900-1999. Conserv Biol. 2001;15(4):1175–8. https://doi.org/10.1046/j.1523-1739.2001.0150041175.x.
Gray BS, Lethaby M. The amphibians and reptiles of Erie County, Pennsylvania. Bull Maryland Herpetol Soc. 2008;44(2):49–69.
Gu L, Hanson PJ, Post WM, Kaiser DP, Yang B, Nemani R, et al. The 2007 eastern US spring freeze: increased cold damage in a warming world? BioScience. 2008;58(3):253–62. https://doi.org/10.1641/B580311.
Hall EM, Brady SP, Mattheus NM, Earley RL, Diamond M, Crespi EJ. Physiological consequences of exposure to salinized roadside ponds on wood frog larvae and adults. Biol Conserv (Elsevier). 2017;209:98–106.
Hau M, Casagrande S, Ouyang JQ, Baugh AT. Glucocorticoid-mediated phenotypes in vertebrates: multilevel variation and evolution. Adv Study Behav. 2016;48:41–115. https://doi.org/10.1016/bs.asb.2016.01.002.
Hintz WD, Relyea RA. A review of the species, community, and ecosystem impacts of road salt salinisation in fresh waters. Freshw Biol. 2019;64(6):1081–97. https://doi.org/10.1111/fwb.13286.
Hopkins GR, Brodie ED Jr. Occurrence of amphibians in saline habitats: a review and evolutionary perspective. Herpetol Monogr. 2015;29(1):1–27. https://doi.org/10.1655/HERPMONOGRAPHS-D-14-00006.
Hsu W-T, Wu C-S, Lai J-C, Chiao Y-K, Hsu C-H, Kam Y-C. Salinity acclimation affects survival and metamorphosis of crab-eating frog tadpoles. Herpetologica. 2012;68(1):14–21. https://doi.org/10.1655/HERPETOLOGICA-D-11-00018.1.
Hua J, Jones DK, Mattes BM, Cothran RD, Relyea RA, Hoverman JT. Evolved pesticide tolerance in amphibians: predicting mechanisms based on pesticide novelty and mode of action. Environ Pollut. 2015;206:56–63. https://doi.org/10.1016/j.envpol.2015.06.030.
Hua J, Jones DK, Relyea RA. Induced tolerance from a sublethal insecticide leads to cross-tolerance to other insecticides. Environ Sci Technol (ACS Publications). 2014;48:4078–85.
Hua J, Morehouse NI, Relyea R. Pesticide tolerance in amphibians: induced tolerance in susceptible populations, constitutive tolerance in tolerant populations. Evol Appl. 2013;6(7):1028–40. https://doi.org/10.1111/eva.12083.
Hua J, Wuerthner VP, Jones DK, Mattes B, Cothran RD, Relyea RA, et al. Evolved pesticide tolerance influences susceptibility to parasites in amphibians. Evol Appl. 2017;10(8):802–12. https://doi.org/10.1111/eva.12500.
Hufkens K, Friedl MA, Keenan TF, Sonnentag O, Bailey A, O’Keefe J, et al. Ecological impacts of a widespread frost event following early spring leaf-out. Glob Chang Biol. 2012;18(7):2365–77. https://doi.org/10.1111/j.1365-2486.2012.02712.x.
Inouye DW. The ecological and evolutionary significance of frost in the context of climate change. Ecol Lett. 2000;3(5):457–63. https://doi.org/10.1046/j.1461-0248.2000.00165.x.
Karraker NE, Gibbs JP. Road deicing salt irreversibly disrupts osmoregulation of salamander egg clutches. Environ Pollut. 2011;159(3):833–5. https://doi.org/10.1016/j.envpol.2010.11.019.
Karraker NE, Gibbs JP, Vonesh JR. Impacts of road deicing salt on the demography of vernal Pool-breeding amphibians. Ecol Appl. 2008;18(3):724–34. https://doi.org/10.1890/07-1644.1.
Karraker NE, Gibbs JP. Amphibian production in forested landscapes in relation to wetland hydroperiod: a case study of vernal pools and beaver ponds. Biol Conserv (Elsevier). 2009;142:2293–302.
Kaushal SS, Groffman PM, Likens GE, Belt KT, Stack WP, Kelly VR, et al. Increased salinization of fresh water in the northeastern United States. Proc Natl Acad Sci. 2005;102(38):13517–20. https://doi.org/10.1073/pnas.0506414102.
Kemker C. Conductivity, Salinity and Total Dissolved Solids. In: Fundamentals of Environmental Measurements. Fairborn: Fondriest Environmental; 2014.
Kimberly DA, Salice CJ. Complex interactions between climate change and toxicants: evidence that temperature variability increases sensitivity to cadmium. Ecotoxicology. 2014;23(5):809–17. https://doi.org/10.1007/s10646-014-1221-y.
Klaus SP. Correlates and temporal variation in call phenology of Eastern Ontario frogs. Queen’s University (Canada); 2012.
Masson-Delmotte TWV, Zhai P, Pörtner HO, Roberts D, Skea J, Shukla PR, et al. Summary for policymakers. In: Global warming of 1.5 C. an IPCC special report on the impacts of global warming of 1.5 C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global. Geneva: IPCC; 2018.
Møller V, Forbes VE, Depledg MH. Influence of acclimation and exposure temperature on the acute toxicity of cadmium to the freshwater snail Potamopyrgus antipodarum (Hydrobiidae). Environ Toxicol Chem. 1994;13(9):1519–24. https://doi.org/10.1002/etc.5620130914.
Newman MC. Fundamentals of ecotoxicology. 3rd ed. Boca Raton: Taylor & Francis; 2010.
Nixon CM, McClain MW. Squirrel population decline following a late spring frost. J Wildl Manag. 1969;33(2):353–7. https://doi.org/10.2307/3799835.
Norberg-King T. Toxicity reduction and toxicity identification evaluations for effluents, ambient waters, and other aqueous media: Society of Environmental Toxicology and Chemistry; 2005.
Parmesan C. Influences of species, latitudes and methodologies on estimates of phenological response to global warming. Glob Chang Biol. 2007;13(9):1860–72. https://doi.org/10.1111/j.1365-2486.2007.01404.x.
Piiroinen S, Lyytinen A, Lindström L. Stress for invasion success? Temperature stress of preceding generations modifies the response to insecticide stress in an invasive pest insect. Evol Appl. 2013;6(2):313–23. https://doi.org/10.1111/eva.12001.
Poloczanska ES, Brown CJ, Sydeman WJ, Kiessling W, Schoeman DS, Moore PJ, et al. Global imprint of climate change on marine life. Nat Clim Chang. 2013;3(10):919–25. https://doi.org/10.1038/nclimate1958.
Renner SS, Zohner CM. Climate change and Phenological mismatch in trophic interactions among plants, insects, and vertebrates. Annu Rev Ecol Evol Syst. 2018;49(1):165–82. https://doi.org/10.1146/annurev-ecolsys-110617-062535.
Rollins-Smith LA, Woodhams DC. Amphibian immunity. New York: Oxford University Press; 2012.
Romero LM. Physiological stress in ecology: lessons from biomedical research. Trends Ecol Evol. 2004;19(5):249–55. https://doi.org/10.1016/j.tree.2004.03.008.
Sanzo D, Hecnar SJ. Effects of road de-icing salt (NaCl) on larval wood frogs (Rana sylvatica). Environ Pollut. 2006;140(2):247–56. https://doi.org/10.1016/j.envpol.2005.07.013.
Searle CL, Belden LK, Bancroft BA, Han BA, Biga LM, Blaustein AR. Experimental examination of the effects of ultraviolet-B radiation in combination with other stressors on frog larvae. Oecologia. 2010;162(1):237–45. https://doi.org/10.1007/s00442-009-1440-8.
Turtle SL. Embryonic survivorship of the spotted salamander (Ambystoma maculatum) in roadside and woodland vernal pools in southeastern New Hampshire. J Herpetol. 2000;34:60–7.
Uchiyama M, Konno N. Hormonal regulation of ion and water transport in anuran amphibians. Gen Comp Endocrinol. 2006;147(1):54–61. https://doi.org/10.1016/j.ygcen.2005.12.018.
Vitt LJ, Caldwell JP. Herpetology. 3rd ed. Amsterdam: Elsevier; 2009.
Walther G-R, Post E, Convey P, Menzel A, Parmesan C, Beebee TJC, et al. Ecological responses to recent climate change. Nature. 2002;416(6879):389–95. https://doi.org/10.1038/416389a.
Warne RW, Kardon A, Crespi EJ. Physiological, behavioral and maternal factors that contribute to size variation in larval amphibian populations. PLoS One. 2013;8(10):e76364. https://doi.org/10.1371/journal.pone.0076364.
Wersebe M, Blackwood P, Guo YT, Jaeger J, May D, Meindl G, et al. The effects of different cold-temperature regimes on development, growth, and susceptibility to an abiotic and biotic stressor. Ecol Evol. 2019;9(6):3355–66. https://doi.org/10.1002/ece3.4957.
Williams CM, Henry HAL, Sinclair BJ. Cold truths: how winter drives responses of terrestrial organisms to climate change. Biol Rev. 2015;90(1):214–35. https://doi.org/10.1111/brv.12105.
Wilson RJ, Maclean IM. Recent evidence for the climate change threat to Lepidoptera and other insects. J Insect Conserv. 2011;15(1–2):259–68. https://doi.org/10.1007/s10841-010-9342-y.
Wingfield JC, Moore MC, Farner DS. Endocrine responses to inclement weather in naturally breeding populations of white-crowned sparrows (Zonotrichia leucophrys pugetensis). Auk. 1983;100(1):56–62. https://doi.org/10.1093/auk/100.1.56.
Wu C-S, Yang W-K, Lee T-H, Gomez-Mestre I, Kam Y-C. Salinity acclimation enhances salinity tolerance in tadpoles living in brackish water through increased Na+, K+-ATPase expression. J Exp Zool A Ecol Genet Physiol. 2014;321(1):57–64. https://doi.org/10.1002/jez.1837.
We thank Devin DiGiacopo, Grascen Shidemantle, Vanessa Wuerthner, Logan Billet and Jared Jaeger for their assistance in collecting egg masses used in the experiments. Thank you to Jessica Barabas and Christopher Davis at Pymatuning Lab of Ecology for help with housing and wood frog monitoring. We additionally thank Devin DiGiacopo, Grascen Shidemantle and Brianna Sander for their help conducting the time-to-death experiments and Sydney Tredo and Rachel Verdi for helping to collect tadpole mass data. Finally, we would like to thank the anonymous reviewers whose comments helped to improve our manuscript.
Funding for this work was provided by a National Science Foundation grant #1655190 to JH.
Ethics approval and consent to participate
This study was approved by Binghamton University’s international animal care and use committee (Protocol 812–18).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Creation of salt solutions. Collection, analysis and results of temperature data. Table 1. Coordinates for wood frog populations. Table 2. Sample size of wood frog tadpoles that were weighed within each respective breeding cohort (early or late), year (2018 or 2019) and NaCl acclimation treatment (0, 0.5 or 1 g L-1 NaCl) combination. Figure A1. Natal pond conductivity for early vs. late breeding cohorts in 2018 and 2019. Different letters represent significance at p <0.05. Figure A2. Survival curves of tadpoles from seven populations of wood frogs in the early breeding cohort in 2018. The blue line represents tadpole tolerance to 8 g L-1 when not placed in an acclimation treatment. The red line represents tadpole tolerance to 8 g L-1 when reared in 0.5 g L-1 NaCl from the egg to tadpole stage. Figure A3. Survival curves of tadpoles from seven populations of wood frogs in the late breeding cohort in 2018. The blue line represents tadpole tolerance to 8 g L-1 when not placed in an acclimation treatment. The red line represents tadpole tolerance to 8 g L-1 when reared in 0.5 g L-1 NaCl from the egg to tadpole stage. Figure A4. Survival curves of tadpoles from seven populations of wood frogs in the early breeding cohort in 2018. The blue line represents tadpole tolerance to 8 g L-1 when not placed in an acclimation treatment. The red line represents tadpole tolerance to 8 g L-1 when reared in 1 g L-1 NaCl from the egg to tadpole stage. Figure A5. Survival curves of tadpoles from seven populations of wood frogs in the late breeding cohort in 2018. The blue line represents tadpole tolerance to 8 g L-1 when not placed in an acclimation treatment. The red line represents tadpole tolerance to 8 g L-1 when reared in 1 g L-1 NaCl from the egg to tadpole stage. Figure A6. Survival curves of tadpoles from five populations of wood frogs in the early breeding cohort in 2019. The blue line represents tadpole tolerance to 8 g L-1 when not placed in an acclimation treatment. The red line represents tadpole tolerance to 8 g L-1 when reared in 0.5 g L-1 NaCl from the egg to tadpole stage. Figure A7. Survival curves of tadpoles from five populations of wood frogs in the late breeding cohort in 2019. The blue line represents tadpole tolerance to 8 g L-1 when not placed in an acclimation treatment. The red line represents tadpole tolerance to 8 g L-1 when reared in 0.5 g L-1 NaCl from the egg to tadpole stage. Figure A8. Survival curves of tadpoles from five populations of wood frogs in the early breeding cohort in 2019. The blue line represents tadpole tolerance to 8 g L-1 when not placed in an acclimation treatment. The red line represents tadpole tolerance to 8 g L-1 when reared in 1 g L-1 NaCl from the egg to tadpole stage. Figure A9. Survival curves of tadpoles from five populations of wood frogs in the late breeding cohort in 2019. The blue line represents tadpole tolerance to 8 g L-1 when not placed in an acclimation treatment. The red line represents tadpole tolerance to 8 g L-1 when reared in 1 g L-1 NaCl from the egg to tadpole stage. Figure A10. Survival curves of tadpoles from three populations of wood frogs. The blue line represents the tolerance of tadpoles in the early-breeding cohort to NaCl. The red line represents the tolerance of tadpoles in the late-breeding cohort to NaCl.
About this article
Cite this article
Buss, N., Swierk, L. & Hua, J. Amphibian breeding phenology influences offspring size and response to a common wetland contaminant. Front Zool 18, 31 (2021). https://doi.org/10.1186/s12983-021-00413-0
- Anthropogenic climate change
- Contaminant tolerance
- Road salt