Phylogeography of Diadophis punctatus: extensive lineage diversity and repeated patterns of historical demography in a trans-continental snake

Dynamic climatic oscillations during the Pleistocene had profound effects on the distributions of species across North America. Although the role of historical climate change on speciation remains controversial, the impact on genetic variation within species has been well documented. We examined mtDNA sequences from the cytochrome b gene (1117 bp) and a portion of the NADH-4 gene (659 bp) for 286 individuals of Diadophis punctatus to infer phylogeographic patterns and population structure and to examine historical demographic patterns in both glaciated and unglaciated regions of North America. We inferred 14 lineages that replace each other geographically across the United States. Several of these lineages appear to be confined to specific habitats (floodplains, grasslands, montane environments) and traverse previously identified genetic barriers for terrestrial vertebrates including the Mississippi and Apalachicola Rivers, the Appalachian Mountains, and the western continental divide. We also observed overlapping ranges between some haplotype groups and several instances of secondary contact associated with ecological transition zones in eastern South Carolina, southern Oklahoma and central California. Within the US, diversification began during the late Miocene and continued into the mid-Pleistocene, suggesting these lineages pre-dated the last glacial maximum. Coalescent and non-coalescent demographic analyses indicate that independent lineages currently occupying previously glaciated or unsuitable areas in eastern, central and western US underwent post-glacial population expansion likely from southern refugia during the late Pleistocene/early Holocene. Conversely, southern lineages display patterns consistent with long-term population stability. Such long-term persistence of genetic structure may be due to the competitive effects between lineages or ecosystem stability in more southern latitudes.     

Identification of glacial refugia in south-eastern North America by phylogeographical analyses of a forest understorey plant, Trillium cuneatum

Aim We examine several hypotheses emerging from biogeographical and fossil records regarding glacial refugia of a southern thermophilic plant species. Specifically, we investigated the glacial history and post‐glacial colonization of a forest understorey species, Trillium cuneatum. We focused on the following questions: (1) Did T. cuneatum survive the Last Glacial Maximum (LGM) in multiple refugia, and (if so) where were they located, and is the modern genetic structure congruent with the fossil record‐based reconstruction of refugia for mesic deciduous forests? (2) What are the post‐glacial colonization patterns in the present geographical range? Location South‐eastern North America. Methods We sampled 45 populations of T. cuneatum throughout its current range. We conducted phylogeographical analyses based on maternally inherited chloroplast DNA (cpDNA haplotypes) and used TCS software to reconstruct intraspecific phylogeny. Results We detected six cpDNA haplotypes, geographically highly structured into non‐overlapping areas. With one exception, none of the populations had mixed haplotype composition. TCS analysis resulted in two intraspecific cpDNA lineages, with one clade subdivided further by shallower diversification. Main conclusions Our investigation revealed that T. cuneatum survived the LGM in multiple refugia, belonging to two (western, eastern) genealogical lineages geographically structured across south‐eastern North America. The western clade is confined to the south‐western corner of T. cuneatum’s modern range along the Lower Mississippi Valley, where fossil records document a major refugium of mesic deciduous forest. For the eastern clade, modern patterns of cpDNA haplotype distribution suggest cryptic vicariance, in the form of forest contractions and subsequent expansions associated with Pleistocene glacial cycles, rather than simple southern survival and subsequent northward colonization. The north–south partitioning of cpDNA haplotypes was unexpected, suggesting that populations of this rather southern thermophilic species may have survived in more northern locations than initially expected based on LGM climate reconstruction, and that the Appalachian Mountains functioned as a barrier to the dispersal of propagules originating in more southern refugia. Furthermore, our results reveal south‐west to north‐east directionality in historical migration through the Valley and Ridge region of north‐west Georgia.     

The demographic history of the New Zealand short-tailed bat Mystacina tuberculata inferred from modified control region sequences

