African genomes illuminate the early history and transition to selfing in Arabidopsis thaliana

The Roles of Drift and Selection on Short Stamen Loss in Arabidopsis thaliana along an Elevational Gradient in the Spanish Pyrenees

Traits that have lost function sometimes persist through evolutionary time. These traits may be maintained by a lack of standing genetic variation for the trait, if selection against the trait is weak relative to drift, or if they have a residual function. To determine the evolutionary processes shaping whether nonfunctional traits are retained or lost, we investigated short stamens in 16 populations of Arabidopsis thaliana along an elevational cline in the Spanish Pyrenees. We found a cline in short stamen number from retention of short stamens in high elevation populations to incomplete loss in low elevation populations. We did not find evidence that limited genetic variation constrains the loss of short stamens at high elevations nor evidence for divergent selection on short stamens between high and low elevations. Finally, we identified loci associated with short stamens in the Spanish Pyrenees that are different from loci associated with variation in short stamen number across latitudes from a previous study. Overall, we did not identify the evolutionary mechanisms maintaining an elevational cline in short stamen number but did identify different genetic loci underlying the variation in short stamen along similar phenotypic clines.

Teaser text

The evolutionary mechanisms underlying loss or retention of traits that have lost function are poorly understood. Short stamens in Arabidopsis thaliana provide a compelling system to investigate the roles of genetic drift and selection in trait loss across latitudinal and elevational clines. This study investigates the role of drift and selection in short stamen loss in 16 populations of A. thaliana along an elevational gradient in the Spanish Pyrenees. An investigation of the genetic loci underlying variation in short stamen number suggests mutations in different genes may cause trait loss in similar phenotypic clines within a species.

Environment-induced heritable variations are common in Arabidopsis thaliana

Parental or ancestral environments can induce heritable phenotypic changes, but whether such environment-induced heritable changes are a common phenomenon remains unexplored. Here, we subject 14 genotypes of Arabidopsis thaliana to 10 different environmental treatments and observe phenotypic and genome-wide gene expression changes over four successive generations. We find that all treatments caused heritable phenotypic and gene expression changes, with a substantial proportion stably transmitted over all observed generations. Intriguingly, the susceptibility of a genotype to environmental inductions could be predicted based on the transposon abundance in the genome. Our study thus challenges the classic view that the environment only participates in the selection of heritable variation and suggests that the environment can play a significant role in generating of heritable variations.

Current status of community resources and priorities for weed genomics research

Weeds are attractive models for basic and applied research due to their impacts on agricultural systems and capacity to swiftly adapt in response to anthropogenic selection pressures. Currently, a lack of genomic information precludes research to elucidate the genetic basis of rapid adaptation for important traits like herbicide resistance and stress tolerance and the effect of evolutionary mechanisms on wild populations. The International Weed Genomics Consortium is a collaborative group of scientists focused on developing genomic resources to impact research into sustainable, effective weed control methods and to provide insights about stress tolerance and adaptation to assist crop breeding.

A pan-genome of 69 Arabidopsis thaliana accessions reveals a conserved genome structure throughout the global species range

Although originally primarily a system for functional biology, Arabidopsis thaliana has, owing to its broad geographical distribution and adaptation to diverse environments, developed into a powerful model in population genomics. Here we present chromosome-level genome assemblies of 69 accessions from a global species range. We found that genomic colinearity is very conserved, even among geographically and genetically distant accessions. Along chromosome arms, megabase-scale rearrangements are rare and typically present only in a single accession. This indicates that the karyotype is quasi-fixed and that rearrangements in chromosome arms are counter-selected. Centromeric regions display higher structural dynamics, and divergences in core centromeres account for most of the genome size variations. Pan-genome analyses uncovered 32,986 distinct gene families, 60% being present in all accessions and 40% appearing to be dispensable, including 18% private to a single accession, indicating unexplored genic diversity. These 69 new Arabidopsis thaliana genome assemblies will empower future genetic research.

Climate biogeography of Arabidopsis thaliana: linking distribution models and individual variation

AIM

Patterns of individual variation are key to testing hypotheses about the mechanisms underlying biogeographic patterns. If species distributions are determined by environmental constraints, then populations near range margins may have reduced performance and be adapted to harsher environments. Model organisms are potentially important systems for biogeographical studies, given the available range-wide natural history collections, and the importance of providing biogeographical context to their genetic and phenotypic diversity.

LOCATION

Global.

TAXON

Arabidopsis thaliana ("Arabidopsis").

METHODS

We fit occurrence records to climate data, and then projected the distribution of Arabidopsis under last glacial maximum, current, and future climates. We confronted model predictions with individual performance measured on 2,194 herbarium specimens, and we asked whether predicted suitability was associated with life-history and genomic variation measured on ~900 natural accessions.

RESULTS

The most important climate variables constraining the Arabidopsis distribution were winter cold in northern and high elevation regions and summer heat in southern regions. Herbarium specimens from regions with lower habitat suitability in both northern and southern regions were smaller, supporting the hypothesis that the distribution of Arabidopsis is constrained by climate-associated factors. Climate anomalies partly explained interannual variation in herbarium specimen size, but these did not closely correspond to local limiting factors identified in the distribution model. Late-flowering genotypes were absent from the lowest suitability regions, suggesting slower life histories are only viable closer to the center of the realized niche. We identified glacial refugia farther north than previously recognized, as well as refugia concordant with previous population genetic findings. Lower latitude populations, known to be genetically distinct, are most threatened by future climate change. The recently colonized range of Arabidopsis was well-predicted by our native-range model applied to certain regions but not others, suggesting it has colonized novel climates.

MAIN CONCLUSIONS

Integration of distribution models with performance data from vast natural history collections is a route forward for testing biogeographical hypotheses about species distributions and their relationship with evolutionary fitness across large scales.

Forces driving transposable element load variation during Arabidopsis range expansion

Genetic load refers to the accumulated and potentially life-threatening deleterious mutations in populations. Understanding the mechanisms underlying genetic load variation of transposable element (TE) insertion, a major large-effect mutation, during range expansion is an intriguing question in biology. Here, we used 1,115 global natural accessions of Arabidopsis (Arabidopsis thaliana) to study the driving forces of TE load variation during its range expansion. TE load increased with range expansion, especially in the recently established Yangtze River basin population. Effective population size, which explains 62.0% of the variance in TE load, high transposition rate, and selective sweeps contributed to TE accumulation in the expanded populations. We genetically mapped and identified multiple candidate causal genes and TEs, and revealed the genetic architecture of TE load variation. Overall, this study reveals the variation in TE genetic load during Arabidopsis expansion and highlights the causes of TE load variation from the perspectives of both population genetics and quantitative genetics.

Incorrect recombination partner associations contribute to meiotic instability of neo-allopolyploid Arabidopsis suecica

