Northern and Southern expansions of Atlantic brown trout (Salmo trutta) populations during the Pleistocene

The phylogeography of Atlantic brown trout ( Salmo trutta ) was analysed using mitochondrial DNA control region complete sequences of 774 individuals from 57 locations. Additionally, the available haplotype information from 100 published populations was incorporated in the analysis. Combined information from nested clade analysis, haplotype trees, mismatch distributions, and coalescent simulations was used to characterize population groups in the Atlantic basin. A major clade involved haplotypes assigned to the Atlantic (AT) lineage, but another major clade should be considered as a distinct endemic lineage restricted to the Iberian Peninsula. The phylogeography of the Atlantic populations showed the mixed distribution of several Atlantic clades in glaciated areas of Northern Europe, whereas diverged haplotypes dominated the coastal Iberian rivers. Populations inhabiting the Atlantic rivers of southern France apparently contributed to postglacial colonization of northern basins, but also comprised the source of southern expansions during the Pleistocene.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 904–917.     

Phylogeographic analysis reveals a deep lineage split within North Atlantic Littorina saxatilis

Phylogeographic studies provide critical insight into the evolutionary histories of model organisms; yet, to date, range-wide data are lacking for the rough periwinkle Littorina saxatilis, a classic example of marine sympatric speciation. Here, we use mitochondrial DNA (mtDNA) sequence data to demonstrate that L. saxatilis is not monophyletic for this marker, but is composed of two distinct mtDNA lineages (I and II) that are shared with sister species Littorina arcana and Littorina compressa. Bayesian coalescent dating and phylogeographic patterns indicate that both L. saxatilis lineages originated in the eastern North Atlantic, around the British Isles, at approximately 0.64 Ma. Both lineages are now distributed broadly across the eastern, central and western North Atlantic, and show strong phylogeographic structure among regions. The Iberian Peninsula is genetically distinct, suggesting prolonged isolation from northeastern North Atlantic populations. Western North Atlantic populations of L. saxatilis lineages I and II predate the last glacial maximum and have been isolated from eastern North Atlantic populations since that time. This identification of two distinct, broadly distributed mtDNA lineages further complicates observed patterns of repeated incipient ecological speciation in L. saxatilis, because the sympatric origins of distinct ecotype pairs on eastern North Atlantic shores may be confounded by admixture of divergent lineages.     

Refugia within refugia – patterns in endemism and genetic divergence are linked to Late Quaternary climate stability in the Iberian Peninsula

In Europe, southern peninsulas served as major refugia during Pleistocene cold periods. However, growing evidence has revealed complex patterns of glacial survival within these southern regions, with multiple glacial refugia within each larger refugial area. We investigated the extent to which patterns of endemism and phylogeographic are concordant across animal species in the Iberian Peninsula, one of the most important unglaciated areas in Europe during the Pleistocene, can be explained in terms of climatic stability. We found that historical climatic stability (notably climate velocity measures integrating macroclimatic shifts with local spatial topoclimate gradients) was often among the most important predictors of endemic species richness for different taxonomic groups using models that also incorporated measures of modern climate. Furthermore, for some taxonomic groups, climatic stability was also correlated with patterns of spatial concordance in interpopulation genetic divergence across multiple taxa, and private haplotypes were more frequently found in relatively stable areas. Overall, our results suggest that both endemism patterns and cross‐taxa concordant phylogeographic patterns across the Iberian Peninsula to some extent are linked to spatial variation in Late Quaternary climate stability, in agreement with the proposed ‘refugia‐within‐refugia’ scenario. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 13–28.     

Variation in mitochondrial DNA reveals multiple distant glacial refugia in black spruce (Picea mariana), a transcontinental North American conifer

Range-wide genetic variation of black spruce (Picea mariana) was studied using polymerase chain reaction-random fragment length polymorphism markers of the mitochondrial genome. Four polymorphic mitochondrial DNA (mtDNA) loci were surveyed and two or three alleles were detected at each locus, resulting in 10 multilocus mtDNA types or mitotypes. A significant subdivision of population genetic diversity was detected (GST = 0.671; NST = 0.726), suggesting low levels of gene flow among populations. The distribution of mitotypes was not random (NST > GST; P < 0.05) and revealed four partially overlapping zones, presumably representative of different glacial populations. Comparison of the genetic structure derived from mtDNA markers and the colonization paths previously deduced from the fossil and pollen records allow us to infer at least three southern and one northeastern glacial populations for black spruce. The patterns revealed in this study suggest that black spruce shares its biogeographical history with other forest-associated North American species.     

