Molecular signatures of divergence and selection in closely related pine taxa

Efforts to detect loci under selection in plants have mostly focussed on single species. However, assuming that intraspecific divergence may lead to speciation, comparisons of genetic variation within and among recently diverged taxa can help to locate such genes. In this study, coalescent and outlier detection methods were used to assess nucleotide polymorphism and divergence at 79 nuclear gene fragments (1212 SNPs) in 16 populations (153 individuals) of the closely related, but phenotypically and ecologically distinct, pine taxa Pinus mugo, P. uliginosa and P. uncinata across their European distributions. Simultaneously, mitochondrial DNA markers, which are maternally inherited in pines and distributed by seeds at short geographic distance, were used to assess genetic relationships of the focal populations and taxa. The majority of nuclear loci showed homogenous patterns of variation between the taxa due to a high number of shared SNPs and haplotypes, similar levels of polymorphism, and low net divergence. However, against this common genetic background and an overall low population structure within taxa at mitochondrial markers, we identified several genes showing signatures of selection, accompanied by significant intra- and interspecific divergence. Our results indicate that loci involved in species divergence may be involved in intraspecific local adaptation.     

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.     

Phylogeography of a stream-dwelling frog (Pseudacris cadaverina) in southern California

Recent phylogeographic studies of animal taxa in California have revealed common geographic patterns of evolutionary divergence and genetic diversity that are generally attributable to landscape influences. However, there remains a paucity of knowledge on the evolution of freshwater taxa in southern California. Here, we investigate phylogeographic patterns in a stream-dwelling frog (Pseudacris cadaverina). Two hundred and twenty-one individuals were collected from 46 populations across the species’ range in southern California. Using 1100 bp of sequence data from cytochrome b and tRNA-Glu, we conducted phylogenetic analyses, analysis of molecular variance, and nested clade phylogeographic analysis to gain insight into the factors contributing to the distribution of genetic diversity in P. cadaverina. We tested for evidence of two putative phylogeographic breaks and tested hypotheses that genetic diversity in this species is partitioned into (1) major watersheds, (2) mountain ranges, and (3) coastal and desert regions. Our results suggest that the eastern Transverse Ranges are the center of origin for extant P. cadaverina lineages and that the observed genetic structure in this species was established during the Pleistocene Epoch. There is strong support for three major haplotype groups and a Transverse Range break in P. cadaverina that is concordant with breaks found in numerous other taxa. The distribution of genetic diversity in P. cadaverina is due in large part to the separation of populations into different major watersheds and mountain ranges. Gene flow appears to be generally limited among disjunct populations throughout the region and some desert populations have been isolated by historical habitat fragmentation.     

Hybridization among Arctic white-headed gulls (Larus spp.) obscures the genetic legacy of the Pleistocene

We studied the influence of glacial oscillations on the genetic structure of seven species of white-headed gull that breed at high latitudes (Larus argentatus, L. canus, L. glaucescens, L. glaucoides, L. hyperboreus, L. schistisagus, and L. thayeri). We evaluated localities hypothesized as ice-free areas or glacial refugia in other Arctic vertebrates using molecular data from 11 microsatellite loci, mitochondrial DNA (mtDNA) control region, and six nuclear introns for 32 populations across the Holarctic. Moderate levels of genetic structure were observed for microsatellites (F(ST)= 0.129), introns (Φ(ST)= 0.185), and mtDNA control region (Φ(ST)= 0.461), with among-group variation maximized when populations were grouped based on subspecific classification. Two haplotype and at least two allele groups were observed across all loci. However, no haplotype/allele group was composed solely of individuals of a single species, a pattern consistent with recent divergence. Furthermore, northernmost populations were not well differentiated and among-group variation was maximized when L. argentatus and L. hyberboreus populations were grouped by locality rather than species, indicating recent hybridization. Four populations are located in putative Pleistocene glacial refugia and had larger τ estimates than the other 28 populations. However, we were unable to substantiate these putative refugia using coalescent theory, as all populations had genetic signatures of stability based on mtDNA. The extent of haplotype and allele sharing among Arctic white-headed gull species is noteworthy. Studies of other Arctic taxa have generally revealed species-specific clusters as well as genetic structure within species, usually correlated with geography. Aspects of white-headed gull behavioral biology, such as colonization ability and propensity to hybridize, as well as their recent evolutionary history, have likely played a large role in the limited genetic structure observed.     

