Ambystoma Opacum Research Articles

Hatching plasticity is associated with a more advanced stage at hatching in an Ambystoma with terrestrial eggs.

Hatching plasticity allows animals to initiate hatching in response to environmental cues including predation, flooding, and hypoxia. In species with terrestrial eggs but aquatic larvae, hatching plasticity often manifests as extended development of embryos when water is not available. Although these effects are taxonomically widespread, little attention has focused on differences in plasticity across closely related species with terrestrial and aquatic embryos. We propose that the terrestrial embryonic environment favors slower and prolonged development and, consequently, that we should see differences in development between closely related species that differ in where they lay their eggs. We test this hypothesis by comparing embryonic development between two mole salamanders, Ambystoma opacum and A. annulatum. Most Ambystoma lay eggs submerged in ponds but A. opacum lays its eggs on land, where hatching is triggered when eggs are submerged by rising pond levels. Embryos of both species were reared under common laboratory conditions simulating both aquatic and terrestrial nest sites. Consistent with our hypothesis, we found that A. opacum embryos exhibited slower development and took longer to hatch than A. annulatum embryos in both rearing environments. Furthermore, we observed in A. opacum a plasticity in hatching stage that was absent in A. annulatum. Our results indicate that the terrestrial-laying A. opacum has evolved slower and prolonged development relative to its aquatic-laying congener and suggest that embryonic survival in the unpredictable terrestrial environment may be facilitated by developmental plasticity.

Local scale population risk determines the adaptive responses of larval salamanders to predator kairomones

Abstract Organisms can exhibit plasticity in phenotypic responses to environmental stimuli but the response can vary based on the evolutionary history of a population. Numerous studies have found that prey will express adaptive responses when subjected to various predators in experimental settings. For species with a large distribution, such as Ambystoma maculatum (spotted salamander), it is not surprising that phenotypic variance can be high across their range. This variance has been hypothesized to be because of different predator regimes of populations. While A. maculatum preferentially oviposit eggs in low-risk ponds (temporary fishless) to decrease egg and larval mortality, Missouri populations still use high-risk ponds (permanent with fish) for reproduction. Using a series of experiments, we investigated how A. maculatum population risk influenced adaptive responses to both native and novel predator kairomones. For natural predators, we used larval A. opacum (marbled salamanders) and Lepomis macrochirus (bluegill) and for novel predators we used Siren intermedia (lesser siren). We found that larval salamanders generally responded with no differences in morphological traits to all three predators. However, head width was larger for low-risk populations. One hypothesized benefit of larger heads (increased foraging efficiency) was not fully supported. All larval salamanders also increased refuge use with predators but this decreased over time. Our results suggest that predation risk of a population may influence the degree of phenotypic expression in response to larval predators. Overall, local adaptation may dictate the ability of prey to respond to environmental conditions within a life stage.

Predicting prey diversity with multiple predator effects

Emergent multiple predator effects can result in either risk enhancement or reduction for prey populations. However, whether emergent multiple predator effects influence survival of entire prey assemblages is unknown, which could alter community diversity. We tested whether emergent multiple predator effects would lead to changes in survival of multiple prey species that would subsequently alter prey diversity. We experimentally tested whether larval ringed (Ambystoma annulatum) and marbled salamanders (A. opacum) had additive or multiplicative effects on hatchlings of an assemblage of five amphibian prey species using outdoor mesocosms. We also examined this question using data from a previously published experiment. We extended the multiplicative risk model to estimate prey diversity based on individual prey survival probabilities from each mesocosm. We found that predators generally had additive effects on all prey species, as the expected survival of prey was not different from the observed data, indicating a lack of emergent multiple predator effects. Consistent with this result, observed diversity was not different from expected diversity based on our extension of the multiplicative risk model, further indicating predator species had additive effects. Re-analysis of previously published data showed that when emergent multiple predator effects resulted in risk reduction, diversity changed in a manner consistent that outcome. Overall, we found that additive and multiplicative effects of multiple predators contributed to structuring entire prey communities, including influencing community metrics like diversity. Future work should consider further under what conditions diversity enhancement or reduction would occur.

Using mechanistic insights to predict the climate‐induced expansion of a key aquatic predator

