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

Confounding of linkage disequilibrium patterns in large scale DNA based gene-gene interaction studies

Background

In Genome-Wide Association Studies (GWAS), the concept of linkage disequilibrium is important as it allows identifying genetic markers that tag the actual causal variants. In Genome-Wide Association Interaction Studies (GWAIS), similar principles hold for pairs of causal variants. However, Linkage Disequilibrium (LD) may also interfere with the detection of genuine epistasis signals in that there may be complete confounding between Gametic Phase Disequilibrium (GPD) and interaction. GPD may involve unlinked genetic markers, even residing on different chromosomes. Often GPD is eliminated in GWAIS, via feature selection schemes or so-called pruning algorithms, to obtain unconfounded epistasis results. However, little is known about the optimal degree of GPD/LD-pruning that gives a balance between false positive control and sufficient power of epistasis detection statistics. Here, we focus on Model-Based Multifactor Dimensionality Reduction as one large-scale epistasis detection tool. Its performance has been thoroughly investigated in terms of false positive control and power, under a variety of scenarios involving different trait types and study designs, as well as error-free and noisy data, but never with respect to multicollinear SNPs.

Results

Using real-life human LD patterns from a homogeneous subpopulation of British ancestry, we investigated the impact of LD-pruning on the statistical sensitivity of MB-MDR. We considered three different non-fully penetrant epistasis models with varying effect sizes. There is a clear advantage in pre-analysis pruning using sliding windows at r² of 0.75 or lower, but using a threshold of 0.20 has a detrimental effect on the power to detect a functional interactive SNP pair (power < 25%). Signal sensitivity, directly using LD-block information to determine whether an epistasis signal is present or not, benefits from LD-pruning as well (average power across scenarios: 87%), but is largely hampered by functional loci residing at the boundaries of an LD-block.

Conclusions

Our results confirm that LD patterns and the position of causal variants in LD blocks do have an impact on epistasis detection, and that pruning strategies and LD-blocks definitions combined need careful attention, if we wish to maximize the power of large-scale epistasis screenings.

Archaic lineages broaden our view on the history of Arabidopsis thaliana

Contents Summary 1194 I. Introduction 1194 II. Origin of the A. thaliana species 1194 III. The classic model of the history of A. thaliana 1195 IV. New genomic data from outside Eurasia challenge our view of A. thaliana history 1195 V. Conclusions 1197 Acknowledgements 1197 References 1197 SUMMARY: Natural variation in Arabidopsis thaliana has contributed to discoveries in diverse areas of plant biology. While A. thaliana has typically been considered a weed associated primarily with human-mediated environments, including agricultural and urban sites and railways, it has recently been shown that it is also native in remote natural areas, including high altitude sites in Eurasia and Africa, from the Atlas mountains in Morocco to the afro-alpine regions in Eastern and South Africa to Yunnan in China, the Himalayas and the Tibetan Plateau. This finding suggests that while A. thaliana has been extensively studied in Europe and Western Asia there are still many open questions about its population history, genotype-phenotype relationships and mechanisms of adaptation.

Natural variation in life history strategy of Arabidopsis thaliana determines stress responses to drought and insects of different feeding guilds

Plants are sessile organisms and, consequently, are exposed to a plethora of stresses in their local habitat. As a result, different populations of a species are subject to different selection pressures leading to adaptation to local conditions and intraspecific divergence. The annual brassicaceous plant Arabidopsis thaliana is an attractive model for ecologists and evolutionary biologists due to the availability of a large collection of resequenced natural accessions. Accessions of A. thaliana display one of two different life cycle strategies: summer and winter annuals. We exposed a collection of 308 European Arabidopsis accessions, that have been genotyped for 250K SNPs, to a range of stresses: one abiotic stress (drought), four biotic stresses (Pieris rapae caterpillars, Plutella xylostella caterpillars, Frankliniella occidentalis thrips and Myzus persicae aphids) and two combined stresses (drought plus P. rapae and Botrytis cinerea fungus plus P. rapae). We identified heritable genetic variation for responses to the different stresses, estimated by narrow-sense heritability. We found that accessions displaying different life cycle strategies differ in their response to stresses. Winter annuals are more resistant to drought, aphids and thrips and summer annuals are more resistant to P. rapae and P. xylostella caterpillars. Summer annuals are also more resistant to the combined stresses of drought plus P. rapae and infection by the fungus Botryris cinerea plus herbivory by P. rapae. Adaptation to drought displayed a longitudinal gradient. Finally, trade-offs were recorded between the response to drought and responses to herbivory by caterpillars of the specialist herbivore P. rapae.

QTL analysis of the developmental response to L-glutamate in Arabidopsis roots and its genotype-by-environment interactions

Primary root growth in Arabidopsis and a number of other species has previously been shown to be remarkably sensitive to the presence of external glutamate, with glutamate signalling eliciting major changes in root architecture. Using two recombinant inbred lines from reciprocal crosses between Arabidopsis accessions C24 and Col-0, we have identified one large-effect quantitative trait locus (QTL), GluS1, and two minor QTLs, GluS2 and GluS3, which together accounted for 41% of the phenotypic variance in glutamate sensitivity. The presence of the GluS1 locus on chromosome 3 was confirmed using a set of C24/Col-0 isogenic lines. GluS1 was mapped to an interval between genes At3g44830-At3g46880. When QTL mapping was repeated under a range of environmental conditions, including temperature, shading and nitrate supply, a strong genotype-by-environment interaction in the controls for the glutamate response was identified. Major differences in the loci controlling this trait were found under different environmental conditions. Here we present evidence for the existence of loci on chromosomes 1 and 5 epistatically controlling the response of the GluS1 locus to variations in ambient temperature, between 20°C and 26°C. In addition, a locus on the long arm of chromosome 1 was found to play a major role in controlling the ability of external nitrate signals to antagonize the glutamate effect. We conclude that there are multiple loci controlling natural variation in glutamate sensitivity in Arabidopsis roots and that epistatic interactions play an important role in modulating glutamate sensitivity in response to changes in environmental conditions.

Natural haplotypes of FLM non-coding sequences fine-tune flowering time in ambient spring temperatures in Arabidopsis

Cool ambient temperatures are major cues determining flowering time in spring. The mechanisms promoting or delaying flowering in response to ambient temperature changes are only beginning to be understood. In Arabidopsis thaliana, FLOWERING LOCUS M (FLM) regulates flowering in the ambient temperature range and FLM is transcribed and alternatively spliced in a temperature-dependent manner. We identify polymorphic promoter and intronic sequences required for FLM expression and splicing. In transgenic experiments covering 69% of the available sequence variation in two distinct sites, we show that variation in the abundance of the FLM-ß splice form strictly correlate (R² = 0.94) with flowering time over an extended vegetative period. The FLM polymorphisms lead to changes in FLM expression (PRO2+) but may also affect FLM intron 1 splicing (INT6+). This information could serve to buffer the anticipated negative effects on agricultural systems and flowering that may occur during climate change.

DOI:http://dx.doi.org/10.7554/eLife.22114.001

Field measurements of genotype by environment interaction for fitness caused by spontaneous mutations in Arabidopsis thaliana

As the ultimate source of genetic diversity, spontaneous mutation is critical to the evolutionary process. The fitness effects of spontaneous mutations are almost always studied under controlled laboratory conditions rather than under the evolutionarily relevant conditions of the field. Of particular interest is the conditionality of new mutations-that is, is a new mutation harmful regardless of the environment in which it is found? In other words, what is the extent of genotype-environment interaction for spontaneous mutations? We studied the fitness effects of 25 generations of accumulated spontaneous mutations in Arabidopsis thaliana in two geographically widely separated field environments, in Michigan and Virginia. At both sites, mean total fitness of mutation accumulation lines exceeded that of the ancestors, contrary to the expected decrease in the mean due to new mutations but in accord with prior work on these MA lines. We observed genotype-environment interactions in the fitness effects of new mutations, such that the effects of mutations in Michigan were a poor predictor of their effects in Virginia and vice versa. In particular, mutational variance for fitness was much larger in Virginia compared to Michigan. This strong genotype-environment interaction would increase the amount of genetic variation maintained by mutation-selection balance.