Combining two or more related homoeologous genomes in a single nucleus, newly formed allopolyploids must rapidly adapt meiosis to restore balanced chromosome segregation, production of euploid gametes and fertility. The poor fertility of such neo-allopolyploids thus strongly selects for the limitation or avoidance of genetic crossover formation between homoeologous chromosomes. In this study, we have reproduced the interspecific hybridization between Arabidopsis thaliana and Arabidopsis arenosa leading to the allotetraploid Arabidopsis suecica and have characterized the first allopolyploid meioses. First-generation neo-allopolyploid siblings vary considerably in fertility, meiotic behavior and levels of homoeologous recombination. We show that centromere dynamics at early meiosis is altered in synthetic neo-allopolyploids compared with evolved A. suecica, with a significant increase in homoeologous centromere interactions at zygotene. At metaphase I, the presence of multivalents involving homoeologous chromosomes confirms that homoeologous recombination occurs in the first-generation synthetic allopolyploid plants and this is associated with a significant reduction in homologous recombination, compared to evolved A. suecica. Together, these data strongly suggest that the fidelity of recombination partner choice, likely during the DNA invasion step, is strongly impaired during the first meiosis of neo-allopolyploids and requires subsequent adaptation.

Molecular insights into self-incompatibility systems: From evolution to breeding

Plants have evolved diverse self-incompatibility (SI) systems for outcrossing. Since Darwin's time, considerable progress has been made toward elucidating this unrivaled reproductive innovation. Recent advances in interdisciplinary studies and applications of biotechnology have given rise to major breakthroughs in understanding the molecular pathways that lead to SI, particularly the strikingly different SI mechanisms that operate in Solanaceae, Papaveraceae, Brassicaceae, and Primulaceae. These best-understood SI systems, together with discoveries in other "nonmodel" SI taxa such as Poaceae, suggest a complex evolutionary trajectory of SI, with multiple independent origins and frequent and irreversible losses. Extensive exploration of self-/nonself-discrimination signaling cascades has revealed a comprehensive catalog of male and female identity genes and modifier factors that control SI. These findings also enable the characterization, validation, and manipulation of SI-related factors for crop improvement, helping to address the challenges associated with development of inbred lines. Here, we review current knowledge about the evolution of SI systems, summarize key achievements in the molecular basis of pollen‒pistil interactions, discuss potential prospects for breeding of SI crops, and raise several unresolved questions that require further investigation.

A large-effect fitness trade-off across environments is explained by a single mutation affecting cold acclimation

Identifying the genetic basis of local adaptation and fitness trade-offs across environments is a central goal of evolutionary biology. Cold acclimation is an adaptive plastic response for surviving seasonal freezing, and costs of acclimation may be a general mechanism for fitness trade-offs across environments in temperate zone species. Starting with locally adapted ecotypes of Arabidopsis thaliana from Italy and Sweden, we examined the fitness consequences of a naturally occurring functional polymorphism in CBF2. This gene encodes a transcription factor that is a major regulator of cold-acclimated freezing tolerance and resides within a locus responsible for a genetic trade-off for long-term mean fitness. We estimated the consequences of alternate genotypes of CBF2 on 5-y mean fitness and fitness components at the native field sites by comparing near-isogenic lines with alternate genotypes of CBF2 to their genetic background ecotypes. The effects of CBF2 were validated at the nucleotide level using gene-edited lines in the native genetic backgrounds grown in simulated parental environments. The foreign CBF2 genotype in the local genetic background reduced long-term mean fitness in Sweden by more than 10%, primarily via effects on survival. In Italy, fitness was reduced by more than 20%, primarily via effects on fecundity. At both sites, the effects were temporally variable and much stronger in some years. The gene-edited lines confirmed that CBF2 encodes the causal variant underlying this genetic trade-off. Additionally, we demonstrated a substantial fitness cost of cold acclimation, which has broad implications for potential maladaptive responses to climate change.

The pan-genome and local adaptation of Arabidopsis thaliana

Arabidopsis thaliana serves as a model species for investigating various aspects of plant biology. However, the contribution of genomic structural variations (SVs) and their associate genes to the local adaptation of this widely distribute species remains unclear. Here, we de novo assemble chromosome-level genomes of 32 A. thaliana ecotypes and determine that variable genes expand the gene pool in different ecotypes and thus assist local adaptation. We develop a graph-based pan-genome and identify 61,332 SVs that overlap with 18,883 genes, some of which are highly involved in ecological adaptation of this species. For instance, we observe a specific 332 bp insertion in the promoter region of the HPCA1 gene in the Tibet-0 ecotype that enhances gene expression, thereby promotes adaptation to alpine environments. These findings augment our understanding of the molecular mechanisms underlying the local adaptation of A. thaliana across diverse habitats.

Genetic Architecture of Flowering Time Differs Between Populations With Contrasting Demographic and Selective Histories

Understanding the evolutionary factors that impact the genetic architecture of traits is a central goal of evolutionary genetics. Here, we investigate how quantitative trait variation accumulated over time in populations that colonized a novel environment. We compare the genetic architecture of flowering time in Arabidopsis populations from the drought-prone Cape Verde Islands and their closest outgroup population from North Africa. We find that trait polygenicity is severely reduced in the island populations compared to the continental North African population. Further, trait architectures and reconstructed allelic histories best fit a model of strong directional selection in the islands in accord with a Fisher-Orr adaptive walk. Consistent with this, we find that large-effect variants that disrupt major flowering time genes (FRI and FLC) arose first, followed by smaller effect variants, including ATX2 L125F, which is associated with a 4-day reduction in flowering time. The most recently arising flowering time-associated loci are not known to be directly involved in flowering time, consistent with an omnigenic signature developing as the population approaches its trait optimum. Surprisingly, we find no effect in the natural population of EDI-Cvi-0 (CRY2 V367M), an allele for which an effect was previously validated by introgression into a Eurasian line. Instead, our results suggest the previously observed effect of the EDI-Cvi-0 allele on flowering time likely depends on genetic background, due to an epistatic interaction. Altogether, our results provide an empirical example of the effects demographic history and selection has on trait architecture.

Ecological genetics of local adaptation in Arabidopsis: An 8-year field experiment

There is considerable evidence for local adaptation in nature, yet important questions remain regarding its genetic basis. How many loci are involved? What are their effect sizes? What is the relative importance of conditional neutrality versus genetic trade-offs? Here we address these questions in the self-pollinating, annual plant Arabidopsis thaliana. We used 400 recombinant inbred lines (RILs) derived from two locally adapted populations in Italy and Sweden, grew the RILs and parents at the parental locations, and mapped quantitative trait loci (QTL) for mean fitness (fruits/seedling planted). We previously published results from the first 3 years of the study, and here add five additional years, providing a unique opportunity to assess how temporal variation in selection might affect QTL detection and classification. We found 10 adaptive and one maladaptive QTL in Italy, and six adaptive and four maladaptive QTL in Sweden. The discovery of maladaptive QTL at both sites suggests that even locally adapted populations are not always at their genotypic optimum. Mean effect sizes for adaptive QTL, 0.97 and 0.55 fruits in Italy and Sweden, respectively, were large relative to the mean fitness of the RILs (approximately 8 fruits/seedling planted at both sites). Both genetic trade-offs (four cases) and conditional neutrality (seven cases) contribute to local adaptation in this system. The 8-year dataset provided greater power to detect QTL and to estimate their locations compared to our previous 3-year study, identifying one new genetic trade-off and resolving one genetic trade-off into two conditionally adaptive QTL.

Transition to Self-compatibility Associated With Dominant S-allele in a Diploid Siberian Progenitor of Allotetraploid Arabidopsis kamchatica Revealed by Arabidopsis lyrata Genomes

