Open Access
How to translate text using browser tools
21 July 2020 Density-Dependent Fitness Attributes and Carry-Over Effects in Crawfish Frogs (Rana areolata), a Species of Conservation Concern
Rochelle M. Stiles, Vanessa C. K. Terrell, John C. Maerz, Michael J. Lannoo
Author Affiliations +
Abstract

Fitness attributes acquired in aquatic habitats by amphibians exhibiting complex life histories have been shown to cascade through terrestrial juveniles into adulthood, a phenomenon termed carry-over effects. We explored density-dependent fitness attributes and carry-over effects in Crawfish Frogs (Rana areolata) using a set of field enclosure experiments and a series of field data. Using field enclosures, we hypothesized that 1) at high densities, intraspecific competition would produce smaller Crawfish Frog juveniles that took longer to metamorphose; 2) at high densities, interspecific competition would also produce smaller Crawfish Frog juveniles that took longer to metamorphose; and 3) vertebrate (ambystomatid salamander) predation on Crawfish Frog larvae would reduce survivorship, but by releasing competition pressure would produce relatively larger tadpoles that metamorphosed earlier. Further, we hypothesized 4) that these enclosure results would apply to field data, and that fitness attributes in newly metamorphosed Crawfish Frogs would carry over to first-time breeding adults. Our results confirmed all four hypotheses. Specifically, in Crawfish Frogs, at high densities, both intra- and interspecific competition reduced size (length and mass) at metamorphosis (hypotheses 1 and 2), and predation reduced survivorship and increased size at metamorphosis (hypothesis 3). Finally, we observed density-dependent fitness effects on newly metamorphosed Crawfish Frog juvenile size (length and mass), and carry-over effects from the larval stage on juvenile survival, adult size, and breeding adult numbers (hypothesis 4). In the absence of predators, high densities of intra- and interspecific competitors had no effect on Crawfish Frog larval survivorship. We also present suggestive evidence for compensatory effects. We discuss the potential mechanisms underlying the patterns of these interactions, as well as the role of these relationships in informing management guidelines intended to ensure the future of this species of conservation concern.

In north America north of Mexico, about two-thirds of the roughly 356 recognized species of amphibians exhibit some form of complex life history, characterized by aquatic eggs and larvae, and, following an abrupt metamorphosis, terrestrial juveniles and adults (Lannoo, 2005; Green et al., 2014; AmphibiaWeb, 2018). In amphibians, fitness attributes (sensu Hunt and Hodgson, 2010) acquired by aquatic larvae have been shown to cascade through terrestrial juveniles to influence adult phenotypes (Semlitsch et al., 1988; Beck and Congdon, 2000; Morey and Reznick, 2001), a phenomenon called carry-over effects (Van Allen et al., 2010; Harrison et al., 2011; Earl and Semlitsch, 2013). Berven (1990) described carry-over effects on amphibian survival and fitness as follows. The negative effects of egg or tadpole density on time of metamorphosis and metamorphic size carry over into juvenile success; larger juveniles and those metamorphosing early enjoy higher survival, earlier age at first reproduction, and are larger as adults. Larger adult size subsequently affects fitness by increasing the numbers and size of eggs produced by females (Berven, 1990). Carryover effects can also influence post-metamorphic survival (Goater, 1994; Beck and Congdon, 2000; Chelgren et al., 2006), activity patterns (Yagi and Green, 2018), and locomotion (John-Alder and Morin, 1990; Álvarez and Nicieza, 2002; Cabrera-Guzmán et al., 2013).

Wilbur (1997) summarized the effects of this density dependence in aquatic amphibian larvae by noting that size at metamorphosis is an exponentially decreasing function of initial density of the population. He observed that at low density, many individuals metamorphose at a large size, while at high densities, most individuals metamorphose at what may be the minimum metamorphic size threshold. Further, he suggested that survival is affected, because as density increases, reduced growth leads to a decreasing probability that a small, less-developed animal will achieve the minimum metamorphic size before the pond dries. Subsequent studies have generally confirmed these effects of larval density on size at and timing to metamorphosis (Semlitsch and Caldwell, 1982; Smith, 1983; Altwegg, 2003; Williams et al., 2012).

Crawfish Frogs (Rana areolata) are large-bodied ranids that breed in the early spring and whose tadpoles metamorphose by mid- to late-summer (Thompson, 1915; Gloyd, 1928; Cagle, 1942). They are considered a species of conservation concern in each of the 13 states they inhabit (Lannoo et al., 2017; Lannoo and Stiles, 2020). As the conservation status of R. areolata becomes more tenuous (of the 243 counties known to have been inhabited across their historic range, today Crawfish Frog populations remain in 159 counties, a 35% decline in occupancy; Lannoo et al., 2017), there is an urgency to understand the population dynamics of this species in order to enact management practices that will prevent further extirpations of these highly specialized animals.

Towards this end, we conducted a set of field enclosure experiments and collected a series of field data to test four hypotheses. We predicted that:

MATERIALS AND METHODS

Study site.—Our study site was located on the 729-hectare Hillenbrand Fish and Wildlife Area-West (HFWA-W), approximately 5 km south of Jasonville, in Greene County, Indiana. Prior to EuroAmerican settlement, this region was eastern deciduous forest punctuated by pocket prairies (Transeau, 1935). With settlement, this landscape was converted to agriculture then farmed for over a century before being surface mined for coal from 1976–1982. Post-mining, the soil was re-contoured and seeded to herbaceous vegetation in accordance with the federal Surface Mining Control and Reclamation Act (SMCRA; P.L. 95-87, enacted on August 3, 1977). The state purchased this property in 1988, and since then the Indiana Department of Natural Resources' (IDNR) Division of Fish and Wildlife has managed this land as native prairie (Lannoo et al., 2009; Terrell et al., 2014).

