- Research article
- Open Access
The evolutionary history of the Arabidopsis lyrata complex: a hybrid in the amphi-Beringian area closes a large distribution gap and builds up a genetic barrier
© Schmickl et al; licensee BioMed Central Ltd. 2010
- Received: 27 November 2009
- Accepted: 8 April 2010
- Published: 8 April 2010
The genomes of higher plants are, on the majority, polyploid, and hybridisation is more frequent in plants than in animals. Both polyploidisation and hybridisation contribute to increased variability within species, and may transfer adaptations between species in a changing environment. Studying these aspects of evolution within a diversified species complex could help to clarify overall spatial and temporal patterns of plant speciation. The Arabidopsis lyrata complex, which is closely related to the model plant Arabidopsis thaliana, is a perennial, outcrossing, herbaceous species complex with a circumpolar distribution in the Northern Hemisphere as well as a disjunct Central European distribution in relictual habitats. This species complex comprises three species and four subspecies, mainly diploids but also several tetraploids, including one natural hybrid. The complex is ecologically, but not fully geographically, separated from members of the closely related species complex of Arabidopsis halleri, and the evolutionary histories of both species compexes have largely been influenced by Pleistocene climate oscillations.
Using DNA sequence data from the nuclear encoded cytosolic phosphoglucoisomerase and Internal Transcribed Spacers 1 and 2 of the ribosomal DNA, as well as the trnL/F region from the chloroplast genome, we unravelled the phylogeography of the various taxonomic units of the A. lyrata complex. We demonstrate the existence of two major gene pools in Central Europe and Northern America. These two major gene pools are constructed from different taxonomic units. We also confirmed that A. kamchatica is the allotetraploid hybrid between A. lyrata and A. halleri, occupying the amphi-Beringian area in Eastern Asia and Northern America. This species closes the large distribution gap of the various other A. lyrata segregates. Furthermore, we revealed a threefold independent allopolyploid origin of this hybrid species in Japan, China, and Kamchatka.
Unglaciated parts of the Eastern Austrian Alps and arctic Eurasia, including Beringia, served as major glacial refugia of the Eurasian A. lyrata lineage, whereas A. halleri and its various subspecies probably survived in refuges in Central Europe and Eastern Asia with a large distribution gap in between. The North American A. lyrata lineage probably survived the glaciation in the southeast of North America. The dramatic climatic changes during glaciation and deglaciation cycles promoted not only secondary contact and formation of the allopolyploid hybrid A. kamchatica, but also provided the environment that allowed this species to fill a large geographic gap separating the two genetically different A. lyrata lineages from Eurasia and North America. With our example focusing on the evolutionary history of the A. lyrata species complex, we add substantial information to a broad evolutionary framework for future investigations within this emerging model system in molecular and evolutionary biology.
- Internal Transcribe Space
- Species Complex
- trnL Intron
- Parental Taxon
- Eurasian Lineage
Molecular biological research during the last decade has largely focussed on model organisms such as Drosophila melanogaster, Caenorhabditis elegans, and Arabidopsis thaliana. Now that knowledge in molecular genetics, cell and developmental biology of these organisms has greatly increased, closely related organisms emerge as promising for studying characteristics not possible to elucidate with model and/or single organisms [1, 2].
Arabidopsis lyrata L. is a close relative of A. thaliana, from which it diverged approximately five million years ago [3, 4]. Arabidopsis lyrata s.l. represents a small species complex of four species and several putative subspecies with a circumpolar arctic-alpine distribution (Additional file 1, Table S1). Populations have been adapted to various ecological conditions, including the harsh environment of the arctic tundra, cryptic warm-stage refugia (exposed rocks, rocky slopes) in Central Europe, and different edaphic conditions with substrates such as dolomite, silicious bedrocks, and even heavy metal rich serpentine soil in Central Europe (Lower Austria, personal observation) and the USA [Maryland; ] . Most members of the species complex are perennial diploid outbreeders (2n = 2x = 16), but also tetraploid cytotypes occur [7, 6, 8]. There are numerous aspects of the biology of the A. lyrata complex that differ from and cannot be addressed in A. thaliana, like self-incompatibility and perennial life cycle.
The A. lyrata complex has already proven to be a suitable study system for the analysis of character traits such as flowering time [9, 10] or pathogen defense . Additionally, molecular mechanisms for the function of sporophytic self-incompatibility have been investigated [12–20], and comparative approaches to analyse sporophytic self-incompatibility in diploids versus polyploids are underway (Jørgensen, unpublished data). Whole genome sequencing of A. lyrata was finished last year, and data have been available for a few months (The A. lyrata genome sequence assembly v1.0, http://genome.jgi-psf.org/Araly1/Araly1.info.html), enabling direct comparisons with the A. thaliana genome. However, in contrast to A. thaliana, where the evolutionary history has been analysed in more detail [e.g. [21, 22]], evolutionary studies on the A. lyrata complex have so far been largely restricted to a small number of populations from Central Europe [23–26] or larger sample sizes with a more general genus-wide perspective [27, 28]. A detailed European study has revealed that population structure is dominated by regional genetic bottlenecks, and genetic structure exists within continents [Ansell, personal communication]. This suggests a comprehensive global study is necessary to resolve the evolutionary history of this complex.
In this study we present the first worldwide evolutionary history of the A. lyrata complex, covering its whole range of distribution and all taxonomically defined units. We use a widely applied nuclear encoded marker system (ITS, internal transcribed spacer region of nuclear encoded ribosomal DNA) to study gene flow between populations, a maternally inherited chloroplast genome marker (trnL intron (trnL) and trnL/F intergenic spacer (trnL/F-IGS) of tRNASer and tRNAThr, respectively) to investigate migrational movements due to seed dispersal, and the nuclear encoded housekeeping gene PgiC (cytosolic phosphoglucoisomerase), a single copy gene, to discriminate between hybridising taxa. We aim to focus on the following four aspects: (1) Unravelling general phylogeographic patterns of the A. lyrata complex by identifying the main genetic lineages, and interpreting genetic variation in space and time , in the context of both climatic and geological events throughout Pleistocene glaciation cycles; (2) Evaluating the role of hybridisation and polyploidisation in the origin of A. kamchatica (Fisch. ex DC.) K. Shimizu & Kudoh, an amphi-Beringian member of the A. lyrata complex; (3) Explaining Pleistocene and postglacial migration routes by analysing genetic diversity statistics: The arctic-alpine A. lyrata complex is one of the rare examples among higher plants with a distribution in both Central Europe and North America and, additionally, a circumpolar distribution - other examples are Cassiope tetragona , and Saxifraga oppositifolia ; and (4) Studying the role of Beringia as a refugia for populations of arctic A. lyrata, since Beringia is assumed to be one of the major refugia for arctic plants during Pleistocene glaciations [30–36].
The following short overview of A. lyrata and A. halleri taxonomy will introduce the concept of Al-Shehbaz and O'Kane . Elven (ed.)  persues a different taxonomic concept, which is summarised in the supplementary material (Additional file 1, Table S1), but will not be discussed here. Apart from A. thaliana, the A. lyrata complex is one of three major species complexes in the genus Arabidopsis , the other two being the A. arenosa (L.) Lawalrée and A. halleri complexes [39, 40, 6, 27, 41]. The A. lyrata complex is considered by different authors to include different numbers of taxa of various ranking, distribution areas, and ploidy levels (Additional file 1, Table S1). Al-Shehbaz and O'Kane  treat the complex as one species, A. lyrata, with three subspecies: ssp. lyrata, ssp. petraea (L.) O'Kane & Al-Shehbaz, and ssp. kamchatica (Fisch. ex DC.) O'Kane & Al-Shehbaz. Subspecies lyrata is considered to be broadly amphi-Pacific and of two ploidy levels (2n = 16, 32), ssp. petraea to be Northern Eurasian and Central European with the same two ploidy levels, and kamchatica to be amphi-Pacific and tetraploid (2n = 32). A comprehensive quantitative morphological analysis using multivariate statistics is still not available for the A. lyrata complex. Moreover, information on the ploidy level is largely lacking throughout the distribution range, and further cytological investigations might be important for taxon delimitation [7, 8]. Arabidopsis arenicola (Richardson) Al-Shehbaz, Elven, D.F. Murray & Warwick has only recently been included as part of the A. lyrata complex . It has, for a long time, been placed within the genus Arabis L. Elven (ed.)  considers this taxon to be diploid (2n = 16) and distributed in north-eastern North America, and they are supported by Al-Shehbaz in the upcoming Flora of North America [Al-Shehbaz, personal communication]. Many of the mentioned taxa include two ploidy levels, suggesting frequent polyploidisation events within the A. lyrata complex. An allopolyploid origin of kamchatica has already been confirmed based on nuclear DNA sequences, with A. lyrata and A. halleri ssp. gemmifera (Matsum.) O'Kane & Al-Shehbaz as possible parental taxa [[43, 27, 28, 8]; Jørgensen et al., unpublished data]. Otherwise, little is known with regard to the number of polyploid units and their origins.
Five subspecies have been recognised in A. halleri: ssp. halleri, ssp. ovirensis (Wulfen) O'Kane & Al-Shehbaz, ssp. dacica (Heuff.) Kolník, comb. nov., ssp. tatrica (Pawl.) Kolník, comb. nov., all distributed in Central Europe, and ssp. gemmifera in Eastern Asia, supported by both morphometric analysis (Kolnik, unpublished data) and genetic AFLP data (Marhold, unpublished data).
DNA isolation, amplification and sequencing
Total DNA was obtained from dried leaf material and extracted according to the CTAB protocol of Doyle and Doyle  with the following modifications: 50-75 mg of dry leaf tissue were ground in 2 ml tubes using a Retsch swing mill (MM 200), 2 units of RNase A per extraction were added to the isolation buffer, and the DNA pellets were washed twice with 70% ethanol. DNA was dissolved in 50 μl TE-buffer for storage and diluted 1:3 in TE-buffer before use.
For the cpDNA markers trnL intron and trnL/F intergenic spacer (trnL/F-IGS), primers and PCR cycling scheme followed the protocol of Dobeš et al. , using a PTC200 (MJ Research, Waltham, USA) thermal cycler. The PCR reaction volume of 50 μl contained 1× PCR buffer (10 mM TRIS/50 mM KCl buffer, pH 8.0), 3 mM MgCl2, 0.4 μM of each primer, 0.2 mM of each dNTP, 1 U Taq DNA polymerase (Amersham Biosciences, Chalfont St Giles, England), and approximately 1 ng of template DNA. Amplified sequences of trnL/F-IGS included the complete trnL/F-IGS and the first 18 bases of the trnF gene. Amplification of the nuclear marker internal transcribed spacer region (ITS) was performed according to Dobeš et al. . PCR reaction conditions were the same as for the two cpDNA markers described above, and PCR cycling scheme was 5 min at 95°C, 35 cycles of 1 min at 95°C, 1 min at 48°C, and 1 min at 72°C, 10 min extension at 72°C, and a final hold at 4°C. PCR products spanned the entire ITS1, 5.8 S rDNA, and ITS2 region.
Before sequencing PCR products were checked for length and concentrations on 1.5% agarose gels and purified with the NucleoFast Kit (Macherey-Nagel, Düren, Germany). Cycle sequencing was performed using the DYEnamic ET Terminator Cycle Sequencing Kit (Amersham Biosciences, Chalfont St Giles, England) and the original primers. Both strands were amplified in order to gain the complete sequence. PCR products were resolved in 10 μl loading solution and run on a MegaBace 500 sequencer.
trnL/F and ITS sequence definition and map reconstruction
Plastidic trnL/F sequences were defined as haplotypes and suprahaplotypes following our previous studies [27, 41, 28]: Haplotypes are characterised by multiple trnF pseudogenes in the 3'-region of the trnL/F-IGS close to the functional trnF gene [47, 48, 27, 49, 41, 50]. When defining trnL/F suprahaplotypes, we excluded the pseudogene-rich region. Pseudogenes evolve with a 10 times higher mutation rate than single nucleotide polymorphisms, which makes them non-applicable for phylogenetic reconstructions on the species level . Additionally, their boundaries are under ongoing discussion. An alternative interpretation is provided by Ansell et al. . In summary, haplotypes belonging to one suprahaplotype share the same base order throughout the whole sequence except for the pseudogene-rich region, where they vary in both length and base content. Suprahaplotypes differ from each other only by single point mutations and/or indels. Newly defined trnL/F haplotypes were assigned to GenBank numbers [GQ922894-GQ922903] (Additional file 2, Table S2).
ITS sequences were obtained from direct sequencing of PCR products in order to detect hybrid individuals according to ambiguous base positions. Sequences were defined as ITS types and supratypes following our previous studies [27, 41, 28]: Most ITS types contained ambiguous sites as a result of the direct sequencing approach. They were combined to ITS supratypes by replacing the ambiguous sites by the bases with higher fluorescence intensity in the electropherogramm. Sequences with equal fluorescence intensity of the two bases at the ambiguous positions were only found between ITS supratypes b and e and labelled b/e ambiguous. Koch and Matschinger  already showed that analyses with either ITS types or ITS supratypes produce the same results. However, working with a limited number of ITS supratypes in contrast to a vast number of ITS types contributes to a clearer data display. Multiple ITS types within a single individual indicate the natural variability among ITS loci within an individual as result of either mutations or gene flow between individuals. In these individuals concerted evolution, which denotes the process of DNA sequence homogenisation among members of multigene families by gene conversion and/or unequal crossing over [52–55], is not yet completed. Newly defined ITS types were assigned to GenBank numbers [GQ922904-GQ922910] (Additional file 2, Table S2).
Maps were constructed using BioOffice version 2.0.6 to create shapefiles and drawn with ArcView version 8.2. Shapefiles for visualising the maximum extent of the ice sheets during the LGM (Last Glacial Maximum) were provided by Ehlers and Gibbard .
Network analyses and genetic diversity statistics
Network analyses and genetic diversity statistics were exclusively performed using the trnL/F suprahaplotypes, as the pseudogene-rich region is not applicable for phylogenetic reconstructions on the species level . The network was constructed using TCS version 1.21  using the statistical parsimony algorithm . Gaps (except polyT stretches) were coded as single additional binary characters. Newly characterised suprahaplotypes were added to the trnL/F network of the A. lyrata complex published earlier . Genetic diversity statistics were performed with Arlequin version 3.11 . Pairwise genetic differentiation was calculated among the following nine taxonomic and regional groups: A. lyrata ssp. petraea from (1) unglaciated Central Europe, (2) glaciated northern Europe, (3) northern Russia and western Beringia, (4) eastern Beringia; A. lyrata ssp. lyrata from (5) unglaciated North America and the glaciated Great Lakes region; (6) Arabidopsis arenicola from glaciated North America and Greenland; and A. kamchatica from (7) Japan, (8) north-eastern Russia, (9) Alaska and western Canada. FST values, regarding haplotype frequencies only, and ΦST values, which take into account the genetic relationship between haplotypes as pairwise character differences, were calculated. The significance of differentiation was examined using a permutation test with 1000 permutations. Additionally, genetic diversity was estimated as nucleotide diversity π.
ITS parsimony analysis
ITS data were analysed based on a 652 bp alignment of ITS supratypes . The total number of variable sites was 44, with 16 of them parsimony informative. From altogether 10 supratypes of both the A. lyrata complex and A. halleri ssp. gemmifera, a strict consensus tree was constructed using maximum parsimony with MEGA version 4.1 . Heuristic searches were performed with 10 random addition sequences and Closest Neighbour Interchange (CNI) branch swapping. Bootstrap values were calculated based on 500 replicates. Length of the most parsimonious trees was 48 mutational steps with a consistency index (CI) = 0.89 and a retention index (RI) = 0.96 (autapomorphies excluded). Arabidopsis thaliana was used as outgroup.
Primer design for the nuclear marker PgiC
Various higher plants are known to have a duplicated locus of the cytosolic enzyme phosphoglucoisomerase [61–63], and the loci are normally unlinked. Also within the genus Arabidopsis, extensive sequencing of the PgiC locus using a cloning strategy revealed a duplication, both in the A. lyrata and A. halleri complex, and this duplication must have predated the evolutionary split between these two species complexes (Additional file 3, Figure S1; Jørgensen, unpublished data). Both loci were initially simultaneously amplified with the general forward primer 5'-TGCTGTSAGCACTAATCTTGCG-3' and the general reverse primer 5'-TCGAACCCGGGAGAGGTAGACCA-3', following the protocol of Wright et al. . The resulting sequence data showed that a group of alleles at the PgiC1 locus were exclusively found in A. halleri and the allopolyploid A. kamchatica. It was thus possible to design A. halleri-specific primers that worked as a high-throughput and simple PCR-based screening marker to discriminate between genomes of the A. lyrata and A. halleri complex. Unfortunately, it was not possible to develop a PCR-based reciprocal marker system characterising alleles from the A. lyrata gene pool because of a lack of appropriate DNA sequence variation. In general, the alleles from both duplicated PgiC loci were only weakly differentiated among and within species. However, within the PgiC1 locus we found a deletion of 7 bp length in the A. lyrata complex compared to the A. halleri complex (Additional file 4, Figure S2). Both groups of alleles were also substantially differentiated by various SNPs. A primer pair with the forward primer located partly within this 7 bp indel (5'-CATTCAACAGATTGTG-3') and the reverse primer 5'-CCAGTAAACATCATGT-3' was developed to amplify a 92 bp fragment within the PgiC1 locus (Additional file 5, Figure S3). The PCR reaction volume of 50 μl contained 1× PCR buffer (10 mM TRIS/50 mM KCl buffer, pH 8.0), 2.5 mM MgCl2, 0.13 μM of each primer, 0.2 mM of each dNTP, 1 U Taq DNA polymerase (Amersham Biosciences, Chalfont St Giles, England), and approximately 1 ng of template DNA. The PCR cycling scheme was 3 min at 94°C, 35 cycles of 20 sec at 94°C, 30 sec at 56°C, and 20 sec at 68°C, 20 sec extension at 68°C, and a final hold at 4°C. By screening the absence/presence of the PCR product of this taxon-specific primer pair, we were able to follow the genetic footprint of the A. halleri complex in its allopolyploid hybrids throughout its distribution range.
Chloroplast sequence data indicate three main genetic lineages: Eurasia, North America, and the amphi-Pacific region
Although all three lineages were characterised both by lineage-specific central and "tip" suprahaplotypes, sharing of central suprahaplotypes was observed, e.g. suprahaplotypes A and B were detected in a few accessions of the Eurasian lineage (Figure 2, Additional file 6, Table S3). This finding is congruent with the observation of central suprahaplotype sharing between the three main species complexes of the genus, A. lyrata, A. halleri, and A. arenosa . This observation has been explained by ancestral cpDNA polymorphism predating the radiation of the genus approximately two million years ago .
Cytosolic phosphoglucose isomerase identifies Arabidopsis halleri ssp. gemmifera as one parent of the allopolyploid amphi-Pacific Arabidopsis kamchatica
Gene diversity statistics show highest genetic diversity in the Eurasian lineage, strongly reduced diversity in the North American lineage, and extremely low diversity in the allopolyploid amphi-Pacific lineage
Taxonomic and regional genetic differentiation based on cpDNA suprahaplotypes.
A. lyrata ssp. petraea: Unglaciated Central Europe
0.603 +/- 0.329
A. lyrata ssp. petraea: Glaciated N Europe
0.764 +/- 0.425
A. lyrata ssp. petraea: N Russia, W Beringia
0.774 +/- 0.422
A. lyrata ssp. petraea: E Beringia
0.339 +/- 0.236
A. lyrata ssp. lyrata: Unglaciated N America, glaciated Great Lakes region
0.080 +/- 0.072
A. arenicola: Glaciated N America, Greenland
0.031 +/- 0.042
A. kamchatica: Japan
0.000 +/- 0.000
A. kamchatica: Far Eastern Federal District of Russia
0.000 +/- 0.000
A. kamchatica: Alaska, W Canada
0.007 +/- 0.019
Pairwise genetic differentiation (FST and ΦST) among taxonomic and regional groups.
A. lyr. ssp. petraea: Unglac. Central Europe
A. lyr. ssp. petraea: Glac.
A. lyr. ssp. petraea:
N Russia, W Beringia
A. lyr. ssp. petraea:
A. lyr. ssp. lyrata: Unglac.
N America, glac.
Great Lakes region
A. arenicola: Glac.
N America, Greenland
A. kamch.: Japan
A. kamch.: Far Eastern Federal District of Russia
A. kamch.: Alaska,
A. lyr. ssp. petraea: Unglac. Central Europe
A. lyr. ssp. petraea: Glac.
A. lyr. ssp. petraea:
N Russia, W Beringia
A. lyr. ssp. petraea:
A. lyr. ssp. lyrata: Unglac.
N America, glac.
Great Lakes region
N America, Greenland
A. kamch.: Japan
A. kamch.: Far Eastern Federal District of Russia
A. kamch.: Alaska,
In contrast to the Eurasian lineage of A. lyrata ssp. petraea, the North American lineage showed an approximately tenfold reduction in nucleotide diversity and strong differentiation according to pairwise ΦST and FST (Table 2). Nucleotide diversity of accessions from predominantly unglaciated southeastern North America (π = 0.0008) was higher than from North America and Greenland (π = 0.0003) (Table 1), which had been under the Laurentide ice sheet during the LGM, possibly indicating genetic bottlenecks with subsequent rapid postglacial immigration.
Extremely reduced nucleotide diversity was reported from amphi-Pacific A. kamchatica (π = 0.0000) (Table 1). When all three genetic groups of A. kamchatica were treated separately, only a single trnL/F suprahaplotype was found in each group (Japan: AD; China: C; Far Eastern Federal District of Russia/Alaska/western Canada: B) (Additional file 6, Table S3). Additionally, haplotype diversity was low, especially in A. kamchatica with suprahaplotype B. Only one haplotype (no. 84) was detected over a vast amphi-Pacific area from Kamchatka to western Canada (Additional file 7, Table S4).
Refugia as areas of secondary contact of formerly allopatric populations: Beringia as an example
High genetic diversity of the Arabidopsis lyratacomplex in Eurasia - postglacial migration from Central European and northern Russianrefugia
For arctic-alpine taxa, centres of species and genetic diversity are, in most cases, considered to concur with Pleistocene refugia [64, 66, 67]. In the genus Arabidopsis, both the number of accepted taxa, trnL/F suprahaplotype and ITS supratype diversity are highest in Central Europe, indicating this area as a centre of diversity [37, 27]. Especially the Eastern Alps could have served as a refugium of A. lyrata ssp. petraea, since both diploid and cytogenetically stabilised tetraploid populations occur there (Schmickl and Koch, unpublished data). For the A. lyrata complex, however, Koch and Matschinger  concluded that Central Europe is not the only centre of diversity. Our results support this conclusion, showing that the Eurasian lineage of the A. lyrata complex is genetically diverse in both Central Europe and arctic Eurasia, including Beringia. To our knowledge such a pattern has not been observed for any other arctic plant with additional distribution in Central Europe. Either genetic diversity is high in the Arctic compared to the Alps, as in Saxifraga oppositifolia L., suggesting long-term evolution in Beringia and more recent colonisation of the Alps [33, 31], or genetic diversity is highest in the Alps and decreasing towards the Arctic, as observed in Arabis alpina L. [68, 69], indicating recent and rapid colonisation from Central Europe. Recent migration into the Arctic, associated with a loss of genetic diversity due to repeated bottlenecks, was also found in Ranunculus glacialis L.  and Rubus chamaemorus L. . The overall high genetic diversity in the Eurasian A. lyrata lineage can probably be explained by long-term glacial survival in multiple refugia (Central Europe, northern Russia, western and eastern Beringia).
So far, periglacial survival of A. lyrata was reported only from unglaciated Central Europe, based on microsatellite data  and cpDNA sequences . This study showed high genetic diversity also in formerly glaciated northern Europe, caused by numerous unique, locally distributed trnL/F suprahaplotypes. A possible explanation might be geographic isolation of populations, which either periglacially survived along the coastline or postglacially migrated into northern Europe, as reported from Ansell [personal communication]. Although several authors suggested periglacial survival on nunataks along the Norwegian coastline, mainly based on geomorphological investigations (reviewed by Brochmann et al. ), we assume postglacial colonisation of formerly glaciated northern Europe from both unglaciated Central Europe and northwestern Russia. Subsequent geographic isolation of populations during Holocene warming probably led to the fixation of local sequence types. Although in our study the major trnL/F suprahaplotypes C and AG were similarly frequent both in Central and northern Europe, a comparative microsatellite study between populations from Central Europe and Iceland revealed significant differences in marker polymorphism . Riihimäki and Savolainen  even found divergent Central and northern European physiological morphotypes, with earlier and more frequent flowering plants in the south, and later and rarer flowering plants in the north.
We observed high genetic diversity within the Eurasian A. lyrata lineage, but also some genetic homogenisation. The trnL/F suprahaplotypes C and AG were found all over the distribution range, which could be explained by repeated gene flow between populations during glacial periods, additionally facilitated by rapid long-distance dispersal of the small Arabidopsis seeds across the smooth snow surface of the tundra and tundra-steppe by strong winds. This mode of long-distance dispersal could have bridged distances up to 2000 km  or even 4000 km (Westergaard et al., in prep.) and seems to be frequent in the Arctic [34, 73].
Ancient split of the Eurasian and North American lineage
The strong genetic differentiation we observed between the Eurasian and North American lineages is probably the result of long-term geographic isolation during Pleistocene glaciations. Muller et al.  detected genetic divergence between North American and western European populations in a comparative microsatellite study, and found that the former had a lower diversity. However, they focused on a few populations only and included only one North American population. The strong reduction of genetic diversity we observed in the North American lineage compared to the Eurasian lineage is congruent with nuclear and plastidic marker data of Wright et al. , but contradicting Balañα-Alcaide et al. , who reported similar genetic diversity based on two nuclear markers. However, both studies included only two populations from each of North America and Europe, and they may not have been representative. Because of the strongly reduced genetic diversity of the North American lineage, we assume this lineage to be derived from the Eurasian lineage. Future studies of more southern populations from the USA will test this assumption. Migration into North America apparently was associated with a strong genetic bottleneck, which restricted populations to ITS supratype e and trnL/F suprahaplotype A. Colonisation of North America is possible from different directions, either from Europe, or from East Asia. Both possibilities have been discussed for various circumboreal species (reviewed by Abbott and Brochmann ). In the case of A. lyrata, colonisation of North America from Russia seems, however, most likely, as the North American genotype was also rarely detected in western Beringia, but not at all in Iceland or Scandinavia.
An amphi-Beringian Arabidopsishybrid zone - due to allopolyploid success?
The majority of arctic polyploids have a history as postglacial colonisers . It has frequently been assumed that polyploids have a broader adaptive potential for recolonising formerly glaciated areas [75, 76]. As one of the major refugia for arctic plants, Beringia has a higher proportion of diploids compared to areas glaciated during Pleistocene . New adaptations in polyploids may evolve either by genome rearrangements  and/or epigenetic changes [78–80] within the first generations after polyploid formation, as reported from rapid gene silencing in the allopolyploid A. suecica (Fr.) Norrl. ex O.E. Schulz . Hybridisation is frequently involved in polyploidisation, leading to the formation of allopolyploids with one each of the parental genomes. Allopolyploidisations have been reported for several arctic species complexes, such as the high polyploid Cerastium alpinum L. complex , high polyploid Primula sect. Aleuritia [82–84], tetra-/hexaploid Silene L. [85, 86], tetra- to octoploid Saxifraga section Mesogyne Sternb. [87, 88], and tetra- to dodecaploid Cardamine digitata Richardson . Introgression, the integration of genetic material from one species into another through repeated backcrossing, was observed between polyploid Saxifraga cernua L. and diploid Saxifraga sibirica L. . The first known Arabidopsis allopolyploid was A. suecica with the maternal parent A. thaliana and the paternal parent A. arenosa, confirmed by artificial crosses . This species probably developed around 20 000 years ago  or between 20 000 and 300 000 years ago , with a single origin in Fennoscandia [91, 92]. The distribution range of this mainly outcrossing  allopolyploid species is rather small. The second natural Arabidopsis allopolyploid, A. kamchatica, has A. halleri ssp. gemmifera and a member of the A. lyrata complex as parental taxa [[43, 27, 8]; Jørgensen, unpublished data]. This allopolyploid origin could be confirmed not only for Japanese, but for all A. kamchatica accessions. According to chloroplast trnL/F data, three different genetic groups were found, which are geographically isolated from each other: (1) accessions of a widespread distribution range from Kamchatka, western and eastern Beringia, to pacific western Canada (trnL/F suprahaplotype B), (2) Japanese accessions (trnL/F suprahaplotype AD), and (3) accessions from pacific eastern China (trnL/F suprahaplotype C). In a more detailed comparison of chloroplast trnL/F and nuclear encoded ITS sequence data, two directions of gene flow could be observed: Either the paternal genome originated from a member of the A. lyrata species complex (ITS supratype b) and the maternal genome from A. halleri ssp. gemmifera (trnL/F suprahaplotype AD), as already reported for A. kamchatica from Japan . Or A. halleri ssp. gemmifera represented the paternal genome (ITS supratype z, derived from ITS supratype r exclusively found in A. halleri ssp. gemmifera), and a member of the A. lyrata species complex served as donor of the maternal genome (trnL/F suprahaplotype C), as in A. kamchatica from China. The North American lineage of the A. lyrata complex could be excluded as a parent, as neither ITS supratype e nor trnL/F suprahaplotype A were detected in A. kamchatica. According to these data we suggest at least three independent origins of A. kamchatica, first with maternal A. halleri ssp. gemmifera in Japan, second with paternal A. halleri ssp. gemmifera in China, and third with an unknown direction of gene flow in Kamchatka, but in all cases with a member of the Eurasian A. lyrata lineage as hybridisation partner.
The most profound change in A. kamchatica in contrast to its parental species is the switch from outcrossing (with sporophytic self-incompatibility system) to selfing [ssp. kamchatica: [15, 8]; ssp. kawasakiana: ]. Selfing increases the possibility of rapid range expansion, as a population can arise from a single individual independent of pollinators and pollen donors. Such a switch is already well known from A. thaliana, dated from around 413 000  to one million years ago , and is therefore not necessarily correlated with hybridisation and polyploidisation. Change of mating systems is one of the major driving forces for speciation, initiating reproductive isolation of populations [97, 98]. However, it is still unclear if there is a correlation between mating systems and hybridisation and/or polyploidisation. The breakdown of self-incompatibility in an artificial cross between A. thaliana and A. lyrata  could indicate a correlation between hybridisation and a switch in mating system, and it could further indicate possible dominance of the selfer A. thaliana over A. lyrata. Otherwise, breakdown of the sporophytic self-incompatibility system has been reported mainly from diploid individuals [16, 100].
In contrast to A. suecica, allopolyploid A. kamchatica has a vast distribution range spanning the whole amphi-Beringian region. Our data suggest that this large distribution is partly due to postglacial colonisation of formerly glaciated areas in eastern Beringia and western Canada. The extremely reduced genetic diversity, particularly of the group with suprahaplotype B, suggests that postglacial immigration may have been rapid, possibly facilitated by de novo adaptations as a result of hybridisation and polyploidisation. Moreover, the success of A. kamchatica as a rapid coloniser may have been enhanced by the availability of large, open landscapes, where habitats were frequently disturbed by glacial and/or permafrost activity. However, the change in mating systems may have had a strong impact on the success of A. kamchatica as a postglacial coloniser, and also on the establishment of the genetic barrier between A. kamchatica and the Eurasian and North American lineages.
Beringia as contact zone between the Eurasian and North American lineage of the Arabidopsis lyratacomplex
Beringia served as a glacial refugium for numerous arctic plant taxa such as Dryas integrifolia Vahl. , Saxifraga hirculus L. , Saxifraga oppositifolia [33, 31], and Vaccinium uliginosum L. [35, 101]. Periglacial survival in Beringia can also be assumed for the A. lyrata complex, in particular, for the Eurasian lineage in arctic western and eastern Beringia north of Brooks Range. The nuclear sequence data indicate that gene exchange with populations of the North American lineage (A. lyrata ssp. lyrata, A. arenicola) occurred inter- and postglacially. However, the data do not support glacial survival of the North American lineage in Beringia, as no plastidic sequence types of this lineage were found in Beringia.
By presenting a worldwide evolutionary history of the Arabidopsis lyrata species complex, we provide solid knowledge about centres of genetic diversity, different genetic lineages, their contact zones, and hybrid speciation. We could clearly differentiate three genetic lineages, a Eurasian, a North American, and an amphi-Pacific one. The latter is constituted of the allopolyploid A. kamchatica, a hybrid between A. lyrata and A. halleri. Further investigations of the population dynamics and the role of selfing within this hybrid species should be conducted to gain a deeper understanding of hybrid establishment in the wild.
We thank Ishan Al-Shehbaz (Missouri, USA), Galina Gusarova (Oslo, Norway), Gu Hongya (Beijing, P. R. China), Barbara Mable (Glasgow, Scotland), David L. Remington (Greensboro, USA), Outi Savolainen (Oulu, Finland), and the curators of the Herbariums of the Natural History Museums London and Vienna for providing plant material, Susanne Ball and Stefan Wötzel for laboratory assistance, and Jürgen Ehlers and Phil Gibbard for shapefiles of the LGM. This research was supported by DFG grant KO 2302/5-2 to Marcus A. Koch.
- Mitchell-Olds T: Arabidopsis thaliana and its wild relatives: a model system for ecology and evolution. Trends Ecol Evol. 2001, 16: 693-700. 10.1016/S0169-5347(01)02291-1.Google Scholar
- Benfey PN, Mitchell-Olds T: From genotype to phenotype: systems biology meets natural variation. Science. 2008, 320: 495-497. 10.1126/science.1153716.PubMed CentralPubMedGoogle Scholar
- Koch MA, Haubold B, Mitchell-Olds T: Comparative evolutionary analysis of chalcone synthase and alcohol dehydrogenase loci in Arabidopsis, Arabis, and related genera (Brassicaceae). Mol Biol Evol. 2000, 17: 1483-1498.PubMedGoogle Scholar
- Koch MA, Haubold B, Mitchell-Olds T: Molecular systematics of the Brassicaceae: evidence from coding plastidic MATK and nuclear CHS sequences. Am J Bot. 2001, 88: 534-544. 10.2307/2657117.PubMedGoogle Scholar
- Turner TL, Von Wettberg EJ, Nuzhdin SV: Genomic analysis of differentiation between soil types reveals candidate genes for local adaptation in Arabidopsis lyrata. PLoS One. 2008, 3: e3183-10.1371/journal.pone.0003183.PubMed CentralPubMedGoogle Scholar
- Clauss MJ, Koch MA: Poorly known relatives of Arabidopsis thaliana. Trends Plant Sci. 2006, 11: 449-459. 10.1016/j.tplants.2006.07.005.PubMedGoogle Scholar
- Zhukova PG, Petrovsky VV, Plieva TV: The chromosome numbers and taxonomy of some plant species from Siberia and Far East. Botanicheskiy Zhurnal. 1973, 58: 1331-1342.Google Scholar
- Shimizu-Inatsugi R, Lihová J, Iwanaga H, Kudoh H, Marhold K, Savolainen O, Watanabe K, Yakubov VV, Shimizu KK: The allopolyploid Arabidopsis kamchatica originated from multiple individuals of Arabidopsis lyrata and Arabidopsis halleri. Mol Ecol. 2009, 18: 4024-4048. 10.1111/j.1365-294X.2009.04329.x.PubMedGoogle Scholar
- Riihimäki M, Savolainen O: Environmental and genetic effects on flowering differences between northern and southern populations of Arabidopsis lyrata (Brassicaceae). Am J Bot. 2004, 91: 1036-1045. 10.3732/ajb.91.7.1036.PubMedGoogle Scholar
- Riihimäki M, Podolsky R, Kuittinen H, Koelewijn H, Savolainen O: Studying genetics of adaptive variation in model organisms: flowering time variation in Arabidopsis lyrata. Genetica. 2005, 123: 63-74. 10.1007/s10709-003-2711-7.PubMedGoogle Scholar
- Clauss MJ, Dietel S, Schubert G, Mitchell-Olds T: Glucosinolate and trichome defenses in a natural Arabidopsis lyrata population. J Chem Ecol. 2006, 32: 2351-2373. 10.1007/s10886-006-9150-8.PubMedGoogle Scholar
- Kusaba M, Dwyer K, Hendershot J, Vrebalov J, Nasrallah JB, Nasrallah ME: Self-incompatibility in the genus Arabidopsis: characterization of the S locus in the outcrossing A. lyrata and its autogamous relative A. thaliana. Plant Cell. 2001, 13: 627-643. 10.1105/tpc.13.3.627.PubMed CentralPubMedGoogle Scholar
- Schierup MH, Mable BK, Awadalla P, Charlesworth D: Identification and characterization of a polymorphic receptor kinase gene linked to the self-incompatibility locus of Arabidopsis lyrata. Genetics. 2001, 158: 387-399.PubMed CentralPubMedGoogle Scholar
- Charlesworth D, Mable BK, Schierup MH, Bartolomé C, Awadalla P: Diversity and linkage of genes in the self-incompatibility gene family in Arabidopsis lyrata. Genetics. 2003, 164: 1519-1535.PubMed CentralPubMedGoogle Scholar
- Mable BK, Beland J, di Berardo C: Inheritance and dominance of self-incompatibility alleles in polyploid Arabidopsis lyrata. Heredity. 2004, 93: 476-486. 10.1038/sj.hdy.6800526.PubMedGoogle Scholar
- Mable BK, Robertson AV, Dart S, diBerardo C, Witham L: Breakdown of self-incompatibility in the perennial Arabidopsis lyrata (Brassicaceae) and its genetic consequences. Evolution. 2005, 59: 1437-1448.PubMedGoogle Scholar
- Charlesworth D, Kamau E, Hagenblad J, Tang C: Trans-specificity at loci near the self-incompatibility loci in Arabidopsis. Genetics. 2006, 172: 2699-2704. 10.1534/genetics.105.051938.PubMed CentralPubMedGoogle Scholar
- Hagenblad J, Bechsgaard J, Charlesworth D: Linkage disequilibrium between incompatibility locus region genes in the plant Arabidopsis lyrata. Genetics. 2006, 173: 1057-1073. 10.1534/genetics.106.055780.PubMed CentralPubMedGoogle Scholar
- Schierup MH, Bechsgaard JS, Nielsen LH, Christiansen FB: Selection at work in self-incompatible Arabidopsis lyrata: mating patterns in a natural population. Genetics. 2006, 172: 477-484. 10.1534/genetics.105.045682.PubMed CentralPubMedGoogle Scholar
- Schierup MH, Bechsgaard JS, Christiansen FB: Selection at work in self-incompatible Arabidopsis lyrata. II. Spatial distribution of S haplotypes in Iceland. Genetics. 2008, 180: 1051-1059. 10.1534/genetics.108.088799.PubMed CentralPubMedGoogle Scholar
- Sharbel TF, Haubold B, Mitchell-Olds T: Genetic isolation by distance in Arabidopsis thaliana: biogeography and postglacial colonization of Europe. Mol Ecol. 2000, 9: 2109-2118. 10.1046/j.1365-294X.2000.01122.x.PubMedGoogle Scholar
- Beck JB, Schmuths H, Schaal BA: Native range genetic variation in Arabidopsis thaliana is strongly geographically structured and reflects Pleistocene glacial dynamics. Mol Ecol. 2008, 17: 902-915. 10.1111/j.1365-294X.2008.03906.x.PubMedGoogle Scholar
- Wright SI, Lauga B, Charlesworth D: Subdivision and haplotype structure in natural populations of Arabidopsis lyrata. Mol Ecol. 2003, 12: 1247-1263. 10.1046/j.1365-294X.2003.01743.x.PubMedGoogle Scholar
- Balañá-Alcaide D, Ramos-Onsins SE, Boone Q, Aguadé M: Highly structured nucleotide variation within and among Arabidopsis lyrata populations at the FAH1 and DFR gene regions. Mol Ecol. 2006, 15: 2059-2068. 10.1111/j.1365-294X.2006.02918.x.PubMedGoogle Scholar
- Clauss MJ, Mitchell-Olds T: Population genetic structure of Arabidopsis lyrata in Europe. Mol Ecol. 2006, 15: 2753-2766.PubMedGoogle Scholar
- Muller MH, Leppälä J, Savolainen O: Genome-wide effects of postglacial colonization in Arabidopsis lyrata. Heredity. 2008, 100: 47-58. 10.1038/sj.hdy.6801057.PubMedGoogle Scholar
- Koch MA, Matschinger M: Evolution and genetic differentiation among relatives of Arabidopsis thaliana. Proc Natl Acad Sci USA. 2007, 104: 6272-6277. 10.1073/pnas.0701338104.PubMed CentralPubMedGoogle Scholar
- Schmickl R, Jørgensen MH, Brysting AK, Koch MA: Phylogeographic implications for the North American boreal-arctic Arabidopsis lyrata complex. Plant Ecol Div. 2008, 1: 245-254. 10.1080/17550870802349138.Google Scholar
- Avise JC, Arnold J, Ball RM, Bermingham E, Lamb T, Neigel JE, Reeb CA, Saunders NC: Intraspecific phylogeography: The mitochondrial DNA bridge between population genetics and systematics. Annu Rev Ecol Syst. 1987, 18: 489-522.Google Scholar
- Eidesen PB, Carlsen T, Molau U, Brochmann C: Repeatedly out of Beringia: Cassiope tetragona embraces the Arctic. J Biogeogr. 2007, 34: 1559-1574. 10.1111/j.1365-2699.2007.01719.x.Google Scholar
- Abott RJ, Comes HP: Evolution in the Arctic: a phylogeographic analysis of the circumarctic plant, Saxifraga oppositifolia (Purple saxifrage). New Phytol. 2004, 161: 211-224. 10.1046/j.1469-8137.2003.00953.x.Google Scholar
- Tremblay NO, Schoen DJ: Molecular phylogeography of Dryas integrifolia: glacial refugia and postglacial recolonization. Mol Ecol. 1999, 8: 1187-1198. 10.1046/j.1365-294x.1999.00680.x.Google Scholar
- Abbott RJ, Smith LC, Milne RI, Crawford RM, 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.PubMedGoogle Scholar
- Brochmann C, Gabrielsen TM, Nordal I, Landvik JY, Elven R: Glacial survival or tabula rasa? The history of North Atlantic biota revisited. Taxon. 2003, 52: 417-450. 10.2307/3647444.Google Scholar
- Alsos IG, Engelskjøn T, Gielly L, Taberlet P, Brochmann C: Impact of ice ages on circumpolar molecular diversity: insights from an ecological key species. Mol Ecol. 2005, 14: 2739-2753. 10.1111/j.1365-294X.2005.02621.x.PubMedGoogle Scholar
- Oliver C, Hollingsworth PM, Gornall RJ: Chloroplast DNA phylogeography of the arctic-montane species Saxifraga hirculus (Saxifragaceae). Heredity. 2006, 96: 222-231. 10.1038/sj.hdy.6800785.PubMedGoogle Scholar
- Al-Shehbaz IA, O'Kane SL: Taxonomy and phylogeny of Arabidopsis (Brassicaceae). The Arabidopsis Book. 2002, Rockville, MD: American Society of Plant BiologistsGoogle Scholar
- Elven R, (Ed): Checklist of the Panarctic Flora (PAF) Vascular Plants. Version. 2007, [http://www.binran.ru/infsys/paflist/index.htm]
- O'Kane SL, Al-Shehbaz IA: A synopsis of Arabidopsis (Brassicaceae). Novon. 1997, 7: 323-327. 10.2307/3391949.Google Scholar
- Al-Shehbaz IA, O'Kane SL, Price RA: Generic placement of species excluded from Arabidopsis (Brassicaceae). Novon. 1999, 9: 296-307. 10.2307/3391724.Google Scholar
- Koch MA, Wernisch M, Schmickl R: Arabidopsis thaliana 's wild relatives: an updated overview on systematics, taxonomy and evolution. Taxon. 2008, 57: 933-943.Google Scholar
- Warwick SI, Al-Shehbaz IA, Sauder CA: Phylogenetic position of Arabis arenicola and generic limits of Aphragmus and Eutrema (Brassicaceae) based on sequences of nuclear ribosomal DNA. Can J Bot. 2006, 84: 269-281. 10.1139/B05-161.Google Scholar
- Shimizu KK, Fujii S, Marhold K, Watanabe K, Kudoh H: Arabidopsis kamchatica (Fisch. ex DC.) K. Shimizu & Kudoh and A. kamchatica subsp. kawasakiana (Makino) K. Shimizu & Kudoh, new combinations. Acta Phytotax Geobot. 2005, 56: 163-172.Google 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
- Dobeš CH, Mitchell-Olds T, Koch MA: Extensive chloroplast haplotype variation indicates Pleistocene hybridization and radiation of North American Arabis drummondii, A. × divaricarpa, and A. holboellii (Brassicaceae). Mol Ecol. 2004, 13: 349-370. 10.1046/j.1365-294X.2003.02064.x.PubMedGoogle Scholar
- Dobeš C, Mitchell-Olds T, Koch MA: Intraspecific diversification in North American Boechera stricta (= Arabis drummondii), Boechera × divaricarpa, and Boechera holboellii (Brassicaceae) inferred from nuclear and chloroplast molecular markers - an integrative approach. Am J Bot. 2004, 91: 2087-2101. 10.3732/ajb.91.12.2087.PubMedGoogle Scholar
- Koch MA, Dobeš C, Matschinger M, Bleeker W, Vogel J, Kiefer M, Mitchell-Olds T: Evolution of the trn F(GAA) gene in Arabidopsis relatives and the Brassicaceae family: monophyletic origin and subsequent diversification of a plastidic pseudogene. Mol Biol Evol. 2005, 22: 1032-1043. 10.1093/molbev/msi092.PubMedGoogle Scholar
- Dobeš C, Kiefer C, Kiefer M, Koch MA: Plastidic trnF UUC pseudogenes in North American genus Boechera (Brassicaceae): mechanistic aspects of evolution. Plant Biol. 2007, 9: 502-515. 10.1055/s-2006-955978.PubMedGoogle Scholar
- Koch MA, Dobeš C, Kiefer C, Schmickl R, Klimeš L, Lysak MA: Supernetwork identifies multiple events of plastid trn F(GAA) pseudogene evolution in the Brassicaceae. Mol Biol Evol. 2007, 24: 63-73. 10.1093/molbev/msl130.PubMedGoogle Scholar
- Schmickl R, Kiefer C, Dobeš C, Koch MA: Evolution of trn F(GAA) pseudogenes in cruciferous plants. Plant Syst Evol. 2009, 282: 229-240. 10.1007/s00606-008-0030-2.Google Scholar
- Ansell SW, Schneider H, Pedersen N, Grundmann M, Russell SJ, Vogel JC: Recombination diversifies chloroplast trnF pseudogenes in Arabidopsis lyrata. J Evol Biol. 2007, 20: 2400-2411. 10.1111/j.1420-9101.2007.01397.x.PubMedGoogle Scholar
- Arnheim N, Krystal M, Schmickel R, Wilson G, Ryder O, Zimmer E: Molecular evidence for genetic exchanges among ribosomal genes on nonhomologous chromosomes in man and apes. Proc Natl Acad Sci USA. 1980, 77: 7323-7327. 10.1073/pnas.77.12.7323.PubMed CentralPubMedGoogle Scholar
- Liao D: Concerted evolution: molecular mechanism and biological implications. Am J Hum Genet. 1999, 64: 24-30. 10.1086/302221.PubMed CentralPubMedGoogle Scholar
- Koch MA, Dobeš C, Mitchell-Olds T: Multiple hybrid formation in natural populations: Concerted evolution of the internal transcribed spacer of nuclear ribosomal DNA (ITS) in North American Arabis divaricarpa (Brassicaceae). Mol Biol Evol. 2003, 20: 338-350. 10.1093/molbev/msg046.PubMedGoogle Scholar
- Carson AR, Scherer SW: Identifying concerted evolution and gene conversion in mammalian gene pairs lasting over 100 million years. BMC Evol Biol. 2009, 9: 156-10.1186/1471-2148-9-156.PubMed CentralPubMedGoogle Scholar
- Ehlers J, Gibbard PL, (Eds): Quaternary glaciations - extent and chronology. 2004, Amsterdam: ElsevierGoogle Scholar
- Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-1659. 10.1046/j.1365-294x.2000.01020.x.PubMedGoogle Scholar
- Templeton AR, Crandall KA, Sing CF: A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics. 1992, 132: 619-633.PubMed CentralPubMedGoogle Scholar
- Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol Bioinform Online. 2005, 1: 47-50.PubMed CentralGoogle Scholar
- Kumar S, Nei M, Dudley J, Tamura K: MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform. 2008, 9: 299-306. 10.1093/bib/bbn017.PubMed CentralPubMedGoogle Scholar
- Gottlieb LD: Evidence for duplication and divergence of the structural gene for phosphoglucoisomerase in diploid species of Clarkia. Genetics. 1977, 86: 289-308.PubMed CentralPubMedGoogle Scholar
- Ghatnekar L: A polymorphic duplicated locus for cytosolic PGI segregating in sheep's fescue (Festuca ovina L.). Heredity. 1999, 83: 451-459. 10.1038/sj.hdy.6885750.PubMedGoogle Scholar
- Ford VS, Gottlieb LD: Single mutations silence PgiC2 genes in two very recent allotetraploid species of Clarkia. Evolution. 2002, 56: 699-707.PubMedGoogle Scholar
- Abbott RT, Brochmann C: History and evolution of the arctic flora: in the footsteps of Eric Hulten. Mol Ecol. 2003, 12: 299-313. 10.1046/j.1365-294X.2003.01731.x.PubMedGoogle Scholar
- DeChain EG: A bridge or a barrier? Beringia's influence on the distribution and diversity of tundra plants. Plant Ecol Div. 2008, 1: 197-207. 10.1080/17550870802328660.Google Scholar
- Tribsch A, Schönswetter P: Patterns of endemism and comparative phylogeography confirm palaeoenvironmental evidence for Pleistocene refugia in the Eastern Alps. Taxon. 2003, 52: 477-497. 10.2307/3647447.Google Scholar
- Hewitt GM: The structure of biodiversity - insights from molecular phylogeography. Front Zool. 2004, 1: 4-10.1186/1742-9994-1-4.PubMed CentralPubMedGoogle Scholar
- Koch MA, Kiefer C, Ehrich D, Vogel J, Brochmann C, Mummenhoff K: Three times out of Asia Minor: the phylogeography of Arabis alpina L. (Brassicaceae). Mol Ecol. 2006, 15: 825-839. 10.1111/j.1365-294X.2005.02848.x.PubMedGoogle Scholar
- Ehrich D, Gaudeul M, Assefa A, Koch MA, Mummenhoff K, Nemomissa S, Brochmann C: Genetic consequences of Pleistocene range shifts: contrast between the Arctic, the Alps and the East African mountains. Mol Ecol. 2007, 16: 2542-2559. 10.1111/j.1365-294X.2007.03299.x.PubMedGoogle Scholar
- Schönswetter P, Paun O, Tribsch A, Niklfeld H: Out of the Alps: colonization of northern Europe by east alpine populations of the glacier buttercup Ranunculus glacialis L. (Ranunculaceae). Mol Ecol. 2003, 12: 3373-3381. 10.1046/j.1365-294X.2003.01984.x.PubMedGoogle Scholar
- Ehrich D, Alsos IG, Brochmann C: Where did the northern peatland species survive the dry glacials: cloudberry (Rubus chamaemorus) as an example. J Biogeogr. 2008, 35: 801-814. 10.1111/j.1365-2699.2007.01864.x.Google Scholar
- Savile DBO: Arctic adaptations in plants. Canada Department of Agriculture Monograph. 1972, 6: 1-81.Google Scholar
- Alsos IG, Eidesen PB, Ehrich D, Skrede I, Westergaard K, Jacobsen GH, Landvik JY, Taberlet P, Brochmann C: Frequent long-distance plant colonization in the changing Arctic. Science. 2007, 316: 1606-1609. 10.1126/science.1139178.PubMedGoogle Scholar
- Brochmann C, Brysting AK, Alsos IG, Borgen L, Grundt HH, Scheen AC, Elven R: Polyploidy in arctic plants. Biol J Linn Soc. 2004, 82: 521-536. 10.1111/j.1095-8312.2004.00337.x.Google Scholar
- Manton I: The problem of Biscutella laevigata L. - II. The evidence from meiosis. Ann Bot. 1937, 1: 439-462.Google Scholar
- Stebbins GL: Polyploidy and the distribution of the arctic-alpine flora: new evidence and a new approach. Bot Helv. 1984, 94: 1-13.Google Scholar
- Soltis PS, Soltis DE: The role of genetic and genomic attributes in the success of polyploids. Proc Natl Acad Sci USA. 2000, 97: 7051-7057. 10.1073/pnas.97.13.7051.PubMed CentralPubMedGoogle Scholar
- Comai L, Tyagi AP, Winter K, Holmes-Davis R, Reynolds SH, Stevens Y, Byers B: Phenotypic instability and rapid gene silencing in newly formed Arabidopsis allotetraploids. Plant Cell. 2000, 12: 1551-1568. 10.1105/tpc.12.9.1551.PubMed CentralPubMedGoogle Scholar
- Wendel JF: Genome evolution in polyploids. Plant Mol Biol. 2000, 42: 225-249. 10.1023/A:1006392424384.PubMedGoogle Scholar
- Liu B, Wendel JF: Epigenetic phenomena and the evolution of plant allopolyploids. Mol Phylogenet Evol. 2003, 29: 365-379. 10.1016/S1055-7903(03)00213-6.PubMedGoogle 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.PubMedGoogle Scholar
- Guggisberg A, Mansion G, Kelso S, Conti E: Evolution of biogeographic patterns, ploidy levels, and breeding systems in a diploid-polyploid species complex of Primula. New Phytol. 2006, 171: 617-632.PubMedGoogle Scholar
- Guggisberg A, Baroux C, Grossniklaus U, Conti E: Genomic origin and organization of the allopolyploid Primula egaliksensis investigated by in situ hybridization. Ann Bot. 2008, 101: 919-927. 10.1093/aob/mcn026.PubMed CentralPubMedGoogle Scholar
- Guggisberg A, Mansion G, Conti E: Disentangling reticulate evolution in an arctic-alpine polyploid complex. Syst Biol. 2009, 58: 55-73. 10.1093/sysbio/syp010.PubMedGoogle Scholar
- Popp M, Erixon P, Eggens F, Oxelman B: Origin and evolution of a circumpolar polyploid species complex in Silene (Caryophyllaceae) inferred from low copy nuclear RNA polymerase introns, rDNA, and chloroplast DNA. Syst Bot. 2005, 30: 302-313. 10.1600/0363644054223648.Google Scholar
- Popp M, Oxelman B: Origin and evolution of North American polyploid Silene (Caryophyllaceae). Am J Bot. 2007, 94: 330-349. 10.3732/ajb.94.3.330.PubMedGoogle Scholar
- Brochmann C, Xiang QY, Brunsfeld SJ, Soltis DE, Soltis PS: Molecular evidence for polyploid origins in Saxifraga (Saxifragaceae): the narrow arctic endemic S. svalbardensis and its widespread allies. Am J Bot. 1998, 85: 135-143. 10.2307/2446562.PubMedGoogle Scholar
- Jørgensen MH, Elven R, Tribsch A, Gabrielsen TM, Stedje B, Brochmann C: Taxonomy and evolutionary relationships in the Saxifraga rivularis complex. Syst Bot. 2006, 31: 702-729. 10.1600/036364406779695988.Google Scholar
- Jørgensen MH, Carlsen T, Skrede I, Elven R: Microsatellites resolve the taxonomy of the polyploid Cardamine digitata aggregate (Brassicaceae). Taxon. 2008, 57: 882-892.Google Scholar
- Kapralov MV, Gabrielsen TM, Sarapultsev IE, Brochmann C: Genetic enrichment of the arctic clonal plant Saxifraga cernua at its southern periphery via the alpine sexual Saxifraga sibirica. Mol Ecol. 2006, 15: 3401-3411. 10.1111/j.1365-294X.2006.03024.x.PubMedGoogle Scholar
- Säll T, Jakobsson M, Lind-Halldén C, Halldén C: Chloroplast DNA indicates a single origin of the allotetraploid Arabidopsis suecica. J Evol Biol. 2003, 16: 1019-1029. 10.1046/j.1420-9101.2003.00554.x.PubMedGoogle Scholar
- Jakobsson M, Hagenblad J, Tavaré S, Säll T, Halldén C, Lind-Halldén C, Nordborg M: A unique recent origin of the allotetraploid species Arabidopsis suecica: Evidence from nuclear DNA markers. Mol Biol Evol. 2006, 23: 1217-1231. 10.1093/molbev/msk006.PubMedGoogle Scholar
- Säll T, Lind-Halldén C, Jakobsson M, Halldén C: Mode of reproduction in Arabidopsis suecica. Hereditas. 2004, 141: 313-317. 10.1111/j.1601-5223.2004.01833.x.PubMedGoogle Scholar
- Sugisaka J, Kudoh H: Breeding system of the annual Cruciferae, Arabidopsis kamchatica subsp. kawasakiana. J Plant Res. 2008, 121: 65-68. 10.1007/s10265-007-0119-7.PubMedGoogle Scholar
- Bechsgaard JS, Castric V, Charlesworth D, Vekemans X, Schierup MH: The transition to self-compatibility in Arabidopsis thaliana and evolution within S-haplotypes over 10 Myr. Mol Biol Evol. 2006, 23: 1741-1750. 10.1093/molbev/msl042.PubMedGoogle Scholar
- Tang C, Toomajian C, Sherman-Broyles S, Plagnol V, Guo YL, Hu TT, Clark RM, Nasrallah JB, Weigel D, Nordborg M: The evolution of selfing in Arabidopsis thaliana. Science. 2007, 317: 1070-1072. 10.1126/science.1143153.PubMedGoogle Scholar
- Schoen DJ, Johnston MO, L'Heureux AM, Marsolais JV: Evolutionary history of the mating system in Amsinckia (Boraginaceae). Evolution. 1997, 51: 1090-1099. 10.2307/2411038.Google Scholar
- Rieseberg LH, Willis JH: Plant speciation. Science. 2007, 317: 910-914. 10.1126/science.1137729.PubMed CentralPubMedGoogle Scholar
- Nasrallah JB, Liu P, Sherman-Broyles S, Schmidt R, Nasrallah ME: Epigenetic mechanisms for breakdown of self-incompatibility in interspecific hybrids. Genetics. 2007, 175: 1965-1973. 10.1534/genetics.106.069393.PubMed CentralPubMedGoogle Scholar
- Mable BK: Genetic causes and consequences of the breakdown of self-incompatibility: case studies in the Brassicaceae. Genet Res. 2008, 90: 47-60.Google Scholar
- Eidesen PB, Alsos IG, Popp M, Stensrud Ø, Suda J, Brochmann C: Nuclear vs. plastid data: complex Pleistocene history of a circumpolar key species. Mol Ecol. 2007, 16: 3902-3925. 10.1111/j.1365-294X.2007.03425.x.PubMedGoogle Scholar
- Goloboff PA, Farris JS, Nixon KC: TNT (Tree analysis using New Technology). 2000, Published by the authors, TucumánGoogle 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.