Differential success in northwards range expansion between ecotypes of the marble gallwasp Andricus kollari: a tale of two lifecycles

The Marble gallwasp Andricus kollari has a native range divided into two geographically separated lifecycles. In Eastern Europe and Turkey, the lifecycle involves a sexual generation on Turkey oak, Quercus cerris, while in Iberia and North Africa the sexual generation host is cork oak, Q. suber. Over the last 500 years, A. kollari has expanded its range into northern Europe, following human planting of Q. cerris from Italy and the Balkans. We ask: (i) what is the genetic relationship between eastern and western distributions of Andricus kollari? Can we determine which lifecycle is ancestral, and how long ago they diverged? (ii) To what extent have eastern and western native ranges contributed to northwards range expansion? (iii) Is there any evidence for hybridization between the two life cycle types? We present analyses of allozyme data for 13 polymorphic loci and of sequence variation for a 433 bp fragment of the mitochondrial cytochrome b gene. These show: (i) that four haplotype lineages (one in Spain, two in Hungary/Italy and one in Turkey) diverged more or less simultaneously between 1 and 2 million years ago, suggesting the existence of at least four refuges through recent ice age cycles. Our data cannot resolve which lifecycle type is ancestral. (ii) Populations north of putative refuges are divided into two sets. Populations in south‐west France are allied to Spain, while all remaining populations in northern Europe have been colonized from Italy and the Balkans. (iii) The transition from one race to another in south‐west France is marked by abrupt transitions in the frequency of refuge‐specific private alleles and corresponds closely to the northern limit of the distribution of cork oak. Although hybrids were detected in north‐west France, none were detected where the two lifecycles meet in south‐western France. The biology of oak gallwasps predicts that any hybrid zone will be narrow, and limited to regions where Q. cerris and Q. suber meet. Our data suggest that eastern and western A. kollari are effectively separate species.     

Longitudinal range expansion and cryptic eastern species in the western Palaearctic oak gallwasp, Andricus coriarius

The oak gallwasp Andricus coriarius is distributed across the Western Palaearctic from Morocco to Iran. It belongs to a clade of host-alternating Andricus species that requires host oaks in two sections of Quercus subgenus Quercus to complete its lifecycle, a requirement that has restricted the historic distribution and dispersal of members of this clade. Here we present nuclear and mitochondrial sequence evidence from the entire geographic range of A. coriarius to investigate the genetic legacy of longitudinal range expansion. We show A. coriarius as currently understood to be para- or polyphyletic, with three evolutionarily independent (but partially sympatric) lineages that diverged c. 10 million years ago (mya). The similarities in gall structure that have justified recognition of single species to date thus represent either strong conservation of an ancestral state or striking convergence. All three lineages originated in areas to the east of Europe, underlining the significance of Turkey, Iran and the Levant as 'cradles' of gallwasp evolution. One of the three lineages gave rise to all European populations, and range expansion from a putative Eastern origin to the present distribution is predicted to have occurred around 1.6 mya.     

Genetic spatial structure of European common hamsters (Cricetus cricetus)--a result of repeated range expansion and demographic bottlenecks

The spatial genetic structure of common hamsters (Cricetus cricetus) was investigated using three partial mitochondrial (mt) genes and 11 nuclear microsatellite loci. All marker systems revealed significant population differentiation across Europe. Hamsters in central and western Europe belong largely to two allopatric mitochondrial lineages south and northwest of the Carpathian and Sudetes. The southern group, 'Pannonia', comprises populations inside the Carpathian basin (Czech Republic, Hungary) while the second group, 'North', includes hamsters from Belgium, the Netherlands, France, and Germany. Isolation of the lineages is maintained by a combination of geographical and ecological barriers. Both main phylogeographical groups show signs of further subdivision. North is separated into highly polymorphic central German and less polymorphic western populations, which most likely split during late glacial expansion (15,000-10,000 bp). Clock estimates based on haplotype distributions predict a divergence of the two major lineages 85,000-147,000 bp. Expansion times fall during the last glaciation (115,000-10,000 bp) corroborating fossil data, which identify Cricetus cricetus as characteristic of colder climatic phases. Despite the allopatry of mt haplotypes, there is an overlap of nuclear microsatellite alleles between phylogeographical units. Although there are strong evidence that Pannonian hamsters have persisted inside the Carpathian basin over the last 50,000 years, genetic differentiation among European hamsters has mainly been caused by immigration from different eastern refugia. Possible source populations are likely to be found in the Ukrainian and the southern Russian plains--core areas of hamster distribution. From there, hamsters have repeatedly expanded during the Quaternary.     