Most of our Crawfish Frog demographic data came from a previously unnamed wetland we called Nate's Pond. This wetland was formed unintentionally by soil slumping due to uneven overburden compaction during the reclamation process. Nate's Pond functions as a semi-permanent wetland, holding water throughout the year in most years, but it dries during droughts. The only time we observed it dry completely when Crawfish Frogs were present was during the 2012 drought, when it dried on 30 May and eliminated that year's cohort of tadpoles.

Study species.—Crawfish Frogs have long been considered the most secretive ranid species in North America (Wright and Myers, 1927) because they occupy abandoned crayfish burrows and typically never range far from them, except to breed (Heemeyer and Lannoo, 2012; Heemeyer et al., 2012). Early field biologists observed that Crawfish Frogs bred in the spring and metamorphosed in July and August, suggesting a larval period of about 90 days (Thompson, 1915; Wright and Meyers, 1927; Gloyd, 1928; Wright and Wright, 1933; Cagle, 1942).

At our study site, peak Crawfish Frog breeding ranged from 14 March to 4 April, but was generally centered around 1 April (Lannoo et al., 2017). Larval period length varied from 71 to 116 days (Table 1). This large variation meant peak metamorphosis occurred during a six-week period ranging from the middle of June through late July.

Table 1

Crawfish Frog field data (collected at Hillenbrand Fish and Wildlife Area-West, Greene County, Indiana) showing effects of larval density on juvenile size, time to metamorphosis, survivorship, and breeding adult size.

img-z3-2_443.gif

Nearby, in southern Illinois, the majority of breeding Crawfish Frogs in one study population were either three or four years old, and none of the 59 breeding adults was more than five years old (Redmer, 2000). Male Crawfish Frogs first breed at two years old, females at three (Redmer, 2000). Over the eight years of our study, we processed 1,102 breeding Crawfish Frog adults consisting of 729 males, 371 females, and 2 unsexed frogs. Confirming Redmer's (2000) results, in our study roughly two-thirds of all Crawfish Frog males first bred in their second year, with the remainder breeding in their third year. Further, only one-third of all Crawfish Frog females first bred in their second year, with the remaining two-thirds breeding in their third year (Lannoo and Stiles, 2020). Breeding adults averaged 92 mm (males) and 97 mm (females) in length (SVL, snout–vent length), and 93 g (males) and 82 g (females). The longest frog we captured was a gravid female at 120 mm (183 g, unknown age). The heaviest frog we captured was a gravid female at 193 g (116 mm, unknown age). The smallest breeding male frog was 61 mm and 33 g (two years old); the smallest breeding female was 71 mm and 42 g gravid (unknown age). Survivorship varied by life history stage, as follows: eggs, ∼98%; larvae, ∼1%; juveniles, ∼3%; adults, ∼70% (Lannoo and Stiles, 2020). From these data, it is clear that the highest mortality in Crawfish Frogs occurs during the larval and juvenile stages.

Size at metamorphosis varied. Wright and Wright (1933) offer that the average size of newly metamorphosed juveniles is ∼30 mm. In our study, in every year except 2014, juveniles were larger than this (Table 1). However, in 2013, a population of Crawfish Frogs on private property, 5 km from our study site, which bred at the same time as frogs in our population, produced juveniles in shallow prairie swale wetlands that metamorphosed a full month earlier and nearly a centimeter shorter (25 mm SVL; n = 51; Lannoo and Stiles, 2020) than juveniles at Nate's Pond.

At our study site, Crawfish Frog tadpoles develop syntopically with three other ranid species—Green Frogs, American Bullfrogs (R. catesbeiana), and Southern Leopard Frogs. Green Frogs and Southern Leopard Frogs are the most abundant of the four ranids (Lannoo et al., 2009), and their tadpoles likely compete with developing Crawfish Frog tadpoles (Terrell et al., unpubl.). At HFWA-W, Green Frogs and Bullfrogs breed during the summer; their tadpoles overwinter and complete metamorphosis the following summer (Minton, 1972). Southern Leopard Frogs may breed in the fall or spring (Lannoo, 2005; Stiles and Lannoo, 2015), but in 2011, when we conducted our enclosure study, adults bred only in the spring, syntopically with Crawfish Frogs. Due to these natural history differences, at the beginning of the enclosure study, second-year Green Frog tadpoles were substantially larger (2.0–2.5 cm SVL) than either recently hatched Crawfish Frog tadpoles or Southern Leopard Frog tadpoles, which were similarly sized (∼0.5 cm SVL).

Enclosure studies.—To test the effects of high density (both intraspecific and interspecific) and vertebrate predation on the growth and survivorship of Crawfish Frog tadpoles, we used 379 L aquatic mesh enclosures (Reptariums®, Apogee, Dallas, TX; filled to ∼340 L). On 27 April 2011 (day 0), we placed 15 enclosures in the shallows of Erosion Control Pond at HFWA-W. We separated tadpoles into five treatments, each with three replicates, as follows: 1) 30 Crawfish Frogs (CF; density of 0.09 tadpoles/L); 2) 90 CF (0.27 tadpoles/L); 3) 45 CF and 45 Leopard Frogs (0.13 tadpoles/L/species); 4) 45 CF and 20 Green Frogs (0.13 CF tadpoles/L; 0.06 GF tadpoles/L); and 5) 90 CF and 5 Small-mouthed Salamander larvae (0.27 CF tadpoles/L; 0.02 SS larvae/L). The enclosures were placed linearly along the bank, blocked by replicates containing treatments. Each replicate had one enclosure of each treatment randomly assigned at the start of the experiment within the block. Each enclosure was situated such that a minimum 7.5 cm gap was present between the water surface inside the enclosure and the mesh top, to ensure tadpoles had access to atmospheric oxygen. To accomplish this, enclosures were slid down-bank as this wetland dried, and up-bank following heavy rains. Enclosure placements were never adjusted more than a meter up or down bank, and never more often than once a week.

