Genetic consequences of climatic oscillations in the Quaternary

Phylogeography and environmental diversification of a highly adaptable marine amphipod, Gammarus duebeni

Genetic diversity and phylogeographic population structure in the gammarid amphipod, Gammarus duebeni, were investigated across its broad latitudinal distribution in the NE and NW Atlantic by analysis of mitochondrial DNA sequence. Gammarus duebeni has exceptional tolerance of salinity change and inhabits environments ranging from marine to freshwater. The longstanding debate on whether there are distinct marine and freshwater subspecies was assessed by sampling populations from sites characterized by different salinities. Our sequence data demonstrates that there are two major lineages, with little internal geographic structuring. Evidence is provided to suggest a pre-glacial divergence of these two clades, involving segregation between a region historically associated with the freshwater form and the majority of the marine localities on both sides of the Atlantic. A modern contact zone between the marine and freshwater forms is proposed in western Britain.

Evolution and genetic differentiation among relatives of Arabidopsis thaliana

Arabidopsis thaliana is one of the most intensively studied plant species. More recently, information is accumulating about its closest relatives, the former genus Cardaminopsis. A. thaliana diverged from these relatives, actually treated within three major lineages (Arabidopsis lyrata, Arabidopsis halleri, and Arabidopsis arenosa), approximately 5 mya. Significant karyotype evolution in A. thaliana with base chromosome number reduction from x=8 to x=5 might indicate and favor effective genetic isolation from these other species, although hybrids are occurring naturally and have been also constituted under controlled conditions. We tested the evolutionary significance to separate the x=5 from the x=8 lineage using DNA sequence data from the plastome and the nuclear ribosomal DNA based on an extensive, representative worldwide sampling of nearly all taxonomic entities. We conclude that (i) A. thaliana is clearly separated phylogenetically from the x=8 lineage, (ii) five major lineages outside A. thaliana can be identified (A. lyrata, A. arenosa, A. halleri, Arabidopsis croatica, and Arabidopsis pedemontana) together with Arabidopsis cebennensis, and (iii) centers of genetic and morphological diversity are mostly in congruence and are located close to the Balkans in Austria and Slovakia outside glaciated and permafrost regions with few notable exceptions.

Genetic structure of the widespread and common Mediterranean bryophyte Pleurochaete squarrosa (Brid.) Lindb. (Pottiaceae) - evidence from nuclear and plastidic DNA sequence variation and allozymes

The Mediterranean Basin as one the world's most biologically diverse regions provides an interesting area for the study of plant evolution and spatial structure in plant populations. The dioecious moss Pleurochaete squarrosa is a widespread and common bryophyte in the Mediterranean Basin. Thirty populations were sampled for a study on molecular diversity and genetic structure, covering most major islands and mainland populations from Europe and Africa. A significant decline in nuclear and chloroplast sequence and allozyme variation within populations from west to east was observed. While DNA sequence data showed patterns of isolation by distance, allozyme markers did not. Instead, their considerable interpopulation genetic differentiation appeared to be unrelated to geographic distance. Similar high values for coefficients of gene diversity (G(ST)) in all data sets provided evidence of geographic isolation and limited gene flow among populations (i) within islands, (ii) within mainland areas, and (iii) between islands and mainland. Notably, populations in continental Spain are strongly genetically isolated from all other investigated areas. Surprisingly, there was no difference in gene diversity and G(ST) between islands and mainland areas. Thus, we conclude that large Mediterranean islands may function as 'mainland' for bryophytes. This hypothesis and its implication for conservation biology of cryptogamic plants warrant further investigation. While sexually reproducing populations were found all over the Mediterranean Basin, high levels of multilocus linkage disequilibrium provide evidence of mainly vegetative propagation even in populations where sexual reproduction was observed.

Molecular evolution of satellite DNA repeats and speciation of lizards of the genus Darevskia (Sauria: Lacertidae)

Satellite DNA repeats were studied in Caucasian populations of 18 rock lizard species of the genus Darevskia. Four subfamilies (Caucasian Lacerta satellites (CLsat)I-IV) were identified, which shared 70%-75% sequence similarity. The distribution of CLsat subfamilies among the species was studied. All the species could be divided into at least 3 clades, depending on the content of CLsat subfamilies in each genome: "saxicola", "rudis", and "mixta" lizards. CLsatI was found in all studied species, but in very different quantities; the "saxicola" group contained this subfamily predominantly. The "rudis" group also contained CLsatIII, and the "mixta" group carried considerable amounts of CLsatII. The highest concentrations of CLsatI and CLsatII were detected in 2 ground lizards--D. derjugini and D. praticola, respectively. D. parvula predominantly carried CLsatIII. CLsatIV was found only in the Crimean species D. lindholmi. The distribution patterns of satellite subfamilies show possible postglacial speciation within the genus Darevskia. A hybrid origin of species that possess 2 or 3 CLsat subfamilies and important clarifications to the systematics of the genus are proposed.

Testing dispersal and cryptic diversity in a widely distributed groundwater amphipod (Niphargus rhenorhodanensis)

Theories about colonization and evolution in groundwater have assumed that the fragmented structure of groundwater strongly limits dispersal. The high number of endemic and allopatric species in groundwater supports this hypothesis, but the occurrence of widely distributed groundwater taxa calls into question its universality. These widely distributed taxa might also be sets of cryptic species because extreme conditions of life in groundwater promote cryptic diversity by inducing convergent morphological evolution. Niphargus rhenorhodanensis is a widely distributed and ubiquitous groundwater amphipod which supposedly colonized the Alps after Quaternary glaciations. We tested the dispersal and the cryptic species hypotheses within this species using a phylogeographic approach based on two mitochondrial genes (COI and 16S) and a nuclear gene (28S). Results support the view that poor dispersal is a main evolutionary factor in groundwater. All genes independently supported the existence of numerous cryptic and mostly allopatric units within N. rhenorhodanensis, indicating that its apparently wide distribution range is an artefact generated by cryptic diversity. We reject the hypothesis of a recent and global colonization of the Alps and argue that some N. rhenorhodanensis lineages probably survived glaciations near or within the Alps.

Phylogeography of the mountain chickadee (Poecile gambeli): diversification, introgression, and expansion in response to Quaternary climate change

Since the late 1990s, molecular techniques have fuelled debate about the role of Pleistocene glacial cycles in structuring contemporary avian diversity in North America. The debate is still heated; however, there is widespread agreement that the Pleistocene glacial cycles forced the repeated contraction, fragmentation, and expansion of the North American biota. These demographic processes should leave genetic 'footprints' in modern descendants, suggesting that detailed population genetic studies of contemporary species provide the key to elucidating the impact of the late Quaternary (late Pleistocene-Holocene). We present an analysis of mitochondrial DNA (mtDNA) variation in the mountain chickadee (Poecile gambeli) in an attempt to examine the genetic evidence of the impact of the late Quaternary glacial cycles. Phylogenetic analyses reveal two strongly supported clades of P. gambeli: an Eastern Clade (Rocky Mountains and Great Basin) and a Western Clade (Sierra Nevada and Cascades). Post-glacial introgression is apparent between these two clades in the Mono Lake region of Central California. Within the Eastern Clade there is evidence of isolation-by-distance in the Rocky Mountain populations, and of limited gene flow into and around the Great Basin. Coalescent analysis of genetic variation in the Western Clade indicates that northern (Sierra Nevada/Cascades) and southern (Transverse/Peninsular Ranges) populations have been isolated and evolving independently for nearly 60,000 years.