Natural selection in a postglacial range expansion: the case of the colour cline in the European barn owl

Gradients of variation—or clines—have always intrigued biologists. Classically, they have been interpreted as the outcomes of antagonistic interactions between selection and gene flow. Alternatively, clines may also establish neutrally with isolation by distance (IBD) or secondary contact between previously isolated populations. The relative importance of natural selection and these two neutral processes in the establishment of clinal variation can be tested by comparing genetic differentiation at neutral genetic markers and at the studied trait. A third neutral process, surfing of a newly arisen mutation during the colonization of a new habitat, is more difficult to test. Here, we designed a spatially explicit approximate Bayesian computation (ABC) simulation framework to evaluate whether the strong cline in the genetically based reddish coloration observed in the European barn owl (Tyto alba) arose as a by‐product of a range expansion or whether selection has to be invoked to explain this colour cline, for which we have previously ruled out the actions of IBD or secondary contact. Using ABC simulations and genetic data on 390 individuals from 20 locations genotyped at 22 microsatellites loci, we first determined how barn owls colonized Europe after the last glaciation. Using these results in new simulations on the evolution of the colour phenotype, and assuming various genetic architectures for the colour trait, we demonstrate that the observed colour cline cannot be due to the surfing of a neutral mutation. Taking advantage of spatially explicit ABC, which proved to be a powerful method to disentangle the respective roles of selection and drift in range expansions, we conclude that the formation of the colour cline observed in the barn owl must be due to natural selection.     

Genetic consequences of glacial survival and postglacial colonization in Norway spruce: combined analysis of mitochondrial DNA and fossil pollen

Norway spruce (Picea abies [L.] Karst.) is a broadly distributed European conifer tree whose history has been intensively studied by means of fossil records to infer the location of full‐glacial refugia and the main routes of postglacial colonization. Here we use recently compiled fossil pollen data as a template to examine how past demographic events have influenced the species’ modern genetic diversity. Variation was assessed in the mitochondrial nad1 gene containing two minisatellite regions. Among the 369 populations (4876 trees) assayed, 28 mitochondrial variants were identified. The patterns of population subdivision superimposed on interpolated fossil pollen distributions indicate that survival in separate refugia and postglacial colonization has led to significant structuring of genetic variation in the southern range of the species. The populations in the northern range, on the other hand, showed a shallow genetic structure consistent with the fossil pollen data, suggesting that the vast northern range was colonized from a single refugium. Although the genetic diversity decreased away from the putative refugia, there were large differences between different colonization routes. In the Alps, the diversity decreased over short distances, probably as a result of population bottlenecks caused by the presence of competing tree species. In northern Europe, the diversity was maintained across large areas, corroborating fossil pollen data in suggesting that colonization took place at high population densities. The genetic diversity increased north of the Carpathians, probably as a result of admixture of expanding populations from two separate refugia.     

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

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

Status of Savi's pipistrelle Hypsugo savii (Chiroptera) and range expansion in Central and south‐eastern Europe: a review

相似文献   

Genetic subdivision, glacial refugia and postglacial recolonization in the golden-striped salamander, Chioglossa lusitanica (Amphibia: urodela)