Short-tailed bats Mystacina tuberculata were widespread throughout the forest that dominated prehuman New Zealand, but extensive deforestation has restricted them to scattered populations in forest fragments. In a previous study, the species' intraspecific phylogeny was investigated using multiple mitochondrial gene sequences. Six phylogroups were identified with estimated divergences of 0.93-0.68 Ma. In the current study, the phylogeographical structure and demographic history of the phylogroups were investigated using control region sequences modified by removing homoplasic sites. Phylogeographical structure in the North Island was generally consistent with an isolation-by-distance dispersal model. Coalescent-based analyses (i.e. mismatch distributions, skyline plots, lineage dispersal analysis and nested clade analysis) indicated that the three phylogroups found in central and southern North Island expanded before the last glacial maximum, presumably during interstadials when Nothofagus forest was most extensive. Genetic structure within a central North Island hybrid zone was consistent with range expansion from separate refugia following reforestation after catastrophic volcanic eruptions. Phylogeographical structure in the South Island was consistent with southern populations originating during rapid southward range expansion from refugia in northern South Island following postglacial reforestation of the South Island 10-9 kya.     

Chloroplast DNA intraspecific phylogeography of plants from the Pacific Northwest of North America

Molecular studies of plants from the Pacific Northwest of North America suggest a recurrent pattern of genetic differentiation and geographic structuring. In each of five angiosperms and one fern species representing diverse life histories, cpDNA data indicate two clades of populations that are geographically structured. A northern group comprises populations from Alaska to central or southern Oregon, whereas populations from central Oregon southward to northern California form a southern group. In several of these species, a few populations having southern genotypes may have survived in glacial refugia further north in the Olympic Peninsula, Queen Charlotte Islands, and Prince of Wales Island. Allozyme data reveal a similar pattern of differentiation in several other plants from the Pacific Northwest. North-south partitioning of genotypes has also been reported for several animal species from this region. On a broader geographic scale, northsouth partitioning of genotypes has also been observed in other plants from western North America having a variety of geographic distributions. Some species also display a reduction of genetic variability in the northern portion of their range compared to the south. The data suggest strongly that past glaciation profoundly influenced the genetic architecture of the flora and fauna of the Pacific Northwest. Two alternative hypotheses are advanced to explain the geographic structuring of genotypes. First, past glaciation may have created discontinuities in the geographic distributions of plant species, with populations surviving in several well-isolated northern and southern refugia. Following glaciation, migration of genetically differentiated, once-isolated populations resulted in the formation of a continuous geographic distribution with a major genetic discontinuity. Alternatively, plants survived and subsequently migrated northward from a southern refugium, and a genotype became fixed in one or a few populations at the leading edge of recolonization. Subsequent long-distance dispersal from this leading edge resulted in a relatively uniform northern genotype that differs from the southern genotype(s). Whatever the underlying mechanism, Pleistocence glaciation may have molded the intraspecific genetic architecture of both plants and animals from the Pacific Northwest in a geographically similar manner. Future studies should seek to obtain a comprehensive phylogeography for regions that includes a diversity of both plants and animals.Dedicated to emer. Univ.-Prof. DrFriedrich Ehrendorfer on the occasion of his 70th birthday     

Phylogeography of Cavernularia hultenii: evidence of slow genetic drift in a widely disjunct lichen

Population structure and history is poorly known in most lichenized ascomycetes. Many species display large-scale infraspecific disjunctions, which have been explained alternately by range fragmentation in species of high age and widespread long-distance dispersal. Using the lichen Cavernularia hultenii, which is widely disjunct across North America and Europe, Pleistocene and Holocene population history was inferred. The internal transcribed spacer (ITS) and part of the the intergenic spacer (IGS) region of the nuclear ribosomal DNA were sequenced in 300 individuals representing 62 populations across the range of the species. While four ancestral haplotypes are found in all areas, none of the observed tip haplotypes is present in more than one of the three part ranges. Although this is evidence for a past fragmentation event, nested clade analysis (NCA) remains equivocal in the choice between allopatric fragmentation and long-distance dispersal. Mismatch distributions indicate exponential population growth, probably during postglacial invasion of C. hultenii into formerly glaciated areas of western North America. The presence of one southern and at least one northern glacial refugium in South Central Alaska is inferred. Evidence for another refugium in the Queen Charlotte Islands or Alexander Archipelago is inconclusive because of sparse sampling. However, a range expansion was not confirmed unambiguously by NCA. The limited power of NCA to infer past range fragmentations and expansions is due apparently to the shallow haplotype network and widespread ancestral haplotypes. This can be explained by slow genetic drift causing incomplete removal of ancestral haplotypes from the postfragmentation and postexpansion areas.     

Chloroplast DNA phylogeography of a distylous shrub (Palicourea padifolia, Rubiaceae) reveals past fragmentation and demographic expansion in Mexican cloud forests