Genetic diversity and population structure of Lamiophlomis rotata (Lamiaceae), an endemic species of Qinghai-Tibet Plateau

Lamiophlomis rotata (Lamiaceae), a perennial medicinal herb, is endemic to the Qinghai-Tibet Plateau. A total of 188 individuals from eight natural populations of L. rotata in Qinghai-Tibet Plateau (four from Tibet, two from Yunnan, and two from Qinghai) were analyzed using intersimple sequence repeats (ISSR) and randomly amplified polymorphic DNA (RAPD) techniques. Our results revealed that the level of genetic variation in L. rotata was relatively high (P = 94.85%, I = 0.440 +/- 0.220, H(T) = 0.289 +/- 0.028). Three genetic groups corresponding to the three geographic regions were detected, suggesting significant geographic structure. Our results suggest that the highly structured geographic pattern found in L. rotata might represent diverging evolutionary processes associated with the uplifting of the Qinghai-Tibet Plateau and Quaternary climatic oscillations. These findings imply that as many populations as possible should be preserved in situ for the conservation of this species. Given their genetic variability and peripheral distribution, Qinghai and Yunnan populations should be assigned priority for conservation. Optimal harvesting strategies, domestication and tissue culture of L. rotata should be developed as soon as possible to guarantee its sustainable use.

Phylogeography of the Western Lyresnake (Trimorphodon biscutatus): testing aridland biogeographical hypotheses across the Nearctic-Neotropical transition

The Western Lyresnake (Trimorphodon biscutatus) is a widespread, polytypic taxon inhabiting arid regions from the warm deserts of the southwestern United States southward along the Pacific versant of Mexico to the tropical deciduous forests of Mesoamerica. This broadly distributed species provides a unique opportunity to evaluate a priori biogeographical hypotheses spanning two major distinct biogeographical realms (the Nearctic and Neotropical) that are usually treated separately in phylogeographical analyses. I investigated the phylogeography of T. biscutatus using maximum likelihood and Bayesian phylogenetic analysis of mitochondrial DNA (mtDNA) from across this species' range. Phylogenetic analyses recovered five well-supported clades whose boundaries are concordant with existing geographical barriers, a pattern consistent with a model of vicariant allopatric divergence. Assuming a vicariance model, divergence times between mitochondrial lineages were estimated using Bayesian relaxed molecular clock methods calibrated using geological information from putative vicariant events. Divergence time point estimates were bounded by broad confidence intervals, and thus these highly conservative estimates should be considered tentative hypotheses at best. Comparison of mtDNA lineages and taxa traditionally recognized as subspecies based on morphology suggest this taxon is comprised of multiple independent lineages at various stages of divergence, ranging from putative secondary contact and hybridization to sympatry of 'subspecies'.

Circumpolar phylogeography of Juncus biglumis (Juncaceae) inferred from AFLP fingerprints, cpDNA sequences, nuclear DNA content and chromosome numbers

We explored the circumpolar phylogeographic history of the arctic-alpine Juncus biglumis using amplified fragment length polymorphisms (AFLPs), sequences of cpDNA, relative nuclear DNA content and chromosome numbers. The analyses of the AFLP and cpDNA data gave congruent results and revealed three distinct clades. One of them, represented by a single population from the Taymyr peninsula in northern Siberia, had approximately fourfold larger genome size than the other samples and produced an AFLP pattern that was too aberrant to be analysed together with the rest of the data set. The two other clades represented different ploidy levels (2n = 60 and 120) as judged from chromosome counts of selected populations but differed only in c. 6% relative DNA content. Based on the AFLP and partly also on the cpDNA data, each of the two main clades was further subdivided into two well-supported subgroups. Three of the subgroups were widespread and exhibited largely overlapping distribution patterns. The fourth subgroup seems to be absent from the North Atlantic region and from western Siberia. We suggest that the four subgroups diverged during isolation in different glacial refugia during the Quaternary. Interestingly, individuals of both main clades were encountered in geographically close populations in eastern Greenland and even within a single population from Svalbard, indicating that both areas were colonised at least twice. The different genome sizes and ploidy levels strongly suggest that the three main clades represent distinct gene pools and act as cryptic species.

Phylogeographic patterns in Drosophila montana

The Drosophila virilis species group offers valuable opportunities for studying the roles of chromosomal re-arrangements and mating signals in speciation. The 13 species are divided into two subgroups, the montana and virilis 'phylads'. There is greater differentiation among species within the montana phylad in both karyotype and acoustic signals than exists among members of the virilis phylad. Drosophila montana is a divergent species which is included in the montana phylad. Here, we analyse the phylogeography of D. montana to provide a framework for understanding divergence of acoustic signals among populations. We analysed mitochondrial sequences corresponding to the cytochrome oxidase I and cytochrome oxidase II genes, as well as 16 microsatellite loci, from 108 lines of D. montana covering most of the species' range. The species shows a clear genetic differentiation between North American and Scandinavian populations. Microsatellite allele frequencies and mitochondrial DNA haplotypes gave significant FST values between populations from Canada, USA and Finland. A Bayesian analysis of population structure based on the microsatellite frequencies showed four genetically distinct groups, corresponding to these three populations plus a small sample from Japan. A network based on mitochondrial haplotypes showed two Finnish clades of very different shape and variability, and another clade with all sequences from North America and Japan. All D. montana populations showed evidence of demographic expansion but the patterns inferred by coalescent analysis differed between populations. The divergence times between Scandinavian and North American clades were estimated to range from 450,000 to 900,000 years with populations in Canada and the USA possibly representing descendants of different refugial populations. Long-term separation of D. montana populations could have provided the opportunity for differentiation observed in male signal traits, especially carrier frequency of the song, but relaxation of sexual selection during population expansion may have been necessary.

Genetic admixture of burbot (Teleostei: Lota lota) in Lake Constance from two European glacial refugia

