text:
- Research article
- Open Access
Taming the wild: resolving the gene pools of non-model Arabidopsislineages
© Hohmann et al.; licensee BioMed Central Ltd. 2014
- Received: 27 June 2014
- Accepted: 15 October 2014
- Published: 27 October 2014
Wild relatives in the genus Arabidopsis are recognized as useful model systems to study traits and evolutionary processes in outcrossing species, which are often difficult or even impossible to investigate in the selfing and annual Arabidopsis thaliana. However, Arabidopsis as a genus is littered with sub-species and ecotypes which make realizing the potential of these non-model Arabidopsis lineages problematic. There are relatively few evolutionary studies which comprehensively characterize the gene pools across all of the Arabidopsis supra-groups and hypothesized evolutionary lineages and none include sampling at a world-wide scale. Here we explore the gene pools of these various taxa using various molecular markers and cytological analyses.
Based on ITS, microsatellite, chloroplast and nuclear DNA content data we demonstrate the presence of three major evolutionary groups broadly characterized as A. lyrata group, A. halleri group and A. arenosa group. All are composed of further species and sub-species forming larger aggregates. Depending on the resolution of the marker, a few closely related taxa such as A. pedemontana, A. cebennensis and A. croatica are also clearly distinct evolutionary lineages. ITS sequences and a population-based screen based on microsatellites were highly concordant. The major gene pools identified by ITS sequences were also significantly differentiated by their homoploid nuclear DNA content estimated by flow cytometry. The chloroplast genome provided less resolution than the nuclear data, and it remains unclear whether the extensive haplotype sharing apparent between taxa results from gene flow or incomplete lineage sorting in this relatively young group of species with Pleistocene origins.
Our study provides a comprehensive overview of the genetic variation within and among the various taxa of the genus Arabidopsis. The resolved gene pools and evolutionary lineages will set the framework for future comparative studies on genetic diversity. Extensive population-based phylogeographic studies will also be required, however, in particular for A. arenosa and their affiliated taxa and cytotypes.
- Chloroplast
- Cytology
- Evolution
- Microsatellites
- Systematics
- Taxonomy
Arabidopsis: life in the fast lane
Less than a decade ago "Arabidopsis and its poorly known relatives" was the title chosen to introduce the closest relatives of Arabidopsis thaliana to a broader readership [1]. This review summarized both the systematics and taxonomy of the genus and also the ecologically important traits to be studied in A. thaliana's "wild" relatives. Its necessity was obvious because until 1999, a huge number of species (60) were recognized in Arabidopsis in the traditional sense. Arabidopsis' taxonomical history was compiled in detail more than 10 years ago [2],[3], and nine Arabidopsis species with several subspecies were recognized by this time. Based on this work and unraveling the evolutionary history of the genus Arabis[4]-[7], which differs morphologically from Arabidopsis only in the position of the cotyledons relative to the radicle in the seeds, a new systematic concept was presented 10-15 years ago [4],[8],[9]. Several species and subspecies have since been added either because molecular studies provided new resolution [10] or because description of new species [11] led to changes in their respective taxonomic rank (species, subspecies, variety) [11]-[17].
Arabidopsis species diversity and taxonomy
Arabidopsis arenosaspecies aggregate
Arabidopsis arenosa (L.) Lawalrée
subsp. arenosa
(2n = 32)
Central and Western Europe, Scandinavia (lower altitudes)
subsp. arenosa var. intermedia (Kovats) Hayek
(2n = 32)
Southeastern Austrian Alps (similar to A. neglecta)
subsp. borbasii (Zapałowicz) O'Kane & Al-Shehbaz
(2n = 32)
Central and Western Europe (mountain ranges, higher altitudes)
Arabidopsis arenosa, unclear taxonomic treatment
(2n = 16)
Arabidopsis carpatica, nom. prov.
(2n = 16)
Carpathians (middle altitudes, calcareous bedrocks)
Arabidopsis neglecta (Schultes) O'Kane & Al-Shehbaz
subsp. neglecta
(2n = 16)
Carpathians (alpine ranges)
subsp. robusta, nom. prov.
(2n = 32)
Carpathians (alpine ranges, only occasionally in lower altitudes)
Arabidopsis nitida, nom. prov.
(2n = 16)
Carpathians (mountain ranges, middle to subalpine altitudes)
Arabidopsis petrogena (A. Kern) V.I. Dorof.
?subsp. petrogena
(2n = 16)
Carpathians, Pannonian lowland (maybe two varieties)
subsp. exoleta, nom. prov.
(2n = 32)
Carpathians (lower altitudes)
Arabidopsis lyrata lineage
Arabidopsis lyrata subsp. lyrata (L.) O'Kane & Al-Shehbaz
(2n = 16)
Alaska, Canada, United States
Arabidopsis lyrata subsp. petraea (L.) O'Kane & Al-Shehbaz
(2n = 16/32)
= A. petraea (L.) V.I. Dorof.
Arabidopsis petraea subsp. umbrosa (Turcz. Ex Steud.) Elven & D.F. Murray
(2n = 16)
Arctic NE Asia, Siberia, Alaska, Canada
Arabidopsis petraea subsp. septentrionalis (N. Busch) Elven & D.F. Murray
(2n = 32)
Arctic NE Europe, European Russia to Siberia
Arabidopsis arenicola (Richardson ex Hook.) Al-Shehbaz et al.
(2n = 16)
Arctic Canada and Greenland
Arabidopsis halleri lineage
Arabidopsis halleri subsp. halleri (L.) O'Kane & Al-Shehbaz
(2n = 16)
Arabidopsis halleri subsp. dacica (Heuff.) Kolník
(2n = 16)
Carpathians, Romania
Arabidopsis halleri subsp. gemmifera (Matsum.) O'Kane & Al-Shehbaz
(2n = 16)
Russia Far East, NE China, Korea, Japan, and Taiwan
Arabidopsis halleri subsp. ovirensis (Wulfen) A. P. Iljinsk.
(2n = 16)
Austria only (all accessions from the Balkans belong to subsp. halleri)
Arabidopsis halleri subsp. tatrica (Pawł.) Kolník
(2n = 16)
Tatra mountains, Slovakia
Arabidopsis umezawana Kadota
(2n = ?)
Japan, Hokkaido (alpine zone of Mt. Rishirizin), annual to biennial
Other diploid taxa
Arabidopsis pedemontana (Boiss.) O'Kane & Al-Shehbaz
(2n = 16)
NW Italy
Arabidopsis cebennensis (DC.) O'Kane & Al-Shehbaz
(2n = 16)
SE France, Massif Central
Arabidopsis croatica (Schott) O'Kane & Al-Shehbaz
(2n = 16)
Allopolyploid taxa
Arabidopsis kamchatica (Fisch. Ex DC.) O'Kane & Al-Shehbaz
(2n = 32)
Boreal Alasca, Canada, E Siberia, Russian Far East, Korea, Japan,
Arabidopsis kamchatica subsp. kawasakiana (Makino) Shimizu & Kudoh
(2n = 32)
Japan, winterannual (coastal, lowland)
Arabidopsis suecica (Fr.) Norrl.
(2n = 26)
Fennoscandinavia and the Baltic region
Monophyly is generally accepted among Arabidopsis taxa by plant scientists at present. However, considering that A. thaliana is a model system taxonomic recognition of new species as Arabidopsis is acknowledged much faster than comparable systematic-taxonomic changes in other genera. One such contrary example from the Brassicaceae family is the genus Noccaea which includes important model species for heavy metal tolerance and hyperaccumulation. Noccaea caerulescens required more than 30 years to be recognized appropriately within the correct evolutionary framework [20],[21]. Systematics and taxonomy in the genus Arabidopsis is thus ever-debatable and in constant need of further improvement.
Developing a comprehensive systematic framework
To date there is limited genetic information across the entire genus which allows for adequate taxonomic and systematic comparison. The first study highlighting centers of genetic variation in Europe for the main evolutionary lineages also provided evidence for extensive shared plastidic variation among species [22]. The female component of nuclear-encoded self-incompatibility genes (SI alleles at the SRK locus) also revealed trans-specific polymorphism among some of the same species [23].
Some major evolutionary lineages have been identified in the Arabidopsis genus [18],[22], namely the following groups: A. halleri, A. lyrata and A. arenosa. Three other genetically isolated diploid species have been identified, A. croatica, A. cebennensis and A. pedemontana. A few allopolyploids are also well studied: A. suecica with A. arenosa and A. thaliana as parental species [24],[25], and A. kamchatica with A. lyrata and A. halleri (subsp. gemmifera) as parents [26]-[28]. Another taxonomically not yet introduced tetraploid taxon (close to A. lyrata) is found in Lower Austria, which is either the result of hybridization and genome doubling between A. arenosa and A. lyrata (allopolyploidy), or genome duplication of diploid A. lyrata (autopolyploid) with subsequent introgression from tetraploid A. arenosa[29].
For some of these major lineages and their subspecies there are more detailed genetic studies available covering either a broader geographic scale or larger sets of taxa. For A. halleri it has been shown that all five subspecies are closely related to each other, and that one major center of genetic diversity is located in the Eastern Austrian Alps [22]. It has also been concluded for A. halleri that metallicolous populations have been founded separately from distinct non-metallicolous populations without suffering from founder effects [30]. The same authors provided a comprehensive phylogeographic scenario [31]; and although the accessions studied were not characterized taxonomically, many helpful comments linking taxonomy with genetic data were provided. For A. lyrata there are several studies available showing general phylogeographic patterns and hybrid speciation on a large scale [26],[27]. Local-scale phylogeographic studies in North America highlighted switches in mating system [32],[33]. Population-based analysis with a few selected populations provided the first evidence for population genetic structure at varying geographic scales [34]-[36]. At a more local scale and focusing on different aspects of adaptation there are numerous contributions covering A. lyrata[37]-[39], and comprehensive reviews have recently been presented to summarize many more aspects [1],[40]. There is very limited information regarding A. arenosa, one of the most diverse evolutionary lineages in Arabidopsis[22], with only one phylogeographic-systematic study at a broad geographic scale [19]. Nevertheless, A. arenosa has proven to be an excellent model to study the formation and evolution of allopolyploids [41],[42] and plant adaptation [43],[44].
A recently published review [45] emphasized the need for all-encompassing evolutionary studies within the genus Arabidopsis that provide a broader framework on genus-wide genetic diversity and differentiation, in order to enable researchers to study molecular mechanisms of speciation-related processes in interspecific comparative approaches. Our goal here was to provide a reliable phylogenetic-systematic base line using ribosomal DNA sequence variation from the internal transcribed spacers 1 and 2 and the trnLF region of the plastid genome [22]. These data were combined with population genetic variation based on a set of nuclear-encoded microsatellite loci shown to be highly sensitive for resolving Arabidopsis lineages [29]. Finally, since genome size and chromosome numbers are important cytological characters that significantly influence various organismal traits, we conducted a comprehensive scan of cytological variation via the homoploid nuclear DNA content within and among the principal gene pools in Arabidopsis.
Here we explore the gene pools of Arabidopsis taxa using a battery of molecular markers and their cytology to identify clearly genetically distinct units over their entire geographic distribution, develop a schematic phylogeographic-systematic scenario based on this data and lastly, comment on any discrepancies between these resolved gene pools and existing taxonomic identifiers.
The number of single nucleotide polymorphisms (SNPs) was not sufficient to resolve taxa below the species level, most likely because the genus' radiation within the last 2.5 million years is too recent.
ITS sequence data recognize major gene pools
Estimation of absolute genome size data of Arabidopsis taxa
1Cx-value (pg) (SD)
Expected genome size
Calculated genome size1
Arabidopsis thaliana (Columbia) (2×)
0.173 (0.010)
135 Mbp(1)
(1)set as reference
Arabidopsis croatica (2×)
0.225 (0.007)
176 Mbp
Arabidopsis arenicola (2×)
0.260 (0.008)
203 Mbp
Arabidopsis lyrata (2×) (Europe)
0.270 (0.005)
210 Mbp
Arabidopsis lyrata (2×) (North America)
0.274 (0.007)
207 Mbp(2)
214 Mbp
Arabidopsis lyrata (4×)
0.238 (0.005)
371 Mbp
Arabidopsis halleri ssp. halleri (2×)
0.266 (0.010)
207 Mbp
Arabidopsis halleri ssp. gemmifera (2×)
0.270 (0.004)
211 Mbp
Arabidopsis kamchatica (4×)
0.235 (0.007)
367 Mbp
Arabidopsis arenosa (4×)
0.237 (0.005)
370 Mbp
Arabidopsis cebennensis (2×)
0.281 (0.010)
219 Mbp
Arabidopsis pedemontana (2×)
0.277 (0.005)
216 Mbp
Nuclear DNA content supports the distinction of major ITS gene pools
In contrast to among-group nuclear DNA content, variation was markedly reduced within the major gene pools and differences among accessions were minimal (the DNA content varied 1.18 fold, 1.1 fold and 1.14 fold, within diploid A. lyrata, A. arenosa, and A. halleri gene pools, respectively; and 1.01 fold and 1.14 fold within tetraploid A. lyrata and A. arenosa, respectively). In tetraploids, the variation was 1.01 fold and 1.14 fold among A. lyrata and A. arenosa accessions, respectively.
Absolute genome size estimates are also provided for all taxa in Table 2 in a Brassicaceae-wide screen. The 1C-value of A. thaliana ecotype Columbia is about 0.17 pg. The estimated physical size of its genome is currently about 135 Mbp, (TAIR, http://www.arabidopsis.org/). Our estimated 1C-value of 0.274 pg for North American A. lyrata indicates a respective genome size of approximately 214 Mbp and is very close to the published physical size (207 Mbp) of the A. lyrata genome [46]. The discrepancy of about 5% could be explained by missing sequence data from centromeric regions.
Chloroplast sequence data recognize some major gene pools but indicates shared polymorphism
Microsatellite analyses characterize distinctive taxa and cytotypes
Arabidopsis microsatellite genetic variation
Mating system
Individuals (populations)
No. of alleles
Mean no. of alleles per locus (SD)2
No. of unique alleles
No. of rare alleles (5%)
A. arenicola
19 (14)
1.04 (0.10)
A. lyrata subsp. lyrata
57 (29)
1.18 (0.16)
A. petraea susbp. umbrosa
25 (16)
1.20 (0.26)
A. petraea subsp. septentrionalis
8 (8)
1.63 (0.31)
A. lyrata subsp. petraea
187 (13)
1.46 (0.19)
A. lyrata subsp. petraea
25 (5)
1.93 (0.37)
A. arenosa subsp. arenosa
76 (13)
2.31 (0.33)
A. arenosa subsp. arenosa var. intermedia
14 (4)
2.19 (0.27)
A. arenosa subsp. borbasii
160 (22)
2.30 (0.33)
A. carpatica
113 (9)
1.59 (0.19)
A. petrogena subsp. petrogena
73 (7)
1.53 (0.20)
A. petrogena subsp. exoleta
56 (11)
2.22 (0.36)
A. neglecta subsp. neglecta
35 (5)
1.51 (0.21)
A. neglecta subsp. robusta
8 (1)
2.24 (0.23)
A. nitida
22 (4)
1.54 (0.22)
A. croatica
15 (3)
1.33 (0.19)
A. halleri subsp. dacica
3 (3)
1.53 (0.31)
A. halleri subsp. halleri
199 (19)
1.37 (0.19)
A. halleri subsp. ovirensis
24 (4)
1.31 (0.16)
A. halleri subsp. tatrica
25 (7)
1.37 (0.14)
A. halleri subsp. gemmifera
8 (5)
1.10 (0.15)
A. cebennensis
153 (11)
1.17 (0.14)
A. pedemontana
40 (9)
1.31 (0.18)
Gene diversity statistics of microsatellites, ITS and cpDNA variation
ITS haplotypes
trnLF haplotypes
Individuals analyzed
Nei's gene diversity (SD)
Individuals analyzed
Nucleotide diversity (πx10-2) (SD)
Nei's gene diversity (SD)
Individuals analyzed
Nucleotide diversity (πx10-2) (SD)
Nei's gene diversity (SD)
A. lyrata group
288 (85)
0.562 (0.312)
126 (92)
0.248 (0.168)
0.797 (0.025)
176 (101)
0.224 (0.150)
0.751 (0.024)
A. arenicola
19 (14)
0.205 (0.156)
17 (16)
0.139 (0.117)
0.323 (0.135)
17 (15)
0.035 (0.047)
0.117 (0.101)
A. lyrata subsp. lyrata
57 (29)
0.378 (0.225)
37 (21)
0.014 (0.030)
0.516 (0.085)
54 (26)
0.058 (0.061)
0.380 (0.065)
A. petraea subsp. umbrosa
25 (16)
0.468 (0.272)
30 (18)
0.154 (0.118)
0.675 (0.061)
28 (18)
0.362 (0.225)
0.738 (0.053)
A. petraea subsp. septentrionalis
8 (8)
not calculated
10 (10)
0.030 (0.046)
0.200 (0.154)
9 (9)
0.323 (0.223)
0.694 (0.147)
A. lyrata subsp. petraea
187 (13)
0.556 (0.309)
31 (26)
0.025 (0.039)
0.898 (0.030)
61 (29)
0.229 (0.154)
0.766 (0.042)
A. lyrata subsp. petraea
25 (5)
not calculated
1 (1)
1.000 (0.000)
7 (4)
0.241 (0.185)
0.714 (0.180)
A. arenosa group
258 (79)
0.560 (0.311)
247 (181)
0.138 (0.107)
0.803 (0.024)
568 (263)
0.171 (0.123)
0.585 (0.023)
A. arenosa subsp. arenosa
76 (13)
not calculated
23 (23)
0.129 (0.106)
0.806 (0.061)
32 (28)
0.163 (0.122)
0.485 (0.107)
A. arenosa var. intermedia
14 (4)
not calculated
4 (3)
1.000 (0.176)
6 (5)
0.149 (0.133)
0.333 (0.215)
A. arenosa subsp. borbasii
160 (22)
not calculated
173 (120)
0.065 (0.066)
0.748 (0.033)
391 (167)
0.142 (0.107)
0.503 (0.029)
A. carpatica
113 (9)
0.554 (0.309)
8 (6)
0.065 (0.075)
0.750 (0.139)
51 (14)
0.282 (0.181)
0.752 (0.046)
A. petrogena
73 (7)
0.512 (0.289)
5 (4)
0.399 (0.298)
1.000 (0.298)
37 (20)
0.052 (0.058)
0.348 (0.077)
A. petrogena subsp. exoleta
56 (11)
not calculated
8 (6)
0.345 (0.241)
0.928 (0.084)
8 (6)
0.341 (0.238)
0.892 (0.085)
A. neglecta subsp. neglecta
35 (5)
0.451 (0.262)
8 (5)
0.323 (0.229)
0.928 (0.084)
16 (5)
0.111 (0.097)
0.450 (0.150)
A. neglecta subsp. robusta
8 (1)
not calculated
6 (3)
0.600 (0.215)
6 (3)
0.099 (0.101)
0.333 (0.215)
A. nitida
22 (4)
0.635 (0.354)
4 (4)
0.384 (0.309)
1.000 (0.176)
11 (6)
0.282 (0.196)
0.709 (0.099)
A. croatica
15 (3)
0.374 (0.228)
8 (7)
0.384 (0.263)
1.000 (0.062)
10 (9)
0.159 (0.129)
0.533 (0.094)
A. halleri group
259 (38)
0.427 (0.254)
103 (90)
0.159 (0.118)
0.901 (0.020)
94 (83)
0.268 (0.173)
0.712 (0.030)
A. halleri subsp. dacica
3 (3)
0.533 (0.380)
8 (7)
0.049 (0.063)
0.928 (0.084)
8 (7)
0.213 (0.165)
0.464 (0.200)
A. halleri subsp. halleri
199 (19)
0.392 (0.237)
67 (61)
0.036 (0.047)
0.858 (0.029)
62 (58)
0.236 (0.158)
0.670 (0.049)
A. halleri subsp. ovirensis
24 (4)
0.330 (0.204)
5 (4)
0.122 (0.122)
0.900 (0.161)
2 (1)
A. halleri subsp. tatrica
25 (7)
0.408 (0.242)
13 (10)
0.023 (0.039)
0.730 (0.096)
12 (10)
0.364 (0.238)
0.727 (0.113)
A. halleri subsp. gemmifera
8 (5)
0.301 (0.205)
9 (7)
0.570 (0.363)
0.916 (0.092)
9 (6)
A. umezawana
not calculated
1 (1)
1.000 (0.000)
1 (1)
1.000 (0.000)
A. cebennensis
153 (11)
0.189 (0.130)
8 (6)
0.185 (0.150)
0.785 (0.150)
148 (12)
0.174 (0.125)
0.702 (0.018)
A. pedemontana
40 (9)
0.259 (0.167)
2 (2)
0.312 (0.382)
1.000 (0.500)
9 (2)
A. arenosa, on the other hand, fell into K = 4 populations; all tetraploid A. arenosa including tetraploid A. neglecta and A. petrogena are set apart from the diploid taxa. Vice versa the diploid A. carpatica and diploid A. petrogena formed distinct groups, respectively. The remaining diploids A. neglecta, A. croatica and A. nitida are combined to form one group. Interestingly, Structure was indifferent to the level of ploidy and the results are total in agreement with the ITS data. However, the initial round of analysis (combining diploids and tetraploids, Figure 5B) produced minor incongruencies such as the occurrence of distinct A. pedemontana genetic variation in other species such as A. umbrosa. This similarity is immediately eliminated when increasing K to the next higher values (data not shown).
When the Structure analysis is strictly confined to each ploidy level separately, the results are more congruent with better resolution (Figure 5A). In the diploid A. lyrata dataset (K = 2/3), A. arenicola was again not separated from A. lyrata subsp. lyrata, but A. umbrosa and A. lyrata subsp. petraea were still distinguished from each other (Figure 5A, upper part). Note that we report multiple K's here and their corresponding barplots in Figure 5 to reflect the fact that delta K was frequently either very similar between two independent runs, or because the less optimal K was more biologically meaningful. Both the optimal K, and the next-best optimal K in the Structure analysis are thus reported.
The diploid A. arenosa group was structured with K = 3/4, and A. carpatica, A. petrogena subsp. petrogena, A. neglecta subsp. neglecta and A. croatica were significantly recognized. Arabidopsis nitida (only few samples analyzed) was less clearly recognized. The five subspecies of A. halleri (K = 4) were grouped with some genetic clusters that distinguished a) A. halleri subsp. ovirensis, b) A. halleri subsp. tatrica/gemmifera, and c) a more complex and mixed cluster of A. halleri subsp. halleri (Figure 5A, upper part). With subsp. dacica the results should be interpreted with caution since only three individuals were analyzed. In summary the structure within A. halleri subsp. halleri is not clear and possibly indicates the need for taxonomic re-evaluation after comprehensive phylogeographic analysis [30],[31].
Analysis of the tetraploid dataset resulted in two unambiguously detectable genetic clusters only (Figure 5C, upper part) and distinguished A. lyrata from A. arenosa. However, it is interesting to see that A. septentrionalis carried approximately 50% genetic admixture from the A. arenosa genetic cluster. When analyzing tetraploid A. lyrata separately (K = 2), again A. septentrionalis was significantly different from tetraploid A. lyrata subsp. petraea from Scotland and Austria. In contrast, the tetraploid A. arenosa group was structured much less clearly (K = 2). At best, two groups could be identified: a) mountainous A. arenosa subsp. borbasii and high-alpine A. arenosa var. intermedia, and b) remaining A. arenosa subsp. arenosa, A. petrogena subsp. exoleta and A. neglecta subsp. robusta. It should be noted that neither ITS data nor plastid DNA data differentiated these groups further.
Genetic diversity is similar in all major groups
Gene diversity and nucleotide diversity is similar among the various groups of taxa (Table 4). Microsatellite gene diversity is highest in A. lyrata and A. arenosa and significantly lower in A. halleri. But this pattern is reverted when considering plastid DNA, where A. arenosa shows significantly lower diversity values compared to the A. lyrata and A. halleri groups. ITS diversity values could be summarized in that A. lyrata comprises the most diverse group. Thus, there is some coincidence with wide distribution ranges (A. lyrata and A. halleri) and high overall genetic variation. However, considering that the A. arenosa group has a much smaller total distribution range compared to the others, it is remarkable that levels of genetic diversity are also high. For the two local endemics, A. cebennensis and A. pedemontana, genetic diversity values are consistently lower.
Mating system affects genetic diversity in the various species and populations. However, detailed information on sporophytic self-incompatibility and mating system is available for A. lyrata and A. halleri only [50]-[52]. Both are self-incompatible with few exceptions (e.g. few populations of A. lyrata subsp. lyrata). Also for the A. arenosa lineage there are only reports of a fully self-incompatibility system [29] and there is only one questionable report of a selfing population, so far [53]. Two of the proven allopolyploids (A. suecica, A. kamchatica, not analyzed herein) are self-compatible [24],[54]. For many of the remaining taxa and cytotypes no data were available, and we added our results from many inbreeding experiments at Heidelberg Botanical Garden (2003-2014) to Table 3. Most of these taxa are also self-incompatible, and only for A. arenicola was self-compatibility shown, which is well-reflected in lowest number of alleles per locus (Table 3) and gene and nucleotide diversity of any marker system used herein (Table 4). Gene diversity (microsatellites, ITS and cpDNA) and nucleotide diversity (ITS and cpDNA) are for both A. arenicola and A. lyrata subsp. lyrata significantly lower than the respective mean values of the whole lyrata group (t-test: P <0.01). The self-incompatible mating system demonstrated for A. cebennensis in our cultivation experiments might not fit with its low values of number of alleles per locus (Table 3) or gene diversity (Table 4). However, these low numbers might be also explained simply by its narrow endemic distribution and small population sizes.
Polyploidy characterizes species groups differently
Mapping of polyploidy levels across the different taxa in our sample reveals that some lineages consist of diploids only (A. halleri, A. cebennensis, A. pedemontana, Table 1). The origins of tetraploid lineages are less clear e.g. tetraploid A. lyrata occurs at low frequency in Great Britain and Austria and there is evidence of introgression from tetraploid close relatives such as A. arenosa[29]. Distinguishing between a simple doubling of diploid A. lyrata genomes within a single ancestral population (autopolyploidy), or the establishment of polyploid lineages as a result of hybridization and genome doubling between two divergent species (A. lyrata/A. arenosa, allopolyploidy) requires further investigation in this system.
For others such as tetraploid A. septentrionalis, no evidence has been presented for a hybrid origin. The most diverse group of taxa with respect to ploidy variation is the tetraploid and diploid lineages within A. arenosa. Of the ten listed taxonomic units within A. arenosa, five are tetraploids (Table 1 and Additional file 5). As a source of raw material for natural selection to shape novel genes, this genome duplication may well have contributed to genomic instability, leading to genome rearrangement and a driver of speciation in this group.
Only for the very rare A. umezawana (from the A. halleri lineage) is no chromosome data available, and unfortunately no leaf material was available for microsatellite analysis. Since sequence data (ITS and chloroplast DNA) do not favor any hybrid origin and the various A. halleri segregates are exclusively diploid, A. umezawana probably also represents a diploid taxon. For A. croatica there are diploid and tetraploid chromosome number reports, but the few reports of tetraploids in the field [53] suggest misidentifications as for A. arenosa (given the geographic origins of the samples).
We have provided some historical evolutionary context for many of the non-model lineages that comprise the Arabidopsis genus. ITS data provided the most robust signature to separate the main evolutionary lineages (Figures 2 and 5): 1) Arabidopsis thaliana, 2) A. cebennensis, 3) A. pedemontana, 4) A. lyrata and its segregates/subspecies, 5) A. arenosa with numerous different species and cytotypes and A. croatica more distinct from the remainder, and 6) A. halleri and its subspecies. This summary excludes two hybrid species, namely Arabidopsis suecica and A. kamchatica "bridging" A. thaliana/A. arenosa and A. halleri/A. lyrata, respectively. These taxa will be discussed subsequently, since there is increasing evidence of substantial gene flow over various species and/or ploidy levels [29],[48].
Taxonomy and systematics of Arabidopsis halleriand its relatives
Delimitation of Arabidopsis halleri is still debated among taxonomists. Up to five subspecies have been recognized [8],[9],[15],[18] though two of these, A. halleri subsp. gemmifera (Matsum.) O'Kane & Al-Shehbaz and A. halleri subsp. ovirensis (Wulfen) O'Kane & Al-Shehbaz, are accepted by some authors as separate species, A. gemmifera (Matsum.) Kadota and A. ovirensis (Wulfen) A. P. Iljinsk., respectively [11],[16].
To date, three predominantly central European subspecies were recognized [15]: subsp. halleri, subsp. tatrica (Pawł.) Kolník, and subsp. dacica (Heuff.) Kolník. The third, Asian A. halleri subsp. gemmifera is geographically separated from the other two subspecies [18]. We did not detect these three subspecies here. A. halleri subsp. gemmifera formed a cluster with A. halleri subsp. tatrica. Arabidopsis halleri subsp. ovirensis was originally described as endemic to the eastern Austrian high mountain range at Mount Hochobir (Carinthia). Reports from other localities are most likely based on misidentifications (e.g. from Romania and Ukraine). Unique sequence types (ITS and cpDNA) in the populations from Mount Hochobir are in agreement with this narrow endemic distribution [18],[22]. Based on microsatellite data, A. halleri subsp. halleri is characterized by different distinct genetic clusters, which is in congruence with the multiple A. halleri gene pools shown earlier [31]: here there were two gene pools with admixture between them, and Arabidopsis halleri subsp. dacica did not form a separate genetic cluster. Limited taxon sampling prohibits further interpretation. Although subsp. tatrica did not show genetic distinctiveness in this study, there is "genetic evidence" for the subspecies A. halleri subsp. tatrica[31]. Based on the data presented here, we suggest recognizing five subspecies within A. halleri: gemmifera, tatrica, halleri, ovirensis, and dacica, of which A. halleri subsp. ovirensis is a genetically distinct local endemic taxon and of which A. halleri subsp. tatrica and subsp. dacica need further and detailed phylogeographic analysis. We had limited access to material from A. umezawana, but based on trnLF and ITS data it is closest to the various subspecies of A. halleri.
The evolutionary history of Arabidopsis halleri can be summarized as follows: It has previously been shown that all five subspecies are closely related to each other, and that one major center of genetic diversity is located in the eastern Austrian Alps [18],[22],[31]. The latter [31] explained this center of genetic diversity by secondary contact and admixture of different European gene pools. Similar to the heavy-metal hyper accumulator N. caerulescens[55], it was concluded that A. halleri metallicolous populations were founded independently from non-metallicolous populations without suffering from founder effects [30]. We think that radiation within A. halleri is likely to have occurred during Pleistocene glaciation and deglaciation cycles [22], which also fits with estimates [56] suggesting it to be 335,000 [272,800-438,200] years ago for subsp. halleri. Note that this study lacks other subspecies, so a deeper evolutionary split is possible. Furthermore, microsatellite data suggest that A. halleri subsp. gemmifera may have originated from A. halleri subsp. tatrica from the Tatra Mountains.
Systematics of Arabidopsis arenosaspp. in relation to resolved gene pools
A. arenosa represents a diploid-tetraploid species complex composed of mainly biennial and predominantly outcrossing individuals [53]. The species complex has a distribution range covering most of Eastern Europe and is found in colline to high-alpine habitats exhibiting wide ecological amplitude, spanning from coastal sand dunes to high-alpine screes. Depending on the author, the A. arenosa complex comprises several taxa at various taxonomic levels The complex has been treated as one species, A. arenosa, with two subspecies of partly overlapping distribution ranges in Central Europe [3]: the tetraploid subsp. arenosa, also occurring in northern Europe, growing mainly on siliceous bedrock and sandy soil, and the tetraploid subsp. borbasii, growing predominantly on calcareous bedrock in mountainous regions. Diploid A. neglecta was described mainly from the Carpathians and rarely from the Alps, where its occurrence is doubtful, since in the Alps this taxon was referred to as Cardaminopsis arenosa var. intermedia[57]. However, we clearly show that this taxon is closer to tetraploid A. arenosa subsp. borbasii. Based on morphological and karyological data, several additional (mainly) diploid Carpathian taxa were proposed at the species and subspecies level, and attributed to the genus Cardaminopsis[53],[58]. Some of these names were never published, however, and kept as "nomina provisoria" (nom. prov.) [59]. Taxonomic concepts in the A. arenosa species complex are still strongly debated [18], and we have endeavored to provide clarification here. The lack of resolution for the slower mutating ITS and trnLF regions suggests that (recent) radiation within the Pleistocene is plausible for this species complex (the presence of shared ancestry notwithstanding). Our Structure results distinguish mountainous-alpine tetraploid A. arenosa subsp. borbasii and A. arenosa subsp. arenosa var. intermedia from the remaining tetraploid taxa. Diploid taxa are resolved into A. neglecta subsp. neglecta, A. carpatica, A. petrogena subsp. petrogena, and Arabidopsis nitida. Diploid A. croatica is also well separated and shows clear affinities with the A. arenosa species group as a whole (see below).
The A. arenosa species complex exhibits the highest levels of genetic diversity within the genus. Only A. lyrata subsp. petraea has comparative values here. In tetraploid A. arenosa subsp. arenosa/subsp. borbasii these levels might be explained by (1) local, periglacial survival, (2) lack of genetic bottlenecks and maintenance of large effective population sizes during postglacial migration into formerly glaciated regions, and (3) gene flow between different taxa and/or ploidy levels [19]. In the cases of A. carpatica and A. petrogena subsp. exoleta, the high levels might be an indicator for past and ongoing speciation within the A. arenosa complex in the Western Carpathians [19].
Taxonomy and systematics of A. lyrataand its close relatives
Worldwide, the phylogeography of A. lyrata largely reflects its recent introduction by humans. Three biogeographically defined groups have been recognized: Eurasia, the amphi-Pacific region, and North America [27]. However, the most widely used taxonomy recognizes only two corresponding subspecies (lyrata and petraea), with a third subspecies representing the allopolyploid A. kamchatica[8]. Additional Eurasian taxa such as A. septentrionalis and A. umbrosa have been treated synonymously under A. lyrata subsp. petraea (A. arenicola was at that time treated as a separate taxon) [3]. Our data clearly shows that the North American taxa A. lyrata subsp. lyrata and A. arenicola are close relatives, and that the self-compatible A. arenicola probably originated postglacially from A. lyrata populations [27].
In accordance with the Panarctic Flora taxonomic concept microsatellites recognized two arctic taxa: A. petraea subsp. umbrosa and A. petraea subsp. septentrionalis (Table 1). Both taxa provide a bridge by connecting the European A. lyrata subsp. petraea with the two North American taxa geographically (and genetically). Remarkably, A. petraea subsp. septentrionalis represents a tetraploid taxon and given the high genetic similarity of subsp. umbrosa with subsp. septentrionalis, the latter is most probably an autotetraploid.
Local endemics and hybrid taxa
These parallel traits support an evolutionary vicariance scenario in potential refugia west of the main distribution area of A. halleri, which itself is distributed along the whole alpine mountain chain. The western and relict occurrence of A. pedemontana and A. cebennensis may reflect adaptation to refugia during warming phases, i.e. high-and sub-alpine spring habitats with cool streaming water. Our data did not provide enough power for divergence time estimates, but it seems likely that speciation in A. pedemontana and A. cebennensis occurred during early Pleistocene glaciation and deglaciation cycles.
A. croatica, on the other hand, is morphologically and ecologically much closer to diploid taxa of A. lyrata and A. arenosa. As such, it could be regarded as a derivative of the ancestral gene pool of these respective diploid species (e.g. A. petraea subsp. lyrata, A. carpatica, A. petrogena) (Figure 2) [19].
We did not consider in detail here hybrid taxa such as A. kamchatica and A. suecica. But it is notable that there is increasing evidence of substantial interspecies and interploidal gene flow [48]. It is accepted that A. kamchatica has a multiple polytopic origin [27],[28], and there is increasing evidence that A. suecica does not result from a single hybridization event [60] but rather, multiple events with genetically distinct parents (Polina Novikova, Magnus Nordborg, personal communication), demonstrated and summarized earlier [18]. The first sightings of diploid A. arenosa from the Baltic Sea area is now documented (Filip Kolář, Karol Marhold, personal communication), and future genomic analyses will highlight the relationships with putative parental populations of A. arenosa and A. thaliana. As noted here, A. petraea subsp. septentrionalis is very likely of hybrid origin, consistent with botanical notes - with limited sampling and one Russian population only - which concluded "… (this population) may have originated from a different refugium probably located more in the East" [36]. Genomic analyses of these endemic and hybrid systems will provide further insight into their evolutionary dynamics.
Genome size variation in Arabidopsis
We did not focus in detail on A. thaliana, but as we saw in this study published estimates of nuclear DNA content size in A. thaliana also show some genomic variation among wild accessions (1.1 fold difference with a mean 1C value of 0.215 pg) [61]. Absolute values in pg are discussed critically in the same work, and differ largely from other estimates [62],[63]. Published results from larger-scale studies were obtained by comparing A. halleri and A. lyrata while focusing on allopolyploid A. kamchatica[64], plants from the latter had a slightly smaller genome size than the sum of its diploid parents. The same study could differentiate genome sizes at the subspecies level, by comparing A. kamchatica subsp. kamchatica and subsp. kawasakiana (larger genome). Data from A. kamchatica[64] showed only small differences compared to data from much smaller sample sizes [65] and focusing on A. lyrata. Another large-scale study [66] focused on European A. lyrata and A. arenosa and demonstrated slightly but significantly larger nuclear DNA content in A. lyrata compared to A. arenosa and its segregates (with both ploidy levels). However, in general there is only a limited number of genome size studies within the genus Arabidopsis. Published genome size with the smallest genome size found in A. arenosa (1Cx of 0.2 pg) and the largest genome size observed in A. cebennensis (1Cx of 0.29 pg) confirm our results [63].
Some discrepancies are apparent among published studies when comparing absolute values of genome sizes either given in pg or in Mbp but this is mostly due to deviations in methodology (e.g. different standards, different fluorescent dyes, sample preparation, diurnal variation within a sample) [64].
Taxonomic remarks
We do not formally propose new taxonomical combinations but rather highlight some changes (below) which need to be implemented pending completion of more detailed morphological and phylogeographic analyses.
Arabidopsis arenosa subsp. arenosa var. intermedia
This taxon is best kept as subsp. of A. arenosa, namely A. arenosa subsp. intermedia. This reflects at best that all tetraploid segregates of the A. arenosa group closely belong to each other, but also considers the morphological distinctiveness and local alpine occurrence of A. arenosa subsp. intermedia.
Arabidopsis umezawana
We have no evidence of a hybrid origin (e.g. close affinities to hybrid A. kamchatica), but instead convincing evidence that it falls into the A. halleri group. Consequently the taxon is at best treated as a subspecies, namely A. halleri subsp. umezawana.
Arabidopsis arenicola
This taxon is closest related to North American A. lyrata subsp. lyrata. Consequently, it should be treated as subspecies under a North American A. lyrata, namely A. lyrata subsp. arenicola. Morphological differences are weak: Compared to A. lyrata subsp. lyrata fruits are terete or only slightly flattened, and cotyledons are incumbent [67].
A. lyrata subsp. petraea
Arabidopsis petraea subsp. septentrionalis and A. petraea subsp. umbrosa have already been described and characterized within the Pan Arctic flora project as subspecies within A. petraea. There are two different options to solve this taxonomic/phylogenetic incongruence. 1) European A. lyrata subsp. petraea is treated as A. petraea subsp. petraea, thereby taking geographical and genetic affinities into account and not changing taxonomy of subsp. septentrionalis and umbrosa. 2) Treating members of the A. lyrata group on subspecies level and establishing the new combinations A. lyrata subsp. septentrionalis and A. lyrata subsp. umbrosa. We prefer the second option, since this would minimize future confusion and misuse of species names [21]. Clearly more and detailed studies of these two taxa are needed. From the herein presented microsatellite analysis it could be hypothesized that the A. halleri genome is also introgressed into both species. The Structure analysis shows some affinities with the purple genetic cluster and linking in particular subsp. septentrionalis and umbrosa with some populations of A. halleri and A. pedemontana (Figure 5B, upper part; Figure 5A, lower part), but see comments given above with A. pedemontana.
We characterized in detail the three main Arabidopsis evolutionary lineages: A. halleri, A. lyrata and A. arenosa, including their respective subspecies in an attempt to present a genus-wide overview on genetic variation and taxon delimitation. The relationship among these three lineages is not completely certain due to the power of resolution across the assays used here, but there is some tendency that the lyrata lineage is more closely related to arenosa than to halleri, consistent with being sister taxa. Three additional well-defined endemic species, A. pedemontana, A. cebennensis and A. croatica do form separate evolutionary lineages, with the latter (croatica) most likely positioned at the base of the A. arenosa lineage. The other two endemics are distantly related to any other lineage, but ecologically and morphologically closer to A. halleri. Aside from these evolutionary lineages, there is a need to characterize some taxa in much more detail, such as the arctic taxa of A. lyrata and members of the A. arenosa species aggregate. One other conclusion which stems from the extensive chloroplast haplotype sharing observed among all major evolutionary lineages is the need to qualify and quantify the extent of gene flow within the entire genus.
Plant material and general sampling strategy
This study was designed to incorporate as much existing data as possible in order to provide a comprehensive perspective on taxon sampling as well as their geographic (spatial) distribution. The internal transcribed spacers (ITS1 and ITS2) separating the small and large rRNA subunits and the plastid trnL intron including the adjacent trnL-F intergenic spacer (hereafter called trnLF region). A single individual of respective accessions rather than population-based sampling is common however from these earlier publications. Consequently, many new species accessions and population level sampling have now been added to the existing sample pool. All new individuals have been genotyped using microsatellites using methods established and optimized for the characterization of an Arabidopsis hybrid zone [29].
The different sampling levels (populations versus individual accessions) is also the main reason why ITS and trnLF data are visualized as trees and/or networks (phylogenetics), and microsatellite data were subject to population based algorithms. World-wide sampling localities are provided in Additional file 5 (including GenBank accession numbers) and illustrated in Figure 1. In brief, we sampled 2909 individuals from 813 populations/accessions representing all taxa and cytotypes of the genus Arabidopsis (see also Table 1). For individual marker sets the sampling is as follows: ITS, 1120 individuals/524 accessions; trnLF, 1777 individuals/632 accessions; microsatellites, 1345 individuals/222 accessions; and cytogenetic analysis, 221 accessions. Note that not all sample material was of sufficient quantity or quality for PCR (material included: voucher, wild, living collections). Information on ploidy level (Table 1) is either based on chromosome counts, genome size measurements or indirectly by the numbers of alleles per locus (based on microsatellite genotyping) (see [19],[29] for cytological methods). Unambiguous detection of polyploids using microsatellite genotyping is, of course, only possible if more than two alleles are present at a given locus.
DNA isolation, amplification, and sequencing
Total DNA was obtained from dried leaf material and extracted according to a CTAB protocol [68] 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 [27],[69], using a PTC200 (MJ Research, Waltham, USA) thermal cycler. The PCR reaction volume of 50 μl contained 1x 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 5 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 ITS region was performed according to [70]. 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.8S, 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). The sequencing was performed by GATC GmbH (Konstanz, Germany) and Eurofins MWG Operon (Ebersberg, Germany). Additionally, cycle-sequencing was performed on the MegaBase500 system using the DYEnamic ET Terminator Cycle Sequencing Kit (Amersham Biosciences, Chalfont St Giles, England).
Microsatellite amplification and allele detection
Microsatellites were chosen from previous studies of A. lyrata[29],[71]. The allopolyploid A. kamchatica, A. suecica and introgressed tetraploid hybrids of A. lyrata subsp. petraea and A. arenosa were excluded from this analysis. Selection criteria, PCR and genotyping conditions are provided in detail together with a list of the seven SSRs finally chosen for the analyses in our previous contribution [29]. Scoring of fragment sizes and fluorescence intensity/peak heights (in tetraploids) was automatically performed with GeneMarker version 1.95 (SoftGenetics, State College PH, USA) using respective panels for each locus with subsequent manual checking of each sample. Allele frequencies within tetraploid individuals could unambiguously be assigned manually for the majority of individuals, based on the fluorescence intensity of the fragment peaks [29].
Estimation of nuclear DNA content
Nuclear DNA content was determined using flow cytometry following a simplified two-step protocol [72]. Approximately 10 mM2 of fresh leaf tissue (or one fresh petal) from each plant was chopped together with an appropriate volume of the internal reference standard (Solanum pseudocapsicum, 2C = 2.59 pg, [73]; an identical individual was used for all measurements) using a razor blade in a Petri-dish containing 0.5 ml of ice-cold Otto I buffer (0.1 m citric acid, 0.5% Tween 20). The suspension was filtered through a 42-μm nylon mesh and incubated for 10 min at room temperature. Isolated nuclei were stained with 1 ml of Otto II buffer (0.4 m Na2HPO4.12H2O) supplemented with propidium iodide and RNase (both in concentration 50 μg mL-1), and ?-mercaptoethanol in concentration 2 μg mL-1. After a few minutes, the relative fluorescence intensity of 5000 particles was recorded using flow cytometer CyFlow SL (Partec GmbH, Germany) equipped with green (532 nm) solid state laser. We applied the following stringent criteria in order to get precise and stable flow cytometric results: (i) only analyses with the coefficient of variation of the sample peak below 3% were taken into account (ii) each sample was measured at least three times on different days to minimize potential random instrumental drift [74], and (iii) the between-day variation was defined to not exceed the 3% threshold; otherwise the most remote value was discarded and the sample was re-analyzed. The histograms were evaluated with FloMax FCS 2.0 program (Partec GmbH, Germany). Differences in homoploid nuclear DNA contents among major gene pools (separately for diploid and tetraploid accessions) were analyzed by one-way ANOVA with TukeyHSD post-hoc comparisons in R v.2.15.2 [75]. The dataset comparing relative genome sizes of taxonomic groups was generated in Prague. A second dataset was generated in Heidelberg to provide some estimates on absolute genome sizes. The second dataset incorporated different standards (Solanum lycopersicum cv. Stupicke, 0.98 pg/1C; and Raphanus sativus cv. Saxa, 0.55 pg/1C) [76] because of comparing to, and integrating into datasets from all over the Brassicaceae. Respective data are deposited in BrassiBase[20],[77]. The two datasets were not merged afterwards and kept separate, because accessions analyzed and standards used were different (as explained above).
ITS and trnLF DNA sequence delimitation
Plastidic trnLF sequences were defined as haplotypes and suprahaplotypes following previous studies [18],[22],[26],[29]: Haplotypes are characterized by multiple trnF pseudogenes in the 3′-region of the trnLF-IGS close to the functional trnF gene [26],[78],[79]. When defining respective trnLF suprahaplotypes, we excluded the pseudogene-rich region and thereby merged sets of haplotypes into suprahaplotypes. The trnF pseudogenes evolve with a mutation rate 10 × higher than single nucleotide polymorphisms, which makes them non-applicable for phylogenetic reconstruction at the species level [26],[80],[81]. 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 trnLF haplotypes were assigned to GenBank [LN610052-LN610063/LN610032-LN610051)] (Additional files 3 and 5). ITS sequences were obtained from direct sequencing of PCR products and defined as previously [18],[22],[26],[29]. A few minor corrections of past ITS type numbering had to be conducted, and codes are indicated in Additional file 5 with new assignments to GenBank [LN610064-610098].
Network, phylogenetic analysis and genetic diversity statistics
Network analyses and genetic diversity statistics were exclusively performed using the trnLF suprahaplotypes, as the pseudogene-rich region is not applicable for phylogenetic reconstruction at the species level [26]. The alignment of the cpDNA sequences was manually made with subsequent adjustment in PhyDE version 0.9971 [82]. The network was constructed using TCS version 1.21 [83] and the statistical parsimony algorithm [84]. Gaps (except polyT stretches) were coded as single additional binary characters. Reliability of certain connections, especially if multiple and internal connections occurred within the network, was tested by analyzing the respective alignment with maximum likelihood-based tree construction methods [85]. Only those connections showing up in both types of analyses were retained. Any unsupported connections are indicated with dashed lines in the respective figure. DNA sequence information from A. thaliana was used to set the root.
ITS sequences were also aligned manually with subsequent adjustment in PhyDE version 0.9971 [82]. Maximum parsimony analysis was performed running PAUP 4.0b10 [86] and using A. thaliana as an outgroup. The parsimony heuristic search was performed with the following settings: gaps were treated as missing data (using the gap-based coded 0/1-matrix), multi-state taxa were interpreted as uncertainty; tree construction was via stepwise addition; tree-bisection-reconnection (TBR) was implemented via the branch-swapping algorithm; MaxTrees limit was set to 10,000; and the MulTrees option was selected (saving all minimal trees found during branch swapping). For bootstrapping, 1000 replicates with a tree maximum of 500 retained trees were run. The resulting phylogenetic hypothesis was used to manually place the root in a reliable way with the subsequently performed network analysis (SplitsTree 4.13.1; [87]). For the network analysis A. thaliana was removed from the dataset to increase resolution of internal splits (removing homoplastic characters).
Genetic diversity statistics were performed with Arlequin version 220.127.116.11 [88] and Nei's genetic diversity and gene diversity was calculated accordingly [89]. Allopolyploids (Arabidopsis kamchatica, A. suecica, introgressed A. lyrata subsp. petraea) and individuals which could only be assigned to a lineage but not to lower taxonomic units were excluded from the analyses.
Genotyping of microsatellite alleles and genetic assignment tests
We obtained comparatively full datasets for diploid and tetraploid microsatellite allele scoring (Additional file 6). Microsatellite genotypes were analyzed using Structure 2.3.4 [90],[91], with ten replicate runs for each K-value, and a burn-in period of 1 × 105 and 2 × 105 iterations. The options 'admixture model' was used in combination with `uncorrelated allele frequencies'. The estimation of the optimal K number of populations (ranging from 1 to 10) was calculated using the R-script Structure-sum [92], which compares the posterior probabilities of the runs [93], the similarity coefficient between the runs, and delta K as defined by [94]. In the visualization of Evanno's delta K, a peak had to appear in the optimal fitting model with consistent results over multiple runs [92],[94]. Input files for CLUMPP were generated with STRUCTURE HARVESTER [95], alignments of replicate runs were conducted in CLUMPP [96] and the mean of 10 runs was visualized [97]. Note that for some of the more complicated groupings (e.g. diploids with all species accessions) the variance between independent runs for K with the highest delta K (optimal K according to the method of Evanno [94]) was high. In these cases we turned to the variance for guidance concerning the correct K, and choose K with the lowest variability across runs. At all times we aimed for the smallest value of K that captured most of the structure in the data with a clear biological interpretation for individual assignments.
To overcome conceptual restrictions in combining diploid and tetraploid data we conducted three separate analyses: I) on the whole ploidy dataset, this combined diploids and tetraploids where diploid alleles were doubled to mimic tetraploid data; II) diploids only; and III) tetraploids only. Following these three analyses Structure was again run on the subsets of accessions detected by analysis I to III.) The whole dataset (I) comprised 24 taxa and 1345 individuals and was subsequently split into two separate runs. The diploid dataset (II) comprised 17 taxa and 998 individuals and was subsequently split into three separate runs. The tetraploid dataset (III) comprised seven taxa and 347 individuals and was subsequently split into two separate runs. For those subsets we also tested for optimal K-values. LocPriors were set with split datasets to optimize search strategies using taxon labels. Two aspects have to be considered regarding the Structure analyses: 1) We excluded a priori any known hybrid taxon from the analyses (e.g. A. suecica, A. kamchatica, A. lyrata from the eastern Austrian Forealps and the Wachau in Austria; for details see: [27],[29],[48], and 2) almost all taxa are obligate outcrossers (known self-compatible exceptions are the few populations of A. lyrata in the Great Lakes region of eastern North America: [33],[51]; A. kamchatica and A. kamchatica subsp. kawasakiana from Japan: [54],[98]; A. suecica[99]; A. arenicola, this study). It can also be assumed that A. kamchatica is self-compatible. Genetic diversity statistics were performed with Arlequin 126.96.36.199 [88] for diploid taxa.
Availability of supporting data
The data sets supporting the results of this article are available online. A complete documentation of the new sequences generated for this study, including GenBank accession numbers, is available from Additional file 5. Further, data files (accession list; ITS alignment; microsatellite dataset) are accessible with the Dryad data repository under doi:10.5061/dryad.497sg.
Authors' contributions
NH and RS carried out the molecular marker studies and statistical analyses, and contributed to drafting the manuscript. ML and FK conducted part of genome size measurements and contributed to drafting the manuscript. TJC and KM helped with plant material and contributed to drafting the manuscript. MAK designed and coordinated the project, analyzed and integrated the results and drafted the manuscript. All authors read and approved the final version. With the exception of part of the genome size analysis, most of the work was done in Heidelberg.
Additional files
We thank Ihsan 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, Liza Kretz and Peter Sack for laboratory assistance. We are very grateful to Graham Muir for countless valuable comments and careful editing of the manuscript. This research was supported by DFG grants KO 2302/5 and KO 2302/14 (priority research program DFG-SPP 1529) to Marcus A. Koch and by the Czech Science Foundation grant no. P506/12/0668 to Karol Marhold.
Authors' Affiliations
- Clauss M, Koch MA: Arabidopsis and its poorly known relatives. Trends Pl Sci. 2006, 11: 449-459. 10.1016/j.tplants.2006.07.005.View ArticleGoogle Scholar
- Al-Shehbaz IA, O'Kane SL, Price RA: Generic placement of species excluded from Arabidopsis. Novon. 1999, 9: 296-307. 10.2307/3391724.View ArticleGoogle Scholar
- Al-Shehbaz IA, O'Kane SL: Taxonomy and phylogeny ofArabidopsis(Brassicaceae). In The Arabidopsis Book 2002, Volume 1. Edited by Torii K. The American Society of Plant Biologists; 2002:e0001. doi:10.1199/tab.0001.,Google Scholar
- Koch M, Bishop J, Mitchell-Olds T: Molecular systematics and evolution of Arabidopsis and Arabis. Pl Biol. 1999, 1: 529-537. 10.1111/j.1438-8677.1999.tb00779.x.View ArticleGoogle 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. 10.1093/oxfordjournals.molbev.a026248.PubMedView ArticleGoogle 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.PubMedView ArticleGoogle Scholar
- Karl R, Koch MA: A world-wide perspective on crucifer speciation and evolution: phylogeny, biogeography and trait evolution in tribe Arabideae. Ann Bot. 2013, 112: 983-1001. 10.1093/aob/mct165.PubMedPubMed CentralView ArticleGoogle Scholar
- O'Kane SL, Al-Shehbaz IA: A synopsis of Arabidopsis (Brassicaceae). Novon. 1997, 7: 323-327. 10.2307/3391949.View ArticleGoogle Scholar
- O'Kane SL, Al-Shehbaz IA: Phylogenetic position and generic limits of Arabidopsis (Brassicaceae) based on sequences of nuclear ribosomal DNA. Ann Missouri Bot Gard. 2003, 90: 603-612. 10.2307/3298545.View ArticleGoogle 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.View ArticleGoogle Scholar
- Kadota Y: Arabidopsis umezawana (Brassicaceae), a new species from Mt. Rishirizan, Rishiri Island, Hokkaido, Northern Japan. J Jpn Bot. 2007, 82: 232-237.Google Scholar
- Dorofeyev VI: Cruciferae of European Russia. Turczaninowia. 2002, 5: 5-114.Google Scholar
- Marhold K, Perný M, Kolník M: Miscellaneous validations in Cruciferae and Crassulaceae. Willdenowia. 2003, 33: 69-70.View ArticleGoogle 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
- Kolnik M, Marhold K: Distribution, chromosome numbers and nomenclature conspect of Arabidopsis halleri (Brassicaceae) in theCarpathians. Biologia (Bratislava). 2006, 61: 41-50. 10.2478/s11756-006-0007-y.View ArticleGoogle Scholar
- Iljinska A, Didukh Y, Burda R, Korotschenko I: Ecoflora of Ukraine. 2007, Phytosociocentre Press, KyivGoogle Scholar
- Elven DR, Murray J: New combinations in the Panarctic vascular plant flora. J Bot Res Inst Texas. 2008, 2: 433-438.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
- Schmickl R, Paule J, Klein J, Marhold K, Koch MA: The evolutionary history of the Arabidopsis arenosa species complex: Highly diverse tetraploids mask that the Western Carpathians are the center of species and genetic diversity. PLoS One. 2012, 7: e42691-10.1371/journal.pone.0042691.PubMedPubMed CentralView ArticleGoogle Scholar
- Koch MA, Kiefer M, German D, Al-Shehbaz IA, Franzke A, Mummenhoff K: BrassiBase: tools and biological resources to study characters and traits in the Brassicaceae - version 1.1. TAXON. 2012, 61: 1001-1009.Google Scholar
- Koch MA, German D: Taxonomy and systematics are key to biological information: Arabidopsis, Eutrema (Thellungiella), Noccaea and Schrenkiella (Brassicaceae) as examples. Frontiers Pl Science. 2013, 4: e267-Google Scholar
- Koch MA, Matschinger M: Evolution and genetic differentiation among relatives of Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2007, 104: 6272-6277. 10.1073/pnas.0701338104.PubMedPubMed CentralView ArticleGoogle Scholar
- Castric V, Bechsgaard J, Schierup MH, Vekemans X: Repeated adaptive introgression at a gene under multiallelic balancing selection. PLoS Genet. 2008, 4: e1000168-10.1371/journal.pgen.1000168.PubMedPubMed CentralView ArticleGoogle 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.PubMedView ArticleGoogle 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.PubMedView ArticleGoogle 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.View ArticleGoogle Scholar
- Schmickl R, Jorgenson M, Brysting A, Koch MA: 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. BMC Evol Biol. 2010, 10: e98-10.1186/1471-2148-10-98.View ArticleGoogle 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.PubMedView ArticleGoogle Scholar
- Schmickl R, Koch MA: Arabidopsis hybrid speciation processes. Proc Natl Acad Sci U S A. 2011, 108: 14192-14197. 10.1073/pnas.1104212108.PubMedPubMed CentralView ArticleGoogle Scholar
- Pauwels M, Saumitou-Laprade P, Holl AC, Petit D, Bonnin I: Multiple origin of metallicolous populations of the pseudometallophyte Arabidopsis halleri (Brassicaceae) in Central Europe: the cpDNA testimony. Molec Ecol. 2005, 14: 4403-4414. 10.1111/j.1365-294X.2005.02739.x.View ArticleGoogle Scholar
- Pauwels M, Vekemans X, Godé C, Frérot H, Castric V, Saimitou-Laprade P: Nuclear and chloroplast DNA phy logeography reveals vicariance among European popula tions of the model species for the study of metal tolerance, Arabidopsis halleri (Brassicaceae). New Phytol. 2012, 193: 916-928. 10.1111/j.1469-8137.2011.04003.x.PubMedView ArticleGoogle Scholar
- Tedder A, Hoebe PN, Ansell SK, Mable BK: Using chloroplast genes for phylogeography in Arabidopsis lyrata. Diversity. 2010, 2: 653-678. 10.3390/d2040653.View ArticleGoogle Scholar
- Hoebe PN, Stift M, Tedder A, Mable BK: Multiple losses of self-incompatibility in North-American Arabidopsis lyrata? Phylogeographic context and population genetic consequences. Mol Ecol. 2009, 18: 4294-4939. 10.1111/j.1365-294X.2009.04400.x.View ArticleGoogle Scholar
- Clauss M, Mitchell-Olds T: Population genetic structure of Arabidopsis lyrata in Europe. Mol Ecol. 2006, 15: 2753-2766. 10.1111/j.1365-294X.2006.02973.x.PubMedView ArticleGoogle Scholar
- Kuittinen H, Niittyvuopio A, Rinne P, Savolainen O: Natural variation in Arabidopsis lyrata vernalization requirement conferred by a FRIGIDA indel polymorphism. Mol Biol Evol. 2008, 25: 319-329. 10.1093/molbev/msm257.PubMedView ArticleGoogle 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.PubMedView ArticleGoogle 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.PubMedView ArticleGoogle Scholar
- Leinonen PH, Sandring S, Quilot B, Clauss MJ, Mitchell-Olds T, Agren J, Savolainen O: Local adaptation in European populations of Arabidopsis lyrata (Brassicaceae). Am J Bot. 2009, 96: 1129-1137. 10.3732/ajb.0800080.PubMedView ArticleGoogle 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Savolainen O, Kuittinen H: Arabidopsis lyratagenetics. In Genetics and Genomics of the Brassicaceae. Edited by Bancroft I, Schmidt R. New York: Springer Verlag; 2011:347-372.,Google 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Madlung A, Tyagi AP, Watson B, Jiang H, Kagochi T, Doerge RW, Martienssen R, Comai L: Genomic changes in synthetic Arabidopsis polyploids. Plant J. 2005, 41: 221-230. 10.1111/j.1365-313X.2004.02297.x.PubMedView ArticleGoogle Scholar
- Hollister J, Arnold B, Svedin E, Xue K, Dilkes B, Bomblies K: Genetic adaptation associated with genome-doubling in autotetraploid Arabidopsis arenosa. PLoS Genet. 2012, 8: e1003093-10.1371/journal.pgen.1003093.PubMedPubMed CentralView ArticleGoogle Scholar
- Yant L, Hollister JD, Wright KM, Arnold BJ, Higgins JD, Franklin FCH, Bomblies K: Meiotic adaptation to genome duplication in Arabidopsis arenosa. Curr Biol. 2013, 23: 2151-2156. 10.1016/j.cub.2013.08.059.PubMedView ArticleGoogle Scholar
- Hunter B, Bomblies K: Progress and promise in usingArabidopsisto study adaptation, divergence and speciation. In The Arabidopsis Book 2010, Volume 8. Edited by Torii K. Rockville, MD: American Society of Plant Biologists; 2010:e0138.,Google Scholar
- Hu TT, Pattyn P, Bakker EG, Cao J, Cheng JF, Clark RM, Fahlgren N, Fawcett JA, Grimwood J, Gundlach H, Haberer G, Hollister JD, Ossowski S, Ottilar RP, Salamov AA, Schneeberger K, Spannagl M, Wang X, Nasrallah ME, Bergelson J, Carrington JC, Gaut BS, Schmutz J, Mayer KFX, Van de Peer Y, Grigoriev IV, Nordborg M, Weigel D, Guo YL: The Arabidopsis lyrata genome sequence and the basis of rapid genome size change. Nat Genet. 2011, 43: 476-481. 10.1038/ng.807.PubMedPubMed CentralView ArticleGoogle Scholar
- Koch M, Dobes 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.PubMedView ArticleGoogle Scholar
- Jorgensen MH, Ehrich D, Schmickl R, Koch MA, Brysting A: Interspecific and interploidal gene flow in Central European Arabidopsis (Brassicaceae). BMC Evol Biol. 2011, 11: e346-10.1186/1471-2148-11-346.View ArticleGoogle Scholar
- Ross-Ibarra J, Wright SI, Foxe JP, Kawabe A, DeRose-Wilson L, Gos G, Charlesworth D, Gaut BS: Patterns of polymorphism and demographic history in natural populations of Arabidopsis lyrata. PLoS One. 2008, 3: e2411-10.1371/journal.pone.0002411.PubMedPubMed CentralView ArticleGoogle Scholar
- Mable BK, Schierup MH, Charlesworth D: Estimating the number, frequency, and dominance of S-alleles in a natural population of Arabidopsis lyrata (Brassicaceae) with sporophytic control of self-incompatibility. Heredity. 2003, 90: 422-431. 10.1038/sj.hdy.6800261.PubMedView ArticleGoogle 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. 10.1111/j.0014-3820.2005.tb01794.x.PubMedView ArticleGoogle Scholar
- Roux C, Pauwels M, Ruggiero MV, Charlesworth D, Castric V, Vekemans X: Recent and ancient signature of balancing selection around the S-locus in Arabidopsis halleri and Arabidopsis lyrata. Mol Biol Evol. 2013, 30: 435-447. 10.1093/molbev/mss246.PubMedPubMed CentralView ArticleGoogle Scholar
- Měsíček J: Chromosome counts in Cardaminopsis arenosa agg. (Cruciferae). Preslia. 1970, 42: 225-248.Google Scholar
- Tsuchimatsu T, Kaiser P, Yew CL, Bachelier JB, Shimizu KK: Recent loss of self-incompatibility by degradation of the male component in allotetraploid Arabidopsis kamchatica. PLoS Genet. 2012, 8: e1002838-10.1371/journal.pgen.1002838.PubMedPubMed CentralView ArticleGoogle Scholar
- Koch M, Mummenhoff K, Hurka H: Systematics and evolutionary history of heavy metal tolerant Thlaspi caerulescens in Western Europe: evidence from genetic studies based on isozyme analysis. Biochem Syst Ecol. 1998, 26: 823-838. 10.1016/S0305-1978(98)00057-X.View ArticleGoogle Scholar
- Roux C, Castric V, Pauwels M, Wright SI, Saumitou-Laprade P, Vekemans X: Does speciation between Arabidopsis halleri and Arabidopsis lyrata coincide with major changes in a molecular target of adaptation?. PLoS One. 2011, 6: e26872-10.1371/journal.pone.0026872.PubMedPubMed CentralView ArticleGoogle Scholar
- Hayek A: Flora von Steiermark. Berlin: Verlag von Gebrüder Bornträger; 1908-1914Google Scholar
- Měsíček J: Cardaminopsis. In Zoznam nižších a vyšších rastlín Slovenska - Checklist of non-vascular and vascular plants of Slovakia. Edited by: Marhold K, Hindák F. 1998, VEDA, Bratislava, 395-396.Google Scholar
- Kolník M: Arabidopsis. Chromosome number Survey of The Ferns and Flowering Plants of Slovakia. Edited by: Marhold K, Mártonfi P, Mereda P Jr, Mráz P. 2012, VEDA, Bratislava, 94-102.Google 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. Molec Biol Evol. 2006, 23: 1217-1231. 10.1093/molbev/msk006.PubMedView ArticleGoogle Scholar
- Schmuths H, Meister A, Horres R, Bachmann K: Genome size variation among accessions of Arabidopsis thaliana. Ann Bot. 2004, 93: 317-321. 10.1093/aob/mch037.PubMedPubMed CentralView ArticleGoogle Scholar
- Johnston SP, Pepper AE, Hall AE, Chen ZF, Hodnett G, Drabek J, Lopez R, Price HJ: Evolution of genome size in Brassicaceae. Ann Bot. 2005, 95: 229-235. 10.1093/aob/mci016.PubMedPubMed CentralView ArticleGoogle Scholar
- Lysak MA, Koch MA, Leitch IJ, Beaulieau JM, Meister A: The dynamic ups and downs of genome size evolution in Brassicaceae. Mol Biol Evol. 2009, 26: 85-98. 10.1093/molbev/msn223.PubMedView ArticleGoogle Scholar
- Wolf DE, Steets JA, Houliston GJ, Takebayashi N: Genome size variation and evolution in a allotetraploidArabidopsis kamchaticaand its parents,Arabidopsis lyrataandArabidopsis halleri.AoB PLANTS 2014, 6: doi:10.1093/aobpla/plu025.,Google Scholar
- Dart S, Kron P, Mable BK: Characterizing polyploidy in Arabidopsis lyrata using chromosome counts and flow cytometry. Canad J Bot. 2004, 82: 185-197. 10.1139/b03-134.View ArticleGoogle Scholar
- Jørgensen MH, Ehrich D, Schmickl R, Koch MA, Brysting AK: Interspecific and interploidal gene flow in central european Arabidopsis (Brassicaceae). BMC Evol Biol. 2011, 11: e346-10.1186/1471-2148-11-346.View ArticleGoogle Scholar
- Al-Shebaz IA: Arabidopsis. Flora of North America. 2010, Oxford University Press, Oxford, 447-449.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. x divaricarpa, and A. holboellii (Brassicaceae). Mol Ecol. 2004, 13: 349-370. 10.1046/j.1365-294X.2003.02064.x.PubMedView ArticleGoogle Scholar
- Dobes C, Mitchell-Olds T, Koch M: Intraspecific diversification in North American Arabis drummondii, A. ×divaricarpa, and A. 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.PubMedView ArticleGoogle Scholar
- Clauss MJ, Cobban H, Mitchell-Olds T: Cross-species microsatellite markers for elucidating population genetic structure in Arabidopsis and Arabis (Brassicaeae). Mol Ecol. 2002, 11: 591-601. 10.1046/j.0962-1083.2002.01465.x.PubMedView ArticleGoogle Scholar
- Doležel J, Greilhuber J, Suda J: Estimation of nuclear DNA content in plants using flow cytometry. Nat Protoc. 2007, 2: 2233-2244. 10.1038/nprot.2007.310.PubMedView ArticleGoogle Scholar
- Temsch EM, Greilhuber J, Krisai R: Genome size in liverworts. Preslia. 2010, 82: 63-80.Google Scholar
- Doležel J, Bartoš J: Plant DNA flow cytometry and estimation of nuclear genome size. Ann Bot. 2005, 95: 99-110. 10.1093/aob/mci005.PubMedPubMed CentralView ArticleGoogle Scholar
- R: A language and Environment for Statistical Computing. 2013, R Foundation for Statistical Computing, Vienna, AustriaGoogle Scholar
- Dolezel J, Sgorbati S, Lucretti S: Comparison of three fluorochromes for flow cytometric estimation of nuclear DNA content in plants. Physiol Plantarum. 1992, 85: 625-631. 10.1111/j.1399-3054.1992.tb04764.x.View ArticleGoogle Scholar
- Kiefer M, Schmickl R, German D, Lysak M, Al-Shehbaz IA, Franzke A, Mummenhoff K, Stamatakis A, Koch MA: BrassiBase: introduction to a novel knowledge database on Brassicaceae evolution. Plant Cell Physiol. 2014, 55: e3-10.1093/pcp/pct158.PubMedView ArticleGoogle Scholar
- Koch MA, Dobeš C, Matschinger M, Bleeker W, Vogel J, Kiefer M, Mitchell-Olds T: Evolution of the trnF(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.PubMedView ArticleGoogle Scholar
- Dobeš C, Kiefer C, Kiefer M, Koch MA: Plastidic trnFUUC pseudogenes in North American genus Boechera (Brassicaceae): mechanistic aspects of evolution. Plant Biol. 2007, 9: 502-515. 10.1055/s-2006-955978.PubMedView ArticleGoogle Scholar
- Koch MA, Dobeš C, Kiefer C, Schmickl R, Klimeš L, Lysak MA: Supernetwork identifies multiple events of plastid trnF(GAA) pseudogene evolution in the Brassicaceae. Mol Biol Evol. 2007, 24: 63-73. 10.1093/molbev/msl130.PubMedView ArticleGoogle Scholar
- Schmickl R, Kiefer C, Dobeš C, Koch MA: Evolution oftrnF(GAA) pseudogenes in cruciferous plants.Plant Syst Evol 2008, [doi:10.1007/s00606-008-0030-2],Google Scholar
- Müller K, Quandt D, Müller J, Neinhuis C: PhyDE, Version 0.92: Phylogenetic Data Editor. 2005Google 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.PubMedView ArticleGoogle 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.PubMedPubMed CentralGoogle Scholar
- Stamatakis A: RAxML Version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies.Bioinformatics 2014, doi:10.1093/bioinformatics/btu033.,Google Scholar
- Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony (*and other methods), Version 4. 2002, Sinauer Associates, Sunderland, MAGoogle Scholar
- Huson DH, Bryant D: Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006, 23: 254-267. 10.1093/molbev/msj030.PubMedView ArticleGoogle Scholar
- Excoffier L, Lischer HEL: Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Eco Res. 2010, 10: 564-567. 10.1111/j.1755-0998.2010.02847.x.View ArticleGoogle Scholar
- Nei M: Molecular Evolutionary Genetics. 1987, Columbia University Press, New YorkGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMedPubMed CentralGoogle Scholar
- Hubisz M, Falush D, Stephens M, Pritchard JK: Inferring weak population structure with the assistance of sample group information. Molec Ecol Res. 2009, 9: 1322-1332. 10.1111/j.1755-0998.2009.02591.x.View ArticleGoogle Scholar
- Ehrich D: AFLPdat: a collection of R functions for convenient handling of AFLP data. Mol Ecol Notes. 2006, 6: 603-604. 10.1111/j.1471-8286.2006.01380.x.View ArticleGoogle Scholar
- Rosenberg NA, Pritchard JK, Weber JL, Cann HM, Kidd KK, Zhivotovsky LA, Feldman MW: Genetic structure of human populations. Science. 2002, 298: 2381-2385. 10.1126/science.1078311.PubMedView ArticleGoogle Scholar
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software structure: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.PubMedView ArticleGoogle Scholar
- Earl DA, vonHoldt BM: Structure harvester: a website and program for visualizing structure output and implementing the Evanno method. Cons Genet Res. 2012, 4: 359-361. 10.1007/s12686-011-9548-7.View ArticleGoogle Scholar
- Jakobsson M, Rosenberg NA: CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics. 2007, 23: 1801-1806. 10.1093/bioinformatics/btm233.PubMedView ArticleGoogle Scholar
- Rosenberg NA: Documentation for Distruct Software: Version 1.1. 2007, University of Michigan, MichiganGoogle 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.PubMedView ArticleGoogle 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.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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.


text:
Twins and their 1 Year Old Sister
Posted on October 16, 2013, in Photography and tagged 1 year old photos, 3 year old photos, affordable photographer, Baby Photographer, Children's Photographer, family portrait photographer, Grand Rapids Photographer, PICS Photography, plainwell photographer, Susan McGowan, West Michigan Photographer. Bookmark the permalink. Leave a comment.

subset: commoncrawl


text:
Holy Comeback! Batgirl Returns After 23 Years
(LONDON) -- Batgirl is back! Twenty-three years after she was shot in the back by The Joker, rendering her a paraplegic, DC Comics has decided to bring her back.
Batgirl, portrayed by actress Yvonne Craig in the 1960s television series, will miraculously regain the ability to walk in an issue to be released in September 2011, according to the Guardian newspaper.
Copyright 2011 ABC News Radio
tagged Batgirl, DC Comics in Entertainment General


text:
Share This:
December 23, 2016
Dangal 1st Day (Friday) Collection Early Trends. Huge Numbers Await.
Hello! Like to say something or having a query regarding the post. Then use the comment section below. Thanks!
10 comments:
30cr is excellent start .. tomorrow should be very huge.. and sunday definitely challenge 40..ReplyDelete
I watched last day... dangal is classic... After 100 years also generation can enjoy the movie. .ReplyDelete
Aamir's performance outstanding... he deserves a national award or more (academic).. in his acting we can find a innocent father and a aggressive coach.. Everyone done fantastic. .ReplyDeleteReplies
100% gurantee dangal will not beat sultanDelete
I want to watch again. . But here in dubai no ticket available till weekend over. . All shows are full. . I know in uae-gcc dangal will be record breaker. . I don't know about india( dangal can break the first weekend record or first day or highest single day)... but one thing am sure. .. dangal definitely cross pk world wide. . But dangal need time... dangal has stamina. . It will run more while doesn't enough power. .ReplyDelete
40 Crore on Sunday .. History will be made.. Only AamirReplyDelete
Thank you Ankur ... I watched this movie today in Bangalore Forum PVR at 11:30 show. Occupancy was more than 90%. This is must watch movie for all kind of audience. Entertainment + Inspiration ...ReplyDelete
Dangal is not a best flim of amir khanReplyDelete
It can't be another 3 idiots or pkReplyDelete
Time will tell . Aamir will break his own recordsReplyDelete


text:
The Public Gardens [Trivandrum]
Photographer: D'Cruz, Zachariah
Medium: Photographic print
Date: 1900
Photograph taken about 1900 by the Government photographer, Zacharias D'Cruz of the Public Gardens in Trivandrum. It is one of 76 prints in an album entitled 'Album of South Indian Views' from the Curzon Collection. George Nathaniel Curzon was Under Secretary of State at the Foreign Office between 1895-98 and Viceroy of India between 1898-1905. The idea of establishing a Public Garden in Trivandrum was sanctioned by the Maharaja in 1859 and was to be accomplished under the direction of John Allan Broun, the Government Astronomer. Several interesting botanical specimens were reared in the garden from seeds and cuttings procured from Bangalore and other cities. In 1879, the Curatorship was abolished and the management of the museum and the gardens was placed under a Committee consisting of the British Resident and three European members. In 1897-98, the institutions were placed under an officer called "Director of the Government Museum and Public Gardens". An artificial lake was constructed in the gardens in the late 19th century.


text:
Hello all. It's been awhile since we last posted on this blog.
Some news just in from the world of BalconyTV.
Firstly ...it's BalconyTV's birthday! We are 11 years old. Crikey.
For a trip down memory lane, here is our first ever test show from 11 years ago today - http://www.balconytv.com/videos/balconytv-test
And here was our official launch show. http://www.balconytv.com/videos/balconytv-launch
18,000+ shows later... here we are.
Every year at this time we like to hand out some awards for the best shows from the past 12 months. We have Shows of the Week, and Shows of the Month... so it just makes sense...
And so with a drumroll -
The Winners of the BalconyTV Music Video Award 2017 as selected by the producers of BalconyTV is Spider Cherry for their performance on BalconyTV Miami. https://youtu.be/2TSAaJzPEJc
They join previous winners of the BalconyTV Music Video Award including RoySeven, The Script, Mumford and Sons, Moon and Maia Vidal.
We must give a shout out to The Runners Up of the BalconyTV Music Video Award 2017 as selected by the producers of BalconyTV - Tha Snakes Crew for their performance on BalconyTV Lyon. https://youtu.be/kxoZHP17vGQ
The Winners of the Peoples Choice BalconyTV Music Video Award 2017 as selected by viewers and fans are Colours of Bubbles for their performance on BalconyTV Poznan. https://youtu.be/4ou-KzEJNso
And the The Runners Up who went really close are Spring Tides for their performance on BalconyTV Cork. https://youtu.be/jkbbSyb-_4c
Congratulations to everyone involved. We have published almost 2000 shows in the past 12 months so to be one of the 12 contenders was a great achievement.
The contenders were -
The Snakes Crew - Lyon
Lauren Ruth Ward - LA
Pavvla - Barcelona
They Must Be Crazy - Lisboa
Spring Tides - Cork
Sculpture Music - Chamonix
Spider Cherry - Miami
Sahara Beck - Sunshine Coast
Lawrence - Brooklyn
Katinka - Copenhagen
Noga Erez - Tel Aviv
Colours of Bubbles - Poznan
You can watch all the shows selected here - http://tinyurl.com/mf3qx7m
And that's it. Here's to another great year of BalconyTV!


text:
Shopping for apparels could be the furthermost unenviable suffer for you especially if your choices are precise fixed. It is because one and only few products fit your waist's size. You might have a dry echelon but not a dry waistline. Don't get the impression down because here are many way to craft your choices turns endless. It includes rhythmic exercises, well-balanced diets, punctual asleep habits, and the supreme working class organic structure property supplements.
Body creation supplements go the new way for achieving a sanguine organic structure. It can distribute your unit near the needful minerals and nutrients so that your muscles can get better after doing a enthusiastic athletics. Thus, agility and force are maintained. In addition, suitable unit structure supplements enclose ingredients such as phosphorus, protein, chloride, calcium, vitamins, zinc, and metal.
People are prone to reflect on that some issue found here that is pertaining to Bodybuilding is incorrect. However, part is assured, all that is textual present is true!Post ads:
software to track cell phone location / 3 signs your wife is cheating / mobile spyware delhi / spy key binds / interscope records publicity contact / surveillance services wa
You essential premiere inquire a physician or a gym pedagogue before mistreatment a natural object structure appendage. However, whatsoever supplements are unarbitrary by doctors to usher a constant medicinal drug to immoderation quite a few diseases. It is also obliging to read reviews about the goods to know its advantages and disadvantages. Take for occasion the Ephedra physical structure place enlarge.
It is singular if you breakthrough whichever utilisation for the matter delineate here on Bodybuilding that we will consciousness the hard work put in writing on Bodybuilding fruitful. So engineer bully employment of it!
Ephedra is one the utmost polemic body place supplements ever released. It was disclosed by the Chinese from a complex Ma Huang more than than 5 a thousand geezerhood ago. It became the base of the Chinese standard medication. Ephedra is utilized to extravagance flu, fever, colds, and some other corresponding symptoms. Recently, whatsoever investigating disclosed that Ephedra can heighten metabolism, advance fat scorching and weight loss, restraint appetite, present excretion for relieving water retention, and unwind the lung airways to assist in treating cough and respiratory disease.Post ads:
affaires zemmour / detective agency cheats / find out if your girlfriend is cheating / cheater dating site / ptr record reverse lookup / record calls app on iphone
Once you are done language what is scripted present on Bodybuilding, have you considered recollecting what has been typewritten and characters them down? This way, you are obligated to have a better-quality considerate on Bodybuilding.
However, the tenderloin personal property of Ephedra specified as occasional beat of the bosom and multiplied humour pressures led the Food and Drug Authority (FDA) to ban its use in April 2004. The reactions of consumers were an assortment of. But Ephedra banning does not affect the old-time flavouring remedies of the Chinese and seasoner tea products containing such atmospheric condition. Likewise, the use of imitative alkaloid and bronchodilator supplements.
This article serves as a sales rep for the implication of Bodybuilding in the room of fluency. Let it equal comprehension all right.
It afraid the FDA because some Ephedra users try to brainwave distance to duck the edict. Since, past one Ephedra somebody products invaded the New York marketplace. The system officials over again called out for prohibition these products. But they were silenced by Judge Tena Campbell who preferred a expand cast in Utah. The FDA ban was challenged by Nutraceutical Corporation claiming that Ephedra is locked to use as tested by its protracted utilisation in the precedent of old-fashioned medicines. FDA never regulates Ephedra as a substance but a genus of remedy. And because FDA failed to typify that use of Ephedra was unsafe, Judge Campbell lifted the Ephedra ban permitting its income for 10 milligrams or less. So, the use of Ephedra was legalized.
A naturopathic physician, Dr. Cathy Wong even welcomes the new Ephedra regnant. Deaths attributed to Ephedra are uncalled-for and rough compared to medical specialty pills fatality taxation.
Learning almost property is what we are flesh and blood present for now. So try to get to cognise as more almost everything, with Bodybuilding whenever mathematical.
Ephedra products can assistance relatives loss unnecessary weights provided that doses are unbendingly followed. Although FDA aerated it as a dangerous remedy that killed lots of ethnic group and transmitted them to poorness due to addiction, much and more grouping are fetching Ephedra for losing weight. And they have established that no remaining herbs are as effective as Ephedra.
We were a bit conditional when embarking on this project on Bodybuilding. However, exploitation the sandstone and force we have, we have make some fine linguistic process bits and pieces on Bodybuilding.
The puzzle of tons society is resolved. They can now have without end choices for apparels that they poorness. With the backing of Ephedra supplements, they achieved the thing improved they poverty.
Please don't nutriment this as an middle hunk of message on Bodybuilding. A lot of energy and thorny pursue has been put to get this end product!


text:
Amritsari Gobi Pakoda
Gobi aka Caulifower Pakoda or Bajo(it's the amchigele version of pakoda) is one of our favourite deep fried snack/starter and when something innovative or say tempting pictures are posted in FB then who can resist. Today's recipe is one of those. Amritsari Gobi Pakoda, name itself is tempting isn't it? Do try it and enjoy.
Recipe Source : Authentic Food Delights
Here is what we need:
2 cups of cauliflower cut into florets
2 onions finely chopped
2 tbsp Coriander leaves finely chopped
Generous pinch of Ajwain
1 tbsp grated ginger
2 tbsp corn flour
3 tbsp Rice flour
5 tbsp Gram flour
4-5 Green chillies minced or finely chopped
1/4 tsp haldi/turmeric powder
Salt to taste
Oil for deep frying
In a large mixing bowl add chopped onion, washed cauliflower florets, salt and mix well add coriander leaves, ajwain, ginger, minced green chillies and mix well.
Now add corn flour, rice flour, gram flour/besan, turmeric powder and give a good mix add very little water and give a good mix it will become sticky. Don't make it like we make for other pakoda or bajji i.e. mixture should be dry and not watery.
Heat oil when it hot then reduce the flame to medium sim and deep fry.
Serve with Dip of your choice and Enjoy!
Kids enjoyed with Ketchup and we had it just like that.
Do try and let me know how you liked it.
The original recipe has used 2 tbsp cornmeal but I haven't because I did not have it in my pantry instead I added 1 tbsp gram flour and 1 tbsp rice flour.
Happy Cooking,


text:
Our product line includes a wide range of nut formers. We also provide operator training in our factory, tooling supplies, special nut tooling designs, technical support, and excellent after sales service.
In the future, we will continue to dedicate our full efforts to improve and meet the needs of customers in the fastener industry. Your satisfaction is our top priority.


text:
My blog friend Moe got me to thinking about the dismal performance of the media when it comes to political and policy coverage. It seems that most media outlets are capable of little more than trying to predict the score — how many Republicans, how many Democrats. Finding thoughtful, well-reasoned, and, particularly, objective analysis of this nation's needs, and how the policy ideas of the different parties might address those needs is hard to come by. TV and other outlets are filled with opinionators, not reporters, and certainly not analysts. Well, other than keeping my subscription to the New Yorker up to date, I don't have any good answers for that.
But I do have an idea on how the media companies can atone for their lack of useful information. They could take a stand against the other great corrupter of our political processes, distortions in campaign and advocacy advertising. Here is my idea. Suppose a media company establishes a baseline fee structure for political advertisements that is, say, three times the usual rate. However, political advertisers could qualify for lower rates by having their ads rated by politifact.com or factcheck.org. The better the rating, the lower the price — all the way back down to the normal rate for an ad that rates completely true.
Sure, there would be lots of details to work through. Obviously such an idea would work much better if all the media in a given market adopted the same policy. Otherwise, a candidate could just forego paying $2000 to run a distorted spot on Channel 2, and run it for $800 on Channel 5 instead. There would have to be some worries that, once they had such power over elections, Politifact and FactCheck might lose their bearings. It happens. But hey, I'm just the big picture guy here. I'm sure there are media-savvy types who could figure out something.
What's appealing about this is that it puts a price on political distortion. Sure, you can still run your ad that implies your opponent once voted to permit snipers to shoot disabled widows on sight — but it will cost you. And that's what we do not have now. Currently, gross distortion in political advertising brings only benefits, but no costs. I heard recently that of the first 58 ads analyzed by Politifact, only two got their best rating. This multi-tiered pricing idea imposes a cost on distortion. Distort all you want, but your budget won't stretch nearly so far. And no one's free speech rights are compromised.
I doubt the various Executive Vice Presidents of Finance for the media companies will rush to embrace my idea. They will worry that this will just cost the company revenue. Maybe. But maybe not. Maybe the compulsion to advertise — and even to distort — is so strong that actually the media companies will generate more revenue. And you know what? If some campaign or some advocacy group decides not to use a media outlet adopting this policy, then the media company should ask itself what it was doing permitting itself to participate in the crass corruption of democracy in the first place. Honor still counts. Or should.