A transition to selfing can be beneficial when mating partners are scarce, for example, due to ploidy changes or at species range edges. Here, we explain how self-compatibility evolved in diploid Siberian Arabidopsis lyrata, and how it contributed to the establishment of allotetraploid Arabidopsis kamchatica. First, we provide chromosome-level genome assemblies for two self-fertilizing diploid A. lyrata accessions, one from North America and one from Siberia, including a fully assembled S-locus for the latter. We then propose a sequence of events leading to the loss of self-incompatibility in Siberian A. lyrata, date this independent transition to ∼90 Kya, and infer evolutionary relationships between Siberian and North American A. lyrata, showing an independent transition to selfing in Siberia. Finally, we provide evidence that this selfing Siberian A. lyrata lineage contributed to the formation of the allotetraploid A. kamchatica and propose that the selfing of the latter is mediated by the loss-of-function mutation in a dominant S-allele inherited from A. lyrata.

Pan-European study of genotypes and phenotypes in the Arabidopsis relative Cardamine hirsuta reveals how adaptation, demography, and development shape diversity patterns

We study natural DNA polymorphisms and associated phenotypes in the Arabidopsis relative Cardamine hirsuta. We observed strong genetic differentiation among several ancestry groups and broader distribution of Iberian relict strains in European C. hirsuta compared to Arabidopsis. We found synchronization between vegetative and reproductive development and a pervasive role for heterochronic pathways in shaping C. hirsuta natural variation. A single, fast-cycling ChFRIGIDA allele evolved adaptively allowing range expansion from glacial refugia, unlike Arabidopsis where multiple FRIGIDA haplotypes were involved. The Azores islands, where Arabidopsis is scarce, are a hotspot for C. hirsuta diversity. We identified a quantitative trait locus (QTL) in the heterochronic SPL9 transcription factor as a determinant of an Azorean morphotype. This QTL shows evidence for positive selection, and its distribution mirrors a climate gradient that broadly shaped the Azorean flora. Overall, we establish a framework to explore how the interplay of adaptation, demography, and development shaped diversity patterns of 2 related plant species.

Cycles of satellite and transposon evolution in Arabidopsis centromeres

Centromeres are critical for cell division, loading CENH3 or CENPA histone variant nucleosomes, directing kinetochore formation and allowing chromosome segregation^1,2. Despite their conserved function, centromere size and structure are diverse across species. To understand this centromere paradox^3,4, it is necessary to know how centromeric diversity is generated and whether it reflects ancient trans-species variation or, instead, rapid post-speciation divergence. To address these questions, we assembled 346 centromeres from 66 Arabidopsis thaliana and 2 Arabidopsis lyrata accessions, which exhibited a remarkable degree of intra- and inter-species diversity. A. thaliana centromere repeat arrays are embedded in linkage blocks, despite ongoing internal satellite turnover, consistent with roles for unidirectional gene conversion or unequal crossover between sister chromatids in sequence diversification. Additionally, centrophilic ATHILA transposons have recently invaded the satellite arrays. To counter ATHILA invasion, chromosome-specific bursts of satellite homogenization generate higher-order repeats and purge transposons, in line with cycles of repeat evolution. Centromeric sequence changes are even more extreme in comparison between A. thaliana and A. lyrata. Together, our findings identify rapid cycles of transposon invasion and purging through satellite homogenization, which drive centromere evolution and ultimately contribute to speciation.

Forecasting trait responses in novel environments to aid seed provenancing under climate change

Revegetation projects face the major challenge of sourcing optimal plant material. This is often done with limited information about plant performance and increasingly requires factoring resilience to climate change. Functional traits can be used as quantitative indices of plant performance and guide seed provenancing, but trait values expected under novel conditions are often unknown. To support climate-resilient provenancing efforts, we develop a trait prediction model that integrates the effect of genetic variation with fine-scale temperature variation. We train our model on multiple field plantings of Arabidopsis thaliana and predict two relevant fitness traits-days-to-bolting and fecundity-across the species' European range. Prediction accuracy was high for days-to-bolting and moderate for fecundity, with the majority of trait variation explained by temperature differences between plantings. Projection under future climate predicted a decline in fecundity, although this response was heterogeneous across the range. In response, we identified novel genotypes that could be introduced to genetically offset the fitness decay. Our study highlights the value of predictive models to aid seed provenancing and improve the success of revegetation projects.

Ancestral self-compatibility facilitates the establishment of allopolyploids in Brassicaceae

Self-incompatibility systems based on self-recognition evolved in hermaphroditic plants to maintain genetic variation of offspring and mitigate inbreeding depression. Despite these benefits in diploid plants, for polyploids who often face a scarcity of mating partners, self-incompatibility can thwart reproduction. In contrast, self-compatibility provides an immediate advantage: a route to reproductive viability. Thus, diploid selfing lineages may facilitate the formation of new allopolyploid species. Here, we describe the mechanism of establishment of at least four allopolyploid species in Brassicaceae (Arabidopsis suecica, Arabidopsis kamchatica, Capsella bursa-pastoris, and Brassica napus), in a manner dependent on the prior loss of the self-incompatibility mechanism in one of the ancestors. In each case, the degraded S-locus from one parental lineage was dominant over the functional S-locus of the outcrossing parental lineage. Such dominant loss-of-function mutations promote an immediate transition to selfing in allopolyploids and may facilitate their establishment.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Unravelling plant diversification: Intraspecific genetic differentiation in hybridizing Anacyclus species in the western Mediterranean Basin

PREMISE

The interfertile species Anacyclus clavatus, A. homogamos, and A. valentinus represent a plant complex coexisting in large anthropic areas of the western Mediterranean Basin with phenotypically mixed populations exhibiting a great floral variation. The goal of this study was to estimate the genetic identity of each species, to infer the role of hybridization in the observed phenotypic diversity, and to explore the effect of climate on the geographic distribution of species and genetic clusters.

METHODS

We used eight nuclear microsatellites to genotype 585 individuals from 31 populations of three Anacyclus species for population genetic analyses by using clustering algorithms based on Bayesian models and ordination methods. In addition, we used ecological niche models and niche overlap analyses for both the species and genetic clusters. We used an expanded data set, including 721 individuals from 129 populations for ecological niche models of the genetic clusters.

RESULTS

We found a clear correspondence between species and genetic clusters, except for A. clavatus that included up to three genetic clusters. We detected individuals with admixed genetic ancestry in A. clavatus and in mixed populations. Ecological niche models predicted similar distributions for species and genetic clusters. For the two specific genetic clusters of A. clavatus, ecological niche models predicted remarkably different areas.

CONCLUSIONS

Gene flow between Anacyclus species likely explains phenotypic diversity in contact areas. In addition, we suggest that introgression could be involved in the origin of one of the two A. clavatus genetic clusters, which also showed ecological differentiation.

Common evolutionary trajectory of short life-cycle in Brassicaceae ruderal weeds

Weed species are detrimental to crop yield. An understanding of how weeds originate and adapt to field environments is needed for successful crop management and reduction of herbicide use. Although early flowering is one of the weed trait syndromes that enable ruderal weeds to overcome frequent disturbances, the underlying genetic basis is poorly understood. Here, we establish Cardamine occulta as a model to study weed ruderality. By genome assembly and QTL mapping, we identify impairment of the vernalization response regulator gene FLC and a subsequent dominant mutation in the blue-light receptor gene CRY2 as genetic drivers for the establishment of short life cycle in ruderal weeds. Population genomics study further suggests that the mutations in these two genes enable individuals to overcome human disturbances through early deposition of seeds into the soil seed bank and quickly dominate local populations, thereby facilitating their spread in East China. Notably, functionally equivalent dominant mutations in CRY2 are shared by another weed species, Rorippa palustris, suggesting a common evolutionary trajectory of early flowering in ruderal weeds in Brassicaceae.