The burbot, Lota lota, is the only freshwater species of the codfish family and has a Holarctic distribution. Pleistocene glaciations caused significant geographical differentiation in the past, but its life history characterized by winter spawning migrations over large distances is likely to homogenize populations by contemporary gene flow. We investigated the population genetic structure of 541 burbots from Lake Constance and adjacent Rhine and Danube tributaries in Europe using the entire mitochondrial DNA (mtDNA) control region and 11 microsatellites. Microsatellites revealed considerable population divergence (F(ST) = 0.26) and evidenced recent bottlenecks in two Central European rivers. In accordance to previous evidence two main phylogeographic lineages (Atlantic and Danubian) were found co-occurring at similar frequencies in Lake Constance, where they currently undergo random mating as indicated by microsatellites. The Danubian lineage contributed only a small proportion to the lake's mtDNA diversity, and probably expanded within the lake shortly after its formation approximately 10,000-15,000 BP. The larger Atlantic haplotype diversity suggested a population expansion older than the lake itself. Levels of admixture at microsatellite loci were less obvious due to their high variability, and coalescence methods were used to estimate past admixture proportions. Our results reinforce a model of a two-step colonization of Europe by burbot from an ancestral Danubian refuge, and confirm the persistence of a secondary Atlantic refuge, as proposed to exist for other freshwater fish. We conclude that the present-day burbot population in Lake Constance bears the genetic signature of both contemporary gene flow and historical separation events.

Microsatellite diversity and broad scale geographic structure in a model legume: building a set of nested core collection for studying naturally occurring variation in Medicago truncatula

BACKGROUND

Exploiting genetic diversity requires previous knowledge of the extent and structure of the variation occurring in a species. Such knowledge can in turn be used to build a core-collection, i.e. a subset of accessions that aim at representing the genetic diversity of this species with a minimum of repetitiveness. We investigate the patterns of genetic diversity and population structure in a collection of 346 inbred lines representing the breadth of naturally occurring diversity in the Legume plant model Medicago truncatula using 13 microsatellite loci distributed throughout the genome.

RESULTS

We confirm the uniqueness of all these genotypes and reveal a large amount of genetic diversity and allelic variation within this autogamous species. Spatial genetic correlation was found only for individuals originating from the same population and between neighbouring populations. Using a model-based clustering algorithm, we identified four main genetic clusters in the set of individuals analyzed. This stratification matches broad geographic regions. We also identified a set of "admixed" individuals that do not fit with this population structure scheme.

CONCLUSION

The stratification inferred is discussed considering potential historical events like expansion, refuge history and admixture between neighbouring groups. Information on the allelic richness and the inferred population structure are used to build a nested core-collection. The set of inbred lines and the core collections are publicly available and will help coordinating efforts for the study of naturally occurring variation in the growing Medicago truncatula community.

Niche breadth, competitive strength and range size of tree species: a trade-off based framework to understand species distribution

Understanding the mechanisms causing latitudinal gradients in species richness and species range size is a central issue in ecology, particularly in the current context of global climate change. Different hypotheses have been put forward to explain these patterns, emphasizing climatic variability, energy availability and competition. Here we show, using a comparative analysis controlling for phylogeny on 234 temperate/boreal tree species, that these hypotheses can be included into a single framework in an attempt to explain latitudinal gradients in species range size. We find that species tend to have larger ranges when (i) closer to the poles, (ii) successionally seral, (iii) having small and light seeds, and (iv) having short generations. The patterns can simply be explained by energy constraints associated with different life-history strategies. Overall, these findings shed a new light on our understanding of species distribution and biodiversity patterns, bringing new insights into underlying large-scale evolutionary processes.

Evolution of mitochondrial relationships and biogeography of Palearctic green toads (Bufo viridis subgroup) with insights in their genomic plasticity

Taxa involving three bisexually reproducing ploidy levels make green toads a unique amphibian system. We put a cytogenetic dataset from Central Asia in a molecular framework and apply phylogenetic and demographic methods to data from the entire Palearctic range. We study the mitochondrial relationships of diploids to infer their phylogeography and the maternal ancestry of polyploids. Control regions (and tRNAs between ND1 and ND2 in representatives) characterize a deeply branched assemblage of twelve haplotype groups, diverged since the Lower Miocene. Polyploidy has evolved several times: Central Asian tetraploids (B. oblongus, B. pewzowi) have at least two maternal origins. Intriguingly, the mitochondrial ancestor of morphologically distinctive, sexually reproducing triploid taxa (B. pseudoraddei) from Karakoram and Hindukush represents a different lineage. We report another potential case of bisexual triploid toads (B. zugmayeri). Identical d-loops in diploids and tetraploids from Iran and Turkmenistan, which differ in morphology, karyotypes and calls, suggest multiple origins and retained polymorphism and/or hybridization. A similar system involves diploids, triploids and tetraploids from Kyrgyzstan and Kazakhstan where green toads exemplify vertebrate genomic plasticity. A new form from Sicily and its African sister species (B. boulengeri) allow internal calibration and divergence time estimates for major clades. The subgroup may have originated in Eurasia rather than Africa since the earliest diverged lineages (B. latastii, B. surdus) and earliest fossils occur in Asia. We delineate ranges, contact and hybrid zones. Phylogeography, including one of the first non-avian datasets from Central Asian high mountains, reflects Quaternary climate and glaciation.

Quaternary diversification in a sexual Holarctic zooplankter, Daphnia galeata

The effects of Quaternary glacial range partitioning on the diversification of Holarctic biota remain unclear. Glacial refugial lineages may form vicariant species, hybrid products, or merge after secondary contact. Here, we assess the effects of Quaternary glaciation on a Holarctic sexual zooplankter, Daphnia galeata, with apparently marked dispersal capacity and a widespread hybrid lineage in the New World. We collected samples of this species from 148 Holarctic lakes, analysed the nuclear and mitochondrial gene sequences, and tested predictions for hypotheses that account for the origin and spread of the New World D. galeata. We detected five nuclear phylogroups and four mitochondrial phylogroups, most of which were restricted to either the New World or the Old World. The oldest mitochondrial phylogroup was restricted to Japan. One major mitochondrial clade was distributed throughout the Holarctic, but only four haplotypes were shared among continents, and analysis of molecular variance indicated significant structure at the continental level. Haplotype sharing among continents could largely be attributed to anthropogenic introductions. Mismatch distributions, haplotype networks, phylogenetic trees, longitudinal haplotype diversity erosion and coalescence analyses are consistent with colonization from an Old World and a New World refugium. Our nuclear and mitochondrial DNA sequence evidence supports the hypothesis that the New World D. galeata underwent introgression with Daphnia dentifera, with dispersal being enhanced by glaciation. We conclude that Quaternary glaciation had a pronounced effect on the diversification of a Holarctic sexual zooplankter.

Differential patterns of morphological and molecular hybridization between Fraxinus excelsior L. and Fraxinus angustifolia Vahl (Oleaceae) in eastern and western France

We examined large-scale patterns of morphology, genetic structure and ecological correlates of Fraxinus excelsior and the closely related species Fraxinus angustifolia in France, in order to determine the degree of hybridization between them. We sampled 24 populations in two putative hybrid zones (Loire and Saône), and five control populations of each species. We measured foliar characteristics of adult trees and used five nuclear microsatellites as molecular markers. Canonical discriminant analysis indicated that the two species differ in morphology, but that intermediate types are common in the Loire region but less frequent in the Saône region. Bayesian population assignment identified one F. angustifolia and two F. excelsior gene pools. Most Loire individuals clustered genetically with the F. angustifolia gene pool. In contrast, the Saône region presented individuals belonging mostly to F. excelsior pools, although the F. angustifolia type was frequent in certain populations. The lowest F(ST) values were found between the Loire and F. angustifolia controls that also exhibited no significant isolation by distance. The proportion of the F. angustifolia gene pool in each locality was negatively correlated with winter temperatures, suggesting that a cold climate may be limiting. Hybridization is probably favoured by the intermediate climatic conditions in the Loire region that allow both species to occur, but is somewhat hampered by the harsher winters in the Saône area where morphological introgression has apparently not yet occurred.

