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

The geography of evolutionary divergence in the highly endemic avifauna from the Sierra Madre del Sur, Mexico

Background

Mesoamerica is a remarkable region with a high geological and ecological complexity. Within northern Mesoamerica, the biotic province of the Sierra Madre del Sur (SMS) in southwestern Mexico harbors exceptionally high avian endemism and diversity. Herein, we searched for spatially and temporally concordant phylogeographic patterns, in four bird genera from three distinct avian orders co-distributed across Mesoamerica and investigated their causes through hypothesis testing regarding historical processes. Selected species include endemic and differentiated populations across the montane forests of Mesoamerica, and particularly within the SMS.

Results

We gathered mitochondrial DNA sequences for at least one locus from 177 individuals across all species. We assessed genetic structure, demographic history, and defined a framework for the coalescent simulations used in biogeographic hypothesis testing temporal and spatial co-variance. Our analyses suggested shared phylogeographic breaks in areas corresponding to the SMS populations, and between the main montane systems in Mesoamerica, with the Central Valley of Oaxaca and the Nicaragua Depression being the most frequently shared breaks among analyzed taxa. Nevertheless, dating analyses and divergence patterns observed were consistent with the hypothesis of broad vicariance across Mesoamerica derived from mechanisms operating at distinct times across taxa in the SMS.

Conclusions

Our study provides a framework for understanding the evolutionary origins and historical factors enhancing speciation in well-defined regions within Mesoamerica, indicating that the evolutionary history of extant biota inhabiting montane forests is complex and often idiosyncratic.

Climate change, extinction, and Sky Island biogeography in a montane lizard

Around the world, many species are confined to "Sky Islands," with different populations in isolated patches of montane habitat. How does this pattern arise? One scenario is that montane species were widespread in lowlands when climates were cooler, and were isolated by local extinction caused by warming conditions. This scenario implies that many montane species may be highly susceptible to anthropogenic warming. Here, we test this scenario in a montane lizard (Sceloporus jarrovii) from the Madrean Sky Islands of southeastern Arizona. We combined data from field surveys, climate, population genomics, and physiology. Overall, our results support the hypothesis that this species' current distribution is explained by local extinction caused by past climate change. However, our results for this species differ from simple expectations in several ways: (a) their absence at lower elevations is related to warm winter temperatures, not hot summer temperatures; (b) they appear to exclude a low-elevation congener from higher elevations, not the converse; (c) they are apparently absent from many climatically suitable but low mountain ranges, seemingly "pushed off the top" by climates even warmer than those today; (d) despite the potential for dispersal among ranges during recent glacial periods (~18,000 years ago), populations in different ranges diverged ~4.5-0.5 million years ago and remained largely distinct; and (e) body temperatures are inversely related to climatic temperatures among sites. These results may have implications for many other Sky Island systems. More broadly, we suggest that Sky Island species may be relevant for predicting responses to future warming.

Different patterns of colonization of Oxalis alpina in the Sky Islands of the Sonoran desert via pollen and seed flow

Historical factors such as climatic oscillations during the Pleistocene epoch have dramatically impacted species distributions. Studies of the patterns of genetic structure in angiosperm species using molecular markers with different modes of inheritance contribute to a better understanding of potential differences in colonization and patterns of gene flow via pollen and seeds. These markers may also provide insights into the evolution of reproductive systems in plants. Oxalis alpina is a tetraploid, herbaceous species inhabiting the Sky Island region of the southwestern United States and northern Mexico. Our main objective in this study was to analyze the influence of climatic oscillations on the genetic structure of O. alpina and the impact of these oscillations on the evolutionary transition from tristylous to distylous reproductive systems. We used microsatellite markers and compared our results to a previous study using chloroplast genetic markers. The phylogeographic structure inferred by both markers was different, suggesting that intrinsic characteristics including the pollination system and seed dispersal have influenced patterns of gene flow. Microsatellites exhibited low genetic structure, showed no significant association between geographic and genetic distances, and all individual genotypes were assigned to two main groups. In contrast, chloroplast markers exhibited a strong association between geographic and genetic distance, had higher levels of genetic differentiation, and were assigned to five groups. Both types of DNA markers showed evidence of a northward expansion as a consequence of climate warming occurring in the last 10,000 years. The data from both types of markers support the hypothesis for several independent transitions from tristyly to distyly.

