How did an annual plant react to Pleistocene glaciations? Postglacial history of Rhinanthus angustifolius in Europe

The impact of climate fluctuations during the Pleistocene on the geographic structure of genetic variation in plant populations is well documented, but there is a lack of studies of annual species at the European scale. The present study aimed to infer the history of the widespread European annual Rhinanthus angustifolius C. C. Gmelin (Orobanchaceae). We explored variation in chloroplast DNA (cpDNA) sequences and amplified fragment length polymorphism (AFLP) in twenty-nine populations covering the entire distribution area of the species. Five AFLP groups were identified, suggesting at least two glacial refugial areas: one area in southwestern Europe and one large eastern area in the Balkan/Caucasus. Recolonization of previously glaciated areas mainly took place from the east of Europe. Despite the difference in life-history traits, the patterns found for the annual R. angustifolius show similarities with those of perennial species in terms of genetic diversity and geographic organization of genetic variation. Although organelle markers have typically been preferred in phylogeographic studies, the cpDNA variation in R. angustifolius did not show any clear geographic structure. The absence of geographic structure in the cpDNA variation may reflect persistence of ancestral polymorphisms or hybridization and introgression with closely-related species.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 98 , 1–13.     

Extensive range persistence in peripheral and interior refugia characterizes Pleistocene range dynamics in a widespread Alpine plant species (Senecio carniolicus, Asteraceae)

Recent evidence suggests that survival of arctic-alpine organisms in peripheral or interior glacial refugia are not mutually exclusive and may both be involved in shaping an organism's Pleistocene history, yet potentially at different time levels. Here, we test this hypothesis in a high-mountain plant (diploid lineage of Senecio carniolicus, Asteraceae) from the Eastern European Alps, in which patterns of morphological variation and current habitat requirements suggest survival in both types of refugia. To this end, we used AFLPs, nuclear and plastid DNA sequences and analysed them, among others, within a graph theoretic framework and using novel Bayesian methods of phylogeographic inference. On the basis of patterns of genetic diversity, occurrence of rare markers, distribution of distinct genetic lineages and patterns of range connectivity both interior refugia in the formerly strongly glaciated central Alps and peripheral refugia along the southern margin of the Alps were identified. The presence of refugia congruently inferred by markers resolving at different time levels suggests that these refugia acted as such throughout several glacial cycles. The high degree of range persistence together with gradual range expansion, which contrasts with the extent of range shifts implied for other Alpine species, is likely responsible for incipient lineage differentiation evident from the genetic data. Replacing a simplistic peripheral vs. interior refugia dualism by more complex models involving both types of refugia and considering different time levels will help identifying common phylogeographic patterns with respect to, for instance, location of refugia and colonization routes and elucidating their underlying genetic and/or ecological causes. DNA sequences have been deposited in GenBank under accession nos. FR796701–FR797793 and nos. HE614296–HE614583.     

Multiple Pleistocene refugia and Holocene range expansion of an abundant southwestern American desert plant species (Melampodium leucanthum,Asteraceae)

Pleistocene climatic fluctuations had major impacts on desert biota in southwestern North America. During cooler and wetter periods, drought‐adapted species were isolated into refugia, in contrast to expansion of their ranges during the massive aridification in the Holocene. Here, we use Melampodium leucanthum (Asteraceae), a species of the North American desert and semi‐desert regions, to investigate the impact of major aridification in southwestern North America on phylogeography and evolution in a widespread and abundant drought‐adapted plant species. The evidence for three separate Pleistocene refugia at different time levels suggests that this species responded to the Quaternary climatic oscillations in a cyclic manner. In the Holocene, once differentiated lineages came into secondary contact and intermixed, but these range expansions did not follow the eastwardly progressing aridification, but instead occurred independently out of separate Pleistocene refugia. As found in other desert biota, the Continental Divide has acted as a major migration barrier for M. leucanthum since the Pleistocene. Despite being geographically restricted to the eastern part of the species’ distribution, autotetraploids in M. leucanthum originated multiple times and do not form a genetically cohesive group.     

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.     

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.     

Chloroplast DNA diversity of the dioecious European tree Ilex aquifolium L. (English holly)

Variation in the chloroplast genome of Ilex aquifolium (English holly), a dioecious evergreen tree native to south, west and central Europe, was analysed using polymerase chain reaction-restriction fragment length polymorphisms (PCR-RFLPs) and microsatellites. Differentiation between populations was high (GST = 0.595) and evidence for phylogeographical structure was detected (NST = 0.697, significantly higher than GST). Two chloroplast lineages were inferred originating from putative glacial refugia in southern Europe (Iberia, Italy and possibly the Balkans). The GST value was higher than reported for endozoochorous hermaphrodite species and for anemochorous dioecious species investigated over a similar geographical scale. It appears that dioecy has contributed to strong differentiation between refugia and that this has been maintained following postglacial recolonization as a result of limited seed flow. Palynological records for I. aquifolium are poor, thus these results give an important insight into patterns of glacial isolation and postglacial recolonization of this species.     

