Patterns of gene flow and population genetic structure in the canyon treefrog, Hyla arenicolor (Cope)

Genetic and morphometric divergence in the Garnet-Throated Hummingbird Lamprolaima rhami (Aves: Trochilidae)

Cloud forests are one of the most endangered ecosystems in the Americas, as well as one of the richest in biological diversity in the world. The species inhabiting these forests are susceptible to environmental changes and characterized by high levels of geographic structure. The Garnet-Throated Hummingbird, Lamprolaima rhami, mainly inhabits cloud forests, but can also be found in other habitats. This species has a highly restricted distribution in Mesoamerica, and five disjunct regions have been delimited within the current geographic distribution of the species from Mexico to Honduras. According to variation in size and color, three subspecies have been described: L. r. rhami restricted to the Mexican highlands and Guatemala, L. r. occidentalis distributed in Guerrero (Mexico), and L. r. saturatior, distributed in the highlands from Honduras and El Salvador. We analyzed the levels of geographic structure in L. rhami and its taxonomic implications. We used mitochondrial and nuclear DNA to analyze genetic variation, demographic history, divergence times, reconstructed a multilocus phylogeny, and performed a species delimitation analyses. We also evaluated morphological variation in 208 specimens. We found high levels of genetic differentiation in three groups, and significant variation in morphological traits corresponding with the disjunct geographic populations. L. rhami presents population stability with the highest genetic variation explained by differences between populations. Divergence time estimates suggest that L. rhami split from its sister group around 10.55 million years ago, and the diversification of the complex was dated ca. 0.207 Mya. The hypotheses tested in the species delimitation analyses validated three independent lineages corresponding to three disjunct populations. This study provides evidence of genetic and/or morphometric differentiation between populations in the L. rhami complex where four separate evolutionary lineages are supported: (1) populations from the Sierra Madre Oriental and the highlands of Oaxaca (rhami), (2) populations from the highlands of Guerrero (occidentalis), (3) populations from the highlands of Chiapas and Guatemala (this is a non-previously proposed potential taxon: tacanensis), and (4) populations from the highlands of Honduras and El Salvador (saturatior). The main promoters of the geographic structure found in the L. rhami complex are likely the Isthmus of Tehuantepec as a geographic barrier, isolation by distance resulting from habitat fragmentation, and climatic conditions during the Pleistocene.

A depth dependent transversely isotropic micromechanic model of articular cartilage

Articular cartilage owing to the variation of collagen fibers orientation through its zones has been indicated to have depth dependent mechanical properties. The aim of this study was to present an innovative micromechanics model to predict the depth dependent mechanical properties of articular cartilage as a function of collagen fibers and proteoglycan matrix mechanical properties, collagen fibers volume fraction as well as angle toward cartilage surface. The variation of collagen fibers angle toward the cartilage surface as a function of cartilage depth was computed using the micromechanics model. This function showed that the collagen fibers parallel to the cartilage surface in the superficial zone have a nonlinear angle variation in the transition zone and become perpendicular to cartilage surface in the deep zone. Depth dependent elastic modulus in perpendicular to cartilage surface plane direction was calculated using presented micromechanics model and variation function of the collagen fibers' angle. The results revealed a suitable agreement with that of the experimental measurements in different samples at different ages and races (R2=0.944). The results also showed that the elastic and aggregate modules perpendicular to the cartilage surface plane in the deep zone were 25.8 and 26.3 times higher than that of the superficial zone, respectively. These findings have implications not only for computing the depth dependent mechanical properties of any type of articular cartilage at different ages and races, but also of potential ability for developing a depth dependent transversely isotropic biphasic model to predict the accurate mechanical behavior of articular cartilage.

The evolution and age of populations of Scaphinotus petersi Roeschke on Arizona Sky Islands (Coleoptera, Carabidae, Cychrini)

Populations of the ground beetle Scaphinotus petersi are isolated in subalpine conifer forest habitats on mountain ranges or Sky Islands in southeastern Arizona. Previous work on this species has suggested these populations have been isolated since the last post-glacial maximum times as warming caused this cool adapted species to retreat to high elevations. To test this hypothesis, we inferred the phylogeny from mitochondrial DNA sequence data from several Arizona Sky Island populations of Scaphinotus petersi and estimated the divergence time of the currently isolated populations. We found two major clades of Scaphinotus petersi, an eastern clade and a western group. Our results indicated most mountain ranges form clades except the Huachucas, which are polyphyletic and the Santa Catalinas, which are paraphyletic. We estimated the Pinaleño population is much older than the last glacial maximum, but the Huachuca and Pinal populations may have been fragmented from the Santa Catalina population since the post-glacial maximum times.