Genetic diversity and population structure of Arabidopsis thaliana along an altitudinal gradient

The natural genetic variation within a plant species is primarily a consequence of its phylogeography and evolutionary history. This variation largely determines its present-day population structure. Arabidopsis thaliana, as a model plant, has been studied in great detail including its probable origin, local as well as global genetic diversity pattern, population structure, adaptation, etc. However, no such studies have so far been reported from the Indian Himalayan region. Here, we describe a comprehensive study on the genetic diversity and population structure of A. thaliana from an altitudinal range of 700-3400 m above mean sea level the highest altitudinal range reported so far. We also compare these populations with previously reported worldwide populations. A total of 48 accessions representing six populations were analysed using 19 microsatellites and 11 chloroplast markers. Genetic diversity analysis indicated populations to be highly diverse and comparable with worldwide populations. STRUCTURE, principal coordinate and isolation by distance (IBD) analyses showed that genetic variation in different populations is structured at geographical and altitudinal level. Further analyses indicate that these populations are genetically distinct from the rest of the world populations. Different parameters of the demographic expansion model support a rapid expansion. Based on mismatch distribution, the initial time of expansion of west Himalayan populations was found to be about 130 000 years. Bayesian analysis of divergence time indicated that these populations have a long evolutionary history in this region. Based on the results of genetic diversity parameters, demographic expansion and divergence time estimation, it appears that west Himalayan populations may be the source of the west-east expansion model.

Global changes in gene expression, assayed by microarray hybridization and quantitative RT-PCR, during acclimation of three Arabidopsis thaliana accessions to sub-zero temperatures after cold acclimation

During cold acclimation plants increase in freezing tolerance in response to low non-freezing temperatures. This is accompanied by many physiological, biochemical and molecular changes that have been extensively investigated. In addition, plants of many species, including Arabidopsis thaliana, become more freezing tolerant during exposure to mild, non-damaging sub-zero temperatures after cold acclimation. There is hardly any information available about the molecular basis of this adaptation. Here, we have used microarrays and a qRT-PCR primer platform covering 1,880 genes encoding transcription factors (TFs) to monitor changes in gene expression in the Arabidopsis accessions Columbia-0, Rschew and Tenela during the first 3 days of sub-zero acclimation at -3 °C. The results indicate that gene expression during sub-zero acclimation follows a tighly controlled time-course. Especially AP2/EREBP and WRKY TFs may be important regulators of sub-zero acclimation, although the CBF signal transduction pathway seems to be less important during sub-zero than during cold acclimation. Globally, we estimate that approximately 5% of all Arabidopsis genes are regulated during sub-zero acclimation. Particularly photosynthesis-related genes are down-regulated and genes belonging to the functional classes of cell wall biosynthesis, hormone metabolism and RNA regulation of transcription are up-regulated. Collectively, these data provide the first global analysis of gene expression during sub-zero acclimation and allow the identification of candidate genes for forward and reverse genetic studies into the molecular mechanisms of sub-zero acclimation.

Lagging adaptation to warming climate in Arabidopsis thaliana

If climate change outpaces the rate of adaptive evolution within a site, populations previously well adapted to local conditions may decline or disappear, and banked seeds from those populations will be unsuitable for restoring them. However, if such adaptational lag has occurred, immigrants from historically warmer climates will outperform natives and may provide genetic potential for evolutionary rescue. We tested for lagging adaptation to warming climate using banked seeds of the annual weed Arabidopsis thaliana in common garden experiments in four sites across the species' native European range: Valencia, Spain; Norwich, United Kingdom; Halle, Germany; and Oulu, Finland. Genotypes originating from geographic regions near the planting site had high relative fitness in each site, direct evidence for broad-scale geographic adaptation in this model species. However, genotypes originating in sites historically warmer than the planting site had higher average relative fitness than local genotypes in every site, especially at the northern range limit in Finland. This result suggests that local adaptive optima have shifted rapidly with recent warming across the species' native range. Climatic optima also differed among seasonal germination cohorts within the Norwich site, suggesting that populations occurring where summer germination is common may have greater evolutionary potential to persist under future warming. If adaptational lag has occurred over just a few decades in banked seeds of an annual species, it may be an important consideration for managing longer-lived species, as well as for attempts to conserve threatened populations through ex situ preservation.

Seed dispersal in time can counteract the effect of gene flow between natural populations of Arabidopsis thaliana

Plants may escape unfavorable environments by dispersing to new sites, or by remaining in an ungerminated state at a given site until environmental conditions become favorable. There is limited evidence regarding the occurrence, interplay and relative importance of dispersal processes in time and space in plant populations. Thirty-six natural populations of the annual ruderal species Arabidopsis thaliana were monitored over five consecutive years, sampling both seed bank and above-ground cohorts. We show that immigration rates are considerably higher than previously inferred, averaging 1.7% per population yr(-1). On the other hand, almost one-third of the individuals in a given above-ground cohort result from seeds shed 2 or 3 yr back in time in 10 of the studied populations. Populations that disappeared one year were recolonized by regeneration from the seed bank the subsequent year. Thus, dispersal in both time and space is an important contributor to the structuring of genetic variability in natural populations of A. thaliana, where a high dispersal rate in time may partly counteract the homogenizing effects of spatial seed and pollen dispersal.

Cardamine hirsuta: a versatile genetic system for comparative studies

A major goal in biology is to identify the genetic basis for phenotypic diversity. This goal underpins research in areas as diverse as evolutionary biology, plant breeding and human genetics. A limitation for this research is no longer the availability of sequence information but the development of functional genetic tools to understand the link between changes in sequence and phenotype. Here we describe Cardamine hirsuta, a close relative of the reference plant Arabidopsis thaliana, as an experimental system in which genetic and transgenic approaches can be deployed effectively for comparative studies. We present high-resolution genetic and cytogenetic maps for C. hirsuta and show that the genome structure of C. hirsuta closely resembles the eight chromosomes of the ancestral crucifer karyotype and provides a good reference point for comparative genome studies across the Brassicaceae. We compared morphological and physiological traits between C. hirsuta and A. thaliana and analysed natural variation in stamen number in which lateral stamen loss is a species characteristic of C. hirsuta. We constructed a set of recombinant inbred lines and detected eight quantitative trait loci that can explain stamen number variation in this population. We found clear phylogeographic structure to the genetic variation in C. hirsuta, thus providing a context within which to address questions about evolutionary changes that link genotype with phenotype and the environment.

The genetic structure of Arabidopsis thaliana in the south-western Mediterranean range reveals a shared history between North Africa and southern Europe

BACKGROUND

Deciphering the genetic structure of Arabidopsis thaliana diversity across its geographic range provides the bases for elucidating the demographic history of this model plant. Despite the unique A. thaliana genomic resources currently available, its history in North Africa, the extreme southern limit in the biodiversity hotspot of the Mediterranean Basin, remains virtually unknown.

RESULTS

To approach A. thaliana evolutionary history in North Africa, we have analysed the genetic diversity and structure of 151 individuals collected from 20 populations distributed across Morocco. Genotyping of 249 genome-wide SNPs indicated that Morocco contains substantially lower diversity than most analyzed world regions. However, IBD, STRUCTURE and PCA clustering analyses showed that genetic variation is strongly geographically structured. We also determined the genetic relationships between Morocco and the closest European region, the Iberian Peninsula, by analyses of 201 populations from both regions genotyped with the same SNPs. These analyses detected four genetic groups, but all Moroccan accessions belonged to a common Iberian/Moroccan cluster that appeared highly differentiated from the remaining groups. Thus, we identified a genetic lineage with an isolated demographic history in the south-western Mediterranean region. The existence of this lineage was further supported by the study of several flowering genes and traits, which also found Moroccan accessions similar to the same Iberian group. Nevertheless, genetic diversity for neutral SNPs and flowering genes was higher in Moroccan than in Iberian populations of this lineage. Furthermore, we analyzed the genetic relationships between Morocco and other world regions by joint analyses of a worldwide collection of 337 accessions, which detected an additional weak relationship between North Africa and Asia.