Several phylogeographic studies in northern Mesoamerica have examined the influence of Pleistocene glaciations on the genetic structure of temperate tree species with their southern limit by the contact zone between species otherwise characteristic of North or South America, but few have featured plant species that presumably colonized northern Mesoamerica from South America. A phylogeographical study of Palicourea padifolia, a fleshy-fruited, bird dispersed distylous shrub, was conducted to investigate genetic variation at two chloroplast regions (trnS-trnG and rpl32-trnL) across cloud forest areas to determine if such patterns are consistent with the presence of Pleistocene refugia and/or with the historical fragmentation of the Mexican cloud forests. We conducted population and spatial genetic analyses as well as phylogenetic and isolation with migration analyses on 122 individuals from 22 populations comprising the distribution of P. padifolia in Mexico to gain insight of the evolutionary history of these populations. Twenty-six haplotypes were identified after sequencing 1389 bp of chloroplast DNA. These haplotypes showed phylogeographic structure (N(ST) = 0.508, G(ST) = 0.337, N(ST) > G(ST), P < 0.05), including a phylogeographic break at the Isthmus of Tehuantepec, with private haplotypes at either side of the isthmus, and a divergence time of the split in the absence of gene flow dating back c. 309,000-103,000 years ago. The patterns of geographic structure found in this study are consistent with past fragmentation and demographic range expansion, supporting the role of the Isthmus of Tehuantepec as a biogeographical barrier in the dispersal of P. padifolia. Our data suggest that P. padifolia populations were isolated throughout glacial cycles by the Isthmus of Tehuantepec, accumulating genetic differences due to the lack of migration across the isthmus in either direction, but the results of our study are not consistent with the existence of the previously proposed Pleistocene refugia for rain forest plant species in the region.     

New microsatellite markers for examining genetic variation in peripheral and core populations of the Coastal Giant Salamander (Dicamptodon tenebrosus)

The Coastal Giant Salamander (Dicamptodon tenebrosus) is classified as threatened at the northern periphery of its range in British Columbia (BC), Canada, primarily due to forestry practices and habitat fragmentation. Characterising dispersal behaviour and population connectivity is therefore a priority for this region, while genetic differentiation in core versus peripheral locations remains unstudied in this wide-ranging species. We present seven new polymorphic microsatellite markers for use in population genetic analyses of D. tenebrosus. We examine locus characteristics and genetic variation in 12 streams at the species' northern range limit in BC, and within two regions representing sub-peripheral (North Cascades) and core localities (South Cascades) in Washington State, United States. In BC, the number of alleles per locus ranged from 2-5 and observed heterozygosity ranged from 0.044-0.825. Genetic differentiation was highest between BC and the South Cascades, and intermediate between BC and the North Cascades. Across loci, mean allelic richness was similar across regions, while private allelic richness was highest in the core locality (corrected for sample size). These new microsatellite loci will be a valuable addition to existing markers for detailed landscape and population genetic analyses of D. tenebrosus across its range.     

Population genetic subdivision in the New Zealand greenshell mussel (Perna canaliculus) inferred from single-strand conformation polymorphism analysis of mitochondrial DNA

Single-strand conformation polymorphism (SSCP) analysis of the NADH IV region of the mitochondrial DNA (mtDNA) molecule in greenshell mussels (Perna canaliculus) indicated strong population genetic structuring in this endemic New Zealand species. A northern and a southern group were differentiated by frequency shifts in common haplotypes and by the occurrence of a unique southern haplotype at approximately 20% frequency. This split occurred south of Cook Strait (the body of water between the North and the South Island) at approximately 42 degrees S latitude. Northern populations were less genetically diverse than southern populations and mussels from the west coast of the South Island were most distinct from northern mussels. We hypothesize that the unique haplotype VIII originated in the lower South Island, and that its spread northwards was obstructed by the opening of Cook Strait approximately 15 000-16 000 years ago and the subsequent establishment of present-day surface water circulation patterns in Greater Cook Strait. We suggest that present-day strong tidal flows and turbulent mixing of water masses in Cook Strait, and intense up-welling on the east and west coasts in this region, represent a barrier to gene flow between mussels located in the North Island and northern South Island vs. mussels in most of the South Island and Stewart Island.     