Variation in heterostylous breeding systems in neighbouring populations of Oxalis alpina (Oxalidaceae)

The heterostylous reproductive system of Oxalis alpina in the Galiuro Mts. of Arizona was investigated using field surveys, controlled crosses in the greenhouse and measurements of reproductive morphs. Although populations in the Pinaleño Mts. to the immediate east and in the Santa Catalina Mts. to the immediate west have derived distylous reproductive systems, tristyly, the ancestral reproductive system in O. alpina, has been retained in the Galiuro Mts.

POPULATION

Tristylous incompatibility relationships in the Galiuro population are modified from the ancestral condition, with significant loss of incompatibility differentiation between stamen whorls of both short- and long-styled morphs. Morphological adjustments of anther positions in the Galiuro population of O. alpina match those expected in light of incompatibility modification, with divergence of the mid-level anthers away from the position of the mid stigmas of the mid-styled morph. The occurrence of tristyly in an area of Arizona where distyly is found in adjacent mountain ranges is particularly remarkable, and indicates both the isolation of populations restricted to the upper elevations of these mountain ranges and variation in the tempo of evolution over short geographic distances.

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.

Molecular patterns of differentiation in canyon treefrogs (Hyla arenicolor): evidence for introgressive hybridization with the Arizona treefrog (H. wrightorum) and correlations with advertisement call differences

Detection of genetic and behavioural diversity within morphologically similar species has led to the discovery of cryptic species complexes. We tested the hypothesis that US populations of the canyon treefrog (Hyla arenicolor) may consist of cryptic species by examining mate-attraction signals among three divergent clades defined by mtDNA. Using a multi-locus approach, we re-analysed phylogenetic relationships among the three clades and a closely related, but morphologically and behaviourally dissimilar species, the Arizona treefrog (H. wrightorum). We found evidence for introgression of H. wrightorum's mitochondrial genome into H. arenicolor. Additionally, the two-clade topology based on nuclear data is more congruent with patterns of call variation than the three-clade topology from the mitochondrial dataset. The magnitude of the call divergence is probably insufficient to promote isolation of the nuclear DNA-defined clades should they become sympatric, but further divergence in call properties significant in species identification could promote speciation in the future.

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.

Integrating paleoecology and genetics of bird populations in two sky island archipelagos

Background

Genetic tests of paleoecological hypotheses have been rare, partly because recent genetic divergence is difficult to detect and time. According to fossil plant data, continuous woodland in the southwestern USA and northern Mexico became fragmented during the last 10,000 years, as warming caused cool-adapted species to retreat to high elevations. Most genetic studies of resulting 'sky islands' have either failed to detect recent divergence or have found discordant evidence for ancient divergence. We test this paleoecological hypothesis for the region with intraspecific mitochondrial DNA and microsatellite data from sky-island populations of a sedentary bird, the Mexican jay (Aphelocoma ultramarina). We predicted that populations on different sky islands would share common, ancestral alleles that existed during the last glaciation, but that populations on each sky island, owing to their isolation, would contain unique variants of postglacial origin. We also predicted that divergence times estimated from corrected genetic distance and a coalescence model would post-date the last glacial maximum.

Results

Our results provide multiple independent lines of support for postglacial divergence, with the predicted pattern of shared and unique mitochondrial DNA haplotypes appearing in two independent sky-island archipelagos, and most estimates of divergence time based on corrected genetic distance post-dating the last glacial maximum. Likewise, an isolation model based on multilocus gene coalescence indicated postglacial divergence of five pairs of sky islands. In contrast to their similar recent histories, the two archipelagos had dissimilar historical patterns in that sky islands in Arizona showed evidence for older divergence, suggesting different responses to the last glaciation.

Conclusion

This study is one of the first to provide explicit support from genetic data for a postglacial divergence scenario predicted by one of the best paleoecological records in the world. Our results demonstrate that sky islands act as generators of genetic diversity at both recent and historical timescales and underscore the importance of thorough sampling and the use of loci with fast mutation rates to studies that test hypotheses concerning recent genetic divergence.

Evolution of Synthyris sect. Dissecta (Plantaginaceae) on sky islands in the Northern Rocky Mountains

