Search

The challenges to which plants are exposed in urban environments represent, in miniature, the challenges plants face as a result of global environmental change. Hence, urban habitats provide a unique opportunity to assess whether processes of local adaptation are taking place despite the short temporal and geographical scales that characterize the anthropoceneWe quantified the ecological diversity of spontaneously occurring urban habitat patches of A. thaliana. Using plant community indicators, we show that these patches differ in their levels of soil nutrient content and disturbance. Accordingly, plants in each patch displayed significantly different flowering time, size, and fitness.Using a deep sampling approach coupled with reduced genome-sequencing, we demonstrate that most individuals can be assigned to a limited set of clonal lineages; the genetic diversity of these lineages represents the diversity observed in western European populations of the species, indicating that established urban populations originate from a broad regional pool of lineages.We assessed the genetic and phenotypic diversity of these lineages in a set of common garden experiments. We report marked genetic differences in life-history traits, including time of primary and secondary dormancy as well as of flowering. These genetic differences in life-history traits are not randomly distributed but sorted out by ecological differences among sites of origin.Synthesis: Our study shows that the genetically diverse phenology of a regional A. thaliana gene pool is not randomly distributed but filtered by urban environmental heterogeneity. This report is the first to show a pattern of local genetic adaptation within urban environments. We conclude that environmental filtering helps maintain functional diversity within species.

Rapid evolution of traits and of plasticity may enable adaptation to climate change, yet solid experimental evidence under natural conditions is scarce. Here, we imposed rainfall manipulations (+30%, control, -30%) for ten years on entire natural plant communities in two Eastern Mediterranean sites. Additional sites along a natural rainfall gradient and selection analyses in a greenhouse assessed whether potential responses were adaptive. In both sites, our annual target species Biscutella didyma consistently evolved earlier phenology and higher reproductive allocation under drought. Multiple arguments suggest that this response was adaptive: it aligned with theory, corresponding trait shifts along the natural rainfall gradient, and selection analyses under differential watering in the greenhouse. However, another seven candidate traits did not evolve, and there was little support for evolution of plasticity. Our results provide compelling evidence for rapid adaptive evolution under climate change. Yet, several non-evolving traits may indicate potential constraints to full adaptation.

Balancing selection is a form of natural selection maintaining diversity at the sites it targets and at linked nucleotide sites. Due to selection favouring heterozygosity, it has the potential to facilitate the accumulation of a “sheltered” load of tightly linked recessive deleterious mutations. However, precisely evaluating the extent of these effects has remained challenging. Taking advantage of plant self-incompatibility as one of the best-understood examples of long-term balancing selection, we provide a highly resolved picture of the genomic extent of balancing selection on the sheltered genetic load. We used targeted genome resequencing to reveal polymorphism of the genomic region flanking the self-incompatibility locus in three sample sets in each of the two closely related plant speciesArabidopsis halleriandA. lyrata, and used 100 control regions from throughout the genome to factor out differences in demographic histories and/or sample structure. Nucleotide polymorphism increased strongly around theS-locus in all sample sets, but only over a limited genomic region, as it became indistinguishable from the genomic background beyond the first 25-30kb. Genes in this chromosomal interval exhibited no excess of mutations at 0-fold degenerated sites relative to putatively neutral sites, hence revealing no detectable weakening of the efficacy of purifying selection even for these most tightly linked genes. Overall, our results are consistent with the predictions of a narrow genomic influence of linkage to theS-locus, and clarify how natural selection in one genomic region affects the evolution of the adjacent genomic regions.

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.

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.

Achieving high intraspecific genetic diversity is a critical goal in ecological restoration as it increases the adaptive potential and long-term resilience of populations. Thus, we investigated genetic diversity within and between pristine sites in a fossil floodplain and compared it to sites restored by hay-transfer between 1997 and 2014. RAD-seq genotyping revealed that the stenoecious flood-plain species Arabis nemorensis is co-occurring with individuals that, based on ploidy, ITS-sequencing and morphology, probably belong to the close relative Arabis sagittata, which has a documented preference for dry calcareous grasslands but has not been reported in floodplain meadows. We show that hay-transfer maintains genetic diversity for both species. Additionally, in A. sagittata, transfer from multiple genetically isolated pristine sites resulted in restored sites with increased diversity and admixed local genotypes. In A. nemorensis, transfer did not create novel admixture dynamics because genetic diversity between pristine sites was less differentiated. Thus, the effects of hay-transfer on genetic diversity also depend on the genetic makeup of the donor communities of each species, especially when local material is mixed. Our results demonstrate the efficiency of hay-transfer for habitat restoration and emphasize the importance of pre-restoration characterization of micro-geographic patterns of intraspecific diversity of the community to guarantee that restoration practices reach their goal, i.e. maximize the adaptive potential of the entire restored plant community. Overlooking these patterns may alter the balance between species in the community. Additionally, our comparison of summary statistics obtained from de novo and reference-based RAD-seq pipelines shows that the genomic impact of restoration can be reliably monitored in species lacking prior genomic knowledge.