AbstractAmeliorating the impacts of climate change on communities requires understanding the mechanisms of change and applying them to predict future responses. One way to prioritize efforts is to identify biotic multipliers, which are species that are sensitive to climate change and disproportionately alter communities. We first evaluate the mechanisms underlying the occupancy dynamics of marbled salamanders, a key predator in temporary ponds in the eastern United States We use long‐term data to evaluate four mechanistic hypotheses proposed to explain occupancy patterns, including autumn flooding, overwintering predation, freezing, and winterkill from oxygen depletion. Results suggest that winterkill and fall flooding best explain marbled salamander occupancy patterns. A field introduction experiment supports the importance of winterkill via hypoxia rather than freezing in determining overwinter survival and rejects dispersal limitation as a mechanism preventing establishment. We build climate‐based correlative models that describe salamander occupancy across ponds and years at two latitudinally divergent sites, a southern and middle site, with and without field‐collected habitat characteristics. Correlative models with climate and habitat variation described occupancy patterns better than climate‐only models for each site, but poorly predicted occupancy patterns at the site not used for model development. We next built hybrid mechanistic metapopulation occupancy models that incorporated flooding and winterkill mechanisms. Although hybrid models did not describe observed site‐specific occupancy dynamics better than correlative models, they better predicted the other site's dynamics, revealing a performance trade‐off between model types. Under future climate scenarios, models predict an increased occupancy of marbled salamanders, especially at the middle site, and expansion at a northern site beyond the northern range boundary. Evidence for the climate sensitivity of marbled salamanders combined with their disproportionate ecological impacts suggests that they might act as biotic multipliers of climate change in temporary ponds. More generally, we predict that top aquatic vertebrate predators will expand into temperate‐boreal lakes as climate change reduces winterkill worldwide. Predaceous species with life histories sensitive to winter temperatures provide good candidates for identifying additional biotic multipliers. Building models that include biological mechanisms for key species such as biotic multipliers could better predict broad changes in communities and design effective conservation actions.

Limited population and community effects of hatching asynchrony in a pond‐breeding salamander

AbstractUnderstanding attributes of phenology beyond the mean date of a life history event, such as variability among individuals within a population, is critical to predict how climate‐induced phenological shifts may alter population dynamics. Identifying how phenological variability impacts organisms is especially needed to better understand how phenological shifts affect trophic dynamics (e.g., shifts in variability of top predators affecting primary production). To better understand the effects of phenological variability on both populations and communities, we examined how variation in egg hatching synchrony of predatory marbled salamanders (Ambystoma opacum) impacted intraspecific interactions at the larval stage, ultimately affecting demographic traits and survival through metamorphosis. We also examined how hatching synchrony affected overall trophic dynamics (e.g., primary consumers and producers) in pond food webs. We experimentally manipulated the degree of hatching synchrony of embryonicA. opacumand subsequently reared larvae in outdoor mesocosms. We monitored demographic traits such as larval growth, size at and time to metamorphosis, and survival. To assess trophic dynamics, we monitored zooplankton abundance and phytoplankton biomass during the experiment. Larvae exhibited greater variability in body size in medium and low hatching synchrony treatments compared to high synchrony treatments. Larval body size variation diminished over time to ultimately result in no differences in most life history traits at metamorphosis or survival among hatching synchrony treatments. We also found no differences among treatments in zooplankton abundance or phytoplankton biomass, likely because of minimal variation inA. opacumsurvival among treatments that would induce top‐down changes. Overall, we found that phenological variation may be context dependent in its influence on demography and overall community structure. Because of concerns for how phenological shifts will affect species interactions, greater scrutiny into conditions that would promote changes in population and community dynamics is needed.

Comparing Genetic and Field-Based Estimates of Population Connectivity in Marbled Salamanders, Ambystoma opacum

Estimating connectivity is key for maintaining population viability for pond-breeding amphibians, especially in areas where habitat alterations occur. Here, we used genetic data (microsatellites) to estimate connectivity of marbled salamanders, Ambystoma opacum, among three focal ponds and compared it to field data (capture-mark-recapture estimates) of movement among the same ponds. In addition, we derived least-cost dispersal paths from genetic data and compared them to field connectivity estimates. We found that genetic and field estimates of dispersal were generally congruent, but field-based paths were more complex than genetic-based paths. While both methods complement each other in identifying important source-sink metapopulation dynamics to inform efficient conservation management plans, field data provide a more biologically accurate understanding of the spatial movement of individual marbled salamanders.

Salamander Demography at Isolated Wetlands within Mature and Regenerating Forests

