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

Cryptic population genetic structure: the number of inferred clusters depends on sample size

Clustering methods have been used extensively to unravel cryptic population genetic structure. We investigated the effect of the number of individuals sampled in each location on the resulting number of clusters. Our study was motivated by recent results in Arabidopsis thaliana: studies in which more than one individual was sampled per location apparently have led to a much higher number of clusters than studies where only one individual was sampled in each location, as is generally done in this species. We show, using computer simulations and microsatellite data in A. thaliana, that the number of sampled individuals indeed has a strong impact on the number of resulting clusters. This effect is smaller if the sampled populations have a hierarchical structure. In most cases, sampling 5-10 individuals per population should be enough. The results argue for abandoning the concept of 'accessions' in partially selfing organisms.

Candidate gene association mapping of Arabidopsis flowering time

The pathways responsible for flowering time in Arabidopsis thaliana comprise one of the best characterized genetic networks in plants. We harness this extensive molecular genetic knowledge to identify potential flowering time quantitative trait genes (QTGs) through candidate gene association mapping using 51 flowering time loci. We genotyped common single nucleotide polymorphisms (SNPs) at these genes in 275 A. thaliana accessions that were also phenotyped for flowering time and rosette leaf number in long and short days. Using structured association techniques, we find that haplotype-tagging SNPs in 27 flowering time genes show significant associations in various trait/environment combinations. After correction for multiple testing, between 2 and 10 genes remain significantly associated with flowering time, with CO arguably possessing the most promising associations. We also genotyped a subset of these flowering time gene SNPs in an independent recombinant inbred line population derived from the intercrossing of 19 accessions. Approximately one-third of significant polymorphisms that were associated with flowering time in the accessions and genotyped in the outbred population were replicated in both mapping populations, including SNPs at the CO, FLC, VIN3, PHYD, and GA1 loci, and coding region deletions at the FRI gene. We conservatively estimate that approximately 4-14% of known flowering time genes may harbor common alleles that contribute to natural variation in this life history trait.

The 1001 genomes project for Arabidopsis thaliana

The 1001 Genomes project for Arabidopsis thaliana could provide an enormous boost for plant research for a modest financial investment.

We advocate here a 1001 Genomes project for Arabidopsis thaliana, the workhorse of plant genetics, which will provide an enormous boost for plant research with a modest financial investment.

Functional and geographical differentiation of candidate balanced polymorphisms in Arabidopsis thaliana

Molecular population genetic analysis of three chromosomal regions in Arabidopsis thaliana suggested that balancing selection might operate to maintain variation at three novel candidate adaptive trait genes, including SOLUBLE STARCH SYNTHASE I (SSI), PLASTID TRANSCRIPTIONALLY ACTIVE 7(PTAC7), and BELL-LIKE HOMEODOMAIN 10 (BLH10). If balanced polymorphisms are indeed maintained at these loci, then we would expect to observe functional variation underlying the previously detected signatures of selection. We observe multiple replacement polymorphisms within and in the 32 amino acids just upstream of the protein-protein interacting BELL domain at the BLH10 locus. While no clear protein sequence differences are found between allele types in SSI and PTAC7, these two genes show evidence for allele-specific variation in expression levels. Geographical patterns of allelic differentiation seem consistent with population stratification in this species and a significant longitudinal cline was observed at all three candidate loci. These data support a hypothesis of balancing selection at all three candidate loci and provide a basis for more detailed functional work by identifying possible functional differences that might be selectively maintained.

Genome sequence comparison of Col and Ler lines reveals the dynamic nature of Arabidopsis chromosomes

Large differences in plant genome sizes are mainly due to numerous events of insertions or deletions (indels). The balance between these events determines the evolutionary direction of genome changes. To address the question of what phenomena trigger these alterations, we compared the genomic sequences of two Arabidopsis thaliana lines, Columbia (Col) and Landsberg erecta (Ler). Based on the resulting alignments large indels (>100 bp) within these two genomes were analysed. There are ∼8500 large indels accounting for the differences between the two genomes. The genetic basis of their origin was distinguished as three main categories: unequal recombination (Urec)-derived, illegitimate recombination (Illrec)-derived and transposable elements (TE)-derived. A detailed study of their distribution and size variation along chromosomes, together with a correlation analyses, allowed us to demonstrate the impact of particular recombination-based mechanisms on the plant genome evolution. The results show that unequal recombination is not efficient in the removal of TEs within the pericentromeric regions. Moreover, we discovered an unexpectedly high influence of large indels on gene evolution pointing out significant differences between the various gene families. For the first time, we present convincing evidence that somatic events do play an important role in plant genome evolution.

Complex rearrangements lead to novel chimeric gene fusion polymorphisms at the Arabidopsis thaliana MAF2-5 flowering time gene cluster

