Evidence for a large-scale population structure among accessions of Arabidopsis thaliana: possible causes and consequences for the distribution of linkage disequilibrium

The existence of a large-scale population structure was investigated in Arabidopsis thaliana by studying patterns of polymorphism in a set of 71 European accessions. We used sequence polymorphism surveyed in 10 fragments of approximately 600 nucleotides and a set of nine microsatellite markers. Population structure was investigated using a model-based inference framework. Among the accessions studied, the presence of four groups was inferred using genetic data, without using prior information on the geographical origin of the accessions. Significant genetic isolation by geographical distance was detected at the group level, together with a geographical gradient in allelic richness across groups. These results are discussed with respect to the previously proposed scenario of postglacial colonization of Europe from putative glacial refugia. Finally, the contribution of the inferred structure to linkage disequilibrium among 171 pairs of essentially unlinked markers was also investigated. Linkage disequilibrium analysis revealed that significant associations detected in the whole sample were mainly due to genetic differentiation among the inferred groups. We discuss the implication of this finding for future association studies in A. thaliana.     

Evidence for a large-scale population structure of Arabidopsis thaliana from genome-wide single nucleotide polymorphism markers

Population-based methods for the genetic mapping of adaptive traits and the analysis of natural selection require that the population structure and demographic history of a species are taken into account. We characterized geographic patterns of genetic variation in the model plant Arabidopsis thaliana by genotyping 115 genome-wide single nucleotide polymorphism (SNP) markers in 351 accessions from the whole species range using a matrix-assisted laser desorption/ionization time-of-flight assay, and by sequencing of nine unlinked short genomic regions in a subset of 64 accessions. The observed frequency distribution of SNPs is not consistent with a constant-size neutral model of sequence polymorphism due to an excess of rare polymorphisms. There is evidence for a significant population structure as indicated by differences in genetic diversity between geographic regions. Accessions from Central Asia have a low level of polymorphism and an increased level of genome-wide linkage disequilibrium (LD) relative to accessions from the Iberian Peninsula and Central Europe. Cluster analysis with the structure program grouped Eurasian accessions into K=6 clusters. Accessions from the Iberian Peninsula and from Central Asia constitute distinct populations, whereas Central and Eastern European accessions represent admixed populations in which genomes were reshuffled by historical recombination events. These patterns likely result from a rapid postglacial recolonization of Eurasia from glacial refugial populations. Our analyses suggest that mapping populations for association or LD mapping should be chosen from regional rather than a species-wide sample or identified genetically as sets of individuals with similar average genetic distances. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Neutral genetic variation among wild North American populations of the weedy plant Arabidopsis thaliana is not geographically structured

We investigated neutral genetic variation within and among 53 wild-collected populations of the weedy annual plant, Arabidopsis thaliana, in North America, using amplified fragment length polymorphism (AFLP) markers. A. thaliana is thought to have been introduced to North America from Eurasia by humans; such an introduction might be expected to leave a clear geographical signal in the genetic data. To detect such patterns, we sampled populations at several hierarchical geographical levels. We collected individuals from populations in two areas of the Southeast and one in the Midwest, as well as individuals from populations in the Pacific Northwest and Northeast. To estimate within-population variation, we sampled eight individuals from each of six populations in the Southeast and Midwest. Among all 95 individuals analysed, we detected 131 polymorphic AFLP fragments. We found no evidence for continental or regional diversification. Individuals sampled from Midwestern and Southeastern populations intermingled in a neighbour-joining tree, and Mantel tests conducted within the Midwestern and Southeastern regions as well as the full data set failed to detect any significant relationship between geographical and genetic distance. These results mirror those found for most global surveys of neutral genetic variability in A. thaliana. Surprisingly, we detected substantial amounts of neutral genetic variability within populations. The levels of genetic variation within populations, coupled with the nongeographical nature of divergence among populations, are consistent with contemporary gene flow and point to a complex and dynamic population history of A. thaliana in North America.     

Genetic diversity in Arabidopsis thaliana L. Heynh. investigated by cleaved amplified polymorphic sequence (CAPS) and inter-simple sequence repeat (ISSR) markers

