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

Can plants build their niche through modulation of soil microbial activities linked with nitrogen cycling? A test with Arabidopsis thaliana

In natural systems, different plant species have been shown to modulate specific nitrogen (N) cycling processes so as to meet their N demand, thereby potentially influencing their own niche. This phenomenon might go beyond plant interactions with symbiotic microorganisms and affect the much less explored plant interactions with free-living microorganisms involved in soil N cycling, such as nitrifiers and denitrifiers. Here, we investigated variability in the modulation of soil nitrifying and denitrifying enzyme activities (NEA and DEA, respectively), and their ratio (NEA : DEA), across 193 Arabidopsis thaliana accessions. We studied the genetic and environmental determinants of such plant-soil interactions, and effects on plant biomass production in the next generation. We found that NEA, DEA, and NEA : DEA varied c. 30-, 15- and 60-fold, respectively, among A. thaliana genotypes and were related to genes linked with stress response, flowering, and nitrate nutrition, as well as to soil parameters at the geographic origin of the analysed genotypes. Moreover, plant-mediated N cycling activities correlated with the aboveground biomass of next-generation plants in home vs away nonautoclaved soil, suggesting a transgenerational impact of soil biotic conditioning on plant performance. Altogether, these findings suggest that nutrient-based plant niche construction may be much more widespread than previously thought.

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.

Transcriptome wide analyses reveal intraspecific diversity in thermal stress responses of a dominant habitat-forming species

The impact of climate change on biodiversity has stimulated the need to understand environmental stress responses, particularly for ecosystem engineers whose responses to climate affect large numbers of associated organisms. Distinct species differ substantially in their resilience to thermal stress but there are also within-species variations in thermal tolerance for which the molecular mechanisms underpinning such variation remain largely unclear. Intertidal mussels are well-known for their role as ecosystem engineers. First, we exposed two genetic lineages of the intertidal mussel Perna perna to heat stress treatments in air and water. Next, we ran a high throughput RNA sequencing experiment to identify differences in gene expression between the thermally resilient eastern lineage and the thermally sensitive western lineage. We highlight different thermal tolerances that concord with their distributional ranges. Critically, we also identified lineage-specific patterns of gene expression under heat stress and revealed intraspecific differences in the underlying transcriptional pathways in response to warmer temperatures that are potentially linked to the within-species differences in thermal tolerance. Beyond the species, we show how unravelling within-species variability in mechanistic responses to heat stress promotes a better understanding of global evolutionary trajectories of the species as a whole in response to changing climate.

Forecasting trait responses in novel environments to aid seed provenancing under climate change

Revegetation projects face the major challenge of sourcing optimal plant material. This is often done with limited information about plant performance and increasingly requires factoring resilience to climate change. Functional traits can be used as quantitative indices of plant performance and guide seed provenancing, but trait values expected under novel conditions are often unknown. To support climate-resilient provenancing efforts, we develop a trait prediction model that integrates the effect of genetic variation with fine-scale temperature variation. We train our model on multiple field plantings of Arabidopsis thaliana and predict two relevant fitness traits-days-to-bolting and fecundity-across the species' European range. Prediction accuracy was high for days-to-bolting and moderate for fecundity, with the majority of trait variation explained by temperature differences between plantings. Projection under future climate predicted a decline in fecundity, although this response was heterogeneous across the range. In response, we identified novel genotypes that could be introduced to genetically offset the fitness decay. Our study highlights the value of predictive models to aid seed provenancing and improve the success of revegetation projects.

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.

Adaptive significance of flowering time variation across natural seasonal environments in Arabidopsis thaliana