We examined the phylogeography of Synthyris sect. Dissecta (Plantaginaceae), which is restricted to the Olympic (S. lanuginosa) and Northern Rocky Mountains (S. dissecta and S. canbyi), to infer effects of sky islands and modes of speciation. Sequences of cpDNA trnT-trnL and psbA-trnH intergenic spacers resolved 22 haplotypes among 185 individuals sampled from 16 populations of the three species. Gene flow in the ancestral lineage and random capture of haplotypes in species lineages of sect. Dissecta have resulted in haplotype clades that are not exclusive to species. Nested clade analysis (NCA) indicates that allopatric fragmentation separated Olympic and Northern Rocky Mountain populations, giving rise to the Olympic endemic S. lanuginosa, which is characterized by unique haplotypes consistent with long temporal isolation. Low haplotype and nucleotide diversity in S. canbyi are consistent with newly founded populations experiencing a bottleneck. Furthermore, we infer S. canbyi evolved as a peripheral isolate of S. dissecta. NCA indicated limited migration in S. dissecta with possible isolation by distance. Both isolation on interglacial sky islands and valley glaciers during at least the last glaciation limited gene flow among populations of S. dissecta in different ranges of the Northern Rocky Mountains.

Phylogeography and demographic history of Chinese black-spotted frog populations (Pelophylax nigromaculata): evidence for independent refugia expansion and secondary contact

Background

Pleistocene glaciations had considerable impact on phylogeographic patterns within and among closely related species of many vertebrates. Compared to Europe and North America, research on the phylogeography of vertebrates in East Asia, particularly in China, remains limited. The black-spotted frog (Pelophylax nigromaculata) is a widespread species in East Asia. The wide distribution of this species in China makes it an ideal model for the study of palaeoclimatic effects on vertebrates in East Asia. Our previous studies of P. nigromaculata revealed significant subdivisions between the northeast China populations and populations in other regions of the mainland. In the present study, we aim to see whether the deepest splits among lineages and perhaps subsequent genealogical divisions are temporally consistent with a Pleistocene origin and whether clade geographic distributions, with insight into expansion patterns, are similarly spatially consistent with this model.

Results

Using 1143 nucleotides of the mitochondrial cytochrome b gene from 262 individuals sampled from 28 localities, two main clades (clade A and clade B) differing by c. 7.72% sequence divergence were defined from parsimony analyses. The corresponding timing of lineage divergence, 0.92 Mya, indicates a most likely Pleistocene split. The A clade is further subdivided into two sub-clades, A1 and A2 with 1.22% sequence divergence. Nested clade phylogeographical and population demographic analyses suggested that the current distribution of this frog species was the result of range expansion from two independent refugia during the last interglacial period. We discovered a population within which haplotype lineages A and B of P. nigromaculata coexist in the Dongliao area of China by nucleotide sequences, PCR-RFLP and ISSR (inter simple sequence repeat) patterns. The ISSR result in particular supported divergence between the mitochondrial clades A and B and implied introgressive gene flow between the two divergent lineages.

Conclusion

Nested clade phylogeographical and population demographic analyses indicate that the current distribution of P. nigromaculata is the result of range expansion from two independent refugia during the last interglacial period in late Pleistocene. One refugium was in east China and the lower elevations of south-western plateau. The distribution of the other mitochondrial clade is consistent with the presence of a refugium in the Korean Peninsula. The gene flow as detected by ISSR markers suggests a range expansion of the two refugia and a secondary contact between the two highly divergent lineages in the Dongliao (DL) area of northeast China.

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.

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.

Phylogeography of the Alpine salamander, Salamandra atra (Salamandridae) and the influence of the Pleistocene climatic oscillations on population divergence

Fifty individuals of the endemic Alpine salamander, Salamandra atra, representing 13 populations throughout the range of the two currently recognized subspecies, atra and aurorae, were examined for sequence variation in a large portion (1050 bp) of the mitochondrial cytochrome b gene. We revealed a large number of mitochondrial DNA (mtDNA) haplotypes (10). Interpopulation sequence divergence was very low, ranging from 0 to 3.1%. The relationships among haplotypes were poorly resolved. The divergence time estimate between several mtDNA haplotypes suggested a pre-Pleistocene differentiation approximately 3 million years ago. Moreover, the impact of the Pleistocene glaciations on the phylogeographical patterns appears to have been secondary, although a somewhat reduced genetic variability was found in populations living in areas that were directly affected by the glaciation.

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.

