- Research article
- Open Access
Age- and environment-dependent changes in chemical defences of larval and post-metamorphic toads
© The Author(s). 2017
Received: 12 February 2017
Accepted: 25 April 2017
Published: 13 June 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–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–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–32]. Moreover, there are only a handful of studies on phenotypic plasticity in chemical defences in amphibians [12, 33–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–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–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
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).
Percent occurrence, retention time and mass signal of bufadienolides in common toad tadpoles
Percent occurrence of bufadienolide compounds
retention time (min)
m/z [M + H]+
unidentified bufadienolide 1
unidentified bufadienolide 2
unidentified bufadienolide 3
unidentified bufadienolide 4
unidentified bufadienolide 5
unidentified bufadienolide 6
unidentified bufadienolide 7
unidentified bufadienolide 8
unidentified bufadienolide 9
unidentified bufadienolide 10
unidentified bufadienolide 11
unidentified bufadienolide 12
unidentified bufadienolide 13
unidentified bufadienolide 14
unidentified bufadienolide 15
unidentified bufadienolide 16
Number of bufadienolide compounds
Effects of ontogeny, treatments, their interactions, and body mass on bufadienolide synthesis of common toads
Number of bufadienolide compounds (NBC)
developmental stage × food level
developmental stage × predation treatment
food level × predation treatment
developmental stage × food level × predation treatment
Total bufadienolide quantity (TBQ, ng)
developmental stage × food level
developmental stage × predation treatment
food level × predation treatment
developmental stage × food level × predation treatment
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–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–54]. Additionally, the ability of predators to distinguish and learn to avoid noxious prey [55–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–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–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–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.
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.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Brodie ED III. Toxins and venoms. Curr Biol. 2009;19(20):R931–5.View ArticlePubMedGoogle Scholar
- Mebs D. Toxicity in animals. Trends in evolution? Toxicon. 2001;39(1):87–96.View ArticlePubMedGoogle Scholar
- Toledo RC, Jared C: Cutaneous granular glands and amphibian venoms. Comp Biochem Physiol A Physiol 1995, 111(1):1-29.Google Scholar
- Blum MS. Chemical defenses of arthropods. New York: Academic Press Inc.; 1981.Google Scholar
- Tachibana K: Chemical defense in fishes. In: Bioorg Mar Chem. Edited by Scheuer PJ, vol. 2. Berlin Heidelberg, Germany: Springer Verlag; 1988: 117-138.Google Scholar
- Pawlik JR. Marine invertebrate chemical defenses. Chem Rev. 1993;93(5):1911–22.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Stankowich T. Armed and dangerous: predicting the presence and function of defensive weaponry in mammals. Adapt Behav. 2012;20(1):32–43.View ArticleGoogle Scholar
- Hettyey A, Tóth Z, Van Buskirk J. Inducible chemical defences in animals. Oikos. 2014;123:1025–8.View ArticleGoogle Scholar
- Benard MF, Fordyce JA. Are induced defenses costly? Consequences of predator-induced defenses in western toads, Bufo boreas. Ecology. 2003;84(1):68–78.View ArticleGoogle Scholar
- 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.Google Scholar
- 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.View ArticleGoogle Scholar
- Morgenstern D, King GF. The venom optimization hypothesis revisited. Toxicon. 2013;63:120–8.View ArticlePubMedGoogle Scholar
- Longson CG, Joss JMP. Optimal toxicity in animals: predicting the optimal level of chemical defences. Funct Ecol. 2006;20(4):731–5.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Brown DD, Cai L. Amphibian metamorphosis. Dev Biol. 2007;306(1):20–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Gosner KL. A simplified table for staging anuran embryos and larvae with notes on identification. Herpetologica. 1960;16(3):183–90.Google Scholar
- Saporito RA, Donnelly MA, Spande TF, Garraffo HM. A review of chemical ecology in poison frogs. Chemoecology. 2012;22(3):159–68.View ArticleGoogle Scholar
- Gunzburger MS, Travis J. Critical literature review of the evidence for unpalatability of amphibian eggs and larvae. J Herpetol. 2005;39(4):547–71.View ArticleGoogle Scholar
- Semlitsch RD, Gibbons JW. Fish predation in size-structured populations of treefrog tadpoles. Oecologia. 1988;75(3):321–6.View ArticlePubMedGoogle Scholar
- Richards SJ, Bull CM. Size-limited predation on tadpoles of three Australian frogs. Copeia. 1990;1990(4):1041–6.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Delfino G, Brizzi R, Feri L. Chemical skin defence in Bufo bufo: an ultrastructural study during ontogenesis. Zool Anz. 1995;234:101–11.Google Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Marion ZH, Fordyce JA, Fitzpatrick BM. Extending the concept of diversity partitioning to characterize phenotypic complexity. Am Nat. 2015;186(3):348–61.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Lardner B. Morphological and life history responses to predators in larvae of seven anurans. Oikos. 2000;88:169–80.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Denton J, Beebee TJC. Palatability of anuran eggs and embryos. Amphibia-Reptilia. 1991;12(1):111–2.View ArticleGoogle Scholar
- Henrikson B-I. Predation on amphibian eggs and tadpoles by common predators in acidified lakes. Holarct Ecol. 1990;13:201–6.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Krenn L, Kopp B. Bufadienolides from animal and plant sources. Phytochemistry. 1998;48(1):1–29.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- APHA. Standard methods for the examination of wastewater. 16th ed. Washington DC: American Public Health Association; 1985.Google Scholar
- Ü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.Google Scholar
- 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.Google Scholar
- Lenth RV. Least-squares means: the R package lsmeans. J Stat Softw. 2016;69(1):1–33.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Camarano S, González A, Rossini C. Chemical defense of the ladybird beetle Epilachna paenulata. Chemoecology. 2006;16(4):179–84.View ArticleGoogle Scholar
- Lindquist N, Hay ME. Palatability and chemical defense of marine invertebrate larvae. Ecol Monogr. 1996;66(4):431–50.View ArticleGoogle Scholar
- 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.Google Scholar
- Kruse KC, Stone BM. Largemouth bass (Micropterus salmoides) learn to avoid feeding on toad (Bufo) tadpoles. Anim Behav. 1984;32(4):1035–9.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Wassersug RJ, Sperry DG. The relationships of locomotion to differential predation on Pseudacris triseriata (Anura: Hylidae). Ecology. 1977;58(4):830–9.View ArticleGoogle Scholar
- Crump ML. Ontogenetic changes in vulnerability to predation in tadpoles of Hyla pseudopuma. Herpetologica. 1984;40(3):265–71.Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Beck CW, Congdon JD. Energetics of metamorphic climax in the southern toad (Bufo terrestris). Oecologia. 2003;137(3):344–51.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.View ArticleGoogle Scholar
- Wells KD. The ecology and behavior of amphibians. Chicago: The University of Chicago Press; 2007.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Steiner UK. Linking antipredator behaviour, ingestion, gut evacuation and costs of predator-induced responses in tadpoles. Anim Behav. 2007;74(5):1473–9.View ArticleGoogle Scholar
- Tollrian R, Harvell CD. The ecology and evolution of inducible defences. New Jersey: Princeton University Press; 1999.Google Scholar
- Benard MF. Predator-induced phenotypic plasticity in organisms with complex life histories. Annu Rev Ecol Evol Syst. 2004;35:651–73.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Reigosa MJ, Pedrol N, González L. Allelopathy: a physiological process with ecological implications. Dordrecht: Springer; 2006.View ArticleGoogle Scholar
- Metlen KL, Aschehoug ET, Callaway RM. Plant behavioural ecology: dynamic plasticity in secondary metabolites. Plant Cell Environ. 2009;32(6):641–53.View ArticlePubMedGoogle Scholar
- Sieg RD, Poulson-Ellestad KL, Kubanek J. Chemical ecology of the marine plankton. Nat Prod Rep. 2011;28(2):388–99.View ArticlePubMedGoogle Scholar
- Jackson JBC, Buss L. Alleopathy and spatial competition among coral reef invertebrates. Proc Natl Acad Sci U S A. 1975;72(12):5160–3.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Crossland MR, Shine R. Embryonic exposure to conspecific chemicals suppresses cane toad growth and survival. Biol Lett. 2012;8(2):226–9.View ArticlePubMedGoogle Scholar
- Licht LE. Growth inhibition in crowded tadpoles: intraspecific and interspecific effects. Ecology. 1967;48(5):736–45.View ArticleGoogle Scholar