The relevance of flowering time variation and plasticity to climate adaptation requires a comprehensive empirical assessment. We investigated natural selection and the genetic architecture of flowering time in Arabidopsis through field experiments in Europe across multiple sites and seasons. We estimated selection for flowering time, plasticity and canalization. Loci associated with flowering time, plasticity and canalization by genome-wide association studies were tested for a geographic signature of climate adaptation. Selection favored early flowering and increased canalization, except at the northernmost site, but was rarely detected for plasticity. Genome-wide association studies revealed significant associations with flowering traits and supported a substantial polygenic inheritance. Alleles associated with late flowering, including functional FRIGIDA variants, were more common in regions experiencing high annual temperature variation. Flowering time plasticity to fall vs spring and summer environments was associated with GIGANTEA SUPPRESSOR 5, which promotes early flowering under decreasing day length and temperature. The finding that late flowering genotypes and alleles are associated with climate is evidence for past adaptation. Real-time phenotypic selection analysis, however, reveals pervasive contemporary selection for rapid flowering in agricultural settings across most of the species range. The response to this selection may involve genetic shifts in environmental cuing compared to the ancestral state.

Genome-wide association mapping of transcriptome variation in Mimulus guttatus indicates differing patterns of selection on cis- versus trans-acting mutations

We measured the floral bud transcriptome of 151 fully sequenced lines of Mimulus guttatus from one natural population. Thousands of single nucleotide polymorphisms (SNPs) are implicated as transcription regulators, but there is a striking difference in the allele frequency spectrum of cis-acting and trans-acting mutations. Cis-SNPs have intermediate frequencies (consistent with balancing selection) while trans-SNPs exhibit a rare-alleles model (consistent with purifying selection). This pattern only becomes clear when transcript variation is normalized on a gene-to-gene basis. If a global normalization is applied, as is typically in RNAseq experiments, asymmetric transcript distributions combined with "rarity disequilibrium" produce a superabundance of false positives for trans-acting SNPs. To explore the cause of purifying selection on trans-acting mutations, we identified gene expression modules as sets of coexpressed genes. The extent to which trans-acting mutations influence modules is a strong predictor of allele frequency. Mutations altering expression of genes with high "connectedness" (those that are highly predictive of the representative module expression value) have the lowest allele frequency. The expression modules can also predict whole-plant traits such as flower size. We find that a substantial portion of the genetic (co)variance among traits can be described as an emergent property of genetic effects on expression modules.

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.

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.

Genetic architecture and adaptation of flowering time among environments

The genetic basis of flowering time changes across environments, and pleiotropy may limit adaptive evolution of populations in response to local conditions. However, little information is known about how genetic architecture changes among environments. We used genome-wide association studies (GWAS) in Boechera stricta (Graham) Al-Shehbaz, a relative of Arabidopsis, to examine flowering variation among environments and associations with climate conditions in home environments. Also, we used molecular population genetics to search for evidence of historical natural selection. GWAS found 47 significant quantitative trait loci (QTLs) that influence flowering time in one or more environments, control plastic changes in phenology between experiments, or show associations with climate in sites of origin. Genetic architecture of flowering varied substantially among environments. We found that some pairs of QTLs showed similar patterns of pleiotropy across environments. A large-effect QTL showed molecular signatures of adaptive evolution and is associated with climate in home environments. The derived allele at this locus causes later flowering and predominates in sites with greater water availability. This work shows that GWAS of climate associations and ecologically important traits across diverse environments can be combined with molecular signatures of natural selection to elucidate ecological genetics of adaptive evolution.

Pearl millet genomic vulnerability to climate change in West Africa highlights the need for regional collaboration

Climate change is already affecting agro-ecosystems and threatening food security by reducing crop productivity and increasing harvest uncertainty. Mobilizing crop diversity could be an efficient way to mitigate its impact. We test this hypothesis in pearl millet, a nutritious staple cereal cultivated in arid and low-fertility soils in sub-Saharan Africa. We analyze the genomic diversity of 173 landraces collected in West Africa together with an extensive climate dataset composed of metrics of agronomic importance. Mapping the pearl millet genomic vulnerability at the 2050 horizon based on the current genomic-climate relationships, we identify the northern edge of the current areas of cultivation of both early and late flowering varieties as being the most vulnerable to climate change. We predict that the most vulnerable areas will benefit from using landraces that already grow in equivalent climate conditions today. However, such seed-exchange scenarios will require long distance and trans-frontier assisted migrations. Leveraging genetic diversity as a climate mitigation strategy in West Africa will thus require regional collaboration.