Pleistocene divergence of Dinaric Drusus endemics (Trichoptera, Limnephilidae) in multiple microrefugia within the Balkan Peninsula

The Balkan Peninsula is one of three major European refugial areas. It has high biodiversity and endemism, but data on the age and origin of its fauna, especially endemics, are limited. Mitochondrial sequence data (COI and 16S genes) were used to study the population structure and phylogeography of the caddisfly Drusus croaticus and the phylogeny and divergence of seven other Drusus species, mostly range-restricted endemics of the Dinaric region of the Balkan Peninsula. The divergence of D. croaticus populations in Croatia and allopatric Drusus species in Bosnia dated to the Pleistocene, showing the importance of this time period for the origin and diversification of Balkan endemic taxa. The divergence of more distantly related species dated to the Late Miocene/Early Pliocene. Population genetic and phylogeographic analysis of 115 individuals from 11 populations of D. croaticus revealed a high level of genetic differentiation and absence of gene flow between populations separated by more than 10 km. The existence of allopatrically fragmented lineages in D. croaticus and the endemic Bosnian species is most likely the result of long-term isolation in multiple microrefugia, probably due to the specific habitat requirements and life-history traits of Drusinae coupled with the topographic complexity and historical changes in geomorphology of the region. Overall, these findings shed light on the processes generating the high genetic complexity of this refugial region that parallels the 'refugia within refugia' pattern widely reported from the Iberian refugium.     

Traces of ancient range shifts in a mountain plant group (Androsace halleri complex, Primulaceae)

Phylogeographical studies frequently detect range shifts, both expansions (including long-distance dispersal) and contractions (including vicariance), in the studied taxa. These processes are usually inferred from the patterns and distribution of genetic variation, with the potential pitfall that different historical processes may result in similar genetic patterns. Using a combination of DNA sequence data from the plastid genome, AFLP fingerprinting, and rigorous phylogenetic and coalescence-based hypothesis testing, we show that Androsace halleri (currently distributed disjunctly in the northwestern Iberian Cordillera Cantábrica, the eastern Pyrenees, and the French Massif Central and Vosges), or its ancestor, was once more widely distributed in the Pyrenees. While there, it hybridized with Androsace laggeri and Androsace pyrenaica, both of which are currently allopatric with A. halleri. The common ancestor of A. halleri and the north Iberian local endemic Androsace rioxana probably existed in the north Iberian mountain ranges with subsequent range expansion (to the French mountain ranges of the Massif Central and the Vosges) and allopatric speciation (A. rioxana, A. halleri in the eastern Pyrenees, A. halleri elsewhere). We have thus been able to use the reticulate evolution in this species group to help elucidate its phylogeographical history, including evidence of range contraction.     

Multiple genetic divergences and population expansions of a Mediterranean sandfly,Phlebotomus ariasi,in Europe during the Pleistocene glacial cycles

Phlebotomus ariasi is one of the two sandflies transmitting the causative agent of zoonotic leishmaniasis, Leishmania infantum, in France and Iberia, and provides a rare case study of the postglacial re-colonization of France by a Mediterranean species. Four DNA sequences were analysed—mitochondrial cytochrome b (cyt b), nuclear elongation factor-1α (EF-1α) and two anonymous nuclear loci—for 14–15 French populations and single populations from northeast Spain, northwest Spain, Portugal and Morocco. The presence of cryptic sibling species was not revealed by phylogenetic analyses and testing for reproductive isolation between sympatric populations defined by the two most divergent cyt b haplogroups. No locus was shown to be under positive directional or balancing selection and, therefore, molecular variation was explained demographically. Each nuclear locus showed shallow isolation by distance from Portugal to the French Pyrenees, but for both cyt b and EF-1α there was then a step change to the upland Massif Central, where leading-edge populations showed low diversity at all loci. Multiple genetic divergences and population expansions were detected by analyses of cyt b and dated to the Pleistocene. Endemicity of one cyt b sub-lineage suggested the presence of a refuge north of the Pyrenees during the last glacial period. Monopolization of the Massif Central by genetically differentiated populations of P. ariasi might possibly hinder the northwards spread of leishmaniasis.     

Pleistocene phylogeography and phylogenetic concordance in cold-adapted spring snails (Bythinella spp.)