CONCLUSIONS

The patterns of genetic diversity and structure of A. thaliana in Morocco show that North Africa is part of the species native range and support the occurrence of a glacial refugium in the Atlas Mountains. In addition, the identification of a genetic lineage specific of Morocco and the Iberian Peninsula indicates that the Strait of Gibraltar has been an A. thaliana migration route between Europe and Africa. Finally, the genetic relationship between Morocco and Asia suggests another migration route connecting north-western Africa and Asia.

Mapping quantitative trait loci for freezing tolerance in a recombinant inbred line population of Arabidopsis thaliana accessions Tenela and C24 reveals REVEILLE1 as negative regulator of cold acclimation

The ability to increase freezing tolerance when exposed to low temperatures is a property of many plant species from temperate climates and involves a wide array of metabolic adjustments and changes in gene expression. In Arabidopsis thaliana, natural accessions show high variation in their acclimation capacity, and freezing tolerance correlates with natural habitat temperatures. To investigate the genetic basis of this variation, a recombinant inbred line population from reciprocal crosses between the accessions C24 and Tenela (Te), showing large variation in tolerance, was established. Over 250 recombinant inbred lines were genotyped for 69 single nucleotide polymorphism markers in a linkage map with 391.9 centimorgans (cM) and phenotyped for their freezing tolerance using the electrolyte leakage method that reports cell damage after a freeze-thaw cycle. Mapping of quantitative trait loci (QTL) for acclimated plants revealed three QTL regions on chromosomes 2, 4 and 5. Based on gene expression data, QTL regions were screened for genes differentially responding to low temperature in C24 and Te. Among the candidate genes, the Myb family transcription factor REVEILLE1 (At5g17300) on chromosome 5 was identified as a novel negative regulator of freezing tolerance in Arabidopsis.

Signatures of demography and recombination at coding genes in naturally-distributed populations of Arabidopsis lyrata subsp. petraea

Demography impacts the observed standing level of genetic diversity present in populations. Distinguishing the relative impacts of demography from selection requires a baseline of expressed gene variation in naturally occurring populations. Six nuclear genes were sequenced to estimate the patterns and levels of genetic diversity in natural Arabidopsis lyrata subsp. petraea populations that differ in demographic histories since the Pleistocene. As expected, northern European populations have genetic signatures of a strong population bottleneck likely due to glaciation during the Pleistocene. Levels of diversity in the northern populations are about half of that in central European populations. Bayesian estimates of historical population size changes indicate that central European populations also have signatures of population size change since the last glacial maxima, suggesting that these populations are not as stable as previously thought. Time since divergence amongst northern European populations is higher than amongst central European populations, suggesting that the northern European populations were established before the Pleistocene and survived glaciation in small separated refugia. Estimates of demography based on expressed genes are complementary to estimates based on microsatellites and transposable elements, elucidating temporal shifts in population dynamics and confirming the importance of marker selection for tests of demography.

Intraspecific sequence variation and differential expression in starch synthase genes of Arabidopsis thaliana

BACKGROUND

Natural accessions of Arabidopsis thaliana are a well-known system to measure levels of intraspecific genetic variation. Leaf starch content correlates negatively with biomass. Starch is synthesized by the coordinated action of many (iso)enzymes. Quantitatively dominant is the repetitive transfer of glucosyl residues to the non-reducing ends of α-glucans as mediated by starch synthases. In the genome of A. thaliana, there are five classes of starch synthases, designated as soluble starch synthases (SSI, SSII, SSIII, and SSIV) and granule-bound synthase (GBSS). Each class is represented by a single gene. The five genes are homologous in functional domains due to their common origin, but have evolved individual features as well. Here, we analyze the extent of genetic variation in these fundamental protein classes as well as possible functional implications on transcript and protein levels.

FINDINGS

Intraspecific sequence variation of the five starch synthases was determined by sequencing the entire loci including promoter regions from 30 worldwide distributed accessions of A. thaliana. In all genes, a considerable number of nucleotide polymorphisms was observed, both in non-coding and coding regions, and several amino acid substitutions were identified in functional domains. Furthermore, promoters possess numerous polymorphisms in potentially regulatory cis-acting regions. By realtime experiments performed with selected accessions, we demonstrate that DNA sequence divergence correlates with significant differences in transcript levels.

CONCLUSIONS

Except for AtSSII, all starch synthase classes clustered into two or three groups of haplotypes, respectively. Significant difference in transcript levels among haplotype clusters in AtSSIV provides evidence for cis-regulation. By contrast, no such correlation was found for AtSSI, AtSSII, AtSSIII, and AtGBSS, suggesting trans-regulation. The expression data presented here point to a regulation by common trans-regulatory transcription factors which ensures a coordinated action of the products of these four genes during starch granule biosynthesis. The apparent cis-regulation of AtSSIV might be related to its role in the initiation of de novo biosynthesis of granules.

Among- and within-population variation in flowering time of Iberian Arabidopsis thaliana estimated in field and glasshouse conditions

The study of the evolutionary and population genetics of quantitative traits requires the assessment of within- and among-population patterns of variation. We carried out experiments including eight Iberian Arabidopsis thaliana populations (10 individuals per population) in glasshouse and field conditions. We quantified among- and within-population variation for flowering time and for several field life-history traits. Individuals were genotyped with microsatellites, single nucleotide polymorphisms and four well-known flowering genes (FRI, FLC, CRY2 and PHYC). Phenotypic and genotypic data were used to conduct Q(ST)-F(ST) comparisons. Life-history traits varied significantly among- and within-populations. Flowering time also showed substantial within- and among-population variation as well as significant genotype × environment interactions among the various conditions. Individuals bearing FRI truncations exhibited reduced recruitment in field conditions and differential flowering time behavior across experimental conditions, suggesting that FRI contributes to the observed significant genotype × environment interactions. Flowering time estimated in field conditions was the only trait showing significantly higher quantitative genetic differentiation than neutral genetic differentiation values. Overall, our results show that these A. thaliana populations are genetically more differentiated for flowering time than for neutral markers, suggesting that flowering time is likely to be under divergent selection.

Polymorphism at genes involved in salt tolerance in Arabidopsis thaliana (Brassicaceae)

PREMISE OF THE STUDY

Genes involved in relevant functions for environmental adaptation can be considered primary candidates for their variation having been shaped by natural selection. Detecting recent selective events through their footprint on nucleotide variation constitutes a challenging task in species with a complex demographic history such as Arabidopsis thaliana. We have surveyed nucleotide variation in this species at nine genes involved in salt tolerance. The available genomewide information for this species has allowed us to contrast the levels and patterns of variation detected at the candidate genes with empirical distributions obtained from noncandidate regions.

METHODS

We sequenced nine genes involved in salt tolerance (~32 kb) in 20 ecotypes of A. thaliana and analyzed polymorphism and divergence at the individual gene and multilocus levels.

KEY RESULTS

Variation at the nine genes studied was characterized by a generalized skew toward polymorphisms with low-frequency variants. Except for genes RCD1 and NHX8, this pattern was similar to that generally detected in the A. thaliana genome and could thus be primarily explained by the species demographic history. The more extreme deviation at the NHX8 gene and its excess of polymorphism relative to divergence points to the recent action of selection on this gene.

CONCLUSIONS

The analysis of nucleotide polymorphism and divergence at nine genes involved in salt tolerance provided little evidence for the recent action of positive selection. Only the signals detected at NHX8 from both polymorphism and divergence were suggestive of the putative contribution of this gene to local adaptation.

Natural variation in the temperature range permissive for vernalization in accessions of Arabidopsis thaliana

