Skip to main content

The impact of elevated temperature and CO2 on growth, physiological and immune responses of Polypedates cruciger (common hourglass tree frog)

Abstract

Background

Amphibians are one of the most susceptible groups to climate change as their development occurs in aquatic environments or in microhabitats with high humidity. Accordingly, our primary objective was to investigate the chronic physiological responses seen in early larval to adult stages of Polypedates cruciger (Common hourglass tree frog) to future climate change based on continuous exposure to elevated temperature and elevated CO2 -induced low water pH. Free-swimming and free-feeding tadpoles were observed until metamorphosis under four experimental treatments; two elevated temperatures, one elevated CO2 (reduced pH) and a control maintained at ambient temperature (29 °C ± 1 °C) and CO2 (pH = 7). Elevated temperature treatments were maintained at 32 °C ± 0.5 °C and 34 °C ± 0.5 °C to represent respectively, the future climate scenarios RCP2.6 (Representative Concentration Pathway 2.6, the ‘base-case’ scenario) and RCP8.5 (‘business-as-usual’ scenario) according to the 5th Assessment Report of the IPCC. Elevated CO2 treatment was maintained within the pH range of 5.5–5.6 representing the range expected between RCP8.5 and RCP2.6.

Results

Compared to the control, elevated CO2 accelerated phenological progression of tadpoles through Gosner stages, thus resulting in lower body size at metamorphosis. Both elevated temperatures significantly delayed the development and reduced the growth of tadpoles. 100% mortality was observed in 34 °C treatment before metamorphosis (before Gosner stage 36) while all the tadpoles died after metamorphosis (at Gosner stage 46) in 32 °C treatment. Elevated CO2 increased tadpole activity, in terms of their swimming speed, while both of the elevated temperatures reduced it compared to the control. Catalase activity increased at elevated CO2. Ammonia excretion by tadpoles was decreased by elevated CO2, but increased under temperature elevation. Both Elevated CO2 and temperature treatments reduced the white blood cell count and its percentage of thrombocytes. Percentages of lymphocytes, monocytes and neutrophils were increased at 32 °C, while lymphocyte percentage and lysozyme activity were increased at elevated CO2. Several deformities were observed in tadpoles at elevated temperature and CO2.

Conclusions

Elevated temperatures and reduced pH due to elevated CO2, being major features of climate change, increase the vulnerability of amphibians, who are already one of the most threatened vertebrate groups. Based on our observations on the model amphibian species P. cruciger, increased vulnerability to climate change occurs by reducing their growth, body size and motility while also reducing their immunity and inducing physical deformities. These impacts are highly-likely to reduce the foraging, competitive and reproductive capabilities in their natural habitats. We conclude further that even the ‘best-case’ scenario of future climate change can impose significant physiological impacts that could threaten amphibian populations on broader spatial and temporal scales.

Graphical abstract

Background

Anthropogenic emission of carbon dioxide (CO2) is widely recognized as the principal cause of rising atmospheric temperatures and long-term climate change [1]. According to the Intergovernmental Panel on Climate Change (IPCC)'s Representative Concentration Pathway 8.5 (RCP8.5), which represents the most fossil fuel-intensive developmental pathway, global atmospheric CO2 (Ca) concentration is projected to increase up to 970 ppm by 2100 [2]. Even under the most environmentally-friendly scenario of RCP2.6, Ca is predicted to increase up to 490 ppm by 2050 before stabilization and subsequent decrease towards the end of this century. Absorption of CO2 by water bodies such as oceans and lakes increases with increasing Ca leading to increased concentrations of dissolved CO2 in water and reduced pH. By 2100, pH of ocean surface water is projected to decrease by 0.065 and 0.31 under RCP2.6 and RCP8.5 respectively [1]. Reductions of similar magnitude have been predicted for the pH of freshwater bodies (e.g. Laurentian Great Lakes) as well when increasing atmospheric partial pressure of CO2 is considered as the only climate forcing [3]. However, pH of freshwater aquatic environments could decrease further due to CO2 generated from decomposing organic matter [4]. Increasing concentrations of CO2 and other greenhouse gases enhance the natural greenhouse effect and accelerate global warming [5]. Consequently for all four RCP scenarios, global mean surface temperature (GMST) is projected to increase by 0.3–0.7 °C during 2016–2035 relative to 1986–2005 [1]. Longer-term projections of warming diverge among RCPs, ranging from 0.3–1.7 °C for RCP2.6 to 2.6–4.8 °C for RCP8.5 by 2081–2100 relative to 1986–2005 [1].

Every biological process shows its optimal performance within a specific range of environmental conditions [6, 7]. Thus, shifts in environmental conditions due to climate change may exert a significant influence on biological systems, at the individual as well as ecosystem level. Although the predicted increase in temperature and decrease in pH due to climate change appear small in magnitude, they could have appreciable impacts on biological systems from species to ecosystem levels [8,9,10]. To assess the impacts of climate change on biological systems, knowledge on three major aspects is required. They are the current climatic conditions and future climatic scenarios, how close organisms are to their tolerance limits in nature; and the degree to which organisms are able to adjust or acclimatize their sensitivity to variations in key climatic parameters [11, 12]. Most research on the response of biological systems to climate change has been based on a few well-studied model organisms such as Drosophila sp. and fish species, particularly temperate fish such as salmons and summer flounder Paralichtys dentatus [13,14,15]. Nevertheless, several studies have focused on amphibians, arguably the most endangered group as one third of all extant amphibian species are threatened with extinction [16].

Inhabitation of temporary aquatic habitats is a major factor that makes amphibians highly susceptible to climate change [17]. Being the link of transition of life from water to land, it is a common characteristic of most amphibians to spend at least one stage of their life cycles in water. Often, it is one of the early life stages that is spent in aquatic habitats. Most amphibians live in close proximity to aquatic habitats and return to water for reproduction. Dramatic fluctuations of temperature and pH occur in the temporary and shallow microhabitats that most amphibians and other ectotherms dwell, causing changes in development rates and development cycles [18, 19]. These fluctuations have detrimental impacts on various biological processes at different levels of biological organization, ranging from specific biochemical reactions within an individual organism to ecological interactions of species and communities [7, 20, 21]. As such loss of conducive habitats due to climate change could be the most probable reason for the rapid decline of amphibian population during recent times [22].

Generally, the rate of physiological processes relate non-linearly to temperature and pH [23, 24]. Although organisms are capable of surviving chronically lethal temperatures and pH for short periods, their growth and development depend on a range of cellular processes that require a specific set of environmental conditions (e.g. temperature, pH etc.) over a longer period for optimum performance [25]. Disruption of key cellular processes due to the absence of favourable environmental conditions appear externally as reduced growth rates and activity, delayed or disrupted development, leading to increased vulnerability to predation, risk of infection and desiccation. Consequently, climate change is regarded as a principal cause of the decline of amphibian and other ectothermic populations worldwide [26, 27].

In living organisms, respiration is a continuous process and reactive oxygen species (ROS) are generated simultaneously. Commonly generated ROS include superoxide and hydroxyl radicals and hydrogen peroxide (H2O2) [28, 29]. Catalase is a universal enzyme present in all aerobic organisms. It is capable of breaking down H2O2 to water and oxygen [30]. All enzymes have an optimal activity range of temperature and pH, Hence, changes in temperature and pH in the external microhabitat could affect enzymatic antioxidant defence systems. Optimum temperature and pH ranges for catalase enzyme activity are species-specific [31, 32]. Accordingly, the activity of antioxidant defense enzymes such as catalase requires investigation in amphibians when they are exposed to increasing temperature and decreasing pH in their microhabitat.

Emergence of new amphibian diseases and increasing severity of prevalent diseases has been reported worldwide and environmental factors are suspected to be their cause [33, 34]. This also has probably contributed to the worldwide decline of amphibian populations [35,36,37,38]. Recent extinction of many tropical frog species could be due to increased infection risk associated with climate change [39]. Accordingly, we investigated the immune response to varying climate change parameters by quantifying the response of white blood cells (WBCs) and lysozyme activity to elevated temperature and CO2. White blood cells play a major role in determining the immunity of living systems [40]. Lysozyme is an enzyme involved in bacterial lysis. Therefore, lysozyme activity is an indicator of immune status of an organism [41]. Lysozyme is known to be stable in temperatures as high as 72 °C and has an optimum activity within the pH range from 6 to 9 [42, 43].

Most research on the impact of environmental variations on immunity of amphibians have involved lower temperatures such as the immune response of amphibians to winter and hibernation [44]. In addition, some research has focused on immune response of ectotherms such as fish to stress caused by handling, transportation, and overcrowding [45]. Apart from the work of Bibi and Christi [46], research on the impacts of elevated temperature and fluctuating acidity on immunity, especially of tropical amphibian species, is limited. Hence, the current study will be important in setting a baseline for tropical species.

Developmental and morphological abnormalities associated with viscera, mouth, eye, and limbs of amphibians are common when exposed to environmental stressors during their developmental stages [47]. Although most research has focused on deformities caused by exposure to toxicants such as pesticides and heavy metals [48, 49], impacts of increasing temperature and acidity on development of deformities is an area that has not been studied extensively.

Sri Lanka is a country rich in amphibian diversity, particularly in anuran diversity [50]. So far, 103 species have been described out of which 87 are endemic to the island [50]. Being a tropical country, it is possible that impacts of climate change on anuran diversity of Sri Lanka would be lower compared to that of temperate regions [51, 52] as the higher latitudes have experienced greater warming than the lower latitudes [53, 54]. On the other hand, tropical species do not usually encounter seasonal changes. Hence they are mostly adapted to constant environmental conditions that are usually close to their physiological optima [55,56,57]. Such species may have limited acclimatizing capacity as they do not encounter seasonal changes [12]. Therefore, tropical ectothermic species such as anurans have limited behavioural and physiological adaptations to adjust to changing climatic factors [58]. This makes tropical species more vulnerable to even small changes of climatic factors [59]. Excessive warming during recent years, with 2016 and 2017 being reported as the warmest years on record, and the tropical region being the warmest part of the world, tropical ectotherms are highly likely to be at risk [59]. Even though Sri Lanka is an amphibian hotspot, research on the impact of climate change on physiology of amphibians is scarce. These facts emphasize the need for research on the potential physiological impacts of climate change on a tropical endemic amphibian species.

Accordingly, our overall objective was to determine the chronic physiological responses of Polypedates cruciger to climate change simulated as long-term, continuous exposure to elevated temperature and elevated CO2 (and consequently reduced pH) from early larval to adult stages. Specifically, we focused on the responses of following processes and parameters: (a) Survival and time taken for metamorphosis; (b) Growth and morphometrics of tadpoles; (c) Activity of tadpoles measured in terms of swimming performance; (d) Immunity as measured by differential white blood cell counts and lysozyme activity; (e) Antioxidant capability as indicated by catalase activity; (f) Excretory metabolism measured as ammonia released.

Material and methods

Organism used in the study

The model organism used for this study was Polypedates cruciger, (Anura: Ranidae), the Common Hourglass Tree Frog. It is an arboreal species that is usually found in forests of both wet and dry zones of Sri Lanka up to an altitude of 1525 m above mean sea level. The geographical selectivity of this species to higher altitudes makes it an excellent model for this study, as such species are considered to be more sensitive to environmental fluctuations, especially in temperature [59]. It usually spawns in a foam nest attached to a surface (e.g. tree branch, pond bank) few inches above a shallow water surface. Hatched tadpoles fall to the water below and carry out their development therein. Recent studies show that P. cruciger has expanded its habitat range, being dominant in human-modified environments [60]. Reproductive behavior of this species also makes it ideal for the current study because the shallow temporary water bodies in which it usually reproduces and tadpoles develop are highly susceptible to temperature and pH fluctuations. The current IUCN conservation status of P. cruciger is Least Concerned (LC) which makes it convenient to use as a model organism for scientific research.

Sample collection

