The rate and potential relevance of new mutations in a colonizing plant lineage

Planting Genomes in the Wild: Arabidopsis from Genetics History to the Ecology and Evolutionary Genomics Era

The genetics model system Arabidopsis thaliana (L.) Heynh. lives across a vast geographic range with contrasting climates, in response to which it has evolved diverse life histories and phenotypic adaptations. In the last decade, the cataloging of worldwide populations, DNA sequencing of whole genomes, and conducting of outdoor field experiments have transformed it into a powerful evolutionary ecology system to understand the genomic basis of adaptation. Here, we summarize new insights on Arabidopsis following the coordinated efforts of the 1001 Genomes Project, the latest reconstruction of biogeographic and demographic history, and the systematic genomic mapping of trait natural variation through 15 years of genome-wide association studies. We then put this in the context of local adaptation across climates by summarizing insights from 73 Arabidopsis outdoor common garden experiments conducted to date. We conclude by highlighting how molecular and genomic knowledge of adaptation can help us to understand species' (mal)adaptation under ongoing climate change.

The genomic secrets of invasive plants

Genomics has revolutionised the study of invasive species, allowing evolutionary biologists to dissect mechanisms of invasion in unprecedented detail. Botanical research has played an important role in these advances, driving much of what we currently know about key determinants of invasion success (e.g. hybridisation, whole-genome duplication). Despite this, a comprehensive review of plant invasion genomics has been lacking. Here, we aim to address this gap, highlighting recent discoveries that have helped progress the field. For example, by leveraging genomics in natural and experimental populations, botanical research has confirmed the importance of large-effect standing variation during adaptation in invasive species. Further, genomic investigations of plants are increasingly revealing that large structural variants, as well as genetic changes induced by whole-genome duplication such as genomic redundancy or the breakdown of dosage-sensitive reproductive barriers, can play an important role during adaptive evolution of invaders. However, numerous questions remain, including when chromosomal inversions might help or hinder invasions, whether adaptive gene reuse is common during invasions, and whether epigenetically induced mutations can underpin the adaptive evolution of plasticity in invasive populations. We conclude by highlighting these and other outstanding questions that genomic studies of invasive plants are poised to help answer.

Museum genomics reveals temporal genetic stasis and global genetic diversity in Arabidopsis thaliana

Global patterns of population genetic variation through time offer a window into evolutionary processes that maintain diversity. Over time, lineages may expand or contract their distribution, causing turnover in population genetic composition. At individual loci, migration, drift, environmental change (among other processes) may affect allele frequencies. Museum specimens of widely distributed species offer a unique window into the genetics of understudied populations and changes over time. Here, we sequenced genomes of 130 herbarium specimens and 91 new field collections of Arabidopsis thaliana and combined these with published genomes. We sought a broader view of genomic diversity across the species, and to test if population genomic composition is changing through time. We documented extensive and previously uncharacterized diversity in a range of populations in Africa, populations that are under threat from anthropogenic climate change. Through time, we did not find dramatic changes in genomic composition of populations. Instead, we found a pattern of genetic change every 100 years of the same magnitude seen when comparing Eurasian populations that are 185 km apart, potentially due to a combination of drift and changing selection. We found only mixed signals of polygenic adaptation at phenology and physiology QTL. We did find that genes conserved across eudicots show altered levels of directional allele frequency change, potentially due to variable purifying and background selection. Our study highlights how museum specimens can reveal new dimensions of population diversity and show how wild populations are evolving in recent history.

TIPPo: A User-Friendly Tool for De Novo Assembly of Organellar Genomes with High-Fidelity Data

Plant cells have two major organelles with their own genomes: chloroplasts and mitochondria. While chloroplast genomes tend to be structurally conserved, the mitochondrial genomes of plants, which are much larger than those of animals, are characterized by complex structural variation. We introduce TIPPo, a user-friendly, reference-free assembly tool that uses PacBio high-fidelity long-read data and that does not rely on genomes from related species or nuclear genome information for the assembly of organellar genomes. TIPPo employs a deep learning model for initial read classification and leverages k-mer counting for further refinement, significantly reducing the impact of nuclear insertions of organellar DNA on the assembly process. We used TIPPo to completely assemble a set of 54 complete chloroplast genomes. No other tool was able to completely assemble this set. TIPPo is comparable with PMAT in assembling mitochondrial genomes from most species but does achieve even higher completeness for several species. We also used the assembled organelle genomes to identify instances of nuclear plastid DNA (NUPTs) and nuclear mitochondrial DNA (NUMTs) insertions. The cumulative length of NUPTs/NUMTs positively correlates with the size of the nuclear genome, suggesting that insertions occur stochastically. NUPTs/NUMTs show predominantly C:G to T:A changes, with the mutated cytosines typically found in CG and CHG contexts, suggesting that degradation of NUPT and NUMT sequences is driven by the known elevated mutation rate of methylated cytosines. Small interfering RNA loci are enriched in NUPTs and NUMTs, consistent with the RdDM pathway mediating DNA methylation in these sequences.

Atlas of telomeric repeat diversity in Arabidopsis thaliana

BACKGROUND

Telomeric repeat arrays at the ends of chromosomes are highly dynamic in composition, but their repetitive nature and technological limitations have made it difficult to assess their true variation in genome diversity surveys.

RESULTS

We have comprehensively characterized the sequence variation immediately adjacent to the canonical telomeric repeat arrays at the very ends of chromosomes in 74 genetically diverse Arabidopsis thaliana accessions. We first describe several types of distinct telomeric repeat units and then identify evolutionary processes such as local homogenization and higher-order repeat formation that shape diversity of chromosome ends. By comparing largely isogenic samples, we also determine repeat number variation of the degenerate and variant telomeric repeat array at both the germline and somatic levels. Finally, our analysis of haplotype structure uncovers chromosome end-specific patterns in the distribution of variant telomeric repeats, and their linkage to the more proximal non-coding region.

CONCLUSIONS

Our findings illustrate the spectrum of telomeric repeat variation at multiple levels in A. thaliana-in germline and soma, across all chromosome ends, and across genetic groups-thereby expanding our knowledge of the evolution of chromosome ends.

Fitness consequences of population bottlenecks in an invasive blowfly