Replacement of local crops with alternative varieties adapted to future conditions may improve food security under climate change. Here the authors apply landscape genomics and ensemble climate modelling to pearl millet in West Africa, supporting the potential of transfrontier assisted seed exchange.

Natural variation in autumn expression is the major adaptive determinant distinguishing Arabidopsis FLC haplotypes

In Arabidopsis thaliana, winter is registered during vernalization through the temperature-dependent repression and epigenetic silencing of floral repressor FLOWERING LOCUS C (FLC). Natural Arabidopsis accessions show considerable variation in vernalization. However, which aspect of the FLC repression mechanism is most important for adaptation to different environments is unclear. By analysing FLC dynamics in natural variants and mutants throughout winter in three field sites, we find that autumnal FLC expression, rather than epigenetic silencing, is the major variable conferred by the distinct Arabidopsis FLChaplotypes. This variation influences flowering responses of Arabidopsis accessions resulting in an interplay between promotion and delay of flowering in different climates to balance survival and, through a post-vernalization effect, reproductive output. These data reveal how expression variation through non-coding cis variation at FLC has enabled Arabidopsis accessions to adapt to different climatic conditions and year-on-year fluctuations.

Populations Are Differentiated in Biological Rhythms without Explicit Elevational Clines in the Plant Mimulus laciniatus

Environmental variation along an elevational gradient can yield phenotypic differentiation resulting from varying selection pressures on plant traits related to seasonal responses. Thus, genetic clines can evolve in a suite of traits, including the circadian clock, that drives daily cycling in varied traits and that shares its genetic background with adaptation to seasonality. We used populations of annual Mimulus laciniatus from different elevations in the Sierra Nevada in California to explore among-population differentiation in the circadian clock, flowering responses to photoperiod, and phenological traits (days to cotyledon emergence, days to flowering, and days to seed ripening) in controlled common-garden conditions. Further, we examined correlations of these traits with environmental variables related to temperature and precipitation. We observed that the circadian period in leaf movement was differentiated among populations sampled within about 100 km, with population means varying by 1.6 h. Significant local genetic variation occurred within 2 populations in which circadian period among families varied by up to 1.8 h. Replicated treatments with variable ecologically relevant photoperiods revealed marked population differentiation in critical day length for flowering that ranged from 11.0 to 14.1 h, corresponding to the time period between late February and mid-May in the wild. Flowering time varied among populations in a 14-h photoperiod. Regardless of this substantial population-level diversity, obvious linear clinality in trait variability across elevations could not be determined based on our genotypic sample; it is possible that more complex spatial patterns of variation arise in complex terrains such as those in the Sierra Nevada. Moreover, we did not find statistically significant bivariate correlations between population means of different traits. Our research contributes to the understanding of genetic variation in the circadian clock and in seasonal responses in natural populations, highlighting the need for more comprehensive investigations on the association between the clock and other adaptive traits in plants.

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.

Plasticity to drought and ecotypic differentiation in populations of a crop wild relative

Populations of widely distributed species often exhibit geographic variation in functional traits in response to environmental heterogeneity. Such trait variation may be the result of different adaptive mechanisms, including genetically based differentiation, phenotypic plasticity or a combination of both. Disentangling the genetic and environmental components of trait variation may be particularly interesting in crop wild relatives, since they may provide unique reservoirs of genetic diversity for crop improvement. In this study, we assessed ecotypic differentiation and patterns of plasticity to drought in populations of Lupinus angustifolius, a Mediterranean crop wild relative, from two climatically distinct regions in the Iberian Peninsula. Using an outdoor common garden, we compared phenotypic responses of inbred maternal families to two ecologically meaningful water availability treatments (drought and high-moisture). We measured 18 different functional traits related to growth, morphology, phenology and reproduction. Plants in the drought treatment grew less, had lower leaf chlorophyll content and photochemical efficiency, but also reproduced faster, produced larger seeds and altered leaflet morphology through increased leaflet thickness, higher leaflet dry matter content and lower specific leaf area. We also found significant differences between regions that likely reflect adaptation to climatically distinct environments, with populations from the south showing a faster onset of reproduction, higher leaf thickness and higher seed size, consistent with the drier conditions experienced in southern sites. Plasticity to drought was in most cases in the same direction as quantitative genetic differentiation (i.e. cogradient variation), providing evidence of the adaptive value of the plastic change. Our results show that both genetic differentiation and plasticity can generate adaptive phenotypic variation in L. angustifolius, and help to identify potentially valuable genetic resources to incorporate into breeding programmes.