Geographic variation in advertisement calls in a tree frog species: gene flow and selection hypotheses

BACKGROUND

In a species with a large distribution relative to its dispersal capacity, geographic variation in traits may be explained by gene flow, selection, or the combined effects of both. Studies of genetic diversity using neutral molecular markers show that patterns of isolation by distance (IBD) or barrier effect may be evident for geographic variation at the molecular level in amphibian species. However, selective factors such as habitat, predator, or interspecific interactions may be critical for geographic variation in sexual traits. We studied geographic variation in advertisement calls in the tree frog Hyla japonica to understand patterns of variation in these traits across Korea and provide clues about the underlying forces for variation.

METHODOLOGY

We recorded calls of H. japonica in three breeding seasons from 17 localities including localities in remote Jeju Island. Call characters analyzed were note repetition rate (NRR), note duration (ND), and dominant frequency (DF), along with snout-to-vent length.

RESULTS

The findings of a barrier effect on DF and a longitudinal variation in NRR seemed to suggest that an open sea between the mainland and Jeju Island and mountain ranges dominated by the north-south Taebaek Mountains were related to geographic variation in call characters. Furthermore, there was a pattern of IBD in mitochondrial DNA sequences. However, no comparable pattern of IBD was found between geographic distance and call characters. We also failed to detect any effects of habitat or interspecific interaction on call characters.

CONCLUSIONS

Geographic variations in call characters as well as mitochondrial DNA sequences were largely stratified by geographic factors such as distance and barriers in Korean populations of H. japonica. Although we did not detect effects of habitat or interspecific interaction, some other selective factors such as sexual selection might still be operating on call characters in conjunction with restricted gene flow.

Elucidation of cryptic diversity in a widespread nearctic treefrog reveals episodes of mitochondrial gene capture as frogs diversified across a dynamic landscape

We investigate the evolutionary history of the wide-ranging Nearctic treefrog Hyla arenicolor through the integration of extensive range-wide sampling, phylogenetic analyses of multilocus genetic data, and divergence dating. Previous phylogeographic studies of this frog documented a potential signature of introgressive hybridization from an ecologically and morphologically divergent sister species. Based on our Bayesian phylogenetic analyses of mitochondrial DNA, we inferred strong phylogeographic structure in H. arenicolor as indicated by seven well-supported clades, five of which correspond to well-defined biogeographic regions. Clades from the Balsas Basin and southwestern Central Mexican Plateau in Mexico, and the Grand Canyon of Arizona, group with the morphologically, behaviorally, and ecologically divergent mountain treefrogs in the H. eximia group, rendering H. arenicolor as paraphyletic. The phylogenetic position of at least two of these three H. arenicolor clades within the H. eximia group, however, is most likely the result of several episodes of introgressive hybridization and subsequent mitochondrial gene capture separated in time and space, as supported by evidence from the nuclear genes. Hyla arenicolor from the Balsas Basin appear to be deeply divergent from other H. arenicolor and represent a distinctly different species. Results suggests that introgressive hybridization events, both ancient and contemporary, coupled with late Neogene vicariance and Pleistocene climate-driven range shifts, have all played a role in the historical diversification of H. arenicolor.

Genetic diversity and structure in Canadian northern leopard frog (Rana pipiens) populations: implications for reintroduction programs

The northern leopard frog ( Rana pipiens Schreber, 1782) underwent a large decline in the western portion of its range and only occurs in 20% of historically occupied sites in Alberta. Its absence may reflect an inability to disperse to these sites because of habitat fragmentation, and human-mediated translocation has been proposed. In this study, we used three criteria to examine the genetic suitability of potential translocation sources: diversity, similarity to area of reintroduction, and evolutionary history. We genotyped 187 samples and sequenced 812 bp of the mitochondrial NADH dehydrogenase 1 gene from 14 Canadian northern leopard frog populations. Nuclear and mitochondrial diversity were highest in Manitoba and western Ontario and declined westward. There was no significant relationship between genetic and geographic distance, suggesting that genetic drift is a driving force affecting the genetic relationships between populations. Regions separated by more than ~50 km were quite differentiated. Therefore, source populations similar to the original inhabitants of an area for reintroduction may be uncommon. Mitochondrial analyses revealed that all populations share a close evolutionary history, belonging to the western haplotype group. While genetic criteria support the use of Manitoba and Ontario as sources, the desirability of environmental similarity to the reintroduction site suggests that ecologically exchangeable Alberta populations should also be considered.