Geographically isolated wetland and surrounding landscape features affect the ecology and life history of amphibian species. We used multistate mark recapture methods and data from over 30,000 captures of adult Ambystoma opacum to explore how survival, breeding, and movement probabilities differed among wetlands surrounded by regenerating 20-year-old clearcuts and mature 100-year-old forest stands. Survival varied among ponds and years but did not differ between regenerating and mature forest habitats. Both sexes at all ponds incurred dramatic mortality during the non-breeding season of a drought year (2001–2002). Females that skipped one or more breeding opportunities had higher breeding probabilities the following year than did successive breeders. Females exiting into regenerating forests had lower breeding probabilities at two of the three ponds. Breeding salamanders tended to make local movements from regenerating to mature forests, particularly when exiting the pond basin. Landscape movements between ponds were generally low, with few individuals moving from mature to regenerating forest habitats. We conclude that clearcuts continue to negatively impact some demographic parameters of salamanders 20 years post-cutting, but other environmental factors may mitigate these effects, and that populations are probably capable of complete recovery, particularly if some mature forest is retained.

Multiple-Predator Effects and Functional Redundancy of Pond Predators

Functional redundancy occurs when different predator species have similar effects on the diversity, abundance, and composition of a prey community. When multiple predators coexist, their interactions can alter prey survival and ultimately diversity through emergent multiple-predator effects (MPEs). MPEs can be exacerbated by differences in predator behavior; however, little is known about the magnitude of MPEs when predators compete for the same prey or have similar behavioral traits. To understand functional redundancy and the interactions of multiple predators in similar niches, as well as their impacts on the composition of a prey community, we conducted two experiments using two ambystomatid salamander predators (Ambystoma opacum and Ambystoma annulatum). We exposed a suite of tadpole prey (Anaxyrus americanus, Rana sphenocephala, Rana sylvatica, Pseudacris feriarum, and Pseudacris crucifer) to different experimental food webs in a mesocosm experiment, and a single prey species (R. sphenocephala) in microcosm experiment with substitutive and additive designs to test for MPEs. We found no evidence for functional redundancy between the two predators: A. annulatum selectively preyed on competitively dominant prey species (A. americanus) and did not alter community diversity. Ambystoma opacum decreased prey diversity relative to the control because of natural phenological mismatch with A. americanus. Interactions between the two predators (e.g., predator inference) were lacking, indicating that predation risk from each predator was independent for the one prey species we tested. A better understanding of community and ecosystem-level effects by A. annulatum can inform future conservation efforts and management decisions regarding this endemic species.

The effects of intraguild interactions (or lack of) on prey diversity in experimental ponds food webs

The interactions between top- and intermediate-level predators can have significant effects on community structure of food webs. The direct and indirect effects of top predators not only affect prey communities both negatively (e.g. predation) and positively (e.g. thinning effects), but can also impact intermediate-level predators (via intraguild predation [IGP]). These IGP interactions can have especially far-reaching implications when predators of increased ecological importance, such as keystone predators, are involved. Yet, the strength of these interactions are often context-dependent with relationships that are highly variable, and little is known about the factors that shape the outcomes. We conducted a mesocosm experiment to test for the potential IGP interactions among top (invertebrate) and intermediate (amphibian) predators and their effects on a diverse prey community. We used larval Boyeria vinosa (Fawn Darner), a rarely studied member of the family Aeshnidae, and Ambystoma opacum (Marbled Salamander), a known keystone predator in parts of its range, in a mesocosm study with an additive experimental design. Each mesocosm received an equal assemblage of five anuran prey species. Total prey survival was lowest when both predators were present, but A. opacum overwhelmingly reduced the survival of three prey species, lowering species diversity in all food webs when present. Species diversity was highest in the B. vinosa food web, with one prey species (Anaxyrus americanus) emerging smaller and sooner than from other food webs. We did not detect any interactions between predator species, indicating that B. vinosa may not be as voracious as other members of Aeshnidae. Our findings suggest that density mediated IGP interactions are non-existent among our focal species. Ultimately, we also demonstrate the importance of examining variability in keystone predators and how abundant, but less-known predator species affect community structure.

Asymmetric density‐dependent competition and predation between larval salamanders

