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

The evolution of selfing from outcrossing ancestors in Brassicaceae: what have we learned from variation at the S-locus?

Evolutionary transitions between mating systems have occurred repetitively and independently in flowering plants. One of the most spectacular advances of the recent empirical literature in the field was the discovery of the underlying genetic machinery, which provides the opportunity to retrospectively document the scenario of the outcrossing to selfing transitions in a phylogenetic perspective. In this review, we explore the literature describing patterns of polymorphism and molecular evolution of the locus controlling self-incompatibility (S-locus) in selfing species of the Brassicaceae family in order to document the transition from outcrossing to selfing, a retrospective approach that we describe as the 'mating system genes approach'. The data point to strikingly contrasted scenarios of transition from outcrossing to selfing. We also perform original analyses of the fully sequenced genomes of four species showing self-compatibility, to compare the orthologous S-locus region with that of functional S-locus haplotypes. Phylogenetic analyses suggest that all species we investigated evolved independently towards loss of self-incompatibility, and in most cases almost intact sequences of either of the two S-locus genes suggest that these transitions occurred relatively recently. The S-locus region in Aethionema arabicum, representing the most basal lineage of Brassicaceae, showed unusual patterns so that our analysis could not determine whether self-incompatibility was lost secondarily, or evolved in the core Brassicaceae after the split with this basal lineage. Although the approach we detail can only be used when mating system genes have been identified in a clade, we suggest that its integration with phylogenetic and population genetic approaches should help determine the main routes of this predominant mating system shift in plants.

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.

Genetic mapping of adaptation reveals fitness tradeoffs in Arabidopsis thaliana

Organisms inhabiting different environments are often locally adapted, and yet despite a considerable body of theory, the genetic basis of local adaptation is poorly understood. Unanswered questions include the number and effect sizes of adaptive loci, whether locally favored loci reduce fitness elsewhere (i.e., fitness tradeoffs), and whether a lack of genetic variation limits adaptation. To address these questions, we mapped quantitative trait loci (QTL) for total fitness in 398 recombinant inbred lines derived from a cross between locally adapted populations of the highly selfing plant Arabidopsis thaliana from Sweden and Italy and grown for 3 consecutive years at the parental sites (>40,000 plants monitored). We show that local adaptation is controlled by relatively few genomic regions of small to modest effect. A third of the 15 fitness QTL we detected showed evidence of tradeoffs, which contrasts with the minimal evidence for fitness tradeoffs found in previous studies. This difference may reflect the power of our multiyear study to distinguish conditionally neutral QTL from those that reflect fitness tradeoffs. In Sweden, but not in Italy, the local genotype underlying fitness QTL was often maladaptive, suggesting that adaptation there is constrained by a lack of adaptive genetic variation, attributable perhaps to genetic bottlenecks during postglacial colonization of Scandinavia or to recent changes in selection regime caused by climate change. Our results suggest that adaptation to markedly different environments can be achieved through changes in relatively few genomic regions, that fitness tradeoffs are common, and that lack of genetic variation can limit adaptation.

Parental genetic distance and patterns in nonrandom mating and seed yield in predominately selfing Arabidopsis thaliana

In this study, we ask two questions: (1) Is reproductive success independent of parental genetic distance in predominately selfing plants? (2) In the absence of early inbreeding depression, is there substantial maternal and/or paternal variation in reproductive success in natural populations? Seed yield in single pollinations and proportion of seeds sired in mixed pollinations were studied in genetically defined accessions of the predominately selfing plant Arabidopsis thaliana by conducting two diallel crosses. The first diallel was a standard, single pollination design that we used to examine variance in seed yield. The second diallel was a mixed pollination design that utilized a standard pollen competitor to examine variance in proportion of seeds sired. We found no correlation between reproductive success and parental genetic distance, and self-pollen does not systematically differ in reproductive success compared to outcross pollen, suggesting that Arabidopsis populations do not experience embryo lethality due to early-acting inbreeding or outbreeding depression. We used these data to partition the contributions to total phenotypic variation from six sources, including maternal contributions, paternal contributions and parental interactions. For seed yield in single pollinations, maternal effects accounted for the most significant source of variance (16.6 %). For proportion of seeds sired in mixed pollinations, the most significant source of variance was paternal effects (17.9 %). Thus, we show that population-level genetic similarities, including selfing, do not correlate with reproductive success, yet there is still significant paternal variance under competition. This suggests two things. First, since these differences are unlikely due to early-acting inbreeding depression or differential pollen viability, this implicates natural variation in pollen germination and tube growth dynamics. Second, this strongly supports a model of fixation of pollen performance genes in populations, offering a focus for future genetic studies in differential reproductive success.

Phylogeography of mugwort (Artemisia indica), a native pioneer herb in Japan

Many phylogeographic studies of various tree species have been conducted to elucidate the locations of refugia and the colonization patterns during the Pleistocene. However, only a few large-scale phylogeographic studies have been conducted on herbaceous plants, especially scarce on herbs that are adapted to disturbance. Artemisia indica is a fast-growing perennial herb found in open habitats. To examine the basic information on the genetic structure of this species, we investigated the chloroplast DNA variation within and among populations across Japan. We detected 26 haplotypes in 604 individuals from 28 Japanese populations. The haplotype A1 had wide geographic distribution, and its close relatives were locally present. Some putative ancestral lineages were found mainly in the Kyushu region. This may be because several lineages migrated from the Eurasian continent to the northern coast in Kyushu via the Korean peninsula during the Pleistocene, and the A1 haplotype expanded northward, whereas others remained in southern areas. Phylogenetic distant haplotypes were present mainly in the Kanto region. Because the geographic distribution pattern of these haplotypes in this region is believed to be unnatural, these haplotypes may be derived from commercial sources for re-vegetation during the last few decades.

Evidence for parallel adaptation to climate across the natural range of Arabidopsis thaliana