In this study, we investigated genetic diversity among 37 accessions in Arabidopsis thaliana from Eurasia, North Africa and North America using morphological traits and two polymerase chain reaction (PCR)-based marker systems: cleaved amplified polymorphic sequences (CAPS) and inter-simple sequence repeats (ISSR). Cluster analysis based on genetic similarities calculated from CAPS data grouped the accessions roughly according to their geographical origin: one large group contained accessions from Western, Northern and Southern Europe as well as North Africa, a second group consisted of Eastern European and Asian continental accessions. North American accessions were interspersed into these groups. Contrary to the CAPS analysis, the dendrogram obtained from the ISSR data did not reflect the geographical origin of the accessions, and the calculated genetic distances did not match the CAPS results. This could be attributable to an uneven genomic distribution of ISSR markers as substantiated by a database search for ISSR binding sites in A. thaliana genomic DNA sequence files, or to the ISSR's different mode of evolution. We recommend CAPS markers for diversity analysis in A. thaliana because a careful selection of markers can ascertain an even representation of the entire genome.     

Clinal variation in freezing tolerance among natural accessions of Arabidopsis thaliana

Low temperature represents a form of abiotic stress that varies predictably with latitude and altitude and to which organisms have evolved multiple physiological responses. Plants provide an especially useful experimental system for investigating the ecological and evolutionary dynamics of tolerance to low temperature because of their sessile lifestyle and inability to escape ambient atmospheric conditions. Here, intraspecific variation in freezing tolerance was investigated in Arabidopsis thaliana by conducting freezing tolerance assays on 71 accessions collected from across the native range of the species. Assays were performed at multiple minimum temperatures and on both cold-acclimated and non-cold-acclimated individuals. Considerable variation in freezing tolerance was observed among accessions both with and without a prior cold-acclimation treatment, suggesting that differences among accessions in cold-acclimation capacity as well as differences in intrinsic physiology contribute to variation in this phenotype. A highly significant positive relationship was observed between freezing tolerance and latitude of origin of accessions, consistent with a major role for natural selection in shaping variation in this phenotype. Clinal variation in freezing tolerance in A. thaliana coupled with considerable knowledge of the underlying genetics and physiology of this phenotype should allow evolutionary genetic analysis at multiple levels.     

Genetic variation in Arabidopsis suecica and its parental species A. arenosa and A. thaliana

Random amplified polymorphic DNA (RAPD) markers were used to estimate the level of genetic variation in Swedish accessions of the allopolyploid Arabidopsis suecica and its parental species A. thaliana and A. arenosa. The results showed clear differences among the three species with respect to the level of variation. A. arenosa was highly variable, A. thaliana showed a moderate level of variation whereas A. suecica was much less variable than the two other species. An extended analysis covering 19 Swedish populations of A. suecica corroborated the low level of variation in this species, yet 16 unique phenotypes were observed. No isolation by distance was observed. When the genetic variation was partitioned among and within populations of A. suecica, the results showed that the majority of the variation (81%) occurred among populations. This result is interpreted as a strong indication that A. suecica is autogamous in nature.     

The genetic architecture of shoot branching in Arabidopsis thaliana: a comparative assessment of candidate gene associations vs. quantitative trait locus mapping

Association mapping focused on 36 genes involved in branch development was used to identify candidate genes for variation in shoot branching in Arabidopsis thaliana. The associations between four branching traits and moderate-frequency haplogroups at the studied genes were tested in a panel of 96 accessions from a restricted geographic range in Central Europe. Using a mixed-model association-mapping method, we identified three loci--MORE AXILLARY GROWTH 2 (MAX2), MORE AXILLARY GROWTH 3 (MAX3), and SUPERSHOOT 1 (SPS1)--that were significantly associated with branching variation. On the basis of a more extensive examination of the MAX2 and MAX3 genomic regions, we find that linkage disequilibrium in these regions decays within approximately 10 kb and trait associations localize to the candidate genes in these regions. When the significant associations are compared to relevant quantitative trait loci (QTL) from previous Ler x Col and Cvi x Ler recombinant inbred line (RIL) mapping studies, no additive QTL overlapping these candidate genes are observed, although epistatic QTL for branching, including one that spans the SPS1, are found. These results suggest that epistasis is prevalent in determining branching variation in A. thaliana and may need to be considered in linkage disequilibrium mapping studies of genetically diverse accessions.     