Phenotypic selection under two contrasting environments in wild sunflower and its crop-wild hybrid

Hybridization is a common phenomenon in plants and can lead to the introgression of alleles from one population into another, generate new hybrid lineages, or cause species extinction. The environmental conditions and the genetic background of the participating populations may influence these outcomes since they can affect the fitness of hybrids, thereby increasing or decreasing the chances of introgression. Thus, it is important to understand the context-dependent prospects for introgression of alleles into diverse populations and under multiple ecological environments. Crop-wild hybridization presents an opportunity to explore these dynamics in agroecosystems. To this end, we used diverse wild and hybrid sunflowers from across the northern United States as a basis for evaluating variation in morphological traits and assessing context-dependent selection. These crop-wild hybrids and their wild counterparts were grown under agricultural conditions in the field with and without wheat competition. Interactions between origin and cross type affected expression of early functional traits, while interactions between competition and cross type acted on reproductive traits. A smattering of early and reproductive traits was affected by interactions between cross type and competition that varied by origin (i.e., 3-way interactions). Seven functional traits, especially number of branches and tertiary head diameter, underwent net and direct directional selection, while six out of these seven traits appear to also be experiencing nonlinear selection dynamics. In general, wild-like traits were favored under both sets of conditions, while, under wheat competition, some crop-like traits related to fast growth and primary head diameter became important. These data reaffirm the hypothesis that stressful conditions establish a scenario more suitable for crop introgression and clarify that nonlinear selection dynamics may play a role in this process.

Large-effect flowering time mutations reveal conditionally adaptive paths through fitness landscapes in Arabidopsis thaliana

Mutations are often assumed to be largely detrimental to fitness, but they may also be beneficial, and mutations with large phenotypic effects can persist in nature. One explanation for these observations is that mutations may be beneficial in specific environments because these conditions shift trait expression toward higher fitness. This hypothesis is rarely tested due to the difficulty of replicating mutants in multiple natural environments and measuring their phenotypes. We did so by planting Arabidopsis thaliana genotypes with large-effect flowering time mutations in field sites across the species’ European climate range. We quantified the adaptive value of mutant traits, finding that certain mutations increased fitness in some environments but not in others.

Contrary to previous assumptions that most mutations are deleterious, there is increasing evidence for persistence of large-effect mutations in natural populations. A possible explanation for these observations is that mutant phenotypes and fitness may depend upon the specific environmental conditions to which a mutant is exposed. Here, we tested this hypothesis by growing large-effect flowering time mutants of Arabidopsis thaliana in multiple field sites and seasons to quantify their fitness effects in realistic natural conditions. By constructing environment-specific fitness landscapes based on flowering time and branching architecture, we observed that a subset of mutations increased fitness, but only in specific environments. These mutations increased fitness via different paths: through shifting flowering time, branching, or both. Branching was under stronger selection, but flowering time was more genetically variable, pointing to the importance of indirect selection on mutations through their pleiotropic effects on multiple phenotypes. Finally, mutations in hub genes with greater connectedness in their regulatory networks had greater effects on both phenotypes and fitness. Together, these findings indicate that large-effect mutations may persist in populations because they influence traits that are adaptive only under specific environmental conditions. Understanding their evolutionary dynamics therefore requires measuring their effects in multiple natural environments.

Population climatic history predicts phenotypic responses in novel environments for Arabidopsis thaliana in North America

PREMISE