How organisms adapt to different climate habitats is a key question in evolutionary ecology and biological conservation. Species distributions are often determined by climate suitability. Consequently, the anthropogenic impact on earth's climate is of key concern to conservation efforts because of our relatively poor understanding of the ability of populations to track and evolve to climate change. Here, we investigate the ability of Arabidopsis thaliana to occupy climate space by quantifying the extent to which different climate regimes are accessible to different A. thaliana genotypes using publicly available data from a large-scale genotyping project and from a worldwide climate database. The genetic distance calculated from 149 single-nucleotide polymorphisms (SNPs) among 60 lineages of A. thaliana was compared to the corresponding climate distance among collection localities calculated from nine different climatic factors. A. thaliana was found to be highly labile when adapting to novel climate space, suggesting that populations may experience few constraints when adapting to changing climates. Our results also provide evidence of a parallel or convergent evolution on the molecular level supporting recent generalizations regarding the genetics of adaptation.

Co-variation between seed dormancy, growth rate and flowering time changes with latitude in Arabidopsis thaliana

Life-history traits controlling the duration and timing of developmental phases in the life cycle jointly determine fitness. Therefore, life-history traits studied in isolation provide an incomplete view on the relevance of life-cycle variation for adaptation. In this study, we examine genetic variation in traits covering the major life history events of the annual species Arabidopsis thaliana: seed dormancy, vegetative growth rate and flowering time. In a sample of 112 genotypes collected throughout the European range of the species, both seed dormancy and flowering time follow a latitudinal gradient independent of the major population structure gradient. This finding confirms previous studies reporting the adaptive evolution of these two traits. Here, however, we further analyze patterns of co-variation among traits. We observe that co-variation between primary dormancy, vegetative growth rate and flowering time also follows a latitudinal cline. At higher latitudes, vegetative growth rate is positively correlated with primary dormancy and negatively with flowering time. In the South, this trend disappears. Patterns of trait co-variation change, presumably because major environmental gradients shift with latitude. This pattern appears unrelated to population structure, suggesting that changes in the coordinated evolution of major life history traits is adaptive. Our data suggest that A. thaliana provides a good model for the evolution of trade-offs and their genetic basis.

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.

Evolution of flux control in the glucosinolate pathway in Arabidopsis thaliana

Network characteristics of biochemical pathways are believed to influence the rate of evolutionary change in constituent enzymes. One characteristic that may affect rate heterogeneity is control of the amount of product produced by a biochemical pathway or flux control. In particular, theoretical analyses suggest that adaptive substitutions should be concentrated in the enzyme(s) that exert the greatest control over flux. Although a handful of studies have found a correlation between position in a pathway and evolutionary rate, these investigations have not examined the relationship between evolutionary rate and flux control. Given that genes with greater control will experience stronger selection and that the probability of fixation is proportional to the selective advantage, we ask the following: 1) do upstream enzymes have majority flux control, 2) do enzymes with majority flux control accumulate adaptive substitutions, and 3) are upstream enzymes under higher selective constraint? First, by perturbing the enzymes in the aliphatic glucosinolate pathway in Arabidopsis thaliana with gene insertion lines, we show that flux control is focused in the first enzyme in the pathway. Next, by analyzing several sequence signatures of selection, we also show that this enzyme is the only one in the pathway that shows convincing evidence of selection. Our results support the hypothesis that natural selection preferentially acts on enzymes with high flux control.

Genome-wide association studies identify heavy metal ATPase3 as the primary determinant of natural variation in leaf cadmium in Arabidopsis thaliana

Understanding the mechanism of cadmium (Cd) accumulation in plants is important to help reduce its potential toxicity to both plants and humans through dietary and environmental exposure. Here, we report on a study to uncover the genetic basis underlying natural variation in Cd accumulation in a world-wide collection of 349 wild collected Arabidopsis thaliana accessions. We identified a 4-fold variation (0.5–2 µg Cd g⁻¹ dry weight) in leaf Cd accumulation when these accessions were grown in a controlled common garden. By combining genome-wide association mapping, linkage mapping in an experimental F2 population, and transgenic complementation, we reveal that HMA3 is the sole major locus responsible for the variation in leaf Cd accumulation we observe in this diverse population of A. thaliana accessions. Analysis of the predicted amino acid sequence of HMA3 from 149 A. thaliana accessions reveals the existence of 10 major natural protein haplotypes. Association of these haplotypes with leaf Cd accumulation and genetics complementation experiments indicate that 5 of these haplotypes are active and 5 are inactive, and that elevated leaf Cd accumulation is associated with the reduced function of HMA3 caused by a nonsense mutation and polymorphisms that change two specific amino acids.

Cadmium (Cd) is a potentially toxic metal pollutant that threatens food quality and human health in many regions of the world. Plants have evolved mechanisms for the acquisition of essential metals such as zinc and iron from the soil. Though often quite specific, such mechanisms can also lead to the accumulation of Cd by plants. Understanding natural variation in the processes that contribute to Cd accumulation in food crops could help minimize the human health risk posed. We have discovered that DNA sequence changes at a single gene, which encodes the Heavy Metal ATPase 3 (HMA3), drives the variation in Cd accumulation we observe in a world-wide sample of Arabidopsis thaliana. We identified 10 major HMA3 protein variants, of which five contribute to reduce Cd accumulation in leaves of A. thaliana.

The evolutionary history of the Arabidopsis arenosa complex: diverse tetraploids mask the Western Carpathian center of species and genetic diversity

