Search

Wheat is an important staple crop since its proteins contribute to human and animal nutrition and are important for its end-use quality. However, wheat proteins can also cause adverse human reactions for a large number of people. We performed a genome wide association study (GWAS) on 114 proteins quantified by LC-MS-based proteomics and expressed in an environmentally stable manner in 148 wheat cultivars with a heritability > 0.6. For 54 proteins, we detected quantitative trait loci (QTL) that exceeded the Bonferroni-corrected significance threshold and explained 17.3–84.5% of the genotypic variance. Proteins in the same family often clustered at a very close chromosomal position or the potential homeolog. Major QTLs were found for four well-known glutenin and gliadin subunits, and the QTL segregation pattern in the protein encoding the high molecular weight glutenin subunit Dx5 could be confirmed by SDS gel-electrophoresis. For nine potential allergenic proteins, large QTLs could be identified, and their measured allele frequencies open the possibility to select for low protein abundance by markers as long as their relevance for human health has been conclusively demonstrated. A potential allergen was introduced in the beginning of 1980s that may be linked to the cluster of resistance genes introgressed on chromosome 2AS from Triticum ventricosum. The reported sequence information for the 54 major QTLs can be used to design efficient markers for future wheat breeding.

Contemporary populations are unlikely to respond to natural selection if much of their genetic variance is non-additive. Understanding the evolutionary and genomic factors that drive amounts of non-additive variance in natural populations is therefore of paramount importance. Here, we use a quantitative genetic breeding design to delineate the additive and non-additive components of expression variance of 17,656 gene transcripts segregating in the outcrossing plantArabidopsis lyrata. We find that most of the genetic variance in gene expression represents non-additive variance, especially among long genes or genes involved in epigenetic gene regulation. Genes with the most non-additive variation not only display a markedly lower rate of synonymous variation, they have also been exposed to stronger purifying selection than genes with high additive variance. This study demonstrates that the potential of standing variation for future evolution is shaped the history of purifying selection.

Polyploidy, or whole-genome duplication, is a common speciation mechanism in plants. An important barrier to polyploid establishment is a lack of compatible mates. Because self-compatibility alleviates this problem, it has long been hypothesized that there should be an association between polyploidy and self-compatibility (SC), but empirical support for this prediction is mixed. Here, we investigate whether the molecular makeup of the Brassicaceae self-incompatibility (SI) system, and specifically dominance relationships among S-haplotypes mediated by small RNAs, could facilitate loss of SI in allopolyploid crucifers. We focus on the allotetraploid species Capsella bursa-pastoris, which formed ~300 kya by hybridization and whole-genome duplication involving progenitors from the lineages of Capsella orientalis and Capsella grandiflora. We conduct targeted long-read sequencing to assemble and analyze eight full-length S-locus haplotypes, representing both homeologous subgenomes of C. bursa-pastoris. We further analyze small RNA (sRNA) sequencing data from flower buds to identify candidate dominance modifiers. We find that C. orientalis-derived S-haplotypes of C. bursa-pastoris harbor truncated versions of the male SI specificity gene SCR and express a conserved sRNA-based candidate dominance modifier with a target in the C. grandiflora-derived S-haplotype. These results suggest that pollen-level dominance may have facilitated loss of SI in C. bursa-pastoris. Finally, we demonstrate that spontaneous somatic tetraploidization after a wide cross between C. orientalis and C. grandiflora can result in production of self-compatible tetraploid offspring. We discuss the implications of this finding on the mode of formation of this widespread weed.