A new cytochrome b phylogroup of the common vole (Microtus arvalis) endemic to the Balkans and its implications for the evolutionary history of the species

The phylogeographic architecture of the common vole, Microtus arvalis, has been well‐studied using mitochondrial DNA and used to test hypotheses relating to glacial refugia. The distribution of the five described cytochrome b (cyt b) lineages in Europe west of Russia has been interpreted as a consequence of postglacial expansion from both southern and central European refugia. A recently proposed competing model suggests that the ‘cradle’ of the M. arvalis lineages is in western central Europe from where they dispersed in different directions after the Last Glacial Maximum. In the present study, we report a new cyt b lineage of the common vole from the Balkans that is not closely related to any other lineage and whose presence might help resolve these issues of glacial refugia. The Balkan phylogroup occurs along the southern distributional border of M. arvalis in central and eastern Bosnia and Herzegovina, Montenegro, and eastern Serbia. Further north and west in Slovenia, Bosnia and Herzegovina, and Serbia, common voles belong to the previously‐described Eastern lineage, whereas both lineages are sympatric in one site in Bosnia and Herzegovina. The Balkan phylogroup most reasonably occupied a glacial refugium already known for various Balkan endemic species, in contrast to the recently proposed model. South‐east Europe is an absolutely crucial area for understanding the postglacial colonization history of small mammals in Europe and the present study adds to the very few previous detailed phylogeographic studies of this region. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 788–796.     

Phylogenetic structure of Neotropical annual fish of the genus Cynopoecilus (Cyprinodontiformes: Rivulidae), with an assessment of taxonomic implications

The definition of species boundaries constitutes an important challenge in biodiversity studies. Cynopoecilus Regan, 1912 encompasses several endangered species of annual fish, occurring in temporary ponds in a restricted area of Southern Brazil and Uruguay. Divergences about the taxonomic status of Cynopoecilus species highlight the importance of species delimitation studies. Therefore, we address here the phylogenetic structure of Cynopoecilus, while assessing its taxonomic implications. For this, fragments of the mitochondrial COI and nuclear RAG1 genes were characterized and analyzed for a set of 275 and 280 specimens, respectively. DNA barcoding and phylogenetic analyses detected subdivision of these specimens in 8–10 clusters, which comprise the six previously described species, and suggest one invalid taxon and at least 3–5 putative new species. The phylogenetic structure also suggests that the Jacuí River and the Patos Lagoon historically acted as effective barriers to gene flow between populations, although some isolated dispersal events across these water bodies could be evidenced, especially for C. melanotaenia Regan, 1912. In general, the results highlight the need of independent conservation strategies within the distribution area of each of the endemic allopatric killifish clusters, while questioning several taxonomic boundaries and distribution data.     

Phylogeography of the montane caddisfly Drusus 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.     

Multiple introgression events and range shifts in Schizocodon (Diapensiaceae) during the Pleistocene

Range shifts during the Pleistocene shaped the unique phylogeographical structures of numerous species. Accompanying species migration, sister taxa may have experienced multiple introgression events. Here, we report the signature of introgression events in multiple areas in Schizocodon, herbs endemic to Japan, using amplified fragment length polymorphism (AFLP) fingerprinting and plastid DNA haplotyping in 48 populations. Although the present distributions of S. soldanelloides and S. ilicifolius are mainly allopatric, the species share plastid DNA haplotypes in each region (north‐eastern, north‐central, south‐central and south‐western Japan); in contrast, the specific groups were highly supported by AFLP analyses. These results support the occurrence of multiple introgression events in Schizocodon. Notably, the disjunct plastid haplotypes found only in S. ilicifolius var. intercedens suggest complete plastid DNA replacement at local areas from S. soldanelloides into S. ilicifolius var. ilicifolius. Furthermore, we found that S. soldanelloides experienced range contraction and expansion during glacial and interglacial cycles based on mismatch distribution analysis and ecological niche modelling. Based on several pieces of evidence, our study supports the idea that historical range shifts associated with Pleistocene climatic oscillations favoured multiple and regional introgression events in Schizocodon. © 2013 The Linnean Society of London, Botanical Journal of the Linnean Society, 2013, 173 , 46–63.     

Beyond the Mediterranean peninsulas: evidence of central European glacial refugia for a temperate forest mammal species, the bank vole (Clethrionomys glareolus)