Tandem gene clusters of multigene families are rearrangement hotspots and may be a major source of novel gene formation. Here, we report on a molecular population genetic analysis of the MAF2-5 gene cluster of the model plant species, Arabidopsis thaliana. The MAF2-5 genes are a MADS-box multigene family cluster spanning approximately 24 kbp on chromosome 5. We find heterogeneous evolutionary dynamics among these genes, all of which are closely related to the floral repressor, FLC, and are believed to play a role in the control of flowering time in A. thaliana. Low levels of nonsynonymous single nucleotide polymorphism (SNP) observed for MAF4 and MAF5 suggest purifying selection and conservation of function. In contrast, high levels of nonsynonymous SNPs, insertion-deletion, and rearrangements are observed for MAF2 and MAF3, including novel gene fusions that persist as a moderate-frequency polymorphism in A. thaliana. These fused genes, involving MAF2 and portions of MAF3, are expressed, resulting in the production of chimeric, alternatively spliced transcripts of MAF2. Association studies support a correlation between the described MAF2-MAF3 gene rearrangements and flowering time variation in the species. The finding that complex rearrangements within gene clusters, such as those observed for MAF2, might play a role in the generation of ecologically important phenotypic variation, emphasize the need for emerging high throughput genotyping and sequencing techniques to correctly reconstruct gene chimeras and other complex polymorphisms.

Flowering time diversification and dispersal in central Eurasian wild wheat Aegilops tauschii Coss.: genealogical and ecological framework

Timing of flowering is a reproductive trait that has significant impact on fitness in plants. In contrast to recent advances in understanding the molecular basis of floral transition, few empirical studies have addressed questions concerning population processes of flowering time diversification within species. We analyzed chloroplast DNA genealogical structure of flowering time variation in central Eurasian wild wheat Aegilops tauschii Coss. using 200 accessions that represent the entire species range. Flowering time measured as days from germination to flowering varied from 144.0 to 190.0 days (average 161.3 days) among accessions in a common garden/greenhouse experiment. Subsequent genealogical and statistical analyses showed that (1) there exist significant longitudinal and latitudinal clines in flowering time at the species level, (2) the early-flowering phenotype evolved in two intraspecific lineages, (3) in Asia, winter temperature was an environmental factor that affected the longitudinal clinal pattern of flowering time variation, and (4) in Transcaucasus-Middle East, some latitudinal factors affected the geographic pattern of flowering time variation. On the basis of palaeoclimatic, biogeographic, and genetic evidence, the northern part of current species' range [which was within the temperate desert vegetation (TDV) zone at the Last Glacial Maximum] is hypothesized to have harbored species refugia. Postglacial southward dispersal from the TDV zone seems to have been driven by lineages that evolved short-flowering-time phenotypes through different genetic mechanisms in Transcaucasus-Middle East and Asia.

Relaxed selection on the CBF/DREB1 regulatory genes and reduced freezing tolerance in the southern range of Arabidopsis thaliana

Elucidating the molecular basis of adaptive phenotypic variation represents a central aim in evolutionary biology. Traits exhibiting patterns of clinal variation represent excellent models for studies of molecular adaptation, especially when variation in phenotype can be linked to organismal fitness in different environments. Natural accessions of the model plant species Arabidopsis thaliana exhibit clinal variation in freezing tolerance that follows a gradient of temperature variability across the species' native range (Zhen Y, Ungerer MC. 2008. Clinal variation in freezing tolerance among natural accessions of A. thaliana. New Phytol. 177:419-427). Here, we report that this pattern of variation is attributable, at least in part, to relaxed purifying selection on members of a small family of transcriptional activators (the CBF/DREB1s) in the species' southern range. These regulatory genes play a critical role in the ability of A. thaliana plants to undergo cold acclimation and thereby achieve maximum freezing tolerance. Relative to accessions from northern regions, accessions of A. thaliana from the southern part of their geographic range exhibit levels of nonsynonymous nucleotide polymorphism that are approximately 2.8-fold higher across this small gene subfamily. Relaxed selection on the CBF/DREB1s in southern accessions also has resulted in multiple mutations in regulatory regions resulting in abrogated expression of particular subfamily members in particular accessions. These coding-region and regulatory mutations compromise the ability of these genes to act as efficient transcriptional activators during the cold acclimation process, as determined by reductions in rates of induction and maximum levels of expression in the downstream genes they regulate. This study highlights the potential role of regulatory genes in underlying adaptive phenotypic variation in nature.

Discordant longitudinal clines in flowering time and phytochrome C in Arabidopsis thaliana

Using seasonal cues to time reproduction appropriately is crucial for many organisms. Plants in particular often use photoperiod to signal the time to transition to flowering. Because seasonality varies latitudinally, adaptation to local climate is expected to result in corresponding clines in photoperiod-related traits. By experimentally manipulating photoperiod cues and measuring the flowering responses and photoperiod plasticity of 138 Eurasian accessions of Arabidopsis thaliana, we detected strong longitudinal but not latitudinal clines in flowering responses. The presence of longitudinal clines suggests that critical photoperiod cues vary among populations occurring at similar latitudes. Haplotypes at PHYC, a locus hypothesized to play a role in adaptation to light cues, were also longitudinally differentiated. Controlling for neutral population structure revealed that PHYC haplotype influenced flowering time; however, the distribution of PHYC haplotypes occurred in the opposite direction to the phenotypic cline, suggesting that loci other than PHYC are responsible for the longitudinal pattern in photoperiod response. Our results provide previously missing empirical support for the importance of PHYC in mediating photoperiod sensitivity in natural populations of A. thaliana. However, they also suggest that other loci and epistatic interactions likely play a role in the determination of flowering time and that the environmental factors influencing photoperiod in plants vary longitudinally as well as latitudinally.

Natural genetic variation of Arabidopsis thaliana is geographically structured in the Iberian peninsula

