Competition between the tadpoles of Japanese toads versus frogs

Abstract

Competition within and among species can play a key role in structuring the assemblages of anuran tadpoles. Previous studies have reported that tadpoles of the invasive cane toad (Rhinella marina) are more strongly disadvantaged by the presence of native frog tadpoles than by the same number of conspecific toad tadpoles. That effect might arise from a lack of coevolution of the invasive toad with its competitors; and/or from a generalized superiority of frog tadpoles over toad tadpoles. To clarify those possibilities, we conducted experimental trials using the larvae of a native rather than invasive toad (Bufo japonicus formosus in Japan) exposed to larvae of native anurans (the sympatric frogs Rana japonica and Rana ornativentris and the parapatric toad Bufo japonicus japonicus). In intraspecific competition trials, higher densities of B. j. formosus prolonged the larval period and reduced size at metamorphosis, but did not affect survival. In interspecific competition trials, the effects of the other anuran species on B. j. formosus were similar to the effects of the same number of conspecific larvae. This similarity in impact of interspecific versus intraspecific competition argues against any overall competitive superiority of frog larvae over toad larvae. Instead, the vulnerability of larval cane toads to frog tadpoles may result from a lack of coevolutionary history.

Download PDF

Introduction

Many ecosystems contain a diverse array of species that overlap considerably in the resources upon which they rely; and ecological theory suggests that such overlap may lead to intense competition, that in turn may favour adaptations to reduce niche overlap^1,2,3. The larvae of anuran amphibians (tadpoles) have been "model organisms" for experimental studies of competitive effects, because a single waterbody often contains larvae of multiple taxa with high niche overlap, that compete for finite food resources^4,5,6. Extensive studies have confirmed that intra- and interspecific competition can affect larval survival rates, larval periods, and body mass at metamorphosis^{6,7,8,9,10,11,12,13}, and lead to the evolution of mechanisms for suppression of competing larvae (e.g.^14,15).

One interesting example of competitive suppression involves the invasive cane toad (Rhinella marina). Experimental studies in two parts of the toads' invasive range (Australia and Ishigaki Island, Japan) have shown that the survival and growth of larval cane toads is strongly reduced by the presence of frog tadpoles^16,17,18,19, via exploitative competition for food²⁰. In these studies, the viability of a cane toad tadpole was reduced more by competition with a frog tadpole than with a conspecific toad tadpole—but why? Possible answers include a generalised competitive superiority of frog larvae over toad larvae, perhaps because the former are often larger than the latter; and/or an effect of coevolution, whereby cane toads are highly vulnerable because they have encountered these frog species only recently (for decades, at most) and hence have not yet evolved mechanisms to reduce that vulnerability.

To explore this question, we can examine the effects of intraspecific and interspecific competition on a toad species that is native rather than invasive—and hence, has had the opportunity to adapt to sympatric frogs over evolutionary time. We can also examine the sensitivity of such a toad to competition from closely related bufonid species that are parapatric to the target taxon, to compare competitive impacts of toads versus frogs. No such studies were possible in earlier studies of cane toads in Australia and on Ishigaki Island (Okinawa, Japan) because neither of those areas contains native toads.

We have conducted such a study using a toad taxon that is native to Japan (Bufo japonicus formosus) as our target species, two ranid frogs as sympatric competitors (Rana japonica, Rana ornativentris) and a closely related subspecies as the parapatric toad competitor (Bufo japonicus japonicus). All four species breed at the same time of year, such that the tadpoles of B. j. formosus co-occur with the two Rana species in temporary waterbodies. We manipulated the numbers and identity of tadpoles in experimental containers to (1) quantify the effects of intraspecific competition, by raising tadpoles of B. j. japonicus at a range of densities; and (2) quantify the effects of interspecific competition, by raising tadpoles at a fixed total larval density but changing the composition of the assemblage in terms of the relative numbers of toad versus frog competitors.

Methods

Study species