Molecular phylogeny and biogeography of Picea (Pinaceae): Implications for phylogeographical studies using cytoplasmic haplotypes

The center of diversity is not necessarily the place of origin, as has been established by many plant molecular phylogenies. Picea is a complicated but very important genus in coniferous forests of the Northern Hemisphere, with a high species diversity in Asia. Its phylogeny and biogeography were investigated here using sequence analysis of the paternally inherited chloroplast trnC-trnD and trnT-trnF regions and the maternally inherited mitochondrial nad5 intron 1. We found that the North American P. breweriana and P. sitchensis were basal to the other spruces that were further divided into three clades in the cpDNA phylogeny, and that the New World species harbored four of five mitotypes detected, including two ancestral ones and three endemics. These results, combined with biogeographic analyses using DIVA and MacClade and fossil evidence, suggest that Picea originated in North America, and that its present distribution could stem from two times of dispersal from North America to Asia by the Beringian land bridge, and then from Asia to Europe. Most of the northeastern Asian species and the European P. abies could arise from a recent radiation given the very low interspecific genetic differentiation and pure mitotype of them. Considering that the ancestral mtDNA polymorphism can be preserved in many descendant species, even distantly related ones, we suggest that more species, at least the closely related ones, should be sampled in the phylogeographical study using cytoplasmic haplotypes if possible. In addition, we also discussed the evolution and phylogenetic utility of morphological characters in Picea.

A tale of two flatties: different responses of two terrestrial flatworms to past environmental climatic fluctuations at Tallaganda in montane southeastern Australia

Comparative phylogeographic studies of animals with low mobility and/or high habitat specificity remain rare, yet such organisms may hold fine-grained palaeoecological signal. Comparisons of multiple, codistributed species can elucidate major historical events. As part of a multitaxon programme, mitochondrial cytochrome oxidase I (COI) variation was analysed in two species of terrestrial flatworm, Artioposthia lucasi and Caenoplana coerulea. We applied coalescent demographic estimators and nested clade analysis to examine responses to past, landscape-scale, cooling-drying events in a model system of montane forest (Tallaganda). Correspondence of haplotype groups in both species to previously proposed microbiogeographic regions indicates at least four refuges from cool, dry conditions. The region predicted to hold the highest quality refuges (the Eastern Slopes Region), is indicated to have been a long-term refuge in both species, but so are several other regions. Coalescent analyses suggest that populations of A. lucasi are declining, while C. coerulea is expanding, although stronger population substructure in the former could yield similar patterns in the data. The differences in spatial and temporal genetic variation in the two species could be explained by differences in ecological attributes: A. lucasi is predicted to have lower dispersal ability but may be better able to withstand cold conditions. Thus, different contemporary population dynamics may reflect different responses to recent (Holocene) climate warming. The two species show highly congruent patterns of catchment-based local genetic endemism with one another and with previously studied slime-mould grazing Collembola.

Population genetic structure of Arabidopsis lyrata in Europe

Population genetic theory predicts that the self-incompatible and perennial herb, Arabidopsis lyrata, will have a genetic structure that differs from the self-fertilizing, annual Arabidopsis thaliana. We quantified the genetic structure for eight populations of A. lyrata ssp. petraea in historically nonglaciated regions of central Europe. Analysis of 20 microsatellite loci for 344 individuals demonstrated that, in accordance with predictions, diploid populations had high genome-wide heterozygosity (H(O) = 0.48; H(E) = 0.52), high within-population diversity (83% of total) compatible with mutation-drift equilibrium, and moderate differentiation among populations (F(ST) = 0.17). Within a single population, the vast majority of genetic variability (92%) was found at the smallest spatial scale (< 3 m). Although there was no evidence of biparental inbreeding or clonal propagation at this scale (F(IS) = 0.003), significant fine-scale spatial autocorrelation indicated localized gene flow presumably due to gravity dispersed seeds (Sp = 0.018). Limited gene flow between isolated population clusters (regions) separated by hundreds of kilometres has given rise to an isolation by distance pattern of diversification, with low, but significant, differentiation among regions (F(ST) = 0.05). The maintenance of geographically widespread polymorphisms and uniformly high diversity throughout central Europe is consistent with periglacial survival of A. lyrata ssp. petraea north of the Alps in steppe-tundra habitats during the last glacial maximum. As expected of northern and previously glaciated localities, A. lyrata in Iceland was genetically less diverse and highly differentiated from central Europe (H(E) = 0.37; F(ST) = 0.27).

A northern glacial refugium for bank voles (Clethrionomys glareolus)

There is controversy and uncertainty on how far north there were glacial refugia for temperate species during the Pleistocene glaciations and in the extent of the contribution of such refugia to present-day populations. We examined these issues using phylogeographic analysis of a European woodland mammal, the bank vole (Clethrionomys glareolus). A Bayesian coalescence analysis indicates that a bank vole population survived the height of the last glaciation (approximately 25,000-10,000 years B.P.) in the vicinity of the Carpathians, a major central European mountain chain well north of the Mediterranean areas typically regarded as glacial refugia for temperate species. Parameter estimates from the fitted isolation with migration model show that the divergence of the Carpathian population started at least 22,000 years ago, and it was likely followed by only negligible immigration from adjacent regions, suggesting the persistence of bank voles in the Carpathians through the height of the last glaciation. On the contrary, there is clear evidence for gene flow out of the Carpathians, demonstrating the contribution of the Carpathian population to the colonization of Europe after the Pleistocene. These findings are consistent with data from animal and plant fossils recovered in the Carpathians and provide the clearest phylogeographic evidence to date of a northern glacial refugium for temperate species in Europe.

European colonization by the spined loach (Cobitis taenia) from Ponto-Caspian refugia based on mitochondrial DNA variation

In the last 20 years, new species, asexual reproduction, polyploidy and hybridization have all been reported within the genus Cobitis. An understanding of the current distribution and baseline phylogeographical history of 'true' nonhybrid Cobitis species is crucial in order to unravel these discoveries. In the present work, we investigated the phylogeography of the spined loach, Cobitis taenia, using 1126 bp of the mitochondrial cytochrome b gene from 174 individuals collected at 47 sites. In total, 51 haplotypes that differed at 49 positions (4.35%) were detected. We deduce that C. taenia survived European glaciations in at least three refugees in the Ponto-Caspian area. Two of these refugees each provided a major lineage that recolonized Europe in separate directions: one westward to England and the other spreading north into Russia before moving west. A third (minor) lineage that contributed little to the recolonization of Europe was also revealed--remaining near its Black Sea refuge. However, more recent history was difficult to resolve with colonization from a more western refugium during the last glacial maximum (LGM) a distinct possibility. Nested clade analysis indicates a pattern of restricted gene flow with isolation by distance at the first two levels and overall. Unlike many other European freshwater fish species, the Danube is not part of the current distribution of C. taenia, nor was it used as either a refuge or a source of colonization of Europe. Low genetic diversity within C. taenia suggests that its colonization of Europe is relatively recent. Demographic analyses revealed a history of recent expansion and isolation by distance.