To understand the demographic history of Arabidopsis thaliana within its native geographical range, we have studied its genetic structure in the Iberian Peninsula region. We have analyzed the amount and spatial distribution of A. thaliana genetic variation by genotyping 268 individuals sampled in 100 natural populations from the Iberian Peninsula. Analyses of 175 individuals from 7 of these populations, with 20 chloroplast and nuclear microsatellite loci and 109 common single nucleotide polymorphisms, show significant population differentiation and isolation by distance. In addition, analyses of one genotype from 100 populations detected significant isolation by distance over the entire Iberian Peninsula, as well as among six Iberian subregions. Analyses of these 100 genotypes with different model-based clustering algorithms inferred four genetic clusters, which show a clear-cut geographical differentiation pattern. On the other hand, clustering analysis of a worldwide sample showed a west-east Eurasian longitudinal spatial gradient of the commonest Iberian genetic cluster. These results indicate that A. thaliana genetic variation displays significant regional structure and consistently support the hypothesis that Iberia has been a glacial refugium for A. thaliana. Furthermore, the Iberian geographical structure indicates a complex regional population dynamics, suggesting that this region contained multiple Pleistocene refugia with a different contribution to the postglacial colonization of Europe.

Patterns of polymorphism and demographic history in natural populations of Arabidopsis lyrata

BACKGROUND

Many of the processes affecting genetic diversity act on local populations. However, studies of plant nucleotide diversity have largely ignored local sampling, making it difficult to infer the demographic history of populations and to assess the importance of local adaptation. Arabidopsis lyrata, a self-incompatible, perennial species with a circumpolar distribution, is an excellent model system in which to study the roles of demographic history and local adaptation in patterning genetic variation.

PRINCIPAL FINDINGS

We studied nucleotide diversity in six natural populations of Arabidopsis lyrata, using 77 loci sampled from 140 chromosomes. The six populations were highly differentiated, with a median FST of 0.52, and structure analysis revealed no evidence of admixed individuals. Average within-population diversity varied among populations, with the highest diversity found in a German population; this population harbors 3-fold higher levels of silent diversity than worldwide samples of A. thaliana. All A. lyrata populations also yielded positive values of Tajima's D. We estimated a demographic model for these populations, finding evidence of population divergence over the past 19,000 to 47,000 years involving non-equilibrium demographic events that reduced the effective size of most populations. Finally, we used the inferred demographic model to perform an initial test for local adaptation and identified several genes, including the flowering time gene FCA and a disease resistance locus, as candidates for local adaptation events.

CONCLUSIONS

Our results underscore the importance of population-specific, non-equilibrium demographic processes in patterning diversity within A. lyrata. Moreover, our extensive dataset provides an important resource for future molecular population genetic studies of local adaptation in A. lyrata.

Clinal distribution of a chromosomal rearrangement: a precursor to chromosomal speciation?

Geographically structured genetic variants provide an effective means to assess sources of natural selection and mechanisms of adaptation to local environments. Correlated selection pressures along environmental gradients favor subdivision of genomes through chromosomal rearrangement. This study examines populations of Drosophila americana to evaluate selection pressures affecting chromosomal forms distinguished by a centromeric fusion. Analyses of chromosomal polymorphism throughout the Mississippi River Valley in the central United States reveal the presence of a distinct latitudinal cline for the chromosomal rearrangement. The cline has a width of 623 km centered at 35.97 degrees N and displays a characteristic sigmoid shape consistent with a balance between selection and dispersal. Extreme low temperature during January, an indicator of winter severity, was identified as the environmental variable that most accurately predicts arrangement frequency. An intriguing relationship identified between the chromosomal cline and operational sex ratio indicates that these alternative arrangements of the X chromosome may influence sex-specific survival. A hypothesis for the cline is presented wherein variation associated with the alternative chromosome forms influences distinct overwintering strategies. The resulting subdivision within the genome embodies a transitory stage of a speciation process in which locally adapted gene complexes provide a foundation for species formation.

Heterosis in the freezing tolerance, and sugar and flavonoid contents of crosses between Arabidopsis thaliana accessions of widely varying freezing tolerance

Heterosis is defined as the increased vigour of hybrids in comparison to their parents. We investigated 24 F(1) hybrid lines of Arabidopsis thaliana generated by reciprocally crossing either C24 or Col with six other parental accessions (Can, Co, Cvi, Ler, Rsch, Te) that differ widely in their freezing tolerance. The crosses differed in the degree of heterosis for freezing tolerance, both in the non-acclimated state and after a 14 d cold acclimation period. Crosses with C24 showed more heterosis than crosses with Col, and heterosis was stronger in acclimated than in non-acclimated plants. Leaf content of soluble sugars and proline showed more deviation from mid-parent values in crosses involving C24 than in those involving Col, and deviations were larger in acclimated than in non-acclimated plants. There were significant correlations between the content of different sugars and leaf freezing tolerance, as well as between heterosis effects in freezing tolerance and sugar content. Flavonoid content and composition varied between accessions, and between non-acclimated and acclimated plants. In the crosses, large deviations from the mid-parent values in the contents of different flavonols occurred, and there were strikingly strong correlations between both flavonol content and freezing tolerance, and between heterosis effects in freezing tolerance and flavonol content.

Natural variation in the degree of autonomous endosperm formation reveals independence and constraints of embryo growth during seed development in Arabidopsis thaliana