Local adaptation across species range is widespread. Yet, much has to be discovered on its environmental drivers, the underlying functional traits and their molecular determinants. Because elevation gradients display continuous environmental changes at a short geographical scale, they provide an exceptional opportunity to investigate these questions. Here, we used two common gardens to phenotype 1664 plants from 11 populations of annual teosintes. These populations were sampled across two elevation gradients in Mexico. Our results point to a syndrome of adaptation to altitude with the production of offspring that flowered earlier, produced less tillers, and larger, longer and heavier grains with increasing elevation. We genotyped these plants for 178 outlier single nucleotide polymorphisms (SNPs), which had been chosen because they displayed excess of allele differentiation and/or correlation with environmental variables in six populations with contrasted altitudes. A high proportion of outlier SNPs associated with the phenotypic variation of at least one trait. We tested phenotypic pairwise correlations between traits, and found that the higher the correlation, the greater the number of common associated SNPs. In addition, allele frequencies at 87 of the outlier SNPs correlated with an environmental component best summarized by altitudinal variation on a broad sample of 28 populations. Chromosomal inversions were enriched for both phenotypically-associated and environmentally-correlated SNPs. Altogether, our results are consistent with the set-up of an altitudinal syndrome promoted by local adaptation of teosinte populations in the face of gene flow. We showed that pleiotropy is pervasive and potentially has constrained the evolution of traits. Finally, we recovered variants underlying phenotypic variation at adaptive traits. Because elevation mimics climate change through space, these variants may be relevant for future maize breeding.Author summaryAcross their native range, species encounter a diversity of habitats promoting local adaptation of geographically distributed populations. While local adaptation is widespread, much has yet to be discovered about the conditions of its emergence, the targeted traits, their molecular determinants and the underlying ecological drivers. Here we employed a reverse ecology approach, combining phenotypes and genotypes, to mine the determinants of local adaptation of teosinte populations distributed along two steep altitudinal gradients in Mexico. Evaluation of 11 populations in two common gardens located at mid-elevation pointed to the set-up of an altitudinal syndrome, in spite of gene flow. We scanned genomes to identify loci with allele frequencies shifts along elevation. Interestingly, variation at these loci was commonly associated to variation of phenotypes. Because elevation mimics climate change through space, these variants may be relevant for future maize breeding.

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

Background and Aims Water limitation is an important determinant of the distribution, abundance and diversity of plant species. Yet, little is known about how the response to limiting water supply changes among closely related plant species with distinct ecological preferences. Comparison of the model annual species A. thaliana to its close perennial relatives A. lyrata and A. halleri, can help disentangle the molecular and physiological changes contributing to tolerance and avoidance mechanisms, because these species must maintain tolerance and avoidance mechanisms to increase long-term survival, but they are exposed to different levels of water stress and competition in their natural habitat.Methods We conducted a dry-down experiment that mimics a period of missing precipitation. We quantified the covariation of progressive decrease in soil water content (SWC) with various physiological and morphological plant traits across a set of representative genotypes in Arabidopsis thaliana, A. lyrata and A. halleri. To quantify the degree of plant stress, transcriptome changes were also monitored.Key Results The analysis of trait co-variation demonstrates that the three species differ in the strategies they deploy to respond to drought stress. A. thaliana showed drought avoidance reaction but failed to survive wilting. A. lyrata efficiently combined avoidance and tolerance mechanisms. By contrast, A. halleri showed some degree of tolerance to wilting but it did not seem to protect itself from the stress imposed by drought. Transcriptome data collected just before plant wilting and after recovery corroborated the phenotypic analysis, with A. lyrata and A. halleri showing a stronger activation of recovery- and stress-related genes, respectively.Conclusions We conclude that these three Arabidopsis species have evolved distinct strategies to face drought stress, and discuss the extent to which these strategic differences reflect their respective ecological priorities.