Genetic differentiation in pikeperch (Sander lucioperca): the relative importance of gene flow, drift and common history

To investigate the relative importance of homogenizing factors, such as gene flow, and diversifying factors, such as drift, genetic variation in pikeperch ( Sander lucioperca ) in two Fennoscandian regions (North and South) was analysed with microsatellites. Allelic richness and the degree of differentiation were significantly higher in the North ( F _ST= 0·20) than in the South ( F _ST= 0·064). In northern areas, assignments of genotypes were almost exclusively to the population of origin, but in southern areas, the proportion of correct assignments was significantly lower. Most samples exhibited significant heterozygote deficits, and the level of relatedness was higher than expected from randomness. These combined results suggest that there has been more gene flow between populations in southern areas than in northern areas, where the importance of genetic drift has been greater. Effective population sizes were small ( c. 100) and did not differ between areas. The effect of a common history appears minor, and thus processes such as genetic drift and gene flow have been more influential in shaping the patterns of genetic diversity in this species.     

Geological barriers and restricted gene flow in the holarctic skipper Hesperia comma (Hesperiidae)

Patterns of genetic variation within a species may be a consequence of historical factors, such as past fragmentation, as well as current barriers to gene flow. Using sequence data from the mitochondrial cytochrome oxidase subunit II region (COII) and the nuclear gene wingless, we conducted a phylogeographical study of the holarctic skipper Hesperia comma to elucidate patterns of genetic diversity and to infer historical and contemporary processes maintaining genetic variation. One hundred and fifty-one individuals were sampled from throughout North America and Eurasia, focusing on California and adjacent regions in the western United States where morphological diversity is highest compared to the rest of the range. Analyses of sequence data obtained from both genes revealed a well-supported division between the Old and New World. Within western North America, wingless shows little geographical structure, while a hierarchical analysis of genetic diversity of COII sequences indicates three major clades: a western clade in Oregon and Northern California, an eastern clade including the Great Basin, Rocky Mountains and British Columbia, and a third clade in southern California. The Sierra Nevada and the Transverse Ranges appear to be the major barriers to gene flow for H. comma in the western United States. Relatively reduced haplotype diversity in Eurasia compared to North America suggests that populations on the two continents have been affected by different historical processes.     

Genetic divergence and phylogeographic relationships among european perch (Perca fluviatilis) populations reflect glacial refugia and postglacial colonization

We used the widely distributed freshwater fish, perch (Perca fluviatilis), to investigate the postglacial colonization routes of freshwater fishes in Europe. Genetic variability within and among drainages was assessed using mitochondrial DNA (mtDNA) D-loop sequencing and RAPD markers from 55 populations all over Europe as well as one Siberian population. High level of structuring for both markers was observed among drainages and regions, while little differentiation was seen within drainages and regions. Phylogeographic relationships among European perch were determined from the distribution of 35 mtDNA haplotypes detected in the samples. In addition to a distinct southern European group, which includes a Greek and a southern Danubian population, three major groups of perch are observed: the western European drainages, the eastern European drainages including the Siberian population, and Norwegian populations from northern Norway, and western side of Oslofjord. Our data suggest that present perch populations in western and northern Europe were colonized from three main refugia, located in southeastern, northeastern and western Europe. In support of this, nested cladistic analysis of mtDNA clade and nested clade distances suggested historical range expansion as the main factor determining geographical distribution of haplotypes. The Baltic Sea has been colonized from all three refugia, and northeastern Europe harbours descendants from both eastern European refugia. In the upper part of the Danube lineages from the western European and the southern European refugia meet. The southern European refugium probably did not contribute to the recolonization of other western and northern European drainages after the last glaciation. However, phylogenetic analyses suggest that the southern European mtDNA lineage is the most ancient, and therefore likely to be the founder of all present perch lineages. The colonization routes used by perch probably also apply to other freshwater species with similar distribution patterns.     

Phylogeography and Pleistocene refugia of the Little Owl Athene noctua inferred from mtDNA sequence data