Differential environmental and genomic architectures shape the natural diversity for trichome patterning and morphology in different Arabidopsis organs

Despite the adaptive and taxonomic relevance of the natural diversity for trichome patterning and morphology, the molecular and evolutionary mechanisms underlying these traits remain mostly unknown, particularly in organs other than leaves. In this study, we address the ecological, genetic and molecular bases of the natural variation for trichome patterning and branching in multiple organs of Arabidopsis (Arabidopsis thaliana). To this end, we characterized a collection of 191 accessions and carried out environmental and genome-wide association (GWA) analyses. Trichome amount in different organs correlated negatively with precipitation in distinct seasons, thus suggesting a precise fit between trichome patterning and climate throughout the Arabidopsis life cycle. In addition, GWA analyses showed small overlapping between the genes associated with different organs, indicating partly independent genetic bases for vegetative and reproductive phases. These analyses identified a complex locus on chromosome 2, where two adjacent MYB genes (ETC2 and TCL1) displayed differential effects on trichome patterning in several organs. Furthermore, analyses of transgenic lines carrying different natural alleles demonstrated that TCL1 accounts for the variation for trichome patterning in all organs, and for stem trichome branching. By contrast, two other MYB genes (TRY and GL1), mainly showed effects on trichome patterning or branching, respectively.

Functional analysis of PHYB polymorphisms in Arabidopsis thaliana collected in Patagonia

Arabidopsis thaliana shows a wide range of natural genetic variation in light responses. Shade avoidance syndrome is a strategy of major adaptive significance that includes seed germination, elongation of vegetative structures, leaf hyponasty, and acceleration of flowering. Previously, we found that the southernmost Arabidopsis accession, collected in the south of Patagonia (Pat), is hyposensitive to light and displays a reduced response to shade light. This work aimed to explore the genetic basis of the shade avoidance response (SAR) for hypocotyl growth by QTL mapping in a recently developed 162 RIL population between Col-0 and Pat. We mapped four QTL for seedling hypocotyl growth: WL1 and WL2 QTL in white light, SHADE1 QTL in shade light, and SAR1 QTL for the SAR. PHYB is the strongest candidate gene for SAR1 QTL. Here we studied the function of two polymorphic indels in the promoter region, a GGGR deletion, and three non-synonymous polymorphisms on the PHYB coding region compared with the Col-0 reference genome. To decipher the contribution and relevance of each PHYB-Pat polymorphism, we constructed transgenic lines with single or double polymorphisms by using Col-0 as a reference genome. We found that single polymorphisms in the coding region of PHYB have discrete functions in seed germination, seedling development, and shade avoidance response. These results suggest distinct functions for each PHYB polymorphism to the adjustment of plant development to variable light conditions.

The selfing syndrome and beyond: diverse evolutionary consequences of mating system transitions in plants

The shift from outcrossing to self-fertilization (selfing) is considered one of the most prevalent evolutionary transitions in flowering plants. Selfing species tend to share similar reproductive traits in morphology and function, and such a set of traits is called the 'selfing syndrome'. Although the genetic basis of the selfing syndrome has been of great interest to evolutionary biologists, knowledge of the causative genes or mutations was limited until recently. Thanks to advances in population genomic methodologies combined with high-throughput sequencing technologies, several studies have successfully unravelled the molecular and genetic basis for evolution of the selfing syndrome in Capsella, Arabidopsis, Solanum and other genera. Here we first introduce recent research examples that have explored the loci, genes and mutations responsible for the selfing syndrome traits, such as reductions in petal size or in pollen production, that are mainly relevant to pre-pollination processes. Second, we review the relationship between the evolution of selfing and interspecific pollen transfer, highlighting the findings of post-pollination reproductive barriers at the molecular level. We then discuss the emerging view of patterns in evolution of the selfing syndrome, such as the pervasive involvement of loss-of-function mutations and the relative importance of selection versus neutral degradation. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

A two-step adaptive walk rewires nutrient transport in a challenging edaphic environment

Most well-characterized cases of adaptation involve single genetic loci. Theory suggests that multilocus adaptive walks should be common, but these are challenging to identify in natural populations. Here, we combine trait mapping with population genetic modeling to show that a two-step process rewired nutrient homeostasis in a population of Arabidopsis as it colonized the base of an active stratovolcano characterized by extremely low soil manganese (Mn). First, a variant that disrupted the primary iron (Fe) uptake transporter gene (IRT1) swept quickly to fixation in a hard selective sweep, increasing Mn but limiting Fe in the leaves. Second, multiple independent tandem duplications occurred at NRAMP1 and together rose to near fixation in the island population, compensating the loss of IRT1 by improving Fe homeostasis. This study provides a clear case of a multilocus adaptive walk and reveals how genetic variants reshaped a phenotype and spread over space and time.

Experimental demonstration and pan-structurome prediction of climate-associated riboSNitches in Arabidopsis

Background

Genome-wide association studies (GWAS) aim to correlate phenotypic changes with genotypic variation. Upon transcription, single nucleotide variants (SNVs) may alter mRNA structure, with potential impacts on transcript stability, macromolecular interactions, and translation. However, plant genomes have not been assessed for the presence of these structure-altering polymorphisms or “riboSNitches.”

Results

We experimentally demonstrate the presence of riboSNitches in transcripts of two Arabidopsis genes, ZINC RIBBON 3 (ZR3) and COTTON GOLGI-RELATED 3 (CGR3), which are associated with continentality and temperature variation in the natural environment. These riboSNitches are also associated with differences in the abundance of their respective transcripts, implying a role in regulating the gene's expression in adaptation to local climate conditions. We then computationally predict riboSNitches transcriptome-wide in mRNAs of 879 naturally inbred Arabidopsis accessions. We characterize correlations between SNPs/riboSNitches in these accessions and 434 climate descriptors of their local environments, suggesting a role of these variants in local adaptation. We integrate this information in CLIMtools V2.0 and provide a new web resource, T-CLIM, that reveals associations between transcript abundance variation and local environmental variation.

Conclusion

We functionally validate two plant riboSNitches and, for the first time, demonstrate riboSNitch conditionality dependent on temperature, coining the term “conditional riboSNitch.” We provide the first pan-genome-wide prediction of riboSNitches in plants. We expand our previous CLIMtools web resource with riboSNitch information and with 1868 additional Arabidopsis genomes and 269 additional climate conditions, which will greatly facilitate in silico studies of natural genetic variation, its phenotypic consequences, and its role in local adaptation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02656-4.

Parallel reduction in flowering time from de novo mutations enable evolutionary rescue in colonizing lineages

Understanding how populations adapt to abrupt environmental change is necessary to predict responses to future challenges, but identifying specific adaptive variants, quantifying their responses to selection and reconstructing their detailed histories is challenging in natural populations. Here, we use Arabidopsis from the Cape Verde Islands as a model to investigate the mechanisms of adaptation after a sudden shift to a more arid climate. We find genome-wide evidence of adaptation after a multivariate change in selection pressures. In particular, time to flowering is reduced in parallel across islands, substantially increasing fitness. This change is mediated by convergent de novo loss of function of two core flowering time genes: FRI on one island and FLC on the other. Evolutionary reconstructions reveal a case where expansion of the new populations coincided with the emergence and proliferation of these variants, consistent with models of rapid adaptation and evolutionary rescue.