Stomata control gas exchanges between the plant and the atmosphere. How natural variation in stomata size and density contributes to resolve trade-offs between carbon uptake and water-loss in response to local climatic variation is not yet understood. We developed an automated confocal microscopy approach to characterize natural genetic variation in stomatal patterning in 330 fully-sequenced Arabidopsis thaliana accessions collected throughout the European range of the species. We compared this to variation in water-use efficiency, measured as carbon isotope discrimination (δ13C). We detect substantial genetic variation for stomata size and density segregating within Arabidopsis thaliana. A positive correlation between stomata size and δ13C further suggests that this variation has consequences on water-use efficiency. Genome-wide association analyses indicate a complex genetic architecture underlying not only variation in stomata patterning but also to its co-variation with carbon uptake parameters. Yet, we report two novel QTL affecting δ13C independently of stomata patterning. This suggests that, in A. thaliana, both morphological and physiological variants contribute to genetic variance in water-use efficiency. Patterns of regional differentiation and co-variation with climatic parameters indicate that natural selection has contributed to shape some of this variation, especially in Southern Sweden, where water availability is more limited in spring relative to summer. These conditions are expected to favor the evolution of drought avoidance mechanisms over drought escape strategies.

Understanding the genetics of lifetime fitness is crucial to understand a species’ ecological preferences and ultimately predict its ability to cope with novel environmental conditions. Yet, there is a dearth of information regarding the impact of the ecological variance experienced by natural populations on expressed phenotypic and fitness differences. Here, we follow the natural dynamics of experimental A. thaliana populations over 5 successive plantings whose timing was determined by the natural progression of the plant’s life cycle and disentangle the environmental and genetic factors that drive plant ecological performance at a given locality. We show that, at the intermediate latitude where the experiment was conducted, a given genotype can experience different life cycles across successive seasons. Lifetime fitness across these seasons varied strongly, with a fall planting yielding 36-fold higher fitness compared to a spring planting. In addition, the actual life-stage at which plant overwinter oscillated across years, depending on the timing of the end of the summer season. We observed a rare but severe fitness differential after inadequate early flowering in one of the five planting. Substrate variation played a comparatively minor role, but also contributed to modulate the magnitude of fitness differentials between genotypes. Finally, reciprocal introgressions on chromosome 4 demonstrated that the fitness effect of a specific chromosomal region is strongly contingent on micro-geographic and seasonal fluctuations. Our study contributes to emphasize the extent to which the fitness impact of phenotypic traits and the genes that encode them in the genome can fluctuate. Experiments aiming at dissecting the molecular basis of local adaptation must apprehend the complexity introduced by temporal fluctuations because they massively affect the expression of phenotype and fitness differences.

Elucidating the molecular basis of natural variation in complex traits is the key for their effective management in crops or natural systems. This review focuses on plant variation. It will first, show that genetic modifications causing major alterations in polygenic phenotypes often hit targets within an array of 'candidate genes', second, present new methods that include mutations of all effect sizes, and help exhaustively describe the molecular systems underlying complex traits, and third, discuss recent findings regarding the role of epigenetic variants, which in plants are often maintained through both mitosis and meiosis. Exploring the whole spectrum of mutations controlling complex traits is made possible by the combination of genetic, genomic and epigenomic approaches.

Local adaptation provides an opportunity to study the genetic basis of adaptation and investigate the allelic architecture of adaptive genes. We study delay of germination 1 (DOG1), a gene controlling natural variation in seed dormancy in Arabidopsis thaliana and investigate evolution of dormancy in 41 populations distributed in four regions separated by natural barriers. Using F(ST) and Q(ST) comparisons, we compare variation at DOG1 with neutral markers and quantitative variation in seed dormancy. Patterns of genetic differentiation among populations suggest that the gene DOG1 contributes to local adaptation. Although Q(ST) for seed dormancy is not different from F(ST) for neutral markers, a correlation with variation in summer precipitation supports that seed dormancy is adaptive. We characterize dormancy variation in several F(2) -populations and show that a series of functionally distinct alleles segregate at the DOG1 locus. Theoretical models have shown that the number and effect of alleles segregatin at quantitative trait loci (QTL) have important consequences for adaptation. Our results provide support to models postulating a large number of alleles at quantitative trait loci involved in adaptation.