Vernalization is an acceleration of flowering in response to chilling, and is normally studied in the laboratory at near-freezing (2-4 °C) temperatures. Many vernalization-requiring species, such as Arabidopsis thaliana, are found in a range of habitats with varying winter temperatures. Natural variation in the temperature range that elicits a vernalization response in Arabidopsis has not been fully explored. We characterized the effect of intermediate temperatures (7-19 °C) on 15 accessions and the well-studied reference line Col-FRI. Although progressively warmer temperatures are gradually less effective at activating expression of the vernalization-specific gene VERNALIZATION-INSENSITIVE 3 (VIN3) and in accelerating flowering, there is substantial natural variation in the upper threshold (T(max) ) of the flowering-time response. VIN3 is required for the T(max) (13 °C) response of Col-FRI. Surprisingly, even 16 °C treatment caused induction of VIN3 in six tested lines, despite the ineffectiveness of this temperature in accelerating flowering for two of them. Finally, we present evidence that mild acceleration of flowering by 19 °C exposure may counterbalance the flowering time delay caused by non-inductive photoperiods in at least one accession, creating an appearance of photoperiod insensitivity.

Family-based linkage and association mapping reveals novel genes affecting Plum pox virus infection in Arabidopsis thaliana

Sharka is a devastating viral disease caused by the Plum pox virus (PPV) in stone fruit trees and few sources of resistance are known in its natural hosts. Since any knowledge gained from Arabidopsis on plant virus susceptibility factors is likely to be transferable to crop species, Arabidopsis's natural variation was searched for host factors essential for PPV infection. To locate regions of the genome associated with susceptibility to PPV, linkage analysis was performed on six biparental populations as well as on multiparental lines. To refine quantitative trait locus (QTL) mapping, a genome-wide association analysis was carried out using 147 Arabidopsis accessions. Evidence was found for linkage on chromosomes 1, 3 and 5 with restriction of PPV long-distance movement. The most relevant signals occurred within a region at the bottom of chromosome 3, which comprises seven RTM3-like TRAF domain-containing genes. Since the resistance mechanism analyzed here is recessive and the rtm3 knockout mutant is susceptible to PPV infection, it suggests that other gene(s) present in the small identified region encompassing RTM3 are necessary for PPV long-distance movement. In consequence, we report here the occurrence of host factor(s) that are indispensable for virus long-distance movement.

Clinal variation in the non-acclimated and cold-acclimated freezing tolerance of Arabidopsis thaliana accessions

Arabidopsis thaliana is a geographically widely spread species consisting of local accessions differing both genetically and phenotypically. These differences may constitute environmental adaptations and a latitudinal cline in freezing tolerance has been shown previously. Many plants, including Arabidopsis, exhibit increased freezing tolerance after cold exposure (cold acclimation). Here we present evidence for geographical clines (both latitudinal and longitudinal) in acclimated (ACC) and non-acclimated (NA) freezing tolerance, estimated from electrolyte leakage measurements on 54 accessions. Leaf Pro contents were not correlated with freezing tolerance, while sugar contents (Glc, Fru, Suc, Raf) were in the ACC, but not the NA state. Expression levels of 14 cold-induced genes were investigated before and after 2 weeks of cold acclimation by quantitative RT-PCR. Expression of the CBF1, 2 and 3 genes was not correlated with freezing tolerance. The expression of some CBF-regulated (COR) genes, however, was correlated specifically with ACC freezing tolerance. A tight correlation between CBF and COR gene expression was only observed under non-acclimating conditions, where CBF and COR expression were also correlated with the expression of PRR5, a component of the circadian clock. Collectively, this study sheds new light on the molecular determinants of plant-freezing tolerance and cold acclimation and their geographical dependence.

Novel natural alleles at FLC and LVR loci account for enhanced vernalization responses in Arabidopsis thaliana

Vernalization, the induction of flowering by low winter temperatures, is likely to be involved in plant climatic adaptation. However, the genetic, molecular and ecological bases underlying the quantitative variation that tunes vernalization sensitivity to natural environments are largely unknown. To address these questions, we have studied the enhanced vernalization response shown by the Ll-0 accession of Arabidopsis thaliana. Quantitative trait locus (QTL) mapping for several flowering initiation traits in relation to vernalization, in a new Ler × Ll-0 recombinant inbred line (RIL) population, identified large effect alleles at FRI, FLC and HUA2, together with two small effect loci named as Llagostera vernalization response (LVR) 1 and 2. Phenotypic analyses of near isogenic lines validated LVR1 effect on flowering vernalization responses. To further characterize the FLC allele from Ll-0, we carried out genetic association analyses using a regional collection of wild genotypes. FLC-Ll-0 appeared as a low-frequency allele that is distinguished by polymorphism Del(-57), a 50-bp-deletion in the 5'-UTR. Del(-57) was significantly associated with enhanced vernalization responses and FLC RNA expression, as well as with altitude and minimum temperatures. These results are consistent with Del(-57) acting as a novel cis-regulatory FLC polymorphism that may confer climatic adaptation by increasing vernalization sensitivity.

Longitudinal trends in climate drive flowering time clines in North American Arabidopsis thaliana

Introduced species frequently show geographic differentiation, and when differentiation mirrors the ancestral range, it is often taken as evidence of adaptive evolution. The mouse-ear cress (Arabidopsis thaliana) was introduced to North America from Eurasia 150–200 years ago, providing an opportunity to study parallel adaptation in a genetic model organism. Here, we test for clinal variation in flowering time using 199 North American (NA) accessions of A. thaliana, and evaluate the contributions of major flowering time genes FRI, FLC, and PHYC as well as potential ecological mechanisms underlying differentiation. We find evidence for substantial within population genetic variation in quantitative traits and flowering time, and putatively adaptive longitudinal differentiation, despite low levels of variation at FRI, FLC, and PHYC and genome-wide reductions in population structure relative to Eurasian (EA) samples. The observed longitudinal cline in flowering time in North America is parallel to an EA cline, robust to the effects of population structure, and associated with geographic variation in winter precipitation and temperature. We detected major effects of FRI on quantitative traits associated with reproductive fitness, although the haplotype associated with higher fitness remains rare in North America. Collectively, our results suggest the evolution of parallel flowering time clines through novel genetic mechanisms.

Characterizing genomic variation of Arabidopsis thaliana: the roles of geography and climate

Arabidopsis thaliana inhabits diverse climates and exhibits varied phenology across its range. Although A. thaliana is an extremely well-studied model species, the relationship between geography, growing season climate and its genetic variation is poorly characterized. We used redundancy analysis (RDA) to quantify the association of genomic variation [214 051 single nucleotide polymorphisms (SNPs)] with geography and climate among 1003 accessions collected from 447 locations in Eurasia. We identified climate variables most correlated with genomic variation, which may be important selective gradients related to local adaptation across the species range. Climate variation among sites of origin explained slightly more genomic variation than geographical distance. Large-scale spatial gradients and early spring temperatures explained the most genomic variation, while growing season and summer conditions explained the most after controlling for spatial structure. SNP variation in Scandinavia showed the greatest climate structure among regions, possibly because of relatively consistent phenology and life history of populations in this region. Climate variation explained more variation among nonsynonymous SNPs than expected by chance, suggesting that much of the climatic structure of SNP correlations is due to changes in coding sequence that may underlie local adaptation.

High throughput screening with chlorophyll fluorescence imaging and its use in crop improvement

Marker assisted plant breeding is a powerful technique for targeted crop improvement in horticulture and agriculture. It depends upon the correlation of desirable phenotypic characteristics with specific genetic markers. This can be determined by statistical models that relate the variation in the value of genetic markers to variation in phenotypic traits. It therefore depends upon the convergence of three technologies; the creation of genetically characterised (and thus marked) populations, high throughput screening procedures, and statistical procedures. While a large number of high throughput screening technologies are available, real-time screening techniques are usually based on some kind of imaging technologies, such as chlorophyll fluorescence imaging, that offers physiological data that are eminently suitable as a quantitative trait for genetic mapping.

Plant sexual reproduction during climate change: gene function in natura studied by ecological and evolutionary systems biology

BACKGROUND

It is essential to understand and predict the effects of changing environments on plants. This review focuses on the sexual reproduction of plants, as previous studies have suggested that this trait is particularly vulnerable to climate change, and because a number of ecologically and evolutionarily relevant genes have been identified.

SCOPE