Detailing how populations adapted to environmental change is needed to predict future responses, but identifying adaptive variants and detailing their fitness effects is rare. Here, the authors show that parallel loss of FRI and FLC function reduces time to flowering and drives adaptation in a drought prone environment.

Parent-of-Origin Effects on Seed Size Modify Heterosis Responses in Arabidopsis thaliana

Parent-of-origin effects arise when a phenotype depends on whether it is inherited maternally or paternally. Parent-of-origin effects can exert a strong influence on F1 seed size in flowering plants, an important agronomic and life-history trait that can contribute to biomass heterosis. Here we investigate the natural variation in the relative contributions of the maternal and paternal genomes to F1 seed size across 71 reciprocal pairs of F1 hybrid diploids and the parental effect on F1 seed size heterosis. We demonstrate that the paternally derived genome influences F1 seed size more significantly than previously appreciated. We further demonstrate (by disruption of parental genome dosage balance in F1 triploid seeds) that hybridity acts as an enhancer of genome dosage effects on F1 seed size, beyond that observed from hybridity or genome dosage effects on their own. Our findings indicate that interactions between genetic hybridity and parental genome dosage can enhance heterosis effects in plants, opening new avenues for boosting heterosis breeding in crop plants.

Phenotypic and Methylome Responses to Salt Stress in Arabidopsis thaliana Natural Accessions

Salt stress threatens plant growth, development and crop yields, and has become a critical global environmental issue. Increasing evidence has suggested that the epigenetic mechanism such as DNA methylation can mediate plant response to salt stress through transcriptional regulation and transposable element (TE) silencing. However, studies exploring genome-wide methylation dynamics under salt stress remain limited, in particular, for studies on multiple genotypes. Here, we adopted four natural accessions of the model species Arabidopsis thaliana and investigated the phenotypic and genome-wide methylation responses to salt stress through whole-genome bisulfite sequencing (WGBS). We found that salt stress significantly changed plant phenotypes, including plant height, rosette diameter, fruit number, and aboveground biomass, and the change in biomass tended to depend on accessions. Methylation analysis revealed that genome-wide methylation patterns depended primarily on accessions, and salt stress caused significant methylation changes in ∼ 0.1% cytosines over the genomes. About 33.5% of these salt-induced differential methylated cytosines (DMCs) were located to transposable elements (TEs). These salt-induced DMCs were mainly hypermethylated and accession-specific. TEs annotated to have DMCs (DMC-TEs) across accessions were found mostly belonged to the superfamily of Gypsy, a type II transposon, indicating a convergent DMC dynamic on TEs across different genetic backgrounds. Moreover, 8.0% of salt-induced DMCs were located in gene bodies and their proximal regulatory regions. These DMCs were also accession-specific, and genes annotated to have DMCs (DMC-genes) appeared to be more accession-specific than DMC-TEs. Intriguingly, both accession-specific DMC-genes and DMC-genes shared by multiple accessions were enriched in similar functions, including methylation, gene silencing, chemical homeostasis, polysaccharide catabolic process, and pathways relating to shifts between vegetative growth and reproduction. These results indicate that, across different genetic backgrounds, methylation changes may have convergent functions in post-transcriptional, physiological, and phenotypic modulation under salt stress. These convergent methylation dynamics across accession may be autonomous from genetic variation or due to convergent genetic changes, which requires further exploration. Our study provides a more comprehensive picture of genome-wide methylation dynamics under salt stress, and highlights the importance of exploring stress response mechanisms from diverse genetic backgrounds.

Seed dormancy varies widely among Arabidopsis thaliana populations both between and within Fennoscandia and Italy

The timing of germination is a key life‐history trait in plants, which is strongly affected by the strength of seed dormancy. Continental‐wide genetic variation in seed dormancy has been related to differences in climate and the timing of conditions suitable for seedling establishment. However, for predictions of adaptive potential and consequences of climatic change, information is needed regarding the extent to which seed dormancy varies within climatic regions and the factors driving such variation.

We quantified dormancy of seeds produced by 17 Italian and 28 Fennoscandian populations of Arabidopsis thaliana when grown in the greenhouse and at two field sites in Italy and Sweden. To identify possible drivers of among‐population variation in seed dormancy, we examined the relationship between seed dormancy and climate at the site of population origin, and between seed dormancy and flowering time.

Seed dormancy was on average stronger in the Italian compared to the Fennoscandian populations, but also varied widely within both regions. Estimates of seed dormancy in the three maternal environments were positively correlated. Among Fennoscandian populations, seed dormancy tended to increase with increasing summer temperature and decreasing precipitation at the site of population origin. In the smaller sample of Italian populations, no significant association was detected between mean seed dormancy and climate at the site of origin. The correlation between population mean seed dormancy and flowering time was weak and not statistically significant within regions.

The correlation between seed dormancy and climatic factors in Fennoscandia suggests that at least some of the among‐population variation is adaptive and that climate change will affect selection on this trait. The weak correlation between population mean seed dormancy and flowering time indicates that the two traits can evolve independently.

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.

Multiple Sources of Introduction of North American Arabidopsis thaliana from across Eurasia

Large-scale movement of organisms across their habitable range, or migration, is an important evolutionary process that can shape genetic diversity and influence the adaptive spread of alleles. Although human migrations have been studied in great detail with modern and ancient genomes, recent anthropogenic influence on reducing the biogeographical constraints on the migration of nonnative species has presented opportunities in several study systems to ask the questions about how repeated introductions shape genetic diversity in the introduced range. We present an extensive overview of population structure of North American Arabidopsis thaliana by studying a set of 500 whole-genome sequenced and over 2,800 RAD-seq genotyped individuals in the context of global diversity represented by Afro-Eurasian genomes. We use methods based on haplotype and rare-allele sharing as well as phylogenetic modeling to identify likely sources of introductions of extant N. American A. thaliana from the native range in Africa and Eurasia. We find evidence of admixture among the introduced lineages having increased haplotype diversity and reduced mutational load. We also detect signals of selection in immune-system-related genes that may impart qualitative disease resistance to pathogens of bacterial and oomycete origin. We conclude that multiple introductions to a nonnative range can rapidly enhance the adaptive potential of a colonizing species by increasing haplotypic diversity through admixture. Our results lay the foundation for further investigations into the functional significance of admixture.

Evidence for extensive hybridisation and past introgression events in feather grasses using genome-wide SNP genotyping

BACKGROUND

The proper identification of feather grasses in nature is often limited due to phenotypic variability and high morphological similarity between many species. Among plausible factors influencing this issue are hybridisation and introgression recently detected in the genus. Nonetheless, to date, only a bounded set of taxa have been investigated using integrative taxonomy combining morphological and molecular data. Here, we report the first large-scale study on five feather grass species across several hybrid zones in Russia and Central Asia. In total, 302 specimens were sampled in the field and classified based on the current descriptions of these taxa. They were then genotyped with high density genome-wide markers and measured based on a set of morphological characters to delimitate species and assess levels of hybridisation and introgression. Moreover, we tested species for past introgression and estimated divergence times between them.