Seed development in flowering plants is a paradigm for the coordination of different tissues during organ growth. It requires a tight interplay between the two typically sexually produced structures: the embryo, developing from the fertilized egg cell, and the endosperm, originating from a fertilized central cell, along with the surrounding maternal tissues. Little is known about the presumptive signal transduction pathways administering and coordinating these different tissues during seed growth and development. Recently, a new signal has been identified emanating from the fertilization of the egg cell that triggers central cell proliferation without prior fertilization. Here, we demonstrate that there exists a large natural genetic variation with respect to the outcome of this signaling process in the model plant Arabidopsis thaliana. By using a recombinant inbred line population between the two Arabidopsis accessions Bayreuth-0 and Shahdara, we have identified two genetic components that influence the development of unfertilized endosperm. Exploiting this natural variation, we could further dissect the interdependence of embryo and endosperm growth during early seed development. Our data show an unexpectedly large degree of independence in embryo growth, but also reveal the embryo's developmental restrictions with respect to endosperm size. This work provides a genetic framework for dissection of the interplay between embryo and endosperm during seed growth in plants.

Demographic history of european populations of Arabidopsis thaliana

The model plant species Arabidopsis thaliana is successful at colonizing land that has recently undergone human-mediated disturbance. To investigate the prehistoric spread of A. thaliana, we applied approximate Bayesian computation and explicit spatial modeling to 76 European accessions sequenced at 876 nuclear loci. We find evidence that a major migration wave occurred from east to west, affecting most of the sampled individuals. The longitudinal gradient appears to result from the plant having spread in Europe from the east approximately 10,000 years ago, with a rate of westward spread of approximately 0.9 km/year. This wave-of-advance model is consistent with a natural colonization from an eastern glacial refugium that overwhelmed ancient western lineages. However, the speed and time frame of the model also suggest that the migration of A. thaliana into Europe may have accompanied the spread of agriculture during the Neolithic transition.

Six new recombinant inbred populations for the study of quantitative traits in Arabidopsis thaliana

Quantitative approaches are now widely used to study the genetic architecture of complex traits. However, most studies have been conducted in single mapping populations, which sample only a fraction of the natural allelic variation available within a gene pool and can identify only a subset of the loci controlling the traits. To enable the progress towards an understanding of the global genetic architecture of a broad range of complex traits, we have developed and characterised six new Arabidopsis thaliana recombinant inbred populations. To evaluate the utility of these populations for integrating analyses from multiple populations, we identified quantitative trait loci (QTL) controlling flowering time in vernalized plants growing in 16 h days. We used the physical positions of markers to align the linkage maps of our populations with those of six existing populations. We identified seven QTL in genomic locations coinciding with those identified in previous studies and in addition a further eight QTL were identified.

Construction and analysis of 2 reciprocal Arabidopsis introgression line populations

Two new large reciprocal sets of introgression lines (ILs) were created between the Arabidopsis accessions Col-0 and C24. In both sets (78 ILs with Col-0 background and 62 ILs with C24 background), the donor segments cover almost the entire genome with an average substitution size of 18.3 cM. In addition to the basic sets of ILs, further subILs were developed for 2 genomic regions allowing better mapping resolution. SubILs carrying donor segments with candidate genes for flowering time and reduced fertility were used to demonstrate the usefulness of the reciprocal ILs for quantitative trait loci detection and fine mapping. For subIL development at high resolution around the reduced fertility locus, we used modified CelI-based assays in one-well format for both marker development and genotyping. This serves as a very flexible and cost-effective approach.

Multilocus analysis of variation using a large empirical data set: phenylpropanoid pathway genes in Arabidopsis thaliana