RAPD variation and population genetic structure of Physalaemus cuvieri (Anura: Leptodactylidae) in Central Brazil

Studies about the organization of the genetic variability and population structure in natural populations are used either to understand microevolutionary processes or the effects of isolation by human-inducted landscape modifications. In this paper, we analyzed patterns of genetic population structure using 126 RAPD loci scored for 214 individuals of Physalaemus cuvieri, sampled from 18 local populations. Around 97% of these loci were polymorphic. The among-population variation component (Phi(ST)) obtained by AMOVA was equal to 0.101 and theta B obtained using a Bayesian approach for dominant markers was 0.103. Genetic divergence, analyzed by Mantel spatial correlogram, revealed only a short-distance significant correlation between genetic and geographic distances. This is expected if low levels of population differentiation, due to high abundance buffering the effect of stochastic processes, are combined with low spatially restricted gene flow. Although this may be consistent with the current knowledge of species' biology, the spatial distribution of local populations observed in this study also suggest that, at least in part, recent human occupation and habitat fragmentation may also explain part of the interpopulational component of the genetic variation.

Phylogeography of the longhorn cactus beetle Moneilema appressum LeConte (Coleoptera: Cerambycidae): was the differentiation of the Madrean sky islands driven by Pleistocene climate changes?

Although it has been suggested that Pleistocene climate changes drove population differentiation and speciation in many groups of organisms, population genetic evidence in support of this scenario has been ambiguous, and it has often been difficult to distinguish putative vicariance from simple isolation by distance. The sky island communities of the American Southwest present an ideal system in which to compare late Pleistocene range fragmentations documented by palaeoenvironmental studies with population genetic data from organisms within these communities. In order to elucidate the impact of Pleistocene climate fluctuations on these environments, biogeographic patterns in the flightless longhorn cactus beetle, Moneilema appressum were examined using mitochondrial DNA sequence data. Gene tree relationships between haplotypes were inferred using parsimony, maximum-likelihood, and Bayesian analysis. Nested clade analysis, Mantel tests, and coalescent modelling were employed to examine alternative biogeographic scenarios, and to test the hypothesis that Pleistocene climate changes drove population differentiation in this species. The program mdiv was used to estimate migration and divergence times between populations, and to measure the statistical support for isolation over ongoing migration. These analyses showed significant geographic structure in genetic relationships, and implicated topography as a key determinant of isolation. However, although the coalescent analyses suggested that a history of past habitat fragmentation underlies the observed geographic patterns, the nested clade analysis indicated that the pattern was consistent with isolation by distance. Estimated divergence times indicated that range fragmentation in M. appressum is considerably older than the end of the most recent glacial, but coincided with earlier interglacial warming events and with documented range expansions in other, desert-dwelling species of Moneilema.

Controlling for the effects of history and nonequilibrium conditions in gene flow estimates in northern bullfrog (Rana catesbeiana) populations

Nonequilibrium conditions due to either allopatry followed by secondary contact or recent range expansion can confound measurements of gene flow among populations in previously glaciated regions. We determined the scale at which gene flow can be estimated among breeding aggregations of bullfrogs (Rana catesbeiana) at the northern limit of their range in Ontario, Canada, using seven highly polymorphic DNA microsatellite loci. We first identified breeding aggregations that likely share a common history, determined from the pattern of allelic richness, factorial correspondence analysis, and a previously published mtDNA phylogeography, and then tested for regional equilibrium by evaluating the association between pairwise F(ST) and geographic distance. Regional breeding aggregations in eastern Ontario separated by <100 km were determined to be at or near equilibrium. High levels of gene flow were measured using traditional F-statistics and likelihood estimates of Nm. Similarly high levels of recent migration (past one to three generations) were estimated among the breeding aggregations using nonequilibrium methods. We also show that, in many cases, breeding aggregations separated by up to tens of kilometers are not genetically distinct enough to be considered separate genetic populations. These results have important implications both for the identification of independent "populations" and in assessing the effect of scale in detecting patterns of genetic equilibrium and gene flow.

High mitochondrial diversity within and among populations of Malagasy poison frogs