Invasive species often undergo demographic bottlenecks that cause a decrease in genetic diversity and associated reductions in population fitness. Despite this, they manage to thrive in novel environments. Investigating the effects of inbreeding and genetic bottlenecks on population fitness for invasive species is, therefore, key to understanding how they may survive in new environments. We used the blowfly Calliphora vicina (Sciences, Mathématiques et Physique, 1830, 2, 1), which is native to Europe and was introduced to Australia and New Zealand, to examine the effects of genetic diversity on population fitness. We first collected 59 samples from 15 populations across New Zealand and one in Australia, and used 20,501 biallelic SNPs to investigate population genomic diversity, structure and admixture. We then explored the impacts of repeated experimental bottlenecks on population fitness by creating inbred and outbred lines of C. vicina and measuring a variety of fitness traits. In wild-caught samples, we found low overall genetic diversity, signals of genetic admixture and limited (<3%) genetic differentiation between North and South Island populations, with genetic links between the South Island and Australia. Following experimental bottlenecks, we found significant reductions in fitness for inbred lines. However, fitness effects were not felt equally across all phenotypic traits. Moreover, they were not enough to cause population collapse in any experimental line, suggesting that C. vicina (when under relaxed selection, as in laboratory settings) may be able to compensate for population bottlenecks even when highly inbred. Our results demonstrate the value of a tractable experimental system for investigating processes that may facilitate or hamper biological invasion.

Century-long timelines of herbarium genomes predict plant stomatal response to climate change

Dissecting plant responses to the environment is key to understanding whether and how plants adapt to anthropogenic climate change. Stomata, plants' pores for gas exchange, are expected to decrease in density following increased CO2 concentrations, a trend already observed in multiple plant species. However, it is unclear whether such responses are based on genetic changes and evolutionary adaptation. Here we make use of extensive knowledge of 43 genes in the stomatal development pathway and newly generated genome information of 191 Arabidopsis thaliana historical herbarium specimens collected over 193 years to directly link genetic variation with climate change. While we find that the essential transcription factors SPCH, MUTE and FAMA, central to stomatal development, are under strong evolutionary constraints, several regulators of stomatal development show signs of local adaptation in contemporary samples from different geographic regions. We then develop a functional score based on known effects of gene knock-out on stomatal development that recovers a classic pattern of stomatal density decrease over the past centuries, suggesting a genetic component contributing to this change. This approach combining historical genomics with functional experimental knowledge could allow further investigations of how different, even in historical samples unmeasurable, cellular plant phenotypes may have already responded to climate change through adaptive evolution.

Rapid phenotypic differentiation in the iconic Japanese knotweed s.l. invading novel habitats

Understanding the mechanisms that underlie plant invasions is critical for management and conservation of biodiversity. At the same time, invasive species also provide a unique opportunity to study rapid adaptation to complex environmental conditions. Using four replicate reciprocal transplant experiments across three habitats, we described patterns of phenotypic response and assessed the degree of local adaptation in knotweed populations. We found plants from beach habitats were generally smaller than plants from marsh and roadside habitats when grown in their home habitat. In the marsh habitat, marsh plants were generally larger than beach plants, but not different from roadside plants. There were no differences among plants grown in the roadside habitat. We found mixed evidence for local adaptation: plants from the marsh habitat had greater biomass in their "home" sites, while plants from beaches and roadsides had greater survival in their "home" sites compared to other plants. In sum, we found phenotypic differentiation and some support for the hypothesis of rapid local adaptation of plants from beach, marsh and roadside habitats. Identifying whether these patterns of differentiation result from genetic or heritable non-genetic mechanisms will require further work.

Predicting species invasiveness with genomic data: Is genomic offset related to establishment probability?

Predicting the risk of establishment and spread of populations outside their native range represents a major challenge in evolutionary biology. Various methods have recently been developed to estimate population (mal)adaptation to a new environment with genomic data via so-called Genomic Offset (GO) statistics. These approaches are particularly promising for studying invasive species but have still rarely been used in this context. Here, we evaluated the relationship between GO and the establishment probability of a population in a new environment using both in silico and empirical data. First, we designed invasion simulations to evaluate the ability to predict establishment probability of two GO computation methods (Geometric GO and Gradient Forest) under several conditions. Additionally, we aimed to evaluate the interpretability of absolute Geometric GO values, which theoretically represent the adaptive genetic distance between populations from distinct environments. Second, utilizing public empirical data from the crop pest species Bactrocera tryoni, a fruit fly native from Northern Australia, we computed GO between "source" populations and a diverse range of locations within invaded areas. This practical application of GO within the context of a biological invasion underscores its potential in providing insights and guiding recommendations for future invasion risk assessment. Overall, our results suggest that GO statistics represent good predictors of the establishment probability and may thus inform invasion risk, although the influence of several factors on prediction performance (e.g., propagule pressure or admixture) will need further investigation.

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.

Exploring the complex pre-adaptations of invasive plants to anthropogenic disturbance: a call for integration of archaeobotanical approaches

Pre-adaptation to anthropogenic disturbance is broadly considered key for plant invasion success. Nevertheless, empirical evidence remains scarce and fragmentary, given the multifaceted nature of anthropogenic disturbance itself and the complexity of other evolutionary forces shaping the (epi)-genomes of recent native and invasive plant populations. Here, we review and critically revisit the existing theory and empirical evidence in the field of evolutionary ecology and highlight novel integrative research avenues that work at the interface with archaeology to solve open questions. The approaches suggested so far focus on contemporary plant populations, although their genomes have rapidly changed since their initial introduction in response to numerous selective and stochastic forces. We elaborate that a role of pre-adaptation to anthropogenic disturbance in plant invasion success should thus additionally be validated based on the analyses of archaeobotanical remains. Such materials, in the light of detailed knowledge on past human societies could highlight fine-scale differences in the type and timing of past disturbances. We propose a combination of archaeobotanical, ancient DNA and morphometric analyses of plant macro- and microremains to assess past community composition, and species' functional traits to unravel the timing of adaptation processes, their drivers and their long-term consequences for invasive species. Although such methodologies have proven to be feasible for numerous crop plants, they have not been yet applied to wild invasive species, which opens a wide array of insights into their evolution.

Unveiling the evolutionary history of lingonberry (Vaccinium vitis-idaea L.) through genome sequencing and assembly of European and North American subspecies

