Nuclear and plastid haplotypes suggest rapid diploid and polyploid speciation in the N Hemisphere Achillea millefoliumcomplex (Asteraceae)
© Guo et al; licensee BioMed Central Ltd. 2012
Received: 8 April 2011
Accepted: 3 January 2012
Published: 3 January 2012
Species complexes or aggregates consist of a set of closely related species often of different ploidy levels, whose relationships are difficult to reconstruct. The N Hemisphere Achillea millefolium aggregate exhibits complex morphological and genetic variation and a broad ecological amplitude. To understand its evolutionary history, we study sequence variation at two nuclear genes and three plastid loci across the natural distribution of this species complex and compare the patterns of such variations to the species tree inferred earlier from AFLP data.
Among the diploid species of A. millefolium agg., gene trees of the two nuclear loci, ncpGS and SBP, and the combined plastid fragments are incongruent with each other and with the AFLP tree likely due to incomplete lineage sorting or secondary introgression. In spite of the large distributional range, no isolation by distance is found. Furthermore, there is evidence for intragenic recombination in the ncpGS gene. An analysis using a probabilistic model for population demographic history indicates large ancestral effective population sizes and short intervals between speciation events. Such a scenario explains the incongruence of the gene trees and species tree we observe. The relationships are particularly complex in the polyploid members of A. millefolium agg.
The present study indicates that the diploid members of A. millefolium agg. share a large part of their molecular genetic variation. The findings of little lineage sorting and lack of isolation by distance is likely due to short intervals between speciation events and close proximity of ancestral populations. While previous AFLP data provide species trees congruent with earlier morphological classification and phylogeographic considerations, the present sequence data are not suited to recover the relationships of diploid species in A. millefolium agg. For the polyploid taxa many hybrid links and introgression from the diploids are suggested.
Species complexes or aggregates consist of a set of closely related species often of different ploidy levels, whose relationships are difficult to reconstruct. Such species complexes are common in angiosperms [1–5]. Rapid genetic, phenotypic and ecological differentiation on one hand, and hybridization/polyploidy on the other, play important roles in their evolutionary bursts [6–9]. The temperate N Hemisphere common yarrow taxa form such a complex, i.e., the Achillea millefolium aggregate. Centered in SE Europe and SW to C Asia, its diploid species are limited to Eurasia, whereas the polyploids have spread throughout the N Hemisphere [10, 11]. In N America, the 4x and 6x cytotypes form a complex of ecological races adapted to many different niches with marked genotypic diversification [6, 12, 13]. By cultivation in experimental gardens, Clausen et al.  documented local adaptation of A. millefolium populations to environments along an altitudinal transect in California from sea level to alpine regions. This has become a classic example of rapid adaptive divergence of plant populations [12–18].
Our earlier AFLP data have characterized the Achillea millefolium complex as a clade . The inferred species relationships of the diploid members conformed to a binary bifurcating tree and generally agreed with traditional species delimitations and taxonomic arrangements. The polyploid members appeared to be polytomic and polyphyletic although geographic patterns can be recognized [10, 11]. Available data show that frequent exchange of genetic materials have been involved in the origins of many polyploidy taxa. Meanwhile, it seems also to have occurred during the divergence of the diploid species [10, 11, 19]. So far, we are still uncertain about the demographic history of A. millefolium agg. and the progenitors of many polyploid taxa.
Due to their dominant nature, AFLP markers are difficult to use for inferring genetic parameters of populations, particularly ancestral population sizes, split times, and migration rates. In principle, this can be accomplished better with DNA sequences either from organelles or from the nuclear genome [20, 21]. Yet, we still meet challenges in practice: Plastid DNA variation is often too low to infer relevant gene trees with confidence. In addition, despite the lower effective population size of plastids, plastid gene trees may still not reflect the species tree due to incomplete lineage sorting. With nuclear genes, frequent birth and death of gene copies, less lineage sorting due to higher effective population size, secondary introgression among split species as well as intragenic recombination tend to hamper interpretation of the patterns of polymorphism [22–25].
Here, we survey DNA sequence variation sampled at two nuclear loci and three plastid fragments from populations across diploid, tetraploid, and hexaploid species throughout the natural distribution range of A. millefolium agg. and from three diploid congeneric species outside the aggregate. On the diploid level, we infer the demographic history among species and populations using the newly generated DNA sequence data in comparison to the relationships inferred from the previous AFLP data. To this end, we also apply a probabilistic model (IMa2) [26, 27] to three widespread diploid species to shed light on the key parameters, i.e., ancestral effective population sizes, time of speciation, and rates of gene flow. For the polyploid populations, we try to untangle their polytomic and polyphyletic relationships which are probably complicated by gene flow on the same and between different ploidy levels using the co-dominant single-gene haplotype data.
Taxa, populations and DNA loci sampled
Geographic locality of populations
Collectors and dates
Nuclear genes: number of indiv./clones analyzed
cpDNA: number of indiv. analyzed
Achillea millefolium agg.
A. asiatica Serg. s.lat. (= A. sergievskiana Shaulo & Shmakov)
China: Daqing Mt., 41°04'52" N, 112°35'56" E; 2010 m
China: Arxan Mt., 47°17'39" N, 120°27'11" E; 1130 m
China: Hebei, 42°26' N, 117°15' E; 1500 m
Russia: Altai, 51°02'52" N, 85°36'47" E; 1100 m
A. asiatica Serg. s.str.
Russia: Altai, 49°32'66" N, 88°13'35" E; 2350 m
Russia: S Siberian lowland near Novosibirsk; 220 m
Uzbekistan: Tashkent, Tschimgan Mt.
A. asplenifolia Vent.
Austria: Burgenland, Zitzmannsdorfer Wiesen
Austria: Burgenland, Rust, near lake "Neusiedler See"
FE, JS, YG, 2003.05.27
Czech Republic: South Moravia, Terezin
FE & LE, 2002.07.12
A. borealis Bong. s.lat. (A. lanulosa Nutt. var. alpicola Rydb.)
USA: Washington, Mt. Rainier National Park; 1350-2070 m
PS & AT, 2002.08.18
A. borealis Bong. s.lat. (A. lanulosa Nutt. var.?)
USA: Connecticut, Hopeville State Park
A. borealis Bong. s.lat. (A. lanulosa Nutt. var. lanulosa)
USA: Utah, ascent from Snowbird Alta to Lake Secret
A. ceretanica Sennen s.str.
France: E Pyrenees
A. ceretanica Sennen s.lat.
France: Massif Central
A. cuspidata Wall.
Kashmir: 34°25.80' N, 75°44.80' E; 3200 m
A. distans Waldst. & Kit. ex Willd.
Austria, Kaltenleutgeben, Flösslberge
A. inundata Kondr.
Ukraine: Kiev, S of Desna mouth into Dnjepr; 100 m
FE & YG, 2003.07.28
A. latiloba Ledeb. ex Nordm.
Georgia: Adjara, 41°29'55″ N, 42°31'46″ E; 2006 m
A. millefolium L. s.lat. (A. apiculata N.I.Orlova)
Russia: Karelia, Louhski region, (a) Kandalaksha Natural Reserve; (b) Kiv bay, near Medvezhij peninsula; (c) near cape Ivanov Navolok
A. millefolium L. s.lat. (A. sudetica Opiz)
Austria: Salzburg, Hohe Tauern; ca. 2300 m
A. roseoalba Ehrend.
Slovenia: Ljubljana, Podpec
Italy: Udine, Kanalta, MalborghettoValbruna
A. schmakovii Kupr.
Russia: Altai, 51°02'52'' N, 85°36'47'' E; 1700 m
A. setacea Waldst. & Kit.
Greece: Thessaloniki, drain from lake Limni Koronia
Ukraine: Kiev, Bald Mt., Lisa Gora
FE & YG 2003.07.22
Austria: Burgenland, E of St. Margarethen; ca 200 m
FE, JS, YG, 2003.05.27
Turkey: Anatolia, Aksaray
A. styriaca J. Saukel & J. Danihelka, ined.
Austria: Styria, Einach, Wald
Species outside A. millefolium agg.
A. acuminata (Ledeb.) Schultz Bip.
China: Jilin, Changbai Mt., Hancong Valley, 680-620 m
YG & GR, 2002.07.24
China: Inner Mongolia, Arxan, N 47°17'39.5", E 120°27'09.9"; 865 m
China: Shanxi, Taibai Mt., N 34°01'17", E 107°18'21"; 1700 m.
A. ligustica Vis.ex Nym.
Italy: Sicily, Nebrodi Mts.
A. nobilis L.
Czech Rep.: Znoimo
LE & FE, 2002.07.13
To check ploidy levels of the populations studied, two methods, either chromosome counting or DNA ploidy level determination, were applied using young flower buds or fresh or silica gel-dried leaves, respectively. Young flower buds were collected in the field and fixed in Carnoy's fluid (ethanol:acetic acid = 3:1). To count the chromosome number, fixed flower buds were stained and squashed in 4% acetocarmine and observed under the microscope. DNA ploidy levels were investigated with propidium iodide flow cytometry [28, 29] from the prepared leaves. Information for ploidy levels obtained with the above two methods are marked in Table 1 by c and f, respectively, while those inferred from previous studies [12 and Ehrendorfer et al., unpubl. Data] are marked by i.
Voucher specimens are deposited in the herbaria of the Institute of Botany (WU) and the Department of Pharmacognosy (HBPh), both at the University of Vienna, Austria, and of the College of Life Sciences, Beijing Normal University (BNU).
Total genomic DNA was extracted from ca. 0.02 g silica gel desiccated leaf materials following the 2 × CTAB protocol  with slight modifications: Before the normal extraction process, sorbitol washing buffer was used to remove polysaccharides in the leaf materials (add 800 μL sorbitol buffer to the ground leaf powder → incubate the sample in ice for 10 min → centrifuge at 10,000 g for 10 min at 4°C and then follow the established 2 × CTAB protocol).
Two nuclear genes were sampled and partially sequenced for this study. They are the chloroplast-expressed Glutamine Synthase gene (ncpGS) and the Sedoheptulose-Bisphosphatase (SBP) gene. The ncpGS gene has been used in many plant phylogenetic studies and shown to be single-copy in all angiosperm species so far studied  and especially in Achillea . We sequenced part of its coding and noncoding regions from exon 7 through to exon 11. The SBP gene has been studied in representative taxa of the family Asteraceae . It was shown to be single-copy by our preliminary analyses in several diploid species of Achillea. Readers are referred to Ma et al.  for primers used for amplifying the ncpGS locus and to Chapman et al.  for the SBP locus (its exon 5 through to 7, i.e. the locus B12 in Chapman et al. 2007).
Three noncoding chloroplast DNA regions, trnH-psbA, trnC-ycf6 (including partial ycf6-psbM) and rpL16 were sequenced. PCR reactions were conducted with universal primer pairs .
The amplification was carried out in a volume of 20 μL with final concentration of 1 × PCR buffer, 0.05 U exTaq (TaKaRa, Shiga, Japan) or HiFi (TransTaq DNA polymerase High Fidelity, TransGen Biotech), 200 μM of each dNTP, 1% DMSO, 0.5 μM of each primer, and with 1 μL template DNA and ddH2O added to the final volume. The amplification was conducted on a Peltier thermocycler (Bio-RAD) initiated with 5 min of pre-denaturing at 94°C followed by 30 cycles of 1 min at 94°C, 30s at 48-55°C, and 1.5 min at 72°C. A final extension was then taken at 72°C for 15 min followed by a hold at 4°C. The PCR products were electrophoresed on and excised from the 1.0% agarose gel in TAE buffer. They were then purified using a DNA Purification kit (TianGen Biotech or TransGen Biotech, Beijing, China). The purified PCR products were either used for direct sequencing (for the cpDNA fragments) or ligated into a pGEM-T Vector (for nuclear genes) with a Promega Kit (Promega Corporation, Madison, USA). For sequencing the nuclear genes, about five to eight positive clones from each diploid and ten to fifteen from each polyploid individual were randomly selected for sequencing. The plasmid was extracted with an Axyprep Kit (Axygene Biotechnology, Hangzhou, China). Cycle Sequencing was conducted using ABI PRISM® BigDye™ Terminator. The same primers used for amplification (for cpDNA fragment) or the vector primers T7/Sp6 (for nuclear genes) were applied here. The sequenced products were run on an ABI PRISM™ 3700 DNA Sequencer (PE Applied Biosystems).
Sequences were assembled with the ContigExpress program (Informax Inc. 2000, North Bethesda, MD), aligned with ClustalX 1.81, and then manually improved with BioEdit version 7.0.1. To prevent possible sequencing errors, single mutations in the nuclear gene data sets likely generated by the cloning sequencing method were excluded from the analyses. Furthermore, unique sequences in the nuclear gene data matrix, which do not fall into any majority-rule consensus sequence group [19, 34] or show inconstant branch positions in trees based on different subsets of data, i.e., with partial characters or randomly selected sequences, during the initial analyses were eliminated to avoid influence of PCR-mediated recombination [19, 35, 36]. The final numbers of individuals/clones analyzed at each locus for each population are listed in Table 1. All the sequences analyzed were submitted to the NCBI GenBank under accession numbers HQ601971-HQ602593 (the nuclear ncpGS and SBP genes) and HQ450864-HQ451071 (the plastid loci).
The allelic data sets of the two nuclear genes, ncpGS and SBP, were analyzed separately, whereas sequences of three cpDNA fragments were combined as one locus.
Gaps in the nuclear data sets were treated as missing data, whereas each indel position (no matter how many nucleotide sites it contained) of the plastid data set was coded as a binary character (0/1 = A/C) using the program GapCoder .
As A. millefolium agg. consists of species with short evolutionary history [10, 11], Neighbor Joining (NJ), Maximum Parsimony (MP) and Median-Joining network were applied to the present data. For the nuclear sequences, Neighbour Joining and Parsimony analyses were performed with MEGA 5.05 and PAUP* 4.0b10a, respectively. All nucleotide substitutions were equally weighted. Gaps were treated as missing data. We first analyzed data of the diploid species to show diversification of the gene lineages, and then of all the taxa to investigate relationships among the polyploids and diploids within A. millefolium agg. The NJ analysis was conducted with Kimura's 2-parameter distances  and bootstrapped with 1000 replicates. For the MP method, heuristic searches were performed using 1000 random taxon addition replicates with ACCTRAN optimization and TBR branch swapping. Up to 10 trees with scores larger than 10 were saved per replicate. The stability of internal nodes of the MP tree was assessed by bootstrapping with 1000 replicates (MulTrees option in effect, TBR branch swapping and simple sequence addition).
Median-Joining network analysis implemented in Network ver. 18.104.22.168 available at http://www.fluxus-engineering.com/sharenet.htm was applied to the cpDNA data set. All variable sites were equally weighted and the homoplasy level parameter (ε) was set to zero given that variation rates of the closely related species is low, especially in their plastid DNA.
To understand the population demography at the time of speciation of the diploid species of A. millefolium agg., we applied a probabilistic model, the Isolation with Migration Model for multiple populations implemented in IMa2 , to three widespread and closely related species A. asplenifolia-2x and A. setacea-2x and A. asiatica-2x. These species are here regarded as three diverged populations which share nuclear sequence variation. Shared alleles could reflect ancestral polymorphism or gene flow after separation of the populations or species. Assuming neutrality, retention of ancestral polymorphisms is likely if speciation is fast relative to drift, which is inverse in intensity to the effective population size. As a rule of thumb, species are well separated with little ancestral polymorphism and thus almost complete lineage sorting, if the time of separation is at least as long as four times the effective population size . Secondary genetic exchange between the diverging species can also lead to shared alleles observed . The multipopulation model IMa2 allows both ancestral polymorphism and gene flow subsequent to divergence. It assumes a known history of the sampled populations, which can be represented by a rooted bifurcating tree. In earlier analysis using AFLP data , we inferred the rooted species tree as: ((A. asiatica, A. asplenifolia), A. setacea). We note that Ima2 provides posterior distributions of parameters, such that the confidence in the inference of each parameter can be obtained from observing the spread of the posterior distribution. The IM model also assumes neutral genetic variation, freely recombining unlinked loci and no intragenic recombination or gene conversion . Sequences of the two nuclear loci, the ncpGS and the SBP genes, and of three plastid fragments were used for this analysis. The polymorphic sites of the sequenced nuclear and plastid loci are mostly of introns or intergenic spacers and thus should fit the neutral variation model. Using the four-gamete criterion , we do not find intragenic recombination in the nuclear sequences among these three species. The data of the three plastid fragments were combined because the chloroplast genes are generally linked and no evidence of recombination between the three regions is found.
To run IMa2, one random haplotype per plant individual was chosen for the nuclear gene data sets, and the plastid data set was composed of sequences from the same plant individuals. This avoids bias but decreases the amount of information and thus leads to broader posterior distributions. The IS (Infinite Sites) model  of sequence evolution was chosen for the plastid locus, whereas, the HKY model  which allows for multiple substitutions was selected for the two nuclear loci because double mutations were found for a few polymorphic sites at both loci. The inheritance scalar was set to 1.0 for the nuclear and to 0.25 for the plastid loci, respectively.
To set upper bounds on the prior distributions of the parameters, we estimated for each of the three species the geometric means of the population mutation rate 4Nu across all three loci using Watterson's estimator θ (per sequence not per site). The largest mean value was found with A. asiatica-2x (an estimate of 4Nu = 9.8205), and this was used to set the upper bound on uniform prior for each of the three population demographic parameters: population size (θ = 4Nu), splitting time (t = Tu, where T is the time in generations since the common ancestry, and it is of the same order of 4N) and migration rates (2NM = 4Nu × m/2). The priors were finally set as follows: the upper bound of population sizes q = 100, splitting times t = 5 and migration rates m = 2.0, respectively. We ran the Markov-chain Monte Carlo (MCMC) simulations with 1,000,000 burn-in steps and 20,000 genealogies sampled per locus. The analysis was done with 10 independent runs in the M mode, each using identical priors and 20 Metropolis-Coupled chains with different random number seeds. The genealogies sampled from the M mode runs were combined in an L mode run to build an estimate of the joint posterior probability of the parameters [26, 27].
Nuclear gene trees with allelic haplotype sequences
Amplification of the partial ncpGS and SBP genes produced a single clear band for each amplification. This and the results from earlier work (see "Methods") suggest sequences of the two nuclear loci obtained here each as belonging to a set of orthologs.
After eliminating some sequences likely containing PCR-recombination (about 10% of the total), 303 sequences (clones) of the ncpGS and 313 of the SBP gene from broadly the same 70 individuals of 29 populations belonging to A. millefolium agg. were used for the data analyses (Table 1). In addition, 56 ncpGS and 36 SBP sequences from five populations of three congeneric species outside the A. millefolium agg., A. nobilis-2x, A. ligustica-2x and A. acuminata-2x, were also analyzed here. The ncpGS alignment contains 918 nucleotide positions with sequence length varying from 795 to 861 bps. The SBP alignment contains 420 nucleotide positions with sequence length varying from 386 to 405 bps.
The diploid-only ncpGS data come from seven species within and three outside A. millefolium agg.. They contain 31 haplotypes with 134 substitution sites from 186 sequences (Additional file 1: S-Figure 1). Out of the 134 polymorphic sites, 122 are in introns. Intragenic recombination among some samples is likely: A discordance in the alignment can be resolved by postulating a recombination in three of the four ncpGS haplotypes of A. nobilis-2x around the 89th polymorphic site (Additional file 1: S-Figure 1). Another discordance can be resolved in the two haplotypes of A. cuspidata-2x (III in Figure 1A) by postulating a recombination around the 26th polymorphic site between haplotype groups of A. millefolium agg. and of A. ligustica-2x (Additional file 1: S-Figure 1).
Figure 1A shows an unrooted Neighbour Joining phylogram based on the diploid-only ncpGS data. The topology of the MP tree on the same data set is broadly comparable, and thus only bootstrap values from the MP analysis are presented in Figure 1A. In this gene tree, haplotypes of A. millefolium agg. fall into two major groups, group I corresponding to A. setacea-2x, and group II those of the other Eurasian diploid species except A. cuspidata-2x. This relationship agrees with that from the previous AFLP analysis (Figure 1B inferred from Guo et al., 2005). In haplotype group II, both the subgroups IIa and IIb harbor sequences of A. asplenifolia-2x and A. asiatiaca-2x. This is not in line with the AFLP tree and is most likely due to retention of ancestral polymorphism or secondary contacts between the two species. A. asplenifolia-2x shares its diverged alleles with A. ceretanica-2x and A. roseoalba-2x, respectively, corresponding to the previous AFLP tree (Figure 1B). Sequences of A. ceretanica-2x (relict in the Eastern Pyrenees) only appear in IIb and those of A. roseoalba-2x (from the meadows in the S-Alps) only in IIa. The sequences of A. cuspidata-2x (in the Himalayas) appear distant to other members of A. millefolium agg. but close to A. ligustica-2x (Figure 1A). The single plant sample of A. cuspidata-2x was collected recently and thus was not included in the previous AFLP study. Clearly more samples of this species should be investigated. Out of the four haplotypes of A. nobilis-2x, one falls into clade IIa together with several members of A. millefolium agg., and three form a group between IIa and IIb. The latter can be explained by the intragenic recombination visible in the alignment (Additional file 1: S-Figure 1). Considering relationships of A. nobilis with A. millefolium agg., the ncpGS gene tree is neither congruent with the AFLP tree (Figure 1B) nor with the SBP tree (see below, Figure 2).
The diploid-only SBP data also include seven species within and three outside A. millefolium agg. They contain 21 haplotypes with 60 substitution sites based on 163 sequences. Out of the 60 polymorphic sites, 48 belong to the intron regions. The alignment of the SBP sequences does not show obvious intragenic recombination (Additional file 2: S-Figure 2). The topology of the NJ tree is broadly comparable with that of the MP tree and thus only bootstrap values from the MP analysis are presented (Figure 2). The SBP gene tree (Figure 2) is remarkably incongruent with the ncpGS tree (Figure 1A) and with the AFLP tree (Figure 1B). In Figure 2, haplotypes belonging to members of A. millefolium agg. do not group together. Some of A. asiatica-2x, all of the Caucasus A. latiloba-2x and the Himalayan A. cuspidata-2x (here defined as Asian types) are distantly related to the others of A. millefolium agg. Surprisingly, sequences of the C European A. nobilis-2x are close to the Asian type, whereas, haplotypes of the E Asian A. acuminata-2x are close to the major haplotype group of A. millefolium agg., which is mostly of the European members. We thus observe little sorting of ancestral polymorphisms of the SBP gene during the speciation processes of the diploid species of Achillea.
The date set of SBP gene of all the diploids and polyploids contains 68 haplotypes generated from 349 sequences with 68 polymorphic substitution sites. Due to the severe conflicts between the SBP gene tree and the species tree inferred from the AFLP and morphological data, sequences of this gene are not suitable for the phylogenetic inference, but could provide some clues about the progenitors of the polyploid taxa. We therefore only present the SBP gene tree of all the diploid and polyploid samples in the supplementary materials as Additional file 3: S-Figure 3.
Phylogenetic networks based on plastid haplotypes
Sequence characters of the analyzed cpDNA non-coding regions
trnC-ycf6(incl. partial ycf6-psbM)
rpL16(for diploid taxa only)
All three fragments combined for diploid taxa only
trnH-psbAand trnC-ycf6combined for all taxa
Number of haplotypes
Number of variable sites
Number of indels (length in base pairs) *
3 (5; 1; 21)
5 (1; 6; 1; 21)
4 (5; 5; 1; 22)
8 (1; 6; 5; 1; 5; 5; 1; 22)
8 (1; 6; 1; 6; 21; 5; 1; 21)
H1 is the most frequent haplotype in Eurasia. It is shared by three out of four European diploid species and is spread among most of the Eurasian polyploid taxa at all ploidy levels (Figure 5). The rare H11 in Europe (only in A asplenifolia-2x) is related to the Eurasian H5, to most of the Asian types H 6-9 &15, and even to H14 which is the most frequent in N America. The European polyploid specific H3 (in A. millefolium subsp. apiculata-6x) and H13 (A. millefolium subsp. sudetica-6x & A. inundata-4x) are directly or indirectly related to H4 found in A. asiatica-2x & -4x. The N American frequent H14 has been derived most probably from the E Asian H10 and evidently gave rise within N America to several rare and more local haplotypes (H18-22, 24, 26).
The 26 plastid haplotypes are mapped associated with the general distribution of populations studied here in the Additional file 4: S-Figure 4.
Demography of major diploid lineages
In contrast, all ancestral population sizes are large compared to the current ones: q3, the ancestral effective population size before the split of A. asiatica-2x and A. asplenifolia-2x, and q4, the ancestral population common to all three species are bounded away from zero. The distributions of the migration parameters (panel C) are broad and quite similar to the equal prior distributions, such that we conclude that there is little information for inference of these parameters in the data. Setting wider maximum priors for m did not result in convergence. We note that differentiation of subdivided populations that exchange few migrants and populations that split and afterwards do not exchange migrants lead to rather similar molecular variation patterns. Therefore, differentiating between migration and drift and temporal population subdivision is difficult. This likely explains the broad posterior distribution of the inferred migration parameters.
Within the Achillea millefolium aggregate, diploid species are usually well separated and their relationships conform to a binary bifurcating tree according to the AFLP data . At the polyploid level, species are difficult to define either by morphology or by molecular data. Previous AFLP data show the polyploid taxa mostly polytomic and polyphyletic [10, 11]. With the haplotype sequence data from both nuclear and plastid genomes available here, we try to infer population history and demography of the diploid species and to untangle the complex relationships among the polyploids.
Gene trees versus species trees of the diploids
The gene trees from the nuclear and the plastid loci of the diploid populations of A. millefolium agg. are incongruent with each other and with the previous AFLP species tree. None of the gene trees corresponds well with the morphological and cytogenetical differentiation of the diploid species, whereas the AFLP tree does (Figures 1, 2 & 4; ). For the ncpGS gene, haplotypes of each of the widespread A. asplenifolia-2x and A. asiatica-2x belong to two groups, IIa and IIb (Figure 1A). Relevant ncpGS sequences indicate intragenic recombination (Additional file 1: S-Figure 1). The SBP gene tree conflicts more severely with the AFLP tree and the species delimitation (Figure 2). Even the plastid sequences show polymorphic haplotypes within each of the three widespread species A. setacea-2x, A. asplenifolia-2x and A. asiatica-2x (Figure 4). As the phylogenetic relationships of these three species indicated by the AFLP tree (Figure 1B) are in line with the morphological and biogeographical information, we suggest that the gene tree incongruence as well as their discordance with the inferred species tree (asserted by the AFLP tree) are due to a lack of sorting of ancestral polymorphic alleles and/or due to introgression after the split of the species. Assuming neutrality, retention of ancestral polymorphism is likely if speciation is fast relative to drift within the populations . The quantitative results from the IM model (Ima2 program ) show large ancestral effective population sizes, short splitting time between them and some migration. This corroborates our inference from the discordant nuclear gene trees and suggests rapid speciation and/or occasional exchange of migrants at the diploid level within A. millefolium agg..
The same pattern apparently also characterizes the more distantly related congeneric species, e.g., A. nobilis-2x, A. ligustica-2x and A. acuminata-2x (Figures 1 & 2). In A. nobilis, a diploid species relatively close to A. millefolim agg., we find intra-locus recombination in the ncpGS gene, indicating its hybrid origin or introgression involving a diploid donor from A. millefolim agg. (Additional file 1: S-Figure 1; Figure 1A).
Phylogeography and rapid speciation
Throughout the N Hemisphere, gene flow between divergent lineages through periods of climate changes has shaped the extant geographical distribution and patterns of genetic variation of many plant species [46–51]. Rapid speciation resulted from post-glacial hybrid contacts and polyploidy during population migration have often been reported and documented. For example, Brysting et al. [34, 52] has untangled the complex history of the polyploid Cerastium alpinum group in the circumpolar arctic flora. Complexity is also evident in the evolutionary radiation of A. millefolium agg [10–12].
The diploid taxa of A. millefolium agg. are limited to Eurasia, following an eco-geographical vicariance pattern but appear to be disjunctive and regressive evidently under the pressure of their more expansive polyploid descendants: A. ceretanica-2x, a relict endemic in subalpine grassland of the E Pyrenees; A. asplenifolia, an endemic relict of humid lowland grasslands in the Pannonian plains from Bulgaria, Hungaria, E Austria and the adjacent Czech Rep.; A. roseoalba-2x, a variable taxon of mesic forest margins and anthropogenous meadows in the geologically recent N Italian plains and foothills with populations in adjacent Switzerland, Austria, and locally even in Bavaria and Slovenia where it is in close contact with A. asplenifolia-2x; A. setacea in steppes from SW Asia and SE Europe to continental areas of C Europe and the Alps; A. latiloba, today a subalpine relict in the mountains of NE Turkey and SW Georgia; A. asiatica-2x in montane to alpine grassland from the Altai to Mongolia and N China; and A. cuspidata-2x, a relict in the W Himalayas.
The previous AFLP data  indicate that A. setacea-2x branches early within the A. millefolium agg. The present sequence data partly show compatible patterns of relationships, i.e., the European species A. ceretanica-2x, A. roseoalba-2x and A. asplenifolia-2x are more closely related to each other and to the Asian A. asiatica-2x than to A. setacea-2x (Figures 1 & 4). The morphology, distribution and habitat preferences suggest that A. roseoalba-2x might have originated from introgression into A. asplenifolia possibly by the A. ceretanica-like populations [15, 53]. Compared to its European sister species, A. asiatica-2x represents a wide geographic extension of A. millefolium agg. into C and E Asia and harbors the richest plastid variation among all the diploid taxa (Figure 4). It could have survived the cold periods of the late Pleistocene in refugia not far from its present occurrence in or near the Asian mountain areas. In spite of the large distributional range, we find little evidence for isolation by distance: The easternmost A. asiatica-2x and the westernmost A. ceretanica-2x share as much genetic variation as each with the geographically intermediate A. asplenifolia-2x and A. roseoalba-2x (Figures 1 & 2).
In contrast to the diploid species that show phylogenetic structuring and are limited to Eurasia, the polyploid taxa of A. millefolium agg. exhibit practically continuous and interrelated relationships and have extended their range to N America. With the available nuclear ncpGS and SPB sequence data alone, it is often hard to decide whether a polyploid taxon studied here is auto- or allopolyploid because even the diploids share polymorphic alleles likely due to incomplete lineage sorting and/or secondary introgression. Only combined with the previous AFLP profiles , the autopolyploid nature of some may become clear, but most are influenced by hybridization.
In Europe, A. ceretanica-4x in central France exhibits some ncpGS correspondence (Figure 3) with A. roseoalba-4x and AFLP affinities with A. ceretanica-2x  but otherwise little affiliation with its diploid cytotype. In populations of A. roseoalba-2x, 4x-individuals occur with a new plastid haplotype (H16 in Figure 5) and corresponding ncpGS alleles (Figure 3). A. styriaca-4x, an endemic from E Austria and the Czech Rep., is ecologically distinct but shares ncpGS and SBP alleles with A. roseoalba-4x and the plastid haplotype H1 with A. asplenifolia-2x and so on. A. setacea-2x and A. asplenifolia-2x have been involved in the origin of the widespread and expansive C to E European allotetraploid A. collina-4x [10, 19]. All the 4x-taxa mentioned above are connected by occasional intermediates, can be easily hybridized in crossing experiments and produce viable progeny with more or less normal meioses [54, 55]. Further to the east in Ukraine, A. inundata-4x is linked to A. setacea-2x and also to A. asplenifolia-2x etc. (Figures 3 & 5).
The diverse 4x-cytotypes from C and E Asia are provisionally called A. asiatica-4x. They are particularly linked to A. asiatica-2x, but also to A. cuspidata-2x (Figure 5; Additional file 3: S-Figure 3). Plastid H1 and H5 also demonstrate links of A. asiatica-4x with the European 2x- and 4x-taxa (Figure 5). How complex the relationships of higher polyploids in A. millefolium agg. are is shown by the recently described Altai endemic A. schmakovii-6x: it shares nuclear haplotypes with A. setacea-2x, A. asplenifolia-2x, A. inundata-4x, A. asiatica-4x and has the rare plastid haplotype H7 from A. asiatica-2x.
The Pleistocene extension of A. millefolium agg. into N America  with populations called A. borealis s.lat. and two main cytotypes (4x and 6x ) is of particular interest. They have developed ecotypes in most habitats from the sand dunes of the Pacific to the peaks of Rocky Mts. and the East Coast forests. Their rapid adaptive differentiation has been well documented by Clausen et al.  and Ramsey et al. . The molecular genetic data from Guo et al. [10, 11] and Ramsey et al.  have shown that all these native N American populations differ from those in Eurasia, are monophyletic and most likely linked to A. asiatica-2x/4x (called A. millefolium var. manshurica Kit. ), but can not be resolved in more details. In the ncpGS and SBP gene trees (Figure 3; Additional file 3: S-Figure 3), they share alleles with A. asiatica-2x and -4x, but also with A. millefolium-6x (A. apiculata) in subarctic Russia and even with A. inundata-4x, A. ceretanica-4x and A. asplenifolia-2x. Considering plastid haplotypes, the most frequent H14 is likely the ancestral, from which H18-H22, 24, 26 have originated (Figure 5). All these and additional fossil data [12, 56, 57] support the assumption that the ancestors of the N American populations might have survived in NE Asian/Alaskan refugia ("Beringia") during the middle or late Pleistocene cold periods [12, 49]. Apart from the westward extension through Siberia to subarctic Europe (A. millefolium-6x/A. apiculata-6x), their main migrations might have been east- and southward into N America. There, without the competition of closely related taxa, they underwent a radiation and formed the 4x and 6x racial complex of A. borealis s.lat.. This can be regarded quasi as a model for the early phase of the eco-geographical radiation of A. millefolium agg. in SE Europe and adjacent SW Asia, and a second phase with A. asiatica-2x+4x, A. alpina-4x and A. schmakovii-6x in C and E Asia [10, 11, 15, 58, 59].
Nuclear and plastid haplotype data analyzed in this study suggest rapid diploid and polyploid speciation in the temperate N Hemisphere Achillea species, especially those within A. millefolium agg.. Hybridization and polyploidy seem to have promoted the recent lineage radiation, shaped the concurrent patterns of genetic variation, and contributed to the wide distribution of this species complex.
The sequence data from two nuclear genes and chloroplast DNA employed in this study result in incongruent trees, obviously due to lack of lineage sorting and/or secondary hybridization, and thus cannot resolve the species phylogeny. This lack of lineage sorting apparently extends to other congeneric taxa. To date, the AFLP tree  is the only molecular tree that can be brought into accordance with species delimitations, morphology and traditional taxonomy of A. millefolium agg.. This is likely due to the averaging effect of the genome-wide sampling of AFLP polymorphism.
With little lineage sorting and frequent gene flow, the species tree can only be recovered using data from many unlinked DNA regions. Despite the development of new techniques such as the RAD tag technology , traditional molecular methods, such as AFLP genome screening, are still useful for non-model species, especially if complemented by likelihood-based Bayesian analytical methods . On the other hand, even if sequence data from the nuclear and chloroplast genomes can not override the AFLP data in resolving recent species divergence, they do help to understand population demography and speciation processes, and to demonstrate that shared ancestral polymorphisms are more common than fixed alleles in young radiating species. Therefore, for the phylogeographic surveys of non-model organisms, we advocate the use of DNA polymorphisms from multiple unlinked loci, e.g., AFLP markers, combined with sequence data from some single genes, as such a combination appears useful and cost and time efficient.
Acknowledgements and funding
This work is supported by the National Natural Science Foundation of China (Grant No. 30770144 and 31121003). We are grateful to E. Temsch, J Li and M. Lambrou for providing technical instructions on DNA ploidy level determinations and chromosome counts. We also thank L. Ehrendorfer-Schratt, A. Tribsch, M. Staudinger, G.-Y. Rao, D.-Y. Tan, and R. Cao for collecting part of the plant samples. Particular thanks are due to J. Ramsey (Rochester, N.Y.) for important additional information on their Achillea research and valuable discussions.
- Grant V: Plant Speciation. 1981, New York, Columbia University PressGoogle Scholar
- Vijverberg K, Kuperus P, Breeuwer JAJ, Bachmann K: Incipient adaptive radiation of New Zealand and Australian Microseris (Asteraceae): an amplified fragment length polymorphism (AFLP) study. J Evol Biol. 2000, 13: 997-1008. 10.1046/j.1420-9101.2000.00241.x.View ArticleGoogle Scholar
- Baldwin BG: A phylogenetic perspective on the origin and evolution of Madiinae. Tarweeds & Silverswords, Evolution of the Madiinae (Asteraceae). Edited by: Carlquist S, Baldwin BG, Carr GD. 2003, St. Louis: Missouri Botanical Garden Press, 193-228.Google Scholar
- Chase MW, Knapp S, Cox AV, Clarkson JJ, Butsko Y, Joseph J, Savolainen V, Parokonny AS: Molecular systematics, GISH and the origin of hybrid taxa in Nicotiana (Solanaceae). Ann Bot. 2003, 92: 107-127. 10.1093/aob/mcg087.PubMedPubMed CentralView ArticleGoogle Scholar
- Soltis DE, Soltis PS, Tate JA: Advances in the study of polyploidy since plant speciation. New Phytol. 2003, 161: 173-191. 10.1046/j.1469-8137.2003.00948.x.View ArticleGoogle Scholar
- Clausen J, Keck D, Hiesey WM: Experimental studies on the nature of species. III. Environmental responses of climatic races of Achillea. Carnegie Institute Washington Publ. 1948, 581-Google Scholar
- Anderson E, Stebbins GL: Hybridization as an evolutionary stimulus. Evolution. 1954, 8: 378-388. 10.2307/2405784.View ArticleGoogle Scholar
- Schemske DW: Population structure and local selection in Impatiens pallida (Balsaminaceae), a selfing annual. Evolution. 1984, 38: 817-832. 10.2307/2408393.View ArticleGoogle Scholar
- Linhart YB, Grant MC: Evolutionary significance of local genetic differentiation in plants. Ann Rev Ecol Syst. 1996, 27: 237-277. 10.1146/annurev.ecolsys.27.1.237.View ArticleGoogle Scholar
- Guo Y-P, Saukel J, Mittermayr R, Ehrendorfer F: AFLP analyses demonstrate genetic divergence, hybridization, and multiple polyploidization in the evolution of Achillea (Asteraceae-Anthemideae). New Phytol. 2005, 166: 273-290. 10.1111/j.1469-8137.2005.01315.x.PubMedView ArticleGoogle Scholar
- Guo Y-P, Saukel J, Ehrendorfer F: AFLP trees versus scatterplots: evolution and phylogeography of the polyploid complex Achillea millefolium agg. (Asteraceae). Taxon. 2008, 57: 153-169.Google Scholar
- Ramsey J, Robertson A, Husband B, Conti E: Rapid adaptive divergence in New World Achillea, an autopolyploid complex of ecological races. Evolution. 2008, 62: 639-653. 10.1111/j.1558-5646.2007.00264.x.PubMedView ArticleGoogle Scholar
- Ramsey J: Polyploidy and ecological adaptation in wild yarrow. Proc Natl Sci USA. 2011, 108: 7096-7101. 10.1073/pnas.1016631108.View ArticleGoogle Scholar
- Ehrendorfer F: Cytology of Achillea hybrids. Carnegie Institute Washington Year Book. 1952, 51: 124-125.Google Scholar
- Ehrendorfer F: Differentiation-hybridization cycles and polyploidy in Achillea. Cold Spring Harbor Symp Quant Biol. 1959, 24: 141-152.PubMedView ArticleGoogle Scholar
- Ehrendorfer F: New chromosome numbers and remarks on the Achillea millefolium polyploid complex in North America. Oesterr Bot Z. 1973, 122: 133-143. 10.1007/BF01624798.View ArticleGoogle Scholar
- Hiesey WM, Nobs MA: Genetic and transplant studies on contrasting species and ecological races of the Achillea millefolium complex. Bot Gaz. 1970, 131: 245-259. 10.1086/336539.View ArticleGoogle Scholar
- Tyrl RJ: Origin and distribution of polyploid Achillea (Compositae) in western North America. Brittonia. 1975, 27: 187-196. 10.2307/2805480.View ArticleGoogle Scholar
- Ma J-X, Li Y-N, Vogl C, Ehrendorfer F, Guo Y-P: Allopolyploid speciation and ongoing backcrossing between diploid progenitor and tetraploid progeny lineages in the Achillea millefolium species complex: Analyses of single-copy nuclear genes and genomic AFLP. BMC Evol Biol. 2010, 10: 100-10.1186/1471-2148-10-100.PubMedPubMed CentralView ArticleGoogle Scholar
- Slotte T, Huang H, Lascoux M, Ceplitis A: Polyploid speciation did not confer instant reproductive isolation in Capsella (Brassicaceae). Mol Biol Evol. 2008, 25: 1472-1481. 10.1093/molbev/msn092.PubMedView ArticleGoogle Scholar
- Strasburg JL, Rieseberg LH: Molecular demographic history of the annual sunflowers Helianthus annuus and H. Petiolaris--large effective population sizes and rates of long-term gene flow. Evolution. 2008, 62: 1936-1950. 10.1111/j.1558-5646.2008.00415.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Slowinski JB, Page RDM: How should species phylogenies be inferred from sequence data?. Syst Biol. 1999, 48: 814-825. 10.1080/106351599260030.PubMedView ArticleGoogle Scholar
- Belfiore NM, Liu L, Moritz C: Multilocus phylogenetics of a rapid radiation in the genus Thomomys (Rodentia: Geomyidae). Syst Biol. 2008, 57: 294-310. 10.1080/10635150802044011.PubMedView ArticleGoogle Scholar
- Edwards SV: Is a new and general theory of molecular systematics emerging?. Evolution. 2009, 63: 1-19. 10.1111/j.1558-5646.2008.00549.x.PubMedView ArticleGoogle Scholar
- Degnan JH, Rosenberg NA: Gene tree discordance, phylogenetic inference and the multispecies coalescent. Trends Ecol Evol. 2009, 24: 332-340. 10.1016/j.tree.2009.01.009.PubMedView ArticleGoogle Scholar
- Hey J: The divergence of chimpanzee species and subspecies as revealed in multipopulation isolation-with-migration analyses. Mol Biol Evol. 2010, 27: 921-933. 10.1093/molbev/msp298.PubMedPubMed CentralView ArticleGoogle Scholar
- Hey J: Isolation with migration models for more than two populations. Mol Biol Evol. 2010, 27: 905-920. 10.1093/molbev/msp296.PubMedPubMed CentralView ArticleGoogle Scholar
- Temsch EM, Greilhuber J: Genome size variation in Arachis hypogaea and A. monticola re-evaluated. Genome. 2000, 43: 449-451.PubMedView ArticleGoogle Scholar
- Suda J, Krahulcová A, Trávnícek P, Krahulec F: Ploidy level versus DNA ploidy level: An appeal for consistent terminology. Taxon. 2006, 55: 447-450. 10.2307/25065591.View ArticleGoogle Scholar
- Doyle JJ, Doyle JL: A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull. 1987, 19: 11-15.Google Scholar
- Emshwiller E, Doyle JJ: Chloroplast-expressed glutamine synthetase (ncpGS): Potential utility for phylogenetic studies with an example from Oxalis (Oxalidaceae). Mol Phylogenet Evol. 1999, 12: 310-319. 10.1006/mpev.1999.0613.PubMedView ArticleGoogle Scholar
- Chapman MA, Chang JC, Weisman D, Kesseli RV, Burke JM: Universal markers for comparative mapping and phylogenetic analysis in the Asteraceae (Compositae). Theor Appl Genet. 2007, 115: 747-755. 10.1007/s00122-007-0605-2.PubMedView ArticleGoogle Scholar
- Shaw J, Lickey EB, Beck JT, Farmer SB, Liu W, Miller J, Siripun KC, Winder CT, Schilling EE, Small RL: The tortoise and the hare II: Relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis. Am J Bot. 2005, 92: 142-166. 10.3732/ajb.92.1.142.PubMedView ArticleGoogle Scholar
- Brysting AK, Oxelman B, Huber KT, Moulton V, Brochmann C: Untangling complex histories of genome mergings in high polyploids. Syst Biol. 2007, 56: 467-476. 10.1080/10635150701424553.PubMedView ArticleGoogle Scholar
- Cronn R, Cedroni M, Haselkorn T, Grover C, Wendel JF: PCR-mediated recombination in amplification products derived from polyploid cotton. Theor Appl Genet. 2002, 104: 482-489. 10.1007/s001220100741.PubMedView ArticleGoogle Scholar
- Wu L, Tang T, Zhou R, Shi S: PCR-mediated recombination of the amplification products of the Hibiscus Tiliaceus cytosolic glyceraldehyde-3-phosphate dehydrogenase gene. J Biochem Mol Biol. 2007, 40: 172-179. 10.5483/BMBRep.2007.40.2.172.PubMedView ArticleGoogle Scholar
- Young ND, Healy J: Gapcoder automates the use of indel characters in phylogenetic analysis. BMC Bioinformatics. 2003, 4: 6-10.1186/1471-2105-4-6.PubMedPubMed CentralView ArticleGoogle Scholar
- Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.PubMedView ArticleGoogle Scholar
- Bandelt HJ, Forster P, Roehl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48.PubMedView ArticleGoogle Scholar
- Clark AG: Neutral behavior of shared polymorphism. Proc Natl Acad Sci USA. 1997, 94: 7730-7734. 10.1073/pnas.94.15.7730.PubMedPubMed CentralView ArticleGoogle Scholar
- Nielsen R, Beaumont MA: Statistical inferences in phylogeography. Mol Ecol. 2009, 18: 1034-1047. 10.1111/j.1365-294X.2008.04059.x.PubMedView ArticleGoogle Scholar
- Hey J, Nielsen R: Multilocus methods for estimating population sizes, migration rates and divergence time, with applications to the divergence of Drosophila pseudoobscura and D. persimilis. Genetics. 2004, 167: 747-760. 10.1534/genetics.103.024182.PubMedPubMed CentralView ArticleGoogle Scholar
- Hudson RR, Kaplan NL: Statistical properties of the number of recombination events in the history of a sample of DNA sequences. Genetics. 1985, 111: 147-164.PubMedPubMed CentralGoogle Scholar
- Kimura M: The number of heterozygous nucleotide sites maintained in a finite population due to steady flux of mutations. Genetics. 1969, 61: 893-903.PubMedPubMed CentralGoogle Scholar
- Hasegawa M, Kishino H, Yano T: Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol. 1985, 22: 160-174. 10.1007/BF02101694.PubMedView ArticleGoogle Scholar
- Comes HP, Kadereit JW: The effect of quaternary climatic changes on plant distribution and evolution. Trends Plant Sci. 1998, 3: 432-438. 10.1016/S1360-1385(98)01327-2.View ArticleGoogle Scholar
- Hewitt GM: Post-glacial re-colonization of European biota. Biol J Linn Soc. 1999, 68: 87-112. 10.1111/j.1095-8312.1999.tb01160.x.View ArticleGoogle Scholar
- Hewitt GM: The genetic legacy of the quaternary ice ages. Nature. 2000, 405: 907-913. 10.1038/35016000.PubMedView ArticleGoogle Scholar
- Hewitt GM: Speciation, hybrid zones and phylogeography--or seeing genes in space and time. Mol Ecol. 2001, 10: 537-549.PubMedView ArticleGoogle Scholar
- Abbott RJ, Smith LC, Milne RI, Crawford RMM, Wolff K, Balfour J: Molecular analysis of plant migration and refugia in the arctic. Science. 2000, 289: 1343-1346. 10.1126/science.289.5483.1343.PubMedView ArticleGoogle Scholar
- Brochmann C, Brysting AK, Alsos IG, Borgen L, Grundt HH, Scheen A-C, Elven R: Polyploidy in arctic plants. Biol Jour Linn Soc. 2004, 82: 521-536. 10.1111/j.1095-8312.2004.00337.x.View ArticleGoogle Scholar
- Brysting AK, Mathiesen C, Marcussen T: Challenges in polyploid phylogenetic reconstruction: A case story from the arctic-alpine Cerastium alpinum complex. TAXON. 2011, 60: 333-347.Google Scholar
- Ehrendorfer F: Achillea roseoalba Ehrendf., spec. nov., eine hybridogene, di- und tetraploide Sippe des Achillea millefolium-Komplexs. Österr Bot Zeitschr. 1959, 106: 363-368. 10.1007/BF01289811.View ArticleGoogle Scholar
- Vetter S, Lambrou M, Franz CH, Ehrendorfer F: Cytogenetics of experimental hybrids within the Achillea millefolium complex (yarrow). Caryologia. 1996, 49: 1-12.View ArticleGoogle Scholar
- Vetter S, Lambrou M, Franz CH, Ehrendorfer F, Saukel J: Chromosome numbers of experimental tetraploid hybrids and selfpollinated progenies within the Achillea millefolium complex (Compositae). Caryologia. 1996, 49: 227-231.View ArticleGoogle Scholar
- Colinvaux PA: The environment of the Bering Land Bridge. Ecol Mono. 1964, 34: 297-329. 10.2307/1948504.View ArticleGoogle Scholar
- Zuzula GD, Froese DG, Westgate JA, La Farge C, Mathewes RW: Paleoecology of Beringian "packrat" middens from central Yukon Territory, Canada. Q Res. 2005, 63: 189-198. 10.1016/j.yqres.2004.11.003.View ArticleGoogle Scholar
- Guo Y-P, Ehrendorfer F, Samuel R: Phylogeny and systematics of Achillea (Asteraceae-Anthemideae) inferred from nrITS and plastid trnL-F DNA sequences. Taxon. 2004, 53: 657-672. 10.2307/4135441.View ArticleGoogle Scholar
- Guo Y-P, Vogl C, Van Loo M, Ehrendorfer F: Hybrid origin and differentiation of two tetraploid Achillea species in East Asia: molecular, morphological and ecogeographical evidence. Mol Ecol. 2006, 15: 133-144.PubMedView ArticleGoogle Scholar
- Emerson KJ, Merz CR, Catchen JM, Hohenlohe PA, Cresko WA, Bradshaw WE, Holzapfel CM: Resolving postglacial phylogeography using high-throughput sequencing. Proc Natl Sci USA. 2010, 107: 16196-16200. 10.1073/pnas.1006538107.View ArticleGoogle Scholar
- Foll M, Fischer MC, Heckel G, Excoffier L: Estimating population structure from AFLP amplification intensity. Mol Ecol. 2010, 19: 4638-4647. 10.1111/j.1365-294X.2010.04820.x.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.