It is proposed that studying gene functions in naturally fluctuating conditions, or gene functions in natura, is important to predict responses to changing environments. First, we discuss flowering time, an extensively studied example of phenotypic plasticity. The quantitative approaches of ecological and evolutionary systems biology have been used to analyse the expression of a key flowering gene, FLC, of Arabidopsis halleri in naturally fluctuating environments. Modelling showed that FLC acts as a quantitative tracer of the temperature over the preceding 6 weeks. The predictions of this model were verified experimentally, confirming its applicability to future climate changes. Second, the evolution of self-compatibility as exemplifying an evolutionary response is discussed. Evolutionary genomic and functional analyses have indicated that A. thaliana became self-compatible via a loss-of-function mutation in the male specificity gene, SCR/SP11. Self-compatibility evolved during glacial-interglacial cycles, suggesting its association with mate limitation during migration. Although the evolution of self-compatibility may confer short-term advantages, it is predicted to increase the risk of extinction in the long term because loss-of-function mutations are virtually irreversible.

CONCLUSIONS

Recent studies of FLC and SCR have identified gene functions in natura that are unlikely to be found in laboratory experiments. The significance of epigenetic changes and the study of non-model species with next-generation DNA sequencers is also discussed.

Whole-genome sequencing of multiple Arabidopsis thaliana populations

The plant Arabidopsis thaliana occurs naturally in many different habitats throughout Eurasia. As a foundation for identifying genetic variation contributing to adaptation to diverse environments, a 1001 Genomes Project to sequence geographically diverse A. thaliana strains has been initiated. Here we present the first phase of this project, based on population-scale sequencing of 80 strains drawn from eight regions throughout the species' native range. We describe the majority of common small-scale polymorphisms as well as many larger insertions and deletions in the A. thaliana pan-genome, their effects on gene function, and the patterns of local and global linkage among these variants. The action of processes other than spontaneous mutation is identified by comparing the spectrum of mutations that have accumulated since A. thaliana diverged from its closest relative 10 million years ago with the spectrum observed in the laboratory. Recent species-wide selective sweeps are rare, and potentially deleterious mutations are more common in marginal populations.

Transcriptome profiling and molecular marker discovery in red pepper, Capsicum annuum L. TF68

Transcriptome from high throughput sequencing-by-synthesis is a good resource of molecular markers. In this study, we present utility of massively parallel sequencing by synthesis for profiling the transcriptome of red pepper (Capsicum annuum L. TF68) using 454 GS-FLX pyrosequencing. Through the generation of approximately 30.63 megabases (Mb) of expressed sequence tag (EST) data with the average length of 375 base pairs (bp), 9,818 contigs and 23,712 singletons were obtained by raw reads assembly. Using BLAST alignment against NCBI non-redundant and a UniProt protein database, 30% of the tentative consensus sequences were assigned to specific function annotation, while 24% returned alignments of unknown function, leaving up to 46% with no alignment. Functional classification using FunCat revealed that sequences with putative known function were distributed cross 18 categories. All unigenes have an approximately equal distribution on chromosomes by aligning with tomato (Solanum lycopersicum) pseudomolecules. Furthermore, 1,536 high quality single nucleotide discrepancies were discovered using the Bukang mature fruit cDNA collection (dbEST ID: 23667) as a reference. Moreover, 758 simple sequence repeat (SSR) motif loci were mined from 614 contigs, from which 572 primer sets were designed. The SSR motifs corresponded to di- and tri- nucleotide motifs (27.03 and 61.92%, respectively). These molecular markers may be of great value for application in linkage mapping and association mapping research.

Unusual signatures of highly adaptable R-loci in closely-related Arabidopsis species

Plant resistance genes (R-genes) evolve rapidly in response to changing environments. What are the most remarkable signatures of fast adaptive genes, besides the commonly revealed rapid divergence and high non-synonymous substitution rate? Here we investigated these changes in five R-loci following recent differentiation between Arabidopsis thaliana and Arabidopsis lyrata. Extreme differences in evolutionary rates were observed: e.g., an overall 5.46-9.83-fold different nucleotide diversity at two R-loci between species, ten-fold higher non-synonymous substitution rates within one species versus the other, significantly different Ka/Ks ratios between species for the same R-gene, and high interspecific divergence at one R-locus. Particularly, we observed an elevated level of trans-specific polymorphism at one R-locus and a differentially maintained presence/absence polymorphism at another. The high frequency of ancestral polymorphisms amongst R-genes suggests that the persistence of some functional variation is an important evolutionary mechanism shaping genetic variation in R-genes, while the variation of presence/absence polymorphisms provides a potential mechanism for malleable activation of adaptive resistance response pathways. The distinct patterns among R-genes suggest that the same R-gene ortholog can be quickly shaped by different evolutionary processes, e.g., purifying selection in one species but positive selection in a closely-related species.

A gene duplication/loss event in the ribulose-1,5-bisphosphate-carboxylase/oxygenase (rubisco) small subunit gene family among accessions of Arabidopsis thaliana

Rubisco (ribulose-1,5-bisphosphate carboxylase/oxygenase; EC 4.1.1.39), the most abundant protein in nature, catalyzes the assimilation of CO(2) (worldwide about 10(11) t each year) by carboxylation of ribulose-1,5-bisphosphate. It is a hexadecamer consisting of eight large and eight small subunits. Although the Rubisco large subunit (rbcL) is encoded by a single gene on the multicopy chloroplast genome, the Rubisco small subunits (rbcS) are encoded by a family of nuclear genes. In Arabidopsis thaliana, the rbcS gene family comprises four members, that is, rbcS-1a, rbcS-1b, rbcS-2b, and rbcS-3b. We sequenced all Rubisco genes in 26 worldwide distributed A. thaliana accessions. In three of these accessions, we detected a gene duplication/loss event, where rbcS-1b was lost and substituted by a duplicate of rbcS-2b (called rbcS-2b*). By screening 74 additional accessions using a specific polymerase chain reaction assay, we detected five additional accessions with this duplication/loss event. In summary, we found the gene duplication/loss in 8 of 100 A. thaliana accessions, namely, Bch, Bu, Bur, Cvi, Fei, Lm, Sha, and Sorbo. We sequenced an about 1-kb promoter region for all Rubisco genes as well. This analysis revealed that the gene duplication/loss event was associated with promoter alterations (two insertions of 450 and 850 bp, one deletion of 730 bp) in rbcS-2b and a promoter deletion (2.3 kb) in rbcS-2b* in all eight affected accessions. The substitution of rbcS-1b by a duplicate of rbcS-2b (i.e., rbcS-2b*) might be caused by gene conversion. All four Rubisco genes evolve under purifying selection, as expected for central genes of the highly conserved photosystem of green plants. We inferred a single positive selected site, a tyrosine to aspartic acid substitution at position 72 in rbcS-1b. Exactly the same substitution compromises carboxylase activity in the cyanobacterium Anacystis nidulans. In A. thaliana, this substitution is associated with an inferred recombination. Functional implications of the substitution remain to be evaluated.

Population structure and genetic differentiation associated with breeding history and selection in tomato (Solanum lycopersicum L.)

Tomato (Solanum lycopersicum L.) has undergone intensive selection during and following domestication. We investigated population structure and genetic differentiation within a collection of 70 tomato lines representing contemporary (processing and fresh-market) varieties, vintage varieties and landraces. The model-based Bayesian clustering software, STRUCTURE, was used to detect subpopulations. Six independent analyses were conducted using all marker data (173 markers) and five subsets of markers based on marker type (single-nucleotide polymorphisms, simple sequence repeats and insertion/deletions) and location (exon and intron sequences) within genes. All of these analyses consistently separated four groups predefined by market niche and age into distinct subpopulations. Furthermore, we detected at least two subpopulations within the processing varieties. These subpopulations correspond to historical patterns of breeding conducted for specific production environments. We found no subpopulation within fresh-market varieties, vintage varieties and landraces when using all marker data. High levels of admixture were shown in several varieties representing a transition in the demarcation between processing and fresh-market breeding. The genetic clustering detected by using the STRUCTURE software was confirmed by two statistics, pairwise F(st) (θ) and Nei's standard genetic distance. We also identified a total of 19 loci under positive selection between processing, fresh-market and vintage germplasm by using an F(st)-outlier method based on the deviation from the expected distribution of F(st) and heterozygosity. The markers and genome locations we identified are consistent with known patterns of selection and linkage to traits that differentiate the market classes. These results demonstrate how human selection through breeding has shaped genetic variation within cultivated tomato.