Our laboratory studies included four anuran taxa, all of which are terrestrial and breed in wide array of freshwater habitats including ponds, marshes and swamps. Two species belong to the "true toads" (family Bufonidae). Our target species was the eastern-Japanese common toad (Bufo japonicus formosus; total length [= TL] of tadpoles up to 30 mm²¹), and the parapatric bufonid was the western-Japanese common toad (Bufo japonicus japonicus; tadpole TL generally up to 35 mm²¹). The other two taxa are members of the family Ranidae, both of which are broadly sympatric with B. j. formosus: the Japanese brown frog (Rana japonica; mean tadpole TL = 38 mm²¹) and the montane brown frog (Rana ornativentris; mean tadpole TL = 43 mm²¹). Larval body sizes in the two groups used in the experiment spanned a similar range (Table 1).

Table 1 Body sizes (mean ± standard errors, and range) and Gosner stages for tadpoles as measured at the beginning of the experiments.

Full size table

Tadpoles of all four species were derived from eggs collected in natural waterbodies from two sites (B. j. formosus and R. japonica—Tochigi prefecture, B. j. japonicus and R. ornativentris—Okayama prefecture) during the period 15–31 March 2014. Tadpoles of B. j. formosus, R. japonica and R. ornativentris were found in the same waterbodies (Haramura, personal observation), confirming that competition is likely to occur in nature. To equalize developmental stage of tadpoles at the onset of the experiment as much as possible, eggs or embryos of early-laid clutches were kept in cool conditions (12 °C) prior to the main experiment. Tadpoles of all four species were maintained in groups in 120 L plastic containers (66 × 86 × 34 cm). Tadpoles were fed algal pellets (Hikari Algae Wafers, Kyorin) ad libitum, and water was changed weekly. The tadpoles used in the experiment were haphazardly selected from these containers and added to experimental bins as described below.

Laboratory experiments

Experiments were conducted using plastic tanks (26 × 38 × 23 cm), each filled with 23 L water and located in a covered building exposed to ambient temperatures. At the start of the experiment, we added a 2 cm layer of soil substrate and 3 g of algal pellets to each bin. We did not provide additional food for the remainder of the experiment. Tadpoles varied in sizes and developmental stages at the beginning of the experiment (see Table 1).

Our experiment consisted of six treatments: (1) 5 larvae of B. j. formosus, (2) 15 B. j. formosus, (3) 50 B. j. formosus, (4) 25 B. j. formosus plus 25 B. j. japonicus, (5) 25 B. j. formosus plus 25 R. japonica, and (6) 25 B. j. formosus plus 25 R. ornativentris. The experiment was a complete randomised block design, with 5 replicate tanks per treatment. We recorded the number of B. j. formosus to metamorphose from each tank (survival), as well as the larval period, and length (snout to urostyle length = SUL) and mass of each B. j. formosus metamorph from each tank. Treatments 1, 2 and 3 allowed us to assess the effect of intraspecific competition on B. j. formosus, whereas treatments 3, 4, 5 and 6 allowed us to assess the strength of interspecific versus intraspecific competition at standardised density. We also measured water temperature and pH in each tank every 4 days.

Because it is not possible to visually distinguish between metamorphs of B. j. formosus and B. j. japonicus, we used a Loop-Mediated Isothermal Amplification (LAMP) assay to distinguish between these two subspecies in the interspecific competition experiments. LAMP is a genetic method which detects the presence/absence of a specific DNA sequence in the tested sample²². Total DNA of each metamorph was extracted from frozen tissue using the DNeasy Blood and Tissue Kit (QIAGEN) with standard protocols. Following extraction, each sample was tested by LAMP assay in two independent systems—assays with B. j. japonicus-positive and B. j. formosus-positive primer sets. For the primer set used and the experimental conditions, we followed the methods with slight modification²³. The reaction mixtures were incubated at 63–65 °C for 90 min and then heated at 95 °C for 2 min to terminate the reaction.

Statistical analyses