The golden-striped salamander (Chioglossa lusitanica) is an ecologically specialized species, endemic to north-western Iberia. Patterns of genetic variation were assessed at seven polymorphic enzyme loci and one mitochondrial DNA (mtDNA) marker (cytochrome b) in 17 populations across its range. Estimates of enzyme genetic diversity revealed a high degree of genetic subdivision (FST = 0.68), mainly attributable to the existence of two groups of populations. The groups were located, respectively, north and south of the Mondego River, indicating that this river coincided with a major historical barrier to gene flow. A significant decrease in genetic variability from the Mondego northwards was associated with the Douro and Minho rivers. mtDNA sequence variation revealed a congruent pattern of two haplotype groups (d = 2.2%), with a geographical distribution resembling that of allozymes. The pattern and depth of genetic variation is consistent with the following hypotheses: (i) subdivision of an ancestral range of the species prior to the middle Pleistocene; (ii) secondary contact between populations representing historical refugia; (iii) relatively recent range expansion giving rise to the northern part of the species range; and (iv) loss of genetic variation through founder effects during range expansion across major rivers.     

Inferring the geographic origin of a range expansion: Latitudinal and longitudinal coordinates inferred from genomic data in an ABC framework with the program x‐origin

Climatic or environmental change is not only driving distributional shifts in species today, but it has also caused distributions to expand and contract in the past. Inferences about the geographic locations of past populations especially regions that served as refugia (i.e., source populations) and migratory routes are a challenging endeavour. Refugial areas may be evidenced from fossil records or regions of temporal stability inferred from ecological niche models. Genomic data offer an alternative and broadly applicable source of information about the locality of refugial areas, especially relative to fossil data, which are either unavailable or incomplete for most species. Here, we present a pipeline we developed (called x ‐origin ) for statistically inferring the geographic origin of range expansion using a spatially explicit coalescent model and an approximate Bayesian computation testing framework. In addition to assessing the probability of specific latitudinal and longitudinal coordinates of refugial or source populations, such inferences can also be made accounting for the effects of temporal and spatial environmental heterogeneity, which may impact migration routes. We demonstrate x ‐origin with an analysis of genomic data collected in the Collared pika that underwent postglacial expansion across Alaska, as well as present an assessment of its accuracy under a known model of expansion to validate the approach.     

Genome‐wide set of SNPs reveals evidence for two glacial refugia and admixture from postglacial recolonization in an alpine ungulate

Past glaciation events have played a major role in shaping the genetic diversity and distribution of wild sheep in North America. The advancement of glaciers can isolate populations in ice‐free refugia, where they can survive until the recession of ice sheets. The major Beringian refugium is thought to have held thinhorn sheep (Ovis dalli) populations during times of glacial advance. While isolation in the major refugium can account for much of the genetic and morphological diversity seen in extant thinhorn sheep populations, mounting evidence suggests the persistence of populations in smaller minor refugia. We investigated the refugial origins of thinhorn sheep using ~10 000 SNPs obtained via a cross‐species application of the domestic sheep ovine HD BeadChip to genotype 52 thinhorn sheep and five bighorn sheep (O. canadensis) samples. Phylogenetic inference revealed a distinct lineage of thinhorn sheep inhabiting British Columbia, which is consistent with the survival of a group of thinhorn sheep in a minor refugium separate from the Beringian refugium. Isolation in separate glacial refugia probably mediated the evolution of the two thinhorn sheep subspecies, the white Dall's sheep (O. d. dalli), which persisted in Beringia, and the dark Stone's sheep (O. d. stonei), which utilized the minor refugium. We also found the first genetic evidence for admixture between sheep from different glacial refugia in south‐central Yukon as a consequence of post glacial expansion and recolonization. These results show that glaciation events can have a major role in the evolution of species inhabiting previously glaciated habitats and the need to look beyond established refugia when examining the evolutionary history of such species.     

Cold-hardiness in several species of land snails

1. 1.|Mean supercooling points (SCP) of Anguispira alternata (Say) gradually declined in autumn and increased in spring. SCP also declined in autumn in Discus conrkhitei (Newcomb) and Gastrocopta armifera (Say).