Lingonberry (Vaccinium vitis-idaea L.) produces tiny red berries that are tart and nutty in flavor. It grows widely in the circumpolar region, including Scandinavia, northern parts of Eurasia, Alaska, and Canada. Although cultivation is currently limited, the plant has a long history of cultural use among indigenous communities. Given its potential as a food source, genomic resources for lingonberry are significantly lacking. To advance genomic knowledge, the genomes for 2 subspecies of lingonberry (V. vitis-idaea ssp. minus and ssp. vitis-idaea var. 'Red Candy') were sequenced and de novo assembled into contig-level assemblies. The assemblies were scaffolded using the bilberry genome (Vaccinium myrtillus) to generate a chromosome-anchored reference genome consisting of 12 chromosomes each with a total length of 548.07 Mb [contig N50 = 1.17 Mb, BUSCO (C%) = 96.5%] for ssp. vitis-idaea and 518.70 Mb [contig N50 = 1.40 Mb, BUSCO (C%) = 96.9%] for ssp. minus. RNA-seq-based gene annotation identified 27,243 and 25,718 genes on the respective assembly, and transposable element detection methods found that 45.82 and 44.58% of the genome were repeats. Phylogenetic analysis confirmed that lingonberry was most closely related to bilberry and was more closely related to blueberries than cranberries. Estimates of past effective population size suggested a continuous decline over the past 1-3 MYA, possibly due to the impacts of repeated glacial cycles during the Pleistocene leading to frequent population fragmentation. The genomic resource created in this study can be used to identify industry-relevant genes (e.g. anthocyanin production), infer phylogeny, and call sequence-level variants (e.g. SNPs) in future research.

Genomics for monitoring and understanding species responses to global climate change

All life forms across the globe are experiencing drastic changes in environmental conditions as a result of global climate change. These environmental changes are happening rapidly, incur substantial socioeconomic costs, pose threats to biodiversity and diminish a species' potential to adapt to future environments. Understanding and monitoring how organisms respond to human-driven climate change is therefore a major priority for the conservation of biodiversity in a rapidly changing environment. Recent developments in genomic, transcriptomic and epigenomic technologies are enabling unprecedented insights into the evolutionary processes and molecular bases of adaptation. This Review summarizes methods that apply and integrate omics tools to experimentally investigate, monitor and predict how species and communities in the wild cope with global climate change, which is by genetically adapting to new environmental conditions, through range shifts or through phenotypic plasticity. We identify advantages and limitations of each method and discuss future research avenues that would improve our understanding of species' evolutionary responses to global climate change, highlighting the need for holistic, multi-omics approaches to ecosystem monitoring during global climate change.

Genomic Tools in Biological Invasions: Current State and Future Frontiers

Human activities are accelerating rates of biological invasions and climate-driven range expansions globally, yet we understand little of how genomic processes facilitate the invasion process. Although most of the literature has focused on underlying phenotypic correlates of invasiveness, advances in genomic technologies are showing a strong link between genomic variation and invasion success. Here, we consider the ability of genomic tools and technologies to (i) inform mechanistic understanding of biological invasions and (ii) solve real-world issues in predicting and managing biological invasions. For both, we examine the current state of the field and discuss how genomics can be leveraged in the future. In addition, we make recommendations pertinent to broader research issues, such as data sovereignty, metadata standards, collaboration, and science communication best practices that will require concerted efforts from the global invasion genomics community.

Two chromosome-level genome assemblies of Rhodiola shed new light on genome evolution in rapid radiation and evolution of the biosynthetic pathway of salidroside

Rhodiola L. is a genus that has undergone rapid radiation in the mid-Miocene and may represent a typic case of adaptive radiation. Many species of Rhodiola have also been widely used as an important adaptogen in traditional medicines for centuries. However, a lack of high-quality chromosome-level genomes hinders in-depth study of its evolution and biosynthetic pathway of secondary metabolites. Here, we assembled two chromosome-level genomes for two Rhodiola species with different chromosome number and sexual system. The assembled genome size of R. chrysanthemifolia (2n = 14; hermaphrodite) and R. kirilowii (2n = 22; dioecious) were of 402.67 and 653.62 Mb, respectively, with approximately 57.60% and 69.22% of transposable elements (TEs). The size difference between the two genomes was mostly due to proliferation of long terminal repeat-retrotransposons (LTR-RTs) in the R. kirilowii genome. Comparative genomic analysis revealed possible gene families responsible for high-altitude adaptation of Rhodiola, including a homolog of plant cysteine oxidase 2 gene of Arabidopsis thaliana (AtPCO2), which is part of the core molecular reaction to hypoxia and contributes to the stability of Group VII ethylene response factors (ERF-VII). We found extensive chromosome fusion/fission events and structural variations between the two genomes, which might have facilitated the initial rapid radiation of Rhodiola. We also identified candidate genes in the biosynthetic pathway of salidroside. Overall, our results provide important insights into genome evolution in plant rapid radiations, and possible roles of chromosome fusion/fission and structure variation played in rapid speciation.

Temporal collections to study invasion biology

Biological invasions represent an extraordinary opportunity to study evolution. This is because accidental or deliberate species introductions have taken place for centuries across large geographical scales, frequently prompting rapid evolutionary transitions in invasive populations. Until recently, however, the utility of invasions as evolutionary experiments has been hampered by limited information on the makeup of populations that were part of earlier invasion stages. Now, developments in ancient and historical DNA technologies, as well as the quickening pace of digitization for millions of specimens that are housed in herbaria and museums globally, promise to help overcome this obstacle. In this review, we first introduce the types of temporal data that can be used to study invasions, highlighting the timescale captured by each approach and their respective limitations. We then discuss how ancient and historical specimens as well as data available from prior invasion studies can be used to answer questions on mechanisms of (mal)adaptation, rates of evolution, or community-level changes during invasions. By bridging the gap between contemporary and historical invasive populations, temporal data can help us connect pattern to process in invasion science. These data will become increasingly important if invasions are to achieve their full potential as experiments of evolution in nature.

Retrieval of long DNA reads from herbarium specimens

High-throughput sequencing of herbarium specimens' DNA with short-read platforms has helped explore many biological questions. Here, for the first time, we investigate the potential of using herbarium specimens as a resource for long-read DNA sequencing technologies. We use target capture of 48 low-copy nuclear loci in 12 herbarium specimens of Silene as a basis for long-read sequencing using SMRT PacBio Sequel. The samples were collected between 1932 and 2019. A simple optimization of size selection protocol enabled the retrieval of both long DNA fragments (>1 kb) and long on-target reads for nine of them. The limited sampling size does not enable statistical evaluation of the influence of specimen age to the DNA fragmentation, but our results confirm that younger samples, that is, collected after 1990, are less fragmented and have better sequencing success than specimens collected before this date. Specimens collected between 1990 and 2019 yield between 167 and 3403 on-target reads > 1 kb. They enabled recovering between 34 loci and 48 (i.e. all loci recovered). Three samples from specimens collected before 1990 did not yield on-target reads > 1 kb. The four other samples collected before this date yielded up to 144 reads and recovered up to 25 loci. Young herbarium specimens seem promising for long-read sequencing. However, older ones have partly failed. Further exploration would be necessary to statistically test and understand the potential of older material in the quest for long reads. We would encourage greatly expanding the sampling size and comparing different taxonomic groups.

