The behavioral phenotypes of hybrids vary in degree of similarity to their parent species. Unisexual salamanders (Ambystoma laterale sp.), the result of ancient hybridization, contain nuclear DNA of multiple sperm-host species whose habitat preferences differ from one another. We radio tracked unisexual salamanders from four vernal pools to quantify migration distances and post-breeding habitat selection and compared these to published accounts for Blue-Spotted Salamanders (A. laterale) and Jefferson Salamanders (Ambystoma jeffersonianum). Unisexual salamanders used sites with higher numbers of small mammal burrows, lower substrate temperatures, and lower cover by forest floor vegetation than available sites, similar to the sperm-hosts. Unisexual salamanders also migrated distances within the range reported for these sperm-hosts. Even so, individual migration distances were context specific. We implore managers to use caution when designating management zones around breeding pools by considering that some populations may move farther than those reported in published accounts.
Unisexual taxa are almost entirely female and reproduce either without sperm or with sperm contributed by males of bisexual species (with separate males and females; Dawley, 1989). Unisexual teleost fish, unisexual amphibians, and many unisexual lizards are the result of past hybridizations between bisexual species (Neaves and Baumann, 2011). Genetic variation influences habitat selection and can result in some hybrids being less selective than either parent species, easing their fit into hybrid zones between allopatric parent populations, whereas others remain sympatric and compete with parent species (Jaenike and Holt, 1991; Saino, 1992; Wood et al., 2016). Unisexuals, likewise, may use a variety of niches, allowing some to thrive in different habitats, persist in changing environments, and reduce competition with parent species (Bullini, 1994; Mee and Rowe, 2010).
The Blue-Spotted Salamander Complex is the result of a 5-million-year-old hybridization event creating a lineage of modern salamanders carrying combinations of the genomes of Blue-Spotted Salamanders (Ambystoma laterale), Jefferson Salamanders (Ambystoma jeffersonianum), Tiger Salamanders (Ambystoma tigrinum), Small-Mouthed Salamanders (Ambystoma texanum), and, rarely, Streamside Salamanders (Ambystoma barbouri; Uzzell, 1964; Morris and Brandon, 1984; Bogart et al., 2009; Bi and Bogart, 2010). Unisexual salamanders have nuclear DNA from two or more of these species and are almost always polyploid (Lowcock and Murphy, 1991; Bogart and Klemens, 1997). We use the convention of abbreviating the genetic composition (genomotype rather than genotype) of individuals by how many replicates of each genome they contain (e.g., LL for A. laterale, LLJ for A. (2) laterale–jeffersonianum, and LLLJ for A. (3) laterale–jeffersonianum, Lowcock et al. 1987). Although size and/or mass can give some indication of genomotype and our unpublished data indicate that salamanders over 7 g in our area are unisexuals, they are similar in appearance to sperm-hosts and genetic methods often are needed for identification.
Unisexual salamander habitat studies have been limited to establishing geographic and climatic niche (Greenwald et al., 2016), examining habitat relations to the subcanopy (Belasen et al., 2013), and modeling breeding site characteristics (Hoffmann, 2017). The post-breeding habitat selection and migration distances of unisexual salamanders are critical to informing management decisions but have not been quantified. In Maine, unisexual salamanders contain more A. laterale than A. jeffersonianum DNA (e.g., are LLJ and LLLJ; Bogart and Klemens, 1997); thus, we hypothesized that their terrestrial habitat preferences and migration distances would be similar to published accounts of A. laterale. Both sperm-hosts occupy burrows in forests and migrate to seasonal wetlands to breed (Petranka, 1998); however, other aspects of habitat use vary. They partition habitat by altitude with A. jeffersonianum typically in well-drained uplands and A. laterale in lowlands (Nyman et al., 1988; Downs, 1989; Klemens, 1993). Ambystoma jeffersonianum can migrate farther than A. laterale can (Williams, 1973; Douglas and Monroe, 1981; Ryan and Calhoun, 2014). Ambystoma jeffersonianum have been documented using forested landscapes with low disturbance (Porej et al., 2004; Rubbo and Kiesecker, 2005; Greenwald et al., 2016), whereas some researchers have documented A. laterale in more open habitat that may have more anthropogenic disturbance (Weller et al., 1978; Downs, 1989; Klemens, 1993; Windmiller et al., 2008).
Our goal was to understand the post-breeding habitat selection of unisexual salamanders in comparison to published accounts of sperm-host species. Specifically, we 1) quantified the emigration distances of LLJ and LLLJ, 2) examined microhabitat selection in late-spring and summer, and 3) compared these to published studies of A. jeffersonianum and A. laterale. We expected the unisexuals to be similar to A. laterale by migrating about 70 m and selecting forested habitat with moist soil, deep leaf litter, and woody debris (Ryan and Calhoun, 2014).
Materials and Methods
Site Selection and Capture
We encircled four vernal pools in Penobscot County, Maine, with drift fence for a related study (Hoffmann, 2017). Pools 1 and 2 were located in Old Town on a parcel managed for forestry by the University of Maine and abutting the Stillwater River (>200 m across). Pools 3 and 4 were located on privately owned residential parcels in the town of Orono. Their forest matrix was penetrated by residential neighborhoods and fields. We constructed the drift fences out of silt fence buried 20 cm into the ground with a pair of #10 aluminum cans every 5 m to act as pitfall traps (Shoop, 1965). We checked the traps every morning from April to September and alternating mornings in October and November (Hoffmann, 2017).
Surgery and Genetic Testing
We implanted 14 unisexual salamanders with radio transmitters in 2013 following the methods of Madison et al. (2010) and implanted another 29 in 2014. We monitored only large emigrating unisexual salamanders because A. laterale (LL) are generally too small to safely implant (Ryan and Calhoun, 2014), and we were outside the geographic range of A. jeffersonianum. We anesthetized salamanders in 3.1 mM tricaine methane sulfonate (MS-222) neutralized to a pH of 7.0 with aqueous NaOH until loss of pain response (toe pinching). We used surgical scissors to remove a 0.5 by 0.3 cm tissue sample (Nöel et al., 2011) from the tip of the tail, which we shipped in 70% ethanol to the University of Guelph to determine genomotype using microsatellite DNA analyses at six loci (AjeD75, AjeD94, AjeD283, AjeD346, AjeD378, and AjeD422; Julian et al., 2003). Our microsatellite DNA methods are described elsewhere (Bogart et al., 2007, 2009). We inserted ATS A2415 transmitters (0.33 g, Advanced Telemetry Systems, Isanti, MN) with the antennas removed and PIT tags (0.09 g, HPT12, 134.2kHz ISO FDXB tag; Biomark, Boise, ID) into coelomic cavities (Ryan and Calhoun, 2014) using 10-mm longitudinal incisions in the left ventrolateral abdominal walls (Ryan and Calhoun, 2014), then closed the wounds with absorbable sutures (Model PDS Plus, RB-1 taper, Size 5-0, Ethicon Inc, Somerville, NJ). Body mass ranged from 7.74–13.38 g, such that the transmitters represented ≤4.26% of body mass. The salamanders recovered overnight and were released under wet leaves outside the drift fence in 2013. After realizing that we might bias movement by placing the animals on one side of the drift fence, we released the salamanders back into the vernal pools in 2014. We extended our 2013 telemetry season from 42 days (battery life of one transmitter) to 92 days by replacing transmitters in six animals (after McDonough and Paton, 2007; Titus et al., 2014). Unfortunately, the skin was weak at the site of the original incision, and we found two animals with an open incision 7 and 12 days after the second reimplant surgery. Therefore, in 2014, we tracked each salamander for the life of only one transmitter.
We relocated 13 salamanders daily in 2013. In 2014, we tracked 27 salamanders and, because of the increase in sample size, we tracked each individual only once every three days. We used a Model R-1000 receiver (Advanced Telemetry Systems) and Model RA-2AK VHF antennae (Telonics, Inc. Mesa, AZ) for direct overhead localization (10-cm accuracy). In 2014 when a transmitter expired, we scanned an ∼20 m radius using a Reader 03 PIT tag reader (West Fork Environmental, Turnwater, Washington, USA) with a custom-built antenna (Blomquist et al., 2008) before dismissing a salamander as lost, extending tracking for 11 salamanders for a mean of 15 days (range: 3–36 days). We recorded locations with a handheld GPS (GPSMAP 62stc and eTrex 10, Garmin International, Inc. Olathe, KS). We considered an animal not to have moved if it was within 3 m on subsequent visits in 2013 and within 0.5 m in 2014 attributable to new information in Ryan and Calhoun (2014), who found that A. laterale was selecting habitat at a small scale.
We measured microhabitat variables at paired used and available plots in succession such that meteorological conditions and vegetation phenology were comparable within pairs. We measured variables at two scales: when the salamanders moved ≥3 m (larger scale), we compared 3-m radius used and available plots (Faccio, 2003); and for movements < 3 m (smaller scale), we compared 0.5-m radius used and available plots (2014 only, based on Ryan and Calhoun, 2014).Three available plots evenly surrounded each used plot with one available plot located along a random bearing and the others 120° to either side (Fig. 1). The center of available plots for the smaller movement was always located 6 m from the center of the used plot. Distances between the centers of plots at the larger scale varied by week and were determined by the median distance between sequential locations during the previous two weeks in 2013 (range: 56–6 m), except for during the initial two weeks when we used the first weeks' median distance. We used the distances measured in 2013 to space plots in 2014.
We measured 22 microhabitat variables within plots based on studies of A. jeffersonianum, A. laterale, and other amphibians (Faccio, 2003; Rittenhouse and Semlitsch, 2007; Ryan and Calhoun, 2014; Groff et al., 2017). We recorded land use as forest, yard/field, or wetland. We visually estimated the percent cover of bare soil, all leaf litter, coniferous leaf litter, water, lawn/hay, moss, rock, vegetation < 1 m tall, vegetation between 1 and 3 m tall, and vegetation > 3 m tall but < 10 cm DBH (diameter at breast height), allowing cover amounts to sum to >100%. We counted the number of stumps. We used a spherical convex densitometer (Model-A, Forestry Suppliers Inc. Jackson, MS) to quantify angular canopy cover because we assumed light interception by tree canopy surrounding the plot rather than just directly above it may be biologically important as a cue for cover (Nuttle, 1997). We measured leaf litter depth to the nearest 0.5 cm, soil moisture (FieldScout TDR200; Spectrum Technologies, Inc. Aurora, IL), and soil temperature in °C (Model 9841; Taylor Precision Products, Oak Brook, IL) near the center of each plot. We measured the diameter and length of coarse woody debris (>10 cm) to calculate the total area covered (cm2) and recorded the maximum decay stage for 3-m plots only (Monti, 1997). Once a salamander had vacated the plot, we brushed away the leaf litter and recorded the number of horizontal and vertical small mammal burrow openings (after Faccio, 2003).
We plotted locations in ArcGIS 10.3 (ERSI, Redlands, CA) and used the “adehabitatLT” package (Calenge, 2006) in Program R (R Core Team, 2016) to determine step length, cumulative distance, and maximum straight line distance from the vernal pool for each salamander. Unless otherwise noted, summary statistics are reported as mean ± SD. We used a Kruskal–Wallace test to determine whether maximum straight-line distances varied by pool and Spearnman's rank to determine whether distance was correlated with salamander mass (Jehle and Arntzen, 2000). We found the farthest distances from the pool for each unisexual salamander and compared the mean of these distances and 95% life zones for each pool to those calculated from published data for the parent species. Faccio (2003) tracked six A. jeffersonianum and followed Semlitsch (1998) in using a 95% confidence interval to determine the radius of a life zone that would include 95% of a study population's distances, suggesting managers keep habitat within this buffer intact to conserve populations. In contrast, we argue the use of confidence intervals to establish 95% life zones is incorrect and that quantiles are more appropriate. Confidence intervals are intended to give precision of estimation of the population mean (μ) such that the maximum of this interval is a possible population mean (P[lower CI ≤ μ ≤ upper CI] = 95%), not the area that includes 95% of the salamanders. We think Semlitsch may have intended to use quantiles, as we simply use the standard deviation rather than the standard error, such that the distance is not inflated by low sample size. Ryan and Calhoun (2014) sorted the distances traveled by A. laterale in ascending order and determined which distance included 95%. We recalculated life zones using the 95% quantile for these published data, the radio-isotope tracked A. jeffersonianum of Williams (1973), and our own observations using t-scores. Neither of these A. jeffersonianum populations has been genomotyped. However, Williams' animals were outside the geographic range of unisexual salamanders (Petranka, 1998; Charney, 2011) and Faccio's salamanders were likely mostly JJ based on sex ratios of the breeding population (S. Faccio, personal communication). Also, half of Faccio's telemetered salamanders were males, and, therefore, almost certainly JJ. Ryan and Calhoun's (2014) A. laterale population was known to contain no unisexual salamanders.
We used conditional logistic regression models to examine variables related to microhabitat selection by comparing plots used by individuals to their own available plots at a given time (i.e., study design IV; Erickson et al., 2001). We used the themes of land use, shelter, ground cover, microclimate, and vegetation to develop 22 specific a priori models (Table 1) about how salamanders may be using the landscape (i.e., if salamanders seek rotten wood when choosing where to settle rather than just any cover object, then a model containing just stumps and coarse woody debris would rank better than one also containing rocks) and 3 composite models (literature models in Table 1) based on the top models to predict plots selected by A. laterale from Ryan and Calhoun (2014). These models were not exclusive, such that variables in our shelter models may also appear in our ground cover models. We Z-standardized continuous variables, checked for collinearity (Pearson r ≥ 0.7) and ranked all models seperately for the 3-m plot and 0.5-m plot scales. We weighed each plot by the proportion of days the salamander spent there, such that the experimental units were animals rather than relocations (Aebischer et al., 1993; Thomas and Taylor, 2006). We used the “survival” package (Therneau, 2015) in program R to conduct conditional logistic regression to examine selection of microhabitat features. We used Akaike Information Criterion adjusted for small sample size (AICc; Burnham and Anderson, 2002) to rank models separately for each scale with package “AICcmodavg” (Mazerolle, 2016). We used the “support.Ces” package (Aizaki, 2012) to determine McFadden's R2 (ρ). Rather than make inference based solely on the fit of individual models, we found the model averaged estimates of the odds ratios of the variables from models with a cumulative model weight of ≤0.9 and considered these variables important if the confidence intervals of their odds ratios did not include 1.
Conditional logistic regression models used to compare each 3-m and 0.5-m plot used by unisexual salamanders to three paired random plots in central Maine. The variables CWD (coarse woody debris) area, CWD decay stage, and forest were not included in models for 0.5-m plots because of lack of field measurements and lack of convergence caused by low variation.
Of the 42 unisexual salamanders sampled in this study, 39 were LLJ, 2 were LLLJ, and 1 was unidentified but greater than the maximum observed mass of local A. laterale (LL; unpubl. data). Seven clones were represented in this sample and included 21 salamanders (Appendix 1). Body mass did not vary by pool (Kruskal-Wallace χ2 = 2.91, df = 3, P = 0.406).
We tracked unisexual salamanders from 5 to 94 days (mean = 51 days), during which they moved an average straight-line distance of 172 m (range = 6–403 m) from the breeding pool (Table 1, Appendix 1). In 2013, the mean cumulative distance was 191 ± 76 m (range = 6–410 m). The mean cumulative distance in 2014 was 209 ± 140 m (range = 47–463 m). Maximum straight-line distance from the pool was not related to the body mass of the salamander (r = 0.090, P = 0.581) but varied by pool (Kruskal-Wallis χ2 = 18.45, df = 3, P = 0.004). Salamanders in Old Town remained closer to the pools (mean = 112 ± 44 m and 36 ± 19 m) than did salamanders in Orono (mean = 244 ± 76 and 214 ± 113 m). Distance to the pool generally reached an asymptote within a week, as salamanders made large initial migrations with few short subsequent movements (Fig. 2).
Ninety-five percent life zones for unisexual salamanders in our study also varied by pool (195 m, 74 m, 383 m, 415 m for Pools 1, 2, 3, and 4, respectively) and extended 362 m with all pools combined (Fig. 3). The 95% quantile for distance traveled by A. jeffersonianum in Indiana, was 478 m (mean = 252 ±136 m, N = 86, based on Williams, 1973), whereas the zone for A. jeffersonianum/unisexual salamanders in Vermont, was 143 m (mean = 92 ± 25, N = 6, based on Faccio, 2003). The life zone for A. laterale was 149 m (mean = 64.9 ± 50.1, N = 43, based on Ryan and Calhoun, 2014). These life zones are slightly smaller than those calculated by the authors (e.g., Faccio reported a 157-m zone based on confidence intervals, Ryan and Calhoun reported a 152-m life zone based on sorting distances by size, and Williams did not calculate a zone).
Most unisexual salamanders (36 of 40) remained within the forest matrix for the entire study, although 12 of these emigrated to post-breeding home ranges that were within about 20 m of forest-lawn or forest-hay field edges (Fig. 4). Seven salamanders occupied swamps dominated by Alder (Alnus incana) and Highbush Blueberry (Vaccinium corymbosum). The remaining four salamanders crossed lawns during emigration and spent the majority of the season underneath buildings (two salamanders under separate sheds with wooden floors and two salamanders under the same garage on a concrete slab).
At the larger movement scale (≥3 m movements), only the “All tunnels” model had substantial support (ΔAICc = 0, Table 2). Seven other models related to shelter, vegetation, and microclimate, including important models for A. laterale based in Ryan and Calhoun (2014) had some support (ΔAICc < 7), though McFadden's adjusted pseudo R2 was low for the five models that did not include “Mammal burrows” as a covariate. We model averaged the β estimates and found the 95% confidence intervals for the odds ratio of each variable that appeared in larger movement scale models with a cumulative 0.9 model weight. Selection of these top models appears to be driven by three important covariates (Table 3), with unisexual salamanders more likely to use plots with more horizontal burrows, lower substrate temperatures, and less vegetation <1 m tall than available plots.
Top-ranked unisexual salamander paired logistic regression models for used and random locations in central Maine. Only models with ΔAICc < 7 are shown. K is the number of parameters, adjusted ρ2 is McFadden's adjusted pseudo R2, and w is the model weight.
Important model averaged parameter estimates (β), odds ratios, and descriptive statistics for unisexual salamander paired logistic regression models for plots of used and random locations in central Maine. Only covariates from models included in 90% of the cumulative model weight and with 95% confidence intervals or odds ratios that did not include one are shown.
At the smaller movement scale (<3 m), both the “All tunnels” model and “Mammal burrows” model had substantial support, and the global shelter and soil moisture models had some support, although all models had low McFadden's adjusted pseudo R2. At the small movement scale, only horizontal burrows had an odds ratio with confidence intervals that did not include 1.
Unisexual salamander post-breeding movement patterns were similar to other ambystomatids, characterized by long movements during emigration over a few nights followed by infrequent and shorter movements in their post-breeding home range (Fig. 2; Williams, 1973; Madison, 1997; Titus et al., 2014). Unisexual salamanders moved as far as 463 m from the pool. Although individuals moved about five times on average during the study, five salamanders (of those tracked > 20 days) moved only once (from the pool to the summer location) and remained within the same 3-m plot for the season. This stationary behavior has been reported for other ambystomatids and also directly observed in Wood Frogs (Lithobates sylvaticus; Douglas and Monroe, 1981; Rittenhouse and Semlitsch, 2007). This may be a “sit-and-wait” predatory strategy.
Mean, median, and 95th percentiles of amphibian migration distances are used to justify the conservation of terrestrial habitat through regulatory or management zones (Semlitsch, 1998); however, these distances have not been widely quantified across and within species, which may be problematic for managers. Managing Orono pools based solely on Old Town Pools could jeopardize the majority of terrestrial habitat. For example, using the 95% life zone from Pool 2 would conserve none of our salamanders at Pool 3 and protect only 12.5% of our salamanders at Pool 4, whereas using Pool 1's life zone would conserve only 25% and 30%.
Unisexual salamanders in our study generally migrated within the range of distances reported for A. laterale and A. jeffersonianum in other studies (Fig. 3). Mean distances from the pool and 95% life zones for unisexual salamanders at three of our four pools were greater than those of A. laterale in Connecticut (Ryan and Calhoun, 2014) and A. jeffersonianum in Vermont (Faccio, 2003), but A. jeffersonianum in Indiana had a larger mean and life zone distances (Williams, 1973). Other references also list mean distances of unisexual salamanders (presumably LLJ) found under cover objects in Michigan as intermediate (110 m) and radioisotope tagged A. jeffersonianum in Kentucky as farther (250 m; Douglas and Monroe, 1981; Belasen et al., 2013). Although we found no relationship between unisexual salamander body mass and distance, many of our populations may have moved farther on average than A. laterale because we selected individuals larger than this parent species. This comparison does not consider variation attributable to geographic location, and we recommend that future work directly compare taxa at the same site.
Ninety percent of tracked unisexual salamanders stayed in the forest, but 13 of 22 salamanders from our Orono pools had post-breeding home ranges near or within residential neighborhoods (i.e., within about 20 m of lawns or fields). We are unsure whether salamanders settled in these areas because they interpreted the neighborhoods as less suitable or whether they sought out locations adjacent to open lands. Ambystomatids are known to cross open areas; however, they also avoid forest edges (deMaynadier and Hunter, 1998; Gibbs, 1998; Regosin et al., 2005; Pittman and Semlitsch, 2013). Forest edges are associated with reduced soil moisture, angular canopy cover, and coarse woody debris and increased forest floor disruption, predation, and temperature (reviewed in Lindenmayer and Fischer, 2006). Pesticide and herbicides may contaminate lawns, but these areas also have high primary production and may have high plant and invertebrate diversity (Falk, 1976; Frankie and Ehler, 1978; McKinney, 2008). Buildings may act as large cover objects to reduce fluctuation in temperature and moisture. We cautiously suggest further study to determine whether neighborhoods are filters, whether salamanders in more urban areas behave similarly, and whether salamanders residing near lawns have lower survival than those in forest interior.
Our top unisexual salamander microhabitat selection models included those based on shelter, vegetation, and microclimate, including a model based on microhabitat selection of A. laterale (from Ryan and Calhoun, 2014; containing soil temperature, leaf litter depth, and soil moisture). Ground cover and land use covariates were not supported, presumably because of the homogeneity of the landscapes in our study area. Only three variables were important: horizontal burrows; forest floor vegetation (herbaceous and woody plants within a meter of the ground); and soil temperature. The most important feature for predicting use by unisexual salamanders was also important in previous studies of parent species and other ambystomatids (Williams, 1973; Douglas and Monroe, 1981; Madison, 1997; Osbourn et al., 2014). Horizontal small mammal burrows were selected by unisexual salamanders both during large movements (migration) and during shorter movements within their post-breeding home range. Horizontal burrows are particularly important and are selected over vertical by both A. jeffersonianum and Ambystoma maculatum in Vermont (Faccio, 2003). Our finding that minimal forest floor vegetation and low temperatures are important may be because shaded areas remain moist and, therefore, indicate conditions conducive to thermoregulation and hydroregulation. Ambystoma jeffersonianum likewise select areas shaded by shrubs (Faccio, 2003). Salamanders, in general, are thought to behaviorally thermoregulate by selecting cool refugia, but rarely are temperature relations observed in the field (Feder and Pough, 1975; Stebbin and Cohen, 1995; Welsh and Lind, 1995). Stebbin and Cohen (1995) also suggest that selection of low temperatures may aid in recovery from high metabolic demands, such as migration and breeding. Our sites were relatively homogenous; hence, we cannot rule out other variables that may be important. Leaf litter, shrubs, logs, soil moisture, and canopy are important to sperm-hosts (Faccio, 2003; Ryan and Calhoun, 2014) but were consistently high in our area.
Our sample of unisexual salamanders selected microhabitat and migrated distances within the known parameters of both sperm-hosts. The overlap in habitat features may allow unisexual salamanders to colonize landscapes wherever sperm-hosts are present; however, unisexuals may require larger forest patches than sperm-hosts if their life zone is larger than the sperm-host's life zone (Mee and Rowe 2010). We suggest future work to track sperm-hosts and unisexual salamanders from the same breeding pool to directly compare habitat selection for differences that might allow coexistence of the taxa.
Unisexual salamanders are generally more abundant than their sympatric sperm host but are unusual among vertebrates in their reproductive system and, therefore, warrant conservation. We recommend maintaining forest conditions that support small mammal populations to provide burrows, avoiding use of lawn chemicals because some salamanders resided near lawns, and we urge further studies to examine the use of rural and suburban/exurban neighborhoods by ambystomatids. We emphasize that migration distances are context specific, and we caution resource managers to be conservative in designating management zones.
D. Harrison and M. Kinnison assisted in study design and early review of this document. Field support was provided by E. D'Urso, A. Feuka, C. Herr, S. Hutchinson, I. Lookabaugh, and K. Sypher. T. Hastings assisted with ArcGIS. We thank B. Cline, L. Groff, and K. Ryan for their influence on our methodology. This work was supported by National Science Foundation Grants 0239915 and EPS-0904155 to the Maine Experimental Program to Stimulate Competitive Research (EPSCoR), the Maine Sustainability Solutions Initiative, and the Maine Agricultural and Forest Experiment Station (MAFES). Fieldwork was conducted under the University of Maine Institutional Animal Care and Use Committee protocol 2013-03-02. This is Maine Agriculture and Forest Experiment Station Publication 3602. This project was supported by the USDA National Institute of Food and Agriculture, Hatch, and McIntire-Stennis project number ME0-21705 through the Maine Agricultural and Forest Experiment Station.
- Aebischer, N. J., P. A. Robinson, and R. E. Kenward. 1993. Compositional analysis of habitat use from animal radio-tracking data. Ecology 74:1313–1325. Google Scholar
- Aizaki, H. 2012. Basic functions for supporting and implementation of choice experiments in R. Journal of Statistical Software 50:1–24. Google Scholar
- Belasen, A., E. Burkell, A. Injaian, K. Li, and I. Perfecto. 2013. Effects of sub-cannopy on haibtat selection in the blue-spotted salamanders (Ambystoma laterale–jeffersonianum unisexual complex). Copeia 2013:254–261. Google Scholar
- Bi, K., and J. P. Bogart. 2010. Time and time again: unisexual salamanders (genus Ambystoma) are the oldest unisexual vertebrates. BMC Evolutionary Biology 10:238. Google Scholar
- Blomquist, S. M., J. D. Zydlewski, and M. L. J. Hunter. 2008. Efficacy of PIT tags for tracking the terrestrial anurans Rana pipiens and Rana sylvatica. Herpetological Review 39:174–179. Google Scholar
- Bogart, J. P., and M. W. Klemens. 1997. Hybrids and genetic interactions of mole salamanders (Ambystoma jeffersonianum and A. laterale) (Amphibia: Caudata) in New York and New England. American Museum Novitates 3218:1–78. Google Scholar
- Bogart, J. P., K. Bi, J. Fu, D. W. A. Noble, and J. Niedzwiecki. 2007. Unisexual salamanders (genus Ambystoma) present a new reproductive mode for eukaryotes. Genome 136:119–136. Google Scholar
- Bogart, J. P., J. Bartoszek, D. W. A. Nöble, and K. Bi. 2009. Sex in unisexual salamanders: discovery of a new sperm donor with ancient affinities. Heredity 103:483–493. Google Scholar
- Bullini, L. 1994. Origion and evolution of animal hybrid species. Trends in Ecology and Evolution 9:422–426. Google Scholar
- Burnham, K. P., and D. R. Anderson. 2002. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. 2nd ed. Springer, USA. Google Scholar
- Calenge, C. 2006. The package “adehabitat” for the R software: a tool for the analysis of space and habitat use by animals. Ecological Modelling 197:516–519. Google Scholar
- Charney, N. D. 2011. Relating hybrid advantage and genome replacement in unisexual salamanders. Evolution 66:1387–1397. Google Scholar
- Dawley, R. M. 1989. An introduction to unisexual vertebrates. Pp. 1–18in R. M. Dawlwyand J. P. Bogart (eds.), Evolution and Ecology of Unisexual Vertebrates. New York State Museum, USA. Google Scholar
- deMaynadier, P. G., and M. L. Hunter. 1998. Effects of silvicultural edges on the distribution and abundance of amphibians in Maine. Conservation Biology 12:340–352. Google Scholar
- Douglas, M. E., and B. L. J. Monroe. 1981. A Comparative study of topographical orientation in Ambystoma (Amphibia: Caudata). Copeia 1981:460–463. Google Scholar
- Downs, F. L. 1989. Ambystoma laterale Hallowell blue-spotted salamander. Pp. 102–107in R. A. Pfingstenand F. L. Downs (eds.), Salamanders of Ohio. Ohio Biological Survey Bulletin 7. Google Scholar
- Erickson, W. P., T. L. McDonald, K. G. Gerow, S. Howlin, and J. W. Kern. 2001. Statistical issues in resource selection studies with radio-marked animals. Pp. 209–242in J. J. Millspaughand J. M. Marzluff (eds.), Radio Tracking and Animal Populations. Academic Press, USA. Google Scholar
- Faccio, S. D. 2003. Postbreeding emigration and habitat use by Jefferson and spotted salamanders in Vermont. Journal of Herpetology 37:479–489. Google Scholar
- Falk, J. H. 1976. Engergetics of a suburban lawn ecosystem. Ecology 57:141–150. Google Scholar
- Feder, M. E., and F. H. Pough. 1975. Temperature selection by the red-backed salamanders, Plethodon c. cinereus (Green) (Caudata: Plethodontidae). Copmarative Biochemistry and Physiology 50:91–98. Google Scholar
- Frankie, G. W., and L. E. Ehler. 1978. Ecology of insects in urban environments. Annual Review of Entomology 23:367–387. Google Scholar
- Gibbs, J. P. 1998. Amphibian movements in response to forest edges, roads, and streambeds in southern New England. Journal of Wildlife Management 62:584–589. Google Scholar
- Groff, L., A. J. K. Calhoun, and C. S. Loftin. 2017. Amphibian terrestrial habitat selection and movenement patterns vary with annual life-hisory period. Canadian Journal of Zoology 95:433–442. Google Scholar
- Hoffmann, K. E. 2017. Breeding Ecology and Habitat Use of Unisexuals Salamanders and Their Sperm-Hosts, Blue-Spotted Salamanders (Ambystoma laterale). Unpubl. Ph.D. diss., Wildlife Ecology, University of Maine, Orono. Google Scholar
- Jaenike, J., and R. D. Holt. 1991. Genetic variation for habitat preference: evidence and explanations. American Naturalist 137:567–590. Google Scholar
- Jehle, R., and J. W. Arntzen. 2000. Post-breeding migrations of newts (Triturus cristatus and T. marmoratus) with contrasting ecological requirements. Journal of Zoology 2000:297–306. Google Scholar
- Julian, S. E., T. L. King, and W. K. Savage. 2003. Novel Jefferson salamander, Ambstoma jeffersonianum, microsatellite DNA markers detect population structure and hybrid complexes. Molecular Ecology Notes 3:95–97. Google Scholar
- Klemens, M. W. 1993. Amphibians and Reptiles of Connecticut and Adjacent Regions. State Geological and Natural History Survey of Connecticut, Bulletin 122. Connecticut Department of Environmental Protection, USA. Google Scholar
- Lindenmayer, D. B., and J. Fischer. 2006. Habitat Fragmentation and Landscape Change. Island Press, USA. Google Scholar
- Lowcock, L. A., and R. W. Murphy. 1991. Pentaploidy in hybrid salamanders demonstrates enhanced tolerance of multiple chromosome sets. Experientia 47:490–493. Google Scholar
- Lowcock, L. A., L. E. Licht, and J. P. Bogart. 1987. Nomenclature in hybrid complexes of Ambystoma (Urodela: Ambystomatidae): no case for the erection of hybrid “ species.”Systematic Zoology 36:328–336. Google Scholar
- Madison, D. M. 1997. The emigration of radio-implanted spotted salamanders, Ambystoma maculatum. Journal of Herpetology 31:542–551. Google Scholar
- Madison, D. M., V. Titus, and V. S. Lamoureux. 2010. Movement patterns and radiotelemetry. Pp. 185–202in C. K. Dodd (ed.), Amphibian Ecology and Conservation: A Handbook of Techniques. Oxford University Press, USA. Google Scholar
- McDonough, C., and P. W. C. Paton. 2007. Salamander dispersal across a forested landscape fragmented by a golf course. Journal of Wildlife Management 71:1163–1169. Google Scholar
- McKinney, M. L. 2008. Effects of urbanization on species richness: a review of plants and animals. Urban Ecosystems 11:161–176. Google Scholar
- Mee, J. A., and L. Rowe. 2010. Distribution of Phoxinus eos, Phoxinus neogaeus, and their asexually-reproducing hybrids (Pisces: Cyprinidae) in Algonquin Provincial Park, Ontario. PLoSOne 5:e13185. Google Scholar
- Monti, L. 1997. Redback Salamander Habitat Preferences in a Maine Forest. Unpubl. master's thesis, University of Maine, Orono. Google Scholar
- Morris, M. A., and R. A. Brandon. 1984. Gynogenesis and hybridization between Ambystoma platineum and Ambystoma texanum in Illinois. Copeia 1984:324–337. Google Scholar
- Neaves, W. B., and P. Baumann. 2011. Unisexual reproduction among vertebrates. Trends in Genetics 27:81–88. Google Scholar
- Nöel, S., P. A. Labonte, and F.-J. Lapointe. 2011. Genomotype frequencies and genetic diversity in urban and protected populations of blue-spotted salamanders (Ambystoma laterale) and related unisexuals. Jounal of Herpetology 45:294–299. Google Scholar
- Nuttle, T. 1997. Densiometer bias? Are we measuring the forest or the trees?Wildlife Society Bulletin 25:610–611. Google Scholar
- Nyman, S., M. J. Ryan, and J. D. Anderson. 1988. The distribution of the Ambystoma jeffersonianum complex in New Jersey. Jounal of Herpetology 22:224–228. Google Scholar
- Osbourn, M. S., G. M. Connette, and R. D. Semlitsch. 2014. Effects of fine-scale forest habitat quality on movement and settling decisions in juvenile pond-breeding salamanders. Ecological Applications 24:1719–1729. Google Scholar
- Petranka, J. W. 1998. Salamanders of the United States and Canada. Smithsonian Institution Press, USA. Google Scholar
- Pittman, S. E., and R. D. Semlitsch. 2013. Habitat type and distance to edge affect movement behavior of juvenile pond-breeding salamanders. Jounral of Zoology 291:154–162. Google Scholar
- Porej, D., M. Micacchion, and T. E. Hetherington. 2004. Core terrestrial habitat for conservation of local populations of salamanders and wood frogs in agricultural landscapes. Biological Conservation 120:399–409. Google Scholar
- Regosin, J. V., B. S. Windmiller, R. N. Homan, and J. M. Reed. 2005. Variation in terrestrial habitat use by four pool-breeding amphibian species. Journal of Wildlife Management 69:1481–1493. Google Scholar
- Rittenhouse, T. A. G., and R. D. Semlitsch. 2007. Postbreeding habitat use of wood frogs in a Missouri oak-hickory forest. Journal of Herpetology 41:645–653. Google Scholar
- Rubbo, M. J., and J. M. Kiesecker. 2005. Amphibian breeding distribution in an urbanized landscape. Conservation Biology 19:504–511. Google Scholar
- Ryan, K. J., and A. J. K. Calhoun. 2014. Postbreeding habitat use of the rare, pure-diploid blue-spotted salamander (Ambystoma laterale). Journal of Herpetology 48:556–566. Google Scholar
- Saino, N. 1992. Selection of foraging habitat and flocking by crow Corvus corone phenotypes in a hybrid zone. Ornis Scandinavica 23:111–120. Google Scholar
- Semlitsch, R. D. 1998. Biological delineation of terrestrial buffer zones for pond-breeding salamanders. Conservation Biology 12:1113–1119. Google Scholar
- Shoop, C. R. 1965. Orientation of Ambystoma maculatum: movements to and from breeding ponds. Science 149:558–559. Google Scholar
- Stebbin, R. C., and N. W. Cohen. 1995. A Natural History of Amphibians. Princeton University Press, USA. Google Scholar
- Thomas, D. L., and E. J. Taylor. 2006. Study designs and tests for comparing resource use and availability II. Journal of Wildlife Management 70:324–336. Google Scholar
- Titus, V., D. Madison, and T. Green. 2014. The importance of maintaining upland forest habitat surrounding salamander breeding ponds: case study of the eastern tiger salamander in New York, USA. Forests 5:3070–3086. Google Scholar
- Uzzell, T. M. 1964. Relations of the diploid and triploid species of the Ambystoma jeffersonianum complex (Amphibia, Caudata). Copeia 1964:257–300. Google Scholar
- Weller, W. F., W. G. Sprules, and T. P. Lamarre. 1978. Distribution of salamanders of the Ambystoma jeffersonianum complex in Ontario. Canadian Field-Naturalist 92:174–181. Google Scholar
- Welsh, H. H., and A. J. Lind. 1995. Habitat correlates of the Del Norte salamander, Plethodon elongatus (Caudata: Plethodontidae), in Northwestern California. Journal of Herpetology 29:198–210. Google Scholar
- Williams, P. K. 1973. Seasonal Movements and Population Dynamics of Four Sympatric Mole Salamanders, Genus Ambystoma. Unpubl. Ph.D. diss., Indiana University, USA. Google Scholar
- Windmiller, B., R. N. Homan, J. V. Regosin, L. A. Willitts, D. L. Wells, and J. M. Reed. 2008. Breeding amphibian population declines following loss of upland forest habitat around vernal pools in Massachusetts, USA. Pp. 41–51in J. C. Mitchell, R. E. J. Brown, and B. Bartholomew (eds.), Urban Herpetology. Society for the Study of Amphibians and Reptiles, USA. Google Scholar
- Wood, E. M., S. E. Barker Swarthout, W. M. Hochachka, J. L. Larkin, R. W. Rohrbaugh, K. V. Rosenberg, and A. D. Rodewald. 2016. Intermediate habitat associations by hybrids may facilitate genetic introgression in a songbird. Journal of Avian Biology 47:508–520. Google Scholar
Summary data for 40 unisexual salamanders radio tracked from 4 vernal pools. ID includes the pool of origin followed by the identification number of each animal. Genomotype indicates both ploidy and how many Ambystoma laterale and Ambystoma jeffersonianum genomes each individual contains. Clone indicates which animals were identical at 3 loci. Body mass and SVL (snout–vent length) were measured under anesthesia prior to transmitter implant surgery. The number of 3-m plots represents the amount of movements > 6 m for which habitat data was recorded, whereas 0.5-m plots represent 1–6 m movements in 2014. Max step indicates the maximum distance moved between successive relocations in 2013 when salamanders were tracked daily, and this distance divided by three in 2014 when salamanders were tracked every 3 days, and estimates the maximum distance traveled in one night. Pool Dist is the maximum Euclidean distance each salamander traveled from the breeding site. The fates of each salamander include mortality events related to the implanted transmitters (MT), mortality events that were unrelated to the transmitters (MU), premature transmitter failure (TF), and battery expiration (BE).