2. 2.|Supercooling ability varied among the species, being the greatest in G. armifera.

3. 3.|Prolonged exposure of winter animals to high temperatures resulted in only a modest increase in SCP.

4. 4.|Low temperature survival of D. cronkhitei and G. armifera was slightly poorer than predicted by the mean SCP.

Does vicariance shape biotas? Biogeographical tests of the vicariance model in the north‐west European land snail fauna

Aim To investigate the importance of vicariance in shaping the north‐west European land snail fauna. Location North‐west Europe. Methods We tested whether there is a non‐random congruence, i.e. a clustering of species ranges, using a Monte Carlo procedure with a null model that generates range data sets such that their range size distribution, the species richness distribution of the geographic cells and the spatial autocorrelation of the occurrences of a taxon approach the parameters in the real data set. Biotic elements, groups of species with similar ranges, were delimited with Model based Gaussian clustering. The prediction that closely related species belong to different biotic elements, has been tested with a chi‐square test. Results The distribution areas of the north‐west European land snail species are significantly clustered as predicted by the vicariance model. One widespread and seven regional biotic elements were identified. Contrary to the predictions of the vicariance model, closely related north‐west European land snail species belong significantly more often to the same biotic element than should be expected by chance. Main conclusions The clustering of closely related north‐west European land snail species within the same biotic element indicates that speciation modes other than vicariance were frequent or that the imprint of vicariance on the ranges was obscured by extensive post‐speciational dispersal. Extensive dispersal may have been caused by Pleistocene climatic fluctuations. The core areas of the regionally restricted biotic elements might indicate the positions of glacial refugia of land snails.     

Unravelling the effects of contemporary and historical range expansion on the distribution of genetic diversity in the damselfly Coenagrion scitulum

Although genetic diversity provides the basic substrate for evolution, there are a limited number of studies that assess the impact of recent climate change on intraspecific genetic variation. This study aims to unravel the degree to which historical and contemporary factors shape genetic diversity and structure across a large part of the range of the range‐expanding damselfly Coenagrion scitulum (Rambur, 1842). A total of 525 individuals from 31 populations were genotyped at nine microsatellites, and a subset was sequenced at two mitochondrial genes. We inferred the importance of geography, environmental factors, and recent range expansion on genetic diversity and structure. Genetic diversity decreased going westwards, suggesting a signature of historical post‐glacial expansion from east to west and the presence of eastern refugia. Although genetic differentiation decreased going northwards, it increased in the northern edge populations, suggesting a role of contemporary range expansion on the genetic make‐up of populations. The phylogeographical context was proven to be essential in understanding and identifying the genetic signatures of local contemporary processes. Within this framework, our results highlight that recent range expansion of a good disperser can decrease genetic diversity and increase genetic differentiation which should be considered when devising suitable conservation strategies.     

Out of Anatolia: longitudinal gradients in genetic diversity support an eastern origin for a circum-Mediterranean oak gallwasp Andricus quercustozae

Many studies have addressed the latitudinal gradients in intraspecific genetic diversity of European taxa generated during postglacial range expansion from southern refugia. Although Asia Minor is known to be a centre of diversity for many taxa, relatively few studies have considered its potential role as a Pleistocene refugium or a potential source for more ancient westward range expansion into Europe. Here we address these issues for an oak gallwasp, Andricus quercustozae (Hymenoptera: Cynipidae), whose distribution extends from Morocco along the northern coast of the Mediterranean through Turkey to Iran. We use sequence data for a fragment of the mitochondrial gene cytochrome b and allele frequency data for 12 polymorphic allozyme loci to answer the following questions: (1) which regions represent current centres of genetic diversity for A. quercustozae? Do eastern populations represent one refuge or several discrete glacial refugia? (2) Can we infer the timescale and sequence of the colonization processes linking current centres of diversity? Our results suggest that A. quercustozae was present in five distinct refugia (Iberia, Italy, the Balkans, southwestern Turkey and northeastern Turkey) with recent genetic exchange between Italy and Hungary. Genetic diversity is greatest in the Turkish refugia, suggesting that European populations are either (a) derived from Asia Minor, or (b) subject to more frequent population bottlenecks. Although Iberian populations show the lowest diversity for putatively selectively neutral markers, they have colonized a new oak host and represent a genetically and biologically discrete entity within the species.     

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.     