Much is known about the evolution of plant immunity components directed against specific pathogen strains: They show pervasive functional variation and have the potential to coevolve with pathogen populations. However, plants are effectively protected against most microbes by generalist immunity components that detect conserved pathogen-associated molecular patterns (PAMPs) and control the onset of PAMP-triggered immunity. In Arabidopsis thaliana, the receptor kinase flagellin sensing 2 (FLS2) confers recognition of bacterial flagellin (flg22) and activates a manifold defense response. To decipher the evolution of this system, we performed functional assays across a large set of A. thaliana genotypes and Brassicaceae relatives. We reveal extensive variation in flg22 perception, most of which results from changes in protein abundance. The observed variation correlates with both the severity of elicited defense responses and bacterial proliferation. We analyzed nucleotide variation segregating at FLS2 in A. thaliana and detected a pattern of variation suggestive of the rapid fixation of a novel adaptive allele. However, our study also shows that evolution at the receptor locus alone does not explain the evolution of flagellin perception; instead, components common to pathways downstream of PAMP perception likely contribute to the observed quantitative variation. Within and among close relatives, PAMP perception evolves quantitatively, which contrasts with the changes in recognition typically associated with the evolution of R genes.

We investigate how different rates of environmental change affect adaptive outcomes and dynamics by selecting Chlamydomonas populations for over 200 generations in environments where the rate of change varies. We find that slower rates of environmental change result in end populations that grow faster and pay a lower cost of adaptation than populations that adapt to a sudden change of the same magnitude. We detected partial selective sweeps in adapting populations by monitoring changes in marker frequency in each population. Although populations adapting to a sudden environmental change showed evidence of mutations of large effect segregating early on, populations adapting to slow rates of change showed patterns that were consistent with mutations of relatively small effect occurring at less predictable times. This work suggests that rates of environmental change may fundamentally alter adaptive dynamics and outcomes of adaptation by changing the size and timing of fitness increases. We suggest that using mutations of smaller effect during adaptation may result in lower levels of pleiotropy and historical constraints, which could in turn result in higher fitness by the end of the experiment.

In plants and animals, gene expression can be down-regulated at the posttranscriptional level by microRNAs (miRNAs), a class of small endogenous RNA. Comparative analysis of miRNA content across species indicates continuous birth and death of these loci in the course of evolution. However, little is known about the microevolutionary dynamics of these genetic elements, especially in plants. In this article we examine polymorphism at two miRNA-encoding loci in Arabidopsis thaliana , miR856 and miR824, which are not found in rice or poplar. We compare their diversity to other miRNA-encoding loci conserved across distant taxa. We find that levels of variation vary significantly across loci and that the two recently derived loci harbor patterns of diversity deviating from neutrality. miRNA miR856 shows a weak signature of a selective sweep whereas miR824 displays signs of balancing selection. A detailed examination of structural variation among alleles found at the miR824-encoding locus suggests nonrandom evolution of a thermoresistant substructure in the precursor. Expression analysis of pre-miR824 and its target, AGL16, indicates that these structural differences likely impact the processing of mature miR824. Our work highlights the relevance of RNA structure in precursor sequence evolution, suggesting that the evolutionary dynamics of miRNA-encoding loci is more complex than suggested by the constraints exerted on the interaction between mature miRNA fragments and their target exon.

Adaptation of a complex trait often requires the accumulation of many modifications to finely tune its underpinning molecular components to novel environmental requirements. The investigation of cis-acting regulatory modifications can be used to pinpoint molecular systems partaking in such complex adaptations. Here, we identify cis-acting modifications with the help of an interspecific crossing scheme designed to distinguish modifications derived in each of the two sister species, Arabidopsis halleri and A. lyrata. Allele-specific expression levels were assessed in three environmental conditions chosen to reflect interspecific ecological differences: cold exposure, dehydration, and standard conditions. The functions described by Gene Ontology categories enriched in cis-acting mutations are markedly different in A. halleri and A. lyrata, suggesting that polygenic adaptation reshaped distinct polygenic molecular functions in the two species. In the A. halleri lineage, an excess of cis-acting changes affecting metal transport and homeostasis was observed, confirming that the well-known heavy metal tolerance of this species is the result of polygenic selection. In A. lyrata, we find a marked excess of cis-acting changes among genes showing a transcriptional response to cold stress in the outgroup species A. thaliana. The adaptive relevance of these changes will have to be validated. We finally observed that polygenic molecular functions enriched in derived cis-acting changes are more constrained at the amino acid level. Using the distribution of cis-acting variation to tackle the polygenic basis of adaptation thus reveals the contribution of mutations of small effect to Darwinian adaptation. Fil: He, Fei. University Of Cologne; Alemania Fil: Arce, Agustín Lucas. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. University Of Cologne; Alemania Fil: Schmitz, Gregor. University Of Cologne; Alemania Fil: Koornneef, Maarten. Max Planck Institute For Plant Breeding Research; Alemania Fil: Novikova, Polina. Gregor Mendel Institute; Austria Fil: Beyer, Andreas. University Of Cologne; Alemania Fil: De Meaux, Juliette. University Of Cologne; Alemania