Previous studies on Pleistocene phylogeography of European taxa are biased towards (i) vertebrates, (ii) terrestrial taxa, (iii) single species, and (iv) taxa that survived the Pleistocene in southern refugia. Relatively little is known about whether evolutionary patterns of vertebrate and terrestrial taxa are also applicable to freshwater invertebrates, whether cold-adapted freshwater species could survive in extensive permafrost areas without retreating into refugia, and whether Pleistocene phylogeographical patterns are influenced by phylogeny. Here, the widespread and species-rich European spring snail genus Bythinella Moquin-Tandon, 1856 is utilized in an attempt to mitigate this bias. These strongly cold-adapted freshwater animals mostly occur in springs — highly isolated habitats that are relatively unaffected by anthropogenic influences. Phylogenetic and phylogeographical analyses based on mitochondrial DNA and nuclear DNA sequence data were conducted in 458 specimens from 142 populations occurring throughout Europe. The study provides evidence that most Bythinella spp. survived the Pleistocene in restricted northern glacial refugia that largely correspond to refugia previously recognized for other European biota. However, survival of Bythinella spp. in extensive permafrost areas outside of refugia can likely be rejected. Low dispersal ability and the isolation and fragmentation of spring habitats, as well as the distribution of perennial springs within permafrost regions, may account for this result. Tests involving a total of 29 nominal species showed that phylogenetically closely related Bythinella species did not occupy similar refugia. This lack of phylogenetic concordance could possibly be explained by the stochasticity of survival and dispersal in spring snails.     

Late Pleistocene origin of the entire circumarctic range of the arctic‐alpine plant Kalmia procumbens

The circumarctic ranges of arctic‐alpine plants are thought to have been established in the late Pliocene/early Pleistocene, when the modern arctic tundra was formed in response to climate cooling. Previous findings of range‐wide genetic structure in arctic‐alpine plants have been thought to support this hypothesis, but few studies have explicitly addressed the temporal framework of the genetic structure. Here, we estimated the demographic history of the genetic structure in the circumarctic Kalmia procumbens using sequences of multiple nuclear loci and examined whether its genetic structure reflects prolonged isolation throughout the Pleistocene. Both Bayesian clustering and phylogenetic analyses revealed genetic distinction between alpine and arctic regions, whereas detailed groupings were somewhat discordant between the analyses. By assuming a population grouping based on the phylogenetic analyses, which likely reflects a deeper intraspecific divergence, we conducted model‐based analyses and demonstrated that the intraspecific genetic divergence in K. procumbens likely originated during the last glacial period. Thus, there is no need to postulate range separation throughout the Pleistocene to explain the current genetic structure in this species. This study demonstrates that range‐wide genetic structure in arctic‐alpine plants does not necessarily result from the late Pliocene/early Pleistocene origin of their circumarctic ranges and emphasizes the importance of a temporal framework of the current genetic structure for understanding the biogeographic history of the arctic flora.     

Chloroplast DNA supports a hypothesis of glacial refugia over postglacial recolonization in disjunct populations of black pine (Pinus nigra) in western Europe

European black pine (Pinus nigra Arn.) is a widely distributed Mediterranean conifer. To test the hypothesis that fragmented populations in western Europe survived in situ during the last glacial rather than having been re-colonized in the postglacial period, genetic variation was assessed using a suite of 10 chloroplast DNA microsatellites. Among 311 individuals analysed, 235 haplotypes were detected revealing high levels of chloroplast haplotype diversity in most populations. Bayesian analysis using a model of linked loci, with no prior assumption of population structure, assigned individuals to 10 clusters that corresponded well with the six predefined sampling regions, while an analysis carried out at the population level and assuming unlinked loci, recovered the original six sampling regions. This regional structure was supported by a biogeographical analysis that detected five barriers, with the two most significant separating Alps from Corsica and southern Italy, and southern Spain from the Pyrenees. No signals of demographic expansion were detected, and comparisons of R(ST) with pR(ST) suggested that a stepwise mutational model was important in regional differentiation, but not in population-within-region differentiation. These tests support long-term persistence of the species within the six regions. The temporal depth estimate, assuming a high mutation rate in coalescent modelling, placed the deepest split between the Alps and the other regions at about 150 000 years ago, and the most recent split of Pyrenees from southern France at about 30 000 years ago. Taken together, the data suggest that chloroplast DNA is structured in black pine and disjunct populations in western Europe are likely to have been present during the Last Glacial Maximum.     

Land masses and oceanic currents drive population structure of Heritiera littoralis,a widespread mangrove in the Indo‐West Pacific

