Nested clade analysis and the genetic evidence for population expansion in the phylogeography of the golden-striped salamander, Chioglossa lusitanica (Amphibia: Urodela)

Genetic traces of hybrid zone movement across a fragmented habitat

Theoretical and empirical studies suggest that the structure and position of hybrid zones can change over time. Evidence for moving hybrid zones has been directly inferred by repeated sampling over time, or indirectly through the detection of genetic footprints left by the receding species and the resulting asymmetric patterns of introgression across markers. We here investigate a hybrid zone formed by two subspecies of the Iberian golden-striped salamander, Chioglossa lusitanica, using a panel of 35 nuclear loci (31 SNPs and 4 allozymes) and one mitochondrial locus in a transect in central Portugal. We found concordant and coincident clines for most of the nuclear loci (n=22, 63%), defining a narrow hybrid zone of ca. 6 km wide, with the centre positioned ca. 15 km south of the Mondego river. Asymmetric introgression was observed at another 14 loci. Their clines are displaced towards the north, with positions located either close to the Mondego river (n=6), or further northwards (n=8). We interpret these profiles as genetic traces of the southward displacement of C. lusitanica lusitanica by C. l. longipes over the wider Mondego river valley. We noted the absence of significant linkage disequilibrium and we inferred low levels of effective selection per locus against hybrids, suggesting that introgression in the area of species replacement occurred under a neutral diffusion process. A species distribution model suggests that the C. lusitanica hybrid zone coincides with a narrow corridor of fragmented habitat. From the position of the displaced clines, we infer that patches of locally suitable habitat trapped some genetic variants that became disassociated from the southward moving hybrid zone. This study highlights the influence of habitat availability on hybrid zone movement.

The Asian plethodontid salamander preserves historical genetic imprints of recent northern expansion

The Korean Peninsula, located at the southern tip of Northeast Asia, has never been covered by ice sheets and was a temperate refugium during the Pleistocene. Karsenia koreana, the sole Asian plethodontid salamander species, occurs only on the southern half of the Korean Peninsula and is thought to have found various climatic refugia. Despite its phylogenetic and biogeographic importance, no population-level genetic analysis has been performed on this species. Here we study the population genetic structure of K. koreana using mitochondrial and microsatellite loci to understand the recent historical dispersion process that shaped its current distribution. Overall, the genetic distance between populations correlated well with the spatial distance, and the genetic structure among populations showed signs of a unilateral northward expansion from a southernmost refugium population. Given the distinct genetic structure formed among the populations, the level of historical gene flow among populations appears to have been very low. As the estimated effective population size of K. koreana was also small, these results suggest that the small, restricted populations of K. koreana are extremely vulnerable to environmental changes that may require high levels of genetic diversity to cope with. Thus, special management strategies are needed to preserve these remnant populations.

Is Malaysia's banded langur, Presbytis femoralis femoralis, actually Presbytis neglectus neglectus? Taxonomic revision with new insights on the radiation history of the Presbytis species group in Southeast Asia

