Age- and environment-dependent changes in chemical defences of larval and post-metamorphic toads
BMC Evolutionary Biology volume 17, Article number: 137 (2017)
Chemical defences are widespread in animals, but how their production is adjusted to ecological conditions is poorly known. Optimal defence theory predicts that inducible defences are favoured over constitutive defences when toxin production is costly and the need for it varies across environments. However, if some environmental changes occur predictably (e.g. coupled to transitions during ontogeny), whereas others are unpredictable (e.g. predation, food availability), changes in defences may have constitutive as well as plastic elements. To investigate this phenomenon, we raised common toad (Bufo bufo) tadpoles with ad libitum or limited food and in the presence or absence of chemical cues on predation risk, and measured their toxin content on 5 occasions during early ontogeny.
The number of compounds showed limited variation with age in tadpoles and was unaffected by food limitation and predator cues. The total amount of bufadienolides first increased and later decreased during development, and it was elevated in young and mid-aged tadpoles with limited food availability compared to their ad libitum fed conspecifics, whereas it did not change in response to cues on predation risk. We provide the first evidence for the active synthesis of defensive toxin compounds this early during ontogeny in amphibians. Furthermore, the observation of increased quantities of bufadienolides in food-restricted tadpoles is the first experimental demonstration of resource-dependent induction of elevated de novo toxin production, suggesting a role for bufadienolides in allelopathy.
Our study shows that the evolution of phenotypic plasticity in chemical defences may depend on the ecological context (i.e. predation vs. competition). Our results furthermore suggest that the age-dependent changes in the diversity of toxin compounds in developing toads may be fixed (i.e., constitutive), timed for the developmental stages in which they are most reliant on their chemical arsenal, whereas inducible plasticity may prevail in the amount of synthesized compounds.
Chemical defences are widespread across the animal kingdom [1, 2] and can serve for deterring predators, parasites, competitors, and pathogens [1,2,3,4,5,6]. Some species sequester toxic compounds from food or symbionts [4, 6, 7], or obtain them from ambiguous sources [5, 8, 9], while others are capable of de novo synthesizing toxins [3, 4, 10]. However, in species that synthesise toxic compounds themselves, it has remained largely unknown if chemical defences are inducible, i.e. if their production can vary plastically in response to changing environmental conditions  and how inducible chemical defences change during ontogeny .
Plastic responses are known to evolve under variable environmental conditions and to come with inherent costs [13, 14]. Therefore, induced chemical defences are especially likely to occur in animals that encounter unpredictably varying environments during their lifetime and synthesise toxins de novo, since such synthesis relies on a specialized biochemical pathway and associated physiological and anatomical structures and therefore is considered to be costly [15, 16]. On the other hand, optimal defence theory predicts that changes in chemical defences may become constitutive when environmental differences are predictable ; for example, if individuals predictably encounter new environments during their life as a consequence of their development. Consequently, in animals that undergo large, predictable changes in their life-history, and thereby become exposed to drastically different environments that also can unpredictably vary in ecologically important characteristics, chemical defences may both have constitutive as well as inducible components, similarly to other types of defences (e.g. [14, 17, 18]).
Among vertebrates, amphibians undergo the most dramatic changes during their post-embryonic development when they metamorphose and leave the aquatic environment to embark on a terrestrial life [19, 20]. Therefore, amphibians are ideal for studies on ontogenetic changes in toxin production and on the inducibility vs. constitutive nature of chemical defences. Also, chemical defences of vertebrates have been most extensively studied in amphibians. While toxin composition of many amphibian species is well known [3, 7, 21], and experiments documenting age-dependent changes in susceptibility to predators are prevalent in the literature (e.g. [22,23,24,25,26]), in-depth studies on ontogenetic changes in the quantity and composition of toxins utilized in chemical defence and on the underlying secretory apparatus are relatively rare [12, 27,28,29,30,31,32]. Moreover, there are only a handful of studies on phenotypic plasticity in chemical defences in amphibians [12, 33,34,35], and only in two of these were larvae sampled for toxin content [12, 35]. Direct evidence for inducible chemical defences in larvae is lacking, and the ability of tadpoles to synthesize toxic compounds has not been confirmed . Also, the studies that so far reported plastic changes in toxin composition in amphibians, and in fact in any vertebrate [12, 33], only documented effects of predators experienced in the larval environment on post-metamorphic animals, while the metamorphic transition from the fully aquatic larval stage to the terrestrial form disrupts selective forces acting during the two life-stages and makes these largely independent of each other. Therefore, evidence for adaptive phenotypic plasticity in chemical defences in species synthesising toxins de novo is lacking.
Here we present a study on the ontogenetic changes and environmental dependence of toxin content in early life-stages of the common toad (Bufo bufo). We aimed to (1) examine ontogenetic variation in chemical defences in larval and post-metamorphic common toads and (2) investigate if ontogenetic changes in toxin production may be constitutive or induced by environmental conditions that may affect the pay-off of chemical defence. We experimentally manipulated the presence of chemical cues on predation risk (i.e. the need for toxin production) and food availability (i.e. the costliness of toxin production) and repeatedly assessed the toxin content of individuals during early ontogeny. We predicted that if cues on predation risk are present during tadpole development and tadpoles are able to synthesize toxins themselves, they would start producing such compounds earlier on during their ontogeny and in higher quantities compared to predator-naïve conspecifics. Given that de novo toxin synthesis is considered to be costly [15, 16], we also predicted that food restriction would lead to decreased production of defensive chemicals, manifesting in lowered quantities and decreased numbers of compounds compared to well-fed conspecifics. We chose the common toad as the study species, because it displays relatively weak inducible defences during the larval stage in terms of morphology and behaviour [36,37,38] and appears to be unpalatable to several predator species [39, 40], suggesting heavy reliance on chemical defence. Also, the chemical composition of Bufo skin secretions is relatively well known, their main defensive chemicals being bufadienolides and biogenic amines [41,42,43,44], and B. bufo are known to contain toxins already in the larval stages [41, 45].
In early spring 2013 we collected 10 common toad pairs from a lake in the Pilis Mountains, Hungary (47° 37′ 24.78″ N, 18° 48′ 27.20″ E) and transported them to the experimental station of the Plant Protection Institute (Centre for Agricultural Research, Hungarian Academy of Sciences) in Budapest. We let the pairs spawn in 200-L containers placed outdoors, filled with 60 L of aged tap water and containing twigs as egg-deposition substrates. After egg-laying, we transferred eggs from each clutch to the laboratory, and placed them into dishpans filled with reconstituted soft water (RSW; ) to a depth of 2 cm. Temperature was set to 17 °C at the beginning and was allowed to gradually increase to 22 °C by the end of the experiment. We set the lighting to a 13: 11 h light: dark cycle.
Upon hatching, we haphazardly selected four hatchlings of each family and stored them in 70% HPLC-grade methanol, resulting in 40 samples collected at the start of the experiment. Hatchlings were at developmental stage 19 (, Additional file 1: Figure S1). We used this sampling to estimate the baseline of bufadienolide content at the start of larval development. We further assigned randomly selected hatchlings in groups of three to 2-L containers filled with 1.5 L RSW, distributed randomly among treatments. We employed a three-factorial design with two predator-cue treatments (control vs. chemical cues on predation risk), two food level treatments (ad libitum vs. limited food), and four sampling occasions during the larval and early metamorphic life-stages (for details see below). We replicated each combination of predator treatment × food level treatment × sampling occasion 20 times, resulting in a total of 320 experimental units at the start of the experiment. Each family was represented twice in each treatment combination. Containers were arranged in groups of four in a randomized block design, where each block contained tadpoles from one family.
As predators we used five 4th instar larvae of the southern hawker, Aeshna cyanea, and five adult, male smooth newts, Lissotriton vulgaris. We kept individuals of both species grouped in 5-L containers filled with 3 L of RSW, and fed them daily with 800 mg B. bufo and 800 mg agile frog, Rana dalmatina, tadpoles each. We prepared stimulus water by mixing the water taken from the tanks in which we housed and fed the predators, and simulated predation risk by transferring 30 ml stimulus water daily to the rearing containers of tadpoles in the predator-cue treatment group, while adding equal amounts of RSW to the containers of the control group.
We fed tadpoles with a finely ground 4:1 mixture of rabbit chow and fish flakes. Tadpoles assigned to the ad libitum food treatment group received a food amount of ca. 12% of their body mass/individual/day; tadpoles in the limited food treatment group received one-third of that amount. We adjusted food quantity by weighing tadpoles to the nearest mg at the sampling occasions (see below, Additional file 1: Figure S2). We changed the water in the tadpoles’ rearing containers every third day. Whenever we observed a dead individual, we removed it, but spontaneous mortality remained relatively low during the experiment (110 tadpoles out of 960, 11.46%).
When tadpoles were approaching metamorphosis, we monitored the rearing containers twice a day. When the first tadpole in a container started to metamorphose (appearance of at least one forelimb, developmental stage 42 according to ), we removed the other individuals from that container, decreased the water level to 1 dl and slightly tilted the containers to allow the metamorph to leave the water.
After the sampling of hatchlings at the start of the experiment, we took samples for toxin analysis four more times, preserving 40 individuals at each occasion . The second and third samplings took place after 14 and 21 days, when tadpoles reached the median developmental stages of 28 (range = 28–30) and 34 (31–36), respectively. We took a fourth sample when tadpoles reached a median developmental stage of 38 (37–41). The date of this sampling occasion was not specified a priori, but was rather determined based on how developed tadpoles were (the presence of well-formed hind limbs), to account for potential differences in growth rate between treatment groups. We performed a final sampling when individuals completed metamorphosis (complete disappearance of the tail at stage 46, Additional file 1: Figure S1). Each container was sampled once during the entire experiment, by haphazardly selecting and conserving one individual from it. From the 320 experimental containers we therefore collected 320 samples, half of which we later analysed for toxin content. In each treatment × family combination, one container was a priori designated to chemical analysis while the other container was used as a backup; the latter samples were analysed only if we encountered problems during sample preparation for HPLC of the respective a priori sample (21 instances out of 160 samples, 13.13%, ). We released adults, unused eggs and all remaining tadpoles, metamorphs and toadlets at the site of collection.
Analysis of toxin content
We used high-performance liquid chromatography with diode-array detection and mass spectrometry (HPLC-DAD-MS) to identify and quantify bufadienolide compounds. We homogenized specimens using a homogenizer (VWR VDI 12) with a dispersing tool (IKA S12 N-7S). After drying samples in vacuo at 45 °C using a rotary evaporator (Büchi Rotavapor R-134), we measured dry weight of samples using an analytical balance (Ohaus Pioneer PA-114) to the nearest 0.1 mg and subsequently re-dissolved samples in 1 ml HPLC-grade absolute methanol, which was further aided by exposing the samples briefly to ultrasound in a bath sonicator (Tesla UC005AJ1). We filtered the samples using nylon syringe filters (FilterBio, pore size = 0.22 μm). We identified compounds as bufadienolides by inspecting the UV spectrum of peaks [27, 33, 45] and by using commercially acquired bufalin, bufotalin, resibufogenin, gamabufotalin, areno- and telocinobufagin (Biopurify Phytochemicals, Chengdu, China), cinobufagin (Chembest, Shanghai, China), cinobufotalin (Quality Phytochemicals, New Jersey, USA) and digitoxigenin (Santa Cruz Biotechnology, Dallas, TX, USA) as standards (Fig. 1). Identification of compounds present in low quantities was further aided by the analysis of a sample obtained from an adult male common toad by gently massaging the parotoid glands.
A single-quadrupole HPLC-MS system (Shimadzu LC-MS-2020) equipped with a binary gradient solvent pump, a vacuum degasser, a thermostated autosampler, a column oven, a diode array detector and a mass analyser with electrospray ionization (ESI-MS) was used. Chromatographic separations were carried out at 35 °C on a C18 2.6 μm column (Kinetex, 100 mm × 3 mm i.d.) in series with a C18 guard column (4 mm × 3 mm i.d.) using 10 μL injections. The mobile phase consisted of water containing 0.05% formic acid (solvent A) and acetonitrile containing 0.05% formic acid (solvent B). The flow rate was 0.8 mL/min and the gradient was as follows: 0–2 min, 15–25% B; 2–15 min, 25–35% B; 15–24 min, 30–50% B; 24–25 min, 50–90% B; 25–30 min 90% B; 30–35 min 15% B. ESI worked under the following conditions: desolvation line (DL) temperature, 250 °C; heat block temperature, 400 °C; drying N2 gas flow, 15 L/min; nebulizer N2 gas flow, 1.5 L/min; positive ionization mode. Data was acquired and processed using the programme LabSolutions 5.42v (Shimadzu).
To calculate the number of bufadienolide compounds (NBC) present in each animal, we assumed a compound to be present when its area value was larger than zero in the chromatogram (Fig. 1). We estimated the quantity of each compound from the area values of chromatogram peaks (Fig. 1) based on the calibration curve of the bufotalin standard, and summed up these values to obtain an estimate of total bufadienolide quantity (TBQ) for each individual. We used the calibration curve of the bufotalin standard, because this was the most ubiquitously identified compound in our samples (Table 1). This approach yields approximate estimates of bufadienolide quantities, but it has been successfully used before in similar studies [12, 33, 45].
We analysed the effects of predator-cue, food treatments and developmental stage on toxin content using linear mixed-effects models (LMM). We entered NBC or TBQ as the dependent variable; we used the log10-transformed values of TBQ to ensure normality of model residuals and homogeneity of variances. Initial models included food level, predator-cue treatment, and developmental stage of tadpoles as fixed factors, their two-way and three-way interactions; and block nested within family as random factors. In the analyses of TBQ, we also entered the log10-transformed dry mass of tadpoles as a covariate, but without interactions with the other explanatory variables. Note that the data obtained from the first sampling occasion (developmental stage 19) were not included in the LMM analyses, because treatments were only applied after this stage. With each initial model, we performed a backward model-simplification procedure based on P-values, with α = 0.05. To calculate relevant statistics for non-significant terms that were dropped during model selection, we re-entered the removed variables one by one into the final models. We ran all analyses in R 3.1.3  using the ‘lme’ function in the ‘nlme’ package . P-values were calculated using ‘anova’ in ‘nlme’, using type-3 sums of squares. We conducted pairwise comparisons among treatment groups and samplings by calculating linear contrasts corrected for false discovery rate  using the ‘lsmeans’ package . Model residuals of NBC and TBQ showed considerable heteroscedasticity between samplings when developmental stage 19 was included in the analysis (Figs. 2 and 3), therefore in these instances we allowed for different within-sampling variances using ‘weights’ with ‘varIdent’ in ‘nlme’ . We also analysed the quantity of each bufadienolide compound separately; the final models of these analyses are presented in the supplementary information (Additional file 1: Table S1). We had to discard one sample from the analysis on NBC, and two samples from the analysis on TBQ due to missing data .
Number of bufadienolide compounds
Toad hatchlings (developmental stage 19) contained a small number of bufadienolides or none at all (median: 0, range: 0-4, N = 40; 25 hatchlings (62.5%) did not contain any compounds in detectable quantities). In contrast, bufadienolides had high diversity in all other age categories (developmental stages 28–46, median NBC: 14, range: 8–20, N = 159; Table 1, Fig. 2). One compound (unidentified bufadienolide 1) was found in all tadpoles and post-metamorphic individuals (Table 1). After hatching, the effect of developmental stage on NBC was marginally non-significant (Table 2), as tadpoles in any stage did not differ from each other significantly, while there was a small but significant difference between post-metamorphic toads and metamorphosing individuals such that the post-metamorphs had slightly fewer (ca. 1 compound less) bufadienolides (Fig. 2, Additional file 1: Table S2). Predation risk and food limitation did not have a significant effect on NBC (Table 2, Additional file 1: Figure S3).
Total bufadienolide quantity
Toad hatchlings contained only minute amounts of bufadienolides (mean ± SE: 12.603 ± 4.065 ng / tadpole, N = 40) compared to all other age categories (mean ± SE: 1555.864 ± 97.796 ng / tadpole, N = 158, Fig. 3). The compound that was present in all individuals after the hatchling stage had the highest mean quantity (unidentified bufadienolide 1; mean ± SE: 313.738 ± 30.424 ng/tadpole). Across the tadpole stages, the total quantity of bufadienolides increased significantly to developmental stage 34 but decreased afterwards (Table 2, Fig. 3, Additional file 1: Table S2). Furthermore, during the first half of tadpole development (developmental stages 28 and 34), tadpoles that received reduced amounts of food contained significantly more bufadienolides than their ad libitum fed conspecifics (Table 2, Fig. 3, Additional file 1: Table S3), despite having significantly lower body mass (LMM of body mass, age: F3, 214 = 160.694, P = <0.0001; food level: F1, 214 = 26.831, P = <0.0001; age × food level: F3, 214 = 34.735, P = <0.0001; N = 303; Additional file 1: Figure S2), whereas this difference in TBQ disappeared in more developed tadpoles and post-metamorphic individuals (Fig. 3, Additional file 1: Table S3). Presence or absence of chemical cues on predation risk did not influence TBQ (Table 2, Additional file 1: Figure S4). Analysing the quantity of each bufadienolide compound separately corroborated our findings that toxin content varied with age and food level but not with predation-cue treatment (Additional file 1: Table S1).
Our study is the first to unequivocally demonstrate de novo production of toxic compounds in amphibian larvae, as indicated by the steep increase in both the number and quantity of bufadienolide compounds after hatching. This finding shows that common toad tadpoles synthesize their toxins de novo, as has been suggested by histological and ultrastructural studies that demonstrated the presence of the underlying secretory cells and glands already during larval life [30, 31]. This contrasts with other toad species [12, 27] in which tadpoles were found not to produce bufadienolides, relying instead on maternal provisioning of these toxins. For example, in the cane toad (Rhinella marina)  the diversity and amount of bufadienolides are highest in eggs and gradually decrease until developmental stage 25 . Because we were primarily interested in the phenotypic plasticity of toxin production, we did not investigate eggs, so it remains possible that the same decrease from the egg stage to hatching occurs in common toads. This is supported by observations that common toad eggs are repulsive to many different predator species [39, 40], although compounds other than bufadienolides (e.g. biogenic amines) may also be responsible for the unpalatability of common toad eggs (and hatchlings), as suggested for larvae of A. boreas . Nonetheless, because the majority of hatchlings in our study contained no bufadienolides at all, the importance of maternal provisioning of these toxins appears to be limited in common toads. Clearly, maternal toxin provisioning and the temporal changes in toxin content during embryonic development of common toads demands further investigation.
We found that bufadienolides accumulated quickly in young tadpoles and, after reaching a peak in mid-aged larvae, decreased to lower quantities as metamorphosis was approaching. This pattern mirrors ontogenetic changes in tadpole vulnerability: young tadpoles are more vulnerable to predators, thus early toxin production may be strongly favoured. Later, when tadpoles grow larger, they reach a size refuge against several predators or are more difficult to capture [23,24,25,26], therefore they may have to rely less on chemical defences. Such an adjustment of toxin dosage to vulnerability to predators may be common in chemically defended organisms [52,53,54]. Additionally, the ability of predators to distinguish and learn to avoid noxious prey [55,56,57] may relieve tadpoles from synthesising large amounts of toxins in later stages. Despite a decrease in total bufadienolide quantity, metamorphosing and post-metamorphic individuals in our study still contained considerable amounts of bufadienolides, most likely providing them with effective defences against certain predators and perhaps also pathogens and parasites [58, 59], although metamorphosing anurans (between developmental stages 42 and 46) are more susceptible to predation than late tadpole stages or already metamorphosed animals [60,61,62,63,64]. The decreased bufadienolide quantity we observed in these later stages may be attributed to proximate constraints associated with metamorphosis, when a complete re-organization of several physiological systems occurs [19, 65, 66]. We would expect a similar pattern also in other animal species that actively synthesise toxins and undergo substantial morphological and physiological changes during their ontogeny, such as many insects .
The observation that predation risk did not induce the production of larger quantities of bufadienolides in tadpoles is surprising, although it agrees with our earlier finding that, in natural ponds, the toxin content of toad tadpoles did not correlate with the density of predators . A previous experimental study did not observe inducible changes in chemical defences of tadpoles either, but this was attributable to the lack of toxin synthesis in tadpoles of the study species . One possible explanation in our case is that predator-induced changes in chemical defence exist in common toad tadpoles, but not in response to the specific predators we used. However, previous studies did document plastic changes in life-history traits, behaviour and morphology of common toad tadpoles to chemical cues on the presence of A. cyanea [38, 68], suggesting that they can detect their presence based on olfactory cues, perceive these predators as dangerous, and respond to them by changes in life-history traits. Because fishes are considered to be the most voracious predators of amphibian larvae , it is possible that in certain populations that live in permanent ponds, such as the one used in the current study, a relatively high baseline level of bufadienolide synthesis becomes fixed via selection . Alternatively, predator-induced plasticity of toxin production may be lacking in toad tadpoles in general, perhaps because of the toxins’ apparent low production cost , and because very high spatiotemporal variability of predator communities may favour constitutive defences . Finally, although bufadienolides may be effective in repelling several predators, it is possible that the evolution of plasticity in toxin production is driven by other factors, such as pathogens [72, 73] or competitors (, see below).
Our results demonstrated inducible changes in toxin production in response to food availability: during early larval life, food-deprived tadpoles contained significantly more bufadienolides than their ad libitum-fed conspecifics. This result, combined with the fact that toxin content was not related to body mass, corroborates our earlier finding that the energetic costs of toxin production in toad larvae may be low . It seems contradictory that an inducible defence may be cheap to produce, but detecting associated costs of expressed plastic traits may be problematic in species with a complex life-history, such as anurans, because costs may not appear synchronously with the displayed trait [74,75,76,77]. Nonetheless, enhanced toxin production in food-limited tadpoles concords with results of our field study showing that common toad tadpoles in ponds with high density of competitors (mainly amphibian larvae) contained more bufadienolide compounds and slightly larger total quantities of bufadienolides than tadpoles coexisting with fewer competitors . Thus, in the current experiment, reduced food level might have acted as an indicator of high competitor density, inducing the synthesis of larger amounts of bufadienolides against competitors or the pathogens and parasites they carry. Allelopathy, which is intra- or interspecific competition mediated by chemical substances , is a phenomenon of fundamental importance in algae and plants [79, 80], but for animals it has been rarely reported so far [81,82,83]. The existence of chemical interference between amphibian larvae was proposed long ago, but the mediating agents involved in the process have not been identified [69, 84]. Bufadienolides have been suggested to act as allelochemicals , but it remains to be tested directly whether the synthesis of these compounds benefits toad tadpoles by negatively affecting competitors or naturally occurring pathogens and parasites. Nonetheless, our results suggest that allelopathy may be a significant factor in the ecology of a wider variety of animals than currently thought.
In conclusion, our results are the first to document plastic changes in chemical defences in response to food availability in any vertebrate capable of de novo toxin synthesis. The observation that tadpoles produced more toxins at low food availability than when food was present ad libitum indicated that bufadienolides may be relatively cheap to produce but their production may respond plastically to the perceived intensity of competition for food. Our results furthermore suggest that the ontogenetic timing of the production of various toxin components may be fixed as a constitutive defence in toad tadpoles, whereas inducible plasticity prevails in how much of these components is produced. These results, coupled with those of previous studies, highlight the existence of surprisingly diverse strategies of toxin provisioning and synthesis even among as closely related taxa as the species of the Bufonidae family and, thus, caution against premature generalization of observed strategies among species of other chemically defended groups of organisms. Our findings also suggest that ontogenetic changes in toxin production may have resulted from adaptation to predictable variation in predation risk over development, and, thus, represent constitutive age-dependent changes in anti-predator defence rather than a phenotypically plastic response. Therefore, the same trait can show different degrees of phenotypic plasticity depending on evolutionary history (i.e. different species) and ecological context (e.g. predators or other enemies, such as competitors). Studies scrutinizing the costs of toxin production, clarifying the role of toxins in competitive interactions and immune defence, and identifying the environmental factors promoting fixation of the rate of toxin synthesis appear to be especially promising avenues of future research and will provide important insights into the evolution and ecology of chemical defences.
Brodie ED III. Toxins and venoms. Curr Biol. 2009;19(20):R931–5.
Mebs D. Toxicity in animals. Trends in evolution? Toxicon. 2001;39(1):87–96.
Toledo RC, Jared C: Cutaneous granular glands and amphibian venoms. Comp Biochem Physiol A Physiol 1995, 111(1):1-29.
Blum MS. Chemical defenses of arthropods. New York: Academic Press Inc.; 1981.
Tachibana K: Chemical defense in fishes. In: Bioorg Mar Chem. Edited by Scheuer PJ, vol. 2. Berlin Heidelberg, Germany: Springer Verlag; 1988: 117-138.
Pawlik JR. Marine invertebrate chemical defenses. Chem Rev. 1993;93(5):1911–22.
Savitzky AH, Mori A, Hutchinson DA, Saporito RA, Burghardt GM, Lillywhite HB, Meinwald J. Sequestered defensive toxins in tetrapod vertebrates: principles, patterns, and prospects for future studies. Chemoecology. 2012;22(3):141–58.
Hanifin CT, Brodie ED III, Brodie ED Jr. Tetrodotoxin levels of the rough-skinned newt, Taricha granulosa, increase in long-term captivity. Toxicon. 2002;40:1142–53.
Williams BL, Hanifin CT, Brodie ED, Caldwell RL. Ontogeny of tetrodotoxin levels in blue-ringed octopuses: maternal investment and apparent independent production in offspring of Hapalochlaena lunulata. J Chem Ecol. 2011;37(1):10–7.
Stankowich T. Armed and dangerous: predicting the presence and function of defensive weaponry in mammals. Adapt Behav. 2012;20(1):32–43.
Hettyey A, Tóth Z, Van Buskirk J. Inducible chemical defences in animals. Oikos. 2014;123:1025–8.
Benard MF, Fordyce JA. Are induced defenses costly? Consequences of predator-induced defenses in western toads, Bufo boreas. Ecology. 2003;84(1):68–78.
DeWitt TJ, Scheiner SM: Phenotypic variation from single genotypes: a primer. In: Phenotypic plasticity: functional and conceptual approaches. Edited by DeWitt TJ, Scheiner SM. New York, USA: Oxford University Press; 2004: 1-9.
McCall AC, Fordyce JA. Can optimal defence theory be used to predict the distribution of plant chemical defences? J Ecol. 2010;98(5):985–92.
Morgenstern D, King GF. The venom optimization hypothesis revisited. Toxicon. 2013;63:120–8.
Longson CG, Joss JMP. Optimal toxicity in animals: predicting the optimal level of chemical defences. Funct Ecol. 2006;20(4):731–5.
Hamilton R, Siva-Jothy M, Boots M. Two arms are better than one: parasite variation leads to combined inducible and constitutive innate immune responses. Proc R Soc B Biol Sci. 2008;275(1637):937–45.
Bourdeau PE. Intraspecific trait cospecialization of constitutive and inducible morphological defences in a marine snail from habitats with different predation risk. J Anim Ecol. 2012;81(4):849–58.
Brown DD, Cai L. Amphibian metamorphosis. Dev Biol. 2007;306(1):20–33.
Gosner KL. A simplified table for staging anuran embryos and larvae with notes on identification. Herpetologica. 1960;16(3):183–90.
Saporito RA, Donnelly MA, Spande TF, Garraffo HM. A review of chemical ecology in poison frogs. Chemoecology. 2012;22(3):159–68.
Gunzburger MS, Travis J. Critical literature review of the evidence for unpalatability of amphibian eggs and larvae. J Herpetol. 2005;39(4):547–71.
Semlitsch RD, Gibbons JW. Fish predation in size-structured populations of treefrog tadpoles. Oecologia. 1988;75(3):321–6.
Richards SJ, Bull CM. Size-limited predation on tadpoles of three Australian frogs. Copeia. 1990;1990(4):1041–6.
Eklöv P, Werner EE. Multiple predator effects on size-dependent behavior and mortality of two species of anuran larvae. Oikos. 2000;88(2):250–8.
Wilson RS, Franklin CE. Effect of ontogenetic increases in body size on burst swimming performance in tadpoles of the striped marsh frog, Limnodynastes peronii. Physiol Biochem Zool. 2000;73(2):142–52.
Hayes RA, Crossland MR, Hagman M, Capon RJ, Shine R. Ontogenetic variation in the chemical defenses of cane toads (Bufo marinus): toxin profiles and effects on predators. J Chem Ecol. 2009;35(4):391–9.
Saporito RA, Isola M, Maccachero VC, Condon K, Donnelly MA. Ontogenetic scaling of poison glands in a dendrobatid poison frog. J Zool. 2010;282:238–45.
Stynoski JL, Torres-Mendoza Y, Sasa-Marin M, Saporito RA. Evidence of maternal provisioning of alkaloid-based chemical defenses in the strawberry poison frog Oophaga pumilio. Ecology. 2014;95(3):587–93.
Delfino G, Brizzi R, Feri L. Chemical skin defence in Bufo bufo: an ultrastructural study during ontogenesis. Zool Anz. 1995;234:101–11.
Chammas SM, Carneiro SM, Ferro RS, Antoniazzi MM, Jared C. Development of integument and cutaneous glands in larval, juvenile and adult toads (Rhinella granulosa): a morphological and morphometric study. Acta Zool-Stockholm. 2014;96(4):460–77.
Gall BG, Stokes AN, French SS, Schlepphorst EA, Brodie ED III, Brodie ED Jr. Tetrodotoxin levels in larval and metamorphosed newts (Taricha granulosa) and palatability to predatory dragonflies. Toxicon. 2011;57:978–83.
Hagman M, Hayes RA, Capon RJ, Shine R. Alarm cues experienced by cane toad tadpoles affect post-metamorphic morphology and chemical defences. Funct Ecol. 2009;23(1):126–32.
Marion ZH, Fordyce JA, Fitzpatrick BM. Extending the concept of diversity partitioning to characterize phenotypic complexity. Am Nat. 2015;186(3):348–61.
Brossman KH, Carlson BE, Stokes AN, Langkilde T. Eastern newt (Notophthalmus viridescens) larvae alter morphological but not chemical defenses in response to predator cues. Can J Zool. 2014;92(4):279–83.
Laurila A, Kujasalo J, Ranta E. Predator-induced changes in life history in two anuran tadpoles: effects of predator diet. Oikos. 1998;83:307–17.
Lardner B. Morphological and life history responses to predators in larvae of seven anurans. Oikos. 2000;88:169–80.
Van Buskirk J. A comparative test of the adaptive plasticity hypothesis: relationships between habitat and phenotype in anuran larvae. Am Nat. 2002;160(1):87–102.
Denton J, Beebee TJC. Palatability of anuran eggs and embryos. Amphibia-Reptilia. 1991;12(1):111–2.
Henrikson B-I. Predation on amphibian eggs and tadpoles by common predators in acidified lakes. Holarct Ecol. 1990;13:201–6.
Mebs D, Wagner MG, Pogoda W, Maneyro R, Kwet A, Kauert G. Lack of bufadienolides in the skin secretion of red bellied toads, Melanophryniscus spp. (Anura, Bufonidae), from Uruguay. Comp Biochem Physiol C Toxicol Pharmacol. 2007;144:398–402.
Flier J, Edwards MW, Daly JW, Myers CW. Widespread occurrence in frogs and toads of skin compounds interacting with the ouabain site of na+,K+-ATPase. Science. 1980;208(4443):503–5.
Gao H, Zehl M, Leitner A, Wu X, Wang Z, Kopp B. Comparison of toad venoms from different Bufo species by HPLC and LC-DAD-MS/MS. J Ethnopharmacol. 2010;131(2):368–76.
Krenn L, Kopp B. Bufadienolides from animal and plant sources. Phytochemistry. 1998;48(1):1–29.
Bókony V, Móricz ÁM, Tóth Z, Gál Z, Kurali A, Mikó Z, Pásztor K, Szederkényi M, Tóth Z, Ujszegi J, et al. Variation in chemical defense among natural populations of common toad, Bufo bufo, tadpoles: the role of environmental factors. J Chem Ecol. 2016;42(4):329–38.
APHA. Standard methods for the examination of wastewater. 16th ed. Washington DC: American Public Health Association; 1985.
Üveges B, Fera G, Móricz ÁM, Krüzselyi D, Bókony V, Hettyey A (2017) Data from: Age- and environment-dependent changes in chemical defences of larval and post-metamorphic toads. figshare. http://dx.doi.org/10.6084/m9.figshare.4635310
R Development Core Team: A language and environment for statistical computing. R package version 3.1.3. Vienna, Austria: R Foundation for Statistical Computing; 2015.
Pinheiro J, Bates D, DebRoy S, Sarkar D, R Development Core Team: nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.1-128. 2015. http://CRAN.R-project.org/package=nlme. Accessed 12 Feb 2017.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B Stat Method. 1995;57(1):289–300.
Lenth RV. Least-squares means: the R package lsmeans. J Stat Softw. 2016;69(1):1–33.
Lopanik NB, Targett NM, Lindquist N. Ontogeny of a symbiont-produced chemical defense in Bugula neritina (bryozoa). Mar Ecol Prog Ser. 2006;327:183–91.
Camarano S, González A, Rossini C. Chemical defense of the ladybird beetle Epilachna paenulata. Chemoecology. 2006;16(4):179–84.
Lindquist N, Hay ME. Palatability and chemical defense of marine invertebrate larvae. Ecol Monogr. 1996;66(4):431–50.
Brodie ED, Formanowicz DR. Larvae of the predaceous diving beetle Dytiscus verticalis acquire an avoidance response to skin secretions of the newt Notophthalmus viridescens. Herpetologica. 1981;37(3):172–6.
Kruse KC, Stone BM. Largemouth bass (Micropterus salmoides) learn to avoid feeding on toad (Bufo) tadpoles. Anim Behav. 1984;32(4):1035–9.
Nelson DWM, Crossland MR, Shine R. Foraging responses of predators to novel toxic prey: effects of predator learning and relative prey abundance. Oikos. 2011;120(1):152–8.
Cunha Filho GA, Schwartz CA, Resck IS, Murta MM, Lemos SS, Castro MS, Kyaw C, Pires OR Jr, Leite JR, Bloch C Jr, et al. Antimicrobial activity of the bufadienolides marinobufagin and telocinobufagin isolated as major components from skin secretion of the toad Bufo rubescens. Toxicon. 2005;45(6):777–82.
Tempone AG, Pimenta DC, Lebrun I, Sartorelli P, Taniwaki NN, de Andrade Jr HF, Antoniazzi MM, Jared C. Antileishmanial and antitrypanosomal activity of bufadienolides isolated from the toad Rhinella jimi parotoid macrogland secretion. Toxicon. 2008;52(1):13–21.
Arnold SJ, Wassersug RJ. Differential predation on metamorphic anurans by garter snakes (Thamnophis): social behavior as a possible defense. Ecology. 1978;59(5):1014–22.
Wassersug RJ, Sperry DG. The relationships of locomotion to differential predation on Pseudacris triseriata (Anura: Hylidae). Ecology. 1977;58(4):830–9.
Crump ML. Ontogenetic changes in vulnerability to predation in tadpoles of Hyla pseudopuma. Herpetologica. 1984;40(3):265–71.
Calsbeek R, Kuchta S. Predator mediated selection and the impact of developmental stage on viability in wood frog tadpoles (Rana sylvatica). BMC Evol Biol. 2011;11:353.
Touchon JC, Jiménez RR, Abinette SH, Vonesh JR, Warkentin KM. Behavioral plasticity mitigates risk across environments and predators during anuran metamorphosis. Oecologia. 2013;173(3):801–11.
Beck CW, Congdon JD. Energetics of metamorphic climax in the southern toad (Bufo terrestris). Oecologia. 2003;137(3):344–51.
Orlofske SA, Hopkins WA. Energetics of metamorphic climax in the pickerel frog (Lithobates palustris). Comp Biochem Physiol A Mol Integr Physiol. 2009;154(2):191–6.
Highnam KC: A survey of invertebrate metamorphosis. In: Metamorphosis: a problem in developmental biology. Edited by Gilbert LI, Frieden E, 2nd edn. New York, USA: Plenum Press; 1981: 43-73.
Almeida E, Nunes A, Andrade P, Alves S, Guerreiro C, Rebelo R. Antipredator responses of two anurans towards native and exotic predators. Amphibia-Reptilia. 2011;32(3):341–50.
Wells KD. The ecology and behavior of amphibians. Chicago: The University of Chicago Press; 2007.
Pfennig DW, Wund MA, Snell-Rood EC, Cruickshank T, Schlichting CD, Moczek AP. Phenotypic plasticity’s impacts on diversification and speciation. Trends in Ecology and Evolution. 2010;25(8):459–67.
Kurali A, Pásztor K, Hettyey A, Tóth Z. Toxin depletion has no effect on antipredator responses in common toad (Bufo bufo) tadpoles. Biol J Linn Soc. 2016;119(4):1000–10.
Miele R, Ponti D, Boman HG, Barra D, Simmaco M. Molecular cloning of a bombinin gene from Bombina orientalis: detection of NF-κB and NF-IL6 binding sites in its promoter. FEBS Lett. 1998;431:23–8.
Mangoni ML, Miele R, Renda TG, Barra D, Simmaco M. The synthesis of antimicrobial peptides in the skin of Rana esculenta is stimulated by microorganisms. FASEB J. 2001;15:1431–2.
Steiner UK. Linking antipredator behaviour, ingestion, gut evacuation and costs of predator-induced responses in tadpoles. Anim Behav. 2007;74(5):1473–9.
Tollrian R, Harvell CD. The ecology and evolution of inducible defences. New Jersey: Princeton University Press; 1999.
Benard MF. Predator-induced phenotypic plasticity in organisms with complex life histories. Annu Rev Ecol Evol Syst. 2004;35:651–73.
Van Buskirk J, Saxer G. Delayed costs of an induced defense in tadpoles? Morphology, hopping, and development rate at metamorphosis. Evolution. 2001;55(4):821–9.
Reigosa MJ, Pedrol N, González L. Allelopathy: a physiological process with ecological implications. Dordrecht: Springer; 2006.
Metlen KL, Aschehoug ET, Callaway RM. Plant behavioural ecology: dynamic plasticity in secondary metabolites. Plant Cell Environ. 2009;32(6):641–53.
Sieg RD, Poulson-Ellestad KL, Kubanek J. Chemical ecology of the marine plankton. Nat Prod Rep. 2011;28(2):388–99.
Jackson JBC, Buss L. Alleopathy and spatial competition among coral reef invertebrates. Proc Natl Acad Sci U S A. 1975;72(12):5160–3.
Kubanek J, Whalen K, Engel S, Kelly S, Henkel T, Fenical W, Pawlik J. Multiple defensive roles for triterpene glycosides from two Caribbean sponges. Oecologia. 2002;131(1):125–36.
Crossland MR, Shine R. Embryonic exposure to conspecific chemicals suppresses cane toad growth and survival. Biol Lett. 2012;8(2):226–9.
Licht LE. Growth inhibition in crowded tadpoles: intraspecific and interspecific effects. Ecology. 1967;48(5):736–45.
The authors would like to thank B. Vági for capturing adult toads, D. Koska for his cheerful assistance during the experiment, J. Van Buskirk, Z. Tóth and all members of the MTA ATK NÖVI Evolutionary Ecology Research Group for their useful comments during the planning phase and V. Verebélyi for drawing Additional file 1: Figure S1.
Funding was provided by a Sparkling Science Project of the Federal Ministry of Science and Research, Austria (BMWF, SPA 04/171), an FP7 Marie Curie Career Integration Grant (PCIG13-GA-2013-631,722), and the “Lendület” programme of the Hungarian Academy of Sciences (MTA, LP2012-24/2012). VB was further supported by the National Research, Development and Innovation Office (NKFIH) of Hungary (grant no. 115402) and the János Bolyai Scholarship of the Hungarian Academy of Sciences.
Availability of data and materials
BÜ and AH designed the study, GF and AH conducted the experiment and took samples, BÜ prepared samples for chemical analysis, ÁMM and DK analysed samples, BÜ and VB conducted statistical analyses, BÜ, ÁMM, VB and AH wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
All experimental procedures were carried out according to the permits issued by the Közép-Duna-Völgyi KTVF (KTVF: 5192-7/2013) and the Government Agency of Pest County, Hungary (PEI/001/389-4/2013). The experiment was further approved by the Ethical Commission of MTA ATK NÖVI. Consent to participate not applicable.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Üveges, B., Fera, G., Móricz, Á.M. et al. Age- and environment-dependent changes in chemical defences of larval and post-metamorphic toads. BMC Evol Biol 17, 137 (2017). https://doi.org/10.1186/s12862-017-0956-5