We investigate how the dynamics and outcomes of adaptation by natural selection are affected by environmental stability by simulating adaptive walks in response to an environmental change of fixed magnitude but variable speed. Here we consider monomorphic lineages that adapt by the sequential fixation of beneficial mutations. This is modeled by selecting short RNA sequences for folding stability and secondary structure conservation at increasing temperatures. Using short RNA sequences allows us to describe adaptive outcomes in terms of genotype (sequence) and phenotype (secondary structure) and to follow the dynamics of fitness increase. We find that slower rates of environmental change affect the dynamics of adaptive walks by reducing the fitness effect of fixed beneficial mutations, as well as by increasing the range of time in which the substitutions of largest effect are likely to occur. In addition, adaptation to slower rates of environmental change results in fitter endpoints with fewer possible end phenotypes relative to lineages that adapt to a sudden change. This suggests that care should be taken when experiments using sudden environmental changes are used to make predictions about adaptive responses to gradual change.

Arabidopsis thaliana is a model plant species and its molecular dissection has greatly contributed to our understanding of the systems preventing genome invasion by transposable elements (TE). Recent advances suggest that A. thaliana may be more efficient than its congener A. lyrata at controlling TE expression and proliferation. The comparative analysis of TE transcription in A. thaliana and A. lyrata, which differ by 40% in genome size, may help understand how silencing mechanisms contribute to the evolution of transposition rate, an important factor controlling genome size variation in plants and animals.

The timing of flowering is pivotal for maximizing reproductive success under fluctuating environmental conditions. Flowering time is tightly controlled by complex genetic networks that integrate endogenous and exogenous cues, such as light, temperature, photoperiod, and hormones. Here, we show that AGAMOUS-LIKE16 (AGL16) and its negative regulator microRNA824 (miR824) control flowering time in Arabidopsis thaliana. Knockout of AGL16 effectively accelerates flowering in nonvernalized Col-FRI, in which the floral inhibitor FLOWERING LOCUS C (FLC) is strongly expressed, but shows no effect if plants are vernalized or grown in short days. Alteration of AGL16 expression levels by manipulating miR824 abundance influences the timing of flowering quantitatively, depending on the expression level and number of functional FLC alleles. The effect of AGL16 is fully dependent on the presence of FLOWERING LOCUS T (FT). Further experiments show that AGL16 can interact directly with SHORT VEGETATIVE PHASE and indirectly with FLC, two proteins that form a complex to repress expression of FT. Our data reveal that miR824 and AGL16 modulate the extent of flowering time repression in a long-day photoperiod.

Arabidopsis seeds rapidly release hydrophilic polysaccharides from the seed coat on imbibition. These form a heavy mucilage layer around the seed that makes it sink in water. Fourteen natural Arabidopsis variants from central Asia and Scandinavia were identified with seeds that have modified mucilage release and float. Four of these have a novel mucilage phenotype with almost none of the released mucilage adhering to the seed and the absence of cellulose microfibrils. Mucilage release was modified in the variants by ten independent causal mutations in four different loci. Seven distinct mutations affected one locus, coding the MUM2 β-D-galactosidase, and represent a striking example of allelic heterogeneity. The modification of mucilage release has thus evolved a number of times independently in two restricted geographical zones. All the natural mutants identified still accumulated mucilage polysaccharides in seed coat epidermal cells. Using nuclear magnetic resonance (NMR) relaxometry their production and retention was shown to reduce water mobility into internal seed tissues during imbibition, which would help to maintain seed buoyancy. Surprisingly, despite released mucilage being an excellent hydrogel it did not increase the rate of water uptake by internal seed tissues and is more likely to play a role in retaining water around the seed., Author Summary Seeds of the model plant Arabidopsis release sticky mucilage on imbibition that is constituted of complex polysaccharides. In this study, we have identified and characterised natural Arabidopsis variants that do not release mucilage and found that their seeds float. The accumulation of unreleased polysaccharides in the seed coat reduced water uptake rates on imbibition and would maintain buoyancy. We subsequently identified additional floating natural variants where mucilage is released, but is not attached to the seed, apparently due to defective cellulose production. The different variants arise from at least ten independent unique mutations and were collected from two discrete geographical areas. Arabidopsis seed flotation has thus evolved several times due to modifications in mucilage release. Released mucilage was found to retain water, but did not improve imbibition of internal seed tissues, indicating a role in maintaining seeds hydrated. These findings highlight the physical and potential physiological effects of mucilage production by the seed coat.