Microbe-associated molecular pattern recognition receptors have little effect on endophytic Arabidopsis thaliana microbiome assembly in the field

Plant microbiome structure affects plant health and productivity. A limited subset of environmental microbes successfully establishes within plant tissues, but the forces underlying this selectivity remain poorly characterized. Transmembrane pattern recognition receptors (PRRs), used by plants to detect microbe-associated molecular patterns (MAMPs), are strong candidates for achieving this selectivity because PRRs can potentially interact with many members of the microbiome. Indeed, MAMPs found in many microbial taxa, including beneficials and commensals, can instigate a robust immune response that affects microbial growth. Surprisingly, we found that MAMP-detecting PRRs have little effect on endophytic bacterial and fungal microbiome structure in the field. We compared the microbiomes of four PRR knockout lines of Arabidopsis thaliana to wild-type plants in multiple tissue types over several developmental stages and detected only subtle shifts in fungal, but not bacterial, β-diversity in one of the four PRR mutants. In one developmental stage, lore mutants had slightly altered fungal β-diversity, indicating that LORE may be involved in plant-fungal interactions in addition to its known role in detecting certain bacterial lipids. No other effects of PRRs on α-diversity, microbiome variability, within-individual homogeneity, or microbial load were found. The general lack of effect suggests that individual MAMP-detecting PRRs are not critical in shaping the endophytic plant microbiome. Rather, we suggest that MAMP-detecting PRRs must either act in concert and/or are individually maintained through pleiotropic effects or interactions with coevolved mutualists or pathogens. Although unexpected, these results offer insights into the role of MAMP-detecting PRRs in plant-microbe interactions and help direct future efforts to uncover host genetic elements that control plant microbiome assembly.

Ancient DNA genomics and the renaissance of herbaria

Herbaria are undergoing a renaissance as valuable sources of genomic data for exploring plant evolution, ecology, and diversity. Ancient DNA retrieved from herbarium specimens can provide unprecedented glimpses into past plant communities, their interactions with biotic and abiotic factors, and the genetic changes that have occurred over time. Here, we highlight recent advances in the field of herbarium genomics and discuss the challenges and opportunities of combining data from modern and time-stamped historical specimens. We also describe how integrating herbarium genomics data with other data types can yield substantial insights into the evolutionary and ecological processes that shape plant communities. Herbarium genomic analysis is a tool for understanding plant life and informing conservation efforts in the face of dire environmental challenges.

An evolutionary epigenetic clock in plants

Molecular clocks are the basis for dating the divergence between lineages over macroevolutionary timescales (~105 to 108 years). However, classical DNA-based clocks tick too slowly to inform us about the recent past. Here, we demonstrate that stochastic DNA methylation changes at a subset of cytosines in plant genomes display a clocklike behavior. This "epimutation clock" is orders of magnitude faster than DNA-based clocks and enables phylogenetic explorations on a scale of years to centuries. We show experimentally that epimutation clocks recapitulate known topologies and branching times of intraspecies phylogenetic trees in the self-fertilizing plant Arabidopsis thaliana and the clonal seagrass Zostera marina, which represent two major modes of plant reproduction. This discovery will open new possibilities for high-resolution temporal studies of plant biodiversity.

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.

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.

Characterization of nuclear DNA diversity in an individual Leymus chinensis

Intraorganismal genetic heterogeneity (IGH) exists when an individual organism harbors more than one genotype among its cells. In general, intercellular DNA diversity occurs at a very low frequency and cannot be directly detected by DNA sequencing from bulk tissue. In this study, based on Sanger and high-throughput sequencing, different species, different organs, different DNA segments and a single cell were employed to characterize nucleotide mutations in Leymus chinensis. The results demonstrated that 1) the nuclear DNA showed excessive genetic heterogeneity among cells of an individual leaf or seed but the chloroplast genes remained consistent; 2) a high density of SNPs was found in the variants of the unique DNA sequence, and the similar SNP profile shared between the leaf and seed suggested that nucleotide mutation followed a certain rule and was not random; and 3) the mutation rate decreased from the genomic DNA sequence to the corresponding protein sequence. Our results suggested that Leymus chinensis seemed to consist of a collection of cells with different genetic backgrounds.

An evolutionary epigenetic clock in plants

Molecular clocks are the basis for dating the divergence between lineages over macro-evolutionary timescales (~10 5 -10 8 years). However, classical DNA-based clocks tick too slowly to inform us about the recent past. Here, we demonstrate that stochastic DNA methylation changes at a subset of cytosines in plant genomes possess a clock-like behavior. This 'epimutation-clock' is orders of magnitude faster than DNA-based clocks and enables phylogenetic explorations on a scale of years to centuries. We show experimentally that epimutation-clocks recapitulate known topologies and branching times of intra-species phylogenetic trees in the selfing plant A. thaliana and the clonal seagrass Z. marina , which represent two major modes of plant reproduction. This discovery will open new possibilities for high-resolution temporal studies of plant biodiversity.

Differences in pseudogene evolution contributed to the contrasting flavors of turnip and Chiifu, two Brassica rapa subspecies

Pseudogenes are important resources for investigation of genome evolution and genomic diversity because they are nonfunctional but have regulatory effects that influence plant adaptation and diversification. However, few systematic comparative analyses of pseudogenes in closely related species have been conducted. Here, we present a turnip (Brassica rapa ssp. rapa) genome sequence and characterize pseudogenes among diploid Brassica species/subspecies. The results revealed that the number of pseudogenes was greatest in Brassica oleracea (CC genome), followed by B. rapa (AA genome) and then Brassica nigra (BB genome), implying that pseudogene differences emerged after species differentiation. In Brassica AA genomes, pseudogenes were distributed asymmetrically on chromosomes because of numerous chromosomal insertions/rearrangements, which contributed to the diversity among subspecies. Pseudogene differences among subspecies were reflected in the flavor-related glucosinolate (GSL) pathway. Specifically, turnip had the highest content of pungent substances, probably because of expansion of the methylthioalkylmalate synthase-encoding gene family in turnips; these genes were converted into pseudogenes in B. rapa ssp. pekinensis (Chiifu). RNA interference-based silencing of the gene encoding 2-oxoglutarate-dependent dioxygenase 2, which is also associated with flavor and anticancer substances in the GSL pathway, resulted in increased abundance of anticancer compounds and decreased pungency of turnip and Chiifu. These findings revealed that pseudogene differences between turnip and Chiifu influenced the evolution of flavor-associated GSL metabolism-related genes, ultimately resulting in the different flavors of turnip and Chiifu.