Constitutive salicylic acid defences do not compromise seed yield, drought tolerance and water productivity in the Arabidopsis accession C24

Plants that constitutively express otherwise inducible disease resistance traits often suffer a depressed seed yield in the absence of a challenge by pathogens. This has led to the view that inducible disease resistance is indispensable, ensuring that minimal resources are diverted from growth, reproduction and abiotic stress tolerance. The Arabidopsis genotype C24 has enhanced basal resistance, which was shown to be caused by permanent expression of normally inducible salicylic acid (SA)-regulated defences. However, the seed yield of C24 was greatly enhanced in comparison to disease-resistant mutants that display identical expression of SA defences. Under both water-replete and -limited conditions, C24 showed no difference and increased seed yield, respectively, in comparison with pathogen-susceptible genotypes. C24 was the most drought-tolerant genotype and showed elevated water productivity, defined as seed yield per plant per millilitre water consumed, and achieved this by displaying adjustments to both its development and transpiration efficiency (TE). Therefore, constitutive high levels of disease resistance in C24 do not affect drought tolerance, seed yield and seed viability. This study demonstrates that it will be possible to combine traits that elevate basal disease resistance and improve water productivity in crop species, and such traits need not be mutually exclusive.

Genetic and physiological bases for phenological responses to current and predicted climates

We are now reaching the stage at which specific genetic factors with known physiological effects can be tied directly and quantitatively to variation in phenology. With such a mechanistic understanding, scientists can better predict phenological responses to novel seasonal climates. Using the widespread model species Arabidopsis thaliana, we explore how variation in different genetic pathways can be linked to phenology and life-history variation across geographical regions and seasons. We show that the expression of phenological traits including flowering depends critically on the growth season, and we outline an integrated life-history approach to phenology in which the timing of later life-history events can be contingent on the environmental cues regulating earlier life stages. As flowering time in many plants is determined by the integration of multiple environmentally sensitive gene pathways, the novel combinations of important seasonal cues in projected future climates will alter how phenology responds to variation in the flowering time gene network with important consequences for plant life history. We discuss how phenology models in other systems--both natural and agricultural--could employ a similar framework to explore the potential contribution of genetic variation to the physiological integration of cues determining phenology.

Extremely low genetic variability and highly structured local populations of Arabidopsis thaliana at higher latitudes

The genetic diversity and population structure of Arabidopsis thaliana populations from Norway were studied and compared to a worldwide sample of A. thaliana to investigate the demographic history and elucidate possible colonization routes of populations at the northernmost species limit. We genotyped 282 individuals from 31 local populations using 149 single nucleotide polymorphism markers. A high level of population subdivision (F(ST) = 0.85 ± 0.007) was found indicating that A. thaliana is highly structured at the regional level. Significant relationships between genetic and geographical distances were found, suggesting an isolation by distance mode of evolution. Genetic diversity was much lower, and the level of linkage disequilibrium was higher in populations from the north (65-68°N) compared to populations from the south (59-62°N); this is consistent with a northward expansion pattern. A neighbour-joining tree showed that populations from northern Norway form a separate cluster, while the remaining populations are distributed over a few minor clusters. Minimal gene flow seems to have occurred between populations in different regions, especially between the geographically distant northern and southern populations. Our data suggest that northern populations represent a homogenous group that may have been established from a few founders during northward expansions, while populations in the central part of Norway constitute an admixed group established by founders of different origins, most probably as a result of human-mediated gene flow. Moreover, Norwegian populations appeared to be homogenous and isolated compared to a worldwide sample of A. thaliana, but they are still grouped with Swedish populations, which may indicate common colonization histories.

Carotenoid biosynthesis genes provide evidence of geographical subdivision and extensive linkage disequilibrium in the carrot

According to the history of the cultivated carrot, root colour can be considered as a structural factor of carrot germplasm. Therefore, molecular variations of carotenoid biosynthesis genes, these being involved in colour traits, represent a good putative source of polymorphism related to diversity structure. Seven candidate genes involved in the carotenoid biosynthesis pathway have been analysed from a sample of 48 individual plants, each one from a different cultivar of carrot (Daucus carota L. ssp. sativus). The cultivars were chosen to represent a large diversity and a wide range of root colour. A high single nucleotide polymorphism (SNP) frequency of 1 SNP per 22 bp (mean pi (sil) = 0.020) was found on average within these genes. The analysis of genetic structure from carotenoid biosynthesis gene sequences and 17 putatively neutral microsatellites showed moderate genetic differentiation between cultivars originating from the West and the East (F (ST) = 0.072), this being consistent with breeding history, but not previously evidenced by molecular tools. Surprisingly, carotenoid biosynthesis genes did not exhibit decay of LD (mean r (2) = 0.635) within the 700-1,000 bp analysed, even though a fast decay level of LD is expected in outcrossing species. The high level of intralocus LD found for carotenoid biosynthesis genes implies that candidate-gene association mapping for carrot root colour should be useful to validate gene function, but may be unable to identify precisely the causative variations involved in trait determinism. Finally this study affords the first molecular evidence of a genetic structure in cultivated carrot germplasm related to phylogeography.

Genes underlying quantitative variation in ecologically important traits: PIF4 (phytochrome interacting factor 4) is associated with variation in internode length, flowering time, and fruit set in Arabidopsis thaliana

Association studies utilize the action of recombination over numerous generations to identify loci that underlie quantitative traits. We use a candidate-gene association approach, segregation analyses and analyses of local linkage disequilibrium (LD) to evaluate the potentially causal effects of molecular variation at PIF4 (PHYTOCROME INTERACTING FACTOR 4) on ecologically important traits in Arabidopsis thaliana. A preliminary analysis of sequence diversity in 14 natural genotypes revealed one intermediate-frequency replacement polymorphism at PIF4. A sample of 161 natural accessions was genotyped at PIF4 and screened for average length of early internodes, inflorescence length, days to flowering and flowering interval (days between bolting and flowering) under high- and low-density environments to test for genotype-phenotype associations. PIF4 was associated with early internode lengths, while the PIF4x treatment interaction was associated with flowering interval in the panel of 161 accessions. Further, in a set of recombinant inbred lines that segregate for the PIF4 polymorphism, nucleotide substitutions at PIF4 co-segregated with early internode lengths, days to flowering and fruit set, suggesting that cryptic population structure in the association-mapping panel and attendant LD with a physically distant locus do not account for the observed association. Finally, in a panel of pseudochromosomes from 20 re-sequenced genotypes, LD appeared to decay rapidly in the immediate vicinity of PIF4, suggesting that flanking loci contribute little to the observed association. In sum, the results suggest that PIF4 causally affects early internode lengths on the primary inflorescence, potentially via effects on reproductive timing and that these traits in turn affect fitness.

Adaptation of tobacco etch potyvirus to a susceptible ecotype of Arabidopsis thaliana capacitates it for systemic infection of resistant ecotypes

Viral pathogens continue to emerge among humans, domesticated animals and cultivated crops. The existence of genetic variance for resistance in the host population is crucial to the spread of an emerging virus. Models predict that rapid spread decreases with the frequency and diversity of resistance alleles in the host population. However, empirical tests of this hypothesis are scarce. Arabiodpsis thaliana--tobacco etch potyvirus (TEV) provides an experimentally suitable pathosystem to explore the interplay between genetic variation in host's susceptibility and virus diversity. Systemic infection of A. thaliana with TEV is controlled by three dominant loci, with different ecotypes varying in susceptibility depending on the genetic constitution at these three loci. Here, we show that the TEV adaptation to a susceptible ecotype allowed the virus to successfully infect, replicate and induce symptoms in ecotypes that were fully resistant to the ancestral virus. The value of these results is twofold. First, we showed that the existence of partially susceptible individuals allows for the emerging virus to bypass resistance alleles that the virus has never encountered. Second, the concept of resistance genes may only be valid for a well-defined viral genotype but not for polymorphic viral populations.