RNA molecules follow a succession of enzyme-mediated processing steps from transcription to maturation. The participating enzymes, for example the spliceosome for mRNAs and Drosha and Dicer for microRNAs, are also produced in the cell and their copy-numbers fluctuate over time. Enzyme copy-number changes affect the processing rate of the substrate molecules; high enzyme numbers increase the processing rate, while low enzyme numbers decrease it. We study different RNA-processing cascades where enzyme copy-numbers are either fixed or fluctuate. We find that for the fixed enzyme copy-numbers, the substrates at steady-state are Poisson-distributed, and the whole RNA cascade dynamics can be understood as a single birth-death process of the mature RNA product. In this case, solely fluctuations in the timing of RNA processing lead to variation in the number of RNA molecules. However, we show analytically and numerically that when enzyme copy-numbers fluctuate, the strength of RNA fluctuations increases linearly with the RNA transcription rate. This linear effect becomes stronger as the speed of enzyme dynamics decreases relative to the speed of RNA dynamics. Interestingly, we find that under certain conditions, the RNA cascade can reduce the strength of fluctuations in the expression level of the mature RNA product. Finally, by investigating the effects of processing polymorphisms, we show that it is possible for the effects of transcriptional polymorphisms to be enhanced, reduced or even reversed. Our results provide a framework to understand the dynamics of RNA processing.

Local adaptation provides an opportunity to study the genetic basis of adaptation and investigate the allelic architecture of adaptive genes. We study delay of germination 1 (DOG1), a gene controlling natural variation in seed dormancy in Arabidopsis thaliana and investigate evolution of dormancy in 41 populations distributed in four regions separated by natural barriers. Using F(ST) and Q(ST) comparisons, we compare variation at DOG1 with neutral markers and quantitative variation in seed dormancy. Patterns of genetic differentiation among populations suggest that the gene DOG1 contributes to local adaptation. Although Q(ST) for seed dormancy is not different from F(ST) for neutral markers, a correlation with variation in summer precipitation supports that seed dormancy is adaptive. We characterize dormancy variation in several F(2) -populations and show that a series of functionally distinct alleles segregate at the DOG1 locus. Theoretical models have shown that the number and effect of alleles segregatin at quantitative trait loci (QTL) have important consequences for adaptation. Our results provide support to models postulating a large number of alleles at quantitative trait loci involved in adaptation.

Transposable elements (TEs) are so abundant and variable that they count among the most important mutational sources in genomes. Nonetheless, little is known about the genetics of their variation in activity or silencing across closely related species. Here, we demonstrate that regulation of TE genes can differ dramatically between the two closely related Arabidopsis species A. thaliana and A. lyrata. In leaf and floral tissues of F1 interspecific hybrids, about 47% of TEs show allele-specific expression, with the A. lyrata copy being generally expressed at higher level. We confirm that TEs are generally expressed in A. lyrata but not in A. thaliana. Allele-specific differences in TE expression are associated with divergence in epigenetic modifications like DNA and histone methylation between species as well as with sequence divergence. Our data demonstrate that A. thaliana silences TEs much better than A. lyrata. For long terminal repeat retrotransposons, these differences are more pronounced for younger insertions. Interspecific differences in TE silencing may have a great impact on genome size changes.

Seasonal germination timing of Arabidopsis thaliana strongly influences overall life history expression and is the target of intense natural selection. This seasonal germination timing depends strongly on the interaction between genetics and seasonal environments both before and after seed dispersal. DELAY OF GERMINATION 1 (DOG1) is the first gene that has been identified to be associated with natural variation in primary dormancy in A. thaliana. Here, we report interaccession variation in DOG1 expression and document that DOG1 expression is associated with seed-maturation temperature effects on germination; DOG1 expression increased when seeds were matured at low temperature, and this increased expression was associated with increased dormancy of those seeds. Variation in DOG1 expression suggests a geographical structure such that southern accessions, which are more dormant, tend to initiate DOG1 expression earlier during seed maturation and achieved higher expression levels at the end of silique development than did northern accessions. Although elimination of the synthesis of phytohormone abscisic acid (ABA) results in the elimination of maternal temperature effects on dormancy, DOG1 expression predicted dormancy better than expression of genes involved in ABA metabolism.

There is now ample evidence that plant development, responses to abiotic environments, and immune responses are tightly intertwined in their physiology. Thus optimization of the immune system during evolution will occur in coordination with that of plant development. Two alternative and possibly complementary forces are at play: genetic constraints due to the pleiotropic action of players in both systems, and coevolution, if developmental changes modulate the cost-benefit balance of immunity. A current challenge is to elucidate the ecological forces driving evolution of quantitative variation for defense at molecular level. The analysis of natural co-variation for developmental and immunity traits in Arabidopsis thaliana promises to bring important insights.