RESULTS

Our findings demonstrated that 250 specimens represent five distinct species: S. baicalensis, S. capillata, S. glareosa, S. grandis and S. krylovii. The remaining 52 individuals provided evidence for extensive hybridisation between S. capillata and S. baicalensis, S. capillata and S. krylovii, S. baicalensis and S. krylovii, as well as to a lesser extent between S. grandis and S. krylovii, S. grandis and S. baicalensis. We detected past reticulation events between S. baicalensis, S. krylovii, S. grandis and inferred that diversification within species S. capillata, S. baicalensis, S. krylovii and S. grandis started ca. 130-96 kya. In addition, the assessment of genetic population structure revealed signs of contemporary gene flow between populations across species from the section Leiostipa, despite significant geographical distances between some of them. Lastly, we concluded that only 5 out of 52 hybrid taxa were properly identified solely based on morphology.

CONCLUSIONS

Our results support the hypothesis that hybridisation is an important mechanism driving evolution in Stipa. As an outcome, this phenomenon complicates identification of hybrid taxa in the field using morphological characters alone. Thus, integrative taxonomy seems to be the only reliable way to properly resolve the phylogenetic issue of Stipa. Moreover, we believe that feather grasses may be a suitable genus to study hybridisation and introgression events in nature.

Arabis alpina: A perennial model plant for ecological genomics and life-history evolution

Many model organisms were chosen and achieved prominence because of an advantageous combination of their life‐history characteristics, genetic properties and also practical considerations. Discoveries made in Arabidopsis thaliana, the most renowned noncrop plant model species, have markedly stimulated studies in other species with different biology. Within the family Brassicaceae, the arctic–alpine Arabis alpina has become a model complementary to Arabidopsis thaliana to study the evolution of life‐history traits, such as perenniality, and ecological genomics in harsh environments. In this review, we provide an overview of the properties that facilitated the rapid emergence of A. alpina as a plant model. We summarize the evolutionary history of A. alpina, including genomic aspects, the diversification of its mating system and demographic properties, and we discuss recent progress in the molecular dissection of developmental traits that are related to its perennial life history and environmental adaptation. From this published knowledge, we derive open questions that might inspire future research in A. alpina, other Brassicaceae species or more distantly related plant families.

The Evolutionary Dynamics of Genetic Incompatibilities Introduced by Duplicated Genes in Arabidopsis thaliana

Although gene duplications provide genetic backup and allow genomic changes under relaxed selection, they may potentially limit gene flow. When different copies of a duplicated gene are pseudofunctionalized in different genotypes, genetic incompatibilities can arise in their hybrid offspring. Although such cases have been reported after manual crosses, it remains unclear whether they occur in nature and how they affect natural populations. Here, we identified four duplicated-gene based incompatibilities including one previously not reported within an artificial Arabidopsis intercross population. Unexpectedly, however, for each of the genetic incompatibilities we also identified the incompatible alleles in natural populations based on the genomes of 1,135 Arabidopsis accessions published by the 1001 Genomes Project. Using the presence of incompatible allele combinations as phenotypes for GWAS, we mapped genomic regions that included additional gene copies which likely rescue the genetic incompatibility. Reconstructing the geographic origins and evolutionary trajectories of the individual alleles suggested that incompatible alleles frequently coexist, even in geographically closed regions, and that their effects can be overcome by additional gene copies collectively shaping the evolutionary dynamics of duplicated genes during population history.

Jumpstarting evolution: How transposition can facilitate adaptation to rapid environmental changes

Because of their ability to replicate across genomes, transposable elements (TEs) represent major generators of large-effect mutations. As a result, chromatin-based mechanisms have evolved to control the mutational potential of TEs at multiple levels, from the epigenetic silencing of TE sequences, through the modulation of their integration space, up to the alleviation of the impact of new insertions. Although most TE insertions are highly deleterious, some can provide key adaptive variation. Together with their remarkable sensitivity to the environment and precise integration preferences, the unique characteristics of TEs place them as potent genomic engines of adaptive innovation. Herein, we review recent works exploring the regulation and impact of transposition in nature and discuss their implications for the evolutionary response of species to drastic environmental changes.

Anno genominis XX: 20 years of Arabidopsis genomics

Twenty years ago, the Arabidopsis thaliana genome sequence was published. This was an important moment as it was the first sequenced plant genome and explicitly brought plant science into the genomics era. At the time, this was not only an outstanding technological achievement, but it was characterized by a superb global collaboration. The Arabidopsis genome was the seed for plant genomic research. Here, we review the development of numerous resources based on the genome that have enabled discoveries across plant species, which has enhanced our understanding of how plants function and interact with their environments.

MYB transcription factors drive evolutionary innovations in Arabidopsis fruit trichome patterning

Both inter- and intra-specific diversity has been described for trichome patterning in fruits, which is presumably involved in plant adaptation. However, the mechanisms underlying this developmental trait have been hardly addressed. Here we examined natural populations of Arabidopsis (Arabidopsis thaliana) that develop trichomes in fruits and pedicels, phenotypes previously not reported in the Arabidopsis genus. Genetic analyses identified five loci, MALAMBRUNO 1-5 (MAU1-5), with MAU2, MAU3, and MAU5 showing strong epistatic interactions that are necessary and sufficient to display these traits. Functional characterization of these three loci revealed cis-regulatory mutations in TRICHOMELESS1 and TRIPTYCHON, as well as a structural mutation in GLABRA1. Therefore, the multiple mechanisms controlled by three MYB transcription factors of the core regulatory network for trichome patterning have jointly been modulated to trigger trichome development in fruits. Furthermore, analyses of worldwide accessions showed that these traits and mutations only occur in a highly differentiated relict lineage from the Iberian Peninsula. In addition, these traits and alleles were associated with low spring precipitation, which suggests that trichome development in fruits and pedicels might be involved in climatic adaptation. Thus, we show that the combination of synergistic mutations in a gene regulatory circuit has driven evolutionary innovations in fruit trichome patterning in Arabidopsis.

Genetic variation, environment and demography intersect to shape Arabidopsis defense metabolite variation across Europe

Plants produce diverse metabolites to cope with the challenges presented by complex and ever-changing environments. These challenges drive the diversification of specialized metabolites within and between plant species. However, we are just beginning to understand how frequently new alleles arise controlling specialized metabolite diversity and how the geographic distribution of these alleles may be structured by ecological and demographic pressures. Here, we measure the variation in specialized metabolites across a population of 797 natural Arabidopsis thaliana accessions. We show that a combination of geography, environmental parameters, demography and different genetic processes all combine to influence the specific chemotypes and their distribution. This showed that causal loci in specialized metabolism contain frequent independently generated alleles with patterns suggesting potential within-species convergence. This provides a new perspective about the complexity of the selective forces and mechanisms that shape the generation and distribution of allelic variation that may influence local adaptation.

Since plants cannot move, they have evolved chemical defenses to help them respond to changes in their surroundings. For example, where animals run from predators, plants may produce toxins to put predators off. This approach is why plants are such a rich source of drugs, poisons, dyes and other useful substances. The chemicals plants produce are known as specialized metabolites, and they can change a lot between, and even within, plant species. The variety of specialized metabolites is a result of genetic changes and evolution over millions of years.