The herbarium of the future

The ~400 million specimens deposited across ~3000 herbaria are essential for: (i) understanding where plants have lived in the past, (ii) forecasting where they may live in the future, and (iii) delineating their conservation status. An open access 'global metaherbarium' is emerging as these specimens are digitized, mobilized, and interlinked online. This virtual biodiversity resource is attracting new users who are accelerating traditional applications of herbaria and generating basic and applied scientific innovations, including e-monographs and floras produced by diverse, interdisciplinary, and inclusive teams; robust machine-learning algorithms for species identification and phenotyping; collection and synthesis of ecological trait data at large spatiotemporal and phylogenetic scales; and exhibitions and installations that convey the beauty of plants and the value of herbaria in addressing broader societal issues.

Population genetics reveal multiple independent invasions of Spodoptera frugiperda (Lepidoptera: Noctuidae) in China

The fall armyworm (Spodoptera frugiperda), a destructive pest that originated in South and North America, spread to China in early 2019. Controlling this invasive pest requires an understanding of its population structure and migration patterns, yet the invasion genetics of Chinese S. frugiperda is not clear. Here, using the mitochondrial cytochrome oxidase subunit I (COI) gene, triose phosphate isomerase (Tpi) gene and eight microsatellite loci, we investigated genetic structure and genetic diversity of 16 S. frugiperda populations in China. The Tpi locus identified most S. frugiperda populations as the corn-strains, and a few were heterozygous strains. The microsatellite loci revealed that the genetic diversity of this pest in China was lower than that in South America. Furthermore, we found moderate differentiation among the populations, distinct genetic structures between adjacent populations and abundant genetic resources in the S. frugiperda populations from China sampled across 2 years. The survival rate of S. frugiperda was significantly higher when it was fed on corn leaves than on rice leaves, and the larval stage mortality rate was the highest under both treatments. Our results showed that S. frugiperda probably invaded China via multiple independent introductions and careful pesticide control, continuous monitoring and further studies will be needed to minimize its potential future outbreak.

Genetic diversity loss in the Anthropocene

Anthropogenic habitat loss and climate change are reducing species' geographic ranges, increasing extinction risk and losses of species' genetic diversity. Although preserving genetic diversity is key to maintaining species' adaptability, we lack predictive tools and global estimates of genetic diversity loss across ecosystems. We introduce a mathematical framework that bridges biodiversity theory and population genetics to understand the loss of naturally occurring DNA mutations with decreasing habitat. By analyzing genomic variation of 10,095 georeferenced individuals from 20 plant and animal species, we show that genome-wide diversity follows a mutations-area relationship power law with geographic area, which can predict genetic diversity loss from local population extinctions. We estimate that more than 10% of genetic diversity may already be lost for many threatened and nonthreatened species, surpassing the United Nations' post-2020 targets for genetic preservation.

Uncovering natural variation in root system architecture and growth dynamics using a robotics-assisted phenomics platform

The plant kingdom contains a stunning array of complex morphologies easily observed above-ground, but more challenging to visualize below-ground. Understanding the magnitude of diversity in root distribution within the soil, termed root system architecture (RSA), is fundamental in determining how this trait contributes to species adaptation in local environments. Roots are the interface between the soil environment and the shoot system and therefore play a key role in anchorage, resource uptake, and stress resilience. Previously, we presented the GLO-Roots (Growth and Luminescence Observatory for Roots) system to study the RSA of soil-grown Arabidopsis thaliana plants from germination to maturity (Rellán-Álvarez et al., 2015). In this study, we present the automation of GLO-Roots using robotics and the development of image analysis pipelines in order to examine the temporal dynamic regulation of RSA and the broader natural variation of RSA in Arabidopsis, over time. These datasets describe the developmental dynamics of two independent panels of accessions and reveal highly complex and polygenic RSA traits that show significant correlation with climate variables of the accessions' respective origins.

Locally adaptive temperature response of vegetative growth in Arabidopsis thaliana

We investigated early vegetative growth of natural Arabidopsis thaliana accessions in cold, nonfreezing temperatures, similar to temperatures these plants naturally encounter in fall at northern latitudes. We found that accessions from northern latitudes produced larger seedlings than accessions from southern latitudes, partly as a result of larger seed size. However, their subsequent vegetative growth when exposed to colder temperatures was slower. The difference was too large to be explained by random population differentiation, and is thus suggestive of local adaptation, a notion that is further supported by substantial transcriptome and metabolome changes in northern accessions. We hypothesize that the reduced growth of northern accessions is an adaptive response and a consequence of reallocating resources toward cold acclimation and winter survival.

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.

Mutation bias reflects natural selection in Arabidopsis thaliana

Since the first half of the twentieth century, evolutionary theory has been dominated by the idea that mutations occur randomly with respect to their consequences1. Here we test this assumption with large surveys of de novo mutations in the plant Arabidopsis thaliana. In contrast to expectations, we find that mutations occur less often in functionally constrained regions of the genome-mutation frequency is reduced by half inside gene bodies and by two-thirds in essential genes. With independent genomic mutation datasets, including from the largest Arabidopsis mutation accumulation experiment conducted to date, we demonstrate that epigenomic and physical features explain over 90% of variance in the genome-wide pattern of mutation bias surrounding genes. Observed mutation frequencies around genes in turn accurately predict patterns of genetic polymorphisms in natural Arabidopsis accessions (r = 0.96). That mutation bias is the primary force behind patterns of sequence evolution around genes in natural accessions is supported by analyses of allele frequencies. Finally, we find that genes subject to stronger purifying selection have a lower mutation rate. We conclude that epigenome-associated mutation bias2 reduces the occurrence of deleterious mutations in Arabidopsis, challenging the prevailing paradigm that mutation is a directionless force in evolution.

Gene targeting in polymerase theta-deficient Arabidopsis thaliana