Pleistocene glaciations greatly affected the distribution of genetic diversity in animal populations. The Little Owl is widely distributed in temperate regions and could have survived the last glaciations in southern refugia. To describe the phylogeographical structure of European populations, we sequenced the mitochondrial cytochrome c oxidase I (COI) and control region (CR1) in 326 individuals sampled from 22 locations. Phylogenetic analyses of COI identified two deeply divergent clades: a western haplogroup distributed in western and northwestern Europe, and an eastern haplogroup distributed in southeastern Europe. Faster evolving CR1 sequences supported the divergence between these two main clades, and identified three subgroups within the eastern clade: Balkan, southern Italian and Sardinian. Divergence times estimated from COI with fossil calibrations indicate that the western and eastern haplogroups split 2.01–1.71 Mya. Slightly different times for splits were found using the standard 2% rate and 7.3% mtDNA neutral substitution rate. CR1 sequences dated the origin of endemic Sardinian haplotypes at 1.04–0.26 Mya and the split between southern Italian and Balkan haplogroups at 0.72–0.21 Mya, coincident with the onset of two Pleistocene glaciations. Admixture of mtDNA haplotypes was detected in northern Italy and in central Europe. These findings support a model of southern Mediterranean and Balkan refugia, with postglacial expansion and secondary contacts for Little Owl populations. Central and northern Europe was predominantly recolonized by Little Owls from Iberia, whereas expansion out of the Balkans was more limited. Northward expansion of the Italian haplogroup was probably prevented by the Alps, and the Sardinian haplotypes remained confined to the island. Results showed a clear genetic pattern differentiating putative subspecies. Genetic distances between haplogroups were comparable with those recorded between different avian species.     

Evolutionary history of a widespread tree species Acer mono in East Asia

East Asia has the most diverse temperate flora in the world primarily due to the lack of Pleistocene glaciation and the geographic heterogeneity. Although increasing phylogeography studies in this region provided more proofs in this issue, discrepancies and uncertainty still exist, especially in northern temperate deciduous broad‐leaved and coniferous mixed forest region (II). And a widespread plant species could reduce the complexity to infer the relationship between diversity and physiographical pattern. Hence, we studied the evolution history of a widespread temperate tree, Acer mono, populations in region II and the influence of physiographic patterns on intraspecific genetic diversity. Analyses of chloroplast sequences and nuclear microsatellites indicated high levels of genetic diversity. The diversity distribution was spatially heterogeneous and a latitudinal cline existed in both markers. The spatial distribution pattern between genetic diversity within A. mono and the diversity at species level was generally consistent. Western subtropical evergreen broad‐leaved forest subregion (IVb) had a unique ancient chloroplast clade (CP3) and a nuclear gene pool (GP5) with dominance indicating the critical role of this area in species diversification. Genetic data and ecological niche model results both suggested that populations in region II disappeared during the last glacial maximum (LGM) and recovered from south of Changbai Mt. and the Korean Peninsula. Two distribution centers were likely during the LGM, one in the north edge of warm temperate deciduous broad‐leaved forest region (III) and another in the south edge of region III. This was reflected by the genetic pattern with two spatially independent genetic groups. This study highlights the key role of region III in sustaining genetic diversity in the northern range and connecting diversity between southern and northern range. We elucidated the diversity relationship between vegetation regions which could facilitate the understanding of biodiversity origin and maintenance in East Asia.     

The role of disjunction and postglacial population expansion on phylogeographical history and genetic diversity of the circumboreal plant Chamaedaphne calyculata

In the last decade a number of studies has illustrated quite different phylogeographical patterns amongst plants with a northern present‐day geographical distribution, spanning the entire circumboreal region and/or circumarctic region and southern mountains. These works, employing several marker systems, have brought to light the complex evolutionary histories of this group. Here I focus on one circumboreal plant species, Chamaedaphne calyculata (leatherleaf), to unravel its phylogeographical history and patterns of genetic diversity across its geographical range. A survey of 29 populations with combined analyses of chloroplast DNA (cpDNA), internal transcribed spacer (ITS) and AFLP markers revealed structuring into two groups: Eurasian/north‐western North American, and north‐eastern North American. The present geographical distribution of C. calyculata has resulted from colonization from two putative refugial areas: east Beringia and south‐eastern North America. The variation of chloroplast DNA (cpDNA) and ITS sequences strongly indicated that the evolutionary histories of the Eurasian/north‐western North American and the north‐eastern North American populations were independent of each other because of a geographical disjunction in the distribution area and ice‐sheet history between north‐eastern and north‐western North America. Mismatch analysis using ITS confirmed that the present‐day population structure is the result of rapid expansion, probably since the last glacial maximum. The AFLP data revealed low genetic diversity of C. calyculata (P = 19.5%, H = 0.085) over the whole geographical range, and there was no evidence of loss of genetic diversity within populations in the continuous range, either at the margins or in formerly glaciated and nonglaciated regions. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 761–775.     