Detecting the signature of adaptation on nucleotide variation is often difficult in species that like Arabidopsis thaliana might have a complex demographic history. Recent re-sequencing surveys in this species provided genome-wide information that would mainly reflect its demographic history. We have used a large empirical data set (LED) as well as multilocus coalescent simulations to analyse sequence variation at loci involved in the phenylpropanoid pathway of this species. We surveyed and examined DNA sequence variation at nine of these loci (about 19.7 kb) in 23 accessions of A. thaliana and one accession of its closely related species Arabidopsis lyrata. Nucleotide variation was lower at nonsynonymous sites than at silent sites in all loci, indicating generalized functional constraint at the protein level. No association between variation and position in the metabolic pathway was detected. When the data were contrasted against the standard neutral model, significant deviations for silent variation were detected with Tajima's D, Fu's F(S) and Fay and Wu's H multilocus test statistics. These deviations were in the same direction than in previous large-scale multilocus analyses, suggesting a genome-wide effect. When the nine-locus data set was contrasted against the large empirical data set, the level (Watterson's theta) and pattern of variation (Tajima's D) detected in these loci did not deviate either at the single-locus or multilocus level from the corresponding empirical distributions. These results would support an important role of the demographic history of A. thaliana in shaping nucleotide variation at the nine studied phenylpropanoid loci. The potential and limitations of the empirical distribution approach are discussed.

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.

Whole-plant mineral partitioning throughout the life cycle in Arabidopsis thaliana ecotypes Columbia, Landsberg erecta, Cape Verde Islands, and the mutant line ysl1ysl3

Minimal information exists on whole-plant dynamics of mineral flow through Arabidopsis thaliana or on the source tissues responsible for mineral export to developing seeds. Understanding these phenomena in a model plant could help in the development of nutritionally enhanced crop cultivars. A whole-plant partitioning study, using sequential harvests, was conducted to characterize growth and mineral concentrations and contents of rosettes, cauline leaves, stems, immature fruit, mature fruit hulls, and seeds of three WT lines (Col-0, Ler, and Cvi) and one mutant line (Col-0::ysl1ysl3). Shoot mineral content increased throughout the life cycle for all minerals, although tissue-specific mineral partitioning differed between genotypes. In particular, iron (Fe), zinc (Zn), and copper (Cu) were aberrantly distributed in ysl1ysl3. Remobilization was observed for several minerals from various tissues, including cauline leaves and silique hulls, but the amounts were generally far below the total mineral accretion observed in seeds. When YSL1 and YSL3 are nonfunctional, Cu, Fe, and Zn are not effectively remobilized from, or do not effectively pass through, leaf and maternal fruit tissues. With respect to seed mineral accretion in Arabidopsis, continued uptake and translocation of minerals to source tissues during seed fill are as important, if not more important, than remobilization of previously stored minerals.

A latitudinal cline and response to vernalization in leaf angle and morphology in Arabidopsis thaliana (Brassicaceae)

Adaptation to latitudinal patterns of environmental variation is predicted to result in clinal variation in leaf traits. Therefore, this study tested for geographic differentiation and plastic responses to vernalization in leaf angle and leaf morphology in Arabidopsis thaliana. Twenty-one European ecotypes were grown in a common growth chamber environment. Replicates of each ecotype were exposed to one of four treatments: 0, 10, 20 or 30 d of vernalization. Ecotypes from lower latitudes had more erect leaves, as predicted from functional arguments about selection to maximize photosynthesis. Lower-latitude ecotypes also had more elongated petioles as predicted by a biomechanical constraint hypothesis. In addition, extended vernalization resulted in shorter and more erect leaves. As predicted by functional and adaptive hypotheses, our results show genetically based clinal variation as well as environmentally induced variation in leaf traits.

Demographic history has influenced nucleotide diversity in European Pinus sylvestris populations

To infer the role of natural selection in shaping standing genetic diversity, it is necessary to assess the genomewide impact of demographic history on nucleotide diversity. In this study we analyzed sequence diversity of 16 nuclear loci in eight Pinus sylvestris populations. Populations were divided into four geographical groups on the basis of their current location and the geographical history of the region: northern Europe, central Europe, Spain, and Turkey. There were no among-group differences in the level of silent nucleotide diversity, which was approximately 0.005/bp in all groups. There was some evidence that linkage disequilibrium extended further in northern Europe than in central Europe: the estimates of the population recombination rate parameter, rho, were 0.0064 and 0.0294, respectively. The summary statistics of nucleotide diversity in central and northern European populations were compatible with an ancient bottleneck rather than the standard neutral model.

Antagonistic pleiotropic effects reduce the potential adaptive value of the FRIGIDA locus

Although the occurrence of epistasis and pleiotropy is widely accepted at the molecular level, its effect on the adaptive value of fitness-related genes is rarely investigated in plants. Knowledge of these features of a gene is critical to understand the molecular basis of adaptive evolution. Here we investigate the importance of pleiotropy and epistasis in determining the adaptive value of a candidate gene using the gene FRI (FRIGIDA), which is thought to be the major gene controlling flowering time variation in Arabidopsis thaliana. The effect of FRI on flowering time was analyzed in an outbred population created by randomly mating 19 natural accessions of A. thaliana. This unique population allows the estimation of FRI effects independent of any linkage association with other loci due to demographic processes or to coadapted genes. It also allows for the estimation of pleiotropic effects of FRI on fitness and inflorescence architecture. We found that FRI explains less variation in flowering time than previously observed among natural accessions, and interacts epistatically with the FLC locus. Although early flowering plants produce more fruits under spring conditions, and nonfunctional alleles of FRI were associated with early flowering, variation at FRI was not associated with fitness. We show that nonfunctional FRI alleles have negative pleiotropic effects on fitness by reducing the numbers of nodes and branches on the inflorescence. We propose that these antagonistic pleiotropic effects reduce the adaptive value of FRI, and helps explain the maintenance of alternative life history strategies across natural populations of A. thaliana.

EVIDENCE OF LOCAL ADAPTATION TO COARSE-GRAINED ENVIRONMENTAL VARIATION IN ARABIDOPSIS THALIANA

Plants can achieve an appropriate phenotype in particular conditions either constitutively or plastically, depending in part on the grain size of the environmental conditions being considered. Coarse-grained environmental variation should result in selection for local adaptation and no selection on plasticity to novel levels of the coarse-grained environmental factors. We tested the hypotheses that natural populations of the well-studied model system Arabidopsis thaliana are locally adapted to spatially coarse-grained environmental variation, and that the photoperiodic regime per se is at least partially responsible for that local adaptation, by exposing natural populations to photoperiodic regimes characteristic of their native and foreign (novel) environments. We also tested the hypothesis that plasticity to novel photoperiodic regimes should appear random. We found that populations showed evidence of local adaptation at a spatially coarse grain, although not to photoperiodic regime per se. We also found that the plasticities to novel photoperiodic regimes appeared random and did not generally show evidence of adaptive divergence. Our study highlights the need for caution in extrapolating from the finding of local adaptation to the causes of local adaptation.

Using multilocus sequence data to assess population structure, natural selection, and linkage disequilibrium in wild tomatoes

We employed a multilocus approach to examine the effects of population subdivision and natural selection on DNA polymorphism in 2 closely related wild tomato species (Solanum peruvianum and Solanum chilense), using sequence data for 8 nuclear loci from populations across much of the species' range. Both species exhibit substantial levels of nucleotide variation. The species-wide level of silent nucleotide diversity is 18% higher in S. peruvianum (pi(sil) approximately 2.50%) than in S. chilense (pi(sil) approximately 2.12%). One of the loci deviates from neutral expectations, showing a clinal pattern of nucleotide diversity and haplotype structure in S. chilense. This geographic pattern of variation is suggestive of an incomplete (ongoing) selective sweep, but neutral explanations cannot be entirely dismissed. Both wild tomato species exhibit moderate levels of population differentiation (average F(ST) approximately 0.20). Interestingly, the pooled samples (across different demes) exhibit more negative Tajima's D and Fu and Li's D values; this marked excess of low-frequency polymorphism can only be explained by population (or range) expansion and is unlikely to be due to population structure per se. We thus propose that population structure and population/range expansion are among the most important evolutionary forces shaping patterns of nucleotide diversity within and among demes in these wild tomatoes. Patterns of population differentiation may also be impacted by soil seed banks and historical associations mediated by climatic cycles. Intragenic linkage disequilibrium (LD) decays very rapidly with physical distance, suggesting high recombination rates and effective population sizes in both species. The rapid decline of LD seems very promising for future association studies with the purpose of mapping functional variation in wild tomatoes.

Natural variation in innate immunity of a pioneer species

By 2010, we will have detailed knowledge about the genome of Arabidopsis thaliana from a Linnean-like effort by an international research community to identify nearly all of the genes in the species and to classify the products that these genes encode according to a primary function in a generic plant cell. To know the wild species, however, we will require knowledge of which genes provide the raw material for phenotypic variation and natural selection, and consequently affect the adaptability of individual plants and local populations across their geographic range, and ultimately survival of the species. Natural variation in innate immunity will be at the forefront of this exciting research frontier as a model for the molecular ecology of plant-microbe interactions.

Genetic diversity of the natural populations of Arabidopsis thaliana in China

Although extensive studies have been conducted on the genetic structure of Arabidopsis thaliana (A. thaliana) populations worldwide, the populations from China have never been studied. In this study, we collected 560 individuals from 19 natural populations of A. thaliana distributed in East China along the lower reaches of the Yangtze River, and two populations from northwest China (Xinjiang Province). We adopted two kinds of molecular marker, inter-simple sequence repeats (ISSRs) and random amplified polymorphic DNA (RAPDs) to investigate the genetic diversity within and among populations, and the correlation between the genetic and geographic distances. Thirteen ISSR primers produced 165 polymorphic bands (PPB) (96%) and 11 RAPD primers produced 162 polymorphic bands (98%) in about 560 individuals. The two marker systems generated similar patterns of genetic diversity in these natural populations. The AMOVA analysis indicated about 42-45% of the total genetic variation existed within populations, and found possible geographic structure. The Mantel test revealed a significant correlation between the geographic distance and the genetic distance of these populations in general. A close genetic relationship was found among four populations in the Jiangxi Province, and these always appeared clustered together as a monophyletic group in unweighted pair-group method with arithmetic averages dendrograms based on both ISSR and RAPD data sets. Based on the observation of recolonization and extinction of naturally distributed populations of A. thaliana, and the pattern of their genetic differentiation, the distribution of this species in China might be a result of natural dispersal under the strong influence of human activity.

Temporal dynamics and linkage disequilibrium in natural Caenorhabditis elegans populations

Caenorhabditis elegans is a major laboratory model system yet a newcomer to the field of population genetics, and relatively little is known of its biology in the wild. Recent studies of natural populations at a single time point revealed strong spatial population structure and suggested that these populations may be very dynamic. We have therefore studied several natural C. elegans populations over time and genotyped them at polymorphic microsatellite loci. While some populations appear to be genetically stable over the course of observation, others seem to go extinct, with full replacement of multilocus genotypes upon regrowth. The frequency of heterozygotes indicates that outcrossing occurs at a mean frequency of 1.7% and is variable between populations. However, in genetically stable populations, linkage disequilibrium between different chromosomes can be maintained over several years at a level much higher than expected from the heterozygote frequency. C. elegans seems to follow metapopulation dynamics, and the maintenance of linkage disequilibrium despite a low yet significant level of outcrossing suggests that selection may act against the progeny of outcrossings.

Fitness Effects Associated with the Major Flowering Time GeneFRIGIDAinArabidopsis thalianain the Field

To date, the effect of natural selection on candidate genes underlying complex traits has rarely been studied experimentally, especially under ecologically realistic conditions. Here we report that the effect of selection on the flowering time gene FRIGIDA (FRI) reverses depending on the season of germination and allelic variation at the interacting gene FLOWERING LOCUS C (FLC). In field studies of 136 European accessions of Arabidopsis thaliana, accessions with putatively functional FRI alleles had higher winter survival in one FLC background in a fall-germinating cohort, but accessions with deletion null FRI alleles had greater seed production in the other FLC background in a spring-germinating cohort. Consistent with FRI's role in flowering, selection analyses suggest that the difference in winter survival can be attributed to time to bolting. However, in the spring cohort, the fitness difference was associated with rosette size. Our analyses also reveal that controlling for population structure with estimates of inferred ancestry and a geographical restriction was essential for detecting fitness associations. Overall, our results suggest that the combined effects of seasonally varying selection and epistasis could explain the maintenance of variation at FRI and, more generally, may be important in the evolution of genes underlying complex traits.

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.

Combining population genomics and quantitative genetics: finding the genes underlying ecologically important traits

A central challenge in evolutionary biology is to identify genes underlying ecologically important traits and describe the fitness consequences of naturally occurring variation at these loci. To address this goal, several novel approaches have been developed, including 'population genomics,' where a large number of molecular markers are scored in individuals from different environments with the goal of identifying markers showing unusual patterns of variation, potentially due to selection at linked sites. Such approaches are appealing because of (1) the increasing ease of generating large numbers of genetic markers, (2) the ability to scan the genome without measuring phenotypes and (3) the simplicity of sampling individuals without knowledge of their breeding history. Although such approaches are inherently applicable to non-model systems, to date these studies have been limited in their ability to uncover functionally relevant genes. By contrast, quantitative genetics has a rich history, and more recently, quantitative trait locus (QTL) mapping has had some success in identifying genes underlying ecologically relevant variation even in novel systems. QTL mapping, however, requires (1) genetic markers that specifically differentiate parental forms, (2) a focus on a particular measurable phenotype and (3) controlled breeding and maintenance of large numbers of progeny. Here we present current advances and suggest future directions that take advantage of population genomics and quantitative genetic approaches - in both model and non-model systems. Specifically, we discuss advantages and limitations of each method and argue that a combination of the two provides a powerful approach to uncovering the molecular mechanisms responsible for adaptation.

Variation of enzyme activities and metabolite levels in 24 Arabidopsis accessions growing in carbon-limited conditions

Our understanding of the interaction of carbon (C) metabolism with nitrogen (N) metabolism and growth is based mainly on studies of responses to environmental treatments, and studies of mutants and transformants. Here, we investigate which metabolic parameters vary and which parameters change in a coordinated manner in 24 genetically diverse Arabidopsis (Arabidopsis thaliana) accessions, grown in C-limited conditions. The accessions were grown in short days, moderate light, and high nitrate, and analyzed for rosette biomass, levels of structural components (protein, chlorophyll), total phenols and major metabolic intermediates (sugars, starch, nitrate, amino acids), and the activities of seven representative enzymes from central C and N metabolism. The largest variation was found for plant weight, reducing sugars, starch at the end of the night, and several enzyme activities. High levels of one sugar correlated with high levels of other sugars and starch, and a trend to increased amino acids, slightly lower nitrate, and higher protein. The activities of enzymes at the interface of C and N metabolism correlated with each other, but were unrelated to carbohydrates, amino acid levels, and total protein. Rosette weight was unrelated or showed a weak negative trend to sugar and amino acid contents at the end of the day in most of the accessions, and was negatively correlated with starch at the end of the night. Rosette weight was positively correlated with several enzyme activities. We propose that growth is not related to the absolute levels of starch, sugars, and amino acids; instead, it is related to flux, which is indicated by the enzymatic capacity to use these central resources.

Comparative phylogeography of unglaciated eastern North America

Regional phylogeographical studies involving co-distributed animal and plant species have been conducted for several areas, most notably for Europe and the Pacific Northwest of North America. Until recently, phylogeographical studies in unglaciated eastern North America have been largely limited to animals. As more studies emerge for diverse lineages (including plants), it seems timely to assess the phylogeography across this region: (i) comparing and contrasting the patterns seen in plants and animals; (ii) assessing the extent of pseudocongruence; and (iii) discussing the potential applications of regional phylogeography to issues in ecology, such as response to climatic change. Unglaciated eastern North America is a large, geologically and topographically complex area with the species examined having diverse distributions. Nonetheless, some recurrent patterns emerge: (i) maritime - Atlantic vs. Gulf Coast; (ii) Apalachicola River discontinuity; (iii) Tombigbee River discontinuity; (iv) the Appalachian Mountain discontinuity; (v) the Mississippi River discontinuity; and (vi) the Apalachicola River and Mississippi River discontinuities. Although initially documented in animals, most of these patterns are also apparent in plants, providing support for phylogeographical generalizations. These patterns may generally be attributable to isolation and differentiation during Pleistocene glaciation, but in some cases may be older (Pliocene). Molecular studies sometimes agree with longstanding hypotheses of glacial refugia, but also suggest additional possible refugia, such as the southern Appalachian Mountains and areas close to the Laurentide Ice Sheet. Many species exhibit distinct patterns that reflect the unique, rather than the shared, aspects of species' phylogeographical histories. Furthermore, similar modern phylogeographical patterns can result from different underlying causal factors operating at different times (i.e. pseudocongruence). One underemphasized component of pseudocongruence may result from the efforts of researchers to categorize patterns visually - similar patterns may, in fact, not fully coincide, and inferring agreement may obscure the actual patterns and lead to erroneous conclusions. Our modelling analyses indicate no clear spatial patterning and support the hypothesis that phylogeographical structure in diverse temperate taxa is complex and was not shaped by just a few barriers.

Segregation distortion in Arabidopsis C24/Col-0 and Col-0/C24 recombinant inbred line populations is due to reduced fertility caused by epistatic interaction of two loci

A new large set of reciprocal recombinant inbred lines (RILs) was created between the Arabidopsis accessions Col-0 and C24 for quantitative trait mapping approaches, consisting of 209 Col-0 x C24 and 214 C24 x Col-0 F(7 )RI lines. Genotyping was performed using 110 evenly distributed framework single nucleotide polymorphism markers, yielding a genetic map of 425.70 cM, with an average interval of 3.87 cM. Segregation distortion (SD) was observed in several genomic regions during the construction of the genetic map. Linkage disequilibrium analysis revealed an association between a distorted region at the bottom of chromosome V and a non-distorted region on chromosome IV. A detailed analysis of the RILs for these two regions showed that an SD occurred when homozygous Col-0 alleles on chromosome IV coincided with homozygous C24 alleles at the bottom of chromosome V. Using nearly isogenic lines segregating for the distorted region we confirmed that this genotypic composition leads to reduced fertility and fitness.

Population genetic structure of Arabidopsis lyrata in Europe

Population genetic theory predicts that the self-incompatible and perennial herb, Arabidopsis lyrata, will have a genetic structure that differs from the self-fertilizing, annual Arabidopsis thaliana. We quantified the genetic structure for eight populations of A. lyrata ssp. petraea in historically nonglaciated regions of central Europe. Analysis of 20 microsatellite loci for 344 individuals demonstrated that, in accordance with predictions, diploid populations had high genome-wide heterozygosity (H(O) = 0.48; H(E) = 0.52), high within-population diversity (83% of total) compatible with mutation-drift equilibrium, and moderate differentiation among populations (F(ST) = 0.17). Within a single population, the vast majority of genetic variability (92%) was found at the smallest spatial scale (< 3 m). Although there was no evidence of biparental inbreeding or clonal propagation at this scale (F(IS) = 0.003), significant fine-scale spatial autocorrelation indicated localized gene flow presumably due to gravity dispersed seeds (Sp = 0.018). Limited gene flow between isolated population clusters (regions) separated by hundreds of kilometres has given rise to an isolation by distance pattern of diversification, with low, but significant, differentiation among regions (F(ST) = 0.05). The maintenance of geographically widespread polymorphisms and uniformly high diversity throughout central Europe is consistent with periglacial survival of A. lyrata ssp. petraea north of the Alps in steppe-tundra habitats during the last glacial maximum. As expected of northern and previously glaciated localities, A. lyrata in Iceland was genetically less diverse and highly differentiated from central Europe (H(E) = 0.37; F(ST) = 0.27).

The PHYTOCHROME C photoreceptor gene mediates natural variation in flowering and growth responses of Arabidopsis thaliana

Light has an important role in modulating seedling growth and flowering time. We show that allelic variation at the PHYTOCHROME C (PHYC) photoreceptor locus affects both traits in natural populations of A. thaliana. Two functionally distinct PHYC haplotype groups are distributed in a latitudinal cline dependent on FRIGIDA, a locus that together with FLOWERING LOCUS C explains a large portion of the variation in A. thaliana flowering time. In a genome-wide scan for association of 65 loci with latitude, there was an excess of significant P values, indicative of population structure. Nevertheless, PHYC was the most strongly associated locus across 163 strains, suggesting that PHYC alleles are under diversifying selection in A. thaliana. Our work, together with previous findings, suggests that photoreceptor genes are major agents of natural variation in plant flowering and growth response.

Challenges of detecting directional selection after a bottleneck: lessons from Sorghum bicolor

Multilocus surveys of sequence variation can be used to identify targets of directional selection, which are expected to have reduced levels of variation. Following a population bottleneck, the signal of directional selection may be hard to detect because many loci may have low variation by chance and the frequency spectrum of variation may be perturbed in ways that resemble the effects of selection. Cultivated Sorghum bicolor contains a subset of the genetic diversity found in its wild ancestor(s) due to the combined effects of a domestication bottleneck and human selection on traits associated with agriculture. As a framework for distinguishing between the effects of demography and selection, we sequenced 204 loci in a diverse panel of 17 cultivated S. bicolor accessions. Genomewide patterns of diversity depart strongly from equilibrium expectations with regard to the variance of the number of segregating sites, the site frequency spectrum, and haplotype configuration. Furthermore, gene genealogies of most loci with an excess of low frequency variants and/or an excess of segregating sites do not show the characteristic signatures of directional and diversifying selection, respectively. A simple bottleneck model provides an improved but inadequate fit to the data, suggesting the action of other population-level factors, such as population structure and migration. Despite a known history of recent selection, we find little evidence for directional selection, likely due to low statistical power and lack of an appropriate null model.

Genetic mechanisms and evolutionary significance of natural variation in Arabidopsis

Genomic studies of natural variation in model organisms provide a bridge between molecular analyses of gene function and evolutionary investigations of adaptation and natural selection. In the model plant species Arabidopsis thaliana, recent studies of natural variation have led to the identification of genes underlying ecologically important complex traits, and provided new insights about the processes of genome evolution, geographic population structure, and the selective mechanisms shaping complex trait variation in natural populations. These advances illustrate the potential for a new synthesis to elucidate mechanisms for the adaptive evolution of complex traits from nucleotide sequences to real-world environments.