The diurnal, brightly colored, and toxic frogs of the genus Mantella are among the most prominent representatives of the endemic anuran fauna of Madagascar. Especially three closely related species, M. aurantiaca, M. crocea, and M. milotympanum, are intensively collected for the pet trade, although basic data on their natural history and genetic diversity are still lacking. Our phylogenetic analyses based on 2.8 kbp of partial 16S rRNA, 12S rRNA, cytochrome b, and rhodopsin DNA sequences confirmed that these species belong to one of the five major clades in Mantella, the M. madagascariensis group. A haplotype network constructed using 830 bp of cytochrome b in 49 individuals from seven populations revealed that M. milotympanum and M. crocea have largely similar haplotypes sharing, confirming doubts about the species validity of M. milotympanum and indicating independent evolution of bright orange pattern in M. milotympanum and M. aurantiaca. Further, clustering of four individuals of M. aurantiaca from Andranomena with M. crocea suggests incomplete lineage sorting or introgression resulting from secondary contact of refugial populations. AMOVA confirmed significant intrapopulation nucleotide diversity (>20%). These diversity patterns and our field observations indicate relatively large population sizes. Hence, overcollecting is probably a minor problem and conservation efforts should rather focus on saving some large populations from habitat destruction through logging and forest fires.

Fine‐scale genetic pattern and evidence for sex‐biased dispersal in the túngara frog,Physalaemus pustulosus

Túngara frogs (Physalaemus pustulosus) are a model system for sexual selection and communication. Population dynamics and gene flow are of major interest in this species because they influence speciation processes and microevolution, and could consequently provide a deeper understanding of the evolutionary processes involved in mate recognition. Although earlier studies have documented genetic variation across the species' range, attempts to investigate dispersal on a local level have been limited to mark-recapture studies. These behavioural studies indicated high mobility at a scale of several hundred metres. In this study we used seven highly polymorphic microsatellite loci to investigate fine-scaled genetic variation in the túngara frog. We analysed the influence of geographical distance on observed genetic patterns, examined the influence of a river on gene flow, and tested for sex-biased dispersal. Data for 668 individuals from 17 populations ranging in distance from 0.26 to 11.8 km revealed significant levels of genetic differentiation among populations. Genetic differentiation was significantly correlated with geographic distance. A river acted as an efficient barrier to gene flow. Several tests of sex-biased dispersal were conducted. Most of them showed no difference between the sexes, but variance of Assignment Indices exhibited a statistically significant male bias in dispersal.

Does advertisement call variation coincide with genetic variation in the genetically diverse frog taxon currently known as Leptodactylus fuscus (Amphibia: Leptodactylidae)?

The frog Leptodactylus fuscus is found throughout much of South America in open and disturbed habitats. Previous study of genetic differentiation in L. fuscus demonstrated that there was lack of genetic exchange among population units consistent with multiple species, rather than a single species. We examine advertisement vocalizations of L. fuscus to determine whether call variation coincides with genetic differentiation. Calls were analyzed for 32 individual frogs from 25 localities throughout the distributional range of L. fuscus. Although there is variation in calls among geographic samples, call variation is not concordant with genetic variation or geographic distance and the call variation observed is less than that typically found among other closely related species of Leptodactylus. This study is an example of the rare pattern of strong genetic differentiation unaccompanied by salient differences in advertisement calls. The relative infrequency of this pattern as currently understood may only reflect the lack of detailed analyses of genetic and acoustic differentiation within population systems currently understood as single species with substantial geographic distributions.

Population genetics of shortnose sturgeon Acipenser brevirostrum based on mitochondrial DNA control region sequences