The Arabidopsis arenosa complex is closely related to the model plant Arabidopsis thaliana. Species and subspecies in the complex are mainly biennial, predominantly outcrossing, herbaceous, and with a distribution range covering most parts of latitudes and the eastern reaches of Europe. In this study we present the first comprehensive evolutionary history of the A. arenosa species complex, covering its natural range, by using chromosome counts, nuclear AFLP data, and a maternally inherited marker from the chloroplast genome [trnL intron (trnL) and trnL/F intergenic spacer (trnL/F-IGS) of tRNA(Leu) and tRNA(Phe), respectively]. We unravel the broad-scale cytogeographic and phylogeographic patterns of diploids and tetraploids. Diploid cytotypes were exclusively found on the Balkan Peninsula and in the Carpathians while tetraploid cytotypes were found throughout the remaining distribution range of the A. arenosa complex. Three centers of genetic diversity were identified: the Balkan Peninsula, the Carpathians, and the unglaciated Eastern and Southeastern Alps. All three could have served as long-term refugia during Pleistocene climate oscillations. We hypothesize that the Western Carpathians were and still are the cradle of speciation within the A. arenosa complex due to the high species number and genetic diversity and the concurrence of both cytotypes there.

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.

Recent loss of self-incompatibility by degradation of the male component in allotetraploid Arabidopsis kamchatica

The evolutionary transition from outcrossing to self-fertilization (selfing) through the loss of self-incompatibility (SI) is one of the most prevalent events in flowering plants, and its genetic basis has been a major focus in evolutionary biology. In the Brassicaceae, the SI system consists of male and female specificity genes at the S-locus and of genes involved in the female downstream signaling pathway. During recent decades, much attention has been paid in particular to clarifying the genes responsible for the loss of SI. Here, we investigated the pattern of polymorphism and functionality of the female specificity gene, the S-locus receptor kinase (SRK), in allotetraploid Arabidopsis kamchatica. While its parental species, A. lyrata and A. halleri, are reported to be diploid and mainly self-incompatible, A. kamchatica is self-compatible. We identified five highly diverged SRK haplogroups, found their disomic inheritance and, for the first time in a wild allotetraploid species, surveyed the geographic distribution of SRK at the two homeologous S-loci across the species range. We found intact full-length SRK sequences in many accessions. Through interspecific crosses with the self-incompatible and diploid congener A. halleri, we found that the female components of the SI system, including SRK and the female downstream signaling pathway, are still functional in these accessions. Given the tight linkage and very rare recombination of the male and female components on the S-locus, this result suggests that the degradation of male components was responsible for the loss of SI in A. kamchatica. Recent extensive studies in multiple Brassicaceae species demonstrate that the loss of SI is often derived from mutations in the male component in wild populations, in contrast to cultivated populations. This is consistent with theoretical predictions that mutations disabling male specificity are expected to be more strongly selected than mutations disabling female specificity, or the female downstream signaling pathway.

Allelic heterogeneity and trade-off shape natural variation for response to soil micronutrient

As sessile organisms, plants have to cope with diverse environmental constraints that may vary through time and space, eventually leading to changes in the phenotype of populations through fixation of adaptive genetic variation. To fully comprehend the mechanisms of evolution and make sense of the extensive genotypic diversity currently revealed by new sequencing technologies, we are challenged with identifying the molecular basis of such adaptive variation. Here, we have identified a new variant of a molybdenum (Mo) transporter, MOT1, which is causal for fitness changes under artificial conditions of both Mo-deficiency and Mo-toxicity and in which allelic variation among West-Asian populations is strictly correlated with the concentration of available Mo in native soils. In addition, this association is accompanied at different scales with patterns of polymorphisms that are not consistent with neutral evolution and show signs of diversifying selection. Resolving such a case of allelic heterogeneity helps explain species-wide phenotypic variation for Mo homeostasis and potentially reveals trade-off effects, a finding still rarely linked to fitness.

Plants are studied for their ability to adapt to their environment and especially to the physical constraints to which they are subjected. It is expected that they evolve in promoting genetic variants favorable under their native conditions, which could lead to negative consequences in other conditions. One approach to study the mechanisms and dynamics of these adaptations is to discover genetic variants that control potentially adaptive traits, and to study directly these variants in wild populations to try to reveal their evolutionary trajectory. We have identified a new polymorphism in a gene coding for a transporter of molybdenum (an essential micronutrient for the plant) in Arabidopsis; we show that this variant has strong phenotypic consequences at the level of plant growth and reproductive value in specific conditions, and that it explains a lot of the species diversity for these traits. Especially, the variant is associated with a clear negative effect under molybdenum-deficient conditions (caused by soil acidity) and with a subtle positive effect under molybdenum-plethoric conditions. Interestingly, the landscape distribution of the variant is not random among Asian populations and correlates well with the availability of molybdenum in the soil at the precise location where the plants are growing in the wild.

Reciprocal transplants demonstrate strong adaptive differentiation of the model organism Arabidopsis thaliana in its native range

To quantify adaptive differentiation in the model plant Arabidopsis thaliana, we conducted reciprocal transplant experiments for five years between two European populations, one near the northern edge of the native range (Sweden) and one near the southern edge (Italy). We planted seeds (years 1-3) and seedlings (years 4-5), and estimated fitness as the number of fruits produced per seed or seedling planted. In eight of the 10 possible site × year comparisons, the fitness of the local population was significantly higher than that of the nonlocal population (3.1-22.2 times higher at the southern site, and 1.7-3.6 times higher at the northern site); in the remaining two comparisons no significant difference was recorded. At both sites, the local genotype had higher survival than the nonlocal genotype, and at the Italian site, the local genotype also had higher fecundity. Across years, the relative survival of the Italian genotype at the northern site decreased with decreasing winter soil temperature. The results provide evidence of strong adaptive differentiation between natural populations of A. thaliana and indicate that differences in tolerance to freezing contributed to fitness variation at the northern site. In ongoing work, we explore the functional and genetic basis of this adaptive differentiation.

Natural variation in GL1 and its effects on trichome density in Arabidopsis thaliana