For the density treatments, we obtained Crawfish Frog and Southern Leopard Frog tadpoles from Cattail Pond. The wetland is a population sink. Each year of our study, from 2009–2016, between 20 and 40 Crawfish Frogs bred in this wetland, but the egg masses deposited (all of which were viable) produced few Crawfish Frog juveniles (46 total over eight years), which in turn generated no breeding adults. We collected tadpoles of both species from egg masses laid on the same night. We trapped Green Frog tadpoles at the same HFWA-W wetland (Erosion Control) where we conducted the treatments. Erosion Control Pond was a recently constructed (2008), permanent wetland naturally colonized by Southern Leopard Frogs, Green Frogs, and Bullfrogs in 2009, and Crawfish Frogs in 2010 (Lannoo et al., 2009; Lannoo and Stiles, 2020). The western end of this wetland is deep, and always held water, while the eastern end, where we placed our enclosures, was shallow (<0.5 m) and would often dry by late summer. Before placing tadpoles in enclosures, we reared Crawfish Frog and Southern Leopard Frog tadpoles in small, floating aquaria until the animals were free swimming. We fed these tadpoles frozen, chopped spinach ad libitum prior to assigning them to treatments. For the predation treatment, we trapped or seined HFWA-W wetlands for Small-mouthed Salamander (Ambystoma texanum) larvae.

We supplied each enclosure with 454 g per week of frozen, chopped spinach, a rate V.C.K.T. considered, based on her undergraduate experiences with mesocosms, sufficient for growth but creating competition. (Indeed, we found that, typically, throughout most of the experiment, spinach would persist from one feeding to the next in the low-density intraspecific treatment but be fully consumed in the high-density treatment.) We examined each enclosure daily for metamorphosed frogs and salamanders, as well as for any predators such as odonate naiads that might have negotiated the mesh or the zippered openings (we never found any). As we collected metamorphosing juveniles, we measured their SVL (mm) and mass (g), and provided an identification mark (toe-clip) prior to release near their natal wetland. We ended the experiment after all the tadpoles had metamorphosed.

Field studies.—From 2009 through 2016, we set up a drift fence/pitfall trap array around Nate's Pond (details in Kinney, 2011; Lannoo et al., 2017), which allowed us to trap breeding Crawfish Frog adults and newly metamorphosed juveniles. Given the low trespass rate of breeding adults (<1%; Kinney, 2011), and the small size and restricted locomotor capacity of newly metamorphosed juveniles, we consider our drift fence data to be a census rather than a sampling of this population (Alford et al., 2001; Lannoo et al., 2017). We weighed and measured each captured adult and juvenile, except in 2013, where the large number of juveniles and time constraints forced us to process every tenth animal. We inserted a passive integrated transponder (PIT) tag into each captured Crawfish Frog adult, which allowed us to identify individuals. We toe-clipped juveniles to indicate year and pond, which allowed us to identify cohorts. We considered the number of days to metamorphosis to be from the date of peak breeding to the date tadpoles completed metamorphosis (Stage 46; Gosner, 1960). We calculated estimated number of eggs laid using the following regression equation (Kinney, 2011), based on data presented in Redmer (2000):

e01_443.gif

