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

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.

Host developmental stages shape the evolution of a plant RNA virus

Viruses are obligate pathogens that entirely rely on their hosts to complete their infectious cycle. The outcome of viral infections depends on the status of the host. Host developmental stage is an important but sometimes overlooked factor impacting host-virus interactions. This impact is especially relevant in a context where climate change and human activities are altering plant development. To better understand how different host developmental stages shape virus evolution, we experimentally evolved turnip mosaic virus (TuMV) on Arabidopsis thaliana at three different developmental stages: vegetative (juvenile), bolting (transition) and reproductive (mature). After infecting plants with an Arabidopsis-naive or an Arabidopsis-well-adapted TuMV isolate, we observed that hosts in later developmental stages were prone to faster and more severe infections. This observation was extended to viruses belonging to different genera. Thereafter, we experimentally evolved lineages of the naive and the well-adapted TuMV isolates in plants from each of the three developmental stages. All evolved viruses enhanced their infection traits, but this increase was more intense in viruses evolved in younger hosts. The genomic changes of the evolved viral lineages revealed mutation patterns that strongly depended on the founder viral isolate as well as on the developmental stage of the host wherein the lineages were evolved. This article is part of the theme issue 'Infectious disease ecology and evolution in a changing world'.

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.

Environmental patterns of adaptation after range expansion in Leontodon longirostris: The effect of phenological events on fitness-related traits

PREMISE

Because of expected range shifts associated with climate change, there is a renewed interest in the evolutionary factors constraining adaptation, among which are genetic bottlenecks, drift, and increased mutational load after range expansion. Here we study adaptation in the short-lived species Leontodon longirostris showing reduced genetic diversity and increased genetic load along an expansion route.

METHODS

We assessed the phenological patterns of variation, and their effect on fitness-related traits, on 42 L. longirostris populations and six populations of the sister taxa L. saxatilis in a common garden located within the current range of both species. The comparison among L. longirostris populations allowed us to test for genetic clines consistent with local adaptation, whereas the comparison between taxa provided evidence for common adaptive features at the species level.

RESULTS

We found significant within-species variability for most traits, as well as differences with its close relative L. saxatilis. In general, seeds from drier, warmer, and unpredictable habitats showed overall lower and more restricted conditions for germination, seedlings emerged later and plants flowered earlier. Consequently, genotypes from arid and unpredictable environments attained smaller reproductive sizes and allocated more biomass to reproduction. Flowering time had the strongest direct effect on total plant size, but seedling emergence also showed an important indirect effect.

CONCLUSIONS

Our results show the crucial role of phenological patterns in shaping adaptive clines for major life-history stage transitions. Furthermore, the genetic load observed in L. longirostris does not seem to preclude adaptation to the climatic variability encountered along the expansion route.

Local adaptation and spatiotemporal patterns of genetic diversity revealed by repeated sampling of Caenorhabditis elegans across the Hawaiian Islands

The nematode Caenorhabditis elegans is among the most widely studied organisms, but relatively little is known about its natural ecology. Genetic diversity is low across much of the globe but high in the Hawaiian Islands and across the Pacific Rim. To characterize the niche and genetic diversity of C. elegans on the Hawaiian Islands and to explore how genetic diversity might be influenced by local adaptation, we repeatedly sampled nematodes over a three‐year period, measured various environmental parameters at each sampling site, and whole‐genome sequenced the C. elegans isolates that we identified. We found that the typical Hawaiian C. elegans niche comprises moderately moist native forests at high elevations (500–1,500 m) where ambient air temperatures are cool (15–20°C). Compared to other Caenorhabditis species found on the Hawaiian Islands (e.g., Caenorhabditis briggsae and Caenorhabditis tropicalis), we found that C. elegans were enriched in native habitats. We measured levels of genetic diversity and differentiation among Hawaiian C. elegans and found evidence of seven genetically distinct groups distributed across the islands. Then, we scanned these genomes for signatures of local adaptation and identified 18 distinct regions that overlap with hyper‐divergent regions, which may be maintained by balancing selection and are enriched for genes related to environmental sensing, xenobiotic detoxification, and pathogen resistance. These results provide strong evidence of local adaptation among Hawaiian C. elegans and contribute to our understanding of the forces that shape genetic diversity on the most remote volcanic archipelago in the world.