Determining how species perform in novel climatic environments is essential for understanding (1) responses to climate change and (2) evolutionary consequences of biological invasions. For the vast majority of species, the number of population characteristics that will predict performance and patterns of natural selection in novel locations in the wild remains limited.

METHODS

We evaluated phenological, vegetative, architectural, and fitness-related traits in experimental gardens in contrasting climates (Ontario, Canada, and South Carolina, USA) in the North American non-native distribution of Arabidopsis thaliana. We assessed the effects of climatic distance, geographic distance, and genetic features of history on performance and patterns of natural selection in the novel garden settings.

RESULTS

We found that plants had greater survivorship, flowered earlier, were larger, and produced more fruit in the south, and that genotype-by-environment interactions were significant between gardens. However, our analyses revealed similar patterns of natural selection between gardens in distinct climate zones. After accounting for genetic ancestry, we also detected that population climatic distance best predicted performance within gardens.

CONCLUSIONS

These data suggest that colonization success in novel, non-native environments is determined by a combination of climate and genetic history. When performance at novel sites was assessed with seed sources from geographically and genetically disparate, established non-native populations, proximity to the garden alone was insufficient to predict performance. Our study highlights the need to evaluate seed sources from diverse origins to describe comprehensively phenotypic responses to novel environments, particularly for taxa in which many source populations may contribute to colonization.

The fitness benefits of genetic variation in circadian clock regulation

Functional circadian clocks are essential for fitness in diverse ecosystems, facilitating detection of predictable light-dark and temperature cycles. The molecular basis of endogenous clocks is variable across the tree of life, but it has one omnipresent attribute: natural genetic diversity that manifests as variation for instance in circadian period length around the hypothesised optimum of 24 hours. Latitudinal variation in photoperiod alone is unlikely to account for the vast diversity documented in varied organisms, but we have yet to achieve a solid understanding of the interplay between clock variability and natural selection. Recent circadian studies sampling populations have drawn attention to the hierarchical structure of genetic diversity in the wild, unveiling pronounced genetic variation even on a scale of metres.

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.

Linking genes with ecological strategies in Arabidopsis thaliana

Arabidopsis thaliana is the most prominent model system in plant molecular biology and genetics. Although its ecology was initially neglected, collections of various genotypes revealed a complex population structure, with high levels of genetic diversity and substantial levels of phenotypic variation. This helped identify the genes and gene pathways mediating phenotypic change. Population genetics studies further demonstrated that this variation generally contributes to local adaptation. Here, we review evidence showing that traits affecting plant life history, growth rate, and stress reactions are not only locally adapted, they also often co-vary. Co-variation between these traits indicates that they evolve as trait syndromes, and reveals the ecological diversification that took place within A. thaliana. We argue that examining traits and the gene that control them within the context of global summary schemes that describe major ecological strategies will contribute to resolve important questions in both molecular biology and ecology.

Variation in wild pea (Pisum sativum subsp. elatius) seed dormancy and its relationship to the environment and seed coat traits

BACKGROUND

Seed germination is one of the earliest key events in the plant life cycle. The timing of transition from seed to seedling is an important developmental stage determining the survival of individuals that influences the status of populations and species. Because of wide geographical distribution and occurrence in diverse habitats, wild pea (Pisum sativum subsp. elatius) offers an excellent model to study physical type of seed dormancy in an ecological context. This study addresses the gap in knowledge of association between the seed dormancy, seed properties and environmental factors, experimentally testing oscillating temperature as dormancy release clue.

METHODS

Seeds of 97 pea accessions were subjected to two germination treatments (oscillating temperatures of 25/15 °C and 35/15 °C) over 28 days. Germination pattern was described using B-spline coefficients that aggregate both final germination and germination speed. Relationships between germination pattern and environmental conditions at the site of origin (soil and bioclimatic variables extracted from WorldClim 2.0 and SoilGrids databases) were studied using principal component analysis, redundancy analysis and ecological niche modelling. Seeds were analyzed for the seed coat thickness, seed morphology, weight and content of proanthocyanidins (PA).