We collected foamy egg clutches from home gardens in Kandy, Sri Lanka (Latitude 7.2906oN; Longitude 80.6337°E). They were transported to the Department of Zoology and Environmental Management, University of Kelaniya, Sri Lanka. The egg clutches were placed in 15 Litre glass aquaria containing dechlorinated tap water. The tanks were placed in a naturally-ventilated room within a building. The diurnal variation of air temperature followed a pattern similar to that of ambient temperature in the external environment. In addition to natural sunlight, artificial light was provided during the day time by fluorescent bulbs. Standard keys were used in species identification of tadpoles [61]. Once they reached the free-feeding and free-swimming stage, we pooled all tadpoles from all egg clutches and assigned fifteen tadpoles randomly to each treatment tank using a pasteur pipette. Each tank contained dechlorinated tap water (5 L) and tadpoles were acclimated to tanks for 1 week. Tadpoles were fed twice daily alternatively with formulated fish feed and chopped spinach leaves ad libitum [62].

Experimental setup

Our experimental setup included four treatments consisting of two elevated temperature treatments, one elevated CO2 treatment and a control treatment at ambient temperature and CO2. In the elevated CO2 (ECO2) treatment, we bubbled CO2 into each assigned tank until pH of water was in the range of 5.5–5.6. The pH of freshwater should be 5.5 and 5.6 respectively corresponding to the atmospheric CO2 concentrations (Ca) projected under RCP8.5 and RCP2.6 pathways. However, one pH (i.e. higher dissolved CO2) treatment in the range of 5.5–5.6 was maintained as precise maintenance of two separate pH levels in close range was not possible. Sodium citrate-citric acid biological buffer was used to maintain pH within the required range [45] and to control the solubility of CO2 in water.

We maintained our control (i.e. ambient temperature and CO2) treatment at an ambient temperature of 29 °C ± 1 °C and a pH of 7.0. In the two elevated temperature treatments, we increased water temperature by fitting Atman-100 W (Atman, China) aquarium heaters. We maintained water temperature at 32 °C (ETem32) and 34 °C (ETem34), which represented predicted temperature increases of 2 °C and 4 °C by 2100 under the IPCC’s ‘best-case’ scenario (RCP2.6) and the ‘business-as-usual’ scenario (RCP8.5) respectively. In both elevated temperature treatments, we increased water temperatures gradually at the rate of 1 °C per day up to the respective temperatures. The daily temperature fluctuation was 0.5 °C in both treatments. We triplicated each treatment. We fixed aerators bubbling ambient air containing O2 to all tanks, except those containing the elevated CO2 treatment. Despite not receiving bubbled ambient air, we expected that the elevated CO2 treatment received adequate oxygen via dissolution from ambient air. We replaced water every 48 h by siphoning and refilling with an equal volume of de-chlorinated water while maintaining all tanks at natural photoperiod.

Measurements

Morphometrics of tadpoles

We measured morphometrics of a sample of three tadpoles selected randomly from each tank once a week by placing a graph sheet underneath the tank. Length was measured by taking a freeze frame photograph and measuring it using ImageJ software (Version 1.51 k, National Institute of Health, USA). Total body length, snout-vent length, tail length, and body width were measured in tadpoles.

Activity of tadpoles

We used swimming speed to assess the activity of tadpoles using the method adopted by Jung and Jagoe [63]. We selected three tadpoles randomly from each tank for this measurement, which was performed fortnightly. We used a narrow channel (60 cm) filled with water taken from the tank in which a selected tadpole was present. Next we placed the tadpole at one end of the channel and prodded its tail gently with a pipette tip. Then we measured the time taken for the tadpole to swim the channel distance using a stopwatch.

Ammonia excretion

We measured the excretory metabolism of tadpoles in terms of their ammonia excretion to the water by the Phenate method at weekly intervals [64]. Briefly, to a water sample (25 ml), phenol solution (1 ml), sodium nitroprusside solution (1 ml) and oxidizing solution (2.5 ml) was added and after standing for 1 h for colour development, absorption was measured at 640 nm using UV and visible absorption spectrophotometer (Cecil, Great Britain). Ammonia concentration of each sample was determined using ammonium chloride as the standard.

Metamorphosis and mortality

We observed the tanks daily for dead tadpoles. Time taken for 50% of tadpoles to reach the stages of hind limb emergence (Gosner Stage 36–39) and fore limb emergence (Gosner Stage 42–46) were determined by daily observation. The experiment was terminated once 50% of tadpoles in a tank reached Gosner Stage 46.

Hematological analysis

Once 50% of tadpoles reached Gosner Stage 46, three tadpoles were randomly selected from each tank. They were euthanized using benzocaine (Sigma-Aldrich) (1.21 mM) and blood was drawn by heart puncture using a micropipette [65]. A blood smear was prepared and was stained using the Leishman-Wright’s stain [66]. Then we performed a differential white blood cell count (WBC) for each slide by determining the percentage of neutrophils, eosinophils, basophils, monocytes, lymphocytes, and thrombocytes in 50 WBC counted [65, 67,68,69]. Also we counted the ratio of WBC per 2000 erythrocytes in each slide as in Schermer [68].

Catalase enzyme activity

We measured the catalase activity in terms of the rate of H2O2 degradation by catalase [70]. We stored the liver extracted from each euthanized adult in 20 mM Tri-HCl buffer (Tris-HCl 20 mM, EDTA 1 mM, DL-dithiothreitol (DTT) 1 mM, sucrose 0.5 M, KCL 0.15 M, phenylmethylsulfonyl fluoride (PMSF) 1 mM), pH 7.4, in -80 °C. The preserved liver tissues were weighed and homogenized in 50 mM phosphate buffer (PB) (pH 7.4). Homogenate was centrifuged at 5000G in 4 °C for 15 min and the supernatant was separated. Lysate samples were diluted to 1:50 with PB and mixed with H2O2 (10 mM, 1 mL) initial absorbance at the wavelength 240 nm was measured immediately by spectrophotometer (Nanospec, Shimadzu, Japan). Decrease of absorbance was measured after 4 min. We used the formula of Cuellar-Cruz et al. to calculate the catalase activity [71]. Protein concentration of the sample was measured by Bradford assay [72]. Catalase activity was quantified relative to the protein concentration and expressed in terms of units per mg of protein.

Stomach enzyme activity

Stomachs from dissected tadpoles were preserved in phosphate-buffered saline (PBS) (pH 7.4) in -20 °C. Weight of each stomach sample was measured and was homogenized in PBS (100 μL). Using the homogenate, activity of stomach lysozyme activity was assessed using the lysoplate method [73]. Briefly, 20 μl of the homogenate was placed into wells (3.5 mm diameter and 4 mm deep) cut on nutrient agar in petri dishes with 100 mm diameter. The nutrient agar was amended with a Micrococcus luteus culture (500 μl of the bacterial culture having a cell concentration of 1 × 106 cfu/ml was added to 100 ml of nutrient agar). The diameter of the lytic zones was measured 48 h after incubation at 34 °C in nine replicates.

Statistical analyses

We used two graphical methods, namely, the normal probability plot (P-P plot) and the quantile-quantile plot (Q-Q plot) to test normality of continuous response variables. As the data points did not deviate appreciably from the fitted straight lines in the respective P-P and Q-Q plots, we proceeded with statistical analyses of all continuous variables using analysis of variance.

We determined the significance of treatment effects on the times taken to reach Gosner stages 36–39 and 42–46 by analysis of variance with the effects of elevated CO2 and temperature considered as fixed effects. Means were compared using the Duncan’s multiple range test.

We determined the time courses of the variation of morphometrics (i.e. total body length, snout-vent length, tail length and body width) by fitting growth curves using second-order polynomial functions. Growth rates of morphometric characters at Gosner stages 36–39 and 42–46 were computed as the first-derivative of the fitted polynomial functions at the time points of reaching the above stages. Initial growth rates were estimated by extrapolating time to zero. We tested the significance of treatment effects on morphometrics by repeated measures analysis of variance (RM ANOVA) using PROC MIXED of the Statistical Analysis System (SAS) [74]. A fixed effects model having a variance-covariance structure with compound symmetry was used after testing several alternative variance-covariance structures (i.e. unstructured, autoregressive and autoregressive with heterogeneous variances). We tested the significance of individual effects of elevated CO2 and elevated temperatures relative to the control (i.e. ambient CO2 and temperature) using separate contrasts within the PROC MIXED procedure. Effects of different treatments were compared among themselves using the same procedure. We compared the morphometrics at the two specific Gosner stages 36–39 and 42–46 by selecting the morphometric data at equivalent stages in different treatments and subjecting them to analysis of variance. Least squares means (LSMEANS) were used to test significance of treatment effects on morphometrics at equivalent stages. We tested the effects of elevated CO2 and temperature treatments on the activity of tadpoles by applying RM ANOVA to swimming speed data following the same procedure that was used for analysis of morphometric data.

We used categorical data analysis with PROC CATMOD in SAS to determine the significance of treatment effects on the counts of white blood cells (expressed as a percentage out of 2000 red blood cells) and on their different types (expressed as a percentage out of 50 white blood cells). Maximum likelihood estimates of frequencies in a log-linear model were used in PROC CATMOD. Significance of individual treatment comparisons was tested using a Z-statistic calculated from the respective maximum likelihood frequencies and their standard errors.

Results

Phenology

Tadpoles in all treatments reached Gosner Stages 26–30 (free feeding and free swimming stage) 2 Weeks After Hatching (2 WAH). In elevated water temperature of 34 °C (ETem34), none of the tadpoles survived until Gosner Stage 36 (hind limb emergence) and all tadpoles died 11 WAH. In the rest of the treatments (i.e. ECO2, ETem32 and the Control), durations to reach Gosner Stages 36–39 and 42–46 were significantly different among treatments (p < 0.05) (Fig. 1). Elevation of water temperature to 32 °C (ETem32) delayed development of tadpoles, with tadpoles in ETem32 taking a substantially longer time to reach the two stages than the control.

Fig. 1
figure 1

Time to reach specific Gosner stages for 50% of tadpoles experiencing elevated temperature and CO2 treatments. Control – Ambient CO2 (water pH = 7) and water temperature at 29 ± 1 °C; ETem32 – Water temperature elevated to 32 °C; ETem34 – Water temperature elevated to 34 °C. ECO2 – CO2 bubbled to water to maintain pH at 5.5–5.6. Each bar is a least squares mean of three replicate tanks each containing 15 tadpoles. Least squares means with the same letter are not significantly different at p = 0.05

Tadpoles in the elevated CO2 (ECO2) treatment reached both stages earlier than those in the control (p < 0.05).

Growth and morphometrics

Variation of total body length, snout-vent length, tail length, and body width of all treatments showed second order polynomial relationships with time (Fig. 2 and Additional file 1: Table S1). Repeated measures analysis of variance (RM ANOVA) showed highly-significant treatment effects on all morphometric characters (Table 1).

Fig. 2
figure 2

Variation of morphometric characters of tadpoles such as (a) Total body length, (b) Snout-Vent length, (c) Tail length and (d) Body width with time after hatching in different treatments. Control – Ambient CO2 (water pH = 7) and water temperature at 29 ± 1 °C; ETem32 – Water temperature elevated to 32 °C; ETem34 – Water temperature elevated to 34 °C. ECO2 – CO2 bubbled to water to maintain pH at 5.5–5.6. Each data point is a mean of three replicates. M – Metamorphosis; D – Death. Lines indicate second-order polynomial curves. Red circles indicate the Control group. Grey triangles indicate the ETem32 group. Green squares indicate the ETem34 group. Blue diamonds indicate the ECO2 group

Table 1 Significance of treatment effects and time on morphometric characters