The ultimate understanding of how biological diversity arises, is maintained, and lost depends on identifying the genes responsible. Although a good deal has been discovered about gene function over the past few decades, far less is understood about gene effects, that is, how natural variation in a gene contributes to natural variation in phenotypes. Trichome density in Arabidopsis thaliana is an ideal trait for studies of natural molecular and phenotypic variation, as trichome initiation is genetically well-characterized and trichome density is highly variable in and among natural populations. Here, we show that variation at GLABRA1 (GL1), an R2R3 MYB transcription factor gene, which has a role in trichome initiation, has qualitative and likely quantitative effects on trichome density in natural accessions of A. thaliana. Specifically, we characterize four independent loss-of-function alleles for GL1, each of which yields a glabrous phenotype. Further, we find that a pattern of common polymorphisms confined to the GL1 locus is associated with quantitative variation for trichome density. While mutations resulting in a glabrous phenotype are primarily coding changes, the pattern resulting in quantitative variation spans both coding and regulatory regions. These data show that GL1 is an important source of trichome density variation within A. thaliana and, along with recent reports, suggest that the TTG1 epidermal cell fate pathway generally may be the key molecular genetic source of natural trichome density variation and an important model for the study of molecular evolution.

Intricate variation patterns in the diploid-polyploid complex of Alyssum montanum-A. repens (Brassicaceae) in the Apennine Peninsula: evidence for long-term persistence and diversification

PREMISE OF THE STUDY

The Apennine Peninsula, one of the three peninsulas of the European Mediterranean, is an important hotspot of genetic and species diversity, but studies devoted to plant evolution are still very scarce in this region. We studied the diploid-polyploid complex of Alyssum montanum-A. repens, focusing on Apennine and adjacent southwestern Alpine populations from southeastern France, with the aim of examining their taxonomic position and evolutionary patterns.

METHODS

We employed AFLP markers and cpDNA sequences, along with cytotype determination using flow cytometry, and a multivariate morphometric approach.

KEY RESULTS

The Italian and French populations formed two well-delimited groupings within the studied complex that were, in contrast to previous taxonomic treatments, clearly distinct from A. montanum. Populations from southeastern France represent A. orophilum, a previously described but abandoned species. Those from central and southern Italy correspond to A. diffusum, exhibiting high, geographically structured variation (central Apennines, Gargano, and southern Apennines/Calabria). This pattern coincides with hotspot refugial regions, in congruence with the "refugia-within-refugia" hypothesis, and is reflected here in the recognition of three subspecies within A. diffusum.

CONCLUSIONS

We provide evidence for the presence of Mediterranean refugia for the studied Alyssum montanum-A. repens complex located in central and southern Italy, which, however, did not contribute to the postglacial colonization of Central Europe. Past extinctions, genetic bottlenecks, and recent expansion were inferred in Central Europe, while long-term accumulation of diversity as well as polyploidization occurred in the Apennines.

Double-strand break repair processes drive evolution of the mitochondrial genome in Arabidopsis

Background

The mitochondrial genome of higher plants is unusually dynamic, with recombination and nonhomologous end-joining (NHEJ) activities producing variability in size and organization. Plant mitochondrial DNA also generally displays much lower nucleotide substitution rates than mammalian or yeast systems. Arabidopsis displays these features and expedites characterization of the mitochondrial recombination surveillance gene MSH1 (MutS 1 homolog), lending itself to detailed study of de novo mitochondrial genome activity. In the present study, we investigated the underlying basis for unusual plant features as they contribute to rapid mitochondrial genome evolution.

Results

We obtained evidence of double-strand break (DSB) repair, including NHEJ, sequence deletions and mitochondrial asymmetric recombination activity in Arabidopsis wild-type and msh1 mutants on the basis of data generated by Illumina deep sequencing and confirmed by DNA gel blot analysis. On a larger scale, with mitochondrial comparisons across 72 Arabidopsis ecotypes, similar evidence of DSB repair activity differentiated ecotypes. Forty-seven repeat pairs were active in DNA exchange in the msh1 mutant. Recombination sites showed asymmetrical DNA exchange within lengths of 50- to 556-bp sharing sequence identity as low as 85%. De novo asymmetrical recombination involved heteroduplex formation, gene conversion and mismatch repair activities. Substoichiometric shifting by asymmetrical exchange created the appearance of rapid sequence gain and loss in association with particular repeat classes.

Conclusions

Extensive mitochondrial genomic variation within a single plant species derives largely from DSB activity and its repair. Observed gene conversion and mismatch repair activity contribute to the low nucleotide substitution rates seen in these genomes. On a phenotypic level, these patterns of rearrangement likely contribute to the reproductive versatility of higher plants.

Population genetic diversity and structure of a naturally isolated plant species, Rhodiola dumulosa (Crassulaceae)

Aims

Rhodiola dumulosa (Crassulaceae) is a perennial diploid species found in high-montane areas. It is distributed in fragmented populations across northern, central and northwestern China. In this study, we aimed to (i) measure the genetic diversity of this species and that of its populations; (ii) describe the genetic structure of these populations across the entire distribution range in China; and (iii) evaluate the extent of gene flow among the naturally fragmented populations.

Methods

Samples from 1089 individuals within 35 populations of R. dumulosa were collected, covering as much of the entire distribution range of this species within China as possible. Population genetic diversity and structure were analyzed using AFLP molecular markers. Gene flow among populations was estimated according to the level of population differentiation.

Important Findings

The total genetic diversity of R. dumulosa was high but decreased with increasing altitude. Population-structure analysis indicated that the most closely related populations were geographically restricted and occurred in close proximity to each other. A significant isolation-by-distance pattern, caused by the naturally fragmented population distribution, was observed. At least two distinct gene pools were found in the 35 sampled populations, one composed of populations in northern China and the other composed of populations in central and northwestern China. The calculation of Nei's gene diversity index revealed that the genetic diversity in the northern China pool (0.1972) was lower than that in the central and northwestern China pool (0.2216). The populations were significantly isolated, and gene flow was restricted throughout the entire distribution. However, gene flow among populations on the same mountain appears to be unrestricted, as indicated by the weak genetic isolation among these populations.

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.

Reduced efficacy of natural selection on codon usage bias in selfing Arabidopsis and Capsella species

