Variation and plasticity in life-history traits and fitness of wild Arabidopsis thaliana populations are not related to their genotypic and ecological diversity.

BACKGROUND: Despite its implications for population dynamics and evolution, the relationship between genetic and phenotypic variation in wild populations remains unclear. Here, we estimated variation and plasticity in life-history traits and fitness of the annual plant Arabidopsis thaliana in two common garden experiments that differed in environmental conditions. We used up to 306 maternal inbred lines from six Iberian populations characterized by low and high genotypic (based on whole-genome sequences) and ecological (vegetation type) diversity. RESULTS: Low and high genotypic and ecological diversity was found in edge and core Iberian environments, respectively. Given that selection is expected to be stronger in edge environments and that ecological diversity may enhance both phenotypic variation and plasticity, we expected genotypic diversity to be positively associated with phenotypic variation and plasticity. However, maternal lines, irrespective of the genotypic and ecological diversity of their population of origin, exhibited a substantial amount of phenotypic variation and plasticity for all traits. Furthermore, all populations harbored maternal lines with canalization (robustness) or sensitivity in response to harsher environmental conditions in one of the two experiments. CONCLUSIONS: Overall, we conclude that the environmental attributes of each population probably determine their genotypic diversity, but all populations maintain substantial phenotypic variation and plasticity for all traits, which represents an asset to endure in changing environments.

Differential environmental and genomic architectures shape the natural diversity for trichome patterning and morphology in different Arabidopsis organs.

Despite the adaptive and taxonomic relevance of the natural diversity for trichome patterning and morphology, the molecular and evolutionary mechanisms underlying these traits remain mostly unknown, particularly in organs other than leaves. In this study, we address the ecological, genetic and molecular bases of the natural variation for trichome patterning and branching in multiple organs of Arabidopsis (Arabidopsis thaliana). To this end, we characterized a collection of 191 accessions and carried out environmental and genome-wide association (GWA) analyses. Trichome amount in different organs correlated negatively with precipitation in distinct seasons, thus suggesting a precise fit between trichome patterning and climate throughout the Arabidopsis life cycle. In addition, GWA analyses showed small overlapping between the genes associated with different organs, indicating partly independent genetic bases for vegetative and reproductive phases. These analyses identified a complex locus on chromosome 2, where two adjacent MYB genes (ETC2 and TCL1) displayed differential effects on trichome patterning in several organs. Furthermore, analyses of transgenic lines carrying different natural alleles demonstrated that TCL1 accounts for the variation for trichome patterning in all organs, and for stem trichome branching. By contrast, two other MYB genes (TRY and GL1), mainly showed effects on trichome patterning or branching, respectively.

Depicting the phenotypic space of the annual plant Diplotaxis acris in hyperarid deserts.

The phenotypic space encompasses the assemblage of trait combinations yielding well-suited integrated phenotypes. At the population level, understanding the phenotypic space structure requires the quantification of among- and within-population variations in traits and the correlation pattern among them. Here, we studied the phenotypic space of the annual plant Diplotaxis acris occurring in hyperarid deserts. Given the advance of warming and aridity in vast regions occupied by drylands, D. acris can indicate the successful evolutionary trajectory that many other annual plant species may follow in expanding drylands. To this end, we conducted a greenhouse experiment with 176 D. acris individuals from five Saudi populations to quantify the genetic component of variation in architectural and life history traits. We found low among-population divergence but high among-individual variation in all traits. In addition, all traits showed a high degree of genetic determination in our study experimental conditions. We did not find significant effects of recruitment and fecundity on fitness. Finally, all architectural traits exhibited a strong correlation pattern among them, whereas for life history traits, only higher seed germination implied earlier flowering. Seed weight appeared to be an important trait in D. acris as individuals with heavier seeds tended to advance flowering and have a more vigorous branching pattern, which led to higher fecundity. Population divergence in D. acris might be constrained by the severity of the hyperarid environment, but populations maintain high among-individual genetic variation in all traits. Furthermore, D. acris showed phenotypic integration for architectural traits and, to a lesser extent, for life history traits. Overall, we hypothesize that D. acris may be fine-tuned to its demanding extreme environments. Evolutionary speaking, annual plants facing increasing warming, aridity, and environmental seasonality might modify their phenotypic spaces toward new phenotypic configurations strongly dominated by correlated architectural traits enhancing fecundity and seed-related traits advancing flowering time.