The rate at which plants grow is a major functional trait in plant ecology. However, little is known about its evolution in natural populations. Here, we investigate evolutionary and environmental factors shaping variation in the growth rate of Arabidopsis thaliana. We used plant diameter as a proxy to monitor plant growth over time in environments that mimicked latitudinal differences in the intensity of natural light radiation, across a set of 278 genotypes sampled within four broad regions, including an outgroup set of genotypes from China. A field experiment conducted under natural conditions confirmed the ecological relevance of the observed variation. All genotypes markedly expanded their rosette diameter when the light supply was decreased, demonstrating that environmental plasticity is a predominant source of variation to adapt plant size to prevailing light conditions. Yet, we detected significant levels of genetic variation both in growth rate and growth plasticity. Genome-wide association studies revealed that only 2 single nucleotide polymorphisms associate with genetic variation for growth above Bonferroni confidence levels. However, marginally associated variants were significantly enriched among genes with an annotated role in growth and stress reactions. Polygenic scores computed from marginally associated variants confirmed the polygenic basis of growth variation. For both light regimes, phenotypic divergence between the most distantly related population (China) and the various regions in Europe is smaller than the variation observed within Europe, indicating that the evolution of growth rate is likely to be constrained by stabilizing selection. We observed that Spanish genotypes, however, reach a significantly larger size than Northern European genotypes. Tests of adaptive divergence and analysis of the individual burden of deleterious mutations reveal that adaptive processes have played a more important role in shaping regional differences in rosette growth than maladaptive evolution., This research was funded by the European Research Council (ERC) through the “AdaptoSCOPE” grant 648617 to JdM, by the German American Fulbright Commission to JdM, and by grant AGL2016-78709-R (MEC, Spain) to S.E.R.-O. S.E.R.-O also acknowledges the financial support of the Spanish Ministry of Economy and Competitivity for the Center of Excellence Severo Ochoa 2016–2019 (SEV-2015-0533) and by the CERCA Programme/Generalitat de Catalunya. We further acknowledge grants R35 GM127131 and RO1 MH101244 from the National Institutes of Health (NIH) for S.S and E.M.K.

Understanding the causes of cis-regulatory variation is a long-standing aim in evolutionary biology. Although cis-regulatory variation has long been considered important for adaptation, we still have a limited understanding of the selective importance and genomic determinants of standing cis-regulatory variation. To address these questions, we studied the prevalence, genomic determinants, and selective forces shaping cis-regulatory variation in the outcrossing plant Capsella grandiflora We first identified a set of 1,010 genes with common cis-regulatory variation using analyses of allele-specific expression (ASE). Population genomic analyses of whole-genome sequences from 32 individuals showed that genes with common cis-regulatory variation (i) are under weaker purifying selection and (ii) undergo less frequent positive selection than other genes. We further identified genomic determinants of cis-regulatory variation. Gene body methylation (gbM) was a major factor constraining cis-regulatory variation, whereas presence of nearby transposable elements (TEs) and tissue specificity of expression increased the odds of ASE. Our results suggest that most common cis-regulatory variation in C. grandiflora is under weak purifying selection, and that gene-specific functional constraints are more important for the maintenance of cis-regulatory variation than genome-scale variation in the intensity of selection. Our results agree with previous findings that suggest TE silencing affects nearby gene expression, and provide evidence for a link between gbM and cis-regulatory constraint, possibly reflecting greater dosage sensitivity of body-methylated genes. Given the extensive conservation of gbM in flowering plants, this suggests that gbM could be an important predictor of cis-regulatory variation in a wide range of plant species.

Genomic imprinting is an epigenetic phenomenon occurring in mammals and flowering plants that causes genes to adopt a parent-of-origin-specific mode of expression. While the imprinting status of genes is well conserved in mammals, clear estimates for the degree of conservation were lacking in plants. We therefore analyzed the genome-wide imprinting status of Capsella rubella, which shared a common recent ancestor with Arabidopsis thaliana ∼10 to 14 million years ago. However, only ∼14% of maternally expressed genes (MEGs) and ∼29% of paternally expressed genes (PEGs) in C. rubella were commonly imprinted in both species, revealing that genomic imprinting is a rapidly evolving phenomenon in plants. Nevertheless, conserved PEGs exhibited signs of selection, suggesting that a subset of imprinted genes play an important functional role and are therefore maintained in plants. Like in Arabidopsis, PEGs in C. rubella are frequently associated with the presence of transposable elements that preferentially belong to helitron and MuDR families. Our data further reveal that MEGs and PEGs differ in their targeting by 24-nucleotide small RNAs and asymmetric DNA methylation, suggesting different mechanisms establishing DNA methylation at MEGs and PEGs.

The formation of an allopolyploid species involves the merger of genomes with separate evolutionary histories and thereby different genomic legacies. Contrary to expectations from theory, genes from one are often lost preferentially in allopolyploids - there is biased fractionation. Here, we provide an overview of two ways in which the genomic legacies of the progenitors may impact the fate of duplicated genes in allopolyploids. Specifically, we discuss the role of homeolog expression biases in setting the stage for biased fractionation, and the evidence for transposable element silencing as a possible mechanism for homeolog expression biases. Finally, we highlight how differences between the progenitors with respect to accumulation of deleterious variation may affect trajectories of duplicate gene evolution in allopolyploids.