Agrobacterium tumefaciens-mediated transformation has been for decades the preferred tool to generate transgenic plants. During this process, a T-DNA carrying transgenes is transferred from the bacterium to plant cells, where it randomly integrates into the genome via polymerase theta (Polθ)-mediated end joining (TMEJ). Targeting of the T-DNA to a specific genomic locus via homologous recombination (HR) is also possible, but such gene targeting (GT) events occur at low frequency and are almost invariably accompanied by random integration events. An additional complexity is that the product of recombination between T-DNA and target locus may not only map to the target locus (true GT), but also to random positions in the genome (ectopic GT). In this study, we have investigated how TMEJ functionality affects the biology of GT in plants, by using Arabidopsis thaliana mutated for the TEBICHI gene, which encodes for Polθ. Whereas in TMEJ-proficient plants we predominantly found GT events accompanied by random T-DNA integrations, GT events obtained in the teb mutant background lacked additional T-DNA copies, corroborating the essential role of Polθ in T-DNA integration. Polθ deficiency also prevented ectopic GT events, suggesting that the sequence of events leading up to this outcome requires TMEJ. Our findings provide insights that can be used for the development of strategies to obtain high-quality GT events in crop plants.

Is there still evolution in the human population?

It is often claimed that humanity has stopped evolving because modern medicine erased all selection on survival. Even if that would be true, and it is not, there would be other mechanisms of evolution which could still led to changes in allelic frequencies. Here I show, by applying basic evolutionary genetics knowledge, that we expect humanity to evolve. The results from genome sequencing projects have repeatedly affirmed that there are still recent signs of selection in our genomes. I give some examples of such adaptation. Then I briefly discuss what our evolutionary future has in store for us.

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.

The Genomic Processes of Biological Invasions: From Invasive Species to Cancer Metastases and Back Again

The concept of invasion is useful across a broad range of contexts, spanning from the fine scale landscape of cancer tumors up to the broader landscape of ecosystems. Invasion biology provides extraordinary opportunities for studying the mechanistic basis of contemporary evolution at the molecular level. Although the field of invasion genetics was established in ecology and evolution more than 50 years ago, there is still a limited understanding of how genomic level processes translate into invasive phenotypes across different taxa in response to complex environmental conditions. This is largely because the study of most invasive species is limited by information about complex genome level processes. We lack good reference genomes for most species. Rigorous studies to examine genomic processes are generally too costly. On the contrary, cancer studies are fortified with extensive resources for studying genome level dynamics and the interactions among genetic and non-genetic mechanisms. Extensive analysis of primary tumors and metastatic samples have revealed the importance of several genomic mechanisms including higher mutation rates, specific types of mutations, aneuploidy or whole genome doubling and non-genetic effects. Metastatic sites can be directly compared to primary tumor cell counterparts. At the same time, clonal dynamics shape the genomics and evolution of metastatic cancers. Clonal diversity varies by cancer type, and the tumors’ donor and recipient tissues. Still, the cancer research community has been unable to identify any common events that provide a universal predictor of “metastatic potential” which parallels findings in evolutionary ecology. Instead, invasion in cancer studies depends strongly on context, including order of events and clonal composition. The detailed studies of the behavior of a variety of human cancers promises to inform our understanding of genome level dynamics in the diversity of invasive species and provide novel insights for management.

Prospects for the natural distribution of crop wild-relatives with limited adaptability: The case of the wild pea Pisum fulvum

Plant breeders and conservationist depend on knowledge about the genetic variation of their species of interest. Pisum fulvum, a wild relative of domesticated pea, has attracted attention as a genetic resource for crop improvement, yet little information about its diversity in the wild has been published hitherto. We sampled 15 populations of P. fulvum from Israeli natural habitats and conducted genotyping by sequencing to analyse their genetic diversity and adaptive state. We also attempted to evaluate the species past demography and the prospects of its future reaction to environmental changes. The results suggest that genetic diversity of P. fulvum is low to medium and is distributed between well diverged populations. Surprisingly, with 56 % in the total population the selfing rate was found to be significantly lower than expected from a species that is commonly assumed to be a predominant selfer. We found a strong genetic bottleneck during the last glacial period and only limited patterns of isolation by distance and environment, which explained 13 %-18 % of the genetic variation. Despite the weak signatures of genome-wide IBE, 1,354 markers were significantly correlated with environmental factors, 1,233 of which were located within known genes with a nonsynonymous to synonymous ratio of 0.382. Species distribution modelling depicted an ongoing fragmentation and decreased habitable area over the next 80 years under two different socio-economic pathways. Our results suggest that complex interactions of substantial drift and selection shaped the genome of P. fulvum. Climate changeis likely to cause further erosion of genetic diversity in P. fulvum. Systematic ex-situ conservation may be advisable to safeguard genetic variability for future utilization of this species.

Genomic analysis of field pennycress (Thlaspi arvense) provides insights into mechanisms of adaptation to high elevation

Background

Understanding how organisms evolve and adapt to extreme habitats is of crucial importance in evolutionary ecology. Altitude gradients are an important determinant of the distribution pattern and range of organisms due to distinct climate conditions at different altitudes. High-altitude regions often provide extreme environments including low temperature and oxygen concentration, poor soil, and strong levels of ultraviolet radiation, leading to very few plant species being able to populate elevation ranges greater than 4000 m. Field pennycress (Thlaspi arvense) is a valuable oilseed crop and emerging model plant distributed across an elevation range of nearly 4500 m. Here, we generate an improved genome assembly to understand how this species adapts to such different environments.

Results

We sequenced and assembled de novo the chromosome-level pennycress genome of 527.3 Mb encoding 31,596 genes. Phylogenomic analyses based on 2495 single-copy genes revealed that pennycress is closely related to Eutrema salsugineum (estimated divergence 14.32–18.58 Mya), and both species form a sister clade to Schrenkiella parvula and genus Brassica. Field pennycress contains the highest percentage (70.19%) of transposable elements in all reported genomes of Brassicaceae, with the retrotransposon proliferation in the Middle Pleistocene being likely responsible for the expansion of genome size. Moreover, our analysis of 40 field pennycress samples in two high- and two low-elevation populations detected 1,256,971 high-quality single nucleotide polymorphisms. Using three complementary selection tests, we detected 130 candidate naturally selected genes in the Qinghai-Tibet Plateau (QTP) populations, some of which are involved in DNA repair and the ubiquitin system and potential candidates involved in high-altitude adaptation. Notably, we detected a single base mutation causing loss-of-function of the FLOWERING LOCUS C protein, responsible for the transition to early flowering in high-elevation populations.

Conclusions