While all morphometric characters showed highly-significant variation with time (p < 0.0001), the time x treatment interaction effect was not significant (p = 0.05). Significance tests for contrasts involving different treatment comparisons showed that there was no significant difference between ECO2 and the control for any of the morphometric characters (Table 2 and Fig. 2). In contrast, morphometrics of tadpoles in elevated temperature treatments, both combined and taken individually, were significantly different from those in the control (p < 0.01). Elevated temperatures reduced all morphometric characters relative to the control (Fig. 2). Snout-vent length and body width were lower at ETem34 in comparison to ETem32 (p < 0.01). Comparison between the effect of ECO2 and the combined effect of elevated temperatures (both ETem32 and ETem34 included in the contrast) was significant for all morphometric characters except tail length (p < 0.05). When effects of the two elevated temperatures were compared separately, comparison between ECO2 and ETem32 was significant only for total body length. On the other hand, the comparison between ECO2 and ETem34 was significant for all characters except tail length. In all instances, where the comparison between ECO2 and elevated temperatures were significant, elevated temperatures had lower morphometrics relative to ECO2 (Fig. 2).

Table 2 Significance of individual contrasts comparing different treatment combinations on morphometric characters

Comparison of morphometric characters at equivalent developmental stages (i.e. Gosner stages 36–39 and 42–46) showed highly-significant variation among experimental treatments (p < 0.001). At both stages, ETem32 did not cause significant changes in morphometric characters in comparison to the control (Figs. 3 and 4). However at Gosner stage 36–39, ECO2 reduced all morphometric characters significantly below those of the control and ETem32. At Gosner stage 42–46 also, ECO2 caused similar reductions of all morphometric characters except body width. Even though tadpoles developing in ETem34 died before reaching Gosner stage 36–39, their morphometrics 7 and 9 WAH were used to make comparisons with morphometrics of other treatments at Gosner stages 36–39 and 42–46 respectively. All morphometric characters of tadpoles in ETem34 at 7 WAH were lower than the corresponding values of tadpoles in the Control and ETem32 at Gosner stage 42–46, but were higher than those in ECO2 (Fig. 3). Total body length and snout-vent length of tadpoles in ETem34 at 9 WAH (at the time of their death) were lower than the corresponding values in the Control and ETem32 at Gosner stage 42–46 (Fig. 4). In contrast, the corresponding values of tail length and body width of tadpoles in ETem34 were not significantly different from those of the control and ETem32. On the other hand, morphometrics of tadpoles in ETem34 at 9 WAH and those in ECO2 at Gosner stage 42–46 were similar (Fig. 4). This was different from the corresponding comparison at the earlier stage (7 WAH and Gosner stage 42–46) (Fig. 3).

Fig. 3
figure 3

Effects of elevated temperature and CO2 treatments on morphometrics of tadpoles such as (a) Total body length, (b) Snout-vent length, (c) Tail length, (d) Body width at Gosner stage 36–39. Control – Ambient CO2 (water pH = 7) and water temperature at 29 ± 1 °C; ETem32 – Water temperature elevated to 32 °C; ETem34 – Water temperature elevated to 34 °C. ECO2 – CO2 bubbled to water to maintain pH at 5.5–5.6. Each bar is a least squares mean of three replicate measurements. Least squares means with the same letter are not significantly different at p = 0.05

Fig. 4
figure 4

Effects of elevated temperature and CO2 treatments on morphometrics of tadpoles such as (a) Total body length, (b) Snout-vent length, (c) Tail length, (d) Body width at Gosner stage 42–46. Control – Ambient CO2 (water pH = 7) and water temperature at 29 ± 1 °C; ETem32 – Water temperature elevated to 32 °C; ETem34 – Water temperature elevated to 34 °C. ECO2 – CO2 bubbled to water to maintain pH at 5.5–5.6. Each bar is a least squares mean of three replicate measurements. Least squares means with the same letter are not significantly different at p = 0.05

At Gosner stage 36–39, both elevated temperatures reduced growth rates of all morphometric characters in comparison to the control (Additional file 2: Table S2). A similar observation was made at Gosner stage 42–46 also, with the exception of snout-vent length and tail length in ETem32. In contrast, morphometric growth rates in ECO2 were greater than in the control at both stages. Similarly, ECO2 increased estimated initial growth rates of all morphometric traits above those of the control. In contrast, no consistent differences could be observed between initial growth rates of the elevated temperature treatments and the control.

Activity of tadpoles

The comparative variation among treatments of tadpole activity (Fig. 5) differed at different times after hatching. This was shown as a highly-significant treatment x time interaction effect in the RM ANOVA (p < 0.0001). In addition, there were significant treatment (p = 0.0126) and time (p = 0.0038) effects on swimming speed. Because of the significant treatment x time interaction, treatments were compared at each time point separately. During 2nd and 3rd WAH, tadpoles in ECO2 had faster swimming speeds than the control and the elevated temperature treatments (Table 3).

Fig. 5
figure 5

Variation of swimming speed of tadpoles in different treatments with time. Control – Ambient CO2 and water temperature at 29 ± 1 °C; ETem32 – Water temperature elevated to 32 °C; ETem34 – Water temperature elevated to 34 °C. ECO2 – CO2 bubbled to water to maintain pH at 5.5–5.6. M – Metamorphosis; D - Death

Table 3 Significance of contrasts comparing the swimming speed of tadpoles at elevated CO2 with other treatments

However, their swimming speed slowed down substantially during the 4th week as they neared metamorphosis in the 5th week. Consequently in the 4th week, activity of tadpoles in ECO2 was significantly lower than that in all other treatments. In contrast to the variation pattern shown in ECO2, activity of tadpoles in the rest of the treatments increased up to the 4th week (Fig. 5). At the early stage (i.e. 2nd week), swimming speed of tadpoles in the two elevated temperature treatments were significantly lower than in the control (Table 4).

Table 4 Significance of contrasts comparing tadpole activity in elevated temperature treatments with those in other treatments

However, with the increase of swimming speeds during the 3rd and 4th weeks in all three treatments, treatment differences were not significant during this period. After the 4th week, activity of tadpoles in ETem32 and ETem34 decreased up to the 7th week until their death in ETem34. Activity of tadpoles in the control fluctuated after the 4th week, with declines in the 5th and 8th weeks, but increases during the rest of the period up to the 10th week. At the 7th WAH, swimming speeds of tadpoles in the two elevated temperature treatments were significantly lower than in the control (Fig. 5 and Table 4). As metamorphosis of tadpoles in ETem32 was prolonged, their activity beyond the 7th week showed fluctuations. At 10 WAH, swimming rate of tadpoles in ETem32 was significantly lower than in the control. Tadpole activity decreased substantially during the week prior to metamorphosis (in ECO2 and ETem32) or death (in ETem34) in elevated temperature and CO2 treatments. Interestingly, such a decline in activity was not observed prior to metamorphosis in the control. There were slight disagreements between the two methods of significance testing (i.e. PROC GLM and PROC MIXED) in the significance of some contrasts involving elevated treatments (Table 4). Hence, results of both methods are given in Table 4.

Comparison of activity of tadpoles at equivalent developmental stages showed significant variation among treatments at Gosner stage 36–39 (p = 0.032), but no significant variation at Gosner stage 42–46. At Gosner stage 36–39, swimming speeds of tadpoles in ECO2 were significantly greater than in other treatments which did not differ significantly among themselves (Fig. 6).

Fig. 6
figure 6

Effects of elevated temperature and CO2 treatments on swimming speed of tadpoles at Gosner stages 36–39 (a) and 42–46 (b). Control – Ambient CO2 (water pH = 7) and water temperature at 29 ± 1 °C; ETem32 – Water temperature elevated to 32 °C; ETem34 – Water temperature elevated to 34 °C. ECO2 – CO2 bubbled to water to maintain pH at 5.5–5.6. Each bar is the least squares mean of three replicate measurements. Least squares means with the same letter are not significantly different at p = 0.05

Mortality of tadpoles

Cumulative mortality % (Cmort%) of tadpoles showed highly-significant variation among treatments throughout the experimental period (p < 0.0001) (Fig. 7). RM ANOVA showed highly-significant (p < 0.0001) variation in Cmort% with time, while the treatment x time interaction effect was also highly-significant. Tadpoles in the control treatment progressed through successive developmental stages to metamorphosis with 0% mortality. Across the whole experimental period, tadpoles in ETem34 showed significantly greater mortality than in all other treatments (Table 5). When the Cmort% data in different weeks were analyzed separately, a similar result was shown in all weeks except the first week (Additional file 3: Tables S3 and Additional file 4: Table S4). When the entire experimental period was considered in the RM ANOVA, tadpoles in ETem32 showed significantly greater mortality than those in the control (Table 5). Similarly, the two elevated temperature treatments, taken together, showed significantly greater mortality than the control. In contrast, Cmort% in ECO2 did not differ significantly from that in the control or ETem32, either across the first 4 weeks or in different weeks (Table 5 and Additional file 3: Table S3). However, Cmort% in ECO2 was significantly lower than that in ETem34, both overall and in different weeks.

Fig. 7
figure 7

Progression of mortality of tadpoles in different treatments. Control – Ambient CO2 and water temperature at 29 ± 1 °C; ETem32 – Water temperature elevated to 32 °C; ETem34 – Water temperature elevated to 34 °C. ECO2 – CO2 bubbled to water to maintain pH at 5.5–5.6. M – Metamorphosis; D – Death. Red circles indicate the Control group. Grey triangles indicae the ETem32 group. Green squares indicated the ETem34 group. Blue diamonds indicate the ECO2 group

Table 5 Significance of contrasts comparing the cumulative mortality percentages of tadpoles in different experimental treatments

Ammonia excretion

Tadpoles in ECO2 showed a continuous increase in ammonia excretion, measured as ammonia concentration (AmConc) in tank water, up to metamorphosis (Fig. 8). Other treatments showed substantial initial increases which were followed by decreases and further fluctuations of lesser amplitude. Across the whole experimental period, RM ANOVA showed a highly-significant (p = 0.0003) treatment x time interaction effect on AmConc. In addition, the main effects of treatments and time were also highly-significant (p < 0.0001). Significance testing of specific treatment contrasts by RM ANOVA showed that ammonia excretion of tadpoles in elevated temperatures were significantly greater than in the control (p = 0.0369) (Table 6). However, AmConc in the two elevated temperature treatments were not significantly different. When compared separately with the control, AmConc in ETemp32 was significantly greater (p = 0.0066) while that in ETemp34 was not significantly different. This was because of fluctuations in AmConc with time (Fig. 8). Therefore, the treatment differences varied at different times after hatching. For example at 2 WAH, AmConc in ETem32 was significantly greater than that in the control (p = 0.0114), but AmConc in ETem34 was lower (p < 0.0001) (Additional file 5: Table S5). However, during the following week, AmConc in the control and ETem32 decreased while that in ETem34 increased (Fig. 8) resulting in both elevated temperatures having greater AmConc than the control. Even though AmConc in the two elevated temperatures decreased at 4 WAH, they were still greater than that in the control. During the period from 5th to 8th WAH, there were no significant variations in AmConc between the elevated temperature treatments and the control. There was an increase in ammonia excretion in tadpoles in ETem32 at 9 WAH, thus raising its AmConc above those of the control and ETem34 (Fig. 8 and Additional file 5: Table S5).

Fig. 8
figure 8

Time courses of variation of ammonia concentration in tank water in different treatments. Control – Ambient CO2 and water temperature at 29 ± 1 °C; ETem32 – Water temperature elevated to 32 °C; ETem34 – Water temperature elevated to 34 °C. ECO2 – CO2 bubbled to water to maintain pH at 5.5–5.6. M – Metamorphosis; D - Death

Table 6 Significance of contrasts comparing the ammonia excretion of tadpoles in different experimental treatments

Across the 5 weeks that tadpoles in ECO2 took to metamorphose, their ammonia excretion was not significantly different from those in the control (Table 6 and Fig. 8). However, AmConc in ECO2 was significantly lower than in both elevated temperature treatments. When the data from different weeks were analyzed separately, the above differences were most prominent during the second week (Additional file 6: Table S6). Here, AmConc in ECO2 was significantly lower than that in the control as well.

Catalase enzyme activity