A role of flowering genes in the tolerance of Arabidopsis thaliana to cucumber mosaic virus

The genetic basis of plant tolerance to parasites is poorly understood. We have previously shown that tolerance of Arabidopsis thaliana to its pathogen cucumber mosaic virus is achieved through changes in host life-history traits on infection that result in delaying flowering and reallocating resources from vegetative growth to reproduction. In this system we analyse here genetic determinants of tolerance using a recombinant inbred line family derived from a cross of two accessions with extreme phenotypes. Three major quantitative trait loci for tolerance were identified, which co-located with three flowering repressor genes, FLC, FRI, and HUA2. The role of these genes in tolerance was further examined in genotypes carrying functional or nonfunctional alleles. Functional alleles of FLC together with FRI and/or HUA2 were required for both tolerance and resource reallocation from growth to reproduction. Analyses of FLC alleles from wild accessions that differentially modulate flowering time showed that they ranked differently for their effects on tolerance and flowering. These results pinpoint a role of FLC in A. thaliana tolerance to cucmber mosaic virus, which is a novel major finding, as FLC has not been recognized previously to be involved in plant defence. Although tolerance is associated with a delay in flowering that allows resource reallocation, our results indicate that FLC regulates tolerance and flowering initiation by different mechanisms. Thus, we open a new avenue of research on the interplay between defence and development in plants.

Cold stress and freezing tolerance negatively affect the fitness of Arabidopsis thaliana accessions under field and controlled conditions

Higher acclimated freezing tolerance improved winter survival, but reduced reproductive fitness of Arabidopsis thaliana accessions under field and controlled conditions. Low temperature is one of the most important abiotic factors influencing plant fitness and geographical distribution. In addition, cold stress is known to influence crop yield and is therefore of great economic importance. Increased freezing tolerance can be acquired by the process of cold acclimation, but this may be associated with a fitness cost. To assess the influence of cold stress on the fitness of plants, long-term field trials over 5 years were performed with six natural accessions of Arabidopsis thaliana ranging from very tolerant to very sensitive to freezing. Fitness parameters, as seed yield and 1000 seed mass, were measured and correlation analyses with temperature and freezing tolerance data performed. The results were compared with fitness parameters from controlled chamber experiments over 3 years with application of cold priming and triggering conditions. Winter survival and seed yield per plant were positively correlated with temperature in field experiments. In addition, winter survival and 1000 seed mass were correlated with the cold-acclimated freezing tolerance of the selected Arabidopsis accessions. The results provide strong evidence for a trade-off between higher freezing tolerance and reproductive fitness in A. thaliana, which might have ecological impacts in the context of global warming.

Host population structure for tolerance determines the evolution of plant-virus interactions

Heterogeneity for plant defences determines both the capacity of host populations to buffer the effect of infection and the pathogen´s fitness. However, little information is known on how host population structure for tolerance, a major plant defence, impacts the evolution of plant-pathogen interactions. By performing 10 serial passages of Turnip mosaic virus (TuMV) in Arabidopsis thaliana populations with varying proportion of tolerant genotypes simulating different structures for this trait, we analysed how host heterogeneity for this defence shapes the evolution of both virus multiplication, the effect of infection on plant fecundity and mortality, and plant tolerance and resistance. Results indicated that a higher proportion of tolerant genotypes in the host population promotes virus multiplication and reduces the effect of infection on plant mortality, but not on plant fecundity. These changes resulted in more effective plant tolerance to virus infection. Conversely, a lower proportion of tolerant genotypes reduced virus multiplication, boosting plant resistance. Our work for the first time provides evidence of the main role of host population structure for tolerance on pathogen evolution and on the subsequent feedback loops on plant defences.

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.

Arabidopsis thaliana Genes Associated with Cucumber mosaic virus Virulence and Their Link to Virus Seed Transmission