This study details the phylogeographic pattern of the bank vole, Clethrionomys glareolus, a European rodent species strongly associated with forest habitat. We used sequences of 1011 base pairs of the mitochondrial DNA cytochrome b gene from 207 bank voles collected in 62 localities spread throughout its distribution area. Our results reveal the presence of three Mediterranean (Spanish, Italian and Balkan) and three continental (western, eastern and 'Ural') phylogroups. The endemic Mediterranean phylogroups did not contribute to the post-glacial recolonization of much of the Palaearctic range of species. Instead, the major part of this region was apparently recolonized by bank voles that survived in glacial refugia in central Europe. Moreover, our phylogeographic analyses also reveal differentiated populations of bank voles in the Ural mountains and elsewhere, which carry the mitochondrial DNA of another related vole species, the ruddy vole (Clethrionomys rutilus). In conclusion, this study demonstrates a complex phylogeographic history for a forest species in Europe which is sufficiently adaptable that, facing climate change, survives in relict southern and northern habitats. The high level of genetic diversity characterizing vole populations from parts of central Europe also highlights the importance of such regions as a source of intraspecific genetic biodiversity.     

Phylogeographic structure and past history of the circum-Mediterranean species Tomicus destruens Woll. (Coleoptera: Scolytinae)

Phylogeographic studies are often focused on temperate European species with relict footholds in the Mediterranean region. Past climatic oscillations usually induced range contractions and expansions from refugial areas located in southern Europe, and spatial distribution of genetic diversity show that northward expansions were usually pioneer-like. Actually, few studies have focused on circum-Mediterranean species, which probably were not influenced in the same way by climatic oscillations. We present the phylogeography of the bark beetle Tomicus destruens, which is restricted to the whole Mediterranean basin and the Atlantic coasts of North Africa and Portugal. We systematically sequenced 617 bp of the mitochondrial genes COI and COII for 42 populations (N = 219). Analysis revealed 53 haplotypes geographically structured in two clades, namely eastern and western clades, that diverged during the Pleistocene. A contact zone was identified along the Adriatic coast of Italy. Interestingly, we found contrasting levels of genetic structure within each clade. The eastern group was characterized by a significant phylogeographic pattern and low levels of gene flow, whereas the western group barely showed a spatial structure in haplotype distribution. Moreover, the main pine hosts were different between groups, with the Aleppo-brutia complex in the east and the maritime pine in the west. Potential roles of host species, climatic parameters and geographical barriers are discussed and the phylogeographic patterns are compared to classical models of postglacial recolonization in Europe.     

Phylogeography of the rare Balkan endemic Martino's vole, Dinaromys bogdanovi, reveals strong differentiation within the western Balkan Peninsula

The spatial genetic structure of Martino's vole, a rare palaeoendemic species of the western Balkans, was investigated using DNA isolated from archived museum samples. The study was based on partial sequencing (555 bp) of the mitochondrial cytochrome b gene for 63 specimens from 20 different localities throughout the species' range. Three highly divergent allopatric phylogenetic lineages (Northwestern, Central and Southeastern) were recognized among 47 haplotypes, suggesting three independent glacial differentiation centres within the western Balkans. The Northwestern lineage, which showed the highest divergence from all other samples (mean sequence divergence of 6.64% +/- 1.10), comprised samples collected from northwest of the Neretva River in Croatia, western Bosnia and Herzegovina. The Central and Southeastern lineages were separated by an average sequence divergence of 2.95% +/- 0.66 and were geographically divided by the Drim River (the Kosovo Basin in Serbia). Overall, haplotype diversity decreased from the Northwestern lineage to the Central and subsequently the Southeastern lineage, in a geographical pattern consistent with a stepping stone colonization. The observed distribution indicates a gradual southerly expansion with subsequent allopatry across the Neretva River and Drim River approximately 1 and 0.3 million years ago, respectively. Such a scenario is concordant with palaeontological evidence. Several highly divergent sublineages within the Northwestern and Central lineages showed no significant geographical structuring, suggesting secondary contact of allopatrically evolved lineages. We hypothesize that the topographical complexity of the Balkans promoted allopatry and isolation on a small geographical scale during interglacial periods, with secondary contact during glacial maxima. Furthermore, the three main lineages should be regarded as evolutionary significant units with important implications for conservation. Ecological data show that the Northwestern lineage in particular fulfils all criteria for a highly endangered, evolutionarily significant unit.     

Mitochondrial differentiation in a polymorphic land snail: evidence for Pleistocene survival within the boundaries of permafrost