Imprints from genetic drift and mutation imply relative divergence times across marine transition zones in a pan-European small pelagic fish (Sprattus sprattus)

Geographic distributions of most temperate marine fishes are affected by postglacial recolonisation events, which have left complex genetic imprints on populations of marine species. This study investigated population structure and demographic history of European sprat (Sprattus sprattus L.) by combining inference from both mtDNA and microsatellite genetic markers throughout the species' distribution. We compared effects from genetic drift and mutation for both genetic markers in shaping genetic differentiation across four transition zones. Microsatellite markers revealed significant isolation by distance and a complex population structure across the species' distribution (overall θ(ST)=0.038, P<0.01). Across transition zones markers indicated larger effects of genetic drift over mutations in the northern distribution of sprat contrasting a stronger relative impact of mutation in the species' southern distribution in the Mediterranean region. These results were interpreted to reflect more recent divergence times between northern populations in accordance with previous findings. This study demonstrates the usefulness of comparing inference from different markers and estimators of divergence for phylogeographic and population genetic studies in species with weak genetic structure, as is the case in many marine species.     

Vicariance and dispersal across an intermittent barrier: population genetic structure of marine animals across the Torres Strait land bridge

Biogeographic barriers, some transitory in duration, are likely to have been important contributing factors to modern marine biodiversity in the Indo-Pacific region. One such barrier was the Torres Strait land bridge between continental Australia and New Guinea that persisted through much of the late Pleistocene and separated Indian and Pacific Ocean taxa. Here, we examine the patterns of mitochondrial DNA diversity for marine animals with present-day distributions spanning the Torres Strait. Specifically, we investigate whether there are concordant signatures across species, consistent with either vicariance or recent colonization from either ocean basin. We survey four species of reef fishes (Apogon doederleini, Pomacentrus coelestis, Dascyllus trimaculatus, and Acanthurus triostegus) for mtDNA cytochrome oxidase 1 and control region variation and contrast these results to previous mtDNA studies in diverse marine animals with similar distributions. We find substantial genetic partitioning (estimated from F-statistics and coalescent approaches) between Indian and Pacific Ocean populations for many species, consistent with regional persistence through the late Pleistocene in both ocean basins. The species-specific estimates of genetic divergence, however, vary greatly and for reef fishes we estimate substantially different divergence times among species. It is likely that Indian and Pacific Ocean populations have been isolated for multiple glacial cycles for some species, whereas for other species genetic connections have been more recent. Regional estimates of genetic diversity and directionality of gene flow also vary among species. Thus, there is no apparent consistency among historical patterns across the Torres Strait for these co-distributed marine animals.     

The distribution of Quercus suber chloroplast haplotypes matches the palaeogeographical history of the western Mediterranean

Combining molecular analyses with geological and palaeontological data may reveal timing and modes for the divergence of lineages within species. The Mediterranean Basin is particularly appropriate for this kind of multidisciplinary studies, because of its complex geological history and biological diversity. Here, we investigated chloroplast DNA of Quercus suber populations in order to detect possible relationships between their geographical distribution and the palaeogeographical history of the western Mediterranean domain. We analysed 110 cork oak populations, covering the whole distribution range of the species, by 14 chloroplast microsatellite markers, among which eight displayed variation among populations. We identified five haplotypes whose distribution is clearly geographically structured. Results demonstrated that cork oak populations have undergone a genetic drift geographically consistent with the Oligocene and Miocene break-up events of the European-Iberian continental margin and suggested that they have persisted in a number of separate microplates, currently found in Tunisia, Sardinia, Corsica, and Provence, without detectable chloroplast DNA modifications for a time span of over 15 million years. A similar distribution pattern of mitochondrial DNA of Pinus pinaster supports the hypothesis of such long-term persistence, in spite of Quaternary climate oscillations and of isolation due to insularity, and suggests that part of the modern geographical structure of Mediterranean populations may be traced back to the Tertiary history of taxa.     

Evidence for stabilizing selection and slow divergent evolution of male genitalia in a millipede (Antichiropus variabilis)