Virulence, the effect of pathogen infection on progeny production, is a major determinant of host and pathogen fitness as it affects host fecundity and pathogen transmission. In plant-virus interactions, ample evidence indicates that virulence is genetically controlled by both partners. However, the host genetic determinants are poorly understood. Through a genome-wide association study (GWAS) of 154 Arabidopsis thaliana genotypes infected by Cucumber mosaic virus (CMV), we identified eight host genes associated with virulence, most of them involved in response to biotic stresses and in cell wall biogenesis in plant reproductive structures. Given that virulence is a main determinant of the efficiency of plant virus seed transmission, we explored the link between this trait and the genetic regulation of virulence. Our results suggest that the same functions that control virulence are also important for CMV transmission through seeds. In sum, this work provides evidence of a novel role for some previously known plant defense genes and for the cell wall metabolism in plant virus interactions.

Natural selection on traits and trait plasticity in Arabidopsis thaliana varies across competitive environments

Interspecific competition reduces resource availability and can affect evolution. We quantified multivariate selection in the presence and absence of strong interspecific competition using a greenhouse experiment with 35 natural accessions of Arabidopsis thaliana. We assessed selection on nine traits representing plant phenology, growth, and architecture, as well as their plasticities. Competition reduced biomass and fitness by over 98%, and plastic responses to competition varied by genotype (significant G × E) for all traits except specific leaf area (SLA). Competitive treatments altered selection on flowering phenology and plant architecture, with significant selection on all phenology traits and most architecture traits under competition-present conditions but little indication that selection occurred in the absence of competitors. Plasticity affected fitness only in competition-present conditions, where plasticity in flowering time and early internode lengths was adaptive. The competitive environment caused changes in the trait correlation structure and surprisingly reduced phenotypic integration, which helped explain some of the observed selection patterns. Despite this overall shift in the trait correlation matrix, genotypes with delayed flowering had lower SLA (thicker, tougher leaves) regardless of the competitive environment, a pattern we have not seen previously reported in the literature. Overall, our study highlights multiple ways in which interspecific competition can alter selective regimes, contributing to our understanding of variability in selection processes over space and time.

Reproductive Assurance Maintains Red-Flowered Plants of Lysimachia arvensis in Mediterranean Populations Despite Inbreeding Depression

Flower color polymorphism, an infrequent but phylogenetically widespread condition in plants, is captivating because it can only be maintained under a few selective regimes but also because it can drive intra-morph assortative mating and promote speciation. Lysimachia arvensis is a polymorphic species with red or blue flowered morphs. In polymorphic populations, which are mostly Mediterranean, pollinators prefer blue-flowered plants to the red ones, and abiotic factors also favors blue-flowered plants. We hypothesize that the red morph is maintained in Mediterranean areas due to its selfing capacity. We assessed inbreeding depression in both color morphs in two Mediterranean populations and genetic diversity was studied via SSR microsatellites in 20 natural populations. Results showed that only 44-47% of selfed progeny of the red plants reached reproduction while about 72-91% of blue morph progeny did it. Between-morph genetic differentiation was high and the red morph had a lower genetic diversity and a higher inbreeding coefficient, mainly in the Mediterranean. Results suggest that selfing maintaining the red morph in Mediterranean areas despite its inbreeding depression. In addition, genetic differentiation between morphs suggests a low gene flow between them, suggesting reproductive isolation.

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.

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.

Trade-offs between host tolerances to different pathogens in plant-virus interactions