Abstract Body size differences among consumers often lead to asymmetric interactions, with larger individuals typically being stronger competitors and/or predators on small individuals. These types of interaction are particularly exemplified in freshwater pond communities, where substantial size variation exists both within and among species of top consumers. We investigated whether density dependence can modify the outcome of size‐structured interactions between larval stages of two pond‐breeding salamanders, Ambystoma annulatum and Ambystoma opacum. Size structure exists in populations of these species due to variation in the timing of breeding, which we hypothesised would amplify predation rates and competitive asymmetries from the early‐arriving species (A. annulatum) on the later‐arriving species (A. opacum). We manipulated the relative densities of both A. annulatum and A. opacum in outdoor mesocosms. We maintained the experiment through metamorphosis, and analysed size at metamorphosis, larval period length and survival of each species. Ambystoma annulatum imparted a strong density‐dependent effect on A. opacum through a combination of predation and competition. Survival of A. opacum was negatively related to the density of A. annulatum. For the A. opacum that survived, body size was reduced and larval period lengthened at higher A. annulatum densities, indicative of interspecific competition that was partly explained by resource pre‐emption. In contrast, A. annulatum was only affected by intraspecific density‐dependent competition. Our results suggest that density‐dependent effects reinforce asymmetric interactions among larval salamanders. However, the intensity of the asymmetric interactions is mediated by the arrival time and size of conspecifics. Specifically, earlier‐arriving species can negatively affect the later‐arriving species via size‐mediated predation and competition. The interactive effects of density dependence and arrival time of community members are probably a common mechanism generating size variability in ecological communities. Yet, most studies only evaluate one mechanism or the other. By interweaving these two processes, our work displays the importance of understanding context‐dependence in species interactions.

Modern-day record of the Marbled Salamander, Ambystoma opacum (Gravenhorst 1807), from Pinellas County, Florida

Modern-day record of the Marbled Salamander, Ambystoma opacum (Gravenhorst 1807), from Pinellas County, Florida

Macrogeographic variation of a pond predator's top‐down effects in a common garden environment

AbstractGeographic variation in species behavior and life history has been well documented in biology. Species with wide geographic distributions (i.e., across a continent) but small home ranges (i.e., <1 km2) likely experience wide variability in abiotic environments across the entirety of their range, possibly exhibiting strong local adaptation. Understanding variation across a large geographic scale is especially important when considering species that have strong ecological importance, such as keystone species. Yet, few studies have compared the potential cascading ecological effects of a predator with a keystone role in at least part of its range. To understand how keystone ability in pond food webs can vary across a large geographic range, we conducted an artificial pond experiment with a known keystone predator in at least part of its range, the marbled salamander (Ambystoma opacum). To do so, we collected size‐matched salamander larvae from three geographically distant populations (>650 km apart) in Ohio, Mississippi, and North Carolina and placed them in mesocosms with a suite of spring breeding amphibian prey species. We observed differential survival of some prey species leading to differences in spring‐breeding amphibian diversity among the three predator populations, indicating that keystone predation may vary at a geographic scale. Prey diversity was lowest with predators from northern (Ohio) populations of salamanders. Further understanding of large‐scale variability in ecologically important predators and the potential effects of translocating wide‐ranging ambystomatid species is needed to direct future conservation efforts and preserve biodiversity.

Functional Responses of Larval Marbled Salamanders (Ambystoma opacum) and Adult Lesser Sirens (Siren intermedia) on Anuran Tadpole Prey

Predation can have strong effects on the structure of pond-breeding amphibian communities. Many factors can influence the outcome of predator–prey interactions, including differences in densities, identities, and body sizes of both predator and prey. These different mediating factors can impart synergistic impacts on predation rates, though distinguishing such interactions among multiple factors are underexplored. We examined whether different body sizes of two predators, larval Marbled Salamanders (Ambystoma opacum) and adult Lesser Sirens (Siren intermedia), varied in their ability to forage on larval anurans across a range of prey densities. We specifically tested whether attack rates and handling times, the two main parameters of functional response models, varied across three size classes in both predator species. We found that larval Marbled Salamanders exhibited a Type II (saturating) functional response and that larger individuals had higher attack rates and shorter handling times, resulting in greater prey mortality at higher prey densities with larger predators. In contrast, Lesser Sirens were largely ineffective predators despite being an order of magnitude larger in body size than Marbled Salamanders; tadpole mortality was largely unrelated to their own density. Predator body size was a significant predictor of prey mortality for both predator species. Overall, our study shows that species identity could be as important as predator body size when predicting the outcomes of predator–prey interactions.

Multiscale resistant kernel surfaces derived from inferred gene flow: An application with vernal pool breeding salamanders.