Population genetic structure of Plasmopara viticola after 125 years of colonization in European vineyards

To examine the within- and among-population genetic structure of Plamopara viticola oosporic populations in Europe, 8991 lesions from 32 vineyard plots were collected and analysed. Four multi-allelic microsatellite markers were used to genotype the pathogen. All populations had high levels of gene and genotypic diversity. Most populations were in Hardy–Weinberg equilibrium and thus randomly mating. Among P. viticola populations, significant low to moderate genetic differentiation was observed, even between geographically close populations. This genetic differentiation was also evident in the neighbour-joining phylogenetic genetic distance tree, showing clear substructure and distinguishing mainly five clusters based on geographical origin. Significant isolation by distance was found in central European P. viticola populations, suggesting a step-wise migration model. No significant isolation by distance was found within Greek populations, most probably owing to natural geographical barriers such as the sea and mountains, as well as the frequent population bottlenecks occurring in these populations, preventing natural migration among populations. The high variability of P. viticola provides explanation for its successful infestation of the heterogeneous European vineyards in the last 125 years after its introduction.     

Unique evolutionary mechanism in R-genes under the presence/absence polymorphism in Arabidopsis thaliana

While the presence/absence polymorphism is commonly observed in disease resistance (R-) genes in Arabidopsis, only a few R-genes under the presence/absence polymorphism (R-P/A) have been investigated. To understand the mechanism of the molecular evolution of R-P/A, we investigated genetic variation of nine R-P/A in A. thaliana from worldwide populations. The number of possessed R-genes varied widely among accessions (two to nine, on average 4.3 +/- 1.6/accession). No pair of accessions shared the same haplotype, and no clear geographic differentiation was observed with respect to the pattern of presence/absence of the R-genes investigated. Presence allele frequencies also varied among loci (25-70%), and no linkage disequilibrium was detected among them. Although the LRR region in regular R-genes is known to be highly polymorphic and has a high Ka/Ks ratio in A. thaliana, nucleotide sequences of this region in the R-P/A showed a relatively low level of genetic variation (pi = 0.0002-0.016) and low Ka/Ks (0.03-0.70, <1). In contrast, the nucleotide diversities around the deletion junction of R-P/A were constantly high between presence and absence accessions for the R-genes (D(xy) = 0.031-0.103). Our results suggest that R-P/A loci evolved differently from other R-gene loci and that balancing selection plays an important role in molecular evolution of R-P/A.     

Looking for natural variation in chiasma frequency in Arabidopsis thaliana

Information concerning natural variation either in chiasma frequency or in the genetic basis of any such variation is a valuable tool to characterize phenotypic traits and their genetic control. Here meiotic recombination frequencies are analysed in nine geographically and ecologically diverse accessions of Arabidopsis thaliana, and a comparative study was carried out incorporating previous data from another eight accessions. Chiasma frequencies, estimated by counting rod and ring bivalents at metaphase I, varied up to 22% among accessions. However, no differences were found among plants of the same accession. There was a relationship, which does not necessarily imply direct proportionality, between the size of the chromosomes and their mean chiasma frequency. Chiasma frequency and distribution between arms and among chromosomes were not consistent over accessions. These findings indicate the existence of genetic factors controlling meiotic recombination both throughout the whole genome and at the whole chromosome level. The reliability of chiasma scoring as an indicator of reciprocal recombination events is also discussed.     

Mapping salinity tolerance during Arabidopsis thaliana germination and seedling growth

To characterize and dissect genetic variation for salinity tolerance, we assessed variation in salinity tolerance during germination and seedling growth for a worldwide sample of Arabidopsis thaliana accessions. By combining QTL mapping, association mapping and expression data, we identified genomic regions involved in salinity response. Among the worldwide sample, we found germination ability within a moderately saline environment (150 mM NaCl) varied considerable, from >90% among the most tolerant lines to complete inability to germinate among the most susceptible. Our results also demonstrated wide variation in salinity tolerance within A. thaliana RIL populations and identified multiple genomic regions that contribute to this variation. These regions contain known candidate genes, but at least four of the regions contain loci not yet associated with salinity tolerance response phenotypes. Our observations suggest A. thaliana natural variation may be an underutilized resource for investigating salinity stress response.     