Although accumulating evidence indicates that tolerance is a plant defence strategy against pathogens as widespread as resistance, how plants evolve tolerance is poorly understood. Theory predicts that hosts will evolve to maximize tolerance or resistance, but not both. Remarkably, most experimental works failed in finding this trade-off. We tested the hypothesis that the evolution of tolerance to one virus is traded-off against tolerance to others, rather than against resistance and identified the associated mechanisms. To do so, we challenged eighteen Arabidopsis thaliana genotypes with Turnip mosaic virus (TuMV) and Cucumber mosaic virus (CMV). We characterized plant life-history trait modifications associated with reduced effects of TuMV and CMV on plant seed production (fecundity tolerance) and life period (mortality tolerance), both measured as a norm of reaction across viral loads (range tolerance). Also, we analysed resistance-tolerance and tolerance-tolerance trade-offs. Results indicate that tolerance to TuMV is associated with changes in the length of the pre-reproductive and reproductive periods, and tolerance to CMV with resource reallocation from growth to reproduction; and that tolerance to TuMV is traded-off against tolerance to CMV in a virulence-dependent manner. Thus, this work provides novel insights on the mechanisms of plant tolerance and highlights the importance of considering the combined effect of different pathogens to understand how plant defences evolve.

Natural selection on the Arabidopsis thaliana genome in present and future climates

Through the lens of evolution, climate change is an agent of natural selection that forces populations to change and adapt, or face extinction. However, current assessments of the risk of biodiversity associated with climate change¹ do not typically take into account how natural selection influences populations differently depending on their genetic makeup². Here we make use of the extensive genome information that is available for Arabidopsis thaliana and measure how manipulation of the amount of rainfall affected the fitness of 517 natural Arabidopsis lines that were grown in Spain and Germany. This allowed us to directly infer selection along the genome³. Natural selection was particularly strong in the hot-dry location in Spain, where 63% of lines were killed and where natural selection substantially changed the frequency of approximately 5% of all genome-wide variants. A significant portion of this climate-driven natural selection of variants was predictable from signatures of local adaptation (R² = 29-52%), as genetic variants that were found in geographical areas with climates more similar to the experimental sites were positively selected. Field-validated predictions across the species range indicated that Mediterranean and western Siberian populations-at the edges of the environmental limits of this species-currently experience the strongest climate-driven selection. With more frequent droughts and rising temperatures in Europe⁴, we forecast an increase in directional natural selection moving northwards from the southern end of Europe, putting many native A. thaliana populations at evolutionary risk.

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.

Cucumber mosaic virus infection as a potential selective pressure on Arabidopsis thaliana populations

It has been proposed that in wild ecosystems viruses are often plant mutualists, whereas agroecosystems favour pathogenicity. We seek evidence for virus pathogenicity in wild ecosystems through the analysis of plant-virus coevolution, which requires a negative effect of infection on the host fitness. We focus on the interaction between Arabidopsis thaliana and Cucumber mosaic virus (CMV), which is significant in nature. We studied the genetic diversity of A. thaliana for two defence traits, resistance and tolerance, to CMV. A set of 185 individuals collected in 76 A. thaliana Iberian wild populations were inoculated with different CMV strains. Resistance was estimated from the level of virus multiplication in infected plants, and tolerance from the effect of infection on host progeny production. Resistance and tolerance to CMV showed substantial genetic variation within and between host populations, and depended on the virus x host genotype interaction, two conditions for coevolution. Resistance and tolerance were co-occurring independent traits that have evolved independently from related life-history traits involved in adaptation to climate. The comparison of the genetic structure for resistance and tolerance with that for neutral traits (QST/FST analyses) indicated that both defence traits are likely under uniform selection. These results strongly suggest that CMV infection selects for defence on A. thaliana populations, and support plant-virus coevolution. Thus, we propose that CMV infection reduces host fitness under the field conditions of the wild A. thaliana populations studied.

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.

Genetic association with high-resolution climate data reveals selection footprints in the genomes of barley landraces across the Iberian Peninsula