Phylogeography of Middle American gophersnakes: mixed responses to biogeographical barriers across the Mexican Transition Zone

Aim We used inferences of phylogeographical structure and estimates of divergence times for three species of gophersnakes (Colubridae: Pituophis) distributed across the Mexican Transition Zone (MTZ) to evaluate the postulated association of three Neogene geological events (marine seaway inundation of the Isthmus of Tehuantepec, formation of the Transvolcanic Belt across central Mexico, and secondary uplifting of the Sierra Madre Occidental) and of Pleistocene climate change with inter‐ and intraspecific diversification. Location Mexico, Guatemala, and the western United States. Methods We combined range‐wide sampling (67 individuals representing three putative species distributed across northern Middle America and western North America) and phylogenetic analyses of 1637 base pairs of mitochondrial DNA to estimate genealogical relationships and divergence times. The hypothesized concordance of inferred gene trees with geological histories was assessed using topology tests. Results We identified three major lineages of Middle American gophersnakes, and strong phylogeographical structure within each lineage. Gene trees were statistically congruent with hypothesized geological histories for two of the three postulated geological events. Estimated divergence dates and the geographical distribution of genetic variation further support mixed responses to these geological events. Considerable phylogeographical structure appears to have been generated during the Pleistocene. Main conclusions Phylogenetic and phylogeographical structure in gophersnakes distributed across northern Middle America and western North America highlights the influence of both Neogene vicariance events and Pleistocene climate change in shaping genetic diversity in this region. Despite the presence of two major geographical barriers in southern Mexico, extreme geological and environmental heterogeneity in this area may have differentially structured genetic diversity in highland taxa. To the north, co‐distributed taxa may display a more predictable pattern of diversification across the warm desert regions. Future studies should incorporate nuclear data to disentangle inferred lineage boundaries and further elucidate patterns of mitochondrial introgression.     

Comparative phylogeography of the Veronica alpina complex in Europe and North America

The Veronica alpina complex comprises eight species of alpine habitats over a wide range of mountain systems in the Northern Hemisphere. The occurrence of sympatric species in the European and North American mountain systems allowed us not only to investigate the effect of the ice ages on intraspecific phylogeographical patterns and genetic diversity in different continents of the Northern Hemisphere, but also to compare these patterns in closely related species. Plastid DNA trnL-F sequences and AFLP (amplified fragment length polymorphism) fingerprints were used to infer the phylogenetic history of the group and phylogeographical patterns within species. Hybrid origin of tetraploid eastern North American V. wormskjoldii from western North American V. nutans (= V. wormskjoldii s.l.) and Eurasian V. alpina is suggested. A number of phylogeographical groups have been found both in V. alpina from Europe and in V. nutans from western North America. Phylogeographical substructuring in the Alps is inferred for V. alpina but not for V. bellidioides, which is moreover characterized by an overall very low genetic diversity. Western North American V. cusickii is much more genetically diverse than its sympatric relative, V. nutans, an effect that is likely due to differences in the breeding system. Populations of V. nutans are differentiated into three groups, those from the Cascades and from the southern and the northern Rocky Mountains. Genetic diversity seems to be higher in the North American V. nutans than in the morphologically and ecologically similar European V. alpina. A possible scenario to explain this pattern is suggested.     

Phylogeography and population history of the crab-eating fox (Cerdocyon thous)

The crab-eating fox is a medium-sized Neotropical canid with generalist habits and a broad distribution in South America. We have investigated its genetic diversity, population structure and demographic history across most of its geographic range by analysing 512 base pairs (bp) of the mitochondrial DNA (mtDNA) control region, 615 bp of the mtDNA cytochrome b gene and 1573 total nucleotides from three different nuclear fragments. MtDNA data revealed a strong phylogeographic partition between northeastern Brazil and other portions of the species' distribution, with complete separation between southern and northern components of the Atlantic Forest. We estimated that the two groups diverged from each other c. 400,000-600,000 years ago, and have had contrasting population histories. A recent demographic expansion was inferred for the southern group, while northern populations seem to have had a longer history of large population size. Nuclear sequence data did not support this north-south pattern of subdivision, likely due at least in part to secondary male-mediated historical gene flow, inferred from multilocus coalescent-based analyses. We have compared the inferred phylogeographic patterns to those observed for other Neotropical vertebrates, and report evidence for a major north-south demographic discontinuity that seems to have marked the history of the Atlantic Forest biota.     