Deleterious amino acid polymorphisms in Arabidopsis thaliana and rice

Plant genetic diversity has been mainly investigated with neutral markers, but large-scale DNA sequencing projects now enable the identification and analysis of different classes of genetic polymorphisms, such as non-synonymous single nucleotide polymorphisms (nsSNPs) in protein coding sequences. Using the SIFT and MAPP programs to predict whether nsSNPs are tolerated (i.e., effectively neutral) or deleterious for protein function, genome-wide nsSNP data from Arabidopsis thaliana and rice were analyzed. In both species, about 20% of polymorphic sites with nsSNPs were classified as deleterious; they segregate at lower allele frequencies than tolerated nsSNPs due to purifying selection. Furthermore, A. thaliana accessions from marginal populations show a higher relative proportion of deleterious nsSNPs, which likely reflects differential selection or demographic effects in subpopulations. To evaluate the sensitivity of predictions, genes from model and crop plants with known functional effects of nsSNPs were inferred with the algorithms. The programs predicted about 70% of nsSNPs correctly as tolerated or deleterious, i.e., as having a functional effect. Forward-in-time simulations of bottleneck and domestication models indicated a high power to detect demographic effects on nsSNP frequencies in sufficiently large datasets. The results indicate that nsSNPs are useful markers for analyzing genetic diversity in plant genetic resources and breeding populations to infer natural/artificial selection and genetic drift.

Genetic variation in tolerance of competition and neighbour suppression in Arabidopsis thaliana

Intraspecific competitive interactions can profoundly influence phenotypic evolution. However, prior studies have rarely evaluated the evolutionary potential of the two components of competitive ability, tolerance of competition and suppression of neighbours. Here, we grow a set of 20 Arabidopsis thaliana recombinant inbred lines in three competitive treatments (noncompetitive, intra-genotypic competition and inter-genotypic competition) to examine if there is genetic variation for the components of competitive ability and whether neighbour relatedness has an effect on fitness. We find evidence for genetic variation in tolerance of competition and neighbour suppression and that these two competitive strategies are correlated, such that genotypes that tolerate competition will also strongly suppress neighbours. We further observe that the effect of neighbour relatedness on fitness of target individuals depends on neighbour identity, i.e. whether target individuals perform better when competing against self vs. nonself individuals depends on the genotypic identity of the nonself neighbour. The results are particularly relevant to evolutionary responses under multi-level selection.

Local-scale patterns of genetic variability, outcrossing, and spatial structure in natural stands of Arabidopsis thaliana

As Arabidopsis thaliana is increasingly employed in evolutionary and ecological studies, it is essential to understand patterns of natural genetic variation and the forces that shape them. Previous work focusing mostly on global and regional scales has demonstrated the importance of historical events such as long-distance migration and colonization. Far less is known about the role of contemporary factors or environmental heterogeneity in generating diversity patterns at local scales. We sampled 1,005 individuals from 77 closely spaced stands in diverse settings around Tübingen, Germany. A set of 436 SNP markers was used to characterize genome-wide patterns of relatedness and recombination. Neighboring genotypes often shared mosaic blocks of alternating marker identity and divergence. We detected recent outcrossing as well as stretches of residual heterozygosity in largely homozygous recombinants. As has been observed for several other selfing species, there was considerable heterogeneity among sites in diversity and outcrossing, with rural stands exhibiting greater diversity and heterozygosity than urban stands. Fine-scale spatial structure was evident as well. Within stands, spatial structure correlated negatively with observed heterozygosity, suggesting that the high homozygosity of natural A. thaliana may be partially attributable to nearest-neighbor mating of related individuals. The large number of markers and extensive local sampling employed here afforded unusual power to characterize local genetic patterns. Contemporary processes such as ongoing outcrossing play an important role in determining distribution of genetic diversity at this scale. Local "outcrossing hotspots" appear to reshuffle genetic information at surprising rates, while other stands contribute comparatively little. Our findings have important implications for sampling and interpreting diversity among A. thaliana accessions.

Single feature polymorphism (SFP)-based selective sweep identification and association mapping of growth-related metabolic traits in Arabidopsis thaliana

Background

Natural accessions of Arabidopsis thaliana are characterized by a high level of phenotypic variation that can be used to investigate the extent and mode of selection on the primary metabolic traits. A collection of 54 A. thaliana natural accession-derived lines were subjected to deep genotyping through Single Feature Polymorphism (SFP) detection via genomic DNA hybridization to Arabidopsis Tiling 1.0 Arrays for the detection of selective sweeps, and identification of associations between sweep regions and growth-related metabolic traits.

Results

A total of 1,072,557 high-quality SFPs were detected and indications for 3,943 deletions and 1,007 duplications were obtained. A significantly lower than expected SFP frequency was observed in protein-, rRNA-, and tRNA-coding regions and in non-repetitive intergenic regions, while pseudogenes, transposons, and non-coding RNA genes are enriched with SFPs. Gene families involved in plant defence or in signalling were identified as highly polymorphic, while several other families including transcription factors are depleted of SFPs. 198 significant associations between metabolic genes and 9 metabolic and growth-related phenotypic traits were detected with annotation hinting at the nature of the relationship. Five significant selective sweep regions were also detected of which one associated significantly with a metabolic trait.

Conclusions

We generated a high density polymorphism map for 54 A. thaliana accessions that highlights the variability of resistance genes across geographic ranges and used it to identify selective sweeps and associations between metabolic genes and metabolic phenotypes. Several associations show a clear biological relationship, while many remain requiring further investigation.

Genome-wide evidence for efficient positive and purifying selection in Capsella grandiflora, a plant species with a large effective population size

Recent studies comparing genome-wide polymorphism and divergence in Drosophila have found evidence for a surprisingly high proportion of adaptive amino acid fixations, but results for other taxa are mixed. In particular, few studies have found convincing evidence for adaptive amino acid substitution in plants. To assess the generality of this finding, we have sequenced 257 loci in the outcrossing crucifer Capsella grandiflora, which has a large effective population size and low population structure. Using a new method that jointly infers selective and demographic effects, we estimate that 40% of amino acid substitutions were fixed by positive selection in this species, and we also infer a low proportion of slightly deleterious amino acid mutations. We contrast these estimates with those for a similar data set from the closely related Arabidopsis thaliana and find significantly higher rates of adaptive evolution and fewer nearly neutral mutations in C. grandiflora. In agreement with results for other taxa, genes involved in reproduction show the strongest evidence for positive selection in C. grandiflora. Taken together, these results imply that both positive and purifying selection are more effective in C. grandiflora than in A. thaliana, consistent with the contrasting demographic history and effective population sizes of these species.

The evolution of transposable elements in natural populations of self-fertilizing Arabidopsis thaliana and its outcrossing relative Arabidopsis lyrata

Background

Transposable Elements (TEs) make up the majority of plant genomes, and thus understanding TE evolutionary dynamics is key to understanding plant genome evolution. Plant reproductive systems are diverse and mating type variation is one factor among many hypothesized to influence TE evolutionary dynamics. Here, we collected a large TE-display data set in self-fertilizing Arabidopsis thaliana, and compared it to data gathered in outcrossing Arabidopsis lyrata. We analyzed seven TE families in four natural populations of each species to tease apart the effects of mating system, demography, transposition, and selection in determining patterns of TE diversity.

Results

Measures of TE band differentiation were largely consistent across TE families. However, patterns of diversity in A. thaliana Ac elements differed significantly from that other TEs, perhaps signaling a lack of recent transposition. Across TE families, we estimated higher allele frequencies and lower selection coefficients on A. thaliana TE insertions relative to A. lyrata TE insertions.

Conclusions

The differences in TE distributions between the two Arabidopsis species represents a synthesis of evolutionary forces that include the transposition dynamics of individual TE families and the demographic histories of populations. There are also species-specific differences that could be attributed to the effects of mating system, including higher overall allele frequencies in the selfing lineage and a greater proportion of among population TE diversity in the outcrossing lineage.