Landraces are local populations of crop plants adapted to a particular environment. Extant landraces are surviving genetic archives, keeping signatures of the selection processes experienced by them until settling in their current niches. This study intends to establish relationships between genetic diversity of barley (Hordeum vulgare L.) landraces collected in Spain and the climate of their collection sites. A high-resolution climatic data set (5 × 5 km spatial, 1-day temporal grid) was computed from over 2,000 temperature and 7,000 precipitation stations across peninsular Spain. This data set, spanning the period 1981-2010, was used to derive agroclimatic variables meaningful for cereal production at the collection sites of 135 barley landraces. Variables summarize temperature, precipitation, evapotranspiration, potential vernalization and frost probability at different times of the year and time scales (season and month). SNP genotyping of the landraces was carried out combining Illumina Infinium assays and genotyping-by-sequencing, yielding 9,920 biallelic markers (7,479 with position on the barley reference genome). The association of these SNPs with agroclimatic variables was analysed at two levels of genetic diversity, with and without taking into account population structure. The whole data sets and analysis pipelines are documented and available at https://eead-csic-compbio.github.io/barley-agroclimatic-association. We found differential adaptation of the germplasm groups identified to be dominated by reactions to cold temperature and late-season frost occurrence, as well as to water availability. Several significant associations pointing at specific adaptations to agroclimatic features related to temperature and water availability were observed, and candidate genes underlying some of the main regions are proposed.

Adaptation genomics: The angel is in the details

This article comments on: Genome-wide signatures of flowering adaptation to climate temperature: Regional analyses in a highly diverse native range of Arabidopsis thaliana.

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.

Genome-wide signatures of flowering adaptation to climate temperature: Regional analyses in a highly diverse native range of Arabidopsis thaliana

Current global change is fueling an interest to understand the genetic and molecular mechanisms of plant adaptation to climate. In particular, altered flowering time is a common strategy for escape from unfavourable climate temperature. In order to determine the genomic bases underlying flowering time adaptation to this climatic factor, we have systematically analysed a collection of 174 highly diverse Arabidopsis thaliana accessions from the Iberian Peninsula. Analyses of 1.88 million single nucleotide polymorphisms provide evidence for a spatially heterogeneous contribution of demographic and adaptive processes to geographic patterns of genetic variation. Mountains appear to be allele dispersal barriers, whereas the relationship between flowering time and temperature depended on the precise temperature range. Environmental genome-wide associations supported an overall genome adaptation to temperature, with 9.4% of the genes showing significant associations. Furthermore, phenotypic genome-wide associations provided a catalogue of candidate genes underlying flowering time variation. Finally, comparison of environmental and phenotypic genome-wide associations identified known (Twin Sister of FT, FRIGIDA-like 1, and Casein Kinase II Beta chain 1) and new (Epithiospecifer Modifier 1 and Voltage-Dependent Anion Channel 5) genes as candidates for adaptation to climate temperature by altered flowering time. Thus, this regional collection provides an excellent resource to address the spatial complexity of climate adaptation in annual plants.

Phenotypic plasticity in response to temperature fluctuations is genetically variable, and relates to climatic variability of origin, in Arabidopsis thaliana

Under current climate change, increasing mean temperatures are not only causing hotter summers, but temperature variability is increasing as well. Phenotypic plasticity can help plants to overcome negative effects of temperature variability and allow them to rapidly adjust traits to adverse conditions. Moreover, genetic variation in such plasticity could provide potential for adaptive evolution in response to changing climate variability. Here, we conducted an experiment with 11 Arabidopsis thaliana genotypes to investigate intraspecific variation in plant responses to two aspects of variable temperature stress: timing and frequency. We found that the timing but not frequency of temperature stress affected the phenology, growth, reproduction and allocation strategy of plants, and that genotypes differed substantially in their responses. Moreover, trait plasticity was positively related to precipitation variability of origin, suggesting an adaptive role of plasticity. Our results indicate that the developmental stage of a plant during heat stress is a key determinant of its response, and that plasticity to temperature variability is an evolving and possibly adaptive trait in natural populations of A. thaliana. More generally, our study demonstrates the usefulness of studying plant responses to climatic variability per se, given that climatic variability is predicted to increase in the future.

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.

Temperature fine-tunes Mediterranean Arabidopsis thaliana life-cycle phenology geographically

