The origin of populations of Arabidopsis thalianain China, based on the chloroplast DNA sequences
© Yin et al; licensee BioMed Central Ltd. 2010
Received: 30 May 2009
Accepted: 8 February 2010
Published: 8 February 2010
In the studies incorporating worldwide sampling of A. thaliana populations, the samples from East Asia, especially from China, were very scattered; and the studies focused on global patterns of cpDNA genetic variation among accessions of A. thaliana are very few. In this study, chloroplast DNA sequence variability was used to infer phylogenetic relationships among Arabidopsis thaliana accessions from around the world, with the emphasis on samples from China.
A data set comprising 77 accessions of A. thaliana, including 19 field-collected Chinese accessions together with three related species (A. arenosa, A. suecica, and Olimarabidopsis cabulica) as the out-group, was compiled. The analysis of the nucleotide sequences showed that the 77 accessions of A. thaliana were partitioned into two major differentiated haplotype classes (MDHCs). The estimated divergence time of the two MDHCs was about 0.39 mya. Forty-nine haplotypes were detected among the 77 accessions, which exhibited nucleotide diversity (π) of 0.00169. The Chinese populations along the Yangtze River were characterized by five haplotypes, and the two accessions collected from the middle range of the Altai Mountains in China shared six specific variable sites.
The dimorphism in the chloroplast DNA could be due to founder effects during late Pleistocene glaciations and interglacial periods, although introgression cannot be ruled out. The Chinese populations along the Yangtze River may have dispersed eastwards to their present-day locations from the Himalayas. These populations originated from a common ancestor, and a rapid demographic expansion began approximately 90,000 years ago. Two accessions collected from the middle range of the Altai Mountains in China may have survived in a local refugium during late Pleistocene glaciations. The natural populations from China with specific genetic characteristics enriched the gene pools of global A. thaliana collections.
Arabidopsis thaliana (L.) Heynh is an annual weed belonging to the family Brassicaceae (Cruciferae). The species is native to Europe and Central Asia, but is now widely distributed in the Northern Hemisphere ranging from 68°N (northern Scandinavia) to Equator (mountains of Tanzania and Kenya) . Many characteristics, from morphological traits to protein and DNA markers, have been used to evaluate natural genetic variation among populations, and to reconstruct an intraspecific phylogeny, for A. thaliana (for example, [2–9]). It has been found that many nuclear genes comprise two or more major differentiated haplotypes, generally referred to as allelic dimorphism [10–20]. Balancing selection or ancient population subdivision was often invoked to explain the pattern. The major mechanisms for balancing selection are heterozygote advantage, frequency-dependent selection, or environmental heterogeneity. It is well known that A. thaliana has an inbreeding mating system. The estimated outcrossing rate of the species is 1% or less . It seems difficult to imagine that so many loci in A. thaliana have experienced balancing selection via heterozygote advantage . Therefore, frequency-dependent selection and/or diversifying selection might be the driving forces for the dimorphism phenomenon, as in the case of pathogen resistance (R) genes [17, 18, 23]. It is not clear yet if the dimorphism also exists in the chloroplast genome.
The chloroplast genome of A. thaliana is a circular DNA composed of 154,478 bp with a pair of inverted repeats of 26,264 bp separated by small and large single-copy regions of 17,780 bp and 84,170 bp, respectively . The uniparentally inherited chloroplast genome has been utilized in many studies in plant population and evolutionary genetics. However, studies focused on global patterns of cpDNA genetic variation among accessions of A. thaliana are scattered. In an investigation on the maternal origins of A. suecica, 12 cpDNA regions were sequenced for 25 A. thaliana accessions, which were mainly collected from Scandinavia . These authors found considerable variation existed among the non-coding single-copy sequences in the chloroplast genome of A. thaliana. In another study, the trnL-trnF cpDNA intergenic spacer region of 475 individuals from 167 A. thaliana populations in its native range was sequenced and 16 haplotypes were identified . Based on the chloroplast and nuclear DNA sequence data, Beck et al proposed the Caucasian area as the possible ancestral area of A. thaliana, and suggested four possibilities for the origin of East Asian populations. They also found that the maternal components of A. suecica shared a high similarity to those in the Asian metapopulation of A. thaliana, especially to those from China .
In the studies incorporating worldwide sampling of A. thaliana populations, the samples from East Asia were very scattered. He et al conducted a study on the genetic diversity of 19 natural Arabidopsis thaliana populations in China based on ISSR and RAPD makers, and found that about 42-45% of the total genetic variation existed within populations and there was a significant correlation between geographic distance and genetic distance . However, the phylogenetic relationships of Chinese populations with those distributed in other regions of the world, and the history of population dispersal in this region, are not clear.
The goals of the present survey are: (1) to examine global patterns of cpDNA genetic variation in A. thaliana; (2) to infer phylogenetic relationships among A. thaliana accessions from all over the world based on cpDNA sequence data, with particular focus on Chinese populations; and (3) to discuss the possible origin(s) of the Chinese populations. It was found in this study that dimorphism did exist in the chloroplast genome of A. thaliana; the 77 accessions studied were grouped into two major clusters; and the Chinese populations might have two independent origins.
Nucleotide variation in the chloroplast DNA sequences
List of the A. thaliana accessions used in this study
Cape Verde Island
Geographic information for the 19 accessions collected from China
Latitude and longitude
47°46' 72" N 88°20' 64" E
46°48' 72" N 90°20' 39" E
29°47' 41" N 106°28' 64" E
29°49' 40" N 106°03' 38" E
32°43' 30" N 105°07' 21" E
27°56' 64" N 108°36' 49" E
28°31' 23" N 110°43' 13" E
29°24' 45" N 110°26' 33" E
32°55' 93" N 107°12' 65" E
33°08' 82" N 106°44' 72" E
30°27' 80" N 117°17' 79" E
30°42' 89" N 116°15' 33" E
31°16' 70" N 115°01' 11" E
30°31' 32" N 114°28' 16" E
32°03' 12" N 118°49' 93" E
26°44' 78" N 114°17' 98" E
26°59' 18" N 116°14' 51" E
29°05' 02" N 120°25' 65" E
29°32' 13" N 119°29' 61" E
Nucleotide diversity (π) and the results of neutral mutation hypothesis tests for the 11 fragments data sets
Fu and Li's D*
Fu and Li's F*
Fu and Li's D
Fu and Li's F
Less frequent nucleotide polymorphisms (such as singleton or doubleton) were in excess for the sequenced regions. Singletons were found at a very high frequency: 44 among the 95 SNPs and 10 among the 26 indels were singletons (Figure 2). The excess of low-frequency polymorphisms resulted in negative Tajima's D, Fu and Li's D* and F* values for most of the sequenced segments; for example, 10 out of 11 Tajima's D values, nine out of 11 Fu and Li's D* values and 10 out of 11 Fu and Li's F* values were negative (Table 3). The values for the combined data matrix were -1.17234 (Tajima's D, P > 0.10), -2.36692 (Fu and Li's D*, P < 0.05) and -2.25760 (Fu and Li's F*, 0.10 > P > 0.05, critical). When Fu and Li's D and F tests were conducted using the A. arenosa ortholog as the reference sequence, similar results were obtained: eight out of 11 Fu and Li's D values and nine out of 11 Fu and Li's F values were negative and for the combined data matrix; both the values were negative (-2.06178 and -2.00798, respectively, 0.10 > P > 0.05; Table 3).
In total, we identified three types of nucleotide variations among the aligned sequences: SNPs, length polymorphisms (including indels), and two small fragmental inversions.
Phylogenetic relationships among the accessions
Although there was incongruence between the two phylogenies, the topological relationship was relatively stable among a large number of accessions. We identified four stable branches in both trees (A, B, D, and E in Figures 4 and 5). The only major difference was in branch C. It was placed within one of the two clusters in the NJ tree, but formed a deep polytonous branch in the MP tree, containing the same accessions except Cvi-0, an accession from Cape Verde Island. The branches A-E comprised 61 out of the 77 accessions.
Discrimination of two major differentiated haplotypes among A. thalianaaccessions
Estimated divergence time for MDHC-I and MDHC-II, and demographic expansion of a monophyletic group of accessions in Asia
The K value between A. arenosa and 77 A. thaliana accessions was 0.0280 ± 0.0037. Using Equation 1, the substitution rate per nucleotide site per year for the sequenced chloroplast regions was 2.8 × 10-9. The K value between MDHC-I and MDHC-II was 0.0022 ± 0.0005. Therefore, the estimated divergence time for MDHC-I and MDHC-II was estimated (using Equation 2) to be about 0.39 ± 0.09 mya.
Although no significant correlation was detected between geographic origin and genetic distance, the 17 accessions collected along the Yangtze River, China, were always clustered together with Kas-2, an accession from Kashmir (74°E, 34°N) in both NJ and MP trees (Figures 4 and 5). In the network analysis, they congregated closely and formed a distinct cluster (Cluster A) in Haplogroup 1 (Figure 8). The level of nucleotide polymorphism of the 18 accessions was very low (π = 0.00030), only six haplotypes (h) were detected (five specific in the 17 accessions along the Yangtze River) and haplotype diversity (Hd) was 0.778. In comparison, the values of π, h and Hd for the 77 accessions were 0.00169, 49 and 0.977, respectively.
To test the model of demographic population growth in the region from which the 18 accessions were sampled, especially for the 17 accessions along the Yangtze River, a mismatch distribution analysis was conducted. Each small fragment inversion was treated as a SNP in the analysis. The SSD between the observed and expected mismatch distribution was 0.093 (P = 0.062) and HRag was 0.300 (P = 0.022). There was an only marginally significant difference in SSD between the observed and the predicted pairwise difference distribution under the sudden-expansion model. This result provided evidence of rapid population expansion along the Yangtze River. The average τ-value was 4.521 (95% confidence intervals: 0.684~8.197). The initial time when the populations expanded along the Yangtze River were calculated using Equation 4 to obtain u (2.52 × 10-5), and then using Equation 3 to obtain t (0.897 × 105). Therefore, the initial time of expansion was estimated to be about 90,000 years ago.
The level and pattern of nucleotide variation in the sequenced chloroplast regions
The π value of the sequenced chloroplast regions among global samples of A. thaliana accessions was 0.00169, which is about one-quarter of that of the mean nucleotide diversity of the nuclear genes in A. thaliana , but double that in another study by Sall et al , in which 12 non-coding single-copy cpDNA regions were sequenced for 25 A. thaliana accessions (π = 0.00061). The differences may be due to the different sampling strategies. The 25 A. thaliana accessions in the latter study were mainly collected from Scandinavia, whereas the 77 accessions in the present study were sampled worldwide. For a highly self-fertilizing species, geographical structure may play an important role on a smaller scale in the level of polymorphism, at least for the uniparentally inherited chloroplast genome. For example, the π value reduced to 0.00030 if only 18 accessions in branch A (Kas-2 and 17 accessions along the Yangtze River) were considered.
Inversions in the chloroplast genome exist in monocotyledonous plants and the Asteraceae. The length of these inversions range from 0.5 to 28 kb, and all have phylogenetic implications [26, 27]. The length of the inversions found in the present study were much shorter, only about 18-20 bp. The accessions with inversions were found mostly scattered on branches B, D, E in the NJ and MP trees. The exception is in branch A, where all accessions had inversion 2. The mechanism responsible for these inversions is not known, but they might have originated several times during the population expansion process. Therefore, it is advisable not to consider them for phylogenetic analysis.
Dimorphism in the chloroplast DNA of A. thaliana
Two significantly differentiated haplotype classes could be identified in the sequenced chloroplast DNA regions, just as in the allelic dimorphism found in some nuclear DNA sequences of A. thaliana. At least three different interpretations have been proposed to explain the nuclear dimorphism phenomenon. First, balanced polymorphisms were usually the mutations maintained in populations by natural selection through heterozygotic advantage [17, 28]. The chloroplast genome is maternally inherited in A. thaliana, and the DNA regions selected for analysis in this study are intergenic regions. Therefore, the dimorphism found in the chloroplast may not be caused by balancing selection via heterozygotic advantage. Furthermore, in our investigation, the value of the Tajima's D-value was negative. A negative Tajima's D value is a general feature of the Arabidopsis thaliana genome , and is correlated to demographic factors, such as population growth , rather than non-neutral forces such as selection .
A second explanation for the nuclear dimorphism is that introgression might result in the allelic dimorphism. Chloroplast DNA introgression has been widely reported [e.g., [31, 32]]. In this study, we found that two related species, Olimarabidopsis cabulica and A. arenosa, had all five identical nucleotide site variations with MDHC-II of A. thaliana, which were the 'markers' to separate MDHC-II from MDHC-I. However, the K values between O. cabulica and the 77 accessions of A. thaliana, and between A. arenosa and the 77 A. thaliana accessions, are 0.0395 and 0.0280, respectively, whereas that between the two MDHCs of A. thaliana is only 0.0022. The interspecific genetic distance is at least one order of magnitude higher than intraspecific genetic distances. The estimated divergence time between O. cabulica and A. thaliana is about 10~14 mya, and that between A. arenosa and A. thaliana is about 3.0~5.8 mya . The genetic distances based on cpDNA between these species pairs correlate to their nuclear gene-based estimated divergent time. These results indicated that the dimorphism in cpDNA found in this study was not the result of recent introgressive hybridization events, but we cannot rule out the possibility that the dimorphism might be the result of ancient introgression events. Hybridization between A. thaliana and its closely related species does occur in nature. For example, several studies confirmed the allotetraploid species, A. suecica, resulted from a hybridization event between A. thaliana and A. arenosa about 10,000 to 50,000 years ago [e.g., [25, 34, 35]].
The third explanation for genetic dimorphism is demographic factors, such as founder effects. Being a small annual weed, A. thaliana is a poor competitor in dense vegetation whereas the highly self-fertilizing characteristic makes it capable of founding a population even from a single seed. As a result, this species has a tendency for rapid colonization and extinction cycles [1, 7, 36]. Founder effects might have occurred repeatedly in the evolutionary history of A. thaliana. The founder event(s) could enable some rare alleles to spread into additional populations when the founder population expanded rapidly if the unoccupied ecological niches were favourable. The divergence time between MDHC-I and MDHC-II was estimated to be about 0.36 mya based on our cpDNA data. This is earlier than the estimated time of demographic expansion during the Eemian interglacial (about 0.122 mya; ). Therefore, another possible explanation for the cpDNA dimorphism might be a founder effect followed by limited gene flow during late Pleistocene glaciations and interglacial periods.
As the accessions in MDHC-II share five specific variable sites with A. arenosa, the chloroplast genomes in MDHC-II might represent more ancient types than those in MDHC-I. It is also supported by the fact that more haplotypes are found in MDHC-II (28) than in MDHC-I (21).
Origin of Chinese populations
The 26 accessions of A. thaliana from Asia included in this study are scattered compared to those collected from Europe. Six of them belong to MDHC-II and 20 belong to MDHC-I. Of the MDHC-II group, two collections from China (XJalt and XJqhx) and two from Kazakhstan (9481 and Kz10) are within or very close to the Altai Mountains. Although these four accessions were not clustered in the same clade in the phylogenies, the two accessions from China were always on the same branch. One of the Kazakhstan accessions (Kz10) was clustered with XJalt and XJqhx together with a Russian accession (N1, from Europe) in the NJ tree (Figure 4). XJalt and XJqhx are unique in that they share six specific variable sites (Figure 2), the most number of specific variable sites in this study. The provenances of these two accessions are about 115 km apart and located in the middle of the Altai Mountains range. Based on the cpDNA data, the populations on the Altai Mountain range may have dispersed there during one of the late Pleistocene glaciations, and some local habitats along the southern slopes of the Altai Mountains might have served as refugia. In contrast to some refugia in Europe, where A. thaliana populations had contributed the postglacial colonization of western and northern Europe , some populations in the Asian refugia, such as XJalt and XJqhx, became relatively isolated genetically from other populations after glaciers retreated. Therefore, some fixed mutations were accumulated specifically in these populations. It is also noticed by Beck et al  that "differentiation is particularly strong between the Central Asian refugia and all others, suggesting that either small historical population sizes in Central Asia, relatively limited gene flow between Central Asia and other areas, or a combination of the two have produced relatively strong genetic drift in this region". The other two Asian accessions in the MDHC-II group are from Japan (Tsu-0) and India (Kas-1). Although Tsu-0 was clustered with Alc-0 from Spain, the bootstrap support is poor in both NJ and MP trees (68% and 52%, respectively). Tsu-0 has two specific variable sites, but only shares one specific site with Alc-0. Its phylogenetic relationships with other accessions are unclear and more samples are needed to clarify its origin. Kas-1 from India shares four specific variable sites with four accessions from Europe (Figure 2); this accession will be discussed later.
In MDHC-I, 17 accessions collected along the Yangtze River, China, were clustered together with an accession from India (Kas-2). All data suggested that the chloroplast genomes of these 18 accessions originated from a single common ancestor. The initial time of expansion of the populations along the Yangtze River was estimated to be about 90,000 years ago based on the cpDNA sequences. At least we can rule out the possibility that these populations were introduced by recent human activities, like the populations found in the USA. Both highlands of the western Himalayas and Caucasus have been proposed as possible ancestral areas of A. thaliana [8, 37]. The populations along the Yangtze River could have dispersed eastwards from the ancestral area to the present-day locations via the Himalayas or Kunlun Mountains, which have an east-west trend. Kas-2 was collected in Kashmir, India (74°E, 34°N), on the southern slopes of the Himalayas. It shares a most frequent haplotype with seven Chinese accessions, but it can be distinguished from them only by one specific variable site (Figure 2). There are at least two explanations for the observed distribution pattern. The first is that Kas-2 represents the type that is most similar to the ancestor of the Yangtze River populations, and the expansion of the latter was immediately or very shortly after the dispersal event. However, in the phylogenies, Kas-2 is not basal among the 18 accessions. The second explanation is that Kas-2 might have been introduced to Kashmir from China. It is uncertain which explanation is more likely. More samples from the western part of the Himalayas and Kunlun Mountains are definitely needed for identification of the possible ancestral haplotype of the Yangtze River populations.
The other three accessions in MDHC-I are H Tsukaya from Japan, Sorbo from Tajikistan, and Kas-1 from India. The former shares one haplotype with four accessions from Europe and one accession from the USA. It is apparent that it represents a recent introduction to Japan, possibly by human activity. The accession from Tajikistan has a specific haplotype that is characterized by four unique SNPs (Figure 2), and it shares a haplotype with one accession from Sweden, Italy and Libya, respectively. Kas-1 shares the exact haplotype with two accessions from Russia (Wassilewskija, Ms-0) and one from Austria (Pi-0). Although Kas-1 and Kas-2 were collected at the same locality (74°E, 34°N), as stated on the ABRC seed stock centre website, their genetic distinctness has been noted previously . Re-sampling from their original locality might help in clarifying the confusion.
The possible maternal parent of Arabidopsis suecica
Arabidopsis suecica is an allotetraploid species whose maternal parent is A. thaliana. The species originated more than 20,000 years ago from an allotetraploid hybrid between A. thaliana and A. arenosa [25, 35]. Several studies indicate that the ancestral area of A. suecica is in Europe [e.g., ]. However, a recent study by Beck et al  revealed a genetic similarity between A. suecica and Chinese accessions of A. thaliana. In the present study, the cpDNA sequences of A. suecica were more similar to those of some European accessions, especially Chi-0 from Russia (34°E, 54°N). Based on the comparison of cpDNA in this study, it is most likely that the maternal parent of A. suecica was from Europe.
Elucidating the dispersal of A. thaliana within Asia is a very complicated issue. Temperature fluctuations during glaciations and interglacial periods in the late Pleistocene, the complex mountain ranges in the Central and Central-East parts of Asia, and human activities have all contributed to the present-day distribution patterns. It is clear from this study that some populations in East Asia, such as those along the Yangtze River, were dispersed there along the Himalaya or Kunlun Mountain ranges and underwent rapid expansion about 90,000 years ago. Altai Mountains may have provided refugia for A. thaliana during Pleistocene glaciations. However, more samples from the Altai, Kunlun, and western part of the Himalaya mountain ranges are needed to elucidate more fully the dispersal history of A. thaliana populations in Asia.
Plants material and DNA extraction
Seventy-seven A. thaliana accessions were used in the survey, seeds of 50 accessions were obtained from the Nottingham Arabidopsis Stock Centre (NASC; University of Nottingham), eight accessions from the Arabidopsis Biological Resource Center (ABRC; Ohio State University) and 19 accessions were field collected in China (Table 1). All sampling locations in China were separated by at least 50 km, with most of the locations separated by more than 300 km (Figure 1, Table 2). The nomenclature of the Chinese accessions follows that of He et al. . Two taxa closely related to A. thaliana, i.e., A. arenosa and A. suecica, and one taxon closely related to the genus Arabidopsis, i.e., Olimarabidopsis cabulica, were used as out-groups. Seeds of selected plants were sterilized and germinated on MS solid agar medium, and then transplanted into soil, as described previously . The seedlings were grown in a controlled environment at 22°C with a photoperiod of 16 h light/8 h dark. About 1 g of green leaf material was harvested from a single plant for each accession and used for total cellular DNA extraction using the CTAB method .
PCR amplification, DNA sequencing, and sequence alignment
Information about eleven pairs of primers for PCR amplification and the length of their products
GenBank Acc No
Intron in rpl16
Two accessions (AHyxx and Abd-0) were randomly chosen to detect intra-accession polymorphism. The seeds of AHyxx were field collected from at least 30 individual plants in Yuexi County, Anhui Province in southeast China . The seeds of Abd-0 were obtained from the ABRC but were originally collected in the UK. Green leaf material was harvested from 10 individual plants for both accessions and total cellular DNA was extracted separately. PCR amplification, purification and sequencing were conducted as described above.
The DNA sequences were aligned initially using MEGA 4.1 . After alignment, the sequences were edited manually. Gaps were positioned to minimize nucleotide mismatches. The 3'- and 5'- flanking regions of the protein-coding sequences were trimmed off and only the intergenic spacer or intron regions were retained for further analysis.
The nucleotide sequences of 11 fragments from each accession were concatenated sequence by sequence to form a combined sequence set. The combined sequence sets formed a combined data matrix of 77 × 8978 nucleotides.
Population genetics and phylogenetic analyses
Polymorphism analyses were conducted for both intra- and inter-specific comparison using DnaSP4.50 . Nucleotide variation was estimated as nucleotide diversity (π) . The number of haplotypes was denoted by h and haplotype diversity was denoted by Hd. Some statistical methods, such as Tajima's D test  and Fu and Li's D*, D, F* and F test , were conducted to test the neutral mutation hypothesis. The A. arenosa orthologs were used as the out-group sequences when needed.
Phylogenetic analyses were performed on the combined data matrix. Both neighbor joining (NJ) and maximum parsimony (MP) methods were used. The NJ tree was constructed using MEGA 4.1, and Maximum Composite Likelihood distance  was used and gaps were treated as a complete deletion. The analysis was done with 1000 bootstrap replicates of the data. The MP tree was constructed using PAUP*4.0 . Bootstrap values were calculated from 1000 replicates with a heuristic search strategy. Gaps in the alignments of the sequences were treated as missing data.
To visualize the relationship among haplotypes, a network was constructed as described by  using NETWORK version 126.96.36.199. A default weight of 10 was applied to each site. 'Epsilon = 0' was chosen for constructing the network for the 77 A. thaliana accessions using the median joining algorithm.
Estimated divergence time and demographic analysis
where mT is the number of the nucleotides of the entire DNA sequence, μ is the substitution rate per nucleotide site per year, and g is the generation time in years.
We would like to thank Professor Ji Yang at Fudan University, Shanghai, for his kind help in phylogenetic analyses and Li Zhang for her excellent skill in DNA sequencing. We also wish to thank Dr. Ihsan Al-Shehbaz at Missouri Botanical Garden, USA, and Dr. Steve O'Kane Jr at University of Northern Iowa, USA for providing the seeds of some out-group species. This study was supported by the National Basic Research Program of China (Grant No. 2006CB100105) and National Natural Science Foundation of China (30370093).
- Koornneef M, Alonso-Blanco C, Vreugdenhil D: Naturally occurring genetic variation in Arabidopsis thaliana. Annu Rev Plant Biol. 2004, 55: 141-172. 10.1146/annurev.arplant.55.031903.141605.PubMedView ArticleGoogle Scholar
- Bergelson J, Stahl E, Dudek S, Kreitman M: Genetic variation within and among populations of Arabidopsis thaliana. Genetics. 1998, 148: 1311-1323.PubMedPubMed CentralGoogle Scholar
- Loridon K, Cournoyer B, Goubely C, Depeiges A, Picard G: Length polymorphism and allele structure of trinucleotide microsatellites in natural accessions of Arabidopsis thaliana. Theor Appl Genet. 1998, 97: 591-604. 10.1007/s001220050935.View ArticleGoogle Scholar
- Zwan Vander C, Brodie SA, Campanella JJ: The intraspecific phylogenetics of Arabidopsis thaliana in worldwide populations. Syst Bot. 2000, 25: 47-59. 10.2307/2666672.View ArticleGoogle Scholar
- Sharbel TF, Haubold B, Mitchell-Olds T: Genetic isolation by distance in Arabidopsis thaliana: biogeography and postglacial colonization of Europe. Mol Ecol. 2000, 9: 2109-2118. 10.1046/j.1365-294X.2000.01122.x.PubMedView ArticleGoogle Scholar
- Barth S, Melchinger AE, Lubberstedt T: Genetic diversity in Arabidopsis thaliana (L.) Heynh investigated by cleaved amplified polymorphic sequence (CAPS) and inter-simple sequence repeat (ISSR) markers. Mol Ecol. 2002, 11: 495-505. 10.1046/j.0962-1083.2002.01466.x.PubMedView ArticleGoogle Scholar
- He F, Kang D, Ren Y, Qu LJ, Zhen Y, Gu H: Genetic diversity of the natural populations of Arabidopsis thaliana in China. Heredity. 2007, 99: 423-431. 10.1038/sj.hdy.6801020.PubMedView ArticleGoogle Scholar
- Beck JB, Schmuths H, Schaal BA: Native range genetic variation in Arabidopsis thaliana is strongly geographically structured and reflects Pleistocene glacial dynamics. Mol Ecol. 2008, 17: 902-915. 10.1111/j.1365-294X.2008.03906.x.PubMedView ArticleGoogle Scholar
- Pico FX, Mendez-Vigo B, Martinez-Zapater JM, Alonso-Blanco C: Natural genetic variation of Arabidopsis thaliana is geographically structured in the Iberian Peninsula. Genetics. 2008, 180: 1009-1021. 10.1534/genetics.108.089581.PubMedPubMed CentralView ArticleGoogle Scholar
- Hanfstingl U, Berry A, Kellogg EA, Costa JT, Rudiger W, Ausubel FM: Haplotypic divergence coupled with lack of diversity at the Arabidopsis thaliana alcohol-dehydrogenase locus - roles for both balancing and directional selection. Genetics. 1994, 138: 811-828.PubMedPubMed CentralGoogle Scholar
- Kawabe A, Innan H, Terauchi R, Miyashita NT: Nucleotide polymorphism in the acidic chitinase focus (ChiA) region of the wild plant Arabidopsis thaliana. Mol Biol Evol. 1997, 14: 1303-1315.PubMedView ArticleGoogle Scholar
- Miyashita NT, Kawabe A, Innan H: DNA variation in the wild plant Arabidopsis thaliana revealed by amplified fragment length polymorphism analysis. Genetics. 1999, 152: 1723-1731.PubMedPubMed CentralGoogle Scholar
- Purugganan MD, Suddith JI: Molecular population genetics of floral homeotic loci: Departures from the equilibrium-neutral model at the APETALA3 and PISTILLATA genes of Arabidopsis thaliana. Genetics. 1999, 151: 839-848.PubMedPubMed CentralGoogle Scholar
- Kawabe A, Yamane K, Miyashita NT: DNA polymorphism at the cytosolic phosphoglucose isomerase (PgiC) locus of the wild plant Arabidopsis thaliana. Genetics. 2000, 156: 1339-1347.PubMedPubMed CentralGoogle Scholar
- Aguade M: Nucleotide sequence variation at two genes of the phenylpropanoid pathway, the FAH1 and F3H genes, in Arabidopsis thaliana. Mol Biol Evol. 2001, 18: 1-9.PubMedView ArticleGoogle Scholar
- Hauser MT, Harr B, Schlotterer C: Trichome distribution in Arabidopsis thaliana and its close relative Arabidopsis lyrata: Molecular analysis of the candidate gene GLABROUS1. Mol Biol Evol. 2001, 18: 1754-1763.PubMedView ArticleGoogle Scholar
- Tian DC, Araki H, Stahl E, Bergelson J, Kreitman M: Signature of balancing selection in Arabidopsis. Proc Natl Acad Sci USA. 2002, 99: 11525-11530. 10.1073/pnas.172203599.PubMedPubMed CentralView ArticleGoogle Scholar
- Mauricio R, Stahl EA, Korves T, Tian DC, Kreitman M, Bergelson J: Natural selection for polymorphism in the disease resistance gene Rps2 of Arabidopsis thaliana. Genetics. 2003, 163: 735-746.PubMedPubMed CentralGoogle Scholar
- Shepard KA, Purugganan MD: Molecular population genetics of the Arabidopsis CLAVATA2 region: The genomic scale of variation and selection in a selfing species. Genetics. 2003, 163: 1083-1095.PubMedPubMed CentralGoogle Scholar
- Yoshida K, Kamiya T, Kawabe A, Miyashita NT: DNA polymorphism at the ACAULIS5 locus of the wild plant Arabidopsis thaliana. Genes Genet Syst. 2003, 78: 11-21. 10.1266/ggs.78.11.PubMedView ArticleGoogle Scholar
- Abbott RJ, Gomes MF: Population genetic-structure and outcrossing rate of Arabidopsis thaliana (L.) Heynh. Heredity. 1989, 62: 411-418. 10.1038/hdy.1989.56.View ArticleGoogle Scholar
- Innan H, Terauchi R, Miyashita NT: Microsatellite polymorphism in natural populations of the wild plant Arabidopsis thaliana. Genetics. 1997, 146: 1441-1452.PubMedPubMed CentralGoogle Scholar
- Ding J, Zhang W, Jing Z, Chen J-Q, Tian D: Unique pattern of R - gene variation within populations in Arabidopsis. Mol Genet Genomics. 2007, 277: 619-629. 10.1007/s00438-007-0213-5.PubMedView ArticleGoogle Scholar
- Sato S, Nakamura Y, Kaneko T, Asamizu E, Tabata S: Complete structure of the chloroplast genome of Arabidopsis thaliana. DNA Res. 1999, 6: 283-290. 10.1093/dnares/6.5.283.PubMedView ArticleGoogle Scholar
- Sall T, Jakobsson M, Lind-Hallden C, Hallden C: Chloroplast DNA indicates a single origin of the allotetraplold Arabidopsis suecica. J Evol Biol. 2003, 16: 1019-1029. 10.1046/j.1420-9101.2003.00554.x.PubMedView ArticleGoogle Scholar
- Doyle JJ, Davis JI, Soreng RJ, Garvin D, Anderson MJ: Chloroplast DNA inversions and the origin of the grass family (Poaceae). Proc Natl Acad Sci USA. 1992, 89: 7722-7726. 10.1073/pnas.89.16.7722.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim KJ, Choi KS, Jansen RKL: Two chloroplast DNA inversions originated simultaneously during the early evolution of the sunflower family (Asteraceae). Mol Biol Evol. 2005, 22: 1783-1792. 10.1093/molbev/msi174.PubMedView ArticleGoogle Scholar
- Castric V, Bechsgaard J, Schierup MH, Vekemans X: Repeated adaptive introgression at a gene under multiallelic balancing selection. PLoS Genetics. 2008, 4: e1000168.PubMedPubMed CentralView ArticleGoogle Scholar
- Schmid KJ, Ramos-Onsins S, Ringys-Beckstein H, Weisshaar B, Mitchell-Olds T: A multilocus sequence survey in Arabidopsis thaliana reveals a genome-wide departure from a neutral model of DNA sequence polymorphism. Genetics. 2005, 169: 1601-1615. 10.1534/genetics.104.033795.PubMedPubMed CentralView ArticleGoogle Scholar
- Aris-Brosou S, Excoffier L: The impact of population expansion and mutation rate heterogeneity on DNA sequence polymorphism. Mol Biol Evol. 1996, 13: 494-504.PubMedView ArticleGoogle Scholar
- Rieseberg LH, Brunsfeld SJ: Molecular evidence and plant introgression. Molecular systematics of plants. Edited by: Soltis PS, Soltis DE, Doyle JJ. New York: Chapman and Hall, 1992:151-176.View ArticleGoogle Scholar
- Smissen RD, Breitwieser I, Ward JM: Phylogenetic implications of trans-specific chloroplast DNA sequence polymorphism in New Zealand Gnaphalieae (Asteraceae). Plant Syst Evol. 2004, 249: 37-53. 10.1007/s00606-004-0209-0.View ArticleGoogle Scholar
- Clauss MJ, Koch MA: Poorly known relatives of Arabidopsis thaliana. Trends Plant Sci. 2006, 11: 449-459. 10.1016/j.tplants.2006.07.005.PubMedView ArticleGoogle Scholar
- O'Kane SL, Schaal BA, Al-Shehbaz IA: The Origins of Arabidopsis suecica (Brassicaceae) as indicated by nuclear rDNA sequences. Systematic Botany. 1996, 21: 559-566. 10.2307/2419615.View ArticleGoogle Scholar
- Pontes O, Neves N, Silva M, Lewis MS, Madlung A, Comai L, Viegas W, Pikaard CS: Chromosomal locus rearrangements are a rapid response to formation of the allotetraploid Arabidopsis suecica genome. Proc Natl Acad Sci USA. 2004, 101: 18240-18245. 10.1073/pnas.0407258102.PubMedPubMed CentralView ArticleGoogle Scholar
- Ostrowski MF, David J, Santoni S, McKhann H, Reboud X, Le Corre V, Camilleri C, Brunel D, Bouchez D, Faure B, Bataillon T: Evidence for a large-scale population structure among accessions of Arabidopsis thaliana: possible causes and consequences for the distribution of linkage disequilibrium. Mol Ecol. 2006, 15: 1507-1517. 10.1111/j.1365-294X.2006.02865.x.PubMedView ArticleGoogle Scholar
- Berger B: The taxonomic confusion within Arabidopsis and allied genera. Arabidopsis research. Edited by: Robbelen G. University of Gottingen, Gottingen, Germany, 1965:19-25.Google Scholar
- Levey S, Wingler A: Natural variation in the regulation of leaf senescence and relation to other traits in Arabidopsis. Plant, Cell and Envir. 2005, 28: 223-231. 10.1111/j.1365-3040.2004.01266.x.View ArticleGoogle Scholar
- Doyle JJ, Doyle JL: A rapid total DNA preparation procedure for fresh plant tissue. Focus. 1990, 12: 13-15.Google Scholar
- Primer Premier. [http://www.premierbiosoft.com/primerdesign/index.html]
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.PubMedView ArticleGoogle Scholar
- GenBank. [http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/Genbank/]
- Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.PubMedView ArticleGoogle Scholar
- Nei M: Molecular evolutionary genetics. Columbia University Press, New York 1987.Google Scholar
- Tajima F: Statistical-method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123: 585-595.PubMedPubMed CentralGoogle Scholar
- Fu YX, Li WH: Statistical tests of neutrality of mutations. Genetics. 1993, 133: 693-709.PubMedPubMed CentralGoogle Scholar
- Tamura K, Nei M, Kumar S: Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci USA. 2004, 101: 11030-11035. 10.1073/pnas.0404206101.PubMedPubMed CentralView ArticleGoogle Scholar
- Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony (* and Other Methods). Version 4. Sinauer Associates, Sunderland, Mass 1998.Google Scholar
- Bandelt HJ, Forster P, Rohl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48.PubMedView 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.PubMedView ArticleGoogle Scholar
- Koch MA, Matschinger M: Evolution and genetic differentiation among relatives of Arabidopsis thaliana. Proc Natl Acad Sci USA. 2007, 104: 6272-6277. 10.1073/pnas.0701338104.PubMedPubMed CentralView ArticleGoogle Scholar
- Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide-sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.PubMedView ArticleGoogle Scholar
- Slatkin M, Hudson RR: Pairwise comparisons of mitochondrial-DNA sequences in stable and exponentially growing populations. Genetics. 1991, 129: 555-562.PubMedPubMed CentralGoogle Scholar
- Rogers AR, Harpending H: Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992, 9: 552-569.PubMedGoogle Scholar
- Rogers AR: Genetic evidence for a Pleistocene population explosion. Evolution. 1995, 49: 608-615. 10.2307/2410314.View ArticleGoogle Scholar
- Excoffier L, Laval G, Schneider S: ARLEQUIN (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform Online. 2005, 1: 47-50.PubMed CentralGoogle Scholar
- Harpending HC: Signature of ancient population-growth in a low-resolution mitochondrial-DNA mismatch distribution. Hum Biol. 1994, 66: 591-600.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.