The importance of assessing spatial data at multiple scales when modelling species-environment relationships has been highlighted by several empirical studies. However, no landscape genetics studies have optimized landscape resistance surfaces by evaluating relevant spatial predictors at multiple spatial scales. Here, we model multiscale/layer landscape resistance surfaces to estimate resistance to inferred gene flow for two vernal pool breeding salamander species, spotted (Ambystoma maculatum) and marbled (A. opacum) salamanders. Multiscale resistance surface models outperformed spatial layers modelled at their original spatial scale. A resistance surface with forest land cover at a 500-m Gaussian kernel bandwidth and normalized vegetation index at a 100-m Gaussian kernel bandwidth was the top optimized resistance surface for A. maculatum, while a resistance surface with traffic rate and topographic curvature, both at a 500-m Gaussian kernel bandwidth, was the top optimized resistance surface for A. opacum. Species-specific resistant kernels were fit at all vernal pools in our study area with the optimized multiscale/layer resistance surface controlling kernel spread. Vernal pools were then evaluated and scored based on surrounding upland habitat (local score) and connectivity with other vernal pools on the landscape, with resistant kernels driving vernal pool connectivity scores. As expected, vernal pools that scored highest were in areas within forested habitats and with high vernal pool densities and low species-specific landscape resistance. Our findings highlight the success of using a novel analytical approach in a multiscale framework with applications beyond vernal pool amphibian conservation.

Understanding variation in salamander ionomes: A nutrient balance approach

AbstractEcological stoichiometry uses information on a few key biological elements (C, N and P) to explain complex ecological patterns. Although factors driving variation in these elements are well established, expanding stoichiometric principles to explore dynamics of the many other essential elements comprising biological tissues (i.e. the ionome) is needed to determine their metabolic relationships and better understand biological control of elemental flows through ecosystems.In this paper, we report observations of ionomic variation in two species of salamander (Ambystoma opacumandA. talpoideum) across ontogenic stages using specimens from biological collections of two wetlands sampled over a 30‐year period. This unique data set allowed us to explore the extent of ionomic variation between species, among ontogenic stages, between sites and through time.We found species‐ and, to a lesser extent, site‐specific differences in C, N and P along with 13 other elements forming salamander ionomes but saw no evidence of temporal changes. Salamander ionomic composition was most strongly related to ontogeny with relatively higher concentrations of many elements in adult males (i.e. Ca, P, S, Mg, Zn and Cu) compared to metamorphic juveniles, which had greater amounts of C, Fe and Mn.In addition to patterns of individual elements, covariance among elements was used to construct multi‐elemental nutrient balances, which revealed differences in salamander elemental composition between species and sites and changes in elemental proportions across ontogenic development. These multi‐elemental balances distinguished among species‐site‐ontogenic stage groups better than using only C, N and P.Overall, this study highlights the responsiveness of consumer ionomes to life‐history and environmental variation while reflecting underlying relationships among elements tied to biological function. As such, ionomic studies can provide important insights into factors shaping consumer elemental composition and for predicting how these changes might affect higher‐order ecological processes.

Microgeographic Adaptation of Wood Frog Tadpoles to an Apex Predator

Rapid adaptation of defenses can alter ecological dynamics following introduction of a new predator. We tested for local adaptation in Wood Frog (Rana sylvatica) populations that face varying selection from an apex predator, the Marbled Salamander (Ambystoma opacum), which is expanding its distribution in the study region. We performed a reciprocal transplant experiment with Wood Frog eggs and tadpoles and tested survival of tadpoles when exposed to Marbled Salamander larvae in experimental predation trials. We also evaluated life history, behavioral, and morphological trait variation with respect to origin and transplant environments. We found that tadpoles from populations exposed to high risk from Marbled Salamanders survived better when raised in high-risk environments than tadpoles from low-risk populations. However, tadpoles from high-risk environments experienced lower survival than those from low-risk environments when raised in low-risk environments. Development rate, activity, and morphology dif...

Salamanders on the bench – A biocompatibility study of salamander skin secretions in cell cultures