Natural genetic resources of Arabidopsis thaliana reveal a high prevalence and unexpected phenotypic plasticity of RPW8-mediated powdery mildew resistance

Here, an approach based on natural genetic variation was adopted to analyse powdery mildew resistance in Arabidopsis thaliana. Accessions resistant to multiple powdery mildew species were crossed with the susceptible Col-0 ecotype and inheritance of resistance was analysed. Histochemical staining was used to visualize archetypal plant defence responses such as callose deposition, hydrogen peroxide accumulation and host cell death in a subset of these ecotypes. In six accessions, resistance was likely of polygenic origin while 10 accessions exhibited evidence for a single recessively or semi-dominantly inherited resistance locus. Resistance in the latter accessions was mainly manifested at the terminal stage of the fungal life cycle by a failure of abundant conidiophore production. The resistance locus of several of these ecotypes was mapped to a genomic region containing the previously analysed atypical RPW8 powdery mildew resistance genes. Gene silencing revealed that members of the RPW8 locus were responsible for resistance to Golovinomyces orontii in seven accessions. These results suggest that broad-spectrum powdery mildew resistance in A. thaliana is predominantly of polygenic origin or based on RPW8 function. The findings shed new light on the natural variation of inheritance, phenotypic expression and pathogen range of RPW8-conditioned powdery mildew resistance.     

Polyploid speciation did not confer instant reproductive isolation in Capsella (Brassicaceae)

Polyploid formation is a major mode of sympatric speciation in flowering plants. Unlike other speciation processes, polyploidization is often assumed to confer instant reproductive isolation. Shared polymorphism across ploidy levels has therefore often been attributed to multiple polyploid origins, whereas the alternative hypothesis of introgressive hybridization has rarely been rigorously tested. Here, we sequence 12 nuclear loci representing 6 genes duplicated by polyploidy in 92 accessions of the tetraploid Capsella bursa-pastoris together with the corresponding loci in 21 accessions of its close diploid relative Capsella rubella. In C. bursa-pastoris accessions from western Eurasia, where the 2 species occur in partial sympatry, we find higher levels of nucleotide diversity than in accessions from eastern Eurasia, where C. rubella does not grow. Furthermore, haplotypes are shared across ploidy levels at 4 loci in western but not in eastern Eurasia. We test whether haplotype sharing is due to retention of ancestral polymorphism or due to hybridization and introgression using a coalescent-based isolation-with-migration model. In western but not in eastern Eurasia, there is evidence for unidirectional gene flow from C. rubella to C. bursa-pastoris. An independent estimate of the timing of dispersal of C. bursa-pastoris to eastern Eurasia indicates that it probably predated introgression. Our results show that polyploid speciation need not result in immediate and complete reproductive isolation, that postpolyploidization hybridization and introgression can contribute significantly to genetic variation in a newly formed polyploid, and that divergence population genetic analysis constitutes a powerful way of testing hypotheses on polyploid speciation.     

Natural variation of submergence tolerance among Arabidopsis thaliana accessions

? The exploitation of natural variation in Arabidopsis thaliana (Arabidopsis) provides a huge potential for the identification of the molecular mechanisms underlying this variation as a result of the availability of a vast array of genetic and genomic resources for this species. Eighty-six Arabidopsis accessions were screened for natural variation in flooding tolerance. This forms the first step towards the identification and characterization of the role of candidate genes contributing to flooding tolerance. ? Arabidopsis accessions at the 10-leaf stage were subjected to complete submergence in the dark. Survival curves were plotted to estimate median lethal times as a measure of tolerance. Flooding-associated survival parameters, such as root and shoot oxygen content, initial carbohydrate content and petiole elongation under water, were also measured. ? There was a significant variation in submergence tolerance among Arabidopsis accessions. However, the order of tolerance did not correlate with root and shoot oxygen content or initial amounts of shoot starch and total soluble sugars. A negative correlation was observed between submergence tolerance and underwater petiole elongation. ? Arabidopsis accessions show considerable variation in the ability to tolerate complete submergence, making it a good species in which to identify and characterize genes and to study mechanisms that contribute to survival under water.     