Catalase activity of tadpoles in ECO2 was greater than in the control (Fig. 9). However, because of the greater variability among different replicates within each treatment, the above difference was not significant at p = 0.05. Measurement of catalase activity in the elevated temperature treatments was not possible because none of the individuals survived to reach the adult stage to extract the liver.

Fig. 9
figure 9

Mean catalase activity of tadpoles in ambient and elevated CO2 treatments. U mg− 1 – Units of catalase per mg of protein in the sample. Amount of catalase needed to degrade 1 μmol of H2O2 per min. is equivalent to one Unit. Control – Ambient CO2 and water temperature at 29 ± 1 °C; ECO2 – CO2 bubbled to water to maintain pH at 5.5–5.6. Error bars which are not overlapping indicate means are not significantly different at p = 0.05

Stomach lysozyme activity

Determination of lysozyme activity of tadpoles in elevated temperature treatments was not possible because all individuals died before reaching Gosner Stage 36 in ETem34 while, in ETem32 individuals died within 24 h from reaching Gosner Stage 46. The distance of lytic zones resulting from lysozyme extracted from tadpoles of ECO2 was significantly higher than in the control (p = 0.046), thus indicating that elevated CO2 raised their lysozyme activity (Fig. 10 and Additional file 7: Plate S1).

Fig. 10
figure 10

Mean distance of lytic zones of Ambient and Elevated CO2 treatments. Ambient – Ambient CO2 and water temperature at 29 ± 1 °C; ECO2 – CO2 bubbled to water to maintain pH at 5.5–5.6. Means with the same letter are not significantly different at p = 0.05

Haematological analysis

Categorical data analysis showed that counts of white blood cells (expressed as a percentage out of 2000 red blood cells) varied significantly among treatments (p < 0.0001). Both elevated CO2 and elevated temperature at 32 °C reduced the WBC count significantly below that in the control (p < 0.0001) (Fig. 11). However, the frequency of WBCs in ETem32 did not differ significantly from that in ECO2 (p > 0.05). Haematological analysis of tadpoles in temperature elevated to 34 °C was not possible because none of the test animals survived until the Gosner Stage 46 at which blood was extracted.

Fig. 11
figure 11

Counts of white blood cells (WBC) in 2000 red blood cells (RBC) in tadpoles under different treatments. Control – Ambient CO2 (water pH = 7) and water temperature at 29 ± 1 °C; ETem32 – Water temperature elevated to 32 °C; ECO2 – CO2 bubbled to water to maintain pH at 5.5–5.6. Means of WBC counts with the same letter are not significantly different at p = 0.05

Different types of WBCs observed in the blood stains included neutrophils, monocytes, basophils, eosinophils, thrombocytes, and lymphocytes (Plate 1). Significant (p < 0.01) variation among treatments was observed in the percentages of thrombocytes, lymphocytes, monocytes, and neutrophils in WBCs (Table 7).

Plate 1
figure 12

Percentages of different types of white blood cells such as (a) Thrombocytes, (b) Lymphocytes, (c) Eosinophils, (d) Basophils, (e) Monocytes, and (d) Neutrophils in tadpoles under different treatments. Control – Ambient CO2 (water pH = 7) and water temperature at 29 ± 1 °C; ETem32 – Water temperature elevated to 32 °C; ECO2 – CO2 bubbled to water to maintain pH at 5.5–5.6. Bars with the same letter are not significantly different at p = 0.05

Table 7 Significance of treatment effects on percentages of white blood cells and its different types

In contrast, percentages of eosinophils and basophils did not differ significantly among treatments. Elevated temperature (ETem32) and CO2 (ECO2) reduced the percentage of thrombocytes (Fig. 12a) and increased the percentage of lymphocytes (Fig. 12b). However, there was no significant difference between ETem32 and ECO2 in the percentages of thrombocytes and lymphocytes. In contrast, ETem32 significantly increased the percentages of monocytes (Fig. 12e) and neutrophils (Fig. 12f), whereas ECO2 did not cause a significant change.

Fig. 12
figure 13

Blood cells observed in the liver tissue of tadpoles: (a) Thrombocytes, (b) Monocyte, (c) Neutrophil, (d) Lymphocyte, (e) Eosinophil, (f) Basophil, and (g) Erythrocytes

Deformities

In the control treatment, the only deformities observed were tail kinks and pale pigmentation while deformities such as tail kinks, oedema, beaked snout, and pale pigmentation were observed in ETem32 (Plate 2). In ECO2, the only deformity observed was the beaked snout. The highest percentage of deformities was observed in ETem32 (Table 8). Deformities could not be identified in ETem34 as all test animals died at an early stage.

Plate 2
figure 14

Tadpoles with (a) Oedema, (b) Tail kink, (c) Pale pigmentation, (d) Beaked snout and (e) and (f) Normal tadpoles

Table 8 aPercentage prevalence of deformities among test animals in different experimental treatments

Discussion

Experimental treatments viz-à-viz the thermal and climate change reality of P. cruciger

The control treatment of our experiment represents the current thermal regime experienced by P. cruciger in Sri Lanka. Temperature in the control treatment (29 °C) is typical of the lower altitude, humid tropical climate in South-Western Sri Lanka where the experiment was carried out. As Sri Lanka is located at a latitudinal range of 6-11oN, it experiences a relatively-narrow diurnal variation in ambient temperature. In this climatic zone, the mean daytime maximum and the mean night-time minimum are ca. 31o and 27 °C respectively. The day length is around 12 h throughout the year with an amplitude of only 40 min between the longest and the shortest day of the year. Therefore, 29 °C in the control treatment represents a mean temperature that tadpoles of P. cruciger would actually experience in their natural habitats in Sri Lanka. The two elevated temperatures, viz. 32o and 34 °C, are based respectively on projected temperature increases of 3o and 5 °C by mid-twenty-first Century according to the ‘best-case’ (RCP2.6) and ‘business-as-usual’ (RCP8.5) scenarios of the IPCC [1, 2]. We acknowledge that the use of constant temperatures in our experiment, rather than a diurnal variation around the mean temperatures, represents a deviation from the actual thermal regimes that tadpoles experience in their natural habitat. However, we believe that this deviation would not introduce a significant deviation in the observed response of tadpoles because of: (a) the relatively-narrow diurnal variation of temperature in this climate; and (b) the long-term decreasing trend in the diurnal temperature range due to the night-time minimum temperature increasing faster than the daytime maximum temperature [75,76,77]. Therefore, the currently narrow diurnal temperature range in Sri Lanka is likely to narrow down further in the future.

The elevated CO2 treatment in our experiment was imposed based on the projected pH reduction in the future climate change scenarios [1, 2]. Here, the absence of ambient air bubbling in the ECO2 treatment could have caused a difference in oxygen concentrations between ECO2 and the rest of the treatments. However, we believe that ECO2 was receiving adequate oxygen via natural dissolution so that any difference in oxygen concentration in tank water did not introduce a significant deviation in the organism responses.

Growth, development and survival: effects of low pH caused by elevated CO2

Our results demonstrate that reduced pH due to increased dissolved CO2 and increased water temperature cause significant changes in the rates of development and growth of Polypedates cruciger tadpoles. These provide important indicators of how future climate change may influence anuran amphibians. While consequent reduction of water pH accelerated the progression of tadpoles through successive developmental stages, elevated temperature delayed their development (Fig. 1). Our observations on the impact of elevated CO2 contradict those of previous studies, where exposure to lower pH resulted in longer larval periods [78,79,80]. This is not surprising as time to metamorphosis shows phenotypic plasticity, particularly in response to environmental stressors, to ensure the survival of young adults [81]. Furthermore, environmental stressors could influence the hormonal control of metamorphosis [82]. Accordingly, lower pH in the elevated CO2 treatment could activate the hypothalamus-pituitary-interrenal axis and accelerate metamorphosis by stimulating the secretion of thyroid and glucocorticoid hormones. Environmental stressors are known to result in fluctuations in thyroid hormone levels which affect metamorphic, developmental and morphometric traits of anuran tadpoles [83, 84]. Traits such as shorter larval period, reduced morphometrics and reduced survival which negatively affect organisms are driven by disruptions in the thyroid hormone homeostasis caused by environmental stressors. Increased thyroid hormone levels resulting from environmental stressors can lead to undersized, early metamorphosed anuran tadpoles [83, 85, 86]. This is in agreement with the shortened larval stages and lower morphometrics observed in tadpoles raised under elevated CO2 conditions in our study. Hence, disruption in thyroid hormone levels due to acid stress could be a possible explanation for this observation. Future studies on thyroid hormone levels of tadpoles of P. cruciger exposed to elevated CO2 levels would be useful for confirmation of this possibility.

Possible hypoxic conditions in the elevated CO2 treatment due to the absence of air bubbling may also have contributed to the changes in the hormonal control induced by lower pH. Furthermore, elevated CO2 and lower pH increased both developmental and growth rates, measured in terms of their morphometrics, (Fig. 2 and Additional file 2: Table S2). This contrasts with the findings of Freda and Dunson [87] and Böhmer and Rahmann who reported reduced growth rates due to acid stress at the embryonic and early larval stages [88]. One possible explanation of this disagreement is the differential responses among different amphibian species as the previous work has been done on a different frog species (i.e. Rana temporaria). On the other hand, our observation could be an acquired adaptation of P. cruciger to the acidic nature of soils in the wet zone of Sri Lanka [89], as supported by the work of Silva and Munaweera [90]. It has been shown that water in the egg collection region (i.e. Kandy) has lower pH, which is within the pH range used study. Our observation of shorter larval periods and earlier maturity could be an adaptive measure for chronic exposure to low pH levels. In nature, tadpoles mostly occupy temporary ponds that are highly vulnerable to drying followed by high acidity. In such a situation, earlier maturity allows a greater chance of survival and constitutes adaptive plasticity [91,92,93,94]. Furthermore, by shortening the larval period under a stressed condition, the risk of predation [95], desiccation [96] and infection [97] is reduced. Therefore, our observation of accelerated development while maintaining growth rates of P. cruciger tadpoles in response to chronic exposure to elevated CO2 and increased acidity may indicate an adaptive response to environmental change in their natural habitat.

It should be noted that in freshwater aquatic environments, decomposing organic matter constitutes an additional source of CO2. Therefore, in a future climate, pH of freshwater could decrease below the level predicted solely based on elevated atmospheric CO2. Thus, it is likely that freshwater organisms may be subjected to even lower levels of acidity than those in the present experiment. Therefore, adaptive responses as observed in our work could be of greater significance. Furthermore, it is possible that CO2 generated from additional sources such as decomposing organic matter could show seasonal variation depending on climatic (e.g. precipitation, temperature) and vegetation (e.g. litter fall) parameters. Accordingly, ability to adapt to fluctuating pH levels may also become an important trait for future survival of organisms which spend at least part of their life cycle in aquatic environments.

Even though reduced water pH did not affect early growth rates of tadpoles, shorter durations taken for progression through successive developmental stages (i.e. Gosner Stages 36–39 and 42–46) meant that at equivalent developmental stages, their morphometrics were lower than those in the control (Figs. 3 and 4). This is in accordance with previous studies [78, 80], which have shown that tadpoles exposed to low pH were lower in body size at metamorphosis in comparison to those at ambient pH. These metamorphic traits (i.e. size and timing of metamorphosis) are critical factors of fitness [80]. When an individual is smaller at the age of metamorphosis, it may have a lower chance of survival and reduced reproductive success [92, 94, 98]. This can be a contributory factor to the decline of amphibian populations globally as increasing CO2 reduces pH in their habitats.

Reduced growth of tadpoles in acidic conditions could be due to disrupted sodium balance [80] as even moderately acidic water disrupts the ion-regulatory process of larval amphibians [99]. Furthermore, living in a pH stressed condition requires maintenance of homeostasis while being subjected to changes of biochemical, physiological, and molecular processes [100]. This requires a greater energy expenditure, which could also result in decreased body size at metamorphosis.

