Long-Term Balancing Selection and the Genetic Load Linked to the Self-Incompatibility Locus in Arabidopsis halleri and A. lyrata.

Balancing selection is a form of natural selection maintaining diversity at the sites it targets and at linked nucleotide sites. Due to selection favoring 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 species Arabidopsis halleri and Arabidopsis 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 the S-locus in all sample sets, but only over a limited genomic region, as it became indistinguishable from the genomic background beyond the first 25-30 kb. 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 the S-locus and clarify how natural selection in one genomic region affects the evolution of the adjacent genomic regions.

Polygenic adaptation of rosette growth in Arabidopsis thaliana.

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.

Approximate Bayesian Computation Untangles Signatures of Contemporary and Historical Hybridization between Two Endangered Species.

Contemporary gene flow, when resumed after a period of isolation, can have crucial consequences for endangered species, as it can both increase the supply of adaptive alleles and erode local adaptation. Determining the history of gene flow and thus the importance of contemporary hybridization, however, is notoriously difficult. Here, we focus on two endangered plant species, Arabis nemorensis and A. sagittata, which hybridize naturally in a sympatric population located on the banks of the Rhine. Using reduced genome sequencing, we determined the phylogeography of the two taxa but report only a unique sympatric population. Molecular variation in chloroplast DNA indicated that A. sagittata is the principal receiver of gene flow. Applying classical D-statistics and its derivatives to whole-genome data of 35 accessions, we detect gene flow not only in the sympatric population but also among allopatric populations. Using an Approximate Bayesian computation approach, we identify the model that best describes the history of gene flow between these taxa. This model shows that low levels of gene flow have persisted long after speciation. Around 10 000 years ago, gene flow stopped and a period of complete isolation began. Eventually, a hotspot of contemporary hybridization was formed in the unique sympatric population. Occasional sympatry may have helped protect these lineages from extinction in spite of their extremely low diversity.

Common gardens in teosintes reveal the establishment of a syndrome of adaptation to altitude.

In plants, local adaptation across species range is frequent. Yet, much has to be discovered on its environmental drivers, the underlying functional traits and their molecular determinants. Genome scans are popular to uncover outlier loci potentially involved in the genetic architecture of local adaptation, however links between outliers and phenotypic variation are rarely addressed. Here we focused on adaptation of teosinte populations along two elevation gradients in Mexico that display continuous environmental changes at a short geographical scale. We used two common gardens, and phenotyped 18 traits in 1664 plants from 11 populations of annual teosintes. In parallel, we genotyped these plants for 38 microsatellite markers as well as for 171 outlier single nucleotide polymorphisms (SNPs) that displayed excess of allele differentiation between pairs of lowland and highland populations and/or correlation with environmental variables. Our results revealed that phenotypic differentiation at 10 out of the 18 traits was driven by local selection. Trait covariation along the elevation gradient indicated that adaptation to altitude results from the assembly of multiple co-adapted traits into a complex syndrome: as elevation increases, plants flower earlier, produce less tillers, display lower stomata density and carry larger, longer and heavier grains. The proportion of outlier SNPs associating with phenotypic variation, however, largely depended on whether we considered a neutral structure with 5 genetic groups (73.7%) or 11 populations (13.5%), indicating that population stratification greatly affected our results. Finally, chromosomal inversions were enriched for both SNPs whose allele frequencies shifted along elevation as well as phenotypically-associated SNPs. Altogether, our results are consistent with the establishment of an altitudinal syndrome promoted by local selective forces in teosinte populations in spite of detectable gene flow. Because elevation mimics climate change through space, SNPs that we found underlying phenotypic variation at adaptive traits may be relevant for future maize breeding.

Rapid adaptive evolution to drought in a subset of plant traits in a large-scale climate change experiment.

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 10 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.

Strengths and potential pitfalls of hay transfer for ecological restoration revealed by RAD-seq analysis in floodplain Arabis species.

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 floodplain 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 make-up 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 prerestoration characterization of microgeographic patterns of intraspecific diversity of the community to guarantee that restoration practices reach their goal, that is 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.

Linking genes with ecological strategies in Arabidopsis thaliana.

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

Natural variation in stomata size contributes to the local adaptation of water-use efficiency in Arabidopsis thaliana.

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 (δ13 C). We detect substantial genetic variation for stomata size and density segregating within Arabidopsis thaliana. A positive correlation between stomata size and δ13 C 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 stomatal patterning but also to its covariation with carbon uptake parameters. Yet, we report two novel QTL affecting δ13 C independently of stomatal 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 covariation 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 favour the evolution of drought avoidance mechanisms over drought escape strategies.

The Footprint of Polygenic Adaptation on Stress-Responsive Cis-Regulatory Divergence in the Arabidopsis Genus.

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.

miR824-Regulated AGAMOUS-LIKE16 Contributes to Flowering Time Repression in Arabidopsis.

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.

Local evolution of seed flotation in Arabidopsis.

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.

Genome-wide association study of 107 phenotypes in Arabidopsis thaliana inbred lines.

Although pioneered by human geneticists as a potential solution to the challenging problem of finding the genetic basis of common human diseases, 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 traits. 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.

Genome-wide analysis of cis-regulatory divergence between species in the Arabidopsis genus.

Cis-regulatory DNA has been suspected to play a preeminent role in adaptive evolution, but understanding the role of cis-regulatory mutations in gene expression divergence first requires an accurate analysis of the functional differences associated with these regions. We analyzed allele-specific expression (ASE) in leaf and floral tissues of F1 interspecific hybrids generated between the two closely related species Arabidopsis thaliana and A. lyrata with a whole-genome SNP (single nucleotide polymorphism) tiling array. We observed 2,205 genes showing ASE pattern in at least one tissue. Nearly 90% of genes displaying ASE preferentially expressed the allele of A. lyrata. Genome-wide comparison of sequence divergence revealed that genes displaying ASE had a higher ratio of nonsynonymous to synonymous substitutions in coding regions. We further observe that the epigenetic landscape of histone methylation in A. thaliana genome associate with ASE. The asymmetry in the direction of allele-specific expression suggests interspecific differences in the efficiency of gene silencing in F1 hybrids.

Widespread interspecific divergence in cis-regulation of transposable elements in the Arabidopsis genus.

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.

Flagellin perception varies quantitatively in Arabidopsis thaliana and its relatives.

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.

Can we always sweep the details of RNA-processing under the carpet?

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.

DOG1 expression is predicted by the seed-maturation environment and contributes to geographical variation in germination in Arabidopsis thaliana.

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.

Structurally different alleles of the ath-MIR824 microRNA precursor are maintained at high frequency in Arabidopsis thaliana.

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.

Co-expression of neighbouring genes in Arabidopsis: separating chromatin effects from direct interactions.

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%.