Mitochondrial DNA-Based phylogeography of North American rubber boas, Charina bottae (Serpentes: Boidae)

We used 783 bp of mitochondrial DNA sequences to study the phylogeography of Charina bottae (rubber boa) in western North America, with an emphasis on populations from California (U.S.A.). Maximum-parsimony and maximum-likelihood methods identified a basal divergence within C. bottae that corresponds to southern and northern segments of its current distribution. These clades coincide with the ranges of the two recognized subspecies, C. b. umbratica in the south and C. b. bottae to the north. A subsequent cladogenetic event in the C. b. bottae clade resulted in two groupings, which we refer to as the Sierra Nevada and the Northwestern subclades, based on the geographic distribution of their constituent populations. The two subclades have completely allopatric distributions, with a genetic break in the vicinity of Lassen Volcanic National Park in northeastern California, an area that was subjected to glaciation during the Pleistocene and that has been volcanically active in the past 100 years. An earlier genetic study documented fixed differences between populations of bottae and umbratica in four of seven allozymes surveyed, and despite noticeable variation and overlap in the characters that define C. b. bottae and C. b. umbratica, the two forms still can be separated in most cases using a suite of morphological traits. All available evidence thus indicates that C. b. umbratica is a genetically cohesive, allopatric taxon that is morphologically diagnosable, and we conclude that it is an independent evolutionary unit that should be recognized as a distinct species, Charina umbratica.

Mitochondrial DNA phylogeography of Acrossocheilus paradoxus (Cyprinidae) in Taiwan

Nucleotide sequences of 3' end of the cytochrome b gene, tRNA genes, D-loop control region, and the 5' end of the 12S rRNA of mitochondrial DNA (mtDNA) were used to assess the genetic and phylogeographic structure of Acrossocheilus paradoxus populations, a Cyprinidae fish of Taiwan. A hierarchical examination of populations in 12 major streams from three geographical regions using an analysis of molecular variance (AMOVA) indicates high genetic differentiation both among populations (PhiST = 0.511, P < 0.001) and among regions (PhiCT = 0.368, P < 0.001). Limited migration largely contributed to the genetic differentiation. High nucleotide diversity (1.13%) and haplotype diversity (0.80%) were detected among populations. The degree of genetic differentiation was correlated with geographical distance between populations, a result consistent with the one-dimensional stepping stone models. A neighbour-joining tree recovered by (DAMBE) supports the pattern of isolation by distance and reveals a closer relationship between populations of the central and southern regions. A minimum spanning network based on nucleotide substitutions reflected migration routes from populations of the central region to the northern and southern regions, respectively. Postglacial colonization and expansion can explain the phylogeographical pattern. Single and ancient migration events may have allowed the northern region to attain the monophyly of mtDNA alleles. In contrast, most populations within geographical regions are either paraphyletic or polyphyletic due to the relatively shorter time period for coalescence. Both low haplotype number and genetic variability suggest a bottleneck event in the Chingmei population of northern Taiwan. Based on coalescence theory, the monophyly of the Tungkang population of the southern region may be associated with a founder event.

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.

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

Patterns of gene flow and genetic structuring were examined in the canyon treefrog, Hyla arenicolor (Cope). Hierarchical analysis of genetic variation was performed on mitochondrial cytochrome b haplotypes from 323 individuals, representing 32 populations from previously described phylogeographic regions. Results from AMOVA revealed that 60.4-78.9% of the recovered genetic variation was the result of differences in the appointment of genetic variation between subdivisions of the primary phylogeographic regions. In contrast, populations only contained between 13.9 and 30.1% of the observed haplotypic variation. Gene flow estimates based on calculations of phi ST revealed moderate levels of gene flow within phylogeographic regions, but there was no evidence of gene flow between these regions, suggesting that geographical boundaries were probably important in the formation of phylogeographic structure in H. arenicolor. Phylogeographic regions exhibited very different patterns of gene flow. One region showed evidence of recent colonization. Another region exhibited very limited gene flow. Moderate to high estimates of gene flow were obtained for populations from two distinct phylogeographic regions characterized by mesic and xeric environments. Isolation by distance was observed in both regions suggesting that these regions are in genetic equilibrium. Because gene flow is extremely unlikely between the populations in the xeric region, this result is interpreted as historical gene flow. These results indicate that isolation-by-distance effects may still be observed even when population genetic structure and gene flow are the result of historical association.