Comparison of the observed mortalities of tadpoles in the present work (Fig. 7) with those in previous work is complicated by the fact that mortality and survival rates are highly dependent on species [87, 99]. Lower mortality of tadpoles exposed to acidic pH (in comparison to higher mortality levels in elevated temperature treatments) in our study agrees with observations on salamanders [101], but contrasts with 25% higher mortality of Rana temporaria tadpoles raised in acidic pH [80]. Generally, pH levels lower than 4.5 have a lethal effect on amphibian larvae while pH levels lower than 5 cause hatching and growth inhibition [102]. On the other hand, pH levels above 5 are considered sub-lethal with hatching proceeding uninhibited, but with detrimental effects on metabolism. This agrees with our observations as the pH range of current study (i.e. 5.5–5.6) falls within the sub-lethal range. In accordance with previous studies [80], mortalities in low pH treatment were not observed until commencement of metamorphosis (Fig. 7). Physiological stress that tadpoles experience with the onset of metamorphosis could have caused the observed increase in mortality with approaching metamorphosis.

Growth, development and survival: effects of elevated temperatures

Our observations demonstrate that the two elevated temperature treatments (i.e. 32 °C and 34 °C), corresponding to ‘best-case’ (RCP2.6) and the ‘worst-case’ (RCP8.5) scenarios of future climate change [1, 2], had substantial adverse impacts on the growth, development, and survival of P. cruciger tadpoles. Their growth and development show a high degree of sensitivity to the 4 – 5 °C increase in temperature in the present study. In contrast to the response to elevated CO2, elevated temperatures delayed development (Fig. 1) and reduced growth rates (Fig. 2 and Additional file 2: Table S2). Death before metamorphosis at 34 °C shows that the lethal temperature for tadpoles of P. cruciger lies between 32 °C and 34 °C. This shows that even though P. cruciger, being a species endemic and evolved in a tropical climate with a relatively high temperature regime, is adapted to inhabit a higher temperature (e.g. 29 °C in the present situation) than temperate species such as Rana temporaria and Bufo bufo [103], it has a narrow thermal tolerance range. This is in accordance with the postulation of Janzen [56] that organisms that inhabit the relatively less variable thermal regime of tropical climates are acclimated and evolutionary adapted to a narrower fluctuation in their environment than comparable organisms inhabiting the more variable temperate climates. Evidence supporting the applicability of Janzen’s postulation to anuran amphibians has been shown by Ruthsatz et al. [83] who found that larvae of the African clawed frog (Xenopus laevis) reared at higher temperatures had a lower range of thermal tolerance that those reared at lower temperatures, despite having a higher maximum thermal limit. Similarly, Drakulic et al. [104] observed that morphometrics, physiological condition and activity of metamorphs of Rana temporaria populations originating from warmer habitats show adaptation to warmer temperatures than metamorphs originating from cooler habitats. Adaptation to the local environment and its fluctuations is particularly important for tadpoles of anuran because of their limited capability for thermoregulation and movement towards favourable habitats [105]. However, data on local adaptations and tolerance limits of P. cruciger and/or similar local species are lacking. Hence, we suggest this as a future direction of study. Increasing mortality rate of tadpoles exposed to 34 °C from the beginning of the experiment indicated their inability to acclimate to the increased temperature (Fig. 7). The time at which 100% mortality was observed in this treatment was approximately synchronous with the time at which tadpoles in the control treatment initiated metamorphosis. Physiological stress caused with the onset of metamorphosis, exacerbated by the elevated temperature, could also have contributed to tadpole mortality before metamorphosis at 34 °C. This indicates that tadpoles of P. cruciger have little chance of acclimation and survival in the predicted future temperature increases under the ‘business-as-usual’ scenario (RCP 8.5).

On the other hand, tadpoles exposed to elevated temperature at 32 °C showed only 10% mortality during the first 2 weeks (Fig. 7). It remained constant until commencement of metamorphosis indicating the ability of tadpoles of P. cruciger to tolerate to 32 °C. However, this was possible only for the larval stage as mortality was observed in all late-stage larvae (Gosner Stage 42–46) raised at 32 °C within 24 h of metamorphosis. This is in agreement with a similar study on a montane frog species Eleutherodactylus portoricensis which rarely experienced temperatures above 30 °C [106]. Similarly, P. cruciger rarely experiences temperatures as high as 32 °C because of the altitude (300–1525 m) of its natural habitat. Hence the late-stage larvae of P. cruciger may not have developed sufficient adaptations to tolerate dehydration and other physiological changes caused by temperatures above 29 °C. This could be the reason for its mortality at the exposure to elevated temperature. However, survival of tadpoles at 32 °C up to the stage of metamorphosis indicates the possibility of this species developing adaptations for survival under predicted future temperatures in the best-case scenario (RCP2.6). This is supported by the observation of Ruthsatz et al. [83] that tadpoles developed at warmer temperatures show higher maximum thermal limits. Drakulic et al. [105] also provide evidence that anuran amphibians have the capacity to adapt to warmer temperatures after they have experienced a higher temperature for some period.

Growth rates of tadpoles exposed to 32 °C and 34 °C decreased relative to that of tadpoles in ambient temperature from the 4th week onwards (Fig. 2). This is because growth and development would proceed slowly when temperature is supra-optimal relative to the range required by the cellular processes [107]. Tadpoles raised at 32 °C have taken a longer time than tadpoles at ambient temperature to reach the stage of metamorphosis, which is in accordance with similar studies [107]. Reduced morphometrics and delayed metamorphosis could be due to changes in energy allocation with most of it being diverted to maintain a high routine metabolic rate. Furthermore, reduced morphometric growth leads to delays in attaining the minimum size threshold required for metamorphosis [83, 108,109,110]. Reduced morphometrics and longer larval period of tadpoles exposed to 32 °C makes P. cruciger mostly vulnerable to predation, infection and desiccation in future temperature predicted even under the best-case-scenario (RCP 2.6).

Negative traits such as the reduced morphometrics and reduced survival observed in both elevated temperature and elevated CO2 conditions could have yielded from changes of the environmental stressor-induced thyroid hormone levels [83, 84]. Existing studies show that thyroid hormone levels altered due to environmental stressors have resulted in substantially low survival rates [83]. Therefore, we recommend that future studies are carried out on P. cruciger to confirm how thyroid hormone homeostasis influence reduced survival and lower morphometrics observed under elevated CO2 and temperature.

Activity of tadpoles

Elevated CO2 and temperature caused changes in the activity of tadpoles. Elevated CO2 increased swimming speed of tadpoles during their early growth stages (i.e. Gosner stages 26–30), elevated temperatures reduced tadpole activity than that of the control (Fig. 5). This observation agrees with previous studies which demonstrate that high temperatures impair motility of tadpoles [107]. Therefore, future temperature increases could increase vulnerability of P. cruciger tadpoles for predation due to reduced motility. We suggest further studies with predatory pressure incorporated to investigate this aspect. The marked reduction of swimming speed that was observed in the week prior to metamorphosis in elevated CO2 and elevated temperature treatments could be due to reduced metabolism and activity during metamorphosis. However, stress induced by reduced pH and increased temperature could also have contributed to this reduction as it was not observed in the control treatment.

Catalase enzyme activity

Increased catalase enzyme activity in tadpoles raised in low pH in our work is in agreement with similar studies [111]. Amphibians live in small, ephemeral aquatic environments are highly susceptible to large fluctuations of pH and temperature, where the resulting biochemical, physiological and molecular changes require maintenance of homeostasis and a higher rate of metabolism. This increases respiration and increases production of reactive oxygen species (ROS) as a byproduct, leading to oxidative stress. This could increase the activity of enzymes such as catalase which is a key component of antioxidant defence systems [100, 112].

Ammonia excretion

Concentration of released ammonia is an indication of excretory metabolism of tadpoles. In our study, elevated CO2 decreased ammonia excretion of tadpoles (Fig. 8, Table 6 and Additional file 6: Table S6), which contradicts observations on anuran tadpoles and fish exposed to acidity [99, 113,114,115]. One explanation for this discrepancy may be the differences between species and the rates of exposure to low pH. In an acidic medium, more ammonia would be protonated upon excretion as a result of increased external H+ ion concentration. It could affect the overall ammonia excretion [116,117,118]. The increase in ammonia gradient across gills and skin would facilitate greater ammonia excretion. Although acute exposure may give this result, as animals in this study were chronically exposed to low pH, a new equilibrium is likely to establish at lower blood ammonia, resulting in reduced concentrations of excreted ammonia [99]. Furthermore, moderately acidic water disrupts the ion-regulatory process of larval amphibians. Short- and long-term exposure to low pH resulted in amphibians losing 21–62% of body sodium resulting in diffusive loss of ions [87]. This was mostly due to stimulation of sodium efflux. Therefore, a part of the increase in ammonia excretion in previous studies may be due to passive ammonia efflux across a leaky gill or membrane. Decrease of ammonia excretion in our study could be due to the sodium present in the citrate buffer in the medium preventing the increased sodium efflux resulting from low pH. Nevertheless, reduced ammonia excretion in the elevated CO2 (i.e. lower pH) treatment merits further investigation.

In contrast to elevated CO2, elevated temperatures increased ammonia excretion in tadpoles (Fig. 8, Table 6 and Additional file 5: Table S5). This could be due to heat stress induced faster metabolism and consequently higher respiration rates. Protein catabolism is increased with increased respiration and thereby increases the excretion of ammonia nitrogen. Weekly variation in ammonia excretion showed reductions during the weeks that a majority of tadpoles were undergoing metamorphosis (e.g. week 7 onwards in the control, week 9 onwards in elevated 32 °C). This was probably because of reduced metabolism and activity during metamorphosis. Reduced ammonia excretion as metamorphosis approached could have been due to tadpoles changing their excretory metabolism from ammonia to urea. However, we do not have information on whether this transition occurred and if so when in this experiment.

Immunity

Research on the impact of elevated temperature and low pH on the immune response of tadpoles is limited, with most reporting the impact of reduced temperatures (i.e. winter, hibernation) on the immunity of amphibians. Reduction in the white blood cell (WBC) count in the elevated CO2 treatment (Fig. 11) is in agreement with previous studies where a reduction of splenic WBC was observed in Rana pipens tadpoles exposed to pH 5.5 [45]. Increased lymphocyte (Fig. 12b), monocyte (Fig. 12e) and neutrophil (Fig. 12f) counts in tadpoles raised in 32 °C suggest an infection as they are the major WBCs involved in phagocytosis. This is supported by the observed deformities which could possibly be related to trematode infections and increased mortality of adults [62, 119]. Reduction of thrombocytes in tadpoles experiencing elevated temperature and reduced pH (Fig. 12a) indicates reduced immunity as thrombocytes play an important role in hemostasis by plugging damaged blood vessels.

Significantly larger lysis zone in the low pH treatment (Fig. 10) indicated increased lysozyme activity in tadpoles compared to those in ambient pH. This may suggest an infection in tadpoles in this treatment [44]. These observations suggest that environmental fluctuations such as elevation of temperature and pH make tadpoles of P. cruciger more vulnerable to infections, even under the best-case scenario of climate change.

Deformities

Developmental and morphological deformities such as oedema, beaked snout, pale pigmentation, and tail kink were observed (Plate 2). Out of these, only two were observed in tadpoles raised under ambient conditions and those were also in very low percentages. Beaked snout was observed in tadpoles raised in low pH while all deformities were observed in high percentages in tadpoles raised in 32 °C. These deformities usually result from filling of fluid in the body cavity (oedema), pale-colored liver, congestive and hemorrhagic organs, malformations in melanin producing cells, and malformations of the spinal cord [120]. Also, most of these deformities were similar to those induced by trematode infections [62, 119], which indicates an infection in tadpoles exposed to 32 °C, where most deformities were observed. Deformities could also be caused by disruption of thyroid hormone system [120], which is highly-sensitive to environmental stressors such as higher temperatures [83] and toxicants in the environment such as air pollutants and agro-chemicals [85].