Evolution is a slow process, yet plants are able to rapidly develop new specialized metabolites to protect them from new threats. Even different populations of the same species produce many distinct metabolites that help them survive in their surroundings. However, the factors that lead plants to produce new metabolites are not well understood, and it is not known how this affects genetic variation.

To gain a better understanding of this process, Katz et al. studied 797 European variants of a common weed species called Arabidopsis thaliana, which is widely studied. The investigation found that many factors affect the range of specialized metabolites in each variant. These included local geography and environment, as well as genetics and population history (demography). Katz et al. revealed a pattern of relationships between the variants that could mirror their evolutionary history as the species spread and adapted to new locations.

These results highlight the complex network of factors that affect plant evolution. Rapid diversification is key to plant survival in new and changing environments and has resulted in a wide range of specialized metabolites. As such they are of interest both for studying plant evolution and for understanding their ecology. Expanding similar work to more populations and other species will broaden the scope of our ability to understand how plants adapt to their surroundings.

Maintenance of Adaptive Dynamics and No Detectable Load in a Range-Edge Outcrossing Plant Population

During range expansion, edge populations are expected to face increased genetic drift, which in turn can alter and potentially compromise adaptive dynamics, preventing the removal of deleterious mutations and slowing down adaptation. Here, we contrast populations of the European subspecies Arabidopsis lyrata ssp. petraea, which expanded its Northern range after the last glaciation. We document a sharp decline in effective population size in the range-edge population and observe that nonsynonymous variants segregate at higher frequencies. We detect a 4.9% excess of derived nonsynonymous variants per individual in the range-edge population, suggesting an increase of the genomic burden of deleterious mutations. Inference of the fitness effects of mutations and modeling of allele frequencies under the explicit demographic history of each population predicts a depletion of rare deleterious variants in the range-edge population, but an enrichment for fixed ones, consistent with the bottleneck effect. However, the demographic history of the range-edge population predicts a small net decrease in per-individual fitness. Consistent with this prediction, the range-edge population is not impaired in its growth and survival measured in a common garden experiment. We further observe that the allelic diversity at the self-incompatibility locus, which ensures strict outcrossing and evolves under negative frequency-dependent selection, has remained unchanged. Genomic footprints indicative of selective sweeps are broader in the Northern population but not less frequent. We conclude that the outcrossing species A. lyrata ssp. petraea shows a strong resilience to the effect of range expansion.

A New Catalog of Structural Variants in 1,301 A. thaliana Lines from Africa, Eurasia, and North America Reveals a Signature of Balancing Selection at Defense Response Genes

Genomic variation in the model plant Arabidopsis thaliana has been extensively used to understand evolutionary processes in natural populations, mainly focusing on single-nucleotide polymorphisms. Conversely, structural variation has been largely ignored in spite of its potential to dramatically affect phenotype. Here, we identify 155,440 indels and structural variants ranging in size from 1 bp to 10 kb, including presence/absence variants (PAVs), inversions, and tandem duplications in 1,301 A. thaliana natural accessions from Morocco, Madeira, Europe, Asia, and North America. We show evidence for strong purifying selection on PAVs in genes, in particular for housekeeping genes and homeobox genes, and we find that PAVs are concentrated in defense-related genes (R-genes, secondary metabolites) and F-box genes. This implies the presence of a "core" genome underlying basic cellular processes and a "flexible" genome that includes genes that may be important in spatially or temporally varying selection. Further, we find an excess of intermediate frequency PAVs in defense response genes in nearly all populations studied, consistent with a history of balancing selection on this class of genes. Finally, we find that PAVs in genes involved in the cold requirement for flowering (vernalization) and drought response are strongly associated with temperature at the sites of origin.

A story from the Miocene: Clock-dated phylogeny of Sisymbrium L. (Sisymbrieae, Brassicaceae)

Morphological variability and imprecise generic boundaries have hindered systematic, taxonomical, and nomenclatural studies of Sisymbrium L. (Brassicaceae, Sisymbrieae DC.). The members of this almost exclusively Old-World genus grow mostly on highly porous substrates across open steppe, semidesert, or ruderal habitats in the temperate zone of the Northern Hemisphere and African subtropics. The present study placed the biological history of Sisymbrium L. into time and space and rendered the tribus Sisymbrieae as monotypic. Five nuclear-encoded and three chloroplast-encoded loci of approximately 85% of all currently accepted species were investigated. Several accessions per species covering their whole distribution range allowed for a more representative assessment of intraspecific genetic diversity. In the light of fossil absence, the impact of different secondary calibration methods and taxon sets on time spans was tested, and we showed that such a combinatorial nested dating approach is beneficial. Multigene phylogeny accompanied with a time divergence estimation analysis placed the onset and development of this tribus into the western Irano-Turanian floristic region during the Miocene. Continuous increase in continentality and decrease in temperatures promoted the diversity of the Sisymbrieae, which invaded the open grasslands habitats in Eurasia, Mediterranean, and South Africa throughout the Pliocene and Pleistocene. Our results support the assumption of the Irano-Turanian region as a biodiversity reservoir for adjacent regions.

Review: Plant eco-evolutionary responses to climate change: Emerging directions

Contemporary climate change is exposing plant populations to novel combinations of temperatures, drought stress, [CO₂] and other abiotic and biotic conditions. These changes are rapidly disrupting the evolutionary dynamics of plants. Despite the multifactorial nature of climate change, most studies typically manipulate only one climatic factor. In this opinion piece, we explore how climate change factors interact with each other and with biotic pressures to alter evolutionary processes. We evaluate the ramifications of climate change across life history stages,and examine how mating system variation influences population persistence under rapid environmental change. Furthermore, we discuss how spatial and temporal mismatches between plants and their mutualists and antagonists could affect adaptive responses to climate change. For example, plant-virus interactions vary from highly pathogenic to mildly facilitative, and are partly mediated by temperature, moisture availability and [CO₂]. Will host plants exposed to novel, stressful abiotic conditions be more susceptible to viral pathogens? Finally, we propose novel experimental approaches that could illuminate how plants will cope with unprecedented global change, such as resurrection studies combined with experimental evolution, genomics or epigenetics.

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.

The rate at which plants grow is a major functional trait in plant ecology. However, little is known about its genetic variation in natural populations. Here, we investigate genetic and environmental factors shaping variation in the growth rate of Arabidopsis thaliana and ask whether genetic variation in plant growth contributes to adaptation to local environmental conditions. We grew plants under two light regimes that mimic latitudinal differences in the intensity of natural light radiation, and measured plant diameter as it grew over time. When the light supply was decreased, plant diameter grew more slowly but reached a markedly larger final size, confirming that plants can adjust their growth to prevailing light conditions. Yet, we also detected significant levels of genetic variation both in growth rate and in how the growth dynamics is adjusted to the light conditions. We show that this variation is encoded by many loci of small effect that are hard to locate in the genome but overall significantly enriched among genes associated with growth and stress reactions. We further observe that Spanish genotypes tended to reach, on average, a significantly larger rosette size than Northern European genotypes. Tests of adaptive divergence indicate that these differences may reflect adaptation to local environmental conditions.

Isolation, Library Preparation, and Bioinformatic Analysis of Historical and Ancient Plant DNA