It is generally accepted that postcopulatory sexual selection drives rapid divergence of genital morphology among isolated populations. The mode of selection operating upon genitalia can be explored by comparing patterns of population divergence in genetic and genitalic traits. We collected Antichiropus variabilis millipedes from eight localities across the species range. Levels of among-population genetic divergence, at microsatellite loci, and the mitochondrial COI gene were very high. Following geometric morphometric analyses, genital morphology was also found to be highly divergent among the populations surveyed, whereas head morphology had not diverged as markedly. However, pairwise comparisons of F(ST) and P(ST) showed that among-population divergence in both genital and head shape was significantly lower than that experienced by neutral genetic markers. Our results suggest that the genitalia of A. variabilis are currently experiencing a period of stabilizing selection, the mode of selection expected for genitalia that function in species recognition via a "lock-and-key" mechanism. Our results demonstrate that although genital morphology can clearly diverge among genetically isolated populations, divergence is not necessarily as rapid as commonly argued, and continuous directional sexual selection may not always underpin the evolutionary divergence of male genitalia.     

Testing Pleistocene refugia theory: phylogeographical analysis of Desmognathus wrighti,a high-elevation salamander in the southern Appalachians

During the colder climates of the Pleistocene, the ranges of high-elevation species in unglaciated areas may have expanded, leading to increased gene flow among previously isolated populations. The phylogeography of the pygmy salamander, Desmognathus wrighti, an endemic species restricted to the highest mountain peaks of the southern Appalachians, was examined to test the hypothesis that the range of D. wrighti expanded along with other codistributed taxa during the Pleistocene. Analyses of genetic variation at 14 allozymic loci and of the 12S rRNA gene in the mtDNA genome was conducted on individuals sampled from 14 population isolates throughout the range of D. wrighti. In contrast to the genetic patterns of many other high-elevation animals and plants, genetic distances derived from both molecular markers showed significant isolation by distance and genetic structuring of populations, suggesting long-term isolation of populations. Phylogeographical analyses revealed four genetically distinct population clusters that probably remained fragmented during the Pleistocene, although there was also evidence supporting recent gene flow among some population groups. Support for isolation by distance is rare among high-elevation species in unglaciated areas of North and Middle America, although not uncommon among Plethodontid Salamanders, and this pattern suggests that populations of D. wrighti did not expand entirely into suitable habitat during the Pleistocene. We propose that intrinsic barriers to dispersal, such as species interactions with other southern Appalachian plethodontid salamanders, persisted during the Pleistocene to maintain the fragmented distribution of D. wrighti and allow for significant genetic divergence of populations by restricting gene flow.     

Barriers,rather than refugia,underlie the origin of diversity in toads endemic to the Brazilian Atlantic Forest

In this study, we investigated the relative contribution of geographic barriers and Pleistocene refuges in the diversification of the Rhinella crucifer species complex, a group of endemic toads with a widespread distribution in the Brazilian Atlantic Forest (AF). We used intensive sampling and multilocus DNA sequence data to compare nucleotide diversity between refuge and nonrefuge areas, investigate regional demographic patterns, estimate demographic parameters related to genetic breaks and test refuge versus barrier scenarios of diversification using approximate Bayesian computation. We did not find higher levels of genetic diversity in putative refuge areas, either at regional or biome scale. Rather, the demographic history of the species complex supports regional differences with moderate population growth in the north and central regions and stability in southern AF. Genetic breaks were dated to the Plio–Pleistocene; however, our analyses rejected the role of refuges in creating a northern and central divergence, supporting a recent colonization scenario at a smaller scale within the central AF. Overall, our data rule out massive climatically driven fragmentation and large‐scale recolonization events for populations across the biome. We confirmed the importance of geographic barriers in creating main divergences and underscored the importance of searching for cryptic discontinuities in the landscape. Comparison of our results with those of other AF taxa indicates organismal specific responses to moderate shifts in habitat and that multiple refuges may constitute a more realistic model for diversification of Atlantic Forest biota.     

Glacial oceanographic contrasts explain phylogeography of Australian bull kelp