We analysed treatment effects on both water temperature and pH using ANOVA. We analysed treatment effects on larval period, metamorph SUL and metamorph mass using linear models (MANOVA, followed by ANOVA), with treatment and spatial block as fixed effects (JMP 9.0, SAS Institute, Cary, NC, USA). MANOVA analyses were based on tank means to avoid pseudoreplication (the JMP statistical package does not support MANOVA with random effects). ANOVA analyses were based on data for all individuals per tank, using tank as a random effect. When the overall ANOVA gave a significant result, we performed post hoc Tukey's HSD tests for pairwise comparison of treatments. We analysed survival to metamorphosis as a binomial response (alive, dead²⁴) using ANOVA, with treatment and spatial block as fixed effects (package carData^25,26). Survival analyses were based on the quasi-binomial distribution to account for overdispersion of data. Alpha level was set at p = 0.05 in all analyses.

Ethics approval

All procedures were approved by Rakuno Gakuen University Animal Care Committee (permit #DH21D6). The study was carried out in compliance with the ARRIVE guidelines, and all methods were carried out in accordance with relevant guidelines and regulations.

Results

The average water temperature and pH in tanks was 19.29 ± 0.10 °C (SE, range: 17.0–22.5) and 8.59 ± 0.01 (SE, range 8.2–8.9) respectively. There was no significant difference among treatments (water temperature: F = 0.0086, df = 5, p = 1.0000, pH: F = 0.0063, df = 5, p = 1.0000).

Intraspecific competition (density = 5, 15, 50 tadpoles per tank)

The density of conspecifics did not have any significant effect on survival to metamorphosis of B. j. formosus (treatment: Wald chi-square = 3.468, df = 2, p = 0.1766; block: Wald chi-square = 7.770, df = 4, p = 0.1004; Fig. 1a). However, conspecific density had a significant effect on the combined responses of variables (larval period, metamorph SUL, metamorph mass) of B. j. formosus (MANOVA treatment: Wilks' Lambda = 0.0181, F = 10.7224, df = 6, 10, p = 0.0007; block: Wilks' Lambda = 0.2028, F = 0.9326, df = 12, 13.52, p = 0.5441). Higher densities of conspecifics increased the duration of the larval period (treatment: F = 6.678, df = 2, 9.30, p = 0.0159; block: F = 0.817, df = 4, 0.40, p = 0.7574; Fig. 1b), and decreased size at metamorphosis (SUL—treatment: F = 49.729, df = 2, 6.94, p < 0.0001; block: F = 1.154, df = 4, 6.88, p = 0.4074; Fig. 1c; mass—treatment: F = 22.949, df = 2, 6.66, p = 0.0010; block: F = 1.031, df = 4, 6.68, p = 0.4566; Fig. 1 d).

Interspecific competition (density = 50 tadpoles per tank)

There was no significant effect of treatment (competitor species) on survival of B. j. formosus to metamorphosis (treatment: Wald chi-square = 4.076, df = 3, p = 0.2533; block: Wald chi-square = 2.708, df = 4, p = 0.6078; Fig. 2a). MANOVA also showed no significant effect of treatment on overall responses (i.e., including variables of larval period, metamorph SUL, metamorph mass) of B. j. formosus (treatment: Wilks' Lambda = 0.3722, F = 1.2275, df = 9, 22.05, p = 0.3285; block: Wilks' Lambda = 0.3269, F = 1.0565, df = 12, 24.10, p = 0.4344). Although tadpoles of B. j. japonicus tended to impose stronger negative effects on B. j. formosus than did B. j. formosus on itself (Fig. 2b–d), there was also no significant effect of competitor species on the duration of the larval period for B. j. formosus (treatment: F = 2.262, df = 3, 9.83, p = 0.1448; block F = 0.783, df = 4, 9.56, p = 0.5627, Fig. 2b), or size at metamorphosis (SUL—treatment: F = 1.895, df = 3, 10.46, p = 0.1917; block: F = 2.039, df = 4, 10.46, p = 0.1615; Fig. 2c; mass—treatment: F = 2.706, df = 3, 10.69, p = 0.0980; block: F = 0.495, df = 4, 10.68, p = 0.7403; Fig. 2d).

Discussion

In our laboratory experiment, tadpoles of the eastern-Japanese common toad (Bufo japonicus formosus) showed strong intraspecific competitive effects: an increase in the number of toad tadpoles per container generated a substantial reduction in rates of growth and development, and in size at metamorphosis (Fig. 1). The main result from interspecific competition treatments, however, was that these effects did not vary among competitors that were either conspecifics or heterospecifics (Fig. 2). That is, the impacts of frog tadpoles on B. j. formosus were similar to those of the same number of toad tadpoles (Fig. 2).

Our results accord with previous studies that have shown negative consequences for tadpoles raised at high densities (e.g.^7,27,28,29). However, we did not find a stronger competitive effect of frog tadpoles than of toad tadpoles, unlike the results of studies on invasive cane toads in Australia and Okinawa^16,17,18,19. Why, then, are tadpoles of the cane toad more sensitive to the presence of frog tadpoles than conspecific toad tadpoles? At least under the conditions under which we conducted our experiments, the answer does not involve a competitive superiority of frog tadpoles over toad tadpoles: we saw no such effect in our trials (Fig. 2). Instead, the results for cane toads may reflect two aspects of this system. First, most of the Australian frog tadpoles tested were much larger than the toad tadpoles—in some cases, by a 20-fold margin^16,17,18. Larger tadpoles may (in general) be better competitors (e.g.³⁰), and this effect may be stronger if the size disparity is greater. However, we note that the competitive superiority of native frog tadpoles over invasive cane toads on Ishigaki occurred despite a relatively small difference in body size (native tadpoles up to 1.44 times the size of cane toad tadpoles¹⁹). In the present study, interspecific differences in tadpole size (Table 1) did not translate to differential competitive effects (Fig. 2), although size differences were relatively minor. Thus, at least over the size range studied in Ishigaki and in the current study, size effects seem unlikely to explain differential vulnerability of cane toad larvae to frog tadpoles than to conspecifics. And as well as body size, outcomes of competition also may be affected by habitat use. However, the tadpoles of all four species used in the present study are primarily bottom-dwellers, minimizing any effects of differential habitat use among species.

The second aspect of the cane toad system is that this is an invasive species; and hence, these toads have had only a brief window of opportunity to adapt in ways that buffer them against the competitive effects of native anurans. In contrast, the toad species that we targeted in the present study (B. j. formosus) is sympatric with the frogs we studied (R. japonica R. ornativentris) over a broad area, and thus likely has coevolved with those frogs over a long period^31,32. These taxa frequently breed in the same waterbodies, at the same time of year, and thus compete with each other in nature as well as in our laboratory studies. That situation, continuing over long periods, should enable coevolution between the competing taxa, in ways that reduce the negative impacts of competitors. Interestingly, we found no significant competitive effects of parapatric B. j. japonicus on B. j. formosus. This result may be due to the close phylogenetic relationship (i.e., belonging to the same genus), and thus ecological similarity, between these two species.

Future work could usefully examine competitive interactions between adult anurans as well as between larvae; and could assess the impacts of a broader range of species under a wider range of conditions (including, outdoor enclosures that more accurately mimic spawning sites in nature). In particular, it would be of great interest to examine larval competition within the native range of the cane toad. If this species' vulnerability to competition from frog tadpoles results from lack of coevolution in invaded areas, then we expect that trials with the tadpoles of South American frog species would provide different results to those seen in Australia and on Ishigaki. That is, cane toad tadpoles should be resilient to the presence of larvae from sympatric frog species, as seen in our work with B. j. formosus. More generally, it would be instructive to compare ecological interactions between invasive species and other fauna not only in the areas they have invaded, but also within their native range, to clarify the impacts of translocation on the intensity of interspecific competition.

References

Colwell, R. K. & Futuyma, D. J. On the measurement of niche breadth and overlap. Ecology 52, 567–576 (1971).
Article  Google Scholar 
Pianka, E. R. Niche overlap and diffuse competition. Proc. Natl. Acad. Sci. USA 71, 2141–2145 (1974).
ADS  CAS  Article  Google Scholar 
Toft, C. Resource partitioning in amphibians and reptiles. Copeia 1985, 1–21 (1985).
Article  Google Scholar 
Alford, R. A. & Wilbur, H. M. Priority effects in experimental pond communities: Competition between Bufo and Rana. Ecology 66, 1097–1105 (1985).
Article  Google Scholar 
Altig, R. & Johnston, G. F. Guilds of anuran larvae: Relationships among developmental modes, morphologies, and habitats. Herp. Monogr. 3, 81–109 (1989).
Article  Google Scholar 
Alford, R. A. Ecology: Resource use, competition, and predation. In Tadpoles: The Biology of Anuran Larvae (eds McDiarmid, R. W. & Altig, R.) 240–278 (University of Chicago Press, 1999).
Google Scholar 
Wilbur, H. M. Density-dependent aspects of growth and metamorphosis in Bufo americanus. Ecology 58, 196–200 (1977).
Article  Google Scholar 
Petranka, J. W. Density-dependent growth and survival of larval Ambystoma: Evidence from whole-pond manipulations. Ecology 70, 1752–1767 (1987).
Article  Google Scholar 
Goater, C. P. Growth and survival of postmetamorphic toads: Interactions among larval history, density, and parasitism. Ecology 75, 2264–2274 (1994).
Article  Google Scholar 
Altwegg, R. Multistage density dependence in an amphibian. Oecologia 136, 46–50 (2003).
ADS  Article  Google Scholar 
Relyea, R. A. Fine-tuned phenotypes: Tadpole plasticity under 16 combinations of predator and competitors. Ecology 85, 172–179 (2004).
Article  Google Scholar 
Boone, M. D. Juvenile frogs compensate for small metamorph size with terrestrial growth: Overcoming the effects of larval density and insecticide exposure. J. Herpetol. 39, 416–423 (2005).
Article  Google Scholar 
Yagi, K. T. & Green, D. M. Mechanisms of density-dependent growth and survival in tadpoles of Fowler's toad, Anaxyrus fowleri: Volume vs. abundance. Copeia 104, 942–951 (2016).
Article  Google Scholar 
Griffiths, R. A., Edgar, P. W. & Wong, A.L.-C. Interspecific competition in tadpoles: Growth inhibition and growth retrieval in natterjack toads, Bufo calamita. J. Anim. Ecol. 60, 1065–1076 (1991).
Article  Google Scholar 
Griffiths, R. A., Denton, J. & Wong, A.L.-C. The effect of food level on competition in tadpoles: Interference mediated by protothecan algae?. J. Anim. Ecol. 62, 274–279 (1993).
Article  Google Scholar 
Crossland, M. R., Alford, R. A. & Shine, R. Impact of the invasive cane toad (Bufo marinus) on an Australian frog (Opisthodon ornatus) depends on minor variation in reproductive timing. Oecologia 158, 625–632 (2009).
ADS  Article  Google Scholar 
Cabrera-Guzmán, E., Crossland, M. R. & Shine, R. Can we use the tadpoles of Australian frogs to reduce recruitment of invasive cane toads?. J. Anim. Ecol. 62, 462–470 (2011).
Article  Google Scholar 
Cabrera-Guzmán, E., Crossland, M. R. & Shine, R. Competing tadpoles: Australian native frogs affect invasive cane toads (Rhinella marina) in natural waterbodies. Austral. Ecol. 38, 896–904 (2013).
Article  Google Scholar 
Haramura, T., Takeuchi, H., Crossland, M. R. & Shine, R. Biotic resistance to an alien amphibian: Larval competition between Japanese frogs and invasive cane toads. PLoS ONE 11, e0156396 (2016).
Article  Google Scholar 
Cabrera-Guzmán, E., Crossland, M. R. & Shine, R. Mechanisms of competition between tadpoles of Australian frogs (Litoria spp.) and invasive cane toads (Rhinella marina). Freshw. Biol. 58, 2584–2600 (2013).
Article  Google Scholar 
Maeda, N. & Matsui, M. Frogs and Toads of Japan Revised. (Bun-ichi Sogo Shuppan, 1999) ((in Japanese)).
Google Scholar 
Notami, T. et al. Loop-mediated isothermal amplification of DNA. Nucleic Acids Res. 28, e63 (2000).
Article  Google Scholar 
Eto, K., Matsui, M., Nishikawa, K. & Haramura, T. Development and evaluation of loop-mediated isothermal amplification (LAMP) assay for quick identification of three Japanese toads. Curr. Herpetol. 35, 33–37 (2016).
Article  Google Scholar 
Warton, D. & Hui, F. K. C. The arcsine is asinine: The analysis of proportions in ecology. Ecology 92, 3–10 (2011).
Article  Google Scholar 
Fox, J., Weisberg, S. & Price, B. carData: Companion to Applied Regression Data Sets. R package version 3.0-2. https://CRAN.R-project.org/package=carData (2018).
R Core Team R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2019).
Brockelman, Y. W. An analysis of density effects and predation in Bufo americanus tadpoles. Ecology 50, 632–644 (1969).
Article  Google Scholar 
Warner, S. C., Dunson, W. A. & Travis, J. Interaction of pH, density, and priority effects on the survivorship and growth of two species of hylid tadpoles. Oecologia 88, 331–339 (1991).
ADS  Article  Google Scholar 
Relyea, R. A. & Auld, J. R. Predator- and competitor-induced plasticity: How changes in foraging morphology affect phenotypic trade-offs. Ecology 86, 1723–1729 (2005).
Article  Google Scholar 
Bardsley, L. & Beebee, T. J. C. Strength and mechanisms of competition between common and endangered anurans. Ecol. Appl. 11, 453–463 (2001).
Article  Google Scholar 
Igawa, T., Kurabayashi, A., Nishioka, M. & Sumida, M. Molecular phylogenetic relationship of toads distributed in the Far East and Europe inferred from the nucleotides sequences of mitochondrial DNA genes. Mol. Phylogenet. Evol. 38, 250–260 (2006).
CAS  Article  Google Scholar 
Yuan, Z. et al. Spatiotemporal diversification of the True Frogs (Genus Rana): A historical framework for a widely studied group of model organisms. Syst. Biol. 65, 824–842 (2016).
Article  Google Scholar 
Gosner, K. L. A simplified table for staging anuran embryos and larvae with notes on identification. Herpetologica 16, 183–190 (1960).
Google Scholar