These malformations hamper movements of tadpoles and make them vulnerable [119], such as pale pigmentation makes them more vulnerable to predation. This is applicable in the present study because tadpoles of P. cruciger normally develop in temporary, shallow ponds formed on places such as rock crevices and garden ponds which are usually dark in background. Therefore, having a pale pigmentation in such a background can make the tadpoles more conspicuous for predators, thus increasing their vulnerability. On the other hand, motor disorders (tail kink and oedema) and mouth deformities (beaked mouth) impair swimming, (by loss of balance) foraging and feeding [48]. This could have contributed to the reduced growth rate and eventually the reduced survival in tadpoles raised in 32 °C. Thus, it is clear that elevated temperatures and low pH levels, even under the best-case scenario, induce deformities in P. cruciger tadpoles, which affect them negatively by increasing their susceptibility to predation and reducing foraging success in natural habitats. This observation carries significant weight, although the exact reason behind increased occurrence of deformities under elevated temperature and low pH levels cannot be pinpointed. That is because, to the best of our knowledge, this is the first time that occurrence of deformities in P. cruciger has been discussed in the context of climate change. Existing literature has discussed how toxicity of pesticides induces deformities in this species [62, 85, 119]. However, the influence of increased acidity and elevated temperature on occurrence of deformities has not been discussed with respect to this species. Furthermore, there exists a remarkable paucity of research discussing the relationship between elevated temperatures and deformities in anurans, although studies exist on how these elevated acidic conditions induce deformities in other anuran species [121, 122]. Therefore, species-specific research, designed to identify the exact causes behind increased prevalence of deformities in tadpoles reared under elevated temperature and acidic conditions is recommended.

Future directions

To the best of our knowledge, this is the first time that impact of climate change on this species has been assessed, although the impacts of agrochemical toxicity on P. cruciger have been investigated. That makes the findings of this study significant and novel, opening several paths for future research based on detailed seasonal field studies, assessing the impact of climate change on this species. An endemic species such as P. cruciger usually experience constant conditions close to its physiological optima, since Sri Lanka is a tropical country not experiencing annual seasonal variations [55,56,57]. Therefore it is likely for P. cruciger to have evolved very limited adaptations to tolerate fluctuating environmental conditions. That makes this species more vulnerable to climate change [58, 59]. Therefore, assessing genetic aspects determining the tolerance capacity and physiological responses of P. cruciger to future climatic conditions and empirical studies based on computational modeling predicting the fate of this species in future climatic conditions are highly recommended as future research avenues.

Furthermore, P. cruciger is a species prevalent in both ‘wet’ (i.e. humid tropical) and ‘dry’ (i.e. sub-humid tropical) zones of Sri Lanka, across an altitude ranging up to 1525 m above sea level, in a wide and patchy distribution pattern. Physiological and behavioral responses of anurans to environmental stressors such as thermal stress are highly population-specific [104, 105, 123]. That is because populations located in different habitats are likely to have developed variable local adaptations resulting in a discrepancy in the ways that they respond to environmental stressors. Hence, there is a possibility that different populations of P. cruciger distributed in varying climatic zones and altitudes of the island to show different responses to the factors tested in this study. Findings of this study are based on test animals that were collected only from the wet zone. We recommend future studies to be conducted on test animals collected from a wider environmental range to obtain a more generalized insight of how P. cruciger would respond to climate change.

Conclusions

This work provides a comprehensive analysis of the physiological response of Polypedates cruciger (Common hourglass tree frog) to two key aspects of climate change; increased temperature and increased CO2-induced reduced pH in aquatic habitats. Based on the observed results, we conclude that chronic exposure to elevated temperature and CO2 alters many physiological responses of tadpoles of P. cruciger which could increase their susceptibility to predation, infection, desiccation, and mortality. Elevated temperature and low pH-induced effects incorporate reduced morphometrics and body size, increased motility changes in immune cells and deformities. Hence, we conclude that climate change can possibly impose a significant threat to P. cruciger. However, it may be possible for the species to adapt to low pH conditions by accelerating development and shortening the larval period. Increased activity of antioxidant enzymes such as catalase and an immune response via increased production of white blood cells could constitute further adaptive responses to physiological stress induced by altered temperatures and pH. This preliminary study is the first of its kind to report the physiological responses of an endemic amphibian species of Sri Lanka to future climate change. We acknowledge that the constant temperatures used in our experimental tanks deviate from the fluctuating temperatures in a natural habitat. However, in view of the narrow amplitude of the diurnal variation of temperature in tropical climates, our findings provide important insights into the response of tropical anuran amphibians to future climate change.

Availability of data and materials

The data sets generated in the current study are available from the corresponding author on reasonable request.

Abbreviations

AmConc:

Ammonia concentration in tank water

Ca :

Atmospheric CO2 concentration

Cmort%:

Cumulative percentage of mortality

ECO2:

Elevated CO2 treatment

ETem32:

Treatment in which tank water temperature was elevated to 32 °C

ETem34:

Water temperature elevated to 34 °C

IPCC:

Intergovernmental Panel on Climate Change

RCP:

Representative Concentration Pathway

ROS:

Reactive oxygen species

WAH:

Weeks after hatching

WBC:

White blood cells