The evolutionary effects of Southern Hemisphere Pleistocene oceanographic conditions — marked by fluctuations in sea levels and water temperatures, and redirected currents — are poorly understood. The southeastern tip of Australia presents an intriguing model system for studying the biological impacts of palaeoceanography. In particular, contrasting oceanographic conditions that existed on eastern vs. western sides of the Bassian Isthmus during Pleistocene glacial periods allow for natural comparisons between putative refugial vs. re-invading populations. Whereas many western Tasmanian marine taxa were likely eliminated by cold subantarctic water during the last glacial period, eastern Tasmanian populations would have persisted in relatively warm temperatures mediated by the ongoing influence of the East Australian Current (EAC). Here we test for the effects of contrasting palaeoceanographic conditions on endemic bull kelp, Durvillaea potatorum , using DNA sequence analysis (COI; rbc L) of more than 100 individuals from 14 localities in southeastern Australia. Phylogenetic reconstructions reveal a deep (maximum divergence 4.7%) genetic split within D. potatorum , corresponding to the 'eastern' and 'western' geographical regions delimited by the Bassian Isthmus, a vicariant barrier during low Pleistocene sea levels. Concordant with the western region's cold glacial conditions, samples from western Tasmania and western Victoria are genetically monomorphic, suggesting postglacial expansion from a mainland refugium. Eastern samples, in contrast, comprise distinct regional haplogroups, suggesting the species persisted in eastern Tasmania throughout recent glacial periods. The deep east–west divergence seems consistent with earlier reports of morphological differences between 'western' and 'eastern' D. potatorum , and it seems likely that these forms represent reproductively isolated species.     

Comparative phylogeography of pitvipers suggests a consensus of ancient Middle American highland biogeography

Aim  We used inferences of phylogenetic relationships and divergence times for three lineages of highland pitvipers to identify broad-scale historical events that have shaped the evolutionary history of Middle American highland taxa, and to test previous hypotheses of Neotropical speciation.
Location  Middle America (Central America and Mexico).
Methods  We used 2306 base pairs of mitochondrial gene sequences from 178 individuals to estimate the phylogeny and divergence times of New World pitviper lineages, focusing on three genera ( Atropoides , Bothriechis and Cerrophidion ) that are broadly co-distributed across Middle American highlands.
Results  We found strong correspondence across three highland lineages for temporally and geographically coincident divergences in the Miocene and Pliocene, and further identified widespread within-species divergences across multiple lineages that occurred in the early–middle Pleistocene.
Main conclusions  Available data suggest that there were at least three major historical events in Middle America that had broad impacts on species divergence and lineage diversification among highland taxa. In addition, we find widespread within-species genetic structure that may be attributable to the climatic changes that affected gene flow among highland taxa during the middle–late Pleistocene.     

Comparative phylogeography of three codistributed stomatopods: origins and timing of regional lineage diversification in the Coral Triangle

The Indonesian-Australian Archipelago is the center of the world's marine biodiversity. Although many biogeographers have suggested that this region is a "center of origin," criticism of this theory has focused on the absence of processes promoting lineage diversification in the center. In this study we compare patterns of phylogeographic structure and gene flow in three codistributed, ecologically similar Indo-West Pacific stomatopod (mantis shrimp) species. All three taxa show evidence for limited gene flow across the Maluku Sea with deep genetic breaks between populations from Papua and Northern Indonesia, suggesting that limited water transport across the Maluku Sea may limit larval dispersal and gene flow across this region. All three taxa also show moderate to strong genetic structure between populations from Northern and Southern Indonesia, indicating limited gene flow across the Flores and Java Seas. Despite the similarities in phylogeographic structure, results indicate varied ages of the genetic discontinuities, ranging from the middle Pleistocene to the Pliocene. Concordance of genetic structure across multiple taxa combined with temporal discordance suggests that regional genetic structures have arisen from the action of common physical processes operating over extended time periods. The presence in all three species of both intraspecific genetic structure as well as deeply divergent lineages that likely represent cryptic species suggests that these processes may promote lineage diversification within the Indonesian-Australian Archipelago, providing a potential mechanism for the center of origin. Efforts to conserve biodiversity in the Coral Triangle should work to preserve both existing biodiversity as well as the processes creating the biodiversity.     

Is Gene Flow Promoting the Reversal of Pleistocene Divergence in the Mountain Chickadee (Poecile gambeli)?