To understand how adaptive evolution in life-cycle phenology operates in plants, we need to unravel the effects of geographic variation in putative agents of natural selection on life-cycle phenology by considering all key developmental transitions and their co-variation patterns. We address this goal by quantifying the temperature-driven and geographically varying relationship between seed dormancy and flowering time in the annual Arabidopsis thaliana across the Iberian Peninsula. We used data on genetic variation in two major life-cycle traits, seed dormancy (DSDS50) and flowering time (FT), in a collection of 300 A. thaliana accessions from the Iberian Peninsula. The geographically varying relationship between life-cycle traits and minimum temperature, a major driver of variation in DSDS50 and FT, was explored with geographically weighted regressions (GWR). The environmentally varying correlation between DSDS50 and FT was analysed by means of sliding window analysis across a minimum temperature gradient. Maximum local adjustments between minimum temperature and life-cycle traits were obtained in the southwest Iberian Peninsula, an area with the highest minimum temperatures. In contrast, in off-southwest locations, the effects of minimum temperature on DSDS50 were rather constant across the region, whereas those of minimum temperature on FT were more variable, with peaks of strong local adjustments of GWR models in central and northwest Spain. Sliding window analysis identified a minimum temperature turning point in the relationship between DSDS50 and FT around a minimum temperature of 7.2 °C. Above this minimum temperature turning point, the variation in the FT/DSDS50 ratio became rapidly constrained and the negative correlation between FT and DSDS50 did not increase any further with increasing minimum temperatures. The southwest Iberian Peninsula emerges as an area where variation in life-cycle phenology appears to be restricted by the duration and severity of the hot summer drought. The temperature-driven varying relationship between DSDS50 and FT detected environmental boundaries for the co-evolution between FT and DSDS50 in A. thaliana. In the context of global warming, we conclude that A. thaliana phenology from the southwest Iberian Peninsula, determined by early flowering and deep seed dormancy, might become the most common life-cycle phenotype for this annual plant in the region.

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

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

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.

Habitat Temperature and Precipitation of Arabidopsis thaliana Ecotypes Determine the Response of Foliar Vasculature, Photosynthesis, and Transpiration to Growth Temperature

Acclimatory adjustments of foliar vascular architecture, photosynthetic capacity, and transpiration rate in Arabidopsis thaliana ecotypes (Italian, Polish [Col-0], Swedish) were characterized in the context of habitat of origin. Temperatures of the habitat of origin decreased linearly with increasing habitat latitude, but habitat precipitation was greatest in Italy, lowest in Poland, and intermediate in Sweden. Plants of the three ecotypes raised under three different growth temperature regimes (low, moderate, and high) exhibited highest photosynthetic capacities, greatest leaf thickness, highest chlorophyll a/b ratio and levels of β-carotene, and greatest levels of wall ingrowths in phloem transfer cells, and, in the Col-0 and Swedish ecotypes, of phloem per minor vein in plants grown at the low temperature. In contrast, vein density and minor vein tracheary to sieve element ratio increased with increasing growth temperature - most strongly in Col-0 and least strongly in the Italian ecotype - and transpirational water loss correlated with vein density and number of tracheary elements per minor vein. Plotting of these vascular features as functions of climatic conditions in the habitat of origin suggested that temperatures during the evolutionary history of the ecotypes determined acclimatory responses of the foliar phloem and photosynthesis to temperature in this winter annual that upregulates photosynthesis in response to lower temperature, whereas the precipitation experienced during the evolutionary history of the ecotypes determined adjustment of foliar vein density, xylem, and transpiration to temperature. In particular, whereas photosynthetic capacity, leaf thickness, and foliar minor vein phloem features increased linearly with increasing latitude and decreasing temperature of the habitats of origin in response to experimental growth at low temperature, transpiration rate, foliar vein density, and minor vein tracheary element numbers and cross-sectional areas increased linearly with decreasing precipitation level in the habitats of origin in response to experimental growth at high temperature. This represents a situation where temperature acclimation of the apparent capacity for water flux through the xylem and transpiration rate in a winter annual responded differently from that of photosynthetic capacity, in contrast to previous reports of strong relationships between hydraulic conductance and photosynthesis in other studies.

Genomic and phenotypic differentiation of Arabidopsis thaliana along altitudinal gradients in the North Italian Alps