Population genetic theory predicts that the efficacy of natural selection in a self-fertilizing species should be lower than its outcrossing relatives because of the reduction in the effective population size (N(e)) in the former brought about by inbreeding. However, previous analyses comparing Arabidopsis thaliana (selfer) with A. lyrata (outcrosser) have not found conclusive support for this prediction. In this study, we addressed this issue by examining silent site polymorphisms (synonymous and intronic), which are expected to be informative about changes in N(e). Two comparisons were made: A. thaliana versus A. lyrata and Capsella rubella (selfer) versus C. grandiflora (outcrosser). Extensive polymorphism data sets were obtained by compiling published data from the literature and by sequencing 354 exon loci in C. rubella and 89 additional loci in C. grandiflora. To extract information from the data effectively for studying these questions, we extended two recently developed models in order to investigate detailed selective differences between synonymous codons, mutational biases, and biased gene conversion (BGC), taking into account the effects of recent changes in population size. We found evidence that selection on synonymous codons is significantly weaker in the selfers compared with the outcrossers and that this difference cannot be fully accounted for by mutational biases or BGC.

Source verification of mis-identified Arabidopsis thaliana accessions

A major strength of Arabidopsis thaliana as a model lies in the availability of a large number of naturally occurring inbred lines. Recent studies of A. thaliana population structure, using thousands of accessions from stock center and natural collections, have revealed a robust pattern of isolation by distance at several spatial scales, such that genetically identical individuals are generally found close to each other. However, some individual accessions deviate from this pattern. While some of these may be the products of rare long-distance dispersal events, many deviations may be the result of mis-identification, in the sense that the data regarding location of origin data are incorrect. Here, we aim to identify such discrepancies. Of the 5965 accessions examined, we conclude that 286 deserve special attention as being potentially mis-identified. We describe these suspicious accessions and their possible origins, and advise caution with regard to their use in experiments in which accurate information on geographic origin is important. Finally, we discuss possibilities for maintaining the integrity of stock lines.

Effective population size is positively correlated with levels of adaptive divergence among annual sunflowers

The role of adaptation in the divergence of lineages has long been a central question in evolutionary biology, and as multilocus sequence data sets have become available for a wide range of taxa, empirical estimates of levels of adaptive molecular evolution are increasingly common. Estimates vary widely among taxa, with high levels of adaptive evolution in Drosophila, bacteria, and viruses but very little evidence of widespread adaptive evolution in hominids. Although estimates in plants are more limited, some recent work has suggested that rates of adaptive evolution in a range of plant taxa are surprisingly low and that there is little association between adaptive evolution and effective population size in contrast to patterns seen in other taxa. Here, we analyze data from 35 loci for six sunflower species that vary dramatically in effective population size. We find that rates of adaptive evolution are positively correlated with effective population size in these species, with a significant fraction of amino acid substitutions driven by positive selection in the species with the largest effective population sizes but little or no evidence of adaptive evolution in species with smaller effective population sizes. Although other factors likely contribute as well, in sunflowers effective population size appears to be an important determinant of rates of adaptive evolution.

Natural variation in Arabidopsis thaliana revealed a genetic network controlling germination under salt stress

Plant responses to environmental stresses are polygenic and complex traits. In this study quantitative genetics using natural variation in Arabidopsis thaliana was used to investigate the genetic architecture of plant responses to salt stress. Eighty seven A. thaliana accessions were screened and showed a large variation for root development and seed germination under 125 and 200 mM NaCl, respectively. Twenty two quantitative trait loci for these traits have been detected by phenotyping two recombinants inbred line populations, Sha x Col and Sha x Ler. Four QTLs controlling germination under salt were detected in the Sha x Col population. Interestingly, only one allelic combination at these four QTLs inhibits germination under salt stress, implying strong epistatic interactions between them. In this interacting context, we confirmed the effect of one QTL by phenotyping selected heterozygous inbred families. We also showed that this QTL is involved in the control of germination under other stress conditions such as KCl, mannitol, cold, glucose and ABA. Our data highlights the presence of a genetic network which consists of four interacting QTLs and controls germination under limiting environmental conditions.

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.

Cytoplasmic phylogeny and evidence of cyto-nuclear co-adaptation in Arabidopsis thaliana

In recent years Arabidopsis thaliana has become a model species for genomic variability and adaptation studies. Although impressive quantities of data have been gathered on the nuclear genomic diversity of this species, little has been published regarding its cytoplasmic diversity. We analyzed the diversity of plastid (pt) and mitochondrial (mt) genomes among 95 accessions, covering most Arabidopsis geographic origins. Four intergenic regions of the pt genome were sequenced, and a total of 68 polymorphisms and 65 pt haplotypes were identified. Several strategies were developed to identify mt polymorphisms among a subset of 14 accessions. Fifteen polymorphisms were further developed as PCR-based markers and used to analyze the whole set of 95 accessions. Using statistical parsimony, we built pt and mt phylogenetic networks of haplotype groups. To root the pt network, the pt intergenic regions of two related Arabidopsis species, Arabidopsis lyrata and Arabidopsis arenosa, were also sequenced. The mt and pt phylogenies are highly congruent and could be combined into a single cytoplasmic phylogeny. To estimate whether co-adaptation between nuclear and cytoplasmic genomes exists in A. thaliana, we tested the germination capacity in challenging conditions of 27 pairs of reciprocal F(2) families. We found that the cytoplasm donor had a significant effect on the germination capacity of some F(2) families.

Natural diversity in flowering responses of Arabidopsis thaliana caused by variation in a tandem gene array