The Pleistocene glacial cycles left a genetic legacy on taxa throughout the world; however, the persistence of genetic lineages that diverged during these cycles is dependent upon levels of gene flow and introgression. The consequences of secondary contact among taxa may reveal new insights into the history of the Pleistocene’s genetic legacy. Here, we use phylogeographic methods, using 20 nuclear loci from regional populations, to infer the consequences of secondary contact following divergence in the Mountain Chickadee (Poecile gambeli). Analysis of nuclear data identified two geographically-structured genetic groups, largely concordant with results from a previous mitochondrial DNA (mtDNA) study. Additionally, the estimated multilocus divergence times indicate a Pleistocene divergence, and are highly concordant with mtDNA. The previous mtDNA study showed a paucity of sympatry between clades, while nuclear patterns of gene flow show highly varied patterns between populations. The observed pattern of gene flow, from coalescent-based analyses, indicates southern populations in both clades exhibit little gene flow within or between clades, while northern populations are experiencing higher gene flow within and between clades. If this pattern were to persist, it is possible the historical legacy of Pleistocene divergence may be preserved in the southern populations only, and the northern populations would become a genetically diverse hybrid species.     

Phylogeography and species biogeography of montane Great Basin stoneflies

Sky islands are ideal systems for determining the effects of climatic oscillations on species distributions and genetic structure. Our study focused on montane stonefly populations in the Great Basin of western North America. We used niche-based distribution modelling, phylogeography and traditional species-based biogeography to test several hypotheses as follows: (i) genetic differentiation among Doroneuria baumanni populations will be independent of hydrologic connectivity (headwater model); (ii) Sky islands were colonized when habitat was more continuous and populations likely experienced multiple expansions and contractions; (iii) Colonization events were coincident with the late Pleistocene and Holocene; and (iv) Shared topography and climate history will result in concordant patterns of genetic differentiation in D. baumanni and occurrences of 32 stonefly species across the region. Overall, Φ(ST) 's and coalescent-based estimates of migration were consistent with the headwater model. Maximum likelihood and Bayesian gene trees identified three major nonoverlapping east-west clades. Distribution modelling indicated more suitable habitat in the Great Basin during the Last Glacial Maximum than at present, but none during the last interglacial period. Demographic analyses showed evidence of population expansion in one of the three major east-west clades. Intra-clade divergence times (60,000-183,000ybp) were well within the late Pleistocene while among-clade divergence times (499.000-719,000ybp) were deeper. Genetic differentiation in D. baumanni and distributions of stonefly species were significantly concordant. These results imply that climatic oscillations have played major roles in shaping the genetic structure and distributions of Great Basin stoneflies, but that divergence among clades occurred much earlier than our late Pleistocence/early Holocene predictions.     

Long‐term panmixia in a cosmopolitan Indo‐Pacific coral reef fish and a nebulous genetic boundary with its broadly sympatric sister species

Phylogeographical studies have shown that some shallow‐water marine organisms, such as certain coral reef fishes, lack spatial population structure at oceanic scales, despite vast distances of pelagic habitat between reefs and other dispersal barriers. However, whether these dispersive widespread taxa constitute long‐term panmictic populations across their species ranges remains unknown. Conventional phylogeographical inferences frequently fail to distinguish between long‐term panmixia and metapopulations connected by gene flow. Moreover, marine organisms have notoriously large effective population sizes that confound population structure detection. Therefore, at what spatial scale marine populations experience independent evolutionary trajectories and ultimately species divergence is still unclear. Here, we present a phylogeographical study of a cosmopolitan Indo‐Pacific coral reef fish Naso hexacanthus and its sister species Naso caesius, using two mtDNA and two nDNA markers. The purpose of this study was two‐fold: first, to test for broad‐scale panmixia in N. hexacanthus by fitting the data to various phylogeographical models within a Bayesian statistical framework, and second, to explore patterns of genetic divergence between the two broadly sympatric species. We report that N. hexacanthus shows little population structure across the Indo‐Pacific and a range‐wide, long‐term panmictic population model best fit the data. Hence, this species presently comprises a single evolutionary unit across much of the tropical Indian and Pacific Oceans. Naso hexacanthus and N. caesius were not reciprocally monophyletic in the mtDNA markers but showed varying degrees of population level divergence in the two nuclear introns. Overall, patterns are consistent with secondary introgression following a period of isolation, which may be attributed to oceanographic conditions of the mid to late Pleistocene, when these two species appear to have diverged.     

The Pleistocene history of the sheepshead minnow (Cyprinodon variegatus): Non-equilibrium evolutionary dynamics within a diversifying species complex