Shortnose sturgeon is an anadromous North American acipenserid that since 1973 has been designated as federally endangered in US waters. Historically, shortnose sturgeon occurred in as many as 19 rivers from the St. John River, NB, to the St. Johns River, FL, and these populations ranged in census size from 10(1) to 10(4), but little is known of their population structure or levels of gene flow. We used the polymerase chain reaction (PCR) and direct sequence analysis of a 440 bp portion of the mitochondrial DNA (mtDNA) control region to address these issues and to compare haplotype diversity with population size. Twenty-nine mtDNA nucleotide-substitution haplotypes were revealed among 275 specimens from 11 rivers and estuaries. Additionally, mtDNA length variation (6 haplotypes) and heteroplasmy (2-5 haplotypes for some individuals) were found. Significant genetic differentiation (P < 0.05) of mtDNA nucleotide-substitution haplotypes and length-variant haplotypes was observed among populations from all rivers and estuaries surveyed with the exception of the Delaware River and Chesapeake Bay collections. Significant haplotype differentiation was even observed between samples from two rivers (Kennebec and Androscoggin) within the Kennebec River drainage. The absence of haplotype frequency differences between samples from the Delaware River and Chesapeake Bay reflects a probable current absence of spawning within the Chesapeake Bay system and immigration of fish from the adjoining Delaware River. Haplotypic diversity indices ranged between 0.817 and 0.641; no relationship (P > 0.05) was found between haplotype diversity and census size. Gene flow estimates among populations were often low (< 2.0), but were generally higher at the latitudinal extremes of their distribution. A moderate level of haplotype diversity and a high percentage (37.9%) of haplotypes unique to the northern, once-glaciated region suggests that northern populations survived the Pleistocene in a northern refugium. Analysis of molecular variance best supported a five-region hierarchical grouping of populations, but our results indicate that in almost all cases populations of shortnose sturgeon should be managed as separate units.

Rapid recovery of genetic diversity of stomatopod populations on Krakatau: temporal and spatial scales of marine larval dispersal

Although the recovery of terrestrial communities shattered by the massive eruption of Krakatau in 1883 has been well chronicled, the fate of marine populations has been largely ignored. We examined patterns of genetic diversity in populations of two coral reef-dwelling mantis shrimp, Haptosquilla pulchella and Haptosquilla glyptocercus (Stomatopoda: Protosquillidae), on the islands of Anak Krakatau and Rakata. Genetic surveys of mitochondrial cytochrome oxidase c (subunit 1) in these populations revealed remarkably high levels of haplotypic and nucleotide diversity that were comparable with undisturbed populations throughout the Indo-Pacific. Recolonization and rapid recovery of genetic diversity in the Krakatau populations indicates that larval dispersal from multiple and diverse source populations contributes substantially to the demographics of local populations over intermediate temporal (tens to hundreds of years) and spatial scales (tens to hundreds of kilometres). Natural experiments such as Krakatau provide an excellent mechanism to investigate marine larval dispersal and connectivity. Results from stomatopods indicate that marine reserves should be spaced no more than 50-100 km apart to facilitate ecological connectivity via larval dispersal.

Phylogeography and conservation genetics of the Columbia spotted frog (Rana luteiventris; Amphibia, Ranidae)

The Columbia spotted frog (Rana luteiventris) has a widespread distribution in western Canada and the western US, although the southern reach of its range is highly fragmented into several isolated populations. Threats from various factors have raised concerns regarding the long-term survival of many small, isolated populations. Here, we report a study designed to determine the phylogeographic and conservation genetic parameters of R. luteiventris in the western US. Mitochondrial DNA (mtDNA) sequences were examined for phylogeographic structuring using phylogenetic reconstruction methods, coupled with networking and nested clade analyses. These methods permitted a distinction to be made between historic and demographic forces acting to generate geographical patterning of genetic variation. Phylogenetic analysis revealed four geographically correlated monophyletic clades. Three of these clades correspond to well-defined, nonoverlapping geographical locations in the fragmented portion of the range. The other is comprised of all samples collected from the contiguous range and includes one isolate from northern Wyoming. Networking and nested clade analyses confirmed these results and revealed that historical processes, such as range expansion and vicariance, rather than recurrent gene flow are likely responsible for observed patterns of genetic variation. A measure of genetic variation (theta = 4N(e)mu) revealed that R. luteiventris populations in Utah have a relatively low amount of genetic variation compared with populations in the continuous portion of the range.

Microsatellite variation and fine-scale population structure in the wood frog (Rana sylvatica)

We investigated genetic population structure in wood frogs (Rana sylvatica) from a series of Prairie Pothole wetlands in the northern Great Plains. Amphibians are often thought to exist in demographic metapopulations, which require some movement between populations, yet genetic studies have revealed strong subdivision among populations, even at relatively fine scales (several km). Wood frogs are highly philopatric and studies of dispersal suggest that they may exhibit subdivision on a scale of approximately 1-2 km. We used microsatellites to examine population structure among 11 breeding assemblages separated by as little as 50 m up to approximately 5.5 km, plus one population separated from the others by 20 km. We found evidence for differentiation at the largest distances we examined and among a few neighbouring ponds, but most populations were strikingly similar in allele frequencies, suggesting high gene flow among all but the most distant populations. We hypothesize that the few significant differences among neighbouring populations at the finest scale may be a transient effect of extinction-recolonization founder events, driven by periodic drying of wetlands in this hydrologically dynamic landscape.

