We describe two new cryptic species of Bufo within the subgenus Anaxyrus discovered in Central Nevada of the western United States. Our analyses revealed that these two localized endemic toads are genetically divergent and morphologically distinct, yet were concealed under the range of the broadly distributed western toad (Bufo boreas), which occurs throughout Nevada. The newly discovered species are close in geographic proximity to each other (albeit, in different hydrological basins) but have evolved unique morphological characters that are distinct from each other and distinctive from all allied taxa within the B. boreas species complex. The delimiting of these two rare toads emphasizes the link between taxonomic crypsis and inadequate conservation as these newly described species are vulnerable to extinction due to severely restricted geographic ranges, unknown population sizes, and dependency on rare, fragile wetland habitat, which is a limited resource within Nevada, the primary state that makes up the arid Great Basin. These two endemics join the Great Basin B. boreas species complex as imperiled new members, and our study demonstrates that our knowledge of anuran diversity is incomplete and that new discoveries can still be made, even in unlikely settings.
Amphibians are among the rarest vertebrates in the Great Basin Desert, yet the western toad (Bufo [Anaxyrus] boreas) can be found throughout much of the region with a range that extends across the western United States. The B. boreas species complex, which includes the cosmopolitan B. boreas, plus four narrowly distributed endemics confined to the hydrologic Great Basin: B. canorus, B. exsul, B. nelsoni, and the newly described B. williamsi (Gordon et al., 2017) are examples of the unique aquatic-dependent taxa within this arid ecoregion. Previous analyses examining the genetic diversity and endemism within the B. boreas species complex have suggested that cryptic lineages are likely present within the western toad's broad geographic range, and that the current taxonomy does not accurately reflect the diversity of B. boreas (Stephens, 2001), particularly around the Great Basin (Goebel, 2005; Goebel et al., 2009; Switzer et al., 2009). In our recent molecular examination of Great Basin diversity of B. boreas, we confirmed the presence of cryptic species, such as B. williamsi and two other lineages described here, which were all concealed under the broad range of the western toad. Our extensive collection of morphological measurements of live toads allowed us to quantify significant features that further distinguish these new species from B. boreas and allied taxa within the regional B. boreas species complex. Here, we describe two new species and highlight the consequence of taxonomic crypsis of undescribed species, which are constrained to extremely limited ranges, as these newest novel discoveries have the smallest known geographic distributions within the B. boreas species complex. Like both B. williamsi and B. exsul, these new species are restricted to rare spring-fed wetlands, a habitat within the Great Basin that is vulnerable to habitat loss and exploitation, warranting urgent conservation initiatives to protect and preserve these rare bufonids.
MATERIALS AND METHODS
Data collection and morphological analyses.—We recorded morphological measurements from live adult toads from 19 populations within the hydrological Great Basin (Fig. 1A) including B. boreas (n = 289), B. exsul (n = 30), B. nelsoni (n = 31), B. williamsi (n = 76), plus individuals from both Hot Creek (n = 42) and Railroad Valley (n = 50) to comprise a large data set (n = 518). Fourteen morphological features were recorded and are as follows: snout–vent length (SVL; tip of snout to posterior end of urostyle), head length (HL; tip of snout to occiput), head width (HW; at widest part of the head), snout length (SL; tip of snout to anterior corner of eye), internarial distance (IND; distance between nares), eye diameter (ED; at widest part of eye), interorbital space (IOD; shortest distance between medial margin of upper eyelids), tympanum diameter (TYM; at maximum width of tympanum), paratoid length (PL; horizontal length of parotoid gland) and width (PW; maximum width of parotoid), interparotoid distance (IPD; shortest distance between medial margin of parotoid glands), femur length (FL; distance between vent and knee), tibial length (TL; distance between knee and heel), and hind foot length (FTL; distance from anterior margin of heel to distal end of the third toe). All morphological characteristics were measured using Mitutoyo digital calipers to a precision of 0.1 mm. All individuals were measured by ETS with the exception of 13 individuals collected from Hot Creek and 19 individuals collected from Railroad Valley, which were measured by MRG, including the holotypes and paratypic series from each site. Digit length from hands and forelimbs were recorded by MRG. Sex was determined in the field, noting body size, behavior, and secondary sex characteristics, such as the utilization of a release call and presence of nuptial pads on males as identifiers.
Individuals collected in May of 2015 from populations in Hot Creek Canyon and the eastern neighboring basin, Railroad Valley, were selected to represent the holotypes and paratypic series for both species and were euthanized and preserved following the guidelines under the Institutional Animal Care and Use Committee (IACUC) from University of Nevada (IACUC #00068). The extracted tissue samples were preserved and deposited at the University of Nevada, Reno (UNR). Specimens were fixed in 10% buffered formalin and transferred to 70% ethanol and deposited as vouchers in the California Academy of Science (CAS) and University of Nevada, Reno.
To characterize morphological differences among species, we used multivariate analysis of covariance (MANCOVA) and used SVL as the covariate in these analyses to account for body-size variability among species (Dahl and Peckarsky, 2002; McCoy et al., 2006). This analysis results in least squares means generated from regressions for each size corrected variable against SVL, which can identify subtle, but statistically significant differences in fine features examined in these toads. We also log transformed the raw dataset as an additional way to account for differences in allometry of measured toads and analyzed the scaled data set using MANCOVA to quantify morphological differences among the species (Lleonart et al., 2000). We used Tukey HSD post hoc pairwise comparisons to identify significant character state differences among the species examined (Feinberg et al., 2014). A cross-validated discriminant function analysis (DFA) was used to evaluate the variation in multivariate space to identify variables that best discriminated among the species (Feinberg et al., 2014). Data collected from multiple measurers can result in inter-observer error due to variations in character assessment, particularly on the fine features of amphibian anatomy (Hayek et al., 2001). These errors can produce results that may bias biological interpretations from morphological analyses (Hayek et al., 2001). To avoid inter-observer biases, only measurements collected by ETS were used in the morphological analyses with the exception of the means table (Table 1), which includes combined raw, unadjusted measurements from ETS and MRG. All statistics were conducted using JMP Pro v. 10 (SAS, Cary, NC).
Morphological variation of six species that comprise the B. boreas species complex of the Great Basin. Fourteen morphological measurements (in mm) are as follows: snout–vent length (SVL), head length (HL), head width (HW), snout length (SL), internarial distance (IND), eye diameter (ED), interorbital distance (IOD), tympanum diameter (TYM), parotoid width (PW), parotoid length (PL), interparotoid distance (IPD), femur length (FL), tibial length (TL), and foot length (FTL). Data include sample size (n), character mean ± standard deviation, and range. All measurements reflect raw, unadjusted values.
Genetics.—Following the methods described in Gordon et al. (2017), tissue samples were collected from individuals identified as B. boreas (Fig. 1B; n = 308), B. canorus (n = 32), B. exsul (n = 30), B. nelsoni (n = 32), and B. williamsi (n = 7). The control region (CR) of the mitochondrial genome was selected due to the region's high rate of evolution ideal for intraspecific analyses (Avise et al., 1987) and because it has been used in previous phylogenetic studies evaluating the diversity in B. boreas (Stephens, 2001; Goebel et al., 2009). Primers were modified using previously published sequences (Goebel et al., 1999; Stephens, 2001) to amplify approximately 1,622 bp area of the CR (Table 2). PCR products were sequenced using ABI 3730 sequencer, and data were analyzed with Sequencher software (Gene Codes, Ann Arbor, MI). The final alignment of the B. boreas species complex (CR 1,622 bp) was completed using ClustalW (Larkin et al., 2007) within Mega 7.0 (Kumar et al., 2016), resulting in 62 unique haplotypes included in further analyses. Haplotype, locality, species, and GenBank accession numbers are in Table 3. To examine pairwise genetic distances among sequences relative to haplotypes identified, a Jukes-Cantor model (Jukes and Cantor, 1969) was applied in Mega 7.0 (Kumar et al., 2016).
Primers used to amplify the mitochondrial control region (CR) in B. boreas and B. punctatus.
Mitochondrial control region sequence data for B. boreas species complex. Haplotype, locality, species, and GenBank accession number. Sequences from Colorado and Utah were downloaded from GenBank (Goebel et al., 2009) with specimen locality and corresponding accession number for B. boreas DNA data.
Genetic analyses.—Previous molecular studies support evidence of recent divergence of allied taxa within the B. boreas species complex (Feder, 1973; Graybeal, 1993; Shaffer et al., 2000; Stephens, 2001; Pauly et al., 2004; Goebel et al., 2009; Switzer et al., 2009). Due to the close ancestry of this species group, we constructed a TCS haplotype network in PopART to examine population level genealogy (Clement et al., 2002; Leigh and Bryant, 2015). Phylogenetic hypotheses were tested using Bayesian inference (BI) in MrBayes v3.1.2 (Ronquist and Huelsenbeck, 2003) using parameters described in Gordon et al. (2017). The program Tracer v1.6 (Rambaut et al., 2014) confirmed that analyses reached stationarity and trees were constructed using FigTree v1.4.2 (Rambaut, 2014). Maximum likelihood (ML) phylogenies were constructed in Mega 7.0 (Kumar et al., 2016) to examine evolutionary relationships and comparative tree topologies amonh taxa of the B. boreas species complex. The evolutionary history was inferred based on the Hasegawa-Kishino-Yano model (Hasegawa et al., 1985) suggested by Akaike Information Criterion (AIC; Burnham and Anderson, 2002), and support was evaluated using 500 bootstrap replications (Pattengale et al., 2010). A discrete gamma distribution was used to model evolutionary rate differences among sites (5 categories [+G, parameter = 0.6336]). Bufo punctatus was selected as the outgroup since this taxon was included in previous studies that examined the same molecular marker investigating the fine-scale relationship of toads within the B. boreas species complex (Stephens, 2001; Goebel et al., 2009). All positions containing gaps and missing data were eliminated. There were a total of 1,401 positions in the final dataset.
To examine the relationships of populations of B. boreas outside the Great Basin to the variant haplotypes for a broader geographic and historical context, sequences for the control region from Goebel et al. (2009) were downloaded from GenBank and added to the dataset (Table 3) in Mega 7.0 (Kumar et al., 2016). To reconstruct the evolutionary history, 87 unique haplotypes were used in ML based on the Hasegawa-Kishino-Yano model (Hasegawa et al., 1985). All positions containing gaps and missing data were eliminated. There were a total of 628 positions in the final dataset using ML in Mega 7.0 (Kumar et al., 2016).
Morphological evidence.—There were significant differences for all 14 size-corrected morphological characters among B. boreas, B. exsul, B. nelsoni, B. williamsi, and populations of Hot Creek (Bufo monfontanus, new species) and Railroad Valley (Bufo nevadensis, new species). Bufo nevadensis, new species, is significantly smaller overall than larger bufonids B. boreas and B. nelsoni (F5,439 = 62.4, P < 0.0001), yet it has a significant, but modestly larger head (F6,439 = 704.0, P < 0.0001; F6,439 = 899.6, P < 0.0001; Table 4), which is longer than similar body-sized congeners, B. exsul, B. monfontanus, new species, and B. williamsi, with moderate, well-separated eyes (F6,439 = 290.1, P < 0.0001). Significant differences were detected among species for snout length (F6,439 = 145.8, P < 0.0001), and in pairwise comparisons recovered from Tukey HSD post hoc tests (Table 4), B. nevadensis, new species, has a short snout, which, while perceptibly subtle, differs significantly from all other congeners examined in this study except B. exsul. The parotoid glands of B. nevadensis, new species, are significantly reduced; parotoid width is narrower than all congeners examined (F6,439 = 145.1, P < 0.0001) and shorter in length than B. boreas and B. monfontanus, new species (F6,439 = 156.6, P < 0.0001; Table 4). Additionally, the small-sized B. nevadensis, new species, has relatively long legs; the femur is comparatively longer than all other congeners with the exception of the larger toads B. boreas and B. nelsoni (F6,439 = 741.2, P < 0.0001; Table 4), and tibial and foot length are longer than the similar-sized toads B. williamsi and B. monfontanus, new species (F6,439 = 770.8, P < 0.0001; F6,439 = 392.3, P < 0.0001; Table 4).
Table of least squares means generated from regressions of MANCOVA analyses for SVL and 13 size corrected traits (see Materials and Methods) against SVL for each species examined within the B. boreas species complex. Each value includes upper and lower 95% confidence intervals (CI). Bolded values indicate statistically significant different states in pairwise comparison identified in Tukey HSD post hoc tests (P < 0.05) with larger (↑) and smaller character states (↓) identified.
Bufo monfontanus, new species, has a significantly smaller body size than the larger sized toads, B. boreas and B. nelsoni (F5,439 = 62.4, P < 0.0001; Table 4), and outsizes only B. williamsi (Tables 1, 4). This small toad has a wide head similar to B. boreas and B. nevadensis, new species, but it has the shortest relative head length of all species within the complex (F6,439 = 899.6, P < 0.0001; Table 4), with moderate, well-separated eyes (F6,439 = 290.1, P < 0.0001; Table 4). Significant differences in parotoid width were detected among the species complex (F6,439 = 145.1, P < 0.0001) and in pairwise comparisons in Tukey HSD post hoc tests, B. monfontanus, new species, exhibits the longest parotoid gland, which differs significantly, though the physical difference is modest, from all congeners of the complex except B. boreas (F6,439 = 156.6, P < 0.0001; Table 4). The large glands are comparatively closer together, which is significantly different from that presented in B. boreas, B. nelson, and B. nevadensis, new species (F6,439 = 380.2, P < 0.0001; Table 4). The leg features of B. monfontanus, new species, are reduced with significant differences detected; while perceptibly subtle, B. monfontanus, new species, has the shortest femur among all congeners of the complex (F6,439 = 741.2, P < 0.0001; Table 4), shorter tibia differing from lengths observed in B. boreas, B. nelsoni, and the small toad, B. nevadensis, new species (F6,439 = 770.8, P < 0.0001; Table 4), and a small foot length distinct from B. boreas and B. nevadensis, new species (F6,439 = 392.3, P < 0.0001; Table 4).
In multivariate morphospace using DFA, significant differences were detected among species (F70,323 = 3.40, P < 0.0001; Fig. 2). The DFA correctly classified 83.7% of predicted species with some morphological overlap detected among B. boreas, B. exsul, B. nelsoni, and B. nevadensis, new species. Little morphological overlap was detected with B. monfontanus, new species. Canonical 1 explained 61% of the variation with head width loading most heavily, while canonical 2 accounted for 17% of the variation with parotoid length loading more heavily than other characters (Fig. 2).
To assess sexual dimorphism in B. monfontanus, new species (females = 16, males = 13), and B. nevadensis, new species (females = 28, males = 3), intraspecific analyses showed significant differences in SVL detected only in B. monfontanus, new species (F1,28 = 5.4, P < 0.03), with females significantly larger than males (P < 0.05). However, males of B. monfontanus, new species, have relatively longer legs than females (F2,28 = 12.7, P < 0.0001 for femur length; F2,28 = 24.1, P < 0.0001 for tibial length; F2,28 = 19.0, P < 0.0001 for foot length). While sexual dimorphism was not detected in SVL of B. nevadensis, new species (F1,30 = 1.2, P = 0.28), males did have relatively longer feet (F2,30 = 8.4, P < 0.001). The unadjusted data collected for the six species examined for all 14 characters are presented in Table 1.
Genetic evidence and phylogenetic relationships.—Similar to the findings presented by Gordon et al. (2017), the combined analyses for the control region of the mitochondrial genome examined in the B. boreas species complex study supported the existence of divergent lineages of undefined toad populations within the Great Basin, warranting inspection of cryptic populations occurring in Central Nevada and described herein. The TCS haplotype network and phylogenetic reconstructions support four major clades that correspond roughly to geographical regions within the hydrological Great Basin and defined as Western Great Basin (W), Humboldt-Lahontan (HL), Mojave (S), and Eastern Great Basin (E; Figs. 1B, 3, 4, Supp. Figs. 1, 2; see Data Accessibility). The haplotype network highlights the divergent lineages of B. monfontanus, new species, and B. nevadensis, new species, from each other and to all B. boreas within the hydrological Great Basin (Fig. 3A) and illustrates the regional divide in diversity among populations of B. boreas and related taxa within the species complex (Figs. 1B, 3). All phylogenies resulted in minor topological differences, but all results supported the divergence of the eastern Great Basin, which includes the newly described species B. monfontanus, new species, and B. nevadensis, new species, forming terminal taxa, and which are distinct from all congeners of the B. boreas species complex (Figs. 3B, 4, Supp. Figs. 1, 2; see Data Accessibility). In all combined analyses, the Mojave clade represents the most species-rich region of the Great Basin, with B. exsul and B. nelsoni plus haplotypes of both B. canorus and B. boreas. Bufo boreas and B. canorus appear to be polyphyletic or paraphyletic with lineages that occur in the Humboldt-Lahontan and Mojave clades, a result consistent with previous studies (Graybeal, 1993; Goebel, 1996; Shaffer et al., 2000; Stephens, 2001; Goebel et al., 2009; Switzer et al., 2009; Gordon et al., 2017).
Assessment of nucleotide diversity evaluating genetic distances for both B. nevadensis, new species, and B. monfontanus, new species, uncovered higher percentages of differentiation than other congeneric taxa within the B. boreas species complex (Table 5). Bufo nevadensis, new species, average genetic distance compared to B. boreas is 3.5% and B. monfontanus new species is 3.2%, and despite their seemingly close relative proximity (Figs. 1, 5B), these two species are highly differentiated from each other at 2.0% divergence (Table 5).
Estimates of evolutionary divergence between mitochondrial DNA control region sequences from taxa of the B. boreas species complex within the Great Basin. The number of base substitutions per site from between sequences is shown. Analyses were conducted using Jukes-Cantor model (Jukes and Cantor, 1969), which involved 14 nucleotide sequences (Fig. 1B), in Mega7 (Kumar et al., 2016). All positions containing gaps and missing data were eliminated, which resulted in a final dataset of 1429 positions. Pairwise comparisons against congeners of the B. boreas species complex for B. nevadensis and B. monfontanus are indicated in bold.
Comparing Bufo nevadensis, new species, and B. monfontanus, new species, to populations of B. boreas outside of the Great Basin, we used maximum likelihood to test evolutionary hypotheses which yielded support for the new species' divergence and illustrates close ancestry with populations of B. boreas in Utah and Colorado (Fig. 4). Bufo nevadensis, new species, forms a terminal clade and is basal to populations of boreal toads near the northwestern Utah border and Colorado (Table 3, Fig. 4), and genetic distance of this species to B. boreas was 1.7%. Bufo monfontanus, new species, forms a terminal end and represents a divergent lineage, but shares a haplotype for this marker with one boreal toad from the northwestern corner of Utah in Box Elder County (Fig. 4) with an average genetic distance of 1.3%. Additional molecular markers may provide greater insight into the fine-scale relationships of these newly described species and taxa of the B. boreas species complex.
Bufo (Anaxyrus) nevadensis, new species
Railroad Valley Toad
Figures 1–4, 5B, 6; Tables 1, 4–5
Paratypes.—UNR 7905, adult male; UNR 7906, adult female; UNR 7907, adult female; UNR 7908, adult male; UNR 7909, adult female; all individuals collected within the identified home range in Railroad Valley, 38°33′9.1″N, 115°46′12.8″W, M. R. Gordon, K. Guadalupe, and C. Burg, 5 May 2015 (Fig. 5B).
Diagnosis.—Bufo (Anaxyrus) nevadensis is a member of the Great Basin B. boreas species complex (Blair, 1972), but traditionally has been identified as B. boreas due to its occurrence within the Western Toad's geographic range, yet is distinct from B. boreas by a combination of morphological characters (Figs. 2, 6, Tables 1, 4), genetic evidence (Figs. 3, 4, Table 5), and restricted geographic distribution (Fig. 5B). Bufo nevadensis is distinguished from B. boreas due to its small adult body size (SVL is approximately 2 cm smaller than B. boreas; Tables 1, 4); significantly, but modestly longer head with a relatively shorter snout; well-separated, perceptibly short and narrow parotoid glands; significantly, but comparatively long legs, large hind feet (Fig. 6B); and distinctive mottling of venter (Fig. 6B, D).
Bufo nevadensis is among the smallest terrestrial bufonids within the B. boreas species complex (Tables 1, 4). However, this new species exhibits a relatively large head unlike similarly small toads, B. exsul and B. monfontanus, new species, with a significantly, but comparatively shorter snout distinctive from all species within the complex except B. exsul (Table 4). The well-separated and severely reduced parotoid glands exhibited in B. nevadensis is divergent from all taxa within the B. boreas species complex, and the shortened gland length distinguishes B. nevadensis from both B. boreas and B. monfontanus, new species. Bufo nevadensis has statistically significant, relatively long legs; longer femur than exhibited in B. exsul, B. monfontanus, new species, and B. williamsi, and a longer tibia and hind feet, which separate B. nevadensis from B. monfontanus, new species, and B. williamsi (Table 4). In addition to morphological shape differences, B. nevadensis displays a dominantly brown and gray toned dorsum with prominent warts and heavily creased skin, which differs from B. exsul, B. monfontanus, new species, B. nelsoni, and B. williamsi. The venter of B. nevadensis is similar to B. exsul and B. williamsi, exhibiting black mottling contrasted against a white background color on the anterior sides of the limbs and belly. The presence of a dorsal stripe is extremely variable among individuals of B. nevadensis, as is similar to the other members of the B. boreas species complex, with the exception of B. exsul. Small, irregular tibial glands may be present in individuals, but this characteristic is highly variable.
In mature male B. nevadensis, distinct nuptial pads develop on the dorsal side of the first finger, a typical secondary sexual characteristic exhibited among most bufonids. This species lacks an advertisement call, but it emits a release call when males come in contact with one another which sound like the weeping of a chick and is similar among congeners of the B. boreas complex (Stebbins, 2003).
Description of holotype.—Body relatively small and robust (SVL = 62.5 mm); head wider (HW = 20.1 mm) than long (HL = 17.2 mm; 85% head length to head width). Snout subelliptical in dorsal view; snout profile moderately truncate in lateral view (SL = 7.5 mm; 1.5 times longer than eye diameter). Canthus rostralis distinct and cuneate. Loreal region slightly concave. Nostrils slightly protuberant, directed dorsolaterally, and closer to anterior corner of eye than to snout. Internarial distance (IND = 3.8 mm) 75% of interorbital distance (IOD = 4.8 mm). Eyes well separated; interorbital space nearly equivalent to eye diameter (ED = 4.8 mm). Upper eyelids prominent in dorsal view; eyes slightly breach profile margin. Tympanum (TYM = 3.0 mm) distinct, subovoid, relatively small (52% of eye diameter). Supratympanic fold present. Parotoid glands weakly present viewed above; parotoid glands narrow (PW = 3.7 mm), severely tapered at posterior corner of eye in lateral view. Parotoid glands 1.5 times longer (PL = 7.5 mm) than eye diameter; parallel, well separated (IPD = 12.3 mm). Forearms robust. Fingers unwebbed; relative length III > I > II > IV; tips rounded, subarticular tubercles moderate, accessory palmar tubercles small, round. Inner metacarpal tubercle distinct, round. Palmar tubercle prominent, elliptical. Hind limbs long; femur slightly longer (FL = 21.2 mm) than tibia (TL = 18.7 mm). Tibial glands irregular, scarcely defined, half the length of parotoid gland. Tarsal fold present; hind feet webbed proximally (FTL = 37.0 mm). Relative toe lengths IV > III > V > II > I; toe tips rounded. Subarticular tubercles faintly evident, small, round. Plantar tubercles numerous, small. Inner metatarsal tubercle pronounced, elevated, relatively large, elliptical. Outer metatarsal tubercle distinct, ovoid. Skin warty on dorsum; primary warts elevated, irregular; finely granular skin between elevated warts from interorbital space increasing in coarseness toward posterior margin of dorsum at articulation with femur. Hind limbs warty; tubercles small, moderate. Venter coarse, seat patch conspicuous.
Color in life.—Dorsal ground color of holotype light brownish gray, flecked with dark brown, irregular spotting (Fig. 6A). Grayish head dappled with dark brown patches, upper eyelids finely speckled black. Smooth dark olive brown blotches along prefrontal to frontal area of head and interorbital space. Pupil black, horizontal with gold streaked iris. Brownish gray parotoid glands exhibit minor black to dark brown spotting. Cream colored dorsal stripe present, originating posteriorly at nares and terminating at posterior margin of urostyle. Elevated, dark olive brown warts at interparotoid space, along dorsum; some warts set in black blotches; olive streaking connects warts along midline bordering dorsal strip; dorsolaterally, olive to dark brown wart color streaked against brownish gray background color. At midaxillary line, black streaking contrasts against white background. Throat white with minor black spotting near lower lip. Venter mottled black against white background color (Fig. 6B). In dorsal view, forelimbs have minor dark olive banding and olive patches against light grayish brown ground color; hind limbs exhibit dark olive brown banding and patches against light grayish brown ground color down to heavily dappled olive feet atop brownish gray background color. Along medial and ventral sides of hind limbs, black spotting and patches occur against white ground color down to medial edge of feet, which appear dark gray on underside.
Variation.—In paratypes, background color ranges from light gray to light brownish gray; warts may be olive colored to olive brown and may be set in black blotches. Brown to dark brown spotting around front of head; minor black spotting on throat and minor to heavy black mottling occur on venter.
Color in preservative.—Color is similar to that in life (Fig. 6A, B) with some notable differences. Overall background color in the holotype appears gray (Fig. 6C). Along dorsum, dark brown warts flattened and appear as heavy streaks against brownish gray background adjacent to dorsal stripe. Marbling along midaxillary line and venter less vibrant than in life and appear dark gray against white ground color and seat patch muted light gray (Fig. 6D).
Distribution.—Bufo nevadensis is known only to occur near and within the spring-fed wetland areas of Lockes Ranch (1,460 m above sea level), a protected wildlife management area located in Railroad Valley, an east-central desert basin between the Pancake Range and Grand Range of Nye County, Nevada (Fig. 5B). The critical marshland habitat for this endemic toad is solely fed from Big, Reynolds, and Hay Corral springs which results in a severely restricted range with an estimated distribution of 1.8 km2. These outflows are remote and isolated, surrounded by cold desert habitat dominated by Big Sagebrush (Artemisia tridentata ssp. tridentata), Greasewood (Sarcobatus vermiculatus), Rubber Rabbitbrush (Ericameria nauseosa), and saltbush (Atriplex spp.) with limited usable corridors for amphibian dispersal, which likely restricts this species' movement to other spring localities within Railroad Valley. Bufo nevadensis co-occurs with the federally listed threatened Railroad Valley Springfish, Crenichthys nevadae, and the Great Basin Spadefoot, Spea intermontana.
Life history and behavior.—Bufo nevadensis is nocturnal, emerging at dusk, and can be found in shallow water or among the vegetation in the perimeter band that transitions from riparian to sagebrush steppe habitat. Characteristic of cold deserts, Railroad Valley experiences extreme fluctuations in day and nighttime temperatures as well as season-to-season variation. As is common for other members of the B. boreas complex, these toads likely retreat to burrows in the fall, not emerging until spring, when males begin to congregate in shallow water for breeding. Mature males, similar to other members of the B. boreas complex (with the singular exception of B. canorus), do not have an advertisement call, but emit a release call when males come in close contact with one another. Egg masses and tadpoles develop in still, shallow water amid the marshy vegetation of the wetland habitat.
The population size for this species is unknown; however, the extreme isolation and restricted range may indicate that the population numbers may be small. Little is known regarding the dispersal and non-breeding behavior of this rare toad.
Etymology.—The species name is a derivative from the state of Nevada (U.S.A) where this rare toad occurs and pays homage to the unique biodiversity found in the desert landscape of its home state.
Remarks.—Railroad Valley is a geothermally active area within the Range and Basin Province with significant opportunities for anthropogenic energy production, including extraction of its oil reservoirs (Liu et al., 1997) that continue to contribute to ongoing economic interests in the valley which are currently overseen by the Bureau of Land Management. Discovery of this rare new species should elicit high conservation concerns due to its severely restricted range and limitations to dispersal due to isolation and remoteness of the spring-fed habitat upon which B. nevadensis is dependent. Any further human anthropogenic modifications of habitat that may degrade this extremely important habitat would imperil this toad. However, with new species designation, conservation initiatives would provide a platform for ongoing policy and monitoring to allow this this endemic toad to persist.
Bufo (Anaxyrus) monfontanus, new species
Hot Creek Toad
Figures 1–4, 5B, 7; Tables 1, 4–5
Paratypes.—UNR 7910, adult male; UNR 7911, adult male; UNR 7912, adult male; UNR 7913, adult male; UNR 7914, adult male; UNR 7915, adult male; UNR 7916, adult male; UNR 7917, adult male; all individuals collected within the identified home range, Hot Creek Canyon, 38°32′19.32″N, 116°27′32.9″W (Fig. 5B), M. R. Gordon, 6 May 2015.
Diagnosis.—Bufo (Anaxyrus) monfontanus occurs within the range of B. boreas but is distinct from the Western Toad by a combination of diagnostic morphological characters (Figs. 2, 7, Tables 1, 4), genetic evidence (Fig. 3, Table 5), and restricted geographic distribution (Fig. 5B). Bufo monfontanus is distinguishable from B. boreas by having a small adult body size (SVL is 2 cm smaller than B. boreas; Tables 1, 5); significantly, but modestly shorter head; perceptibly large, parotoid glands; significantly, but comparatively shorter legs with small hind feet; and weakly warted body (Fig. 7A).
Bufo monfontanus is among the smallest bufonids within the B. boreas species complex, and it is only larger than B. williamsi (Tables 1, 4). Bufo monfontanus has a significant, but relatively shorter head with a comparatively long snout, with a relative head width more comparable to larger-sized taxa B. boreas and B. nelsoni (Table 4). An important diagnostic feature among B. monfontanus is the presentation of well-defined, relatively large parotoid glands, which distinguishes this small toad from all other small-sized toads within the B. boreas species complex (Table 4). Bufo monfontanus has shorter legs; significantly, but relatively the shortest femur and tibia of all taxa within the B. boreas species complex and relatively small feet distinct from both B. boreas and B. nevadensis (Table 4). The dorsal stripe is extremely variable among individuals of B. monfontanus; a characteristic typical among taxa within the B. boreas species complex with the exception of B. exsul. Small, irregular tibial glands may be present among individuals of B. monfontanus as seen in both B. nevadensis and B. williamsi.
In adult male B. monfontanus, distinct nuptial pads develop on the dorsal side of the first finger which is a typical secondary sexual characteristic exhibited among most bufonids. Akin to congeners of the B. boreas species complex, except B. canorus, males of this species emit a release call that sounds like a weeping chick (Stebbins, 2003).
Description of holotype.—Body small, robust (SVL = 59.6 mm); head wider (HW = 19.3 mm) than long (HL = 16.2 mm). Snout subovoid in dorsal view; snout rounded in lateral view, long (SL = 6.4 mm; 40% of head length). Canthus rostralis distinct, concave, angular, sloping up to medial orbital margins. Loreal region slightly concave. Nostrils protuberant, directed dorsolaterally, closer to anterior corner of eye than end of snout. Internarial distance 75% of eye-to-naris distance (IND = 2.3 mm). Relatively moderate eyes well separated (ED = 4.4 mm; IOD = 3.8 mm); eyes do not breach snout profile in dorsal view. Tympanum distinct, oval, small (TYM = 2.5 mm; 53% of eye diameter). Supratympanic fold weakly present, flat. Parotoid glands longer (PL = 8.8 mm) than wide (PW = 5.3 mm; length nearly twice eye diameter). Glands elevated, slightly convergent at posterior ends in dorsal view. In lateral view, parotoid glands elongated longitudinally from posterior corner of eye, oval, wider than the eye (1.4 times eye diameter). Interparotoid space large (IPD = 9.4 mm; 2.5 times interorbital distance). Forearms robust, smooth. Fingers unwebbed; relative lengths III > I > IV > II; nuptial pads present, raised on dorsal side of digit I; tips rounded, subarticular tubercles moderate, round; accessory palmar tubercles small, round. Inner metacarpal tubercle raised, prominent, round. Palmar tubercle distinct, large, subovoid, borders medial margin of inner metacarpal tubercle. Hind limbs long; femur longer (FL = 23.1 mm) than tibia (TL = 19.7 mm). Tibial gland weakly present in dorsal view and equivalent to the width of the parotoid gland. Tarsal fold present. Hind feet webbed proximally (FTL = 35.0 mm). Relative toe lengths IV > III > V > II > I; tips rounded. Subarticular tubercles distinct, small and round; plantar tubercles numerous, small. Inner metatarsal tubercle conspicuous, elevated, relatively large and elliptical. Outer metatarsal tubercle pronounced, ovoid. Longitudinally along dorsum, dorsal stripe weakly present, originating posterior to interorbital space, terminating at urostyle; irregular, elevated but scattered warts present, increasing in size from interorbital space to posterior margin of urostyle. Skin between warts nearly smooth; forearms smooth; hind limbs exhibit minor warts, tibial gland scarcely present and irregular. Originating posterior to labial commissure, inferior to tympanum, dense, small tubercles occur along posterior axillary line, terminating near anterior articulation of femur. Venter coarse; seat patch granular and conspicuous.
Color in life.—Dorsal background color light olive gray with minor black flecks throughout dorsum; elevated brown warts encircled by narrow, incomplete black border (Fig. 7A). More pronounced black rings border brown warts laterally along posterior axillary line. Dorsal stripe present; fine line, light green, interrupted at parotoid region. Pupils horizontal; iris flecked gold. Parotoid glands flecked with black, minute spots against olive gray. Thick brown stripe present, inferior to eye, occurring from orbit anterior border down to upper lip margin; buffy colored tubercles originate at labial commissure, transition to more olive colored spines moving laterally along body. Along midaxillary line, olive background color transitions to cream color with heavy black blotching. Throat clear, buffy colored. Venter nearly clear, buffy with minor black spots except at seat patch, which is gray and coarse with fine white granules (Fig. 7B). Dorsally, forelimbs exhibit little color variation from olive with some minor dark brown patches; hind limbs have minor black stripes atop white background color along inner shank down to dorsal side of foot.
Variation.—In paratypes, dorsal stripe present or weakly present; dorsum ground color variable shades of olive; venter may exhibit minor black spotting against buffy ventral ground color. Parotoids directed parallel behind eyes or slightly convergent at posterior ends. Brown striping above lip and inferior to eyes may be weakly present.
Color in preservative.—There are notable differences in preservation (Fig. 7C, D) when compared in life (Fig. 7A, B). Overall ground color of preserved holotype is dark and gray. Dorsal stripe thin, white, broken at posterior edges of parotoid glands, continues to vent. Warts along dorsum flattened, dark gray, and similar in shade to ground color and parotoid glands. Head dark gray, striping inferior to eye muted and tubercles posterior to labial commissure appear white. Arms very smooth and dark. Nuptial pads present and brown in preservative. Marbling along midaxillary line black against muted gray ground color. Venter light gray overall; seat patch muted in preservative. Tubercles on feet and hands diminished and dark, tips appear brown.
Distribution.—Bufo monfontanus is only found within the marshes fed by thermal spring outflows in Hot Creek Canyon (1,859 m above sea level) in Hot Creek Range in Central Nevada. The narrow canyon is nestled between Box and Corral Canyons and runs east to west, and toads have only been found in the small Hot Creek stream fed by Upper Warm Springs. This locality is extremely remote and isolated, and the spring flows through the canyon are interrupted, likely restricting toad dispersal from east to west, with few suitable corridors outside the canyon. The estimated range for B. monfontanus is extremely small at 1 km2, as essential wetland habit is critically limited and isolated by cold desert habitat dominated by sagebrush steppe. The population size is unknown for this endemic toad, but is likely small and warrants further examination. Bufo monfontanus co-occurs with introduced Crenichthys nevadae (Railroad Valley Springfish) near Upper Warm Springs, and the introduced Rana catesbeiana (American Bullfrog) can be found within the interior streams of the canyon.
Life history and behavior.—Bufo monfontanus emerges only after dusk, exhibiting typical nocturnal behavior similar among taxa of the B. boreas species complex with the exception of B. canorus, which is diurnal. This species can be found in the marshy water or in the perimeter band that transition from riparian habitat to sagebrush steppe. Typical of other species in the B. boreas species complex, B. monfontanus likely hibernates, using burrows. Little is known about the dispersal behavior of this toad. The high elevation and extreme temperatures, both daily and season-to-season variability, likely prohibit a long breeding season, and further investigation through monitoring and annual surveys could provide insight into the life history strategy of this unique toad.
Etymology.—The species name monfontanus (from the Latin “mons” for mountain and Latin “fons,” a spring or fountain) is descriptive of the high-elevation spring habitat where this toad occurs and pays tribute to the nature of the rare spring habitat and the biodiversity relying on this important resource within the Nevada.
Remarks.—This new toad species, similar to other taxa of the B. boreas species complex, represents another narrow endemic that is found only in small and fragile spring ecosystems of the Great Basin. This new species warrants high conservation concern and urgent initiatives to monitor and to study how to preserve this new toad species. Especially troubling is that little is known about the population size, breeding, and dispersal of this small toad species. Additionally, the predatory generalist Rana catesbeiana co-occurs within the spring-fed stream system of Hot Creek Canyon, and this species is known to outcompete and prey upon much smaller anurans such as bufonids. Rana catesbeiana is also a known vector for potentially lethal anuran diseases, such as chytridiomycosis (Kats and Ferrer, 2003; Daszak et al., 2004), whose effects on the narrowly distributed B. monfontanus is unknown.
Our combined genetic and comparative morphological evidence strongly support recognizing two newly described toads, B. nevadensis and B. monfontanus, increasing the diversity of the Great Basin B. boreas complex to seven species. Taxonomic nomenclature within the Nearctic genus Bufo remains unstable, so we recommend retaining the two new species in the genus Bufo, with Anaxyrus used at the subgeneric rank (Pauly et al., 2009). The descriptions of these two new species bring the number of species of Anaxyrus to 25 (Pauly et al., 2009; Frost, 2015). These new species are morphologically distinct (Table 1) and genetically differentiated (Fig. 3) from the broadly distributed B. boreas, as well as from each other, adding to the intricate phylogeographic story of the B. boreas complex of the arid Great Basin. Bufo nevadensis and B. monfontanus are seemingly close in geographic proximity, separated by approximately 61 km of mountainous desert landscape, yet both have evolved unique phenotypic traits that are nearly opposite of each other in some respects. Bufo nevadensis, found only a small locality within Railroad Valley, Nye County, Nevada, is overall squat and warty, but with dramatically diminished parotoid glands (Fig. 6A, C), while the relatively high elevation toad, B. monfontanus of Hot Creek Canyon, Nye County, Nevada, which is nestled in the mountains of Hot Creek Range, exhibits relatively large, close-set parotoid glands and a weakly warted, comparatively slender body (Figs. 6C, 7A). Additionally, sexual dimorphism was only detected for B. monfontanus. While MANCOVA analyses are robust to uneven sample sizes, further analysis of a larger sample of B. nevadensis males may uncover sexual variation among B. nevadensis. The life histories for both of these newly named species warrants serious attention as the population numbers are unknown, and the geographic ranges for these species are among the smallest known among the taxa of the B. boreas species complex. Breeding and reproduction timing have been documented in B. nevadensis, but not in B. monfontanus. However, this may be due to phenological differences in the species due to extreme dissimilarities in temperature and weather experienced by the mountain-inhabiting B. monfontanus. These parameters may limit the reproductive window, akin to the environmental restrictions experienced by B. canorus (Karlstrom, 1962), a narrowly distributed high elevation relative found only in the Sierra Nevada Mountains (Fig. 5B). Despite a strong difference in habitat elevation, both B. nevadensis and B. monfontanus rely on rare spring-fed wetlands, a habitat that is often a small and isolated aquatic resource within the state of Nevada, a region that comprises much of the interior of the arid Great Basin Desert (Sada and Vinyard, 2002). Although scarce, these riparian habitats are important hubs supporting widespread biodiversity and are often identified as sites rich in endemism, demonstrated by the high taxonomic diversity of spring fish (Hubbs and Miller, 1948; Hewitt, 1996, 2000; Smith et al., 2002; Finger and May, 2015), spring snails (Hershler and Sada, 2002; Sada and Vinyard, 2002), riparian insects (Shepard, 1992) and toads (Wang, 2009; Gordon et al., 2017). While new species within the region have been recognized, undetected and cryptic diversity is still likely within the region as many of these aquatic resources are difficult to detect, severely isolated (Shepard, 1993), and rarely studied.
The Great Basin B. boreas species complex presents an ideal vehicle for evolutionary study demonstrated by the recent discoveries of cryptic species B. williamsi (Gordon et al., 2017) and the two new species described here. While neither species is sympatric with B. boreas, these divergent lineages represent new evolutionary trajectories from a common ancestor shared with the western toad, which is supported by molecular evidence (Fig. 3, Table 5) and emphasized by the variations exhibited in morphology by the region's toad species (Fig. 5B). The Great Basin Desert may appear to be an unlikely setting for new bufonid species, as amphibians are among the rarest animals to occur within this region. However, the occurrence of this species complex, whose localized endemics are confined to extremely restricted ranges within this exceptionally dry ecoregion, continues to spur study examining range-wide diversity of B. boreas to gain insight into the evolutionary relationships among these close relative toads (Graybeal, 1993; Stephens, 2001; Goebel, 2005; Goebel et al., 2009; Switzer et al., 2009). An interesting result from our molecular study suggests that B. nevadensis and B. monfontanus are more closely related to western toads in Colorado (B. nevadensis) and Utah (B. monfontanus). In the work of Goebel et al. (2009), genetic analyses suggested that there were divergent lineages of B. boreas uncovered in the phylogenetic analyses from both states, but morphology was not investigated. Research examining the historic hydrological connections into Railroad Valley suggests a connection from the now disjunct White River to the Colorado River, which may have provided the corridors necessary for toad dispersal into the southern Great Basin (Noles, 2010). This may provide some explanation for the high level of divergence exhibited in B. nevadensis and B. monfontanus to Great Basin B. boreas and allied taxa (Table 5), and elucidate the genetic link to western toads outside the region. Additional study is required to advance our understanding of the fine-scale relationships among the toads within this species complex and work to identify cryptic taxa that may still remain undetected under the broad range of the Western Toad.
As species are fundamental to biological study, accurate taxonomy is critical for proper evaluation of diversity and conservation implementation (Bickford et al., 2007). Bufo nevadensis and B. monfontanus represent novel species concealed within a widely distributed nominal species (Bickford et al., 2007). These species went undetected until our recent molecular study, which demonstrates the increased risk of extinction for cryptic species that are rare due to inadequately resolved taxonomy. With the exception of widely distributed B. boreas, the congeneric taxa of the B. boreas species complex are currently Threatened (IUCN, 2015) and themselves restricted to extremely small ranges (Fig. 5B). Both B. nevadensis and B. monfontanus inhabit severely small geographic distributions and join this complex as critically imperiled new members of the B. boreas species complex. This study exposes the link between taxonomic crypsis and high extinction risk, which can have profound consequences to the preservation of biodiversity due to inaccurate taxonomy that may result in improper conservation initiatives. Inadequate conservation can result in serious ramifications, which further endanger these rare, endemic toads reliant on rare, isolated, and fragile wetland habitat open to mismanagement. Unknown population sizes, limited knowledge of life histories, and small geographic ranges further emphasize the increased risk of extinction for both these newly discovered bufonids. Moreover, these new toad species reveal that our knowledge of North American anuran diversity remains incomplete (Bickford et al., 2007; Crawford et al., 2010), and taxonomic crypsis among frogs poses an important challenge in preservation of anuran diversity and conservation for a class experiencing global declines and extinctions (Collins and Storfer, 2003; Corn, 2005; Köhler et al., 2005; Crawford et al., 2010; Lannoo, 2012). Delimiting both B. monfontanus and B. nevadensis is the first step in refining our knowledge of the diversity within the B. boreas species complex, enriching our understanding of bufonid evolution within the Great Basin, and is the necessary first step in launching critical conservation initiatives to protect these vulnerable, rare toads.
Supplemental material is available at https://www.copeiajournal.org/ch-18-086.
We would like to extend our gratitude to US Fish and Wildlife Service (especially Bob Williams and Todd Gilmore), the Biological Resources Research Center, and the Department of Biology of the University of Nevada, Reno for providing funding for this research. Further thanks are extended for the invaluable field support from Sarah Snyder and Molly Betchel, plus field assistance from Kevin Guadalupe and Chris Burg of Nevada Department of Wildlife (NDOW #223071); laboratory support and genetic modeling from Bridgette Hagerty; Molly Stephens for bufonid tissues used to develop DNA markers; assistance in phylogenetic reconstruction from Joshua Hallas, Mo Beck, and Matthew Forister; and aid from Pete Noles and Erich Purpur on map development and design. We extend thanks to Chris Feldman for his help in preserving the holotype and the paratypic series. Additional thanks to Jens Vindum for assistance with the deposition of the valuable specimens into the Herpetology Department of the California Academy of Sciences. We were grateful for thoughtful reviews and attention to detail on early drafts from Chava Weitzman.