Download references

Acknowledgements

We thank Gregory Brown (Macquarie University) for statistical support and suggestions. This work was supported by an internal grant from The Hakubi Center, Kyoto University, the Environment Research and Technology Develoment Fund (4RF-1402) of the Ministry of the Environment, Japan.

Author information

Affiliations

Department of Environmental Sciences, Rakuno Gakuen University, Hokkaido, Japan
Takashi Haramura
Graduate School of Human and Environmental Studies, Kyoto University, Kyoto, Japan
Koshiro Eto & Kanto Nishikawa
Kitakyushu Museum of Natural History & Human History, Kitakyushu, Fukuoka, Japan
Koshiro Eto
School of Life and Environmental Sciences, University of Sydney, Camperdown, NSW, 2006, Australia
Michael R. Crossland & Richard Shine
Graduate School of Global Environmental Studies, Kyoto University, Kyoto, Japan
Kanto Nishikawa
School of Natural Sciences, Macquarie University, North Ryde, NSW, 2109, Australia
Richard Shine

Authors

Takashi Haramura
View author publications
You can also search for this author in PubMed Google Scholar
Koshiro Eto
View author publications
You can also search for this author in PubMed Google Scholar
Michael R. Crossland
View author publications
You can also search for this author in PubMed Google Scholar
Kanto Nishikawa
View author publications
You can also search for this author in PubMed Google Scholar
Richard Shine
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

T.H., and M.R.C. designed research experiment; T.H. and K.E. performed research; T.H., M.R.C., and R.S. analyzed data; T.H., M.R.C., K.N., and R.S. wrote the paper.

Corresponding author

Correspondence to Takashi Haramura.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Search This Blog

Famous Reptiles