The origin, current diversity and future conservation of the modern lion (Panthera leo)

Understanding the phylogeographic processes affecting endangered species is crucial both to interpreting their evolutionary history and to the establishment of conservation strategies. Lions provide a key opportunity to explore such processes; however, a lack of genetic diversity and shortage of suitable samples has until now hindered such investigation. We used mitochondrial control region DNA (mtDNA) sequences to investigate the phylogeographic history of modern lions, using samples from across their entire range. We find the sub-Saharan African lions are basal among modern lions, supporting a single African origin model of modern lion evolution, equivalent to the 'recent African origin' model of modern human evolution. We also find the greatest variety of mtDNA haplotypes in the centre of Africa, which may be due to the distribution of physical barriers and continental-scale habitat changes caused by Pleistocene glacial oscillations. Our results suggest that the modern lion may currently consist of three geographic populations on the basis of their recent evolutionary history: North African-Asian, southern African and middle African. Future conservation strategies should take these evolutionary subdivisions into consideration.

An unusual phylogeography in the bushcricket Ephippiger ephippiger from Southern France

Pleistocene glaciations have played a major role in species divergence. The bushcricket Ephippiger ephippiger shows unusual patterns of intraspecific variation in multiple traits across Southern Europe. This is centred in Southern France, and evidence implies that it results from secondary contact after differentiation in Pleistocene refugia. However, the possible time scales involved, locations of the refugia and patterns of expansion remain obscure. This study sequenced the COII (507 BP) and cyt b (428 BP) mitochondrial genes to examine the intraspecific phylogeography of Western European samples of E. ephippiger. A minimum evolution tree revealed little resolution between described subspecies of E. ephippiger. Strikingly, populations from the Pyrenees and Mediterranean coastal region contained a complex genetic structure corresponding to major river valleys, independent of the traditional taxonomy. Samples of the subspecies E. e. vitium formed a distinct clade, perhaps supporting their taxonomic status. However, other forms (cruciger and cunii) were not genetically distinct, which is surprising given differences in their morphology and behaviour. The extent of the genetic divergence between Pyreneen valleys is unexpectedly deep, with average Tamura-Nei distances of around 14% (net distances of 11%) separating the main clades of coding COII sequences. Cyt b showed a similar pattern, but was confounded by some non-coding probable pseudogenes. If a conventional insect molecular clock is applied, these cryptic clades must pre-date the Pleistocene, and hypotheses for their history are discussed. However, mtDNA divergence in Ephippiger is not evolving in a clock-like manner, because a likelihood ratio test rejects clock assumptions for the COII sequences.

Genetic diversity and phylogeography of the Apennine yellow-bellied toad Bombina pachypus, with implications for conservation

Genetic variation was investigated in 17 populations of the Italian endemic Apennine yellow-bellied toad using both mitochondrial (598 bp of the cytochrome b gene) and nuclear (21 allozyme loci) markers. Populations from central Calabria (southern Italy) showed the highest levels of intrapopulation genetic variation, whereas samples located north of this region were nearly lacking in variation. This appears to be a typical pattern of 'southern richness and northern purity', usually attributed to the prolonged population stability within southern refugia coupled with the loss of variation during postglacial northward expansion. However, the overall pattern of genetic variation observed has a strong geographical component, suggesting two Calabrian plains, Catanzaro and Crati-Sibari, as historical barriers to dispersal separating three population groups. These findings cannot be explained by the prolonged stability of southern populations alone, and suggest that the southern richness has been at least in part shaped by allopatric differentiation within the refugial range, followed by intermixing of previously differentiated lineages. From a conservation standpoint, Calabria is the major genetic diversity reservoir for this species, thus deserving particular conservation efforts. Furthermore, although the low intrapopulation genetic variation outside Calabria appears to be of clear historical origin, evidence of a current reduction of gene flow suggests that human disturbance has also played a part, particularly in the anthropogenic impacted Volturno river drainage basin.

Phylogeography of the endangeredCathaya argyrophylla(Pinaceae) inferred from sequence variation of mitochondrial and nuclear DNA

Cathaya argyrophylla is an endangered conifer restricted to subtropical mountains of China. To study phylogeographical pattern and demographic history of C. argyrophylla, species-wide genetic variation was investigated using sequences of maternally inherited mtDNA and biparentally inherited nuclear DNA. Of 15 populations sampled from all four distinct regions, only three mitotypes were detected at two loci, without single region having a mixed composition (G(ST) = 1). Average nucleotide diversity (theta(ws) = 0.0024; pi(s) = 0.0029) across eight nuclear loci is significantly lower than those found for other conifers (theta(ws) = 0.003 approximately 0.015; pi(s) = 0.002 approximately 0.012) based on estimates of multiple loci. Because of its highest diversity among the eight nuclear loci and evolving neutrally, one locus (2009) was further used for phylogeographical studies and eight haplotypes resulting from 12 polymorphic sites were obtained from 98 individuals. All the four distinct regions had at least four haplotypes, with the Dalou region (DL) having the highest diversity and the Bamian region (BM) the lowest, paralleling the result of the eight nuclear loci. An AMOVA revealed significant proportion of diversity attributable to differences among regions (13.4%) and among populations within regions (8.9%). F(ST) analysis also indicated significantly high differentiation among populations (F(ST) = 0.22) and between regions (F(ST) = 0.12-0.38). Non-overlapping distribution of mitotypes and high genetic differentiation among the distinct geographical groups suggest the existence of at least four separate glacial refugia. Based on network and mismatch distribution analyses, we do not find evidence of long distance dispersal and population expansion in C. argyrophylla. Ex situ conservation and artificial crossing are recommended for the management of this endangered species.

Phylogeography and Pleistocene refugia of the adder (Vipera berus) as inferred from mitochondrial DNA sequence data

In order to contribute to the debate about southern glacial refugia used by temperate species and more northern refugia used by boreal or cold-temperate species, we examined the phylogeography of a widespread snake species (Vipera berus) inhabiting Europe up to the Arctic Circle. The analysis of the mitochondrial DNA (mtDNA) sequence variation in 1043 bp of the cytochrome b gene and in 918 bp of the noncoding control region was performed with phylogenetic approaches. Our results suggest that both the duplicated control region and cytochrome b evolve at a similar rate in this species. Phylogenetic analysis showed that V. berus is divided into three major mitochondrial lineages, probably resulting from an Italian, a Balkan and a Northern (from France to Russia) refugial area in Eastern Europe, near the Carpathian Mountains. In addition, the Northern clade presents an important substructure, suggesting two sequential colonization events in Europe. First, the continent was colonized from the three main refugial areas mentioned above during the Lower-Mid Pleistocene. Second, recolonization of most of Europe most likely originated from several refugia located outside of the Mediterranean peninsulas (Carpathian region, east of the Carpathians, France and possibly Hungary) during the Mid-Late Pleistocene, while populations within the Italian and Balkan Peninsulas fluctuated only slightly in distribution range, with larger lowland populations during glacial times and with refugial mountain populations during interglacials, as in the present time. The phylogeographical structure revealed in our study suggests complex recolonization dynamics of the European continent by V. berus, characterized by latitudinal as well as altitudinal range shifts, driven by both climatic changes and competition with related species.