Phylogeographic forces driving evolution of sea‐dispersed plants are often influenced by regional and species characteristics, although not yet deciphered at a large spatial scale for many taxa like the mangrove species Heritiera littoralis. This study aimed to assess geographic distribution of genetic variation of this widespread mangrove in the Indo‐West Pacific region and identify the phylogeographic factors influencing its present‐day distribution. Analysis of five chloroplast DNA fragments’ sequences from 37 populations revealed low genetic diversity at the population level and strong genetic structure of H. littoralis in this region. The estimated divergence times between the major genetic lineages indicated that glacial level changes during the Pleistocene epoch induced strong genetic differentiation across the Indian and Pacific Oceans. In comparison to the strong genetic break imposed by the Sunda Shelf toward splitting the lineages of the Indian and Pacific Oceans, the genetic differentiation between Indo‐Malesia and Australasia was not so prominent. Long‐distance dispersal ability of H. littoralis propagules helped the species to attain transoceanic distribution not only across South East Asia and Australia, but also across the Indian Ocean to East Africa. However, oceanic circulation pattern in the South China Sea was found to act as a barrier creating further intraoceanic genetic differentiation. Overall, phylogeographic analysis in this study revealed that glacial vicariance had profound influence on population differentiation in H. littoralis and caused low genetic diversity except for the refugia populations near the equator which might have persisted through glacial maxima. With increasing loss of suitable habitats due to anthropogenic activities, these findings therefore emphasize the urgent need for conservation actions for all populations throughout the distribution range of H. littoralis.     

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

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

A phylogeographic, demographic and historical analysis of the short-tailed pit viper (Gloydius brevicaudus): evidence for early divergence and late expansion during the Pleistocene

The impact of quaternary glaciation in eastern China on local fanua and flora has been a topic of considerable interest. We used mitochondrial DNA (mtDNA) sequence data and coalescent simulations to test two general biogeographic hypothesis related to the effects of Pleistocene climatic fluctuations for a widespread ophidian species (Gloydius brevicaudus) in eastern China and Korean Peninsula. The phylogenetic analysis revealed three major lineages, the southeast Coastal, Yangtze and North Lineages. The latter two are closely related and jointly form a continental lineage. Divergence dating and coalescent simulations indicate a Late Pliocene to Early Pleistocene divergence between lineages from the southeast coast and continental interior, followed by a mid-to-late Pleistocene divergence between lineages from the north and the middle-lower Yangtze Valley across East China, suggesting that all these lineages predated the last glacial maximum. An overlapping range between the two lineages within the continental lineage and a secondary contact associated with ecological transition zones on the margins of the North China Plain were also observed. These results show that vicariance patterns dominated the history of G. brevicaudus. Though the climatic events of the Pleistocene have had a marked effect on the historical distribution and intra-specific divergence of reptiles in China, coalescent and non-coalescent demographic analyses indicate that all lineages of G. brevicaudus seem not to have been adversely affected by glacial cycles during the Late Pleistocene, presumably because of an increase in the amount of climatically mild habitat in East Asia due to a decline in elevation and the development of monsoons since the Mid-End Pleistocene.     

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.     

Multiple Pleistocene refugia and post‐glacial colonization in the European chub (Squalius cephalus) revealed by combined use of nuclear and mitochondrial markers

Aim To analyse patterns of nuclear and mitochondrial genetic variation in the European chub, Squalius cephalus (Linnaeus, 1758), in order to understand the evolutionary history of this species and to test biogeographical hypotheses for the existence of co‐distributed European freshwater fish species. Location Rivers in Europe (Finland, Poland, Czech Republic, France, Bulgaria, Spain, Italy). Methods We genotyped 12 polymorphic microsatellite markers derived from 310 individuals collected from across the distribution of S. cephalus in Europe (including a total of 15 populations) and sequenced mitochondrial DNA (mtDNA) from a subset of 75 individuals. Sequences of mtDNA cytochrome b were analysed using both phylogenetic (median‐joining networks) and population genetic methods (tests for demographic history, mismatch distributions, Bayesian coalescent analysis). Geographical structure in microsatellite loci was examined using a distance method (F_ST), factorial correspondence analysis (FCA) and a Bayesian clustering method (structure ). Results The mtDNA network showed a clear split into four different haplogroup lineages: Western (separated into Atlantic and Danubian sublineages), Eastern, Aegean (occurring in two distinct sublineages in the Balkans and in Spain) and Adriatic. Our results indicate recent population expansion in the Eastern and Western Atlantic lineages and the admixture of two previously separate sublineages (Atlantic and Danubian) in the Western lineage. Bayesian structure analysis as well as FCA results roughly corresponded to the mtDNA‐based structure, separating the sampled individuals into almost non‐overlapping groups. Main conclusions Our results support hypotheses suggesting origins of extant lineages of freshwater fishes in multiple refugia and the subsequent post‐glacial colonization of Europe via different routes. We confirmed the previously proposed two‐step expansion scenario from the Danube refuge, the existence of a secondary (Atlantic) refuge during the last glaciation (probably in the Rhone River) and population expansion of this lineage. Conspicuous divergences among Mediterranean populations reflect their different origin, as well as their low contribution to the recent genetic pool of chub in central Europe.