Introduction to the Special Issue: The ecology and genetics of population differentiation in plants.

Population differentiation is a pervasive process in nature. At present, evolutionary studies on plant population differentiation address key questions by undertaking joint ecological and genetic approaches and employing a combination of molecular and experimental means. In this special issue, we gathered a collection of papers dealing with various ecological and genetic aspects of population differentiation in plants. In particular, this special issue encompasses eight research articles and two reviews covering a wide array of worldwide environments, plant functional types, genetic and genomic approaches, and common garden experiments to quantify molecular and/or quantitative trait differentiation in plant populations. Overall, this special issue stresses the validity of traditional evolutionary studies focused on plant populations, whilst emphasizing the integration of classical biological disciplines and state-of-the-art molecular techniques into a unique toolkit for evolutionary plant research.

MYB transcription factors drive evolutionary innovations in Arabidopsis fruit trichome patterning.

Both inter- and intra-specific diversity has been described for trichome patterning in fruits, which is presumably involved in plant adaptation. However, the mechanisms underlying this developmental trait have been hardly addressed. Here we examined natural populations of Arabidopsis (Arabidopsis thaliana) that develop trichomes in fruits and pedicels, phenotypes previously not reported in the Arabidopsis genus. Genetic analyses identified five loci, MALAMBRUNO 1-5 (MAU1-5), with MAU2, MAU3, and MAU5 showing strong epistatic interactions that are necessary and sufficient to display these traits. Functional characterization of these three loci revealed cis-regulatory mutations in TRICHOMELESS1 and TRIPTYCHON, as well as a structural mutation in GLABRA1. Therefore, the multiple mechanisms controlled by three MYB transcription factors of the core regulatory network for trichome patterning have jointly been modulated to trigger trichome development in fruits. Furthermore, analyses of worldwide accessions showed that these traits and mutations only occur in a highly differentiated relict lineage from the Iberian Peninsula. In addition, these traits and alleles were associated with low spring precipitation, which suggests that trichome development in fruits and pedicels might be involved in climatic adaptation. Thus, we show that the combination of synergistic mutations in a gene regulatory circuit has driven evolutionary innovations in fruit trichome patterning in Arabidopsis.

Temporal migration rates affect the genetic structure of populations in the biennial Erysimum mediohispanicum with reproductive asynchrony.

Migration is a process with important implications for the genetic structure of populations. However, there is an aspect of migration seldom investigated in plants: migration between temporally isolated groups of individuals within the same geographic population. The genetic implications of temporal migration can be particularly relevant for semelparous organisms, which are those that reproduce only once in a lifetime after a certain period of growth. In this case, reproductive asynchrony in individuals of the same population generates demes of individuals differing in their developmental stage (non-reproductive and reproductive). These demes are connected by temporal migrants, that is, individuals that become annually asynchronous with respect to the rest of individuals of their same deme. Here, we investigated the extent of temporal migration and its effects on temporal genetic structure in the biennial plant Erysimum mediohispanicum. To this end, we conducted two independent complementary approaches. First, we empirically estimated temporal migration rates and temporal genetic structure in four populations of E. mediohispanicum during three consecutive years using nuclear microsatellites markers. Second, we developed a demographic genetic simulation model to assess genetic structure for different migration scenarios differing in temporal migration rates and their occurrence probabilities. We hypothesized that genetic structure decreased with increasing temporal migration rates due to the homogenizing effect of migration. Empirical and modelling results were consistent and indicated a U-shape relationship between genetic structure and temporal migration rates. Overall, they indicated the existence of temporal genetic structure and that such genetic structure indeed decreased with increasing temporal migration rates. However, genetic structure increased again at high temporal migration rates. The results shed light into the effects of reproductive asynchrony on important population genetic parameters. Our study contributes to unravel the complexity of some processes that may account for genetic diversity and genetic structure of natural populations.

Ecological, genetic and evolutionary drivers of regional genetic differentiation in Arabidopsis thaliana.