International audience; Hybrid seed lethality is a widespread type of reproductive barrier among angiosperm taxa1,2 that contributes to species divergence by preventing gene flow between natural populations3,4. Besides its ecological importance, it is an important obstacle to plant breeding strategies. Hybrid seed lethality is mostly due to a failure of the nourishing endosperm tissue, resulting in embryo arrest. The cause of this failure is a parental dosage imbalance in the endosperm that can be a consequence of either differences in parental ploidy levels or differences in the 'effective ploidy', also known as the endosperm balance number (EBN)8,9. Hybrid seed defects exhibit a parent-of-origin pattern, suggesting that differences in number or expression strength of parent-of-origin-specific imprinted genes underpin, as the primary or the secondary cause, the molecular basis of the EBN7,10. Here, we have tested this concept in the genus Capsella and show that the effective ploidy of three Capsella species correlates with the number and expression level of paternally expressed genes (PEGs). Importantly, the number of PEGs and the effective ploidy decrease with the selfing history of a species: the obligate outbreeder Capsella grandiflora had the highest effective ploidy, followed by the recent selfer Capsella rubella and the ancient selfer Capsella orientalis. PEGs were associated with the presence of transposable elements and their silencing mark, DNA methylation in CHH context (where H denotes any base except C). This suggests that transposable elements have driven the imprintome divergence between Capsella species. Together, we propose that variation in transposable element insertions, the resulting differences in PEG number and divergence in their expression level form one component of the effective ploidy variation between species of different breeding system histories, and, as a consequence, allow the establishment of endosperm-based hybridization barriers.

Rapid advances in short-read DNA sequencing technologies have revolutionized population genomic studies, but there are genomic regions where this technology reaches its limits. Limitations mostly arise due to the difficulties in assembly or alignment to genomic regions of high sequence divergence and high repeat content, which are typical characteristics for loci under strong long-term balancing selection. Studying genetic diversity at such loci therefore remains challenging. Here, we investigate the feasibility and error rates associated with targeted long-read sequencing of a locus under balancing selection. For this purpose, we generated bacterial artificial chromosomes (BACs) containing the Brassicaceae S-locus, a region under strong negative frequency-dependent selection which has previously proven difficult to assemble in its entirety using short reads. We sequence S-locus BACs with single-molecule long-read sequencing technology and conduct de novo assembly of these S-locus haplotypes. By comparing repeated assemblies resulting from independent long-read sequencing runs on the same BAC clone we do not detect any structural errors, suggesting that reliable assemblies are generated, but we estimate an indel error rate of 5.7×10−5. A similar error rate was estimated based on comparison of Illumina short-read sequences and BAC assemblies. Our results show that, until de novo assembly of multiple individuals using long-read sequencing becomes feasible, targeted long-read sequencing of loci under balancing selection is a viable option with low error rates for single nucleotide polymorphisms or structural variation. We further find that short-read sequencing is a valuable complement, allowing correction of the relatively high rate of indel errors that result from this approach.

Intermediate outcrossing rates are theoretically predicted to maintain effective selection against harmful alleles, but few studies have empirically tested this prediction with the use of genomic data. We used whole-genome resequencing data from alpine rock-cress to study how genetic variation and purifying selection vary with mating system. We find that populations with intermediate outcrossing rates have similar levels of genetic diversity as outcrossing populations, and that purifying selection against harmful alleles is efficient in mixed-mating populations. In contrast, self-fertilizing populations from Scandinavia have strongly reduced genetic diversity and accumulate harmful mutations, likely as a result of demographic effects of postglacial colonization. Our results suggest that mixed-mating populations can avoid some of the negative evolutionary consequences of high self-fertilization rates.