The genetic differentiation of populations having colonized formerly unsuitable habitats after the Pleistocene glaciations depends to a great extent on the speed of expansion. Slow dispersers maintain their refugial diversity whereas fast dispersal leads to a reduction of diversity in the newly colonized areas. During the Pleistocene, almost the entire current range of the land snail Arianta arbustorum has repeatedly been covered with ice or been subjected to permafrost. Owing to the low potential for dispersal of land snails, slow (re)colonization of the wide range from southern refugia can be excluded. Alternatively, fast, passive dispersal from southern refugia or survival in and expansion from multiple refugia within the area subjected to permafrost may account for the current distribution. To distinguish between these scenarios we reconstructed a phylogeography based on the sequences of a fragment of the cytochrome oxidase I from 133 individuals collected at 45 localities and analysed the molecular variance. Seventy-five haplotypes were found that diverged on average at 7.52% of positions. This high degree of diversity suggests that A. arbustorum is an old species in which the population structure, isolation and the hermaphroditic nature have reduced the probability of lineage extinction. The genetic structure was highly significant with the highest variance partition found among regions. Geographic distance and mitochondrial differentiation were not congruent. Lineages had overlapping ranges. The clear genetic differentiation and the patchy pattern of haplotype distribution suggest that colonization of formerly unsuitable habitats was mainly achieved from multiple populations from within the permafrost area.     

Isolation and characterization of microsatellite markers in the caddisfly Drusus discolor (Trichoptera: Limnephilidae)

We describe the development of and amplification conditions for microsatellite primers isolated from the caddisfly Drusus discolor. Eight polymorphic microsatellite loci were developed and screened for variability using 37 individuals from two populations from central Europe. The primers yielded an average of 8.6 alleles per loci. No linkage disequilibrium between loci was detected, while three loci showed deviations from Hardy–Weinberg equilibrium in one of the two tested populations.     

Phylogeography and historical demography of the Italian treefrog, Hyla intermedia, reveals multiple refugia, population expansions and secondary contacts within peninsular Italy

We investigated the geographical patterns of genetic diversity in the Italian treefrog through sequence analysis of a mitochondrial cytochrome b gene fragment. Three main mitochondrial lineages were identified, distributed in northern, central and southern Italy, respectively. Their divergence appears indicative of a split time largely predating Late Pleistocene climatic oscillations, and syntopy between them was only observed in the geographically intermediate populations. The historical demographic reconstructions suggest that in both northern and central Italy, an expansion occurred during the last major glacial phase, when a vast widening of the lowland habitats followed the glaciation-induced fall of the sea level. Instead, in southern Italy an expansion event likely followed the end of the last glaciation, although the inference of expansion appears less reliable for the southern clade than for the others. Within this geographical area, a sharp phylogeographic discontinuity separated peninsular from Sicilian populations, and the overall pattern of diversity suggests that the latter derived from a recent colonization of the island, probably through a Late Pleistocene land bridge. Phylogenetic, phylogeographic and historical demographic analyses thus concur in delineating a scenario of multiple refugia, with four groups of populations which survived the last glacial-interglacial cycles in at least three distinct refugia arranged along peninsular Italy, and have recently come into contact following range expansions. Therefore, these results support the hypothesis that a plethora of microevolutionary processes, rather than the prolonged stability of populations, were mainly responsible for shaping the patterns of diversity within this major biodiversity hotspot.     

Phylogeography of red deer (Cervus elaphus) in Europe

Aim To investigate the phylogeographical patterns of red deer (Cervus elaphus) in Europe, and to disentangle the influence of ancient (e.g. Pleistocene ice ages) from more recent processes (e.g. human translocations). Location Europe. Methods In this study we provide by far the most extensive analysis of genetic structure in European red deer, based on analyses of variation at two mitochondrial markers (cyt b and D‐loop) in a large number of individuals from 39 locations. Relationships of mitochondrial DNA haplotypes were determined using minimum spanning networks and phylogenetic analyses. Population structure was examined by analyses of molecular variance. Historical processes shaping the present patterns were inferred from nested clade analysis and nucleotide diversity statistics. Results Within Europe, we detected three deeply divergent mitochondrial DNA lineages. The three lineages displayed a phylogeographical pattern dividing individuals into western European, eastern European and Mediterranean (Sardinia, Spain and Africa) groups, suggesting contraction into three separate refugia during the last glaciation. Few haplotypes were shared among these three groups, a finding also confirmed by F_ST values. Calculations of divergence times suggest that the groups probably split during the Pleistocene. Main conclusions The observed pattern is interpreted to result from isolation in different refugia during the last glaciation. The western and eastern European lineages could be linked to an Iberian and Balkan refugium, respectively. The third lineage might originate from a Sardinian or African refugium. We link local phylogeographical patterns observed in Europe to the post‐glacial recolonization process, shaped by the geographical localization of refugia and barriers to gene flow. Regardless of the importance of red deer as a game species and the tradition of translocating red deer in Europe, we detected few individuals that did not match the trichotomous pattern, suggesting that translocations have occurred mainly at smaller spatial scales.     