The disjunct distribution of Presbytis femoralis subspecies across Sumatra (P. f. percura), southern (P. f. femoralis) and northern (P. f. robinsoni) Peninsular Malaysia marks the unique vicariance events in the Sunda Shelf. However, the taxonomic positions and evolutionary history of P. f. femoralis are unresolved after decades of research. To elucidate this evolutionary history, we analyzed 501 base pairs of the mitochondrial HVSI gene from 25 individuals representing Malaysia's banded langur, with the addition of 29 sequences of Asian Presbytis from Genbank. Our results revealed closer affinity of P. f. femoralis to P. m. mitrata and P. m. sumatrana while maintaining the monophyletic state of P. f. femoralis as compared to P. f. robinsoni. Two central theses were inferred from the results; (1) P. f. femoralis does not belong in the same species classification as P. f. robinsoni, and (2) P. f. femoralis is the basal lineage of the Presbytis in Peninsular Malaysia. Proving the first hypothesis through genetic analysis, we reassigned P. f. femoralis of Malaysia to Presbytis neglectus (Schlegel's banded langur) (Schlegel in Revue Methodique, Museum d'Histoire Naturelle des Pays-Bas 7:1, 1876) following the International Code of Zoological Nomenclature (article 23.3). The ancestors of P. neglectus are hypothesized to have reached southern Peninsular Malaysia during the Pleistocene and survived in refugium along the western coast. Consequently, they radiated upward, forming P. f. robinsoni and P. siamensis resulting in the highly allopatric distribution in Peninsular Malaysia. This study has successfully resolved the taxonomic position of P. neglectus in Peninsular Malaysia while providing an alternative biogeographic theory for the Asian Presbytis.

Mitochondrial phylogeography of the Iberian endemic frog Rana iberica, with implications for its conservation

Genetic characterization of species using phylogeographic approaches represents a basic reference to understand their evolutionary history as well as to identify conservation priorities to protect areas of particular interest regarding evolutionary potential. Even in well-studied regions such information is lacking for the majority of species, including many endemic species with reduced distribution ranges. We investigate the phylogeographic pattern of the Iberian frog Rana iberica, an endemic amphibian restricted to Central and North-Western Iberian Peninsula. Using mitochondrial sequences, we reconstruct the phylogeographic history of the species to test the effect of Quaternary climate changes on the evolutionary diversification of lineages, that is, the differentiation of mitochondrial lineages and the formation of genetic diversity melting pots, and integrate phylogeographic evidence for future conservation planning. Our results indicate the existence of 3 main mitochondrial lineages differentiated during the Upper Pleistocene. Both historical demographic analyses and climatic niche modeling show a strong effect of glacial climate changes, suggesting recurrent range contractions and expansions. Under such circumstances, differentiation took place most likely by isolation in allopatric interglacial refugia. Secondary lineage admixture in northern Portugal generated a broad mixed zone with highest nucleotide diversity. Given its particular evolutionary potential, its reduced distribution and eventual threats under current climate change scenario, conservation priorities should focus on the isolated lineage from Sierra de Guadalupe.

Three in One--Multiple Faunal Elements within an Endangered European Butterfly Species

Ice ages within Europe forced many species to retreat to refugia, of which three major biogeographic basic types can be distinguished: "Mediterranean", "Continental" and "Alpine / Arctic" species. However, this classification often fails to explain the complex phylogeography of European species with a wide range of latitudinal and altitudinal distribution. Hence, we tested for the possibility that all three mentioned faunal elements are represented within one species. Our data was obtained by scoring 1,307 Euphydryas aurinia individuals (46 European locations) for 17 allozyme loci, and sequencing a subset of 492 individuals (21 sites) for a 626 base pairs COI fragment. Genetic diversity indices, F statistics, hierarchical analyses of molecular variance, individual-based clustering, and networks were used to explore the phylogeographic patterns. The COI fragment represented 18 haplotypes showing a strong geographic structure. All but one allozyme loci analysed were polymorphic with a mean F_ST of 0.20, supporting a pronounced among population structure. Interpretation of both genetic marker systems, using several analytical tools, calls for the recognition of twelve genetic groups. These analyses consistently distinguished different groups in Iberia (2), Italy, Provence, Alps (3), Slovenia, Carpathian Basin, the lowlands of West and Central Europe as well as Estonia, often with considerable additional substructures. The genetic data strongly support the hypothesis that E. aurinia survived the last glaciation in Mediterranean, extra-Mediterranean and perialpine refugia. It is thus a rare example of a model organism that combines attributes of faunal elements from all three of these sources. The observed differences between allozymes and mtDNA most likely result from recent introgression of mtDNA into nuclear allozyme groups. Our results indicate discrepancies with the morphologically-based subspecies models, underlining the need to revise the current taxonomy.

Elevated Genetic Diversity in the Emerging Blueberry Pathogen Exobasidium maculosum

Emerging diseases caused by fungi are increasing at an alarming rate. Exobasidium leaf and fruit spot of blueberry, caused by the fungus Exobasidium maculosum, is an emerging disease that has rapidly increased in prevalence throughout the southeastern USA, severely reducing fruit quality in some plantings. The objectives of this study were to determine the genetic diversity of E. maculosum in the southeastern USA to elucidate the basis of disease emergence and to investigate if populations of E. maculosum are structured by geography, host species, or tissue type. We sequenced three conserved loci from 82 isolates collected from leaves and fruit of rabbiteye blueberry (Vaccinium virgatum), highbush blueberry (V. corymbosum), and southern highbush blueberry (V. corymbosum hybrids) from commercial fields in Georgia and North Carolina, USA, and 6 isolates from lowbush blueberry (V. angustifolium) from Maine, USA, and Nova Scotia, Canada. Populations of E. maculosum from the southeastern USA and from lowbush blueberry in Maine and Nova Scotia are distinct, but do not represent unique species. No difference in genetic structure was detected between different host tissues or among different host species within the southeastern USA; however, differentiation was detected between populations in Georgia and North Carolina. Overall, E. maculosum showed extreme genetic diversity within the conserved loci with 286 segregating sites among the 1,775 sequenced nucleotides and each isolate representing a unique multilocus haplotype. However, 94% of the nucleotide substitutions were silent, so despite the high number of mutations, selective constraints have limited changes to the amino acid sequences of the housekeeping genes. Overall, these results suggest that the emergence of Exobasidium leaf and fruit spot is not due to a recent introduction or host shift, or the recent evolution of aggressive genotypes of E. maculosum, but more likely as a result of an increasing host population or an environmental change.

Continental monophyly and molecular divergence of Peninsular Malaysia's Macaca fascicularis fascicularis

The phylogenetic relationships of long-tailed macaque (Macaca fascicularis fascicularis) populations distributed in Peninsular Malaysia in relation to other regions remain unknown. The aim of this study was to reveal the phylogeography and population genetics of Peninsular Malaysia's M. f. fascicularis based on the D-loop region of mitochondrial DNA. Sixty-five haplotypes were detected in all populations, with only Vietnam and Cambodia sharing four haplotypes. The minimum-spanning network projected a distant relationship between Peninsular Malaysian and insular populations. Genetic differentiation (F_ST, Nst) results suggested that the gene flow among Peninsular Malaysian and the other populations is very low. Phylogenetic tree reconstructions indicated a monophyletic clade of Malaysia's population with continental populations (NJ = 97%, MP = 76%, and Bayesian = 1.00 posterior probabilities). The results demonstrate that Peninsular Malaysia's M. f. fascicularis belonged to Indochinese populations as opposed to the previously claimed Sundaic populations. M. f. fascicularis groups are estimated to have colonized Peninsular Malaysia ~0.47 million years ago (MYA) directly from Indochina through seaways, by means of natural sea rafting, or through terrestrial radiation during continental shelf emersion. Here, the Isthmus of Kra played a central part as biogeographical barriers that then separated it from the remaining continental populations.

A new subspecies identification and population study of the Asian small-clawed otter (Aonyx cinereus) in Malay Peninsula and southern Thailand based on fecal DNA method

Three species of otter can be found throughout Malay Peninsula: Aonyx cinereus, Lutra sumatrana, and Lutrogale perspicillata. In this study, we focused on the A. cinereus population that ranges from the southern and the east coast to the northern regions of Malay Peninsula up to southern Thailand to review the relationships between the populations based on the mitochondrial D-loop region. Forty-eight samples from six populations were recognized as Johor, Perak, Terengganu, Kelantan, Ranong, and Thale Noi. Among the 48 samples, 33 were identified as A. cinereus, seven as L. sumatrana, and eight as L. perspicillata. Phylogenetically, two subclades formed for A. cinereus. The first subclade grouped all Malay Peninsula samples except for samples from Kelantan, and the second subclade grouped Kelantan samples with Thai sample. Genetic distance analysis supported the close relationships between Thai and Kelantan samples compared to the samples from Terengganu and the other Malaysian states. A minimum-spanning network showed that Kelantan and Thailand formed a haplogroup distinct from the other populations. Our results show that Thai subspecies A. cinereus may have migrated to Kelantan from Thai mainland. We also suggest the classification of a new subspecies from Malay Peninsula, the small-clawed otter named A. cinereus kecilensis.

Phylogeography and demographic history of Lacerta lepida in the Iberian Peninsula: multiple refugia, range expansions and secondary contact zones

BACKGROUND

The Iberian Peninsula is recognized as an important refugial area for species survival and diversification during the climatic cycles of the Quaternary. Recent phylogeographic studies have revealed Iberia as a complex of multiple refugia. However, most of these studies have focused either on species with narrow distributions within the region or species groups that, although widely distributed, generally have a genetic structure that relates to pre-Quaternary cladogenetic events. In this study we undertake a detailed phylogeographic analysis of the lizard species, Lacerta lepida, whose distribution encompasses the entire Iberian Peninsula. We attempt to identify refugial areas, recolonization routes, zones of secondary contact and date demographic events within this species.

RESULTS

Results support the existence of 6 evolutionary lineages (phylogroups) with a strong association between genetic variation and geography, suggesting a history of allopatric divergence in different refugia. Diversification within phylogroups is concordant with the onset of the Pleistocene climatic oscillations. The southern regions of several phylogroups show a high incidence of ancestral alleles in contrast with high incidence of recently derived alleles in northern regions. All phylogroups show signs of recent demographic and spatial expansions. We have further identified several zones of secondary contact, with divergent mitochondrial haplotypes occurring in narrow zones of sympatry.

CONCLUSIONS

The concordant patterns of spatial and demographic expansions detected within phylogroups, together with the high incidence of ancestral haplotypes in southern regions of several phylogroups, suggests a pattern of contraction of populations into southern refugia during adverse climatic conditions from which subsequent northern expansions occurred. This study supports the emergent pattern of multiple refugia within Iberia but adds to it by identifying a pattern of refugia coincident with the southern distribution limits of individual evolutionary lineages. These areas are important in terms of long-term species persistence and therefore important areas for conservation.

High levels of population subdivision in a morphologically conserved Mediterranean toad (Alytes cisternasii) result from recent, multiple refugia: evidence from mtDNA, microsatellites and nuclear genealogies

Pleistocene glaciations often resulted in differentiation of taxa in southern European peninsulas, producing the high levels of endemism characteristic of these regions (e.g. the Iberian Peninsula). Despite their small ranges, endemic species often exhibit high levels of intraspecific differentiation as a result of a complex evolutionary history dominated by successive cycles of fragmentation, expansion and subsequent admixture of populations. Most evidence so far has come from the study of species with an Atlantic distribution in northwestern Iberia, and taxa restricted to Mediterranean-type habitats remain poorly studied. The Iberian Midwife toad (Alytes cisternasii) is a morphologically conserved species endemic to southwestern and central Iberia and a typical inhabitant of Mediterranean habitats. Applying highly variable genetic markers from both mitochondrial and nuclear genomes to samples collected across the species' range, we found evidence of high population subdivision within A. cisternasii. Mitochondrial haplotypes and microsatellites show geographically concordant patterns of genetic diversity, suggesting population fragmentation into several refugia during Pleistocene glaciations followed by subsequent events of geographical and demographic expansions with secondary contact. In addition, the absence of variation at the nuclear beta-fibint7 and Ppp3caint4 gene fragments suggests that populations of A. cisternasii have been recurrently affected by episodes of extinction and recolonization, and that documented patterns of population subdivision are the outcome of recent and multiple refugia. We discuss the evolutionary history of the species with particular interest in the increasing relevance of Mediterranean refugia for the survival of genetically differentiated populations during the Pleistocene glaciations as revealed by studies in co-distributed taxa.

The coup de grâce for the nested clade phylogeographic analysis?

Nested clade phylogeographic analysis (NCPA) has become a popular method for reconstructing the history of populations across species ranges. Ever since its invention in 1995, criticisms have been formulated, but the method, which has been regularly updated, continues to attract investigators. Molecular Ecology has published a large fraction of the literature on the topic--both pro and con. A recent study by Panchal and Beaumont (2007) finally allows a precise evaluation of the method by developing software that automates the somewhat complicated NCPA procedure. Using simulations of random-mating populations, Panchal and Beaumont find a high frequency of false-positives with their automated NCPA procedure (over 75%). These findings, which echo and amplify earlier warnings, appear serious enough to suggest to researchers to await further evaluation of the method. Although no other all-encompassing method such as the NCAP currently exists to evaluate phylogeographic data sets, researchers have many alternative methods to test ever more refined hypotheses.

Reading the history of a hybrid fish complex from its molecular record

Squalius alburnoides is a widely distributed intergeneric hybrid complex with fish of both sexes, varying ploidy levels and proportions of the parental genomes. Its dispersal routes were here delineated and framed by the reconstruction of the phylogeny and phylogeography of other Squalius with which it hybridizes, based on the available data on the paleohydrographical history of the Iberian Peninsula. Results based on sequences of cytochrome b and beta-actin genes showed that: proto-Squalius pyrenaicus originated at least five species as it dispersed throughout the Iberian Peninsula in the Mio-Pliocene; the S. alburnoides complex likely had a single origin in the bulk of Iberia, in the Upper Tagus/Guadiana area, when hydrographical rearrangements allowed the contact between its ancestors (around 700,000 years ago); interspecific crosses allowed the introgression of mitochondrial and nuclear genes of S. alburnoides in allopatric species/populations of other Squalius and vice-versa; and reconstituted S. alburnoides non-hybrid males may contribute to the replacement of the typical mtDNA of the complex (in the populations where they occur, crosses with females of other Squalius seem to have been especially frequent). A number of dispersal events and colonization routes are proposed.

Phylogeography of a ground beetle species in La Gomera (Canary Islands): the effects of landscape topology and population history

Paraeutrichopus pecoudi (Coleoptera, Carabidae) is a species endemic to the mountain laurel forests of La Gomera in the Canary archipelago. La Gomera is an island of volcanic origin (9.4 mya), where a well-preserved laurel forest is partly fragmented by valleys and ridges. Historically restricted gene flow with isolation by distance are deduced at different clade levels for P. pecoudi populations, an expected scenario arising from landscape discontinuity, which has caused local isolation of these flightless insects. Methodological considerations were investigated for nested clade phylogeographic analysis (NCPA), comparing the results obtained using either (a) pairwise distances between collecting sites measured as surface (route) distances or (b) standard geographical distances calculated from latitude-longitude coordinates. Some differences were found in the NCPA statistics and associated inferences at the higher clade levels in relation to the assumed distance criteria. Coalescent simulations and posterior automated NCPAs assuming different geographical distances were performed to test the robustness of the method when considering the distances used. These analyses showed no significant differences in NCPA conclusions, and the scenario of restricted gene flow was recovered using both distance calculators at the same rate, although the parameter statistics were slightly different. We suggest that in landscapes with extreme topography, geographical distances separating populations are certainly underestimated by standard techniques. More complex and explicit descriptions of the potential dispersal of terrestrial organisms should be explored for implementation in statistical phylogeography in such cases.

Mitochondrial DNA phylogeography and mating compatibility reveal marked genetic structuring and speciation in the NE Atlantic bryozoan Celleporella hyalina

The marine bryozoan Celleporella hyalina is a species complex composed of many highly divergent and mostly allopatric genetic lineages that are reproductively isolated but share a remarkably similar morphology. One such lineage commonly encrusts macroalgae throughout the NE Atlantic coast. To explore the processes leading to geographical diversification, reproductive isolation and speciation in this taxon, we (i) investigated NE Atlantic C. hyalina mitochondrial DNA phylogeography, and (ii) used breeding trials between geographical isolates to ascertain reproductive isolation. We find that haplotype diversity is geographically variable and there is a strong population structure, with significant isolation by distance. NE Atlantic C. hyalina is structured into two main parapatric lineages that appear to have had independent Pleistocene histories. Range expansions have resulted in two contact zones in Spain and W Ireland. Lineage 1 is found from Ireland to Spain and has low haplotype diversity, with closely related haplotypes, suggesting a recent population expansion into the Irish Sea, S Ireland, S England and Spain. Lineage 2 is found from Iceland to Spain and has high haplotype diversity. Complete reproductive isolation was found between some geographical isolates representing both lineages, whereas it was incomplete or asymmetric between others, suggesting these latter phylogeographical groups probably represent incipient species. The phylogeographical distribution of NE Atlantic C. hyalina does not fall easily into a pattern of southern refugia, and we discuss likely differences between terrestrial and marine system responses to Pleistocene glacial cycles.

Molecular biogeography of Europe: Pleistocene cycles and postglacial trends

The climatic cycles with subsequent glacial and intergalcial periods have had a great impact on the distribution and evolution of species. Using genetic analytical tools considerably increased our understanding of these processes. In this review I therefore give an overview of the molecular biogeography of Europe. For means of simplification, I distinguish between three major biogeographical entities: (i) "Mediterranean" with Mediterranean differentiation and dispersal centres, (ii) "Continental" with extra-Mediterranean centres and (iii) "Alpine" and/or "Arctic" with recent alpine and/or arctic distribution patterns. These different molecular biogeographical patterns are presented using actual examples. Many "Mediterranean" species are differentiated into three major European genetic lineages, which are due to glacial isolation in the three major Mediterranean peninsulas. Postglacial expansion in this group of species is mostly influenced by the barriers of the Pyrenees and the Alps with four resulting main patterns of postglacial range expansions. However, some cases are known with less than one genetic lineage per Mediterranean peninsula on the one hand, and others with a considerable genetic substructure within each of the Mediterranean peninsulas, Asia Minor and the Maghreb. These structures within the Mediterranean sub-centres are often rather strong and in several cases even predate the Pleistocene. For the "Continental" species, it could be shown that the formerly supposed postglacial spread from eastern Palearctic expansion centres is mostly not applicable. Quite the contrary, most of these species apparently had extra-Mediterranean centres of survival in Europe with special importance of the perialpine regions, the Carpathian Basin and parts of the Balkan Peninsula. In the group of "Alpine" and/or "Arctic" species, several molecular biogeographical patterns have been found, which support and improve the postulates based on distribution patterns and pollen records. Thus, genetic studies support the strong linkage between southwestern Alps and Pyrenees, northeastern Alps and Carpathians as well as southeastern Alps and the Dinaric mountain systems, hereby allowing conclusions on the glacial distribution patterns of these species. Furthermore, genetic analyses of arctic-alpine disjunct species support their broad distribution in the periglacial areas at least during the last glacial period. The detailed understanding of the different phylogeographical structures is essential for the management of the different evolutionary significant units of species and the conservation of their entire genetic diversity. Furthermore, the distribution of genetic diversity due to biogeographical reasons helps understanding the differing regional vulnerabilities of extant populations.

Molecular biogeography of Europe: Pleistocene cycles and postglacial trends

The climatic cycles with subsequent glacial and intergalcial periods have had a great impact on the distribution and evolution of species. Using genetic analytical tools considerably increased our understanding of these processes. In this review I therefore give an overview of the molecular biogeography of Europe. For means of simplification, I distinguish between three major biogeographical entities: (i) "Mediterranean" with Mediterranean differentiation and dispersal centres, (ii) "Continental" with extra-Mediterranean centres and (iii) "Alpine" and/or "Arctic" with recent alpine and/or arctic distribution patterns. These different molecular biogeographical patterns are presented using actual examples.

Many "Mediterranean" species are differentiated into three major European genetic lineages, which are due to glacial isolation in the three major Mediterranean peninsulas. Postglacial expansion in this group of species is mostly influenced by the barriers of the Pyrenees and the Alps with four resulting main patterns of postglacial range expansions. However, some cases are known with less than one genetic lineage per Mediterranean peninsula on the one hand, and others with a considerable genetic substructure within each of the Mediterranean peninsulas, Asia Minor and the Maghreb. These structures within the Mediterranean sub-centres are often rather strong and in several cases even predate the Pleistocene.

For the "Continental" species, it could be shown that the formerly supposed postglacial spread from eastern Palearctic expansion centres is mostly not applicable. Quite the contrary, most of these species apparently had extra-Mediterranean centres of survival in Europe with special importance of the perialpine regions, the Carpathian Basin and parts of the Balkan Peninsula.

In the group of "Alpine" and/or "Arctic" species, several molecular biogeographical patterns have been found, which support and improve the postulates based on distribution patterns and pollen records. Thus, genetic studies support the strong linkage between southwestern Alps and Pyrenees, northeastern Alps and Carpathians as well as southeastern Alps and the Dinaric mountain systems, hereby allowing conclusions on the glacial distribution patterns of these species. Furthermore, genetic analyses of arctic-alpine disjunct species support their broad distribution in the periglacial areas at least during the last glacial period.

The detailed understanding of the different phylogeographical structures is essential for the management of the different evolutionary significant units of species and the conservation of their entire genetic diversity. Furthermore, the distribution of genetic diversity due to biogeographical reasons helps understanding the differing regional vulnerabilities of extant populations.

Contrasting patterns of population subdivision and historical demography in three western Mediterranean lizard species inferred from mitochondrial DNA variation

Pleistocene climatic oscillations were a major force shaping genetic variability in many taxa. We analyse the relative effects of the ice ages across a latitudinal gradient in the Western Mediterranean region, testing two main predictions: (i) species with historical distributions in northern latitudes should have experienced greater loss of suitable habitat, resulting in higher extinction of historical lineages than species distributed in southern latitudes, where the effects of the ice ages were not as drastic. This would be reflected in the observation of lower diversity and number of differentiated lineages in northern areas. (ii) a signature of demographic expansion following the climate amelioration should be obvious in northern species, whereas in the south evidence of long-term effective population size stability should be observed. We used as models three species of wall lizards (Podarcis bocagei, Podarcis carbonelli and Podarcis vaucheri) that replace each other along the study area. We investigated the patterns of mitochondrial DNA diversity and subdivision and obtained demographic parameter estimates for each species. Our results suggest that P. bocagei, the northernmost species, bears low genetic diversity, a shallow coalescent history and marks of a demographic expansion. In contrast, P. vaucheri, the species with a southernmost distribution, shows deeper coalescence events, complex geographical substructure and no evidence for population growth. The species with an intermediate distribution, P. carbonelli, shows average levels of diversity, substructure and population growth. Taken together, these results conform to our main predictions and are explained by a differential influence of the ice ages on distinct latitudes.

Mitochondrial DNA phylogeography of Lissotriton boscai (Caudata, Salamandridae): evidence for old, multiple refugia in an Iberian endemic

In Europe, southern peninsulas served as refugia during cold periods in the Pleistocene, acting both as centres of origin of endemisms and as sources from which formerly glaciated areas were recolonized during interglacial periods. Previous studies have revealed that within the main refugial areas, intraspecific lineages often survived in allopatric refugia. We analysed two mitochondrial markers (nad4, control region, approximately 1.4 kb) in 103 individuals representing the entire distribution of Lissotriton boscai, a newt endemic to the western Iberian Peninsula. We inferred the evolutionary history of the species through phylogenetic, phylogeographic and historical demographic analyses. The results revealed unexpected, deep levels of geographically structured genetic variability. We identified two main evolutionary lineages, each containing three well-supported clades. The first historical split involved populations from central-southwestern coastal Portugal and the ancestor of all the remaining populations around 5.8 million years ago. Both lineages were subsequently fragmented into different population groups between 2.5 and 1.2 million years ago. According to nested clade analysis, at lower hierarchical levels the patterns suggest restricted gene flow with isolation by distance, whereas at higher levels the clades exhibit signatures of contiguous range expansion. Bayesian Skyline Plots show recent bottlenecks, followed by demographic expansions in all lineages. The significant genetic structure found is consistent with long-term survival of populations in allopatric refugia, supporting the 'refugia-within-refugia' scenario for southern European peninsulas. The comparison of our results with other co-distributed species highlights the generality of this hypothesis for the Iberian herpetofauna and suggests that Mediterranean refuges had more relevance for the composition and distribution of present biodiversity patterns than currently acknowledged. We briefly discuss the taxonomic and conservation implications of our results.

A molecular assessment of phylogenetic relationships and lineage accumulation rates within the family Salamandridae (Amphibia, Caudata)

We examine phylogenetic relationships among salamanders of the family Salamandridae using approximately 2700 bases of new mtDNA sequence data (the tRNALeu, ND1, tRNAIle, tRNAGln, tRNAMet, ND2, tRNATrp, tRNAAla, tRNAAsn, tRNACys, tRNATyr, and COI genes and the origin for light-strand replication) collected from 96 individuals representing 61 of the 66 recognized salamandrid species and outgroups. Phylogenetic analyses using maximum parsimony and Bayesian analysis are performed on the new data alone and combined with previously reported sequences from other parts of the mitochondrial genome. The basal phylogenetic split is a polytomy of lineages ancestral to (1) the Italian newt Salamandrina terdigitata, (2) a strongly supported clade comprising the "true" salamanders (genera Chioglossa, Mertensiella, Lyciasalamandra, and Salamandra), and (3) a strongly supported clade comprising all newts except S. terdigitata. Strongly supported clades within the true salamanders include monophyly of each genus and grouping Chioglossa and Mertensiella as the sister taxon to a clade comprising Lyciasalamandra and Salamandra. Among newts, genera Echinotriton, Pleurodeles, and Tylototriton form a strongly supported clade whose sister taxon comprises the genera Calotriton, Cynops, Euproctus, Neurergus, Notophthalmus, Pachytriton, Paramesotriton, Taricha, and Triturus. Our results strongly support monophyly of all polytypic newt genera except Paramesotriton and Triturus, which appear paraphyletic, and Calotriton, for which only one of the two species is sampled. Other well-supported clades within newts include (1) Asian genera Cynops, Pachytriton, and Paramesotriton, (2) North American genera Notophthalmus and Taricha, (3) the Triturus vulgaris species group, and (4) the Triturus cristatus species group; some additional groupings appear strong in Bayesian but not parsimony analyses. Rates of lineage accumulation through time are evaluated using this nearly comprehensive sampling of salamandrid species-level lineages. Rate of lineage accumulation appears constant throughout salamandrid evolutionary history with no obvious fluctuations associated with origins of morphological or ecological novelties.

Phylogeography of postglacial range expansion in Nigronia serricornis Say (Megaloptera: Corydalidae)

We have examined the effects of post-Wisconsinan glacial range expansion on the phylogeography of the saw-combed fishfly, Nigronia serricornis Say (Megaloptera: Corydalidae) because aquatic insects are under-represented in postglacial studies (and in phylogeography in general), and because the effects of ecological degradation on the population genetics of environmental indicator species like N. serricornis cannot be measured unless the underlying phylogeography is understood. Sequence data from a 630-base fragment of the mitochondrial cytochrome oxidase I (COI) gene were subjected to amova and nested clade analysis for 30 populations (n = 344) of N. serricornis. Both the amova and nested clade analysis revealed substantial population structure; 44.4% of the variance occured among populations. Three northward migrations are apparent: one from Tennessee into Illinois, Indiana, Wisconsin, Michigan and Ontario, a second that radiated eastward from Pennsylvania, and a third that moved along the coast from North Carolina into Connecticut, Massachusetts, Maine and then into New York. The latter two of these migrations were the result of contiguous range expansions, while the former expansion, out of Tennessee, appears to have been rapid with little gene flow from the source population. Additional clades included a group of haplotypes in central Kentucky that appear to have expanded along preglacial drainages, and clades in North Carolina and Georgia that have remained centrally located. Haplotype diversity decreased from south to north, a pattern that has been widely reported for animal and plant populations that expanded with the retreat of the Wisconsinan glaciation.

Phylogeography of the montane caddisflyDrusus discolor: evidence for multiple refugia and periglacial survival

We studied the genetic population structure and phylogeography of the montane caddisfly Drusus discolor across its entire range in central and southern Europe. The species is restricted to mountain regions and exhibits an insular distribution across the major mountain ranges. Mitochondrial sequence data (COI) of 254 individuals from the entire species range is analysed to reveal population genetic structure. The data show little molecular variation within populations and regions, but distinct genetic differentiation between mountain ranges. Most populations are significantly differentiated based on F(ST) and exact tests of population differentiation and most haplotypes are unique to a single mountain range. Phylogenetic analyses reveal deep divergence between geographically isolated lineages. Combined, these results suggest that past fragmentation is the prominent process structuring the populations across Europe. We use tests of selective neutrality and mismatch distributions, to study the demographic population history of regions with haplotype overlap. The high level of genetic differentiation between mountain ranges and estimates of demographic history provide evidence for the existence of multiple glacial refugia, including several in central Europe. The study shows that these aquatic organisms reacted differently to Pleistocene cooling than many terrestrial species. They persisted in numerous refugia over multiple glacial cycles, allowing many local endemic clades to form.

Phylogeography of ground tit (Pseudopodoces humilis) based on mtDNA: evidence of past fragmentation on the Tibetan Plateau

Pseudopodoces humilis, a long misclassified terrestrial tit, is the only species of parid whose distribution is limited to treeless terrain and endemic to the Tibetan Plateau. We revealed the phylogeographic structure of the species by using mitochondrial control region, as well as comparing morphological characters. The distinct geographic distributions of two major clades suggest spatial and temporal separations that coincide with important climatic and paleogeographic changes following the uplift of the Tibetan Plateau. Population expansion was inferred for the population at the platform of the Plateau 0.17 million years before present (Ma B.P.), and restricted gene flow with isolation by distance was detected within this region, congruent with expansion occurring after the extensive glacial period. A significant decrease in body size with decreasing altitude was found, possibly indicating selection for larger-sized birds at higher altitude.

From descriptive to predictive distribution models: a working example with Iberian amphibians and reptiles

Background

Aim of the study was to identify the conditions under which spatial-environmental models can be used for the improved understanding of species distributions, under the explicit criterion of model predictive performance. I constructed distribution models for 17 amphibian and 21 reptile species in Portugal from atlas data and 13 selected ecological variables with stepwise logistic regression and a geographic information system. Models constructed for Portugal were extrapolated over Spain and tested against range maps and atlas data.

Results

Descriptive model precision ranged from 'fair' to 'very good' for 12 species showing a range border inside Portugal ('edge species', kappa (k) 0.35–0.89, average 0.57) and was at best 'moderate' for 26 species with a countrywide Portuguese distribution ('non-edge species', k = 0.03–0.54, average 0.29). The accuracy of the prediction for Spain was significantly related to the precision of the descriptive model for the group of edge species and not for the countrywide species. In the latter group data were consistently better captured with the single variable search-effort than by the panel of environmental data.

Conclusion

Atlas data in presence-absence format are often inadequate to model the distribution of species if the considered area does not include part of the range border. Conversely, distribution models for edge-species, especially those displaying high precision, may help in the correct identification of parameters underlying the species range and assist with the informed choice of conservation measures.

Genetic differentiation of the marbled white butterfly, Melanargia galathea, accounts for glacial distribution patterns and postglacial range expansion in southeastern Europe

Isolation of Mediterranean species in the southern European peninsulas during the cold glacial phases often resulted in differentiation of several genetic lineages confined to the respective peninsulas. However, whilst there is good genetic evidence for multiple refugia in Iberia, there are only limited data available for the Balkans. Therefore, we wish to examine the hypothesis of a strong genetic structuring within southeastern Europe for the existence of multiple Balkan differentiation centres and/or several leading edges. As a model we use the marbled white butterfly, Melanargia galathea. We studied 18 allozyme loci of 564 individuals from 16 populations distributed over a large part of southeastern Europe. The single populations showed moderately high genetic diversity and no northward decline of genetic diversity was detected. The overall genetic differentiation between populations was considerable (F(ST) 7.0%). Cluster analysis discriminated three genetic groups: (i) a western flank in the former Yugoslavia, parts of eastern Austria and Hungary; (ii) an eastern flank with populations from Bulgaria and Romania (south of the southern Carpathians and eastern Carpathians); and (iii) the eastern Carpathian Basin. Hierarchical variance analysis distributed 53% of the variance among populations between these three groups. One sample from the Greek-Bulgarian border clustered within the eastern flank, but showed some tendency towards the eastern Carpathian Basin populations. Two populations from Carinthia clustered together with the eastern Carpathian Basin ones and a population from Styria showed an intermediate genetic composition between the three groups. Most probably, the eastern and the western flank groups are due to postglacial range expansion from the northeastern and the northwestern edges of the glacial differentiation centre (so-called leading edges). The eastern Carpathian Basin group may have resulted from postglacial expansion from northern Greece through valley systems of the central Balkan peninsula, maybe even expanding westwards north of the Balkan mountains reaching some parts of eastern Austria (e.g. Carinthia). Therefore, the Balkanic refugium of M. galathea may or may not have been continuous along the coastal areas of the Mediterranean, but must have been strongly genetically structured.

Life on the edge: the long-term persistence and contrasting spatial genetic structure of distinct brown trout life histories at their ecological limits

The Atlantic Iberian brown trout is at the southwestern limit of its distribution. At this ecological edge, which was once a glacial refugia, anadromy becomes less common as increased water temperatures restricted populations closer to the headwaters. We examined 847 individuals from 20 populations from throughout this region and assessed spatial genetic structure using 11 protein and four microsatellite loci. The higher levels of heterozygosity and allelic diversity north of the southernmost limit of anadromy (SLA), as well as an isolation-by-distance model of population structure, likely influenced by the anadromous forms, suggest that more stable demographic conditions existed over time in this region. Populations south of the SLA were highly differentiated given the restricted size of the area (protein F(ST) = 0.16 in the north and 0.63 in the south of the SLA; microsatellite F(ST) = 0.18 in the north and 0.70 in the south of the SLA). The low levels of heterozygosity and the pattern of southward allele depletion in resident populations is indicative of fragmentation, caused by stressful ecological conditions that reduced the anadromy (restricting gene flow) and the effective population sizes (higher genetic drift), which, in combination, dramatically decreased within-population genetic variation and increased among-population genetic variation. The higher gene diversity north of the SLA does not reflect ancestry but rather the signature of a population size expansion, as evidence suggest the persistence of older populations (with several private alleles) south of the SLA. These data support a scenario that demonstrates how contemporary events (critical ecological conditions) can moderate historical influences, suggesting that careful interpretation of the evolutionary history of glacial refugia is necessary, especially where populations persisted for a long time but not always with optimal ecological conditions. These peripheral populations are of high conservation value and should be preserved to help conserve the future potential of the species.

Extensive intraspecific polymorphism detected by SSCP at the nuclear C-mos gene in the endemic Iberian lizard Lacerta schreiberi

C-mos is a highly conserved intronless gene that has proved useful in the analysis of ancient phylogenetic relationships within vertebrates. We selected the Iberian endemic Schreiber's green lizard (Lacerta schreiberi) that persisted in allopatric refugia since the late Pliocene to investigate the utility of the C-mos nuclear gene for intraspecific phylogeographic studies. Our combination of DNA sequencing with the high resolving power of single-strand conformational polymorphism (SSCP) effectively discriminated four common alleles showing strong population structuring (F(ST) = 0.46). In addition, reconstruction of allele phylogenetic relationships further improved our understanding of C-mos spatial patterns of variation and allowed a comparison with previously described mitochondrial DNA data. Finally, limited sequencing of an extended C-mos fragment in six additional Lacerta species showed extensive polymorphism, to our knowledge representing a rare example of variation in a highly conserved nuclear gene.

Genetic exchange across a hybrid zone within the Iberian endemic golden-striped salamander, Chioglossa lusitanica

The study of hybrid zones resulting from Pleistocene vicariance is central in examining the potential of genetically diverged evolutionary units either to introgress and merge or to proceed with further isolation. The hybrid zone between two mitochondrial lineages of Chioglossa lusitanica is located near the Mondego River in Central Portugal. We used mitochondrial and nuclear diagnostic markers to conduct a formal statistical analysis of the Chioglossa hybrid zone in the context of tension zone theory. Key results are: (i) cline centres are not coincident for all markers, with average widths of ca. 2-15 km; (ii) heterozygote deficit was not observed across loci near the transect centre; (iii) associations of parental allele combinations ('linkage disequilibrium'R) were not detected either across loci or across the transect. These observations suggest that the Chioglossa hybrid zone is not a tension zone with strong selection against hybrids but instead one shaped mostly by neutral mixing. The patterns uncovered suggest a complex history of populations over a small scale that may be common in southern Pleistocene refugia.

Phylogeographic evidence for the existence of an ancient biogeographic barrier: the Isthmus of Kra Seaway

Biogeographic boundaries are characterised by distinct faunal and floral assemblages restricted on either side, but patterns among groups of taxa often vary and may not be discrete. Historical biogeography as a consequence, while providing crucial insights into the relationship between biological diversity and earth history, has some limitations. Patterns of intraspecific molecular variation, however, may show unambiguous evidence for such historical divides, and can be used to test competing biogeographic hypotheses (often based on the dispersal-vicariance debate). Here, we utilise this method to test the hypothesis that a major biogeographic transition zone between the Sundaic and Indochinese biotas, located just north of the Isthmus of Kra in SE Asia, is the result of Neogene marine transgressions that breached the Isthmus in two locations for prolonged periods of time (>1 million year duration). Phylogeographic analyses of a freshwater decapod crustacean, the giant freshwater prawn Macrobrachium rosenbergii, strongly supports the historical existence of the more northerly postulated seaway. Results presented here highlight the power of utilising intraspecific molecular variation in testing biogeographical hypotheses.

The structure of biodiversity - insights from molecular phylogeography

DNA techniques, analytical methods and palaeoclimatic studies are greatly advancing our knowledge of the global distribution of genetic diversity, and how it evolved. Such phylogeographic studies are reviewed from Arctic, Temperate and Tropical regions, seeking commonalities of cause in the resulting genetic patterns. The genetic diversity is differently patterned within and among regions and biomes, and is related to their histories of climatic changes. This has major implications for conservation science.

Genetic consequences of climatic oscillations in the Quaternary

An appreciation of the scale and frequency of climatic oscillations in the past few million years is modifying our views on how evolution proceeds. Such major events caused extinction and repeated changes in the ranges of those taxa that survived. Their spatial effects depend on latitude and topography, with extensive extinction and recolonization in higher latitudes and altitudinal shifts and complex refugia nearer the tropics. The associated population dynamics varied with life history and geography, and the present genetic constitution of the populations and species carry attenuated signals of these past dynamics. Phylogeographic studies with DNA have burgeoned recently and studies are reviewed from the arctic, temperate and tropical regions, seeking commonalities of cause in the resulting genetic patterns. Arctic species show distinct shallow genetic clades with common geographical boundaries. Thus Beringia is distinct phylogeographically, but its role as a refugial source is complex. Arctic taxa do not show the common genetic pattern of southern richness and northern purity in north-temperate species. Temperate refugial regions in Europe and North America show relatively deep DNA divergence for many taxa, indicating their presence over several Ice Ages, and suggesting a mode of speciation by repeated allopatry. DNA evidence indicates temperate species in Europe had different patterns of postglacial colonization across the same area and different ones in previous oscillations, whereas the northwest region of North America was colonized from the north, east and south. Tropical montane regions contain deeply diverged lineages, often in a relatively small geographical area, suggesting their survival there from the Pliocene. Our poor understanding of refugial biodiversity would benefit from further combined fossil and genetic studies.

Statistical phylogeography: methods of evaluating and minimizing inference errors

Nested clade phylogeographical analysis (NCPA) has become a common tool in intraspecific phylogeography. To evaluate the validity of its inferences, NCPA was applied to actual data sets with 150 strong a priori expectations, the majority of which had not been analysed previously by NCPA. NCPA did well overall, but it sometimes failed to detect an expected event and less commonly resulted in a false positive. An examination of these errors suggested some alterations in the NCPA inference key, and these modifications reduce the incidence of false positives at the cost of a slight reduction in power. Moreover, NCPA does equally well in inferring events regardless of the presence or absence of other, unrelated events. A reanalysis of some recent computer simulations that are seemingly discordant with these results revealed that NCPA performed appropriately in these simulated samples and was not prone to a high rate of false positives under sampling assumptions that typify real data sets. NCPA makes a posteriori use of an explicit inference key for biological interpretation after statistical hypothesis testing. Alternatives to NCPA that claim that biological inference emerges directly from statistical testing are shown in fact to use an a priori inference key, albeit implicitly. It is argued that the a priori and a posteriori approaches to intraspecific phylogeography are complementary, not contradictory. Finally, cross-validation using multiple DNA regions is shown to be a powerful method of minimizing inference errors. A likelihood ratio hypothesis testing framework has been developed that allows testing of phylogeographical hypotheses, extends NCPA to testing specific hypotheses not within the formal inference key (such as the out-of-Africa replacement hypothesis of recent human evolution) and integrates intra- and interspecific phylogeographical inference.

PHYLOGENY AND PHYLOGEOGRAPHY OF THE LIOLAEMUS DARWINII COMPLEX (SQUAMATA: LIOLAEMIDAE): EVIDENCE FOR INTROGRESSION AND INCOMPLETE LINEAGE SORTING

Although mitochondrial DNA markers have several properties that make them suitable for phylogeographic studies, they are not free of difficulties. Phylogeographic inferences within and between closely related species can be mislead by introgression and retention of ancestral polymorphism. Here we combine different phylogenetic, phylogeographic, and population genetic methods to extract the maximum information from the Liolaemus darwinii complex. We estimate the phylogeographic structure of L. darwinii across most of its distributional range, and we then estimate relationships between L. darwinii and the syntopic species L. laurenti and L. grosseorum. Our results suggest that range expansion of these lineages brought them into secondary contact in areas where they are presently in syntopy. Here we present the first evidence for introgression in lizards from temperate South America (of L. danwinii mitochondrial DNA into L. laurenti and L. grosseorum), and for incomplete lineage sorting (between L. darwinii and L. laurenti). We show that a combination of methods can provide additional support for inferences derived from any single method and thus provide more robust interpretations and narrow the range of plausible hypotheses about mechanisms and processes of divergence. Additional studies are needed in this group of lizards and in other codistributed groups to determine if Pleistocene climatic changes could be a general factor influencing the evolutionary history of a regional biota.