Shifts from outcrossing to self-fertilisation have occurred repeatedly in many different lineages of flowering plants, and often involve the breakdown of genetic outcrossing mechanisms. In the Brassicaceae, self-incompatibility (SI) allows plants to ensure outcrossing by recognition and rejection of self-pollen on the stigma. This occurs through the interaction of female and male specificity components, consisting of a pistil based receptor and a pollen-coat protein, both of which are encoded by tightly linked genes at theS-locus. When benefits of selfing are higher than costs of inbreeding, theory predicts that loss-of-function mutations in the male (pollen) SI component should be favoured, especially if they are dominant. However, it remains unclear whether mutations in the male component of SI are predominantly responsible for shifts to self-compatibility, and testing this prediction has been difficult due to the challenges of sequencing the highly polymorphic and repetitive ~100 kbpS-locus. The crucifer genusCapsellaoffers an excellent opportunity to study multiple transitions from outcrossing to self-fertilization, but so far, little is known about the genetic basis and timing of loss of SI in the self-fertilizing diploidCapsella orientalis. Here, we show that loss of SI inC. orientalisoccurred within the past 2.6 Mya and maps as a dominant trait to theS-locus. Using targeted long-read sequencing of multiple complete S-haplotypes, we identify a frameshift deletion in the male specificity geneSCRthat is fixed inC. orientalis, and we confirm loss of male SI specificity. We further analyze RNA sequencing data to identify a conserved,S-linked small RNA (sRNA) that is predicted to cause dominance of self-compatibility. Our results suggest that degeneration of pollen SI specificity in dominantS-alleles is important for shifts to self-fertilization in the Brassicaceae.Author SummaryAlready Darwin was fascinated by the widely varying modes of plant reproduction. The shift from outcrossing to self-fertilization is considered one of the most frequent evolutionary transitions in flowering plants, yet we still know little about the genetic basis of these shifts. In the Brassicaceae, outcrossing is enforced by a self-incompatibility (SI) system that enables the recognition and rejection of self pollen. This occurs through the action of two tightly linked genes at theS-locus, that encode a receptor protein located on the stigma (female component) and a pollen ligand protein (male component), respectively. Nevertheless, SI has frequently been lost, and theory predicts that mutations in the male component should have an advantage during the loss of SI, especially if they are dominant. To test this hypothesis, we mapped the loss of SI in a selfing species from the genusCapsella, a model system for evolutionary genomics. We found that loss of SI mapped to theS-locus, which harbored a dominant loss-of-function mutation in the male SI protein, and as expected, we found that male specificity was indeed lost inC. orientalis. Our results suggest that transitions to selfing often involve parallel genetic changes.

Plant mating systems have profound effects on levels and structuring of genetic variation, and can affect the impact of natural selection. While theory predicts that intermediate outcrossing rates may allow plants to prevent accumulation of deleterious alleles, few studies have empirically tested this prediction using genomic data. Here, we study the effect of mating system on purifying selection by conducting population genomic analyses on whole-genome resequencing data from 38 European individuals of the arctic-alpine cruciferArabis alpina. We find that outcrossing and mixed-mating populations maintain genetic diversity at similar levels, whereas highly self-fertilizing ScandinavianA. alpinashow a strong reduction in genetic diversity, most likely as a result of a postglacial colonization bottleneck. We further find evidence for accumulation of genetic load in highly self-fertilizing populations, whereas the genome-wide impact of purifying selection does not differ greatly between mixed-mating and outcrossing populations. Our results demonstrate that intermediate levels of outcrossing may allow efficient selection against harmful alleles whereas demographic effects can be important for relaxed purifying selection in highly selfing populations. Thus, both mating system and demography shape the impact of purifying selection on genomic variation inA. alpina. These results are important for an improved understanding of the evolutionary consequences of mating system variation and the maintenance of mixed-mating strategies.SignificanceIntermediate outcrossing rates are theoretically predicted to maintain effective selection against harmful alleles, but few studies have empirically tested this prediction using genomic data. We used whole-genome resequencing data from alpine rock-cress to study how genetic variation and purifying selection vary with mating system. We find that populations with intermediate outcrossing rates have similar levels of genetic diversity as outcrossing populations, and that purifying selection against harmful alleles is efficient in mixed-mating populations. In contrast, self-fertilizing populations from Scandinavia have strongly reduced genetic diversity, and accumulate harmful mutations, likely as a result of demographic effects of postglacial colonization. Our results suggest that mixed-mating populations can avoid the negative evolutionary consequences of high self-fertilization rates.