A complex scenario of glacial survival in Mediterranean and continental refugia of a temperate continental vole species (Microtus arvalis) in Europe

The role of glacial refugia in shaping contemporary species distribution is a long-standing question in phylogeography and evolutionary ecology. Recent studies are questioning previous paradigms on glacial refugia and postglacial recolonization pathways in Europe, and more flexible phylogeographic scenarios have been proposed. We used the widespread common vole Microtus arvalis as a model to investigate the origin, locations of glacial refugia, and dispersal pathways, in the group of “Continental” species in Europe. We used a Bayesian spatiotemporal diffusion analysis (relaxed random walk model) of cytochrome b sequences across the species range, including newly collected individuals from 10 Iberian localities and published sequences from 68 localities across 22 European countries. Our data suggest that the species originated in Central Europe, and we revealed the location of multiple refugia (in both southern peninsulas and continental regions) for this continental model species. Our results confirm the monophyly of Iberian voles and the pre-LGM divergence between Iberian and European voles. We found evidence of restricted postglacial dispersal from refugia in Mediterranean peninsulas. We inferred a complex evolutionary and demographic history of M. arvalis in Europe over the last 50,000 years that does not adequately fit previous glacial refugial scenarios. The phylogeography of M. arvalis provides a paradigm of ice-age survival of a temperate continental species in western and eastern Mediterranean peninsulas (sources of endemism) and multiple continental regions (sources of postglacial spread). Our findings also provide support for a major role of large European river systems in shaping geographic boundaries of M. arvalis in Europe.     

Contrasting phylogeographic patterns in Sphagnum fimbriatum and Sphagnum squarrosum (Bryophyta, Sphagnopsida) in Europe

The chloroplast phylogeography of two peat mosses (Sphagnum fimbriatum and Sphagnum squarrosum) with similar distributions but different life history characteristics was investigated in Europe. Our main aim was to test whether similar distributions reflect similar phylogeographic and phylodemographic processes. Accessions covering the European distributions of the species were collected and approx. 2000 bp of the chloroplast genome of each species was sequenced. Maximum parsimony, statistical parsimony and phylodemographic analyses were used to address the question of whether these species with similar distributions show evidence of similar phylogeographic and phylodemographic processes. The chloroplast haplotypes of the currently spreading species S. fimbriatum showed strong geographic structure, whereas those of S. squarrosum, which has stable historical population sizes, showed only very weak geographic affinity and were widely distributed. We hypothesize that S. fimbriatum survived the last glaciations along the Atlantic coast of Europe, whereas S. squarrosum had numerous, scattered refugia in Europe. The dominance of one haplotype of S. fimbriatum across almost all of Europe suggests rapid colonization after the last glacial maximum. We hypothesize that high colonizing ability is an inherent characteristic of the species and its recent expansion in Europe is a response to climate change.     

Multilocus tests of Pleistocene refugia and ancient divergence in a pair of Atlantic Forest antbirds (Myrmeciza)

The Atlantic Forest (AF) harbours one of the most diverse vertebrate faunas of the world, including 199 endemic species of birds. Understanding the evolutionary processes behind such diversity has become the focus of many recent, primarily single locus, phylogeographic studies. These studies suggest that isolation in forest refugia may have been a major mechanism promoting diversification, although there is also support for a role of riverine and geotectonic barriers, two sets of hypotheses that can best be tested with multilocus data. Here we combined multilocus data (one mtDNA marker and eight anonymous nuclear loci) from two species of parapatric antbirds, Myrmeciza loricata and M. squamosa, and Approximate Bayesian Computation to determine whether isolation in refugia explains current patterns of genetic variation and their status as independent evolutionary units. Patterns of population structure, differences in intraspecific levels of divergence and coalescent estimates of historical demography fit the predictions of a recently proposed model of refuge isolation in which climatic stability in the northern AF sustains higher diversity and demographic stability than in the southern AF. However, a pre‐Pleistocene divergence associated with their abutting range limits in a region of past tectonic activity also suggests a role for rivers or geotectonic barriers. Little or no gene flow between these species suggests the development of reproductive barriers or competitive exclusion. Our results suggests that limited marker sampling in recent AF studies may compromise estimates of divergence times and historical demography, and we discuss the effects of such sampling on this and other studies.