SummaryMorphological diversity is often caused by altered gene expression of key developmental regulators. However, the precise developmental trajectories through which morphologies evolved remain poorly understood. It is also unclear to what degree genetic changes contributing to morphological divergence were fixed by natural selection. Here we investigate these problems in the context of evolutionary developmental transitions that produced the simple unlobed leaf of the model species Arabidopsis thaliana. We demonstrate that A. thaliana leaf shape likely derived from a more complex lobed ancestral state that persists in extant Arabidopsis species. We also show that evolution of the unlobed leaf form in A. thaliana involved loss of expression of the knotted1-like homeobox gene SHOOTMERISTEMLESS (STM) in leaves and that cis-regulatory divergence contributed to this process. Further, we provide evidence for a selective sweep at the A. thaliana STM locus, indicating that loss of STM expression in A. thaliana leaves may have been fixed by positive selection. In summary, our data provide key information as to when and how the characteristic leaf form of A. thaliana evolved.

The orientation of flanking genes may influence the evolution of intergenic regions in which cis-regulatory elements are likely to be located: divergently transcribed genes share their 5′ regions, resulting either in smaller “private” spaces or in overlapping regulatory elements. Thus, upstream sequences of divergently transcribed genes (bi-directional upstream regions, or URs) may be more constrained than those of uni-directional gene pairs. We investigated this effect by analyzing nucleotide variation segregating within and between Arabidopsis species. Compared to uni-directional URs, bi-directional URs indeed display lower population mutation rate, as well as more low-frequency polymorphisms. Furthermore, we find that bi-directional regions undergo selection for the maintenance of intergenic distance. Altogether, however, we observe considerable variation in evolutionary rates, with putative signatures of selection on two uni-directional upstream regions., This research was funded by the German Funding Agency (SFB680).

International audience; Although pioneered by human geneticists as a potential solution to the challenging problem of finding the genetic basis of common human diseases1,2, genome-wide association (GWA) studies have, owing to advances in genotyping and sequencing technology, become an obvious general approach for studying the genetics of natural variation and traits of agricultural importance. They are particularly useful when inbred lines are available, because once these lines have been genotyped they can be phenotyped multiple times, making it possible (as well as extremely cost effective) to study many different traits in many different environments, while replicating the phenotypic measurements to reduce environmental noise. Here we demonstrate the power of this approach by carrying out a GWA study of 107 phenotypes in Arabidopsis thaliana, a widely distributed, predominantly self-fertilizing model plant known to harbour considerable genetic variation for many adaptively important traits3. Our results are dramatically different from those of human GWA studies, in that we identify many common alleles of major effect, but they are also, in many cases, harder to interpret because confounding by complex genetics and population structure make it difficult to distinguish true associations from false. However, a-priori candidates are significantly over-represented among these associations as well, making many of them excellent candidates for follow-up experiments. Our study demonstrates the feasibility of GWA studies in A. thaliana and suggests that the approach will be appropriate for many other organisms.

Background In all eukaryotic species examined, genes that are chromosomal neighbours are more similar in their expression than random gene pairs. Currently, it is still unclear how much of this local co-expression is caused by direct transcriptional interactions, and how much is due to shared chromatin environments. Results We analysed neighbouring genes in Arabidopsis thaliana. At large intergenic distances (>400 bp), divergently and convergently transcribed gene pairs show very similar levels of co-expression, mediated most likely by shared chromatin environments. At gene distances below 400 bp, co-expression is strongly enhanced only for divergently transcribed gene pairs, indicating bi-directional transcription from a single promoter. Conversely, co-expression is suppressed for short convergently or uni-directionally transcribed pairs. This suppression points to transcriptional interference concentrated at the 3' end, e.g., in the context of transcription termination. Conclusions Classifying linked gene pairs by their orientation, we are able to partially tease apart the different levels of regional expression modulation. (i) Regional chromatin characteristics modulate the accessibility for regulation and transcription, regardless of gene orientation; the strength of this chromatin effect can be assessed from divergently or convergently transcribed distant neighbours. (ii) Shared promoter regions up to 400 bp in length enhance the co-expression of close bi-directional neighbours. (iii) Transcriptional interference of close neighbours is concentrated at the 3' ends of genes, and reduces co-expression on average by 40%.