Spatial patterns of genetic diversity across European subspecies of the mountain hare, Lepus timidus L

Fossil evidence shows that populations of species that currently inhabit arctic and boreal regions were not isolated in refugia during glacial periods, but instead maintained populations across large areas of central Europe. These species commonly display little reduction in genetic diversity in northern areas of their range, in contrast to many temperate species. The mountain hare currently inhabits both temperate and arctic-boreal regions. We used nuclear microsatellite and mtDNA sequence data to examine population structure and alternate phylogeographic hypotheses for the mountain hare, that is, temperate type (lower genetic diversity in northern areas) and arctic-boreal type (high northern genetic diversity). Both data sets revealed concordant patterns. Highest allelic richness, expected heterozygosity and mtDNA haplotype diversity were identified in the most northerly subspecies, indicating that this species more closely maps to phylogeographic patterns observed in arctic-boreal rather than temperate species. With regard to population structure, the Alpine and Fennoscandian subspecies were most genetically similar (F(ST) approximately 0.1). These subspecies also clustered together on the mtDNA tree and were assigned with highest likelihood to a common Bayesian cluster. This is consistent with fossil evidence for intermediate populations in the central European plain, persisting well into the postglacial period. In contrast, the geographically close Scottish and Irish populations occupied separate Bayesian clusters, distinct clades on the mtDNA maximum likelihood tree and were genetically divergent from each other (F(ST) > 0.4) indicating the influence of genetic drift, long isolation (possibly dating from the late glacial era) and/or separate postglacial colonisation routes.

E Pluribus Unum: A phylogenetic and phylogeographic reassessment of Laevapex (Pulmonata: Ancylidae), a North American genus of freshwater limpets

The North American freshwater limpet genus Laevapex (Walker, 1903) is a ubiquitous inhabitant of lentic and slow-moving lotic habitats east of the Rocky Mountains, but uncertainty clouds its systematic affinities, the phylogenetic validity of its constituent nominal species, and its degree of genetic connectivity among drainages. We addressed these issues by sampling the genus throughout much of its collective range and constructing representative nuclear and mitochondrial (mt) gene trees, in addition to performing morphometric analyses of shell shape variation. Our results identify neotropical Gundlachia and South American Uncancylus as sister lineages for Laevapex and reveal a pronounced sub-familial dichotomy within the Ancylidae, separating these three New World genera from a Holarctic (Ferrissia (Ancylus, Rhodacmea)) sister clade. Five nominal taxa (L. fuscus, L. diaphanus, L. peninsulae, L. sp., and "F."arkansasensis), indistinguishable in our morphometric analyses, were polyphyletic in the mt gene trees, exhibited modest levels (< 3.9%) of genetic divergence in the primary (103 of 109 individuals) mt clade and, with one minor exception, they appeared fixed for a single nuclear ITS-2 genotype. Although complicated by the presence of rare, highly divergent mt lineages (of either introgressive or persistent ancestral polymorphic origin) in some populations, the molecular data were consistent with a taxonomic conclusion that these five nominal taxa represent a single polymorphic lineage of the type species L. fuscus. AMOVA analyses indicated that 56% of the observed mt variation could be attributed to among population differences, only two of 36 haplotypes were detected in more than one sampling location, and estimates of among-population mt gene flow were generally low at both regional and continental scales. Unrooted network analyses revealed a number of mt tip clades, one restricted to the southwestern part of the range, the remainder having overlapping distributions in eastern North America. All of the eastern tip clades occurred in the Mid-Atlantic region, and these samples displayed by far the highest levels of collective mt diversity. However, directional gene flow estimates indicated that this region has been a recipient (especially from Alabama populations), rather than a source of haplotypic diversity, implying that it likely represents a center of overlap, not a primary ice age refugium, for this limpet species.

Phylogeographic analyses of a circumarctic coastal and a boreal lacustrine mysid crustacean, and evidence of fast postglacial mtDNA rates

Phylogeographic structures of two weakly dispersing Mysis sibling species, one with a circumarctic coastal, the other with a boreal lacustrine-Baltic distribution, were studied from mitochondrial cytochrome c oxidase subunit I (COI) gene sequences. Mysis segerstralei showed high overall diversity and little phylogeographic structure across the Arctic, indicating late-glacial dispersal among coastal and lake populations from Alaska, Siberia and the north of Europe. A strongly divergent refugial lineage was however identified in Beringia. The boreal 'glacial relict'Mysis salemaai in turn displayed clear structuring among postglacially isolated Scandinavian lake populations. The inferred pattern of intralake mitochondrial DNA (mtDNA) monophyly in Scandinavia suggested relatively small population sizes and a remarkably fast postglacial mtDNA divergence rate (0.27% per 10 000 years). Nevertheless, the broader phylogeographic pattern did not support distinct eastern and western glacial refugia in Northern Europe, unlike in some other aquatic taxa. In all, the two species comprised three equidistant mitochondrial lineages (approximately 2% divergence), corresponding to M. salemaai, to the bulk of M. segerstralei, and to the Beringian M. segerstralei lineage. The lack of reciprocal monophyly of the two species in respect to their mitochondrial genealogy could indicate postspeciation mitochondrial introgression, also exemplified by an evidently more recent capture of M. segerstralei mitochondria in a Karelian population of M. salemaai. Overall, the data suggest that the continental boreal M. salemaai has a relatively recent ancestry in arctic coastal waters, whereas two other boreal 'glacial relict'Mysis sibling species in Europe (Mysis relicta) and North America (Mysis diluviana) have colonized inland waters much earlier (approximately 8% COI divergence).

Divergence and diversity: lessons from an arctic-alpine distribution (Pardosa saltuaria group, Lycosidae)