RESULTS

Seed total germination ranged from 0% to 100%. Cluster analysis of germination patterns of seeds under two temperature treatments differentiated the accessions into three groups: (1) non-dormant (28 accessions, mean germination of 92%), (2) dormant at both treatments (29 acc., 15%) and (3) responsive to increasing temperature range (41 acc., with germination change from 15 to 80%). Seed coat thickness differed between groups with dormant and responsive accessions having thicker testa (median 138 and 140 µm) than non-dormant ones (median 84 mm). The total PA content showed to be higher in the seed coat of dormant (mean 2.18 mg g-1) than those of non-dormant (mean 1.77 mg g-1) and responsive accessions (mean 1.87 mg g-1). Each soil and bioclimatic variable and also germination responsivity (representing synthetic variable characterizing germination pattern of seeds) was spatially clustered. However, only one environmental variable (BIO7, i.e., annual temperature range) was significantly related to germination responsivity. Non-dormant and responsive accessions covered almost whole range of BIO7 while dormant accessions are found in the environment with higher annual temperature, smaller temperature variation, seasonality and milder winter. Ecological niche modelling showed a more localized potential distribution of dormant group. Seed dormancy in the wild pea might be part of a bet-hedging mechanism for areas of the Mediterranean basin with more unpredictable water availability in an otherwise seasonal environment. This study provides the framework for analysis of environmental aspects of physical seed dormancy.

Variation in wild pea (Pisum sativum subsp. elatius) seed dormancy and its relationship to the environment and seed coat traits

BACKGROUND

Seed germination is one of the earliest key events in the plant life cycle. The timing of transition from seed to seedling is an important developmental stage determining the survival of individuals that influences the status of populations and species. Because of wide geographical distribution and occurrence in diverse habitats, wild pea (Pisum sativum subsp. elatius) offers an excellent model to study physical type of seed dormancy in an ecological context. This study addresses the gap in knowledge of association between the seed dormancy, seed properties and environmental factors, experimentally testing oscillating temperature as dormancy release clue.

METHODS

Seeds of 97 pea accessions were subjected to two germination treatments (oscillating temperatures of 25/15 °C and 35/15 °C) over 28 days. Germination pattern was described using B-spline coefficients that aggregate both final germination and germination speed. Relationships between germination pattern and environmental conditions at the site of origin (soil and bioclimatic variables extracted from WorldClim 2.0 and SoilGrids databases) were studied using principal component analysis, redundancy analysis and ecological niche modelling. Seeds were analyzed for the seed coat thickness, seed morphology, weight and content of proanthocyanidins (PA).

RESULTS

Seed total germination ranged from 0% to 100%. Cluster analysis of germination patterns of seeds under two temperature treatments differentiated the accessions into three groups: (1) non-dormant (28 accessions, mean germination of 92%), (2) dormant at both treatments (29 acc., 15%) and (3) responsive to increasing temperature range (41 acc., with germination change from 15 to 80%). Seed coat thickness differed between groups with dormant and responsive accessions having thicker testa (median 138 and 140 µm) than non-dormant ones (median 84 mm). The total PA content showed to be higher in the seed coat of dormant (mean 2.18 mg g^-1) than those of non-dormant (mean 1.77 mg g^-1) and responsive accessions (mean 1.87 mg g^-1). Each soil and bioclimatic variable and also germination responsivity (representing synthetic variable characterizing germination pattern of seeds) was spatially clustered. However, only one environmental variable (BIO7, i.e., annual temperature range) was significantly related to germination responsivity. Non-dormant and responsive accessions covered almost whole range of BIO7 while dormant accessions are found in the environment with higher annual temperature, smaller temperature variation, seasonality and milder winter. Ecological niche modelling showed a more localized potential distribution of dormant group. Seed dormancy in the wild pea might be part of a bet-hedging mechanism for areas of the Mediterranean basin with more unpredictable water availability in an otherwise seasonal environment. This study provides the framework for analysis of environmental aspects of physical seed dormancy.

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.