The contribution of cis-regulation to adaptive evolutionary change is believed to be essential, yet little is known about the evolutionary rules that govern regulatory sequences. Here, we characterize the short-term evolutionary dynamics of a cis-regulatory region within and among two closely related species, A. lyrata and A. halleri, and compare our findings to A. thaliana. We focused on the cis-regulatory region of chalcone synthase (CHS), a key enzyme involved in the synthesis of plant secondary metabolites. We observed patterns of nucleotide diversity that differ among species but do not depart from neutral expectations. Using intra- and interspecific F1 progeny, we have evaluated functional cis-regulatory variation in response to light and herbivory, environmental cues, which are known to induce CHS expression. We find that substantial cis-regulatory variation segregates within and among populations as well as between species, some of which results from interspecific genetic introgression. We further demonstrate that, in A. thaliana, CHS cis-regulation in response to herbivory is greater than in A. lyrata or A. halleri. Our work indicates that the evolutionary dynamics of a cis-regulatory region is characterized by pervasive functional variation, achieved mostly by modification of response modules to one but not all environmental cues. Our study did not detect the footprint of selection on this variation.

Promoters are thought to play a major role in adaptive evolution, yet little is known about the regulatory diversity within species, where microevolutionary processes take place. To investigate the potential for evolutionary change in the promoter of a gene, we examined nucleotide and functional variation of the Chalcone Synthase (CHS) cis-regulatory region in Arabidopsis thaliana. CHS is the branch point enzyme of a biosynthetic pathway that leads to the production of secondary metabolites influencing the interaction between the plant and its environment. We found that nucleotide diversity in the intergenic region encompassing the CHS promoter (π = 0.003) is compatible with neutral expectations. To quantify functional variation specifically as a result of cis-regulation of CHS mRNA levels, we developed an assay using F1 individuals in which distinct promoter alleles are compared within a common trans-regulatory background. We examined functional cis-regulatory variation in response to different stimuli representing a variety of CHS transcriptional environments (dark, light, and insect feeding). We observed extensive functional variation, some of which appeared to be independent of the trans-regulatory background. Comparison of functional and nucleotide diversity suggested a candidate point mutation that may explain cis-regulatory differences in light response. Our results indicate that functional changes in promoters can arise from a few mutations, pointing to promoter regions as a fundamental determinant of functional genetic variation.

A good understanding of the dynamics of host/pathogen interactions and of the factors that shape the spatial distribution of resistance genes is a prerequisite of metapopulation dynamic management of resistance genes. We studied the diversity and spatial structure of natural populations of common bean (Phaseolus vulgaris) for resistance to Colletotrichum lindemuthianum, the causal agent of anthracnose. This study was carried out in Mexico and Argentina both at the phenotypic level and at the molecular level for two families of resistance gene candidates (RGCs). Using a simulation model, we also investigated the effects of migration and selection on the spatial structure of resistance phenotypes in a metapopulation for two genetic determinisms of the interaction (gene-for-gene and matching allele). Our results showed a differentiation between the countries for all the markers and indicated that the RGC, polymorphic in both countries, do not behave as neutral markers. Comparison of the diversities for resistance to strains isolated from wild or cultivated plants suggested that, although there is local adaptation of C. lindemuthianum between the two countries, the coevolution process seems to occur at a very local scale with the maintenance of resistances to allopatric strains, a result consistent with simulations of the models.

It has been brought to our attention that our paper (Kronholm et al. BMC Genetics 2010, 11:33) may have caused some confusion for readers interested in the correct quantification of population differentiation. We feel that this issue is of some importance and wish to clarify any confusion that might have resulted.

The selective impact of pathogen epidemics on host defenses can be strong but remains transient. By contrast, life-history shifts can durably and continuously modify the balance between costs and benefits of immunity, which arbitrates the evolution of host defenses. Their impact on the evolutionary dynamics of host immunity, however, has seldom been documented. Optimal investment into immunity is expected to decrease with shortening lifespan, because a shorter life decreases the probability to encounter pathogens or enemies. Here, we document that in natural populations of Arabidopsis thaliana, the expression levels of immunity genes correlate positively with flowering time, which in annual species is a proxy for lifespan. Using a novel genetic strategy based on bulk-segregants, we partitioned flowering time-dependent from – independent immunity genes and could demonstrate that this positive co-variation can be genetically separated. It is therefore not explained by the pleiotropic action of some major regulatory genes controlling both immunity and lifespan. Moreover, we find that immunity genes containing variants reported to impact fitness in natural field conditions are among the genes whose expression co-varies most strongly with flowering time. Taken together, these analyses reveal that natural selection has likely assorted alleles promoting lower expression of immunity genes with alleles that decrease the duration of vegetative lifespan in A. thaliana and vice versa. This is the first study documenting a pattern of variation consistent with the impact that selection on flowering time is predicted to have on diversity in host immunity.