Phylogeographic structure, gene flow and species status in blue grouse (Dendragapus obscurus)

We investigated the genetic population structure and species status of a relatively sedentary bird that is a permanent resident of western North American forests, the blue grouse (Dendragapus obscurus). Phylogenetic analysis of complete mitochondrial control region DNA sequences resulted in the identification of three basal clades of haplotypes that were largely congruent with well-known biogeographical regions. These clades corresponded to the parapatric sooty (D. o. fuliginosus) and dusky (D. o. obscurus) subspecies groups of blue grouse plus a previously unrecognized division between northern and southern dusky grouse populations; the latter does not correspond closely to any currently recognized subspecies boundary. Approximately 66% of the total genetic variation was distributed among these three regions. Maximum likelihood estimates of gene flow between the regions were low or asymmetric; gene flow has been insufficient to prevent genetic divergence between dusky and sooty grouse. Estimates of gene flow among populations within sooty grouse were large except across the Columbia River valley. Among populations of dusky grouse, estimates of gene flow were heterogeneous and asymmetrical, reflecting large-scale fragmentation of the distribution due to landscape features and associated vegetation. Genetic, morphological and behavioural evidence suggest that sooty and dusky grouse are species-level taxa; the specific status of a third clade remains ambiguous.     

Postglacial range expansion from northern refugia by the wood frog, Rana sylvatica

Although the range dynamics of North American amphibians during the last glacial cycle are increasingly better understood, the recolonization history of the most northern regions and the impact of southern refugia on patterns of intraspecific genetic diversity and phenotypic variation in these regions are not well reconstructed. Here we present the phylogeographic history of a widespread and primarily northern frog, Rana sylvatica . We surveyed 551 individuals from 116 localities across the species' range for a 650-bp region of the NADH dehydrogenase subunit 2 and tRNA^TRP mitochondrial genes. Our phylogenetic analyses revealed two distinct clades corresponding to eastern and western populations, as well as a Maritime subclade within the eastern lineage. Patterns of genetic diversity support multiple refugia. However, high-latitude refugia in the Appalachian highlands and modern-day Wisconsin appear to have had the biggest impact on northern populations. Clustering analyses based on morphology further support a distinction between eastern and western wood frogs and suggest that postglacial migration has played an important role in generating broad-scale patterns of phenotypic variation in this species.     

Glacial refugia in a maritime temperate climate: Cicada (Kikihia subalpina) mtDNA phylogeography in New Zealand

Understanding the biological significance of Pleistocene glaciations requires knowledge of the nature and extent of habitat refugia during glacial maxima. An opportunity to examine evidence of glacial forest refugia in a maritime, Southern Hemisphere setting is found in New Zealand, where the extent of Pleistocene forests remains controversial. We used the mitochondrial phylogeography of a forest-edge cicada ( Kikihia subalpina ) to test the hypothesis that populations of this species survived throughout South Island during the Last Glacial Maximum. We also compared mitochondrial DNA phylogeographic patterns with male song patterns that suggest allopatric divergence across Cook Strait. Cytochrome oxidase I and II sequences were analyzed using network analysis, maximum-likelihood phylogenetic estimation, Bayesian dating and Bayesian skyline plots. K. subalpina haplotypes from North Island and South Island form monophyletic clades that are concordant with song patterns. Song divergence corresponds to approximately 2% genetic divergence, and Bayesian dating suggests that the North Island and South Island population-lineages became isolated around 761 000 years bp . Almost all South Island genetic variation is found in the north of the island, consistent with refugia in Marlborough Sounds, central Nelson and northwest Nelson. All central and southern South Island and Stewart Island haplotypes are extremely similar to northern South Island haplotypes, a 'northern richness/southern purity' pattern that mirrors genetic patterns observed in many Northern Hemisphere taxa. Proposed southern South Island forest habitat fragments may have been too small to sustain populations of K. subalpina , and/or they may have harboured ecological communities with no modern-day analogues.