Altitudinal gradients in mountain regions are short-range clines of different environmental parameters such as temperature or radiation. We investigated genomic and phenotypic signatures of adaptation to such gradients in five Arabidopsis thaliana populations from the North Italian Alps that originated from 580 to 2350 m altitude by resequencing pools of 19-29 individuals from each population. The sample includes two pairs of low- and high-altitude populations from two different valleys. High-altitude populations showed a lower nucleotide diversity and negative Tajima's D values and were more closely related to each other than to low-altitude populations from the same valley. Despite their close geographic proximity, demographic analysis revealed that low- and high-altitude populations split between 260 000 and 15 000 years before present. Single nucleotide polymorphisms whose allele frequencies were highly differentiated between low- and high-altitude populations identified genomic regions of up to 50 kb length where patterns of genetic diversity are consistent with signatures of local selective sweeps. These regions harbour multiple genes involved in stress response. Variation among populations in two putative adaptive phenotypic traits, frost tolerance and response to light/UV stress was not correlated with altitude. Taken together, the spatial distribution of genetic diversity reflects a potentially adaptive differentiation between low- and high-altitude populations, whereas the phenotypic differentiation in the two traits investigated does not. It may resemble an interaction between adaptation to the local microhabitat and demographic history influenced by historical glaciation cycles, recent seed dispersal and genetic drift in local populations.

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.

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.

The relationship between host lifespan and pathogen reservoir potential: an analysis in the system Arabidopsis thaliana--cucumber mosaic virus

Identification of the determinants of pathogen reservoir potential is central to understand disease emergence. It has been proposed that host lifespan is one such determinant: short-lived hosts will invest less in costly defenses against pathogens, so that they will be more susceptible to infection, more competent as sources of infection and/or will sustain larger vector populations, thus being effective reservoirs for the infection of long-lived hosts. This hypothesis is sustained by analyses of different hosts of multihost pathogens, but not of different genotypes of the same host species. Here we examined this hypothesis by comparing two genotypes of the plant Arabidopsis thaliana that differ largely both in life-span and in tolerance to its natural pathogen Cucumber mosaic virus (CMV). Experiments with the aphid vector Myzus persicae showed that both genotypes were similarly competent as sources for virus transmission, but the short-lived genotype was more susceptible to infection and was able to sustain larger vector populations. To explore how differences in defense against CMV and its vector relate to reservoir potential, we developed a model that was run for a set of experimentally-determined parameters, and for a realistic range of host plant and vector population densities. Model simulations showed that the less efficient defenses of the short-lived genotype resulted in higher reservoir potential, which in heterogeneous host populations may be balanced by the longer infectious period of the long-lived genotype. This balance was modulated by the demography of both host and vector populations, and by the genetic composition of the host population. Thus, within-species genetic diversity for lifespan and defenses against pathogens will result in polymorphisms for pathogen reservoir potential, which will condition within-population infection dynamics. These results are relevant for a better understanding of host-pathogen co-evolution, and of the dynamics of pathogen emergence.

Understanding pathogen emergence is a major goal of pathology, because of the high impact of emerging diseases. Pathogens emerge onto a new host from a reservoir, hence the relevance of identifying the determinants of host's reservoir potential. Host lifespan is considered as one such determinant: short-lived hosts will invest less in defenses, being more susceptible to infection, more competent as infection sources and/or will sustain larger vector populations, and thus, are effective reservoirs for long-lived host infection. Evidence for this hypothesis derives from analyses of different hosts of multihost pathogens, and here we examine whether it holds at the within-species level by comparing two genotypes of the plant Arabidopsis thaliana that differ in life-span and in tolerance to its natural pathogen Cucumber mosaic virus. Experiments showed that defenses to the virus and its aphid vector were less efficient in the short-lived genotype that, according to model simulations, was an effective reservoir under a large range of conditions. Reservoir potential, though, was modulated by the demography of host and vector and by the genetic composition of the host population. Thus, within-species genetic diversity for lifespan and pathogen defense will result in differences in reservoir potential, which will condition infection dynamics and host-pathogen co-evolution.