The relationship of interpopulation genetic divergence and within-population diversity has been studied for many temperate species in Europe, but not for the cold-adapted fauna. Here we present the first European-wide phylogeographical study of an arctic-alpine distribution in invertebrates, focusing on wolf spiders of the Pardosa saltuaria group. One hundred twenty-seven (127) specimens from 14 populations were examined. Within Europe, these populations were distributed among six high mountain ranges and Scandinavia. We sequenced the whole 921 base pair mitochondrial (mt) ND1 gene. The resulting 55 unique haplotypes form three monophyletic phylogroups of deep divergence: a Pyrenean, a Balkan and a 'northern' clade. Genetic distances (3.6-4.0%) between the major clades indicate that the arctic-alpine range disjunction was initiated by vicariance events, which precede the four major Alpine glaciations. However, low divergence and incomplete lineage sorting within the 'northern clade' suggest a late Pleistocene separation of the Alpine, Scandinavian, Carpathian and Sudetian populations. Thus, we provide evidence for a multiglacial origin of arctic-alpine distributions in Europe, i.e. the current disjunction results from range fragmentation in several glacial cycles. The pattern of genetic diversity within populations seems predominantly determined by historical factors, but is modified by contemporary aspects. Overall, diversity and divergence are negatively correlated. We suggest that low diversity values might result from (i) ancient bottlenecking during warm interglacial periods, as seen in the Pyrenees and Balkans; (ii) recent bottlenecking in small modern areas, as seen in the Giant Mountains and Bohemian Forest; and (iii) dispersal bottlenecking in northern Scandinavia.

Pan-European phylogeography of the aquatic snail Theodoxus fluviatilis (Gastropoda: Neritidae)

Investigating the geographical distribution of genetic lineages within species is critical to our understanding of how species evolve. As many species inhabit large and complex ranges, it is important that phylogeographical research take into account the entire range of widespread species to clarify how myriad extrinsic variables have affected their evolutionary history. Using phylogenetic, nested clade, and mismatch distribution analyses on a portion of the mitochondrial COI gene, I demonstrate that the wide-ranging freshwater snail Theodoxus fluviatilis possesses in parallel many of the phylogeographical patterns seen in less widespread freshwater species of Europe. Fragmentary forces play a major part in structuring the range of this species, with 12 of 14 geographically structured nested clades displaying a distribution consistent with fragmentation or restricted dispersal. Certain regions of southern Europe harbour the majority of genetic diversity (total haplotype diversity, H = 0.87), particularly Italy (H = 0.87) and areas surrounding the Black Sea (H = 0.81). Post-Pleistocene range expansion is pronounced, with the majority of northern European populations (95% of sample sites) having arisen from northern Italian individuals that initially colonized northern Germany. Additionally, two highly divergent haplotype lineages present in northern Germany imply that there were at least two postglacial recolonization routes. Estuaries may also provide a means of dispersal given that no genetic differentiation was found between estuarine populations and neighbouring freshwater populations. Taken together, these data reveal a species with a complex genetic history resulting from the fragmentary effects of European geology as well as continuous and discrete range expansion related to their aquatic biology.

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.

Identifying the genome of wood barley Hordelymus europaeus (Poaceae: Triticeae)

Wood barley, Hordelymus europaeus, was compared with other Triticeae species by Southern and fluorescence in situ hybridisation using total genomic DNA and repetitive sequences as probes. On Southern blots, the total genomic probe from H. europaeus hybridised strongly to DNA of its own species and to Leymus and Psathyrostachys, indicating the presence of Ns genome in H. europaeus. Furthermore, the total genomic probe from P. fragilis hybridised to DNA of H. europaeus as much as to all of the Psathyrostachys and Leymus species examined. Ns genome-specific DNA sequences isolated from L. mollis (pLmIs1, pLmIs44 and pLmIs53) hybridised essentially to H. europaeus and all of the species of Leymus and Psathyrostachys. Chromosomal localization of these clones on H. europaeus confirmed the presence of Ns genome-specific DNA on all chromosomes, indiscriminately. Under moderate hybridisation stringency the Ns genome-specific probes, together with repetitive sequences pTa71 and pAesKB7, produced species-specific RFLP banding profiles on Southern blots. A phenetic tree based on these profiles revealed a distinct Ns species cluster within the Triticeae, represented by Leymus and Psathyrostachys species. Hordelymus europaeus belonged to this Ns cluster. Chromosomal mapping of the 18S-25S and the 5S ribosomal genes, together with the repetitive sequence pLrTaiI, corroborated that H. europaeus was most probably related to Leymus, especially the European/Eurasian members of sect. Leymus. In an attempt to identify the genome of H. europaeus, different approaches were employed; the results clearly showed that wood barley had the Ns basic genome and nothing else.

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.

The recent history and population structure of five Mandarina snail species from subtropical Ogasawara (Bonin Islands, Japan)

The effect of Pleistocene climate change on the organisms of tropical and subtropical regions is rather poorly understood. We therefore studied the land snail genus Mandarina (Bradybaenidae) of oceanic Ogasawara (Bonin Islands, Japan), with the aim of using population genetic data to understand their recent history. Our analysis of a mitochondrial 16S ribosomal RNA region from more than 600 snails in five ground-living species suggests that populations on the small islands of Mukoujima, Anejima, Imotojima and Meijima, as well as on the low-lying southern and central parts of Hahajima, have probably undergone recent bottlenecks followed by subsequent expansions. Except between the main island of Hahajima and Mukouijima, there is almost no evidence for gene flow among islands even though the islands were connected repeatedly by land bridges through the Pleistocene. Within islands the population structure is severe, suggestive of a long-term, low level of gene flow (F(ST) is frequently greater than 0.5 among geographically close populations). Finally, there is a marked genetic patchiness, meaning that genetically close populations are sometimes separated by genetically distant populations. These patterns could be a consequence of expansion from bottlenecks, low active dispersal and founder effects caused by rare long-distance migrants. Unfortunately, the exact nature of the refugia and bottlenecks remains unknown because the palaeoclimate of this region is poorly understood. Dating the population size changes is also challenging because the molecular clock is uncertain. We suggest, however, that arid conditions or deforestation induced by decreased atmospheric CO(2) may have been the main factor in determining population size.

Mitochondrial phylogeography and genetic diversity of Tibetan gazelle (Procapra picticaudata): implications for conservation

The Tibetan gazelle (Procapra picticaudata) is a threatened species and distributed on the Qinghai-Tibet Plateau of China (Qinghai Province, Tibet Autonomous Region and the adjacent Gansu Province, Sichuan Province, and Xinjiang Uigur Autonomous Region). Small peripheral populations of Tibetan gazelle were once found in northern Sikkim and Ladakh, but now these are close to extinction. To describe the evolutionary history and to assess the genetic diversity within this monotypic species and population structure among different geographic locations in China, we sequenced mitochondrial DNA from the control region (CR) and cytochrome (cyt) b gene for 46 individuals from 12 geographic localities in Qinghai, Tibet, Xinjiang, Gansu, and Sichuan. A total of 25 CR haplotypes and 16 cyt b haplotypes were identified from these gazelle samples. CR haplotype diversity (0.98+/-0.01) and nucleotide diversity (0.08+/-0.009) were both high. Phylogenetic trees indicate that the Tibetan gazelle in China can be divided into three main clades: Tibet, Sichuan (SCH) and Qinghai-Arjin Shan-Kekexili (QH-ARJ-KKXL). Analysis of molecular variance (AMOVA) and network analysis consistently support this geographic structure in both datasets. Significant differentiation between populations argues for the presence of management units (MUs). Such differentiation may reflect a geographic separation resulting from the uplift of the Qinghai-Tibet Plateau during the Late Pliocene and Pleistocene. Mismatch distribution analysis implies that Tibetan gazelle has undergone complex population changes. We suggest that the present population structure has resulted from habitat fragmentation during the recent glacial period on the Qinghai-Tibet Plateau and population expansion from glacial refugia after the glacial period. It is likely that the present populations of Tibetan gazelle exhibit a pattern reminiscent of several bottlenecks and expansions in the recent past.