The ability to sequence DNA retrieved from ancient and historical material plays a crucial role in reinforcing evolutionary and anthropological inference. While the focus of the field is largely on analyzing DNA from ancient hominids and other animals, we have also learned from plant ancient DNA (aDNA), in particular, about human farming practices, crop domestication, environment management, species invasion, and adaptation to various environmental conditions. In the following protocols, we outline best practices for plant aDNA isolation, preparation for sequencing, bioinformatic processing, and authentication. We describe the process all the way from processing of archaeological or historical plant material to characterizing and authenticating sequencing reads. In alternative protocols, we include modifications to this process that are tailored to strongly degraded DNA. Throughout, we stress the importance of precautionary measures to successfully analyze aDNA. Finally, we discuss the evolution of the archaeogenomics field and the development of new methods, which both shaped this protocol. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Isolation of aDNA Alternate Protocol 1: Isolation of ultra-short DNA (Dabney modification) Support Protocol 1: Preparation of PTB-based mix Support Protocol 2: Preparation of binding buffer Basic Protocol 2: Preparation of genomic libraries Alternate Protocol 2: Preparation of genomic libraries with uracil removal Basic Protocol 3: Bioinformatic processing and authentication of aDNA.

Strong stabilizing selection on timing of germination in a Mediterranean population of Arabidopsis thaliana

PREMISE

Timing of germination can strongly influence plant fitness by affecting seedling survival and by having cascading effects on later life-history traits. In seasonal environments, the period favorable for seedling establishment and growth is limited, and timing of germination is likely to be under stabilizing selection because of conflicting selection through survival and fecundity. Moreover, optimal germination time may vary among genotypes because of inherent differences in later life-history traits.

METHODS

To examine how germination time affects survival, fecundity, and the relative fitness of two genotypes differing in time to first flower, we conducted a field experiment in an Italian population of the winter annual Arabidopsis thaliana, in which seedling establishment occurs mainly in November. We transplanted seedlings of the local genotype and of a Swedish genotype monthly from August to December and monitored survival and fecundity.

RESULTS

Only seedlings transplanted in November and December survived until reproduction, and fitness of the November cohort was 35 times higher than that of the December cohort, indicating strong stabilizing selection on timing of germination. There was no evidence of conflicting selection: seedling survival, adult survival, and fecundity were all highest in the November cohort. Moreover, the relative fitness of the two genotypes did not differ significantly between cohorts.

CONCLUSIONS

The very narrow window of opportunity for seedling establishment was related to rapid seasonal changes in soil moisture and temperature, suggesting that rate of seasonal change is an important aspect to consider for understanding spatiotemporal variation in selection on phenological traits.

A Nested Association Mapping Panel in Arabidopsis thaliana for Mapping and Characterizing Genetic Architecture

Linkage and association mapping populations are crucial public resources that facilitate the characterization of trait genetic architecture in natural and agricultural systems. We define a large nested association mapping panel (NAM) from 14 publicly available recombinant inbred line populations (RILs) of Arabidopsis thaliana, which share a common recurrent parent (Col-0). Using a genotype-by-sequencing approach (GBS), we identified single nucleotide polymorphisms (SNPs; range 563-1525 per population) and subsequently built updated linkage maps in each of the 14 RIL sets. Simulations in individual RIL populations indicate that our GBS markers have improved power to detect small effect QTL and enhanced resolution of QTL support intervals in comparison to original linkage maps. Using these robust linkage maps, we imputed a common set of publicly available parental SNPs into each RIL linkage map, generating overlapping markers across all populations. Though ultimately depending on allele frequencies at causal loci, simulations of the NAM panel suggest that surveying between 4 to 7 of the 14 RIL populations provides high resolution of the genetic architecture of complex traits, relative to a single mapping population.

Genomic analysis reveals the genetic diversity, population structure, evolutionary history and relationships of Chinese pepper

Chinese pepper, mainly including Zanthoxylum bungeanum and Zanthoxylum armatum, is an economically important crop popular in Asian countries due to its unique taste characteristics and potential medical uses. Numerous cultivars of Chinese pepper have been developed in China through long-term domestication. To better understand the population structure, demographic history, and speciation of Chinese pepper, we performed a comprehensive analysis at a genome-wide level by analyzing 38,395 genomic SNPs that were identified in 112 cultivated and wild accessions using a high-throughput genome-wide genotyping-by-sequencing (GBS) approach. Our analysis provides genetic evidence of multiple splitting events occurring between and within species, resulting in at least four clades in Z. bungeanum and two clades in Z. armatum. Despite no evidence of recent admixture between species, we detected substantial gene flow within species. Estimates of demographic dynamics and species distribution modeling suggest that climatic oscillations during the Pleistocene (including the Penultimate Glaciation and the Last Glacial Maximum) and recent domestication events together shaped the demography and evolution of Chinese pepper. Our analyses also suggest that southeastern Gansu province is the most likely origin of Z. bungeanum in China. These findings provide comprehensive insights into genetic diversity, population structure, demography, and adaptation in Zanthoxylum.

Chromosome-Level Reference Genome and Population Genomic Analysis Provide Insights into the Evolution and Improvement of Domesticated Mulberry (Morus alba)

Mulberry (Morus spp.) is the sole plant consumed by the domesticated silkworm. However, the genome of domesticated mulberry has not yet been sequenced, and the ploidy level of this species remains unclear. Here, we report a high-quality, chromosome-level domesticated mulberry (Morus alba) genome. Analysis of genomic data and karyotype analyses confirmed that M. alba is a diploid with 28 chromosomes (2n = 2x = 28). Population genomic analysis based on resequencing of 134 mulberry accessions classified domesticated mulberry into three geographical groups, namely, Taihu Basin of southeastern China (Hu mulberry), northern and southwestern China, and Japan. Hu mulberry had the lowest nucleotide diversity among these accessions and demonstrated obvious signatures of selection associated with environmental adaptation. Further phylogenetic analysis supports a previous proposal that multiple domesticated mulberry accessions previously classified as different species actually belong to one species. This study expands our understanding of genome evolution of the genus Morus and population structure of domesticated mulberry, which would facilitate mulberry breeding and improvement.

Biased Gene Conversion Constrains Adaptation in Arabidopsis thaliana

Reduction of fitness due to deleterious mutations imposes a limit to adaptive evolution. By characterizing features that influence this genetic load we may better understand constraints on responses to both natural and human-mediated selection. Here, using whole-genome, transcriptome, and methylome data from >600 Arabidopsis thaliana individuals, we set out to identify important features influencing selective constraint. Our analyses reveal that multiple factors underlie the accumulation of maladaptive mutations, including gene expression level, gene network connectivity, and gene-body methylation. We then focus on a feature with major effect, nucleotide composition. The ancestral vs. derived status of segregating alleles suggests that GC-biased gene conversion, a recombination-associated process that increases the frequency of G and C nucleotides regardless of their fitness effects, shapes sequence patterns in A. thaliana Through estimation of mutational effects, we present evidence that biased gene conversion hinders the purging of deleterious mutations and contributes to a genome-wide signal of decreased efficacy of selection. By comparing these results to two outcrossing relatives, Arabidopsis lyrata and Capsella grandiflora, we find that protein evolution in A. thaliana is as strongly affected by biased gene conversion as in the outcrossing species. Last, we perform simulations to show that natural levels of outcrossing in A. thaliana are sufficient to facilitate biased gene conversion despite increased homozygosity due to selfing. Together, our results show that even predominantly selfing taxa are susceptible to biased gene conversion, suggesting that it may constitute an important constraint to adaptation among plant species.