Open Access

The change in allele frequencies within a population over time represents a fundamental process of evolution. By monitoring allele frequencies, we can analyze the effects of natural selection and genetic drift on populations. To efficiently track time-resolved genetic change, large experimental or wild populations can be sequenced as pools of individuals sampled over time using high-throughput genome sequencing (called the Evolve & Resequence approach, E&R). Here, we present a set of experiments using hundreds of natural genotypes of the model plant Arabidopsis thaliana to showcase the power of this approach to study rapid evolution at large scale. First, we validate that sequencing DNA directly extracted from pools of flowers from multiple plants -- organs that are relatively consistent in size and easy to sample -- produces comparable results to other, more expensive state-of-the-art approaches such as sampling and sequencing of individual leaves. Sequencing pools of flowers from 25-50 individuals at ∼40X coverage recovers genome-wide frequencies in diverse populations with accuracy r > 0.95. Secondly, to enable analyses of evolutionary adaptation using E&R approaches of plants in highly replicated environments, we provide open source tools that streamline sequencing data curation and calculate various population genetic statistics two orders of magnitude faster than current software. To directly demonstrate the usefulness of our method, we conducted a two-year outdoor evolution experiment with A. thaliana to show signals of rapid evolution in multiple genomic regions. We demonstrate how these laboratory and computational Pool-seq-based methods can be scaled to study hundreds of populations across many climates.

Abstract Natural plant populations often harbour substantial heritable variation in DNA methylation. However, a thorough understanding of the genetic and environmental drivers of this epigenetic variation requires large-scale and high-resolution data, which currently exist only for a few model species. Here, we studied 207 lines of the annual weed Thlaspi arvense (field pennycress), collected across a large latitudinal gradient in Europe and propagated in a common environment. By screening for variation in DNA sequence and DNA methylation using whole-genome (bisulfite) sequencing, we found significant epigenetic population structure across Europe. Average levels of DNA methylation were strongly context-dependent, with highest DNA methylation in CG context, particularly in transposable elements and in intergenic regions. Residual DNA methylation variation within all contexts was associated with genetic variants, which often co-localized with annotated methylation machinery genes but also with new candidates. Variation in DNA methylation was also significantly associated with climate of origin, with methylation levels being higher in warmer regions and lower in more variable climates. Finally, we used variance decomposition to assess genetic versus environmental associations with differentially methylation regions (DMRs). We found that while genetic variation was generally the strongest predictor of DMRs, the strength of environmental associations increased from CG to CHG and CHH, with climate-of-origin as the strongest predictor in about one third of the CHH DMRs. In summary, our data show that natural epigenetic variation in Thlaspi arvense is significantly associated with both DNA sequence and environment of origin, and that the relative importance of the two factors strongly depends on the sequence context of DNA methylation. T. arvense is an emerging biofuel and winter cover crop; our results may hence be relevant for breeding efforts and agricultural practices in the context of rapidly changing environmental conditions. Author Summary: Variation within species is an important level of biodiversity, and it is key for future adaptation. Besides variation in DNA sequence, plants also harbour heritable variation in DNA methylation, and we want to understand the evolutionary significance of this epigenetic variation, in particular how much of it is under genetic control, and how much is associated with the environment. We addressed these questions in a high-resolution molecular analysis of 207 lines of the common plant field pennycress (Thlaspi arvense), which we collected across Europe, propagated under standardized conditions, and sequenced for their genetic and epigenetic variation. We found large geographic variation in DNA methylation, associated with both DNA sequence and climate of origin. Genetic variation was generally the stronger predictor of DNA methylation variation, but the strength of environmental association varied between different sequence contexts. Climate-of-origin was the strongest predictor in about one third of the differentially methylated regions in the CHH context, which suggests that epigenetic variation may play a role in the short-term climate adaptation of pennycress. As pennycress is currently being domesticated as a new biofuel and winter cover crop, our results may be relevant also for agriculture, particularly in changing environments.

Climate forecasts show that in many regions the temporal distribution of precipitation events will become less predictable. Root traits may play key roles in dealing with changes in precipitation predictability, but their functional plastic responses, including transgenerational processes, are scarcely known. We investigated root trait plasticity of Papaver rhoeas with respect to higher versus lower intra-seasonal and inter-seasonal precipitation predictability (i.e., the degree of temporal autocorrelation among precipitation events) during a four-year outdoor multi-generation experiment. We first tested how the simulated predictability regimes affected intra-generational plasticity of root traits and allocation strategies of the ancestors, and investigated the selective forces acting on them. Second, we exposed three descendant generations to the same predictability regime experienced by their mothers or to a different one. We then investigated whether high inter-generational predictability causes root trait differentiation, whether transgenerational root plasticity existed and whether it was affected by the different predictability treatments. We found that the number of secondary roots, root biomass and root allocation strategies of ancestors were affected by changes in precipitation predictability, in line with intra-generational plasticity. Lower predictability induced a root response, possibly reflecting a fast-acquisitive strategy that increases water absorbance from shallow soil layers. Ancestors’ root traits were generally under selection, and the predictability treatments did neither affect the strength nor the direction of selection. Transgenerational effects were detected in root biomass and root weight ratio (RWR). In presence of lower predictability, descendants significantly reduced RWR compared to ancestors, leading to an increase in performance. This points to a change in root allocation in order to maintain or increase the descendants’ fitness. Moreover, transgenerational plasticity existed in maximum rooting depth and root biomass, and the less predictable treatment promoted the lowest coefficient of variation among descendants’ treatments in five out of six root traits. This shows that the level of maternal predictability determines the variation in the descendants’ responses, and suggests that lower phenotypic plasticity evolves in less predictable environments. Overall, our findings show that roots are functional plastic traits that rapidly respond to differences in precipitation predictability, and that the plasticity and adaptation of root traits may crucially determine how climate change will affect plants.