The genetic structure of Drosophila ananassae populations from Asia, Australia and Samoa

Information about genetic structure and historical demography of natural populations is central to understanding how natural selection changes genomes. Drosophila ananassae is a widespread species occurring in geographically isolated or partially isolated populations and provides a unique opportunity to investigate population structure and molecular variation. We assayed microsatellite repeat-length variation among 13 populations of D. ananassae to assess the level of structure among the populations and to make inferences about their ancestry and historic biogeography. High levels of genetic structure are apparent among all populations, particularly in Australasia and the South Pacific, and patterns are consistent with the hypothesis that the ancestral populations are from Southeast Asia. Analysis of population structure and use of F-statistics and Bayesian analysis suggest that the range expansion of the species into the Pacific is complex, with multiple colonization events evident in some populations represented by lineages that show no evidence of recent admixture. The demographic patterns show isolation by distance among populations and population expansion within all populations. A morphologically distinct sister species, D. pallidosa, collected in Malololelei, Samoa, appears to be more closely related to some of the D. ananassae populations than many of the D. ananassae populations are to one another. The patterns of genotypic diversity suggest that many of the individuals that we sampled may be morphologically indistinguishable nascent species.     

Native range genetic variation in Arabidopsis thaliana is strongly geographically structured and reflects Pleistocene glacial dynamics

Despite Arabidopsis thaliana 's pre-eminence as a model organism, major questions remain regarding the geographic structure of its genetic variation due to the geographically incomplete sample set available for previous studies. Many of these questions are addressed here with an analysis of genome-wide variation at 10 loci in 475 individuals from 167 globally distributed populations, including many from critical but previously un-sampled regions. Rooted haplotype networks at three loci suggest that A. thaliana arose in the Caucasus region. Identification of large-scale metapopulations indicates clear east–west genetic structure, both within proposed Pleistocene refugia and post-Pleistocene colonized regions. The refugia themselves are genetically differentiated from one another and display elevated levels of within-population genetic diversity relative to recolonized areas. The timing of an inferred demographic expansion coincides with the Eemian interglacial (approximately 120 000 years ago). Taken together, these patterns are strongly suggestive of Pleistocene range dynamics. Spatial autocorrelation analyses indicate that isolation by distance is pervasive at all hierarchical levels, but that it is reduced in portions of Europe.     

Population genetic structure and hybridization patterns in the Mediterranean endemics Phlomis lychnitis and P. crinita (Lamiaceae)

BACKGROUND AND AIMS: The historical influence of gene flow and genetic drift after the last glacial phase of the Quaternary Period is reflected in current levels of genetic diversity and population structure of plant species. Moreover, hybridization after secondary contact might also affect population genetic diversity and structure. An assessment was made of the genetic variation and hybrid zone structure in Iberian populations of the Mediterranean Phlomis lychnitis and P. crinita, for which phylogenetic relationships are controversial, and hybridization and introgression are common. METHODS: Allozyme variation at 13 loci was analysed in 1723 individual plants sampled from 35 natural locations of P. lychnitis, P. crinita subsp. malacitana and P. crinita subsp. crinita in southern and eastern Spain. Standard genetic diversity parameters were calculated and patterns of genetic structure in each taxon were tested to fit the equilibrium between gene flow and genetic drift. Individual multilocus genotypes were subjected to Bayesian clustering analysis to estimate hybridization and introgression rates for both geographic regions. KEY RESULTS: Contrasting patterns in the distribution of genetic variation among the three taxa were found. Phlomis lychnitis showed no significant inbreeding, low genetic differentiation among populations and no evidence of isolation by distance. Phlomis crinita subsp. malacitana and P. crinita subsp. crinita showed high levels of genetic structure consistent with a pattern of gene flow-drift equilibrium. Higher instances of hybridization and introgression were detected in locations from southern Spain compared with locations from eastern Spain, matching unimodal and bimodal hybrid zones, respectively. CONCLUSIONS: High instances of historical gene flow, range expansion and altitudinal movement during the Quaternary Period, and lineage sorting can explain the diversity of patterns observed. The results suggest that P. lychnitis is the most differentiated lineage in the group; however, the relationship between the three taxa remains unclear.