Inferring the mode of colonization of the rapid range expansion of a solitary bee from multilocus DNA sequence variation

Rapid geographic range expansions can have dramatic effects on the distribution of genetic diversity, both within and among populations. Based on field records collected over the past two decades in Western Europe, we report on the rapid geographic range expansion in Colletes hederae, a solitary bee species. To characterize how this expansion shaped the distribution of genetic diversity within and among populations, we performed a genetic analysis based on the sequencing of three nuclear loci (RNAp, CAD and WgL). We then simulated the evolution of DNA sequences under a spatially explicit model of coalescence to compare different hypotheses regarding the mode of colonization associated with this rapid expansion and to identify those that are most consistent with the observed molecular data. Our genetic analyses indicate that the range expansion was not associated with an important reduction in genetic diversity, even in the most recently colonized area in the United Kingdom. Moreover, little genetic differentiation was observed among populations. Our comparative analysis of simulated data sets indicates that the observed genetic data are more consistent with a demographic scenario involving relatively high migration rates than with a scenario based on a high reproduction rate associated with few migrants. In the light of these results, we discuss the factors that might have contributed to the rapid geographic range expansion of this pollen‐specialist solitary bee species across Western Europe.     

Spatial genetic structure in Beta vulgaris subsp. maritima and Beta macrocarpa reveals the effect of contrasting mating system,influence of marine currents,and footprints of postglacial recolonization routes

Understanding the factors that contribute to population genetic divergence across a species' range is a long‐standing goal in evolutionary biology and ecological genetics. We examined the relative importance of historical and ecological features in shaping the present‐day spatial patterns of genetic structure in two related plant species, Beta vulgaris subsp. maritima and Beta macrocarpa. Using nuclear and mitochondrial markers, we surveyed 93 populations from Brittany (France) to Morocco – the southern limit of their species' range distribution. Whereas B. macrocarpa showed a genotypic structure and a high level of genetic differentiation indicative of selfing, the population genetic structure of B. vulgaris subsp. maritima was consistent with an outcrossing mating system. We further showed (1) a strong geographic clustering in coastal B. vulgaris subsp. maritima populations that highlighted the influence of marine currents in shaping different lineages and (2) a peculiar genetic structure of inland B. vulgaris subsp. maritima populations that could indicate the admixture of distinct evolutionary lineages and recent expansions associated with anthropogenic disturbances. Spatial patterns of nuclear diversity and differentiation also supported a stepwise recolonization of Europe from Atlantic‐Mediterranean refugia after the last glacial period, with leading‐edge expansions. However, cytoplasmic diversity was not impacted by postglacial recolonization: stochastic long‐distance seed dispersal mediated by major oceanic currents may mitigate the common patterns of reduced cytoplasmic diversity observed for edge populations. Overall, the patterns we documented here challenge the general view of reduced genetic diversity at the edge of a species' range distribution and provide clues for understanding how life‐history and major geographic features interact to shape the distribution of genetic diversity.     

Glacial bottleneck and postglacial recolonization of a seed parasitic weevil, Curculio hilgendorfi, inferred from mitochondrial DNA variation