Isolation by distance in the Atlantic cod, Gadus morhua, at large and small geographic scales

Genetic isolation by distance (IBD) has rarely been described in marine species with high potential for dispersal at both the larval and adult life-history stages. Here, we report significant relationships between inferred levels of gene flow and geographic distance in the Atlantic cod, Gadus morhua, at 10 nuclear restriction-fragment-length-polymorphism (RFLP) loci at small regional scales in the western north Atlantic region (< 1,600 km) that mirror those previously detected over its entire geographic range (up to 7,300 km). Highly significant allele frequency differences were observed among eight northwestern Atlantic populations, although the mean FST for all 10 loci was only 0.014. Despite this weak population structuring, the distance separating populations explained between 54% and 62% of the variation in gene flow depending on whether nine or 10 loci were used to estimate Nm. Across the species' entire geographic range, highly significant differences were observed among six regional populations at nine of the 10 loci (mean FST = 0.068) and seven loci exhibited significant negative relationships between gene flow and distance. At this large geographic scale, natural selection acting in the vicinity of one RFLP locus (GM798) had a significant effect on the correlation between gene flow and distance, and eliminating it from the analysis caused the coefficient of determination to increase from 17% to 62%. The role of vicariance was assessed by sequentially removing populations from the analysis and was found to play a minor role in contributing to the relationship between gene flow and distance at either geographic scale. The correlation between gene flow and distance detected in G. morhua at small and large spatial scales suggests that dispersal distances and effective population sizes are much smaller than predicted for the species and that the recent age of populations, rather than extensive gene flow, may be responsible for its weak population structure. Our results suggest that interpreting limited genetic differences among populations as reflecting high levels of ongoing gene flow should be made with caution.

Phylogeography of the jumping spider Habronattus pugillis (araneae: salticidae): recent vicariance of sky island populations?

In island systems with diverging populations, the history of island formation and genealogical estimates of divergence dates can be mutually informative. In the "sky islands" of southeastern Arizona, climate-induced contraction of woodlands appears to have fragmented populations of woodland-dwelling species onto disjunct mountain ranges. Montane populations of the jumping spider, Habronattus pugillis, display striking amounts of phenotypic divergence among ranges. Paleoclimatic estimates date woodland fragmentation at approximately 10,000 years ago, suggesting that phenotypic divergence has been extraordinarily rapid in these spiders. This phylogeographic study of populations of H. pugillis attempts to clarify the species' history of isolation and divergence and to address the suitability of available paleoclimatic data for dating divergences among populations of the region's woodland-dwelling organisms. Mitochondrial sequence data of spiders from 13 mountain ranges was used to reconstruct genealogical relationships. Gene trees show that small mountain ranges tend to have populations whose sequences form monophyletic groups, whereas larger ranges do not. Paraphyly among genes from larger ranges could result from either recent migration or incomplete lineage sorting. I use phylogenetic and geographic information to test these alternatives, and conclude that incomplete lineage sorting best explains the observed paraphyly. Gene trees are concordant with some of the predictions of vegetation history generated by examination of topography. Dates estimated for divergence of populations vary from 30,000 years to more than 2 million years ago, suggesting multiple vicariance events that are older than would be inferred from paleoclimatic studies. These findings illustrate that use of any single paleontological dataset to calibrate molecular clocks can potentially greatly underestimate actual divergence times.

Phylogeography of the canyon treefrog, Hyla arenicolor (Cope) based on mitochondrial DNA sequence data

Patterns of phylogeography and gene flow were examined in the canyon treefrog, Hyla arenicolor. A total of 973 bp of mitochondrial cytochrome b sequence data were obtained for 65 individuals from 53 populations, yielding 50 unique haplotypes. Interpopulation sequence variation ranged from 0 to 13.7%. Phylogenetic analysis revealed three deeply divergent mtDNA lineages. These three Clades were mapped onto geography and found to represent completely concordant, nonoverlapping, geographical regions. Levels of sequence divergence between the three Clades were equal to or greater than published levels of divergence found in other vertebrate species and genera. Furthermore, one Clade of H. arenicolor was found to be more closely related to the outgroup H. eximia than to other H. arenicolor, suggesting that the taxonomy of this species may require revision.