Tandemly arrayed genes that belong to gene families characterize genomes of many organisms. Gene duplication and subsequent relaxation of selection can lead to the establishment of paralogous cluster members that may evolve along different trajectories. Here, we report on the structural variation in MADS AFFECTING FLOWERING 2 (MAF2) gene, one member of the tandemly duplicated cluster of MADS-box-containing transcription factors in Arabidopsis thaliana. The altered gene structure at the MAF2 locus is present as a moderate-frequency polymorphism in Arabidopsis and leads to the extensive diversity in transcript patterns due to alternative splicing. Rearrangements at the MAF2 locus are associated with an early flowering phenotype in BC(5) lines. The lack of suppression of flowering time in a MAF2-insertion line expressing the MAF2-specific artificial miRNA suggests that these MAF2 variants are behaving as loss-of-function alleles. The variation in gene architecture is also associated with segregation distortion, which may have facilitated the spread and the establishment of the corresponding alleles throughout the Eurasian range of the A. thaliana population.

Genetic structure and phylogeography of a temperate-boreal herb, Cardamine scutata (Brassicaceae), in northeast Asia inferred from AFLPs and cpDNA haplotypes

PREMISE OF THE STUDY

Studies on genetic structure of plant populations help us understand the history of local flora and vegetation. In this study, we focus on the temperate-boreal herb Cardamine scutata from northeast Asia, an area with scarce phylogeographic studies. We explore patterns of genetic variation within this species, with an aim to infer its (post-) glacial history with reference to colonization routes and migrations via land bridges. •

METHODS

We analyzed 46 populations sampled in Japan, Kamchatka, and Korea using AFLP and cpDNA sequence data. •

KEY RESULTS

Two intraspecific genetic groups were resolved, distributed in the northeastern and southwestern part of the study area, most likely reflecting lineages isolated from each other during (at least) the last glaciation. A zone of secondary contacts was found in central/northern Honshu, and a few cases of long-distance dispersal were observed. We detected efficient gene flow across the marine straits, supporting the role of land bridges created by sea level decline during the last glacial period. The cpDNA data indicated extensive recent expansion and diversification within both lineages. We inferred recent colonization of Kamchatka from Hokkaido, associated with genetic impoverishment. •

CONCLUSIONS

The pattern of north-south genetic differentiation found in C. scutata is rather common among several other plant species studied in Japan, despite their distinct biological features. We assume that different processes and factors may have brought about this similarity. Overall, this study contributes to better understanding of the biogeography of northeast Asia.

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.

Influence of mutation rate on estimators of genetic differentiation--lessons from Arabidopsis thaliana

BACKGROUND

The analysis of molecular variation within and between populations is crucial to establish strategies for conservation as well as to detect the footprint of spatially heterogeneous selection. The traditional estimator of genetic differentiation (F(ST)) has been shown to be misleading if genetic diversity is high. Alternative estimators of F(ST) have been proposed, but their robustness to variation in mutation rate is not clearly established. We first investigated the effect of mutation and migration rate using computer simulations and examined their joint influence on Q(ST), a measure of genetic differentiation for quantitative traits. We further used experimental data in natural populations of Arabidopsis thaliana to characterize the effect of mutation rate on various estimates of population differentiation. Since natural species exhibit various degrees of self-fertilisation, we also investigated the effect of mating system on the different estimators.

RESULTS

If mutation rate is high and migration rate low, classical measures of genetic differentiation are misleading. Only Phi(ST), an estimator that takes the mutational distances between alleles into account, is independent of mutation rate, for all migration rates. However, the performance of Phi(ST) depends on the underlying mutation model and departures from this model cause its performance to degrade. We further show that Q(ST) has the same bias. We provide evidence that, in A. thaliana, microsatellite variation correlates with mutation rate. We thereby demonstrate that our results on estimators of genetic differentiation have important implications, even for species that are well established models in population genetics and molecular biology.

CONCLUSIONS

We find that alternative measures of differentiation like F'(ST) and D are not suitable for estimating effective migration rate and should not be used in studies of local adaptation. Genetic differentiation should instead be measured using an estimator that takes mutation rate into account, such as Phi(ST). Furthermore, in systems where migration between populations is low, such as A. thaliana, Q(ST) < F(ST) cannot be taken as evidence for homogenising selection as has been traditionally thought.

Historical spatial range expansion and a very recent bottleneck of Cinnamomum kanehirae Hay. (Lauraceae) in Taiwan inferred from nuclear genes

Background

Species in the varied geographic topology of Taiwan underwent obvious demographic changes during glacial periods. Cinnamomum kanehirae has been exploited for timber and to obtain medicinal fungi for the past 100 years. Understanding anthropogenic factors influencing the demography of this species after the last glacial maximum (LGM) is critically important for the conservation of this species.

Results

Populations of C. kanehirae were classified into four geographic regions: northwestern (NW), west-central (WC), southwestern (SW), and southeastern (SE). In total, 113 individuals from 19 localities were sampled, and variations in the chalcone synthase gene (Chs) intron and leafy (Lfy) intron-2 sequences of nuclear DNA were examined in order to assess phylogeographic patterns, the timescales of demographic and evolutionary events, and recent anthropogenic effects. In total, 210 Chs and 170 Lfy sequences, which respectively constituted 36 and 35 haplotypes, were used for the analyses. Estimates of the migration rate (M) through time revealed a pattern of frequent gene flow during previous and the present interglacials. The isolation-by-distance test showed that there generally was no significant correlation between genetic and geographic distances. The level of among-region genetic differentiation was significant when comparing eastern to western populations. However, no significant among-region genetic differentiation was found in comparisons among the four geographic regions. Moreover, essentially no genetic structuring was found for the three regions west of the CMR. A fit of spatial range expansion was found for pooled and regional samples according to the non-significant values of the sum of squared deviations. Using the Bayesian skyline plot (BSP) method, a recent bottleneck after the LGM expansion was detected in both regional and pooled samples.

Conclusions

Common haplotype distributions among geographic regions and the relatively shallow genetic structuring displayed are the result of historical gene flows. Southward dispersals in an earlier time frame from the NW region and in a later time frame from the SE region were inferred. The BSP analysis suggested a postglacial expansion event. Recent trends, however, refer to a bottleneck due to human interventions observed for both pooled and regional C. kanehirae samples.