References

  1. Stocker TF, Qin D, Plattner G-K, Alexander LV, Allen SK, Bindoff NL, Bréon F-M, Church JA, Cubasch U, Emori S, et al. In: Stocker TF, Qin D, Plattner G-K, Tignor M, Allen SK, Boschung J, Nauels A, Xia Y, Bex V, Midgley PM, editors. editors Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the IPCC. Cambridge: Cambridge University Press; 2013. p. 33–115.

    Google Scholar 

  2. van Vuuren DP, Edmonds J, Kainuma M, Riahi K, Thomson A, Hibbard K, Hurtt GC, Kram T, Krey V, Lamarque JF, Masui T. The representative concentration pathways: an overview. Clim Chang. 2011;109:5–31.

    Article  Google Scholar 

  3. Houghton JT. Global warming: the complete briefing. 4th ed. Cambridge: Cambridge University Press; 2009.

    Book  Google Scholar 

  4. Phillips JC, McKinley GA, Bennington V, Bootsma HA, Pilcher DJ, Sterner RW, Urban NR. The potential for CO2-induced acidification in freshwater: a Great Lakes case study. Oceanography. 2015;28:136–45.

    Article  Google Scholar 

  5. Solomon CT, Jones SE, Weidel BC, Buffam I, Fork ML, Karlsson J, Larsen S, Lennon JT, Read JS, Sadro S, Saros JE. Ecosystem consequences of changing inputs of terrestrial dissolved organic matter to lakes: current knowledge and future challenges. Ecosystems. 2015;18:376–89.

    Article  Google Scholar 

  6. Dunson WA, Travis J. The role of abiotic factors in community organization. Am Nat. 1991;138:1067–91.

    Article  Google Scholar 

  7. Hochochka PW, Somero GN. Biochemical adaptation: mechanism and process in biochemical evolution. New York: Oxford University Press; 2002.

    Google Scholar 

  8. Case TJ, Holt RD, McPeek MA, Keitt TH. The community context of species’ borders: ecological and evolutionary perspectives. Oikos. 2005;108:28–46.

    Article  Google Scholar 

  9. Pearson RG, Dawson TP. Predicting the impacts of climate change on the distribution of species: are bioclimate envelope models useful? Glob Ecol Biogeogr. 2003;12:361–71.

    Article  Google Scholar 

  10. Southward AJ, Hawkins SJ, Burrows MT. Seventy years’ observations of changes in distribution and abundance of zooplankton and intertidal organisms in the western English Channel in relation to rising sea temperature. J Therm Biol. 1995;20:127–55.

    Article  Google Scholar 

  11. Gilman EL, Ellison J, Jungblut V, Van Lavieren H, Wilson L, Areki F, Brighouse G, Bungitak J, Dus E, Henry M, Kilman M. Adapting to Pacific Island mangrove responses to sea level rise and climate change. Clim Res. 2006;32:161–76.

    Article  Google Scholar 

  12. Stillman JH. Acclimation capacity underlies susceptibility to climate change. Science. 2003;301(5629):65.

    Article  CAS  PubMed  Google Scholar 

  13. Hoffmann AA, Sørensen JG, Loeschcke V. Adaptation of Drosophila to temperature extremes: bringing together quantitative and molecular approaches. J Therm Biol. 2003;28:175–216.

    Article  Google Scholar 

  14. Johns DM, Howell WH, Klein-MacPhee G. Yolk utilization and growth to yolk-sac absorption in summer flounder (Paralichthys dentatus) larvae at constant and cyclic temperatures. Mar Biol. 1981;63:301–8.

    Article  Google Scholar 

  15. Parson EA, Fisher-Vanden AK. Integrated assessment models of global climate change. Annu Rev Energ Env. 1997;22:589–628.

    Article  Google Scholar 

  16. Stuart SN, Chanson JS, Cox NA, Young BE, Rodrigues AS, Fischman DL, Waller RW. Status and trends of amphibian declines and extinctions worldwide. Science. 2004;306(5702):1783–6.

    Article  CAS  PubMed  Google Scholar 

  17. Wells KD. The ecology and behavior of amphibians. Chicago: University of Chicago Press; 2010.

    Google Scholar 

  18. Dadour IR, Cook DF, Wirth N. Rates of development of Hydrotaea rostrata under summer and winter (cyclic and constant) temperature regimes. Med Vet Entomol. 2001;15:177–82.

    Article  CAS  PubMed  Google Scholar 

  19. Niehaus AC, Wilson RS, Franklin CE. Short-and long-term consequences of thermal variation in the larval environment of anurans. J Anim Ecol. 2006;75:686–92.

    Article  PubMed  Google Scholar 

  20. Angilletta MJ Jr, Niewiarowski PH, Navas CA. The evolution of thermal physiology in ectotherms. J Therm Biol. 2002;27:249–68.

    Article  Google Scholar 

  21. Jiang L, Morin PJ. Temperature-dependent interactions explain unexpected responses to environmental warming in communities of competitors. J Anim Ecol. 2004;73:569–76.

    Article  Google Scholar 

  22. Katzenberger MJ. Thermal tolerance and sensitivity of amphibian larvae from Palearctic and Neotropical communities, vol. 2009: Univesidade de Lisboa; PhD dissertation., 2009.

  23. Christian KA, Tracy CR, Porter WP. The effect of cold exposure during incubation of Sceloporus undulatus eggs. Copeia. 1986;1986:1012–4.

    Article  Google Scholar 

  24. Tingle CC, Copland MJ. Predicting development of the mealybug parasitoids Anagyrus pseudococci, Leptomastix dactylopii and Leptomastidae abnormis under glasshouse conditions. Entomol Exp Appl. 1988;46:19–28.

    Article  Google Scholar 

  25. Beck SD. Insect thermoperiodism. Annu Rev Entomol. 1983;28:91–108.

    Article  Google Scholar 

  26. Foden WB, Butchart SHM, Stuart SN, Vié J-C, Akçakaya HR, Angulo A, et al. Identifying the world’s most climate change vulnerable species: a systematic trait-based assessment of all birds, amphibians and corals. PLoS One. 2013;8(6):e65427. https://doi.org/10.1371/journal.pone.0065427.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Bickford D, Howard SD, Ng DJ, Sheridan JA. Impacts of climate change on the amphibians and reptiles of Southeast Asia. Biodivers Conserv. 2010;19:1043–62.

    Article  Google Scholar 

  28. Ahmad I, Maria VL, Oliveira M, Serafim A, Bebianno MJ, Pacheco M, Santos MA. DNA damage and lipid peroxidation vs. protection responses in the gill of Dicentrarchus labrax L. from a contaminated coastal lagoon (Ria de Aveiro, Portugal). Sci Total Environ. 2008;406:298–307.

    Article  CAS  PubMed  Google Scholar 

  29. Halliwell B, Gutteridge JM. The importance of free radicals and catalytic metal ions in human diseases. Mol Asp Med. 1985;8:89–193.

    Article  CAS  Google Scholar 

  30. Chelikani P, Fita I, Loewen PC. Diversity of structures and properties among catalases. Cell Mol Life Sci. 2004;61:192–208.

    Article  CAS  PubMed  Google Scholar 

  31. Toner K, Sojka G, Ellis R. A quantitative enzyme study; CATALASE. bucknell.edu. Retrieved. 2007;2007:02–11 http://web.archive.org/web.

    Google Scholar 

  32. Schomburg, D., 2015. BRENDA-the Comprehensive enzyme information system. http://www.Brenda-enzymes.org. Accessed on 19 Feb 2018.

    Google Scholar 

  33. Blaustein AR, Romansic JM, Kiesecker JM, Hatch AC. Ultraviolet radiation, toxic chemicals and amphibian population declines. Divers Distrib. 2003;9:123–40.

    Article  Google Scholar 

  34. Kiesecker JM. Synergism between trematode infection and pesticide exposure: a link to amphibian limb deformities in nature? Proc Natl Acad Sci U S A. 2002;99:9900–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Daszak P, Berger L, Cunningham AA, Hyatt AD, Green DE, Speare R. Emerging infectious diseases and amphibian population declines. Emerg Infect Dis. 1999;5:735–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Daszak P, Cunningham AA, Hyatt AD. Infectious disease and amphibian population declines. Divers Distrib. 2003;9:141–50.

    Article  Google Scholar 

  37. Kiesecker JM, Belden LK, Shea K, Rubbo MJ. Amphibian decline and emerging disease: what can sick frogs teach us about new and resurgent diseases in human populations and other species of wildlife? Am Sci. 2004;92:138–47.

    Article  Google Scholar 

  38. Lips KR, Brem F, Brenes R, Reeve JD, Alford RA, Voyles J, Carey C, Livo L, Pessier AP, Collins JP. Emerging infectious disease and the loss of biodiversity in a Neotropical amphibian community. Proc Natl Acad Sci U S A. 2006;103:3165–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Pounds JA, Bustamante MR, Coloma LA, Consuegra JA, Fogden MP, Foster PN, La Marca E, Masters KL, Merino-Viteri A, Puschendorf R, Ron SR. Widespread amphibian extinctions from epidemic disease driven by global warming. Nature. 2006;439(7073):161–7.

    Article  CAS  PubMed  Google Scholar 

  40. Janeway CA, Travers P, Walport M, Shlomchik M. Immunobiology: the immune system in health and disease. London: 5th ed. New York: Garland Science; 1996.

    Google Scholar 

  41. Sahoo NR, Kumar P, Bhusan B, Bhattacharya TK, Dayal S, Sahoo M. Lysozyme in livestock: a guide to selection for disease resistance: a review. J Anim Sci Adv. 2012;2:347–60.

    Google Scholar 

  42. Jollès P, Jollès J. What's new in lysozyme research? Mol Cell Biochem. 1984;63:165–89.

    Article  PubMed  Google Scholar 

  43. Miyazaki T. Influences of pH and temperature on lysozyme activity in the plasma of Japanese flounder and Japanese char. Fish Pathol. 1998;33:7–10.

    Article  Google Scholar 

  44. Raffel TR, Rohr JR, Kiesecker JM, Hudson PJ. Negative effects of changing temperature on amphibian immunity under field conditions. Funct Ecol. 2006;20:819–28.

    Article  Google Scholar 

  45. Brodkin M, Vatnick I, Simon M, Hopey H, Butler-Holston K, Leonard M. Effects of acid stress in adult Rana pipiens. J Exp Zool Part A: Comp Exp Biol. 2003;298:16–22.

    Article  Google Scholar 

  46. Bibi R, Christi K. Effect of carbonate ions (pH) on the growth rate of Rhinella marina tadpoles. Int J Adv Res Rev. 2017;2(3):1–7.

    Google Scholar 

  47. Cooke AS. Tadpoles as indicators of harmful levels of pollution in the field. Environ Pollut A, Ecol Biol. 1981;25:123–33.

    Article  CAS  Google Scholar 

  48. Brunelli E, Bernabò I, Berg C, Lundstedt-Enkel K, Bonacci A, Tripepi S. Environmentally relevant concentrations of endosulfan impair development, metamorphosis and behaviour in Bufo bufo tadpoles. Aquat Toxicol. 2009;91:135–42.

    Article  CAS  PubMed  Google Scholar 

  49. Burger J, Snodgrass J. Heavy metals in bullfrog (Rana catesbeiana) tadpoles: effects of depuration before analysis. Environ Toxicol Chem. 1998;17:2203–9.

    Article  CAS  Google Scholar 

  50. Manamendra-Arachchi K, Pethiyagoda R. Amphibians of Sri Lanka, vol. 2009. Colombo: Wildlife Heritage Trust; 2006.

    Google Scholar 

  51. Parmesan C. Influences of species, latitudes and methodologies on estimates of phenological response to global warming. Glob Chang Biol. 2007;13:1860–72.

    Article  Google Scholar 

  52. Root TL, Price JT, Hall KR, Schneider SH, Rosenzweig C, Pounds JA. Fingerprints of global warming on wild animals and plants. Nature. 2003;421(6918):57–60.

    Article  CAS  PubMed  Google Scholar 

  53. IPCC. Climate change 2007: synthesis report. Contribution of working groups I, II and III to the fourth assessment report of the intergovernmental panel on climate change. Core writing team, Pachauri RK, Reisinger a, editors. Geneva: IPCC; 2007. p. 104.

    Google Scholar 

  54. Hartmann DL, Klein Tank AMG, Rusticucci M, Alexander LV, Bronnimann S, et al. Observations: Atmosphere and Surface. In: Stocker TF, Qin D, Plattner G-K, Tignor M, Allen SK, Boschung J, Nauels A, Xia Y, Bex V, Midgley PM, editors. editors Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press; 2013. p. 159–254.

    Google Scholar 

  55. Deutsch CA, Tewksbury JJ, Huey RB, Sheldon KS, Ghalambor CK, Haak DC, Martin PR. Impacts of climate warming on terrestrial ectotherms across latitude. P Natl Acad Sci USA. 2008;105:6668–72.

    Article  CAS  Google Scholar 

  56. Janzen DH. Why mountain passes are higher in the tropics. Am Nat. 1967;101:233–49.

    Article  Google Scholar 

  57. Perez TM, Stroud JT, Feeley KJ. Thermal trouble in the tropics. Science. 2016;351(6280):1392–3.

    Article  CAS  PubMed  Google Scholar 

  58. Ruibal R. Thermal relations of five species of tropical lizards. Evolution. 1961;15:98–111.

    Article  Google Scholar 

  59. Tewksbury JJ, Huey RB, Deutsch CA. Putting the heat on tropical animals. Science. 2008;320(5881):1296–7.

    Article  CAS  PubMed  Google Scholar 

  60. Pethiyagoda R, Manamendra-Arachchi K. Evaluating Sri Lanka’s amphibian diversity. Occas Pap Wildl Herit Trust. 1998;2:1–2.

    Google Scholar 

  61. Kirthisinghe P. The Amphibia of Ceylon, 1st ed. Charles Circus. Colombo: Charles Circus; 1955.

    Google Scholar 

  62. Jayawardena UA, Rohr JR, Navaratne AN, Amerasinghe PH, Rajakaruna RS. Combined effects of pesticides and trematode infections on hourglass tree frog Polypedates cruciger. Ecohealth. 2016;13:111–22.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Jung RE, Jagoe CH. Effects of low pH and aluminum on body size, swimming performance, and susceptibility to predation of green tree frog (Hyla cinerea) tadpoles. Can J Zool. 1995;73:2171–83.

    Article  CAS  Google Scholar 

  64. Verdouw H, Van Echteld CJ, Dekkers EM. Ammonia determination based on indophenol formation with sodium salicylate. Water Res. 1978;12:399–402.

    Article  CAS  Google Scholar 

  65. Das M, Mahapatra PK. Hematology of wild caught dubois’s tree frog Polypedates teraiensis, Dubois, 1986 (Anura: Rhacophoridae). Sci World J. 2014. https://doi.org/10.1155/2014/491415.

    CAS  Google Scholar 

  66. Gajendra S, Jha B, Goel S, Sahni T, Sharma R, Shariq M, Jaiswal S, Sachdev R. Leishman and Giemsa stain: a new reliable staining technique for blood/bone marrow smears. Int J of Lab Hematol. 2015;37:774–82.

    Article  CAS  Google Scholar 

  67. Hadji-Azimi I. Atlas of adult Xenopus laevis laevis hematology. Dev Comp Immunol. 1987;11:807–74.

    Article  CAS  PubMed  Google Scholar 

  68. Schermer S. Blood morphology of laboratory animals. 3rd ed. Philadelphia: F.A. Davis Co.; 1987.

    Google Scholar 

  69. Ussing AP, Rosenkilde P. Effect of induced metamorphosis on the immune system of the axolotl, Ambystoma mexicanum. Gen Comp Endoc. 1995;97:308–19.

    Article  CAS  Google Scholar 

  70. Aebi H. 1984. [13] Catalase in vitro. Methods Enzymol. 1984;105:121–6.

    Article  CAS  PubMed  Google Scholar 

  71. Cuéllar-Cruz M, Castaño I, Arroyo-Helguera O, De Las Peñas A. Oxidative stress response to menadione and cumene hydroperoxide in the opportunistic fungal pathogen Candida glabrata. Memorias do Instituto Oswaldo Cruz. 2009;104:649–54.

    Article  PubMed  Google Scholar 

  72. Bradford MM. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976;72:248–54.

    Article  CAS  PubMed  Google Scholar 

  73. Yousif AN, Albright LJ, Evelyn TP. Occurrence of lysozyme in the eggs of coho salmon Oncorhynchus kisutch. Dis Aquat Org. 1991;10:45–9.

    Article  Google Scholar 

  74. Littell RC, Henry PR, Ammerman CB. Statistical analysis of repeated measures data using SAS procedures. J Anim Sci. 1998;76:1216–31.

    Article  CAS  PubMed  Google Scholar 

  75. IPCC. In: Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, Tignor M, Miller HL, editors. editors Climate change 2007: the physical science basis. Contribution of working group I to the fourth assessment report of the intergovernmental panel on climate change. Cambridge: Cambridge University Press, Cambridge; 2007. p. 996.

    Google Scholar 

  76. Vose RS, Easterling DR, Gleason B. Maximum and minimum temperature trends for the globe: An update through 2004. Geophys Res Lett. 2005;32:L23822. https://doi.org/10.1029/2005GL024379.

    Article  Google Scholar 

  77. Donat MG, Alexander LV, Yang H, Durre I, Vose R, Dunn RJH, Willett KM, Aguilar E, Brunet M, Caesar J, Hewitson B. Updated analyses of temperature and precipitation extreme indices since the beginning of the twentieth century: the HadEX2 dataset. J Geophys Res-Atmos. 2013;118:2098–118.

    Article  Google Scholar 

  78. Cummins CP. Effects of aluminium and low pH on growth and development in Rana temporaria tadpoles. Oecologia. 1986;69:248–52.

    Article  CAS  PubMed  Google Scholar 

  79. Grant KP, Licht LE. Acid tolerance of anuran embryos and larvae from Central Ontario. J Herpetol. 1993;27:1–6.

    Article  Google Scholar 

  80. Räsänen K. Evolutionary implications of acidification: a frog’s eye view. Doctoral dissertation. Uppasala: Acta Universitatis Upsaliensis; 2002.

    Google Scholar 

  81. Denver RJ. Environmental stress as a developmental cue: corticotropin-releasing hormone is a proximate mediator of adaptive phenotypic plasticity in amphibian metamorphosis. Horm Behav. 1997;31:169–79.

    Article  CAS  PubMed  Google Scholar 

  82. Kikuyama S, Kawamura K, Tanaka S, Yamamoto K. Aspects of amphibian metamorphosis: hormonal control. Int Rev Cytol. 1993;145:105–48.

    Article  CAS  PubMed  Google Scholar 

  83. Ruthsatz K, Dausmann KH, Peck MA, Drees C, Sabatino NM, Becker LI, Reese J, Hartmann L, Glos J. Thyroid hormone levels and temperature during development alter thermal tolerance and energetics of Xenopus laevis larvae. Conserv Physiol. 2018;6:coy059. https://doi.org/10.1093/conphys/coy059.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Sowers AD, Mills MA, Klaine SJ. The developmental effects of a municipal wastewater effluent on the northern leopard frog, Rana pipiens. Aquat Toxicol. 2009;94:145–52.

    Article  CAS  PubMed  Google Scholar 

  85. Mann RM, Hyne RV, Choung CB, Wilson SP. Amphibians and agricultural chemicals: review of the risks in a complex environment. Environ Pollut. 2009;157:2903–27.

    Article  CAS  PubMed  Google Scholar 

  86. Cauble K, Wagner RS. Sublethal effects of the herbicide glyphosate on amphibian metamorphosis and development. B Environ Contam Tox. 2005;75:429–35.

    Article  CAS  Google Scholar 

  87. Freda J, Dunson WA. Field and laboratory studies of ion balance and growth rates of ranid tadpoles chronically exposed to low pH. Copeia. 1985;1985:415–23.

    Article  Google Scholar 

  88. Böhmer J, Rahmann H. Influence of surface water acidification on amphibians. Fort S Zool Biol Physiol Amphibians. 1990;38:287–309.

    Google Scholar 

  89. Panabokke CR. Soils and agro-ecological environments of Sri Lanka. Colombo: Natural Resources, Energy and Science Authority of Sri Lanka; 1996.

    Google Scholar 

  90. Silva EI, Manuweera L. Surface and rainwater chemistry in Sri Lanka–A risk of acidification. Asian J Water Environ Pollut. 2004;1(1, 2):79–86.

    CAS  Google Scholar 

  91. Newman RA. Adaptive plasticity in amphibian metamorphosis. BioSci. 1992;42:671–8.

    Article  Google Scholar 

  92. Berven KA. Factors affecting population fluctuations in larval and adult stages of the wood frog (Rana sylvatica). Ecology. 1990;71:1599–608.

    Article  Google Scholar 

  93. Martof B. Growth and development of the green frog, Rana clamitans, under natural conditions. Am Midl Nat. 1956;55:101–17.

    Article  Google Scholar 

  94. Smith DC. Factors controlling tadpole populations of the chorus frog (Pseudacris triseriata) on isle Royale, Michigan. Ecology. 1983;64:501–10.

    Article  Google Scholar 

  95. Wedekind C, Müller R. Risk-induced early hatching in salmonids. Ecology. 2005;86:2525–9.

    Article  Google Scholar 

  96. Semlitsch RD, Wilbur HM. Effects of pond drying time on metamorphosis and survival in the salamander Ambystoma talpoideum. Copeia. 1988;1988:978–83.

    Article  Google Scholar 

  97. Warkentin KM, Currie CR, Rehner SA. Egg-killing fungus induces early hatching of red-eyed treefrog eggs. Ecology. 2001;82:2860–9.

    Article  Google Scholar 

  98. Semlitsch RD, Scott DE, Pechmann JH. Time and size at metamorphosis related to adult fitness in Ambystoma talpoideum. Ecology. 1988;69:184–92.

    Article  Google Scholar 

  99. Tattersall GJ, Wright PA. The effects of ambient pH on nitrogen excretion in early life stages of the American toad (Bufo americanus). Comp Biochem Phys A: Physiol. 1996;113:369–74.

    Article  CAS  Google Scholar 

  100. Madeira D, Narciso L, Cabral HN, Vinagre C, Diniz MS. Influence of temperature in thermal and oxidative stress responses in estuarine fish. Comp Biochem Phys A: Mol Integr Physiol. 2013;166:237–43.

    Article  CAS  Google Scholar 

  101. Horne MT, Dunson WA. Toxicity of metals and low pH to embryos and larvae of the Jefferson salamander, Ambystoma jeffersoniannum. Arch Environ Con Tox. 1995;29:110–4.

    Article  CAS  Google Scholar 

  102. Clark KL, LaZerte BD. Intraspecific variation in hydrogen ion and aluminum toxicity in Bufo americanus and Ambystoma maculatum. Can J Fish Aquat Sci. 1987;44:1622–8.

    Article  CAS  Google Scholar 

  103. Watson G. An investigation into the development of upland and lowland tadpoles of the common frog (Rana temporaria) and common toad (Bufo bufo) when cultured at different temperatures. Doctoral dissertation. Durham: Durham University; 1992.

    Google Scholar 

  104. Drakulić S, Feldhaar H, Lisičić D, Mioč M, Cizelj I, Seiler M, Spatz T, Rödel M-O. Population-specific effects of developmental temperature on body condition and jumping performance of a widespread European frog. Ecol Evol. 2016;6:3115–28.

    Article  PubMed  PubMed Central  Google Scholar 

  105. Drakulić S, Feldhaar H, Lisičić D, Mioč M, Cizelj I, Seiler M, Spatz T, Rödel M-O. Local differences of thermal preferences in European common frog (Rana temporaria Linnaeus, 1758) tadpoles. Zool Anz. 2017;268:47–54.

    Article  Google Scholar 

  106. Beuchat CA, Pough FH, Stewart MM. Response to simultaneous dehydration and thermal stress in three species of Puerto Rican frogs. J Comp Physiol B. 1984;154:579–85.

    Article  Google Scholar 

  107. Niehaus AC, Wilson RS, Seebacher F, Franklin CE. Striped marsh frog (Limnodynastes peronii) tadpoles do not acclimate metabolic performance to thermal variability. J Exp Biol. 2011;214:1965–70.

    Article  PubMed  Google Scholar 

  108. Day T, Rowe L. Developmental thresholds and the evolution of reaction norms for age and size at life-history transitions. Am Nat. 2002;159:338–50.

    Article  PubMed  Google Scholar 

  109. Hentschel BT. Complex life cycles in a variable environment: predicting when the timing of metamorphosis shifts from resource dependent to developmentally fixed. Am Nat. 1999;154:549–58.

    Article  PubMed  Google Scholar 

  110. Wilbur HM, Collins JP. Ecological aspects of amphibian metamorphosis: nonnormal distributions of competitive ability reflect selection for facultative metamorphosis. Science. 1973;182(4119):1305–14.

    Article  CAS  PubMed  Google Scholar 

  111. Freitas JS, Almeida EA. Antioxidant defense system of tadpoles (Eupemphix nattereri) exposed to changes in temperature and pH. Zool Sci. 2016;33:186–95.

    Article  CAS  Google Scholar 

  112. McDonald DG, Wood CM. Branchial and renal acid and ion fluxes in the rainbow trout, Salmo gairdneri, at low environmental pH. J Exp Biol. 1981;93:101–18.

    CAS  Google Scholar 

  113. Booth JH, Jansz GF, Holeton GF. Cl−, K+, and acid–base balance in rainbow trout during exposure to, and recovery from, sublethal environmental acidification. Can J Zool. 1982;60:1123–30.

    Article  CAS  Google Scholar 

  114. McDonald DG. The interaction of environmental calcium and low pH on the physiology of the rainbow trout, Salmo gairdneri: I. Branchial and renal net ion and H+ fluxes. J Exp Biol. 1983;102:123–40.

    Google Scholar 

  115. Ultsch GR, Ott ME, Heisler N. Acid-base and electrolyte status in carp (Cyprinus carpio) exposed to low environmental pH. J Exp Biol. 1981;93:65–80.

    CAS  PubMed  Google Scholar 

  116. Kormanik GA, Cameron JN. Ammonia excretion in animals that breathe water: a review. Mar Biol Lett. 1981;2:11–23.

    CAS  Google Scholar 

  117. Randall DJ, Wright PA. The interaction between carbon dioxide and ammonia excretion and water pH in fish. Can J Zool. 1989;67:2936–42.

    Article  CAS  Google Scholar 

  118. Wright PA, Randall DJ, Perry SF. Fish gill water boundary layer: a site of linkage between carbon dioxide and ammonia excretion. J Comp Physiol B. 1989;158:627–35.

    Article  Google Scholar 

  119. Rajakaruna RS, Piyatissa PM, Jayawardena UA, Navaratne AN, Amerasinghe PH. Trematode infection induced malformations in the common hourglass treefrogs. J Zool. 2008;275:89–95.

    Article  Google Scholar 

  120. Mazzoni R, Cunningham AA, Daszak P, Apolo A, Perdomo E, Speranza G. Emerging pathogen in wild amphibians and frogs (Rana catesbeiana) farmed for international trade. Emerg Infect Dis. 2003;9:995–8.

    Article  PubMed  PubMed Central  Google Scholar 

  121. Dodd CE. The effects of acidic conditions and tannins on the survival, development, and behavior of Cope’S Gray Treefrog (Hyla chrysoscelis); 2018. p. 379. Electronic Theses and Dissertations. https://egrove.olemiss.edu/etd/379

    Google Scholar 

  122. Farquharson C, Wepener V, Smit NJ. Acute and chronic effects of acidic pH on four subtropical frog species. Water SA. 2016;42:52–62.

    Article  CAS  Google Scholar 

  123. Dittrich C, Drakulić S, Schellenberg M, Thein J, Rödel MO. Some like it hot? Developmental differences in yellow-bellied toad (Bombina variegata) tadpoles from geographically close but different habitats. Can J Zool. 2016;94:69–77.

    Article  Google Scholar 