Our results provide a genome-wide perspective of how plants adapt to distinct environmental conditions across extreme elevation differences and the potential for further follow-up research with extensive data from additional populations and species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01079-0.

De Novo Mutation and Rapid Protein (Co-)evolution during Meiotic Adaptation in Arabidopsis arenosa

A sudden shift in environment or cellular context necessitates rapid adaptation. A dramatic example is genome duplication, which leads to polyploidy. In such situations, the waiting time for new mutations might be prohibitive; theoretical and empirical studies suggest that rapid adaptation will largely rely on standing variation already present in source populations. Here, we investigate the evolution of meiosis proteins in Arabidopsis arenosa, some of which were previously implicated in adaptation to polyploidy, and in a diploid, habitat. A striking and unexplained feature of prior results was the large number of amino acid changes in multiple interacting proteins, especially in the relatively young tetraploid. Here, we investigate whether selection on meiosis genes is found in other lineages, how the polyploid may have accumulated so many differences, and whether derived variants were selected from standing variation. We use a range-wide sample of 145 resequenced genomes of diploid and tetraploid A. arenosa, with new genome assemblies. We confirmed signals of positive selection in the polyploid and diploid lineages they were previously reported in and find additional meiosis genes with evidence of selection. We show that the polyploid lineage stands out both qualitatively and quantitatively. Compared with diploids, meiosis proteins in the polyploid have more amino acid changes and a higher proportion affecting more strongly conserved sites. We find evidence that in tetraploids, positive selection may have commonly acted on de novo mutations. Several tests provide hints that coevolution, and in some cases, multinucleotide mutations, might contribute to rapid accumulation of changes in meiotic proteins.

Genetic and environmental modulation of transposition shapes the evolutionary potential of Arabidopsis thaliana

BACKGROUND

How species can adapt to abrupt environmental changes, particularly in the absence of standing genetic variation, is poorly understood and a pressing question in the face of ongoing climate change. Here we leverage publicly available multi-omic and bio-climatic data for more than 1000 wild Arabidopsis thaliana accessions to determine the rate of transposable element (TE) mobilization and its potential to create adaptive variation in natural settings.

RESULTS

We demonstrate that TE insertions arise at almost the same rate as base substitutions. Mobilization activity of individual TE families varies greatly between accessions, in association with genetic and environmental factors as well as through complex gene-environment interactions. Although the distribution of TE insertions across the genome is ultimately shaped by purifying selection, reflecting their typically strong deleterious effects when located near or within genes, numerous recent TE-containing alleles show signatures of positive selection. Moreover, high rates of transposition appear positively selected at the edge of the species' ecological niche. Based on these findings, we predict through mathematical modeling higher transposition activity in Mediterranean regions within the next decades in response to global warming, which in turn should accelerate the creation of large-effect alleles.

CONCLUSIONS

Our study reveals that TE mobilization is a major generator of genetic variation in A. thaliana that is finely modulated by genetic and environmental factors. These findings and modeling indicate that TEs may be essential genomic players in the demise or rescue of native populations in times of climate crises.

Natural Bacterial Assemblages in Arabidopsis thaliana Tissues Become More Distinguishable and Diverse during Host Development

Developing synthetic microbial communities that can increase plant yield or deter pathogens requires basic research on several fronts, including the efficiency with which microbes colonize plant tissues, how plant genes shape the microbiome, and the microbe-microbe interactions involved in community assembly. Findings on each of these fronts depend upon the spatial and temporal scales at which plant microbiomes are surveyed.

To study the spatial and temporal dynamics of bacterial colonization under field conditions, we planted and sampled Arabidopsis thaliana during 2 years at two Michigan sites and surveyed colonists by sequencing 16S rRNA gene amplicons. Mosaic and dynamic assemblages revealed the plant as a patchwork of tissue habitats that differentiated with age. Although assemblages primarily varied between roots and shoots, amplicon sequence variants (ASVs) also differentiated phyllosphere tissues. Increasing assemblage diversity indicated that variants dispersed more widely over time, decreasing the importance of stochastic variation in early colonization relative to tissue differences. As tissues underwent developmental transitions, the root and phyllosphere assemblages became more distinct. This pattern was driven by common variants rather than those restricted to a particular tissue or transiently present at one developmental stage. Patterns also depended critically on fine phylogenetic resolution: when ASVs were grouped at coarse taxonomic levels, their associations with host tissue and age weakened. Thus, the observed spatial and temporal variation in colonization depended upon bacterial traits that were not broadly shared at the family level. Some colonists were consistently more successful at entering specific tissues, as evidenced by their repeatable spatial prevalence distributions across sites and years. However, these variants did not overtake plant assemblages, which instead became more even over time. Together, these results suggested that the increasing effect of tissue type was related to colonization bottlenecks for specific ASVs rather than to their ability to dominate other colonists once established.

Fitness effects of mutation in natural populations of Arabidopsis thaliana reveal a complex influence of local adaptation

Little is empirically known about the contribution of mutations to fitness in natural environments. However, Fisher's Geometric Model (FGM) provides a conceptual foundation to consider the influence of the environment on mutational effects. To quantify mutational properties in the field, we established eight sets of MA lines (7-10 generations) derived from eight founders collected from natural populations of Arabidopsis thaliana from French and Swedish sites, representing the range margins of the species in Europe. We reciprocally planted the MA lines and their founders at French and Swedish sites, allowing us to test predictions of FGM under naturally occurring environmental conditions. The performance of the MA lines relative to each other and to their respective founders confirmed some and contradicted other predictions of the FGM: the contribution of mutation to fitness variance increased when the genotype was in an environment where its fitness was low, that is, in the away environment, but mutations were more likely to be beneficial when the genotype was in its home environment. Consequently, environmental context plays a large role in the contribution of mutations to the evolutionary process and local adaptation does not guarantee that a genotype is at or close to its optimum.

A heuristic model of the effects of phenotypic robustness in adaptive evolution

A recent theoretical, deterministic model of the effects of phenotypic robustness on adaptive evolutionary dynamics showed that a certain level of phenotypic robustness (critical robustness) is a required condition for adaptation to occur and to be maintained during evolution in most real organismal systems. We built an individual-based heuristic model to verify the soundness of these theoretical results through computer simulation, testing expectations under a range of scenarios for the relevant parameters of the evolutionary dynamics. These include the mutation probability, the presence of stochastic effects, the introduction of environmental influences and the possibility for some features of the population (like selection coefficients and phenotypic robustness) to change themselves during adaptation. Overall, we found a good match between observed and expected results, even for evolutionary parameter values that violate some of the assumptions of the deterministic model, and that robustness can itself evolve. However, from more than one simulation it appears that very high robustness values, higher than the critical value, can limit or slow-down adaptation. This possible trade-off was not predicted by the deterministic model.