The sheepshead minnow, Cyprinodon variegatus, is a widespread fish species that typically inhabits coastal tidal marsh and mangrove swamp environments, ranging from Cape Cod, Massaschusetts to northern Mexico and into the Caribbean. This wide range crosses several biogeographic boundaries which are coincident with genetic structuring within numerous species originating in the Pleistocene. In addition, the more northerly reaches of this species range have been further subject to the evolutionary consequences of Pleistocene glaciation due to local extinction and recolonization of formerly glaciated sites. C. variegatus thus provides an excellent vertebrate model system within which to test the extent of genetic differentiation among populations in a dominant coastal ecosystem and examine patterns of historical demography in populations distributed along a latitudinal gradient. Using mitochondrial control region and ND2 sequence data, we discovered monophyletic clades within C. variegatus with divergence times within the Pleistocene, and very low gene flow between most sites. Intraspecific genetic breaks appear to correspond broadly to biogeographic or oceanic boundaries. Pleistocene climate change appears to have had dramatic impacts on the size and distribution of populations within and near the glacial margins, but has also affected populations far from formerly glaciated regions.     

Genetic relationships among marine and freshwater populations of the European three-spined stickleback (Gasterosteus aculeatus) revealed by microsatellites

To assess the population genetic structure of the three-spined stickleback, Gasterosteus aculeatus, variability at 18 microsatellite loci was examined in 1724 individuals from 74 locations covering most of the species distribution range in Europe. The results revealed high overall degree of differentiation (F(ST) = 0.21) but contrasting level of divergence and genetic variability between habitat types. Marine populations were genetically relatively uniform even across great geographical distances as compared to substantial differentiation among freshwater populations. Analysis of molecular variance indicated low but significant (2.7%) variation in allele frequencies between geographical regions, but a negligible effect of habitat type (0.2%). The phylogenetic pattern was not explained by habitat type, but a weak signal of populations clustering according to geographical or water system origin was found. The results support the view that three-spined stickleback marine ancestors colonized northern European fresh waters during the postglacial marine submergence c. 10,000 years ago, whereas in the Mediterranean region colonization probably dates back to the Pleistocene. The independent origins of river and lake populations indicate that they originate from multiple colonizations rather than sharing common ancestry. In the continuous marine environment, the low degree of differentiation among populations can be explained by gene flow among subpopulations and large effective population size buffering divergence in neutral markers. In contrast, among postglacially established freshwater populations differentiation appears to be driven by genetic drift and isolation. The stepwise mutations appear to have contributed to the population differentiation in the southern part of the three-spined stickleback distribution range.     

Evolution of Bombina bombina and Bombina variegata (Anura: Discoglossidae) in the Carpathian Basin: a history of repeated mt-DNA introgression across species

The structure and geographic location of hybrid zones change through time. Current patterns result from present and historical population-environment interactions that act on each of the hybridizing taxa. This is particularly evident for species involved in complex hybrid zones, such as that formed by the toad species Bombina bombina and Bombina variegata (Anura: Discoglossidae), which interact along extensive areas in Central Europe. We used data on external morphology and partial sequences of the cytochrome oxidase I (cox1) and nicotinamide adenine dinucleotid dehydrogenase subunit 4 (nad4) mitochondrial DNA (mt-DNA) genes to analyze the current patterns of genetic structure shown by both species of Bombina along their contact zone in Hungary. Phylogenetic, phylogeographic, and historical demography analyses were applied to 1.5kb mt-DNA obtained from 119 individuals representing 24 populations from Hungary and additional specimens from Slovakia, Albania, and Bosnia-Herzegovina. We use these data to infer the evolutionary history of the isolated populations of B. variegata in Hungary and to discriminate between competing biogeographic scenarios accounting for the historical interactions between species in this region. Results from the inferred phylogenetic branching pattern and sequence divergence among species and populations support the following: (i) recent population expansion has occurred in Hungarian populations of B. bombina, which are genetically very homogeneous; (ii) the Hungarian populations of B. variegata correspond to two distinct mitochondrial lineages (Carpathian and Alpine, respectively); average maximum-likelihood-corrected sequence divergence between these lineages is 8.96% for cox1 and 10.85% for nad4; (iii) mt-DNA divergence among the three isolated western populations of B. variegata from Transdanubia is low, with four closely related haplotypes, which suggests that the isolation between these populations is the result of a recent process, possibly mediated by the invasion of B. bombina; and (iv) we have detected discordances between morphology and mt-DNA data in the Transdanubia region (Bakony Mountains, Mecsek Mountains, Orség area), suggesting mt-DNA introgression across species in this regions. These results are discussed with reference to previous biogeographic hypotheses.