Salamanders have evolved a wide variety of antipredator mechanisms and behavior patterns, including toxins and noxious or adhesive skin secretions. The high bonding strength of the natural bioadhesives makes these substances interesting for biomimetic research and applications in industrial and medical sectors. Secretions of toxic species may help to understand the direct effect of harmful substances on the cellular level.In the present study, the biocompatibility of adhesive secretions from four salamander species (Plethodon shermani, Plethodon glutinosus, Ambystoma maculatum, Ambystoma opacum) were analyzed using the MTT assay in cell culture and evaluated against toxic secretions of Pleurodeles waltl, Triturus carnifex, Pseudotriton ruber, Tylototriton verrucosus, and Salamandra salamandra. Their effect on cells was tested in direct contact (direct culture) or under the influence of the extract (indirect exposure) in accordance with the protocol of the international standard norm ISO 10993-5. Human dermal fibroblasts (NHDF), umbilical vein endothelial cells (HUVEC), and articular chondrocytes (HAC), as well as the cell lines C2C12 and L929 were used in both culture types.While the adhesive secretions from Plethodon shermani are cytocompatible and those of Ambystoma opacum are even advantageous, those of Plethodon glutinosus and Ambystoma maculatum appear to be cytotoxic to NDHF and HUVEC. Toxic secretions from Salamandra salamandra exhibited harmful effects on all cell types. Pseudotriton ruber and Triturus carnifex secretions affected certain cell types marginally; those from Pleurodeles waltl and Tylototriton verrucosus were generally well tolerated.The study shows for the first time the effect of salamander secretions on the viability of different cell types in culture. Two adhesive secretions appeared to be cell compatible and are therefore promising candidates for future investigations in the field of medical bioadhesives. Among the toxic secretions tested, only two of the five had a harmful effect on cells, indicating different cell toxicity mechanisms.

Multiple Paternity Benefits Female Marbled Salamanders by Increasing Survival of Progeny to Metamorphosis

AbstractMultiple paternity occurs in most species and animal groups that have been studied. Because mating involves fitness costs to individual females, theory predicts that polyandrous females gain greater fitness benefits than costs, allowing the behavior to be maintained. Genetic, rather than material, benefits often occur in species where males provide females with little more than sperm and seminal fluid. We compared fitness correlates of single‐ and double‐sire clutches from female marbled salamanders (Ambystoma opacum) at the egg, hatchling, and metamorph stages of offspring development. Because clutches were collected from experimental breeding groups, strict paternity exclusion of offspring using microsatellite data allowed us to categorize each clutch as having either one or two fathers. Early offspring viability and size of hatchlings were not different between single‐ and multiple‐paternity clutches. Larvae from the two clutch types were allowed to develop together in field enclosures until metamorphosis. Although there was no difference in size at metamorphosis, survival to metamorphosis was significantly higher in multiple‐paternity clutches (44% vs. 40%) suggesting a benefit for females. The results were consistent with genetic benefits, although maternal effects could not be ruled out. The data did not support predictions of the genetic bet‐hedging and good sperm hypotheses for genetic benefits of polyandry.

Agonistic Behavior and Resource Defense among Sympatric Juvenile Pond-Breeding Salamanders

Abstract Interference competition is frequently observed in sexually mature adults as they defend breeding territories. However, it remains unclear in many taxa whether juveniles respond aggressively to other juveniles or if they defend resources. To test whether postmetamorphic juveniles of a pond-breeding amphibian were aggressive towards other juveniles or were defending resources, we staged terrestrial encounters between three species of sympatric Ambystoma salamanders. We observed biting and other aggressive behaviors by juvenile Spotted (Ambystoma maculatum) and Marbled Salamanders (Ambystoma opacum). However, we did not observe aggressive behaviors by Ringed Salamanders (Ambystoma annulatum). In addition to species-level variations in aggression, these three species also differed in whether aggression was targeted primarily intra- or interspecifically. This study suggests that juveniles of pond-breeding amphibians of some species may defend essential habitat with agonistic behavior.

Differences in Larval Allometry among Three Ambystomatid Salamanders

Abstract Ecological and evolutionary processes commonly result in morphological variation among larval amphibians. Variation in head shape plays a critical role in both food capture and predation risk in gape-limited salamanders, yet in situ studies of head shape variability are rare outside of cannibal morph assessments. We examined allometry differences in larval head width (HW) and snout–vent length (SVL) among three sympatric species of ambystomatid salamanders from 166 ponds in Missouri, USA: Ringed Salamander (Ambystoma annulatum), Marbled Salamander (A. opacum), and Spotted Salamander (A. maculatum). We also tested whether several abiotic and biotic factors would predict HW after accounting for SVL. We found that larval HW and SVL were strongly correlated for all species but that the strength of this relationship varied among species. For early-stage larvae, Marbled Salamanders showed isometric scaling relationships, whereas both Spotted Salamanders and Ringed Salamanders were allometric. For late-...