The evolutionary history of the Arabidopsis lyrata complex: a hybrid in the amphi-Beringian area closes a large distribution gap and builds up a genetic barrier

BACKGROUND

The genomes of higher plants are, on the majority, polyploid, and hybridisation is more frequent in plants than in animals. Both polyploidisation and hybridisation contribute to increased variability within species, and may transfer adaptations between species in a changing environment. Studying these aspects of evolution within a diversified species complex could help to clarify overall spatial and temporal patterns of plant speciation. The Arabidopsis lyrata complex, which is closely related to the model plant Arabidopsis thaliana, is a perennial, outcrossing, herbaceous species complex with a circumpolar distribution in the Northern Hemisphere as well as a disjunct Central European distribution in relictual habitats. This species complex comprises three species and four subspecies, mainly diploids but also several tetraploids, including one natural hybrid. The complex is ecologically, but not fully geographically, separated from members of the closely related species complex of Arabidopsis halleri, and the evolutionary histories of both species compexes have largely been influenced by Pleistocene climate oscillations.

RESULTS

Using DNA sequence data from the nuclear encoded cytosolic phosphoglucoisomerase and Internal Transcribed Spacers 1 and 2 of the ribosomal DNA, as well as the trnL/F region from the chloroplast genome, we unravelled the phylogeography of the various taxonomic units of the A. lyrata complex. We demonstrate the existence of two major gene pools in Central Europe and Northern America. These two major gene pools are constructed from different taxonomic units. We also confirmed that A. kamchatica is the allotetraploid hybrid between A. lyrata and A. halleri, occupying the amphi-Beringian area in Eastern Asia and Northern America. This species closes the large distribution gap of the various other A. lyrata segregates. Furthermore, we revealed a threefold independent allopolyploid origin of this hybrid species in Japan, China, and Kamchatka.

CONCLUSIONS

Unglaciated parts of the Eastern Austrian Alps and arctic Eurasia, including Beringia, served as major glacial refugia of the Eurasian A. lyrata lineage, whereas A. halleri and its various subspecies probably survived in refuges in Central Europe and Eastern Asia with a large distribution gap in between. The North American A. lyrata lineage probably survived the glaciation in the southeast of North America. The dramatic climatic changes during glaciation and deglaciation cycles promoted not only secondary contact and formation of the allopolyploid hybrid A. kamchatica, but also provided the environment that allowed this species to fill a large geographic gap separating the two genetically different A. lyrata lineages from Eurasia and North America. With our example focusing on the evolutionary history of the A. lyrata species complex, we add substantial information to a broad evolutionary framework for future investigations within this emerging model system in molecular and evolutionary biology.

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.

Population structure and historical biogeography of European Arabidopsis lyrata

Understanding the natural history of model organisms is important for the effective use of their genomic resources. Arabidopsis lyrata has emerged as a useful plant for studying ecological and evolutionary genetics, based on its extensive natural variation, sequenced genome and close relationship to A. thaliana. We studied genetic diversity across the entire range of European Arabidopsis lyrata ssp. petraea, in order to explore how population history has influenced population structure. We sampled multiple populations from each region, using nuclear and chloroplast genome markers, and combined population genetic and phylogeographic approaches. Within-population diversity is substantial for nuclear allozyme markers (mean P=0.610, A(e)=1.580, H(e)=0.277) and significantly partitioned among populations (F(ST)=0.271). The Northern populations have modestly increased inbreeding (F(IS)=0.163 verses F(IS)=0.093), but retain comparable diversity to central European populations. Bottlenecks are common among central and northern Europe populations, indicating recent demographic history as a dominant factor in structuring the European diversity. Although the genetic structure was detected at all geographic scales, two clear differentiated units covering northern and central European areas (F(CT) =0.155) were identified by Bayesian analysis and supported by regional pairwise F(CT) calculations. A highly similar geographic pattern was observed from the distribution of chloroplast haplotypes, with the dominant northern haplotypes absent from central Europe. We conclude A. l. petraea's cold-tolerance and preference for disturbed habitats enabled glacial survival between the alpine and Nordic glaciers in central Europe and an additional cryptic refugium. While German populations are probable peri-glacial leftovers, Eastern Austrian populations have diversity patterns possibly compatible with longer-term survival.

The scale of population structure in Arabidopsis thaliana

The population structure of an organism reflects its evolutionary history and influences its evolutionary trajectory. It constrains the combination of genetic diversity and reveals patterns of past gene flow. Understanding it is a prerequisite for detecting genomic regions under selection, predicting the effect of population disturbances, or modeling gene flow. This paper examines the detailed global population structure of Arabidopsis thaliana. Using a set of 5,707 plants collected from around the globe and genotyped at 149 SNPs, we show that while A. thaliana as a species self-fertilizes 97% of the time, there is considerable variation among local groups. This level of outcrossing greatly limits observed heterozygosity but is sufficient to generate considerable local haplotypic diversity. We also find that in its native Eurasian range A. thaliana exhibits continuous isolation by distance at every geographic scale without natural breaks corresponding to classical notions of populations. By contrast, in North America, where it exists as an exotic species, A. thaliana exhibits little or no population structure at a continental scale but local isolation by distance that extends hundreds of km. This suggests a pattern for the development of isolation by distance that can establish itself shortly after an organism fills a new habitat range. It also raises questions about the general applicability of many standard population genetics models. Any model based on discrete clusters of interchangeable individuals will be an uneasy fit to organisms like A. thaliana which exhibit continuous isolation by distance on many scales.

Much of the modern field of population genetics is premised on particular models of what an organism's population structure is and how it behaves. The classic models generally start with the idea of a single randomly mating population that has reached an evolutionary equilibrium. Many models relax some of these assumptions, allowing for phenomena such as assortative mating, discrete sub-populations with migration, self-fertilization, and sex-ratio distortion. Virtually all models, however, have as their core premise the notion that there exist classes of exchangeable individuals each of which represents an identical, independent sample from that class' distribution. For certain organisms, such as Drosophila melanogaster, these models do an excellent job of describing how populations work. For other organisms, such as humans, these models can be reasonable approximations but require a great deal of care in assembling samples and can begin to break down as sampling becomes locally dense. For the vast majority of organisms the applicability of these models has never been investigated.