The formation of an allopolyploid species involves the merger of two genomes with separate evolutionary histories. In allopolyploids, genes derived from one progenitor species are often expressed at higher levels than those from the other progenitor. It has been suggested that this could be due to differences in transposable element (TE) content among progenitors, as silencing of TEs can affect expression of nearby genes. Here, we examine the role of TEs for expression biases in the widespread allotetraploid Capsella bursa-pastoris and in diploid F1 hybrids generated by crossing Capsella orientalis and Capsella rubella, two close relatives of the progenitors of C. bursa-pastoris. As C. rubella harbors more TEs than C. orientalis, we expect C. orientalis alleles to be expressed at higher levels if TE content is key for expression biases. To test this hypothesis, we quantified expression biases at approximately 5800 genes in flower buds and leaves, while correcting for read mapping biases using genomic data. While three of four C. bursa-pastoris accessions exhibited a shift toward higher relative expression of C. orientalis alleles, the fourth C. bursa-pastoris accession had the opposite direction of expression bias, as did diploid F1 hybrids. Associations between TE polymorphism and expression bias were weak, and the effect of TEs on expression bias was small. These results suggest that differences in TE content alone cannot fully explain expression biases in these species. Future studies should investigate the role of differences in TE silencing efficacy, as well as a broader set of other factors. Our results are important for a more general understanding of the role of TEs for cis-regulatory evolution in plants.

Understanding the causes of cis-regulatory variation is a long-standing aim in evolutionary biology. Although cis-regulatory variation has long been considered important for adaptation, we still have a limited understanding of the selective importance and genomic determinants of standing cis-regulatory variation. To address these questions, we studied the prevalence, genomic determinants and selective forces shaping cis-regulatory variation in the outcrossing plant Capsella grandiflora. We first identified a set of 1,010 genes with common cis-regulatory variation using analyses of allele-specific expression (ASE). Population genomic analyses of whole-genome sequences from 32 individuals showed that genes with common cis-regulatory variation are 1) under weaker purifying selection and 2) undergo less frequent positive selection than other genes. We further identified genomic determinants of cis-regulatory variation. Gene-body methylation (gbM) was a major factor constraining cis-regulatory variation, whereas presence of nearby TEs and tissue specificity of expression increased the odds of ASE. Our results suggest that most common cis-regulatory variation in C. grandiflora is under weak purifying selection, and that gene-specific functional constraints are more important for the maintenance of cis-regulatory variation than genome-scale variation in the intensity of selection. Our results agree with previous findings that suggest TE silencing affects nearby gene expression, and provide novel evidence for a link between gbM and cis-regulatory constraint, possibly reflecting greater dosage-sensitivity of body-methylated genes. Given the extensive conservation of gene-body methylation in flowering plants, this suggests that gene-body methylation could be an important predictor of cis-regulatory variation in a wide range of plant species.

Whole-genome duplication (WGD) events have occurred repeatedly during flowering plant evolution, and there is growing evidence for predictable patterns of gene retention and loss following polyploidization. Despite these important insights, the rate and processes governing the earliest stages of diploidization remain poorly understood, and the relative importance of genetic drift, positive selection, and relaxed purifying selection in the process of gene degeneration and loss is unclear. Here, we conduct whole-genome resequencing in Capsella bursa-pastoris, a recently formed tetraploid with one of the most widespread species distributions of any angiosperm. Whole-genome data provide strong support for recent hybrid origins of the tetraploid species within the past 100,000–300,000 y from two diploid progenitors in the Capsella genus. Major-effect inactivating mutations are frequent, but many were inherited from the parental species and show no evidence of being fixed by positive selection. Despite a lack of large-scale gene loss, we observe a decrease in the efficacy of natural selection genome-wide due to the combined effects of demography, selfing, and genome redundancy from WGD. Our results suggest that the earliest stages of diploidization are associated with quantitative genome-wide decreases in the strength and efficacy of selection rather than rapid gene loss, and that nonfunctionalization can receive a “head start” through a legacy of deleterious variants and differential expression originating in parental diploid populations.

The selfing syndrome constitutes a suite of floral and reproductive trait changes that have evolved repeatedly across many evolutionary lineages in response to the shift to selfing. Convergent evolution of the selfing syndrome suggests that these changes are adaptive, yet our understanding of the detailed molecular genetic basis of the selfing syndrome remains limited. Here, we investigate the role of cis-regulatory changes during the recent evolution of the selfing syndrome in Capsella rubella, which split from the outcrosser Capsella grandiflora less than 200 kya. We assess allele-specific expression (ASE) in leaves and flower buds at a total of 18,452 genes in three interspecific F1 C. grandiflora x C. rubella hybrids. Using a hierarchical Bayesian approach that accounts for technical variation using genomic reads, we find evidence for extensive cis-regulatory changes. On average, 44% of the assayed genes show evidence of ASE, however only 6% show strong allelic expression biases. Flower buds, but not leaves, show an enrichment of cis-regulatory changes in genomic regions responsible for floral and reproductive trait divergence between C. rubella and C. grandiflora. We further detected an excess of heterozygous transposable element (TE) insertions near genes with ASE, and TE insertions targeted by uniquely mapping 24-nt small RNAs were associated with reduced expression of nearby genes. Our results suggest that cis-regulatory changes have been important during the recent adaptive floral evolution in Capsella and that differences in TE dynamics between selfing and outcrossing species could be important for rapid regulatory divergence in association with mating system shifts.