Climatic changes during glacial periods have had a major influence on the recent evolutionary history of living organisms, even in the warm temperate zone. We investigated phylogeographical patterns of a weevil Curculio hilgendorfi (Curculionidae), a host-specific seed predator of Castanopsis (Fagaceae) growing in the broadleaved evergreen forests in Japan. We examined 2709 bp of mitochondrial DNA for 204 individuals collected from 62 populations of the weevil. Four major haplogroups were detected, in southwestern and northeastern parts of the main islands and in central and southern parts of the Ryukyu Islands. The demographic population expansion was detected for the two groups in the main islands but not for the Ryukyu groups. The beginning time of the expansion was dated to 39 000–59 000 years ago, which is consistent with the end of the last glacial period. Our data also demonstrated that the southwestern population of the main islands has experienced a more severe bottleneck and more rapid population growth after glacial ages than the northeastern population. At least three refugial areas in the main islands were likely to have existed during the last glacial periods, one of which had not previously been recognized by analyses of intraspecific chloroplast DNA variation of several plant species growing in the broadleaved evergreen forests. Our results represent the first phylogeographical and population demographic analysis of an insect species associated with the broadleaved evergreen forests in Japan, and reveal more detailed postglacial history of the forests.     

Genetic evidence for the origins of range disjunctions in the Australian dry sclerophyll plant Hardenbergia violacea

Aim The distribution of genetic variation in the Australian dry sclerophyll plant Hardenbergia violacea (Fabaceae) is examined in the context of Pleistocene climate change in order to identify likely refugia. Particular consideration is given to the origin of range disjunctions in South Australia and Tasmania, and to determining whether the Tasmanian population is indigenous or recently introduced from mainland Australia. Location Southeastern Australian mainland and Tasmania. Methods A combination of chloroplast polymerase chain reaction–restriction fragment length polymorphism and genomic amplified fragment length polymorphism (AFLP) marker systems was used to examine the genetic structure of 292 individuals from 13 populations across the range of H. violacea in southeastern Australia. Results Hardenbergia violacea populations in Tasmania and southern Victoria were characterized by low, almost monotypic chloroplast diversity. New South Wales showed higher haplotype diversity and haplotype sharing among widely distributed populations. Principal coordinates analysis (PCoA) of the AFLP data found a strong latitudinal cline in AFLP variation from northern New South Wales south to Tasmania. The Tasmanian population formed an isolated and somewhat disjunct genetic cluster at one end of this cline. However, the South Australian population was an exception to the clinal variation shown by all other populations, forming a highly disjunct cluster in the PCoA. Within‐population genetic diversity was low in both disjunct populations. Main conclusions The genetic evidence indicates that the Tasmanian population is likely to be indigenous and probably the product of vicariance, which was followed by range contraction at the Last Glacial Maximum or an earlier glacial event. The deep phylogenetic disjunction in South Australia is evidence of a much earlier separation on mainland Australia. The chloroplast structure indicates that, during the Pleistocene, H. violacea underwent broad‐scale recolonization in southern Victoria and Tasmania, possibly from a large continental refugium in eastern New South Wales. We conclude that H. violacea, and presumably the sclerophyll communities in which it occurs, have undergone multiple range contractions to large continental refugia during different Pleistocene glaciations in southeastern Australia.     

Repeated Range Expansion and Glacial Endurance of Potentilla glabra (Rosaceae) in the Qinghai-Tibetan Plateau

To date, little is still known about how alpine species occurring in the Qinghai-Tibetan Plateau (QTP) responded to past climatic oscillations. Here, by using variations of the chloroplast trnJ-L, we examined the genetic distribution pattern of 101 individuals of Potentilla glabra, comprising both the interior QTP and the plateau edge. Phylogenetic and network analyses of 31 recovered haplotypes identified three tentative clades (A, B and C). Analysis of molecular variance (amova) revealed that most of the genetic variability was found within populations (0.693), while differentiations between populations were obviously distinct (Fst = 0.307). Two independent range expansions within clades A and B occurring at approximately 316 and 201 thousand years ago (kya) were recovered from the hierarchical mismatch analysis, and these two expansions were also confirmed by Fu's Fs values and 'g' tests. However, distant distributions of clade C and private haplotypes from clades A and B suggest that they had survived the Last Glacial Maximum (LGM) and previous glaciers in situ since their origins. Our findings based on available limited samples support that multiple refugia of a few cold-enduring species had been maintained in the QTP platform during LGM and/or previous glacial stages.