BACKGROUND: Disentangling the drivers of genetic differentiation is one of the cornerstones in evolution. This is because genetic diversity, and the way in which it is partitioned within and among populations across space, is an important asset for the ability of populations to adapt and persist in changing environments. We tested three major hypotheses accounting for genetic differentiation-isolation-by-distance (IBD), isolation-by-environment (IBE) and isolation-by-resistance (IBR)-in the annual plant Arabidopsis thaliana across the Iberian Peninsula, the region with the largest genomic diversity. To that end, we sampled, genotyped with genome-wide SNPs, and analyzed 1772 individuals from 278 populations distributed across the Iberian Peninsula. RESULTS: IBD, and to a lesser extent IBE, were the most important drivers of genetic differentiation in A. thaliana. In other words, dispersal limitation, genetic drift, and to a lesser extent local adaptation to environmental gradients, accounted for the within- and among-population distribution of genetic diversity. Analyses applied to the four Iberian genetic clusters, which represent the joint outcome of the long demographic and adaptive history of the species in the region, showed similar results except for one cluster, in which IBR (a function of landscape heterogeneity) was the most important driver of genetic differentiation. Using spatial hierarchical Bayesian models, we found that precipitation seasonality and topsoil pH chiefly accounted for the geographic distribution of genetic diversity in Iberian A. thaliana. CONCLUSIONS: Overall, the interplay between the influence of precipitation seasonality on genetic diversity and the effect of restricted dispersal and genetic drift on genetic differentiation emerges as the major forces underlying the evolutionary trajectory of Iberian A. thaliana.

Microsatellite marker development in the crop wild relative Linum bienne using genome skimming.

PREMISE: Nuclear microsatellite markers were developed for Linum bienne, the sister species of the crop L. usitatissimum, to provide molecular genetic tools for the investigation of L. bienne genetic diversity and structure. METHODS AND RESULTS: Fifty microsatellite loci were identified in L. bienne by means of genome skimming, and 44 loci successfully amplified. Of these, 16 loci evenly spread across the L. usitatissimum reference nuclear genome were used for genotyping six L. bienne populations. Excluding one monomorphic locus, the number of alleles per locus ranged from two to 12. Four out of six populations harbored private alleles. The levels of expected and observed heterozygosity were 0.076 to 0.667 and 0.000 to 1.000, respectively. All 16 loci successfully cross-amplified in L. usitatissimum. CONCLUSIONS: The 16 microsatellite loci developed here can be used for population genetic studies in L. bienne, and 28 additional loci that successfully amplified are available for further testing.

Genetic rescue by distant trees mitigates qualitative pollen limitation imposed by fine-scale spatial genetic structure.

Restricted seed dispersal frequently leads to fine-scale spatial genetic structure (i.e., FSGS) within plant populations. Depending on its spatial extent and the mobility of pollinators, this inflated kinship at the immediate neighbourhood can critically impoverish pollen quality. Despite the common occurrence of positive FSGS within plant populations, our knowledge regarding the role of long-distance pollination preventing reproductive failure is still limited. Using microsatellite markers, we examined the existence of positive FSGS in two low-density populations of the tree Pyrus bourgaeana. We also designed controlled crosses among trees differing in their kinship to investigate the effects of increased local kinship on plant reproduction. We used six pollination treatments and fully monitored fruit production, fruit and seed weight, proportion of mature seeds per fruit, and seed germination. Our results revealed positive FSGS in both study populations and lower fruit initiation in flowers pollinated with pollen from highly-genetically related individuals within the neighbourhood, with this trend intensifying as the fruit development progressed. Besides, open-pollinated flowers exhibited lower performance compared to those pollinated by distant pollen donors, suggesting intense qualitative pollen limitation in natural populations. We found positive fine-scale spatial genetic structure is translated into impoverished pollen quality from nearby pollen donors which negatively impacts the reproductive success of trees in low-density populations. Under this scenario of intrapopulation genetic rescue by distant pollen donors, the relevance of highly-mobile pollinators for connecting spatially and genetically distant patches of trees may be crucial to safeguarding population recruitment.

A hierarchical Bayesian Beta regression approach to study the effects of geographical genetic structure and spatial autocorrelation on species distribution range shifts.

Global climate change (GCC) may be causing distribution range shifts in many organisms worldwide. Multiple efforts are currently focused on the development of models to better predict distribution range shifts due to GCC. We addressed this issue by including intraspecific genetic structure and spatial autocorrelation (SAC) of data in distribution range models. Both factors reflect the joint effect of ecoevolutionary processes on the geographical heterogeneity of populations. We used a collection of 301 georeferenced accessions of the annual plant Arabidopsis thaliana in its Iberian Peninsula range, where the species shows strong geographical genetic structure. We developed spatial and nonspatial hierarchical Bayesian models (HBMs) to depict current and future distribution ranges for the four genetic clusters detected. We also compared the performance of HBMs with Maxent (a presence-only model). Maxent and nonspatial HBMs presented some shortcomings, such as the loss of accessions with high genetic admixture in the case of Maxent and the presence of residual SAC for both. As spatial HBMs removed residual SAC, these models showed higher accuracy than nonspatial HBMs and handled the spatial effect on model outcomes. The ease of modelling and the consistency among model outputs for each genetic cluster was conditioned by the sparseness of the populations across the distribution range. Our HBMs enrich the toolbox of software available to evaluate GCC-induced distribution range shifts by considering both genetic heterogeneity and SAC, two inherent properties of any organism that should not be overlooked.