The origin of populations of Arabidopsis thaliana in China, based on the chloroplast DNA sequences

BACKGROUND

In the studies incorporating worldwide sampling of A. thaliana populations, the samples from East Asia, especially from China, were very scattered; and the studies focused on global patterns of cpDNA genetic variation among accessions of A. thaliana are very few. In this study, chloroplast DNA sequence variability was used to infer phylogenetic relationships among Arabidopsis thaliana accessions from around the world, with the emphasis on samples from China.

RESULTS

A data set comprising 77 accessions of A. thaliana, including 19 field-collected Chinese accessions together with three related species (A. arenosa, A. suecica, and Olimarabidopsis cabulica) as the out-group, was compiled. The analysis of the nucleotide sequences showed that the 77 accessions of A. thaliana were partitioned into two major differentiated haplotype classes (MDHCs). The estimated divergence time of the two MDHCs was about 0.39 mya. Forty-nine haplotypes were detected among the 77 accessions, which exhibited nucleotide diversity (pi) of 0.00169. The Chinese populations along the Yangtze River were characterized by five haplotypes, and the two accessions collected from the middle range of the Altai Mountains in China shared six specific variable sites.

CONCLUSIONS

The dimorphism in the chloroplast DNA could be due to founder effects during late Pleistocene glaciations and interglacial periods, although introgression cannot be ruled out. The Chinese populations along the Yangtze River may have dispersed eastwards to their present-day locations from the Himalayas. These populations originated from a common ancestor, and a rapid demographic expansion began approximately 90,000 years ago. Two accessions collected from the middle range of the Altai Mountains in China may have survived in a local refugium during late Pleistocene glaciations. The natural populations from China with specific genetic characteristics enriched the gene pools of global A. thaliana collections.

Dissecting the genetic control of natural variation in salt tolerance of Arabidopsis thaliana accessions

Many accessions (ecotypes) of Arabidopsis have been collected. Although few differences exist among their nucleotide sequences, these subtle differences induce large genetic variation in phenotypic traits such as stress tolerance and flowering time. To understand the natural variability in salt tolerance, large-scale soil pot experiments were performed to evaluate salt tolerance among 350 Arabidopsis thaliana accessions. The evaluation revealed a wide variation in the salt tolerance among accessions. Several accessions, including Bu-5, Bur-0, Ll-1, Wl-0, and Zu-0, exhibited marked stress tolerance compared with a salt-sensitive experimental accession, Col-0. The salt-tolerant accessions were also evaluated by agar plate assays. The data obtained by the large-scale assay correlated well with the results of a salt acclimation (SA) assay, in which plants were transferred to high-salinity medium following placement on moderate-salinity medium for 7 d. Genetic analyses indicated that the salt tolerance without SA is a quantitative trait under polygenic control, whereas salt tolerance with SA is regulated by a single gene located on chromosome 5 that is common among the markedly salt-tolerant accessions. These results provide important information for understanding the mechanisms underlying natural variation of salt tolerance in Arabidopsis.

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.

The allopolyploid Arabidopsis kamchatica originated from multiple individuals of Arabidopsis lyrata and Arabidopsis halleri

Polyploidization, or genome duplication, has played a critical role in the diversification of animals, fungi and plants. Little is known about the population structure and multiple origins of polyploid species because of the difficulty in identifying multiple homeologous nuclear genes. The allotetraploid species Arabidopsis kamchatica is closely related to the model species Arabidopsis thaliana and is distributed in a broader climatic niche than its parental species. Here, we performed direct sequencing of homeologous pairs of the low-copy nuclear genes WER and CHS by designing homeolog-specific primers, and obtained also chloroplast and ribosomal internal transcribed spacer sequences. Phylogenetic analysis showed that 50 individuals covering the distribution range including North America are allopolyploids derived from Arabidopsis lyrata and Arabidopsis halleri. Three major clusters within A. kamchatica were detected using Bayesian clustering. One cluster has widespread distribution. The other two are restricted to the southern part of the distribution range including Japan, where the parent A. lyrata is not currently distributed. This suggests that the mountains in Central Honshu and surrounding areas in Japan served as refugia during glacial-interglacial cycles and retained this diversity. We also found that multiple haplotypes of nuclear and chloroplast sequences of A. kamchatica are identical to those of their parental species. This indicates that multiple diploid individuals contributed to the origin of A. kamchatica. The haplotypes of low-copy nuclear genes in Japan suggest independent polyploidization events rather than introgression. Our findings suggest that self-compatibility and gene silencing occurred independently in different origins.

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.

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.

Pleiotropic effects of environment-specific adaptation in Arabidopsis thaliana

Local adaptation may be important for the preservation of genetic diversity and the promotion of speciation. However, local adaptation may also constrain establishment in different environments. The consequences of local adaptation depend strongly on the pleiotropic effects of the genes involved in adaptation. Here, we investigated the pleiotropic effects of the genetic response to selection in outbred lines of Arabidopsis artificially selected to flower earlier under both winter- and spring-annual simulated conditions. The consequences of adaptation were evaluated by reciprocally transplanting selected and control lines between the two conditions. Selected lines always flower earlier than their controls, independent of growing conditions. However, selected lines, growing in the same condition in which they were selected, flower earlier than plants selected in the alternative environment. Plants selected to flower earlier in spring produce more fruits than controls when growing in the spring, and less fruits when growing in the winter; indicating that local adaptation has negative pleiotropic effects in another environment. Our results indicate that local adaptation can arise even when selection targets the same trait in the same direction. Furthermore, it suggests that adaptation under the two different environments can generate fitness trade-offs that can maintain genetic variation for flowering time.