Recency, range expansion, and unsorted lineages: implications for interpreting neutral genetic variation in the sharp-tailed grouse (Tympanuchus phasianellus)

Both current and historical patterns of variation are relevant to understanding and managing ecological diversity. Recently derived species present a challenge to the reconstruction of historical patterns because neutral molecular data for these taxa are more likely to exhibit effects of recent and ongoing demographic processes. We studied geographical patterns of neutral molecular variation in a species thought to be of relatively recent origin, Tympanuchus phasianellus (sharp-tailed grouse), using mitochondrial control region sequences (CR-I), amplified fragment length polymorphisms (AFLP), and microsatellites. For historical context, we also analysed CR-I in all species of Tympanuchus. Within T. phasianellus, we found evidence for restricted gene flow between eastern and western portions of the species range, generally corresponding with the range boundary of T. p. columbianus and T. p. jamesi. The mismatch distribution and molecular clock estimates from the CR-I data suggested that all Tympanuchus underwent a range expansion prior to sorting of mitotypes among the species, and that sorting may have been delayed as a result of mutation-drift disequilibrium. This study illustrates the challenge of using genetic data to detect historical divergence in groups that are of relatively recent origin, or that have a history dominated by nonequilibrium conditions. We suggest that in such cases, morphological, ecological, and behavioural data may be particularly important adjuncts to molecular data for the recognition of historically or adaptively divergent groups.

The determinant role of temporary proglacial drainages on the genetic structure of fishes

Phylogeographic studies have shed light on Pleistocene glaciations as a key factor in shaping present-day genetic structure of many organisms. In formerly glaciated regions, the combined action of several factors such as refuges origin, physiological capacities and demographic parameters have contributed importantly to this process but specifically for each species. Therefore, a fine-scale genetic structure is not expected to be similar for different species, unless it has been modulated by the action of a strong environmental pressure. The aim of this study is to investigate the effects of postglacial environment on the genetic structure of fishes. To achieve this objective, three fish species (northern pike, lake whitefish and yellow perch) commonly found in sympatry in Laurentian Shield lakes but displaying different ecological and physiological characteristics were analysed. The comparison of these unrelated species was performed to identify the factors determining the organization of their genetic structure. Populations of all species mostly originated from the Mississippian refuge. Low genetic differentiation was observed among populations but significant structures were detected for the three species. Despite marked differences among species, these structures presented common characteristics: a lack of congruence with drainage and a longitudinal organization. This suggested that the dispersion of species occurred independently, leading to a species-specific structure. However, the settling of populations appeared to be mediated by a dynamic system of proglacial meltwater streams associated to the glacial Lake Ojibway-Barlow, providing such similarities among species.

The evolutionary history of Sium sensu lato (Apiaceae): dispersal, vicariance, and domestication as inferred from ITS rDNA phylogeny

The biogeographic history of Sium sensu lato (s.l.) (including Sium sensu stricto [s.s.], Afrocarum, and Berula) was inferred using a nuclear rDNA ITS sequence phylogeny and dispersal-vicariance analysis. One hundred accessions were analyzed, including 86 specimens of Sium s.l., representing all taxa throughout their geographical range. The phylogenetic analyses showed that Afrocarum and the African and St. Helena species of Sium were nested within Berula forming the Berula s.l. clade. The remaining species of Sium constituted a weakly supported Sium s.s. clade. The cultivated S. sisarum is closely related to the western Asian S. sisaroideum. The divergence of Sium s.s. from Berula s.l., estimated under the assumption of a molecular clock and using a substitution rate of 0.61% per million years, took place around 10.4 million years ago. Dispersal-vicariance analysis suggested that the biogeographic pattern of this group resulted from simple vicariance and eventual migration to adjacent areas rather than from long-distance dispersal. The ancestral area of Sium s.l. was probably western Palearctic. The dispersal of Berula to North America occurred ca. 3 million years ago, possibly by island hopping. The migration of S. suave to North America took place ca 0.6 million years ago through Beringia.

Rare arctic-alpine plants of the European Alps have different immigration histories: the snow bed species Minuartia biflora and Ranunculus pygmaeus

Minuartia biflora and Ranunculus pygmaeus are circumarctic plants with a few isolated occurrences in the European Alps. We analysed amplified fragment length polymorphism (AFLP) and chloroplast DNA sequence data to unravel the history of their immigration into the Alps and to provide data on their circumpolar phylogeography. In spite of the similar ecological requirements of the two species, they exhibit strikingly different immigration histories into the Alps. In M. biflora, the Alpine populations are most probably derived from source populations located between the Alpine and Scandinavian ice sheets, in accordance with the traditional biogeographic hypothesis. In contrast, the Alpine populations of R. pygmaeus cluster with those from the Tatra Mountains and the Taymyr region in northern Siberia, indicating that the distant Taymyr area served as source for the Alpine populations. Both species showed different levels of genetic diversity in formerly glaciated areas. In contrast to the considerable AFLP diversity observed in M. biflora, R. pygmaeus was virtually nonvariable over vast areas, with a single phenotype dominating all over the Alps and another, distantly related one dominating the North Atlantic area from Greenland over Svalbard to Scandinavia. The same pattern was observed in chloroplast DNA sequence data. Thus, postglacial colonization of R. pygmaeus was accompanied by extreme founder events.

Inferring phylogeny despite incomplete lineage sorting

It is now well known that incomplete lineage sorting can cause serious difficulties for phylogenetic inference, but little attention has been paid to methods that attempt to overcome these difficulties by explicitly considering the processes that produce them. Here we explore approaches to phylogenetic inference designed to consider retention and sorting of ancestral polymorphism. We examine how the reconstructability of a species (or population) phylogeny is affected by (a) the number of loci used to estimate the phylogeny and (b) the number of individuals sampled per species. Even in difficult cases with considerable incomplete lineage sorting (times between divergences less than 1 N(e) generations), we found the reconstructed species trees matched the "true" species trees in at least three out of five partitions, as long as a reasonable number of individuals per species were sampled. We also studied the tradeoff between sampling more loci versus more individuals. Although increasing the number of loci gives more accurate trees for a given sampling effort with deeper species trees (e.g., total depth of 10 N(e) generations), sampling more individuals often gives better results than sampling more loci with shallower species trees (e.g., depth = 1 N(e)). Taken together, these results demonstrate that gene sequences retain enough signal to achieve an accurate estimate of phylogeny despite widespread incomplete lineage sorting. Continued improvement in our methods to reconstruct phylogeny near the species level will require a shift to a compound model that considers not only nucleotide or character state substitutions, but also the population genetics processes of lineage sorting. [Coalescence; divergence; population; speciation.].