Quantifying temporal change in plant population attributes: insights from a resurrection approach.

Rapid evolution in annual plants can be quantified by comparing phenotypic and genetic changes between past and contemporary individuals from the same populations over several generations. Such knowledge will help understand the response of plants to rapid environmental shifts, such as the ones imposed by global climate change. To that end, we undertook a resurrection approach in Spanish populations of the annual plant Arabidopsis thaliana that were sampled twice over a decade. Annual weather records were compared to their historical records to extract patterns of climatic shifts over time. We evaluated the differences between samplings in flowering time, a key life-history trait with adaptive significance, with a field experiment. We also estimated genetic diversity and differentiation based on neutral nuclear markers and nucleotide diversity in candidate flowering time (FRI and FLC) and seed dormancy (DOG1) genes. The role of genetic drift was estimated by computing effective population sizes with the temporal method. Overall, two climatic scenarios were detected: intense warming with increased precipitation and moderate warming with decreased precipitation. The average flowering time varied little between samplings. Instead, within-population variation in flowering time exhibited a decreasing trend over time. Substantial temporal changes in genetic diversity and differentiation were observed with both nuclear microsatellites and candidate genes in all populations, which were interpreted as the result of natural demographic fluctuations. We conclude that drought stress caused by moderate warming with decreased precipitation may have the potential to reduce within-population variation in key life-cycle traits, perhaps as a result of stabilizing selection on them, and to constrain the genetic differentiation over time. Besides, the demographic behaviour of populations probably accounts for the substantial temporal patterns of genetic variation, while keeping rather constant those of phenotypic variation.

Spatio-temporal variation in fitness responses to contrasting environments in Arabidopsis thaliana.

The evolutionary response of organisms to global climate change is expected to be strongly conditioned by preexisting standing genetic variation. In addition, natural selection imposed by global climate change on fitness-related traits can be heterogeneous over time. We estimated selection of life-history traits of an entire genetic lineage of the plant Arabidopsis thaliana occurring in north-western Iberian Peninsula that were transplanted over multiple years into two environmentally contrasting field sites in southern Spain, as southern environments are expected to move progressively northwards with climate change in the Iberian Peninsula. The results indicated that natural selection on flowering time prevailed over that on recruitment. Selection favored early flowering in six of eight experiments and late flowering in the other two. Such heterogeneity of selection for flowering time might be a powerful mechanism for maintaining genetic diversity in the long run. We also found that north-western A. thaliana accessions from warmer environments exhibited higher fitness and higher phenotypic plasticity for flowering time in southern experimental facilities. Overall, our transplant experiments suggested that north-western Iberian A. thaliana has the means to cope with increasingly warmer environments in the region as predicted by trends in global climate change models.

African genomes illuminate the early history and transition to selfing in Arabidopsis thaliana.

Over the past 20 y, many studies have examined the history of the plant ecological and molecular model, Arabidopsis thaliana, in Europe and North America. Although these studies informed us about the recent history of the species, the early history has remained elusive. In a large-scale genomic analysis of African A. thaliana, we sequenced the genomes of 78 modern and herbarium samples from Africa and analyzed these together with over 1,000 previously sequenced Eurasian samples. In striking contrast to expectations, we find that all African individuals sampled are native to this continent, including those from sub-Saharan Africa. Moreover, we show that Africa harbors the greatest variation and represents the deepest history in the A. thaliana lineage. Our results also reveal evidence that selfing, a major defining characteristic of the species, evolved in a single geographic region, best represented today within Africa. Demographic inference supports a model in which the ancestral A. thaliana population began to split by 120-90 kya, during the last interglacial and Abbassia pluvial, and Eurasian populations subsequently separated from one another at around 40 kya. This bears striking similarities to the patterns observed for diverse species, including humans, implying a key role for climatic events during interglacial and pluvial periods in shaping the histories and current distributions of a wide range of species.