Whole genome duplication events have occurred repeatedly during flowering plant evolution, and there is growing evidence for predictable patterns of gene retention and loss following polyploidization. Despite these important insights, the rate and processes governing the earliest stages of diploidization remain poorly understood, and the relative importance of genetic drift, positive selection and relaxed purifying selection in the process of gene degeneration and loss is unclear. Here, we conduct whole genome resequencing in Capsella bursa-pastoris, a recently formed tetraploid with one of the most widespread species distributions of any angiosperm. Whole genome data provide strong support for recent hybrid origins of the tetraploid species within the last 100-300,000 years from two diploid progenitors in the Capsella genus. Major-effect inactivating mutations are frequent, but many were inherited from the parental species and show no evidence of being fixed by positive selection. Despite a lack of large-scale gene loss, we observe a decrease in the efficacy of natural selection genome-wide, due to the combined effects of demography, selfing and genome redundancy from whole genome duplication. Our results suggest that the earliest stages of diploidization are associated with quantitative genome-wide decreases in the strength and efficacy of selection rather than rapid gene loss, and that non-functionalization can receive a 'head start' through a legacy of deleterious variants and differential expression originating in parental diploid populations.

International audience; The shift from outcrossing to selfing is common in flowering plants, but the genomic consequences and the speed at which they emerge remain poorly understood. An excellent model for understanding the evolution of self fertilization is provided by Capsella rubella, which became self compatible

The wild tomato species Solanum chilense and S. peruvianum are a valuable non-model system for studying plant adaptation since they grow in diverse environments facing many abiotic constraints. Here we investigate the sequence evolution of regulatory regions of drought and cold responsive genes and their expression regulation. The coding regions of these genes were previously shown to exhibit signatures of positive selection. Expression profiles and sequence evolution of regulatory regions of members of the Asr (ABA/water stress/ripening induced) gene family and the dehydrin gene pLC30-15 were analyzed in wild tomato populations from contrasting environments. For S. chilense, we found that Asr4 and pLC30-15 appear to respond much faster to drought conditions in accessions from very dry environments than accessions from more mesic locations. Sequence analysis suggests that the promoter of Asr2 and the downstream region of pLC30-15 are under positive selection in some local populations of S. chilense. By investigating gene expression differences at the population level we provide further support of our previous conclusions that Asr2, Asr4, and pLC30-15 are promising candidates for functional studies of adaptation. Our analysis also demonstrates the power of the candidate gene approach in evolutionary biology research and highlights the importance of wild Solanum species as a genetic resource for their cultivated relatives.

Phenotypic plasticity is the variation in phenotype that a single genotype can produce in different environments and, as such, is an important component of individual fitness. However, whether the effect of new mutations, and hence evolution, depends on the direction of plasticity remains controversial. Here, we identify the cis-acting modifications that have reshaped gene expression in response to dehydration stress in three Arabidopsis species. Our study shows that the direction of effects of most cis-regulatory variants differentiating the response between A. thaliana and the sister species A. lyrata and A. halleri depends on the direction of pre-existing plasticity in gene expression. A comparison of the rate of cis-acting variant accumulation in each lineage indicates that the selective forces driving adaptive evolution in gene expression favors regulatory changes that magnify the stress response in A. lyrata. The evolutionary constraints measured on the amino-acid sequence of these genes support this interpretation. In contrast, regulatory changes that mitigate the plastic response to stress evolved more frequently in A. halleri. Our results demonstrate that pre-existing plasticity may be a stepping stone for adaptation, but its selective remodeling differs between lineages., Plasticity allows organisms to respond to environmental change. Here the authors compare the distribution of cis-regulatory variants in the transcriptomes of Arabidopsis lyrata and A. halleriafter exposure to stress, to trace the role of polygenic selection in the evolution of gene expression plasticity.