Download references

Acknowledgements

We thank the Department of Zoology and Environmental Management for facilitating the work, Gihan Athapaththu and Krishan Rajapakse assisted in experimental set up. Dr. Anupama Halmillawewa, Dr. Renuka Attanayake and Dr. Shashikala Rajapakse provided microbial cultures for assessing lysozyme activity. Prof. D.M. de Costa facilitated the conduction of catalase assay and Dr. Saminda Fernando provided information on frog species in Sri Lanka.

Author information

Authors and Affiliations

Authors

Contributions

WAMTW conceived, designed and conducted the experiment, analysed the results and wrote the manuscript under the supervision and advice of RPGKR. RPGKR designed and supervised the study, reviewed and edited the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Gayani Rajapaksa.

Ethics declarations

Ethics approval

Ethics approval for this work was obtained from the Department of Zoology and Environmental Management of the University of Kelaniya, Sri Lanka.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1:

Table S1. Second-order polynomial functions fitted to variation of tadpole morphometrics with time (wk).

Additional file 2:

Table S2. Growth rates of tadpole morphometric characters at different developmental stages.

Additional file 3:

Table S3. Significance of contrasts comparing cumulative tadpole mortality in elevated CO2 with that of other treatments.

Additional file 4:

Table S4. Significance of contrasts comparing cumulative tadpole mortality in elevated temperatures with that of others.

Additional file 5:

Table S5. Significance of contrasts comparing ammonia excretion of tadpoles in elevated temperatures with that of others.

Additional file 6:

Table S6. Significance of contrasts comparing tadpole ammonia excretion in elevated CO2 with that of other treatments.

Additional file 7:

Plate S1. Lytics zones on gel plates 48 h after innoculation. Con - Ambient CO2 (water pH = 7) and water temperature at 29 ± 1 °C; ECO2 – CO2 bubbled to water to maintain pH at 5.5–5.6. R1 – R4 – Replicates.

Rights and permissions

Open Access This 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Weerathunga, W.A.M.T., Rajapaksa, G. The impact of elevated temperature and CO2 on growth, physiological and immune responses of Polypedates cruciger (common hourglass tree frog). Front Zool 17, 3 (2020). https://doi.org/10.1186/s12983-019-0348-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12983-019-0348-3

Keywords