When assessing adult size across cohorts, we only used data from 2-year-old, first-time breeding animals (to avoid confounding our results by mixing size data from these animals and larger 3-year-old animals). When assessing juvenile survivorship to first reproduction, we used data from all known first-time breeding adults produced by a cohort. These animals included both 2- and 3-year-old first-time breeders, plus any adults collected subsequently (one female bred at Nate's Pond every year of our study [from 2009–2016], which made her at least ten years old; Lannoo et al., 2017; Lannoo and Stiles, 2020). Because we terminated our study in 2016, in order to estimate juvenile survivorship for the 2014 cohort, we calculated the ratio of 2-year-old breeders to all other adults collected for the 2009, 2010, and 2011 cohorts (0.99:1.00), then used this ratio to estimate survivorship from metamorphosis to first reproduction in the 2014 cohort.

Statistics.—In our enclosure studies, to determine whether competition and predation affected survival, length, mass, and days to metamorphosis, we compared treatments using one-way analysis of variance (ANOVA) and post hoc Tukey's honestly significant difference tests in Program R version 3.3.0 (R Core Team, 2013). We transformed survivorship data using an arcsine function to meet the parametric assumption of a normally distributed sample. For these results, we report means with 95% confidence intervals and statistical significance at P < 0.05.

Our field data, representing only four breeding years (cohorts; see below) of useable data, but encompassing 12 variables (Table 1), presented the statistical challenge of a small sample size relative to the number of variables being considered. We eliminated potential explanatory variables based on correlations with other variables (e.g., number of females breeding and female size was positively correlated with estimated eggs deposited, and juvenile mass was positively correlated with juvenile SVL). We chose candidate models a priori based on knowledge acquired from our eight years of field data collected on Crawfish Frogs (as recommended by Burnham and Anderson, 2010).

We divided our analysis into direct density effects (effects of larval factors on newly metamorphosed juvenile characteristics) and carry-over effects (effects of larval factors on older juvenile and breeding adult characteristics; Table 2). To assess direct density effects, we used the explanatory variables of estimated number of eggs laid (a proxy for early density) and cohort size (juvenile number, a proxy for late density), on juvenile SVL, juvenile mass, and days to metamorphosis (Table 2A). To address carry-over effects, we used the explanatory variables of cohort size and juvenile SVL on the response variables of survival to first breeding, adult SVL, adult mass, and adult number (Table 2B). We constructed simple linear regression models in Program R. As indicated above, we excluded data from two years with low juvenile numbers (2012, n = 0; 2013, n = 8), and as a result of small sample sizes, we only included one predictor variable (df = 2) in each model. We present R2 values (i.e., the variation in the dependent variable that can be explained by the explanatory variable) and P-values (with significance indicated by P < 0.05). To compare direct density effects with carry-over effects, we used the common explanatory variable of cohort size.

Table 2

Linear regression models constructed in program R assessing fitness and recruitment. Data are from Nate's Pond cohorts 2009, 2010, 2011, 2014. Abbreviations: #Eggs = estimated number of eggs deposited; Cohort size = number of juveniles; Days = days to metamorphosis; and Survival = survivorship from metamorphosis to first breeding. * P < 0.05.

img-z4-9_443.gif

RESULTS

High intraspecific density effects.—In our enclosure studies, the high-density treatment among Crawfish Frog tadpoles decreased growth, increased the number of days to metamorphosis, and had no effect on larval survivorship (Fig. 1AD). Newly metamorphosed juveniles from the high-density tadpole treatment (0.27 tadpoles/L) were on average smaller and weighed less than juveniles from the low-density tadpole treatment (0.09 tadpoles/L; ANOVA: SVL, F4,10 = 18.73, P < 0.001; mass, F4,10 = 30.90, P < 0.001; Fig. 1B, C). Further, tadpoles raised at the higher density required an average of five days longer to reach metamorphosis than tadpoles raised at the lower density (F4,10 = 17.27, P < 0.001; Fig. 1D). Larval survivorship was unaffected by high intraspecific competition; there was no statistical difference in survivorship between tadpoles raised at high and low densities (F4,10 = 7.35, P = 0.01; Tukey's HSD: P = 0.94; Fig. 1A).

Fig. 1

Results of enclosure experiments to determine effects of intraspecific competition, interspecific competition, and predation treatments on larval survivorship, size (SVL and mass), and time to metamorphosis. In the two Intraspecific treatments, Low = 30 Crawfish Frog tadpoles/enclosure (0.09 tadpoles/L) and High = 90 Crawfish Frog tadpoles/enclosure (0.27 tadpoles/L). In the two Interspecific treatments, Lc = Green Frog tadpoles and Ls = Southern Leopard Frog tadpoles. In these treatments, in each enclosure we placed 45 Crawfish Frog and 45 Leopard Frog tadpoles (0.13 tadpoles/L/species), and 45 Crawfish Frog (0.13 tadpoles/L) and 20 Green Frog tadpoles (0.06 tadpoles/L), respectively. In the Predation treatment (At), we placed 90 Crawfish Frog tadpoles (0.27 tadpoles/L) and 5 Small-mouthed Salamander larvae (0.02 larvae/L) per enclosure. Lowercase letters above each column reflect statistical significance using ANOVA. Columns with different letters are statistically significant (P < 0.05).

img-z5-1_443.jpg

High interspecific density effects.—In our enclosures, Crawfish Frogs reared with Green Frogs and Southern Leopard Frogs were smaller (ANOVA: SVL, F4,10 = 18.73, P < 0.001; weight, F4,10 = 30.90, P < 0.001) and required more time to reach metamorphosis (F4,10 = 17.27, P < 0.001) than Crawfish Frogs reared alone at the low density (Fig. 1BD). Crawfish Frogs reared with Green Frogs were a similar size and weight at metamorphosis as those reared with Southern Leopard Frogs. We found no statistically significant differences across interspecific competition treatments in Crawfish Frog growth (Tukey's HSD: SVL, P = 1.00; mass, P = 1.00) and time to metamorphosis (P = 0.31; Fig. 1). Larval survivorship was high in Crawfish Frogs reared at high densities with both Green Frogs and Southern Leopard Frogs, with no statistically significant species effect (Tukey's HSD: P = 0.95; Fig. 1A).

Predation effect.—In our enclosure study, predation by Small-mouthed Salamanders significantly reduced survivorship of Crawfish Frog tadpoles (ANOVA: F4,10 = 7.35, P < 0.01; Fig. 1A). However, Crawfish Frog tadpoles reared with salamanders produced larger juveniles that weighed more than juveniles raised in the high-density intraspecific treatment (SVL, F4,10 = 18.73, P < 0.001; mass, F4,10 = 30.90, P < 0.001). These juveniles were, however, smaller and weighed less than in the low-density intraspecific treatment (Fig. 1B, C). Crawfish Frog tadpoles reared with salamanders metamorphosed significantly faster than tadpoles reared in either the high-density intraspecific treatment or co-reared with Southern Leopard Frogs (F4,10 = 17.27, P < 0.001; Fig. 1D).

Field data: number, size, and survivorship.—The number of breeding females at Nate's Pond varied from 8 to 52 (Table 1A). On average, the largest females bred in 2012, the smallest in 2013. We estimate the highest number of eggs was laid in 2013, the lowest in 2012. Survivorship to metamorphosis (larval survivorship) ranged from 1.7% of estimated eggs laid to 0 (100% mortality).

The number of newly metamorphosed Crawfish Frog juveniles at Nate's Pond varied from 286 to 3,122 (Table 1B). Juveniles were longest and heaviest in 2009, and shortest and lightest in 2014. The number of days required for metamorphosis was fewest in 2014 and greatest in 2013. As indicated above, we excluded cohorts with low juvenile numbers (2012, n = 0; 2013, n = 8) in our analyses.

The number of breeding adult Crawfish Frogs produced by each cohort also varied, from 8 (2.8% survivorship) to 87 (4.1%; Table 1B, C). Based on the ratio of 2-year-old breeders to all other adults, we estimate an additional 21 frogs from the 2014 cohort (making a total of 42 frogs) would have bred if our study continued beyond 2016. This estimate increases the survivorship of the 2014 cohort from 2.5% to 5.0%. The size of breeding adults produced also varied by cohort. The 2009 cohort produced the longest adults, while 2014 produced the heaviest adults (Table 1C). The 2011 cohort produced both the shortest and the lightest adults.

Field data: direct density dependence and carry-over effects.—From a statistical perspective, there were no significant effects of estimated number of eggs laid or cohort size on juvenile SVL and mass, or days to metamorphosis (Table 2A). Similarly, with one exception (the inverse effects of cohort size on adult SVL; R2 = 0.94, df = 0.94, P = 0.03), there were no significant effects of cohort size and juvenile SVL on survival to first breeding, adult SVL or mass, or adult number (Table 2B). A visualization of these data, however, suggests that the relationships between cohort size and juvenile and adult characteristics are stronger than the power of our statistical tests to detect them (Fig. 2). In particular, in 2009, 2010, and 2011, but not 2014, cohort size had a clear negative effect on juvenile length (Fig. 2A) and mass (Fig. 2C). In 2010, 2011, and 2014, but not 2009, cohort size had a clear positive effect on days to metamorphosis (Fig. 2E). These relationships between cohort size and body size carried over to adult length (Fig. 2B) and mass (Fig. 2D). There was no correlation between cohort size and survivorship from metamorphosis to first breeding (R2 = 0.04, df = 2, P = 0.79; Fig. 2F).

Fig. 2

Field data showing effects of juvenile density on juvenile body size (A, C) and time to metamorphosis (E), and breeding adult size (B, D) and survivorship (F). Lines reflect regression calculations. While suggestive, none of these trends is significant except for the adult SVL data shown in B.

img-z6-1_443.jpg

DISCUSSION

Using cage experiments, we tested three hypotheses addressing the effects of high-density intra- and interspecific competition, and predation, on Crawfish Frog larvae. Using field data, we then tested a fourth hypothesis on cohort size on newly metamorphosed juvenile size, survivorship, and length of larval period (density-dependent effects). As well, we assessed the effects on adult breeding size and survivorship (carry-over effects). We consider these hypotheses, below.

Density dependence and predation.—We hypothesized the high intraspecific density treatment would produce smaller Crawfish Frog juveniles that took longer to metamorphose. Further, we hypothesized high interspecific density would have the same negative effect on Crawfish Frog larvae as high intraspecific density. Our results generally confirmed these hypotheses—high intra- and interspecific densities slowed larval growth and increased time to metamorphosis (Fig. 1BD). Neither the high intraspecific or interspecific density treatments had an effect on Crawfish Frog larval survivorship (Fig. 1A). Our results echoed (Williams et al., 2012) and confirmed for Crawfish Frogs the effects of high larval densities on juvenile size and time to metamorphosis (Semlitsch and Caldwell, 1982; Smith, 1983; Altwegg, 2003). These results may have differed if we had introduced Crawfish Frog eggs into enclosures containing overwintered Green Frog tadpoles. Sours and Petranka (2007) have shown that ranids will prey on newly laid anuran eggs, and reduce survivorship to hatching.

We further hypothesized that ambystomatid salamander predation on Crawfish Frog larvae would be severe (e.g., Sredl and Collins, 1992; Relyea and Hoverman, 2003), but releasing competition pressure would produce larger larvae that metamorphosed early (Relyea and Rosenberger, 2018), which our data confirmed (Fig. 1B, C). However, the magnitude of this effect was less than observed in other studies (see below).

Importantly, we observed similar density effects in our field data. Cohort size negatively influenced newly metamorphosed Crawfish Frog juvenile length and mass (Fig. 2A, C), and positively influenced time to metamorphosis (Fig. 2E). (Outliers from 2014 [Fig. 2A, C] influenced our statistical results [Table 2]. Our suspicion is that food [algae] levels were reduced in 2014 compared to other years, but we did not measure this variable, and therefore can only speculate.) Conversely, cohort size positively influenced time to metamorphosis (in 2010, 2011, and 2014; Fig. 2E), consistent with other field studies (Semlitsch et al., 1988; Berven, 1990; Wilbur, 1997). (The 2009 cohort took much longer to metamorphose—potentially the result of lower than average water temperatures [once again, a variable we did not track].)

Carry-over and compensatory effects.—Our field data suggested both carry-over and compensatory effects from the larval to the juvenile and adult life history stages (Fig. 2). In particular, cohort size negatively affected juvenile size (in three years of our study), which in turn affected adult size and number (Fig. 2B, D). These suggested carry-over effects on Crawfish Frogs generally agree with results from studies on a range of other species (A. talpoideum [Semlitsch et al., 1988]; R. sylvatica [Berven, 1990]; A. opacum [Scott, 1994]; Bufo terrestris [Beck and Congdon, 2000]).

What is perhaps most interesting about these data, however, is that the 2014 outliers we observed in data detailing the effects of cohort size on juvenile SVL and mass (Fig. 2A, C) are no longer outliers when considering the effect cohort size has on 2-year-old breeding adult SVL and mass (Fig. 2B, D); they have been compensated for, similar to the findings of Boone (2005). As further evidence for compensatory effect, we note the slopes describing the relationship between cohort size and 2-year-old breeding adult SVLs and mass (Fig. 2B, D) are shallower than the slopes describing the relationship between cohort size and newly metamorphosed juvenile SVL and mass (Fig. 2A, C), suggesting that the high and low extremes seen in juvenile size metrics were moderated in first-time breeding adults.

How should we be assessing larval density?—Our regression models suggest that cohort size (a proxy for final density) correlated better with juvenile SVL and mass, and days to metamorphosis than did the estimated number of eggs laid (a proxy for initial density; Table 2A). As most aquatic amphibian biologists know, an accurate assessment of larval density is difficult to obtain, or even estimate. While the general trend is for seasonal and semi-permanent wetlands to fill with snowmelt and rain water in the early spring, then gradually dry as the warm season progresses, this is far from predictable on a day-to-day, or year-to-year basis (Semlitsch et al., 1996). Late-spring thunderstorms can refill drying basins, and in a wet spring following a dry winter, there may be more water present at the time of metamorphosis than there was three months earlier when eggs were laid (Euliss and Mushet, 1996). In general, drawdowns concentrate animals and increase density, while heavy predation reduces density (Wilbur et al., 1983; Sredl and Collins, 1992; Relyea and Rosenberger, 2018).

Moreover, most density estimates fail to consider animal size. It is possible to fit 1,000 newly hatched Crawfish Frog tadpoles in a liter of water, but impossible to fit 1,000 premetamorphic tadpoles in the same volume. Rather than number of animals per volume of water, a more accurate assessment of density effects might compare volumes of animals per volume of water. All of these uncertainties are in play when considering density-dependent effects on the animals experiencing them.

Despite Wilbur's (1997) observation that size at metamorphosis is an exponentially decreasing function of initial density of the population, field studies of density-dependent effects in amphibians with complex life histories have used juvenile number (final density) as the proxy for density (Semlitsch et al., 1988). In fact, during years with wetland drawdowns and low predation, it is likely that the time immediately prior to metamorphosis is the time of highest tadpole density. But this may not always be true, and other proxies for density, such as number of breeding females and estimated number of eggs deposited—which assess densities earlier in the larval period—might serve just as well (and be consistent with Wilbur, 1997).

How to define metamorphosis?—There is one final consideration regarding the timing of metamorphosis that has never been fully addressed, which certainly contributes noise to data such as ours. If metamorphosis is defined not when animals actually metamorphose, but when juveniles emerge and appear in drift fence buckets (technically, post-metamorphic dispersal), then metamorphosis as defined in studies such as ours (e.g., Semlitsch et al., 1988) is not some continuous, normally distributed function related to the completion of larval development, but rather the saltatory response of post-metamorphic animals to periodic nighttime rains—an environmental condition permitting distant dispersal across terrestrial landscapes hostile to animals prone to desiccation (Kinney, 2011; Lannoo et al., 2017). True metamorphosis likely occurs days, if not longer, before metamorphosis is signified by dispersing juveniles appearing in pitfall traps.

Conservation implications.—Our experimental data confirm the negative effects of high larval density on juvenile body size and time to metamorphosis in Crawfish Frogs (Williams et al., 2012; Fig. 1). Our field data are consistent with this relationship, and further show that the negative effects of high larval densities on juvenile body size carry over to influence breeding adult body size, although we also show evidence for offsetting compensatory effects (Fig. 2). With these facts in mind, we revisit the findings of Berven (1990). Given that larger adult size positively affects fitness by increasing the numbers of eggs produced by females, our data then suggest that the second consequence of increased adult fitness—that is the consequence of having increased numbers of eggs and therefore larvae—becomes diminished fitness, as density-dependent factors act on this increased number of larvae by making each animal smaller. These small larvae then produce small juveniles, which in turn produce smaller adults, which have comparatively reduced fitness. Following this logic, one consequence of increased adult fitness is, eventually, reduced adult fitness. We suppose by reverse logic the converse is also true—that smaller adults produce fewer eggs that produce larvae that, by not being constrained by density, grow bigger and produce juveniles that go on to produce larger adults that have increased fitness.

The question for our Crawfish Frog population, and pond-breeding amphibian populations in general, then becomes, which situation is better for population persistence: large numbers of small adults, or small numbers of large adults? Typically, amphibian population assessments are made based on body counts (e.g., Pechmann et al., 1991; Green, 2003), suggesting that population size is more important than body size when considering population viability. But, from the standpoint of reproductive potential, might a decline in population size be compensated for by an increase in the body size of the remaining individuals? The data we present here suggest they can, up to a point. In 2013, 52 small (88 mm SVL) Crawfish Frog females at Nate's Pond produced an estimated 211,828 eggs, while in 2009 only 31 large (100 mm SVL) females produced an estimated 196,478 eggs (Table 1A); that is, a reduction of 40% in the number of females breeding (31 vs. 52) produced a reduction of only 7% in the estimated number of eggs deposited (196,478 vs. 211,828). It is beyond the scope of this paper to formally consider this question, but if due to larval density-dependent factors there is an inverse relationship between adult size and adult number, a population decline may become irreversible when the increased reproductive potential of the remaining large adults is insufficient to overcome the reduced reproductive potential caused by the diminished population size.

ACKNOWLEDGMENTS

We thank K. Gremillion-Smith, S. Klueh-Mundy, S. Johnson, R. Ronk, B. Feaster, and T. Stoelting of the Indiana Department of Natural Resources for supporting this research at HFWA-W; P. Moler and D. Scott for collecting belostomatids; J. Swan, L. Sawyer, H. Nesius, S. Ronk, R. Ronk, D. Ronk, M. Goode, D. Ramesh, and D. Schaefer for assisting with experimental setup and data collection; and J. Klemish and S. Lannoo for comments on earlier drafts of this manuscript. This work was funded by a US Fish and Wildlife Service State Wildlife Grant (E2-08-WDS13) from the Indiana Department of Natural Resources, Division of Fish and Wildlife. This project was approved by and conducted in accordance with Indiana State University's Animal Care and Use (IACUC number 245168-1: ML) and the Indiana Scientific Collection (DNR permits 11-0017 and 14-063) guidelines.

LITERATURE CITED

1.

Alford, R. A., P. M. Dixon, and J. H. K. Pechmann. 2001. Global amphibian population declines. Nature 412:499–500. Google Scholar

2.

Alford, R. A., and H. M. Wilbur. 1985. Priority effects in experimental pond communities: competition between Bufo and Rana. Ecology 66:1097–1105. Google Scholar

3.

Altwegg, R. 2003. Multistage density dependence in an amphibian. Oecologia 136:46–50. Google Scholar

4.

Álvarez, D., and A. G. Nicieza. 2002. Effects of induced variation in anuran larval development on postmetamorphic energy reserves and locomotion. Oecologia 131:186–195. Google Scholar

5.

Amphibia Web. 2018. Amphibian data base.  https://amphibiaweb.org/search/index.html (accessed 1 June 2018). Google Scholar

6.

Beck, C. W., and J. D. Congdon. 2000. Effects of age and size at metamorphosis on growth and survivorship of Southern Toad (Bufo terrestris) metamorphs. Canadian Journal of Zoology 77:944–951. Google Scholar

7.

Berven, K. A. 1990. Factors affecting population fluctuations in larval and adult stages of the Wood Frog (Rana sylvatica). Ecology 71:1599–1608. Google Scholar

8.

Boone, M. D. 2005. Juvenile frogs compensate for small metamorph size with terrestrial growth: overcoming the effects of larval density and insecticide exposure. Journal of Herpetology 39:416–423. Google Scholar

9.

Burnham, K. P., and D. R. Anderson. 2010. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. Springer-Verlag, New York. Google Scholar

10.

Cabrera-Guzmán, E., M. R. Crossland, G. P. Brown, and R. Shine. 2013. Larger body size at metamorphosis enhances survival, growth, and performance of young Cane Toads (Rhinella marina). PLoS ONE 8:e70121. Google Scholar

11.

Cagle, F. R. 1942. Herpetological fauna of Jackson and Union counties, Illinois. American Midland Naturalist 28:164–200. Google Scholar

12.

Chelgren, N. D., D. K. Rosenberg, S. S. Heppell, and A. I. Gitelman. 2006. Carryover aquatic effects on survival of metamorphic frogs during pond emigration. Ecological Applications 16:627–636. Google Scholar

13.

Earl, J. E., and R. D. Semlitsch. 2013. Carryover effects in amphibians: are characteristics of the larval habitat needed to predict juvenile survival? Ecological Applications 23: 1429–1442. Google Scholar

14.

Euliss, N. H., Jr., and D. M. Mushet. 1996. Water-level fluctuation in wetlands as a function of landscape condition in the Prairie Pothole Region. Wetlands 16: 587–593. Google Scholar

15.

Gloyd, H. K. 1928. Amphibians and reptiles of Franklin County. Transactions of the Kansas Academy of Science 31: 116–119. Google Scholar

16.

Goater, C. P. 1994. Growth and survival of postmetamorphic toads: interactions among larval history, density, and parasitism. Ecology 75:2264–2274. Google Scholar

17.

Gosner, K. L. 1960. A simplified table for staging anuran embryos and larvae with notes on identification. Herpetologica 16:183–190. Google Scholar

18.

Green, D. M. 2003. The ecology of extinction: population fluctuation and decline in amphibians. Biological Conservation 111:331–343. Google Scholar

19.

Green, D. M., L. Weir, G. S. Casper, and M. J. Lannoo. 2014. North American Amphibians: Distribution and Diversity. University of California Press, Berkeley, California. Google Scholar

20.

Harrison, X. A., J. D. Blount, R. Inger, D. R. Norris, and S. Bearhop. 2011. Carry-over effects as drivers of fitness differences in animals. Journal of Animal Ecology 80:4–18. Google Scholar

21.

Heemeyer, J. L., and M. J. Lannoo. 2012. Breeding migrations in Crawfish Frogs (Lithobates areolatus): long-distance movements, burrow philopatry, and mortality in a near-threatened species. Copeia 2012:440–450. Google Scholar

22.

Heemeyer, J. L., P. J. Williams, and M. J. Lannoo. 2012. Obligate crayfish burrow use and core habitat requirements of Crawfish Frogs. Journal of Wildlife Management 76: 1081–1091. Google Scholar

23.

Hunt, J., and D. Hodgson. 2010. What is fitness, and how do we measure it?, p. 46–70. In : Evolutionary Behavioral Ecology. D. Westneat and C. W. Fox (eds.). Oxford University Press, Oxford. Google Scholar

24.

John-Alder, H. B., and P. J. Morin. 1990. Effects of larval density on jumping ability and stamina in newly metamorphosed Bufo woodhousii fowleri. Copeia 1990:856–860. Google Scholar

25.

Kinney, V. C. 2011. Adult survivorship and juvenile recruitment in populations of Crawfish Frogs (Lithobates areolatus), with additional consideration of the population sizes of associated pond breeding species. Unpubl. M.S. thesis, Indiana State University, Terre Haute, Indiana. Google Scholar

26.

Lannoo, M. J. (Ed.). 2005. Amphibian Declines: The Conservation Status of United States Species. University of California Press, Berkeley, California. Google Scholar

27.

Lannoo, M. J., V. C. Kinney, J. Heemeyer, N. J. Engbrecht, A. L. Gallant, and R. W. Klaver. 2009. Mine spoil prairies expand critical habitat for endangered and threatened amphibian and reptile species. Diversity 1:118–132. Google Scholar

28.

Lannoo, M. J., and R. M. Stiles. 2020. The Call of The Crawfish Frog. Taylor and Francis, London. Google Scholar

29.

Lannoo, M. J., R. M. Stiles, M. A. Sisson, J. W. Swan, V. C. K. Terrell, and K. E. Robinson. 2017. Patch dynamics inform management decisions in a threatened frog species. Copeia 105:53–63. Google Scholar

30.

Minton, S. A. 1972. Amphibians and Reptiles of Indiana. Indiana Academy of Science Monograph, No. 3. Indiana Academy of Science, Indianapolis, Indiana. Google Scholar

31.

Morey, S., and D. Reznick. 2001. Effects of larval density on postmetamorphic Spadefoot Toads (Spea hammondii). Ecology 82:510–522. Google Scholar

32.

Pechmann, J. H. K., D. E. Scott, R. D. Semlitsch, J. P. Caldwell, L. J. Vitt, and J. W. Gibbons. 1991. Declining amphibian populations: the problem of separating human impacts from natural fluctuations. Science 253:892–895. Google Scholar

33.

R Core Team. 2013. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.  https://www.R-project.org/  Google Scholar

34.

Redmer, M. 2000. Demographic and reproductive characteristics of a southern Illinois population of the Crayfish Frog, Rana areolata. Journal of the Iowa Academy of Science 107: 128–133. Google Scholar

35.

Relyea, R. A., and J. T. Hoverman. 2003. The impact of larval predators and competitors on the morphology and fitness of juvenile treefrogs. Oecologia 134:596–604. Google Scholar

36.

Relyea, R. A., and D. Rosenberger. 2018. Predator effects on metamorphosis: the effects of scaring versus thinning at high prey densities. Copeia 106:457–467. Google Scholar

37.

Scott, D. E. 1994. The effect of larval density on adult demographic traits in Ambystoma opacum. Ecology 75:296–306. Google Scholar

38.

Semlitsch, R. D., and J. P. Caldwell. 1982. Effects of density on growth, metamorphosis, and survivorship in tadpoles of Scaphiopus holbrookii. Ecology 63:905–911. Google Scholar

39.

Semlitsch, R. D., D. E. Scott, and J. H. K. Pechmann. 1988. Time and size at metamorphosis related to adult fitness in Ambystoma talpoideum. Ecology 71:184–192. Google Scholar

40.

Semlitsch, R. D., D. E. Scott, J. H. K. Pechmann, and J. W. Gibbons. 1996. Structure and dynamics of an amphibian community: evidence from a 16-year study of a natural pond, p. 217–248. In : Long-Term Studies of Vertebrate Communities. M. L. Cody and J. A. Smallwood (eds.). Academic Press, San Diego, California. Google Scholar

41.

Smith, D. C. 1983. Factors controlling tadpole populations of the Chorus Frog (Pseudacris triseriata) on Isle Royale, Michigan. Ecology 64:501–510. Google Scholar

42.

Sours, G. N., and J. W. Petranka. 2007. Intraguild predation and competition mediate stage-structured interactions between Wood Frog (Rana sylvatica) and Upland Chorus Frog (Pseudacris feriarum) larvae. Copeia 2007:131–139. Google Scholar

43.

Sredl, M. J., and J. P. Collins. 1992. The interaction of predation, competition, and habitat complexity in structuring an amphibian community. Copeia 1992:607–614. Google Scholar

44.

Stiles, R. M., and M. J. Lannoo. 2015. Lithobates sphenocephalus (Southern Leopard Frog). Fall breeding. Herpetological Review 46:414. Google Scholar

45.

Terrell, V. C. K., J. L. Klemish, N. J. Engbrecht, J. A. May, P. J. Lannoo, R. M. Stiles, and M. J. Lannoo. 2014. Amphibian and reptile recolonization of reclaimed coal spoil grasslands. Journal of North American Herpetology 2014:59–68. Google Scholar

46.

Thompson, C. 1915. Notes on the habits of Rana areolata Baird and Girard. Occasional papers of the Museum of Zoology, University of Michigan, Number 9. University of Michigan, Ann Arbor. Google Scholar

47.

Transeau, E. N. 1935. The prairie peninsula. Ecology 16:423–437. Google Scholar

48.

Van Allen, B. G., V. S. Briggs, M. W. McCoy, and J. R. Vonesh. 2010. Carry-over effects of the larval environment on post-metamorphic performance in two hylid frogs. Oecologia 164:891–898. Google Scholar

49.

Werner, E. E. 1992. Competitive interactions between Wood Frog and Northern Leopard Frog larvae. Ecology 77:157–169. Google Scholar

50.

Wilbur, H. M. 1977a. Density-dependent aspects of growth and metamorphosis in Bufo americanus. Ecology 58:196–200. Google Scholar

51.

Wilbur, H. M. 1977b. Interactions of food level and population density in Rana sylvatica. Ecology 58:206–209. Google Scholar

52.

Wilbur, H. M. 1997. Experimental ecology of food webs: complex systems in temporary ponds. Ecology 78:2279–2302. Google Scholar

53.

Wilbur, H. M., and R. A. Alford. 1985. Priority effects in experimental pond communities: responses of Hyla to Bufo and Rana. Ecology 66:1106–1114. Google Scholar

54.

Wilbur, H. M., P. J. Morin, and R. N. Harris. 1983. Salamander predation and the structure of experimental communities: anuran responses. Ecology 64:1423–1429. Google Scholar

55.

Williams, P. J., J. R. Robb, R. H. Kappler, T. E. Piening, and D. R. Karns. 2012. Intraspecific density dependence in larval development of the Crawfish Frog, Lithobates areolatus. Herpetological Review 43:36–38. Google Scholar

56.

Wright, A. H., and A. A. Wright. 1933. Handbook of Frogs and Toads of the United States and Canada. Comstock Publishing Associates, Ithaca, New York. Google Scholar

57.

Wright, H. P., and G. S. Myers. 1927. Rana areolata at Bloomington, Indiana. Copeia 1927:173–175. Google Scholar

58.

Yagi, K. T., and D. M. Green. 2018. Post-metamorphic carryover effects in a complex life history: behavior and growth and two life stages in an amphibian, Anaxyrus fowleri. Copeia 106:77–85. Google Scholar
© 2020 by the American Society of Ichthyologists and Herpetologists
Rochelle M. Stiles, Vanessa C. K. Terrell, John C. Maerz, and Michael J. Lannoo "Density-Dependent Fitness Attributes and Carry-Over Effects in Crawfish Frogs (Rana areolata), a Species of Conservation Concern," Copeia 108(3), 443-452, (21 July 2020). https://doi.org/10.1643/CH-19-246
Received: 3 June 2019; Accepted: 18 March 2020; Published: 21 July 2020
Back to Top