Tackling intraspecific genetic structure in distribution models better reflects species geographical range.

Genetic diversity provides insight into heterogeneous demographic and adaptive history across organisms' distribution ranges. For this reason, decomposing single species into genetic units may represent a powerful tool to better understand biogeographical patterns as well as improve predictions of the effects of GCC (global climate change) on biodiversity loss. Using 279 georeferenced Iberian accessions, we used classes of three intraspecific genetic units of the annual plant Arabidopsis thaliana obtained from the genetic analyses of nuclear SNPs (single nucleotide polymorphisms), chloroplast SNPs, and the vernalization requirement for flowering. We used SDM (species distribution models), including climate, vegetation, and soil data, at the whole-species and genetic-unit levels. We compared model outputs for present environmental conditions and with a particularly severe GCC scenario. SDM accuracy was high for genetic units with smaller distribution ranges. Kernel density plots identified the environmental variables underpinning potential distribution ranges of genetic units. Combinations of environmental variables accounted for potential distribution ranges of genetic units, which shrank dramatically with GCC at almost all levels. Only two genetic clusters increased their potential distribution ranges with GCC. The application of SDM to intraspecific genetic units provides a detailed picture on the biogeographical patterns of distinct genetic groups based on different genetic criteria. Our approach also allowed us to pinpoint the genetic changes, in terms of genetic background and physiological requirements for flowering, that Iberian A. thaliana may experience with a GCC scenario applying SDM to intraspecific genetic units.

Altitudinal and climatic associations of seed dormancy and flowering traits evidence adaptation of annual life cycle timing in Arabidopsis thaliana.

The temporal control or timing of the life cycle of annual plants is presumed to provide adaptive strategies to escape harsh environments for survival and reproduction. This is mainly determined by the timing of germination, which is controlled by the level of seed dormancy, and of flowering initiation. However, the environmental factors driving the evolution of plant life cycles remain largely unknown. To address this question we have analysed nine quantitative life history traits, in a native regional collection of 300 wild accessions of Arabidopsis thaliana. Seed dormancy and flowering time were negatively correlated, indicating that these traits have coevolved. In addition, environmental-phenotypic analyses detected strong altitudinal and climatic clines for most life history traits. Overall, accessions showing life cycles with early flowering, small seeds, high seed dormancy and slow germination rate were associated with locations exposed to high temperature, low summer precipitation and high radiation. Furthermore, we analysed the expression level of the positive regulator of seed dormancy DELAY OF GERMINATION 1 (DOG1), finding similar but weaker altitudinal and climatic patterns than seed dormancy. Therefore, DOG1 regulatory mutations are likely to provide a quantitative molecular mechanism for the adaptation of A. thaliana life cycle to altitude and climate.

Environmental and genetic interactions reveal FLOWERING LOCUS C as a modulator of the natural variation for the plasticity of flowering in Arabidopsis.

The timing of flowering initiation depends strongly on the environment, a property termed as the plasticity of flowering. Such plasticity determines the adaptive potential of plants because it provides phenotypic buffer against environmental changes, and its natural variation contributes to evolutionary adaptation. We addressed the genetic mechanisms of the natural variation for this plasticity in Arabidopsis thaliana by analysing a population of recombinant inbred lines derived from Don-0 and Ler accessions collected from distinct climates. Quantitative trait locus (QTL) mapping in four environmental conditions differing in photoperiod, vernalization treatment and ambient temperature detected the folllowing: (i) FLOWERING LOCUS C (FLC) as a large effect QTL affecting flowering time differentially in all environments; (ii) numerous QTL displaying smaller effects specifically in some conditions; and (iii) significant genetic interactions between FLC and other loci. Hence, the variation for the plasticity of flowering is determined by a combination of environmentally sensitive and specific QTL, and epistasis. Analysis of FLC from Don identified a new and more active allele likely caused by a cis-regulatory deletion covering the non-coding RNA COLDAIR. Further characterization of four FLC natural alleles showed different environmental and genetic interactions. Thus, FLC appears as a major modulator of the natural variation for the plasticity of flowering to multiple environmental factors.

Natural Variation Identifies ICARUS1, a Universal Gene Required for Cell Proliferation and Growth at High Temperatures in Arabidopsis thaliana.