Evidence for the Accumulation of Nonsynonymous Mutations and Favorable Pleiotropic Alleles During Wheat Breeding

Plant breeding leads to the genetic improvement of target traits by selecting a small number of genotypes from among typically large numbers of candidate genotypes after careful evaluation. In this study, we first investigated how mutations at conserved nucleotide sites normally viewed as deleterious, such as nonsynonymous sites, accumulated in a wheat, Triticum aestivum, breeding lineage. By comparing a 150 year old ancestral and modern cultivar, we found recent nucleotide polymorphisms altered amino acids and occurred within conserved genes at frequencies expected in the absence of purifying selection. Mutations that are deleterious in other contexts likely had very small or no effects on target traits within the breeding lineage. Second, we investigated if breeders selected alleles with favorable effects on some traits and unfavorable effects on others and used different alleles to compensate for the latter. An analysis of a segregating population derived from the ancestral and modern parents provided one example of this phenomenon. The recent cultivar contains the Rht-B1b green revolution semi-dwarfing allele and compensatory alleles that reduce its negative effects. However, improvements in traits other than plant height were due to pleiotropic loci with favorable effects on traits and to favorable loci with no detectable pleiotropic effects. Wheat breeding appears to tolerate mutations at conserved nucleotide sites and to only select for alleles with both favorable and unfavorable effects on traits in exceptional situations.

What Do We Really Know About Adaptation at Range Edges?

Recent theory and empirical evidence have provided new insights regarding how evolutionary forces interact to shape adaptation at stable and transient range margins. Predictions regarding trait divergence at leading edges are frequently supported. However, declines in fitness at and beyond edges show that trait divergence has sometimes been insufficient to maintain high fitness, so identifying constraints to adaptation at range edges remains a key challenge. Indirect evidence suggests that range expansion may be limited by adaptive genetic variation, but direct estimates of genetic constraints at and beyond range edges are still scarce. Sequence data suggest increased genetic load in edge populations in several systems, but its causes and fitness consequences are usually poorly understood. The balance between maladaptive and positive effects of gene flow on fitness at range edges deserves further study. It is becoming increasingly clear that characterizations about degree of adaptation based solely on geographical peripherality are unsupported.

Limited divergent adaptation despite a substantial environmental cline in wild pea

Isolation by environment (IBE) is a widespread phenomenon in nature. It is commonly expected that the degree of difference among environments is proportional to the level of divergence between populations in their respective environments. It is therefore assumed that a species' genetic diversity displays a pattern of IBE in the presence of a strong environmental cline if gene flow does not mitigate isolation. We tested this common assumption by analysing the genetic diversity and demographic history of Pisum fulvum, which inhabits contrasting habitats in the southern Levant and is expected to display only minor migration rates between populations, making it an ideal test case. Ecogeographical and subpopulation structure were analysed and compared. The correlation of genetic with environmental distances was calculated to test the effect of isolation by distance and IBE and detect the main drivers of these effects. Historical effective population size was estimated using stairway plot. Limited overlap of ecogeographical and genetic clustering was observed, and correlation between genetic and environmental distances was statistically significant but small. We detected a sharp decline of effective population size during the last glacial period. The low degree of IBE may be the result of genetic drift due to a past bottleneck. Our findings contradict the expectation that strong environmental clines cause IBE in the absence of extensive gene flow.

A genome assembly and the somatic genetic and epigenetic mutation rate in a wild long-lived perennial Populus trichocarpa

BACKGROUND

Plants can transmit somatic mutations and epimutations to offspring, which in turn can affect fitness. Knowledge of the rate at which these variations arise is necessary to understand how plant development contributes to local adaption in an ecoevolutionary context, particularly in long-lived perennials.

RESULTS

Here, we generate a new high-quality reference genome from the oldest branch of a wild Populus trichocarpa tree with two dominant stems which have been evolving independently for 330 years. By sampling multiple, age-estimated branches of this tree, we use a multi-omics approach to quantify age-related somatic changes at the genetic, epigenetic, and transcriptional level. We show that the per-year somatic mutation and epimutation rates are lower than in annuals and that transcriptional variation is mainly independent of age divergence and cytosine methylation. Furthermore, a detailed analysis of the somatic epimutation spectrum indicates that transgenerationally heritable epimutations originate mainly from DNA methylation maintenance errors during mitotic rather than during meiotic cell divisions.

CONCLUSION

Taken together, our study provides unprecedented insights into the origin of nucleotide and functional variation in a long-lived perennial plant.

Adaptation to ionizing radiation of higher plants: From environmental radioactivity to chernobyl disaster

The purpose of this work is to highlight the effects of ionizing radiation on the genetic material in higher plants by assessing both adaptive processes as well as the evolution of plant species. The effects that the ionizing radiation has on greenery following a nuclear accident, was examined by taking the Chernobyl Nuclear Power Plant disaster as a case study. The genetic and evolutionary effects that ionizing radiation had on plants after the Chernobyl accident were highlighted. The response of biota to Chernobyl irradiation was a complex interaction among radiation dose, dose rate, temporal and spatial variation, varying radiation sensitivities of the different plants' species, and indirect effects from other events. Ionizing radiation causes water radiolysis, generating highly reactive oxygen species (ROS). ROS induce the rapid activation of detoxifying enzymes. DeoxyriboNucleic Acid (DNA) is the object of an attack by both, the hydroxyl ions and the radiation itself, thus triggering a mechanism both direct and indirect. The effects on DNA are harmful to the organism and the long-term development of the species. Dose-dependent aberrations in chromosomes are often observed after irradiation. Although multiple DNA repair mechanisms exist, double-strand breaks (DSBs or DNA-DSBs) are often subject to errors. Plants DSBs repair mechanisms mainly involve homologous and non-homologous dependent systems, the latter especially causing a loss of genetic information. Repeated ionizing radiation (acute or chronic) ensures that plants adapt, demonstrating radioresistance. An adaptive response has been suggested for this phenomenon. As a result, ionizing radiation influences the genetic structure, especially during chronic irradiation, reducing genetic variability. This reduction may be associated with the fact that particular plant species are more subject to chronic stress, confirming the adaptive theory. Therefore, the genomic effects of ionizing radiation demonstrate their likely involvement in the evolution of plant species.