Predicting Arabidopsis freezing tolerance and heterosis in freezing tolerance from metabolite composition

Heterosis, or hybrid vigor, is one of the most important tools in plant breeding and has previously been demonstrated for plant freezing tolerance. Freezing tolerance is an important trait because it can limit the geographical distribution of plants and their agricultural yield. Plants from temperate climates increase in freezing tolerance during exposure to low, non-freezing temperatures in a process termed 'cold acclimation'. Metabolite profiling has indicated a major reprogramming of plant metabolism in the cold, but it has remained unclear in previous studies which of these changes are related to freezing tolerance. In the present study, we have used metabolic profiling to discover combinations of metabolites that predict freezing tolerance and its heterosis in Arabidopsis thaliana. We identified compatible solutes and, in particular, the pathway leading to raffinose as crucial statistical predictors for freezing tolerance and its heterosis, while some TCA cycle intermediates contribute only to predicting the heterotic phenotype. This indicates coordinate links between heterosis and metabolic pathways, suggesting that a limited number of regulatory genes may determine the extent of heterosis in this complex trait. In addition, several unidentified metabolites strongly contributed to the prediction of both freezing tolerance and its heterosis and we present an exemplary analysis of one of these, identifying it as a hexose conjugate.

Hormone- and light-mediated regulation of heat-induced differential petiole growth in Arabidopsis

Plants react quickly and profoundly to changes in their environment. A sudden increase in temperature, for example, induces differential petiole growth-driven upward leaf movement (hyponastic growth) in Arabidopsis (Arabidopsis thaliana). We show that accessions that face the strongest fluctuations in diurnal temperature in their natural habitat are least sensitive for heat-induced hyponastic growth. This indicates that heat-induced hyponastic growth is a trait subject to natural selection. The response is induced with kinetics remarkably similar to ethylene- and low light-induced hyponasty in several accessions. Using pharmacological assays, transcript analysis, and mutant analyses, we demonstrate that ethylene and the photoreceptor protein phytochrome B are negative regulators of heat-induced hyponastic growth and that low light, phytochrome A, auxin, polar auxin transport, and abscisic acid are positive regulators of heat-induced hyponastic growth. Furthermore, auxin, auxin polar transport, phytochrome A, phytochrome B, and cryptochromes are required for a fast induction of heat-induced hyponastic growth.

Demographic and genetic patterns of variation among populations of Arabidopsis thaliana from contrasting native environments

BACKGROUND

Understanding the relationship between environment and genetics requires the integration of knowledge on the demographic behavior of natural populations. However, the demographic performance and genetic composition of Arabidopsis thaliana populations in the species' native environments remain largely uncharacterized. This information, in combination with the advances on the study of gene function, will improve our understanding on the genetic mechanisms underlying adaptive evolution in A. thaliana.

METHODOLOGY/PRINCIPAL FINDINGS

We report the extent of environmental, demographic, and genetic variation among 10 A. thaliana populations from Mediterranean (coastal) and Pyrenean (montane) native environments in northeast Spain. Geographic, climatic, landscape, and soil data were compared. Demographic traits, including the dynamics of the soil seed bank and the attributes of aboveground individuals followed over a complete season, were also analyzed. Genetic data based on genome-wide SNP markers were used to describe genetic diversity, differentiation, and structure. Coastal and montane populations significantly differed in terms of environmental, demographic, and genetic characteristics. Montane populations, at higher altitude and farther from the sea, are exposed to colder winters and prolonged spring moisture compared to coastal populations. Montane populations showed stronger secondary seed dormancy, higher seedling/juvenile mortality in winter, and initiated flowering later than coastal populations. Montane and coastal regions were genetically differentiated, montane populations bearing lower genetic diversity than coastal ones. No significant isolation-by-distance pattern and no shared multilocus genotypes among populations were detected.

CONCLUSIONS/SIGNIFICANCE

Between-region variation in climatic patterns can account for differences in demographic traits, such as secondary seed dormancy, plant mortality, and recruitment, between coastal and montane A. thaliana populations. In addition, differences in plant mortality can partly account for differences in the genetic composition of coastal and montane populations. This study shows how the interplay between variation in environmental, demographic, and genetic parameters may operate in natural A. thaliana populations.

Phytochrome B and histone deacetylase 6 control light-induced chromatin compaction in Arabidopsis thaliana

Natural genetic variation in Arabidopsis thaliana exists for many traits and often reflects acclimation to local environments. Studying natural variation has proven valuable in the characterization of phenotypic traits and, in particular, in identifying genetic factors controlling these traits. It has been previously shown that chromatin compaction changes during development and biotic stress. To gain more insight into the genetic control of chromatin compaction, we investigated the nuclear phenotype of 21 selected Arabidopsis accessions from different geographic origins and habitats. We show natural variation in chromatin compaction and demonstrate a positive correlation with latitude of geographic origin. The level of compaction appeared to be dependent on light intensity. A novel approach, combining Quantitative Trait Locus (QTL) mapping and microscopic examination, pointed at PHYTOCHROME-B (PHYB) and HISTONE DEACETYLASE-6 (HDA6) as positive regulators of light-controlled chromatin compaction. Indeed, mutant analyses demonstrate that both factors affect global chromatin organization. HDA6, in addition, strongly promotes the light-mediated compaction of the Nucleolar Organizing Regions (NORs). The accession Cape Verde Islands-0 (Cvi-0), which shows sequence polymorphism in the PHYB gene and in the HDA6 promotor, resembles the hda6 mutant in having reduced chromatin compaction and decreased methylation levels of DNA and histone H3K9 at the NORs. We provide evidence that chromatin organization is controlled by light intensity. We propose that chromatin plasticity is associated with acclimation of Arabidopsis to its environment. The polymorphic alleles such as PHYB and HDA6 control this process.

The habitat of the plant model species Arabidopsis thaliana can be found throughout the Northern hemisphere. As a consequence, individual populations have acclimated to a great diversity of environmental conditions. This is reflected by a wealth of natural genetic variation in many phenotypic traits. We utilized this natural variation via a novel approach, combining microscopic examination, quantitative genetics, and analysis of environmental parameters, to understand the regulation of nuclear chromatin compaction in leaf mesophyll cells. We show that the level of chromatin compaction among natural Arabidopsis thaliana accessions correlates with latitude of origin and depends on local light intensity. Our study provides evidence that the photoreceptor PHYTOCHROME-B (PHYB) and the histone modifier HISTONE DEACETYLASE 6 (HDA6) are positive regulators of global chromatin organization in a light-dependent manner. In addition, HDA6 specifically controls light-mediated chromatin compaction of the Nucleolar Organizing Regions (NORs). We propose that the observed light-controlled plasticity of chromatin plays a role in acclimation and survival of plants in their natural environment.

A focus on natural variation for abiotic constraints response in the model species Arabidopsis thaliana

Plants are particularly subject to environmental stress, as they cannot move from unfavourable surroundings. As a consequence they have to react in situ. In any case, plants have to sense the stress, then the signal has to be transduced to engage the appropriate response. Stress response is effected by regulating genes, by turning on molecular mechanisms to protect the whole organism and its components and/or to repair damage. Reactions vary depending on the type of stress and its intensity, but some are commonly turned on because some responses to different abiotic stresses are shared. In addition, there are multiple ways for plants to respond to environmental stress, depending on the species and life strategy, but also multiple ways within a species depending on plant variety or ecotype. It is regularly accepted that populations of a single species originating from diverse geographic origins and/or that have been subjected to different selective pressure, have evolved retaining the best alleles for completing their life cycle. Therefore, the study of natural variation in response to abiotic stress, can help unravel key genes and alleles for plants to cope with their unfavourable physical and chemical surroundings. This review is focusing on Arabidopsis thaliana which has been largely adopted by the global scientific community as a model organism. Also, tools and data that facilitate investigation of natural variation and abiotic stress encountered in the wild are set out. Characterization of accessions, QTLs detection and cloning of alleles responsible for variation are presented.