Plants are highly sensitive to environmental changes and even small variations in ambient temperature have severe consequences on their growth and development. Temperature affects multiple aspects of plant development, but the processes and mechanisms underlying thermo-sensitive growth responses are mostly unknown. Here we exploit natural variation in Arabidopsis thaliana to identify and characterize novel components and processes mediating thermo-sensitive growth responses in plants. Phenotypic screening of wild accessions identified several strains displaying pleiotropic growth defects, at cellular and organism levels, specifically at high ambient temperatures. Positional cloning and characterization of the underlying gene revealed that ICARUS1 (ICA1), which encodes a protein of the tRNAHis guanylyl transferase (Thg1) superfamily, is required for plant growth at high temperatures. Transcriptome and gene marker analyses together with DNA content measurements show that ICA1 loss-of-function results in down regulation of cell cycle associated genes at high temperatures, which is linked with a block in G2/M transition and endoreduplication. In addition, plants with mutations in ICA1 show enhanced sensitivity to DNA damage. Characterization of additional strains that carry lesions in ICA1, but display normal growth, shows that alternative splicing is likely to alleviate the deleterious effects of some natural mutations. Furthermore, analyses of worldwide and regional collections of natural accessions indicate that ICA1 loss-of-function has arisen several times independently, and that these occur at high frequency in some local populations. Overall our results suggest that ICA1-mediated-modulation of fundamental processes such as tRNAHis maturation, modify plant growth responses to temperature changes in a quantitative and reversible manner, in natural populations.

Range-wide population genetics and variation in morph ratio in style-dimorphic Narcissus papyraceus.

UNLABELLED: • PREMISE OF THE STUDY: Theoretical models state that natural selection and mating patterns account for floral morph ratio in style-polymorphic plants. However, the demographic history of populations can also influence variation in morph ratios. If so, we hypothesize an association between the morph ratios and the genetic structure across populations.• METHODS: We used nuclear microsatellites to assess genetic variation and structure in populations of Narcissus papyraceus, a style-dimorphic plant whose floral morph ratios (L-morph to S-morph) gradually vary throughout its distribution range in the southwestern Mediterranean Basin. We implemented analyses to relate the genetic features of populations with their morph ratios.• KEY RESULTS: We found greater frequencies of the S-morph in central populations and declining frequencies toward the periphery. This geographic pattern was not associated with the genetic structure of populations. Instead, we found two distinct genetic groups, mainly separated by the Strait of Gibraltar, with a mixture of morph ratios within each one. Overall, there was a weak genetic structure. Genetic diversity was greater in central and southern dimorphic populations than in northern L-monomorphic populations.• CONCLUSIONS: Altogether, our results do not support the hypothesis that the demographic history of populations can account for the observed geographical pattern of morph ratios in N. papyraceus. We suggest that adaptive processes shown in previous studies in the species are the main determinant of the existing variation in the morph composition of populations.

Deciphering the adjustment between environment and life history in annuals: lessons from a geographically-explicit approach in Arabidopsis thaliana.

The role that different life-history traits may have in the process of adaptation caused by divergent selection can be assessed by using extensive collections of geographically-explicit populations. This is because adaptive phenotypic variation shifts gradually across space as a result of the geographic patterns of variation in environmental selective pressures. Hence, large-scale experiments are needed to identify relevant adaptive life-history traits as well as their relationships with putative selective agents. We conducted a field experiment with 279 geo-referenced accessions of the annual plant Arabidopsis thaliana collected across a native region of its distribution range, the Iberian Peninsula. We quantified variation in life-history traits throughout the entire life cycle. We built a geographic information system to generate an environmental data set encompassing climate, vegetation and soil data. We analysed the spatial autocorrelation patterns of environmental variables and life-history traits, as well as the relationship between environmental and phenotypic data. Almost all environmental variables were significantly spatially autocorrelated. By contrast, only two life-history traits, seed weight and flowering time, exhibited significant spatial autocorrelation. Flowering time, and to a lower extent seed weight, were the life-history traits with the highest significant correlation coefficients with environmental factors, in particular with annual mean temperature. In general, individual fitness was higher for accessions with more vigorous seed germination, higher recruitment and later flowering times. Variation in flowering time mediated by temperature appears to be the main life-history trait by which A. thaliana adjusts its life history to the varying Iberian environmental conditions. The use of extensive geographically-explicit data sets obtained from field experiments represents a powerful approach to unravel adaptive patterns of variation. In a context of current global warming, geographically-explicit approaches, evaluating the match between organisms and the environments where they live, may contribute to better assess and predict the consequences of global warming.

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.