Adaptation to climate across the Arabidopsis thaliana genome

Within-Generation Polygenic Selection Shapes Fitness-Related Traits across Environments in Juvenile Sea Bream

Understanding the genetic underpinnings of fitness trade-offs across spatially variable environments remains a major challenge in evolutionary biology. In Mediterranean gilthead sea bream, first-year juveniles use various marine and brackish lagoon nursery habitats characterized by a trade-off between food availability and environmental disturbance. Phenotypic differences among juveniles foraging in different habitats rapidly appear after larval settlement, but the relative role of local selection and plasticity in phenotypic variation remains unclear. Here, we combine phenotypic and genetic data to address this question. We first report correlations of opposite signs between growth and condition depending on juvenile habitat type. Then, we use single nucleotide polymorphism (SNP) data obtained by Restriction Associated DNA (RAD) sequencing to search for allele frequency changes caused by a single generation of spatially varying selection between habitats. We found evidence for moderate selection operating at multiple loci showing subtle allele frequency shifts between groups of marine and brackish juveniles. We identified subsets of candidate outlier SNPs that, in interaction with habitat type, additively explain up to 3.8% of the variance in juvenile growth and 8.7% in juvenile condition; these SNPs also explained significant fraction of growth rate in an independent larval sample. Our results indicate that selective mortality across environments during early-life stages involves complex trade-offs between alternative growth strategies.

Root System Depth in Arabidopsis Is Shaped by EXOCYST70A3 via the Dynamic Modulation of Auxin Transport

Root system architecture (RSA), the distribution of roots in soil, plays a major role in plant survival. RSA is shaped by multiple developmental processes that are largely governed by the phytohormone auxin, suggesting that auxin regulates responses of roots that are important for local adaptation. However, auxin has a central role in numerous processes, and it is unclear which molecular mechanisms contribute to the variation in RSA for environmental adaptation. Using natural variation in Arabidopsis, we identify EXOCYST70A3 as a modulator of the auxin system that causes variation in RSA by acting on PIN4 protein distribution. Allelic variation and genetic perturbation of EXOCYST70A3 lead to alteration of root gravitropic responses, resulting in a different RSA depth profile and drought resistance. Overall our findings suggest that the local modulation of the pleiotropic auxin pathway can gives rise to distinct RSAs that can be adaptive in specific environments.

Environmental Genome-Wide Association Reveals Climate Adaptation Is Shaped by Subtle to Moderate Allele Frequency Shifts in Loblolly Pine

Understanding the genomic basis of local adaptation is crucial to determine the potential of long-lived woody species to withstand changes in their natural environment. In the past, efforts to dissect the genomic architecture in gymnosperms species have been limited due to the absence of reference genomes. Recently, the genomes of some commercially important conifers, such as loblolly pine, have become available, allowing whole-genome studies of these species. In this study, we test for associations between 87k SNPs, obtained from whole-genome resequencing of loblolly pine individuals, and 270 environmental variables and combinations of them. We determine the geographic location of significant loci and identify their genomic location using our newly constructed ultradense 26k SNP linkage map. We found that water availability is the main climatic variable shaping local adaptation of the species, and found 821 SNPs showing significant associations with climatic variables or combinations of them based on the consistent results of three different genotype–environment association methods. Our results suggest that adaptation to climate in the species might have occurred by many changes in the frequency of alleles with moderate to small effect sizes, and by the smaller contribution of large effect alleles in genes related to moisture deficit, temperature and precipitation. Genomic regions of low recombination and high population differentiation harbored SNPs associated with groups of environmental variables, suggesting climate adaptation might have evolved as a result of different selection pressures acting on groups of genes associated with an aspect of climate rather than on individual environmental variables.

Evaluation of Seven Different RNA-Seq Alignment Tools Based on Experimental Data from the Model Plant Arabidopsis thaliana

Quantification of gene expression is crucial to connect genome sequences with phenotypic and physiological data. RNA-Sequencing (RNA-Seq) has taken a prominent role in the study of transcriptomic reactions of plants to various environmental and genetic perturbations. However, comparative tests of different tools for RNA-Seq read mapping and quantification have been mainly performed on data from animals or humans, which necessarily neglect, for example, the large genetic variability among natural accessions within plant species. Here, we compared seven computational tools for their ability to map and quantify Illumina single-end reads from the Arabidopsis thaliana accessions Columbia-0 (Col-0) and N14. Between 92.4% and 99.5% of all reads were mapped to the reference genome or transcriptome and the raw count distributions obtained from the different mappers were highly correlated. Using the software DESeq2 to determine differential gene expression (DGE) between plants exposed to 20 °C or 4 °C from these read counts showed a large pairwise overlap between the mappers. Interestingly, when the commercial CLC software was used with its own DGE module instead of DESeq2, strongly diverging results were obtained. All tested mappers provided highly similar results for mapping Illumina reads of two polymorphic Arabidopsis accessions to the reference genome or transcriptome and for the determination of DGE when the same software was used for processing.

Growth under Fluctuating Light Reveals Large Trait Variation in a Panel of Arabidopsis Accessions

The capacity of photoautotrophs to fix carbon depends on the efficiency of the conversion of light energy into chemical potential by photosynthesis. In nature, light input into photosynthesis can change very rapidly and dramatically. To analyze how genetic variation in Arabidopsis thaliana affects photosynthesis and growth under dynamic light conditions, 36 randomly chosen natural accessions were grown under uniform and fluctuating light intensities. After 14 days of growth under uniform or fluctuating light regimes, maximum photosystem II quantum efficiency (Fv/Fm) was determined, photosystem II operating efficiency (ΦPSII) and non-photochemical quenching (NPQ) were measured in low light, and projected leaf area (PLA) as well as the number of visible leaves were estimated. Our data show that ΦPSII and PLA were decreased and NPQ was increased, while Fv/Fm and number of visible leaves were unaffected, in most accessions grown under fluctuating compared to uniform light. There were large changes between accessions for most of these parameters, which, however, were not correlated with genomic variation. Fast growing accessions under uniform light showed the largest growth reductions under fluctuating light, which correlated strongly with a reduction in ΦPSII, suggesting that, under fluctuating light, photosynthesis controls growth and not vice versa.

In the presence of population structure: From genomics to candidate genes underlying local adaptation

Understanding the genomic signatures, genes, and traits underlying local adaptation of organisms to heterogeneous environments is of central importance to the field evolutionary biology. To identify loci underlying local adaptation, models that combine allelic and environmental variation while controlling for the effects of population structure have emerged as the method of choice. Despite being evaluated in simulation studies, there has not been a thorough investigation of empirical evidence supporting local adaptation across these alleles. To evaluate these methods, we use 875 Arabidopsis thaliana Eurasian accessions and two mixed models (GEMMA and LFMM) to identify candidate SNPs underlying local adaptation to climate. Subsequently, to assess evidence of local adaptation and function among significant SNPs, we examine allele frequency differentiation and recent selection across Eurasian populations, in addition to their distribution along quantitative trait loci (QTL) explaining fitness variation between Italy and Sweden populations and cis-regulatory/nonsynonymous sites showing significant selective constraint. Our results indicate that significant LFMM/GEMMA SNPs show low allele frequency differentiation and linkage disequilibrium across locally adapted Italy and Sweden populations, in addition to a poor association with fitness QTL peaks (highest logarithm of odds score). Furthermore, when examining derived allele frequencies across the Eurasian range, we find that these SNPs are enriched in low-frequency variants that show very large climatic differentiation but low levels of linkage disequilibrium. These results suggest that their enrichment along putative functional sites most likely represents deleterious variation that is independent of local adaptation. Among all the genomic signatures examined, only SNPs showing high absolute allele frequency differentiation (AFD) and linkage disequilibrium (LD) between Italy and Sweden populations showed a strong association with fitness QTL peaks and were enriched along selectively constrained cis-regulatory/nonsynonymous sites. Using these SNPs, we find strong evidence linking flowering time, freezing tolerance, and the abscisic-acid pathway to local adaptation.

Evolutionary genomics can improve prediction of species' responses to climate change

Global climate change (GCC) increasingly threatens biodiversity through the loss of species, and the transformation of entire ecosystems. Many species are challenged by the pace of GCC because they might not be able to respond fast enough to changing biotic and abiotic conditions. Species can respond either by shifting their range, or by persisting in their local habitat. If populations persist, they can tolerate climatic changes through phenotypic plasticity, or genetically adapt to changing conditions depending on their genetic variability and census population size to allow for de novo mutations. Otherwise, populations will experience demographic collapses and species may go extinct. Current approaches to predicting species responses to GCC begin to combine ecological and evolutionary information for species distribution modelling. Including an evolutionary dimension will substantially improve species distribution projections which have not accounted for key processes such as dispersal, adaptive genetic change, demography, or species interactions. However, eco-evolutionary models require new data and methods for the estimation of a species' adaptive potential, which have so far only been available for a small number of model species. To represent global biodiversity, we need to devise large-scale data collection strategies to define the ecology and evolutionary potential of a broad range of species, especially of keystone species of ecosystems. We also need standardized and replicable modelling approaches that integrate these new data to account for eco-evolutionary processes when predicting the impact of GCC on species' survival. Here, we discuss different genomic approaches that can be used to investigate and predict species responses to GCC. This can serve as guidance for researchers looking for the appropriate experimental setup for their particular system. We furthermore highlight future directions for moving forward in the field and allocating available resources more effectively, to implement mitigation measures before species go extinct and ecosystems lose important functions.

In the presence of population structure: From genomics to candidate genes underlying local adaptation

Understanding the genomic signatures, genes, and traits underlying local adaptation of organisms to heterogeneous environments is of central importance to the field evolutionary biology. To identify loci underlying local adaptation, models that combine allelic and environmental variation while controlling for the effects of population structure have emerged as the method of choice. Despite being evaluated in simulation studies, there has not been a thorough investigation of empirical evidence supporting local adaptation across these alleles. To evaluate these methods, we use 875 Arabidopsis thaliana Eurasian accessions and two mixed models (GEMMA and LFMM) to identify candidate SNPs underlying local adaptation to climate. Subsequently, to assess evidence of local adaptation and function among significant SNPs, we examine allele frequency differentiation and recent selection across Eurasian populations, in addition to their distribution along quantitative trait loci (QTL) explaining fitness variation between Italy and Sweden populations and cis-regulatory/nonsynonymous sites showing significant selective constraint. Our results indicate that significant LFMM/GEMMA SNPs show low allele frequency differentiation and linkage disequilibrium across locally adapted Italy and Sweden populations, in addition to a poor association with fitness QTL peaks (highest logarithm of odds score). Furthermore, when examining derived allele frequencies across the Eurasian range, we find that these SNPs are enriched in low-frequency variants that show very large climatic differentiation but low levels of linkage disequilibrium. These results suggest that their enrichment along putative functional sites most likely represents deleterious variation that is independent of local adaptation. Among all the genomic signatures examined, only SNPs showing high absolute allele frequency differentiation (AFD) and linkage disequilibrium (LD) between Italy and Sweden populations showed a strong association with fitness QTL peaks and were enriched along selectively constrained cis-regulatory/nonsynonymous sites. Using these SNPs, we find strong evidence linking flowering time, freezing tolerance, and the abscisic-acid pathway to local adaptation.

Natural variation on whole-plant form in the wild is influenced by multivariate soil nutrient characteristics: natural selection acts on root traits

PREMISE

In the complex soil nutrient environments of wild populations of annual plants, in general, low nutrient availability restricts growth and alters root-shoot relationships. However, our knowledge of natural selection on roots in field settings is limited. We sought to determine whether selection acts directly on root traits and to identify which components of the soil environment were potential agents of selection.

METHODS

We studied wild native populations of Arabidopsis thaliana across 4 years, measuring aboveground and belowground traits and analyzing soil nutrients. Using multivariate methods, we examined patterns of natural selection and identified soil attributes that contributed to whole-plant form. In a common garden experiment at two field sites with contrasting soil texture, we examined patterns of selection on root and shoot traits.

RESULTS

In wild populations, we uncovered selection for above- and belowground size and architectural traits. We detected variation through time and identified soil components that influenced fruit production. In the garden experiment, we detected a distinct positive selection for total root length at the site with greater water-holding capacity and negative selection for measures of root architecture at the field site with reduced nutrient availability and water holding capacity.

CONCLUSIONS

Patterns of natural selection on belowground traits varied through time, across field sites and experimental gardens. Simultaneous investigations of above- and belowground traits reveal trait functional relationships on which natural selection can act, highlighting the influence of edaphic features on evolutionary processes in wild annual plant populations.

Identification of loci controlling adaptation in Chinese soya bean landraces via a combination of conventional and bioclimatic GWAS

Landraces often contain genetic diversity that has been lost in modern cultivars, including alleles that confer enhanced local adaptation. To comprehensively identify loci associated with adaptive traits in soya bean landraces, for example flowering time, a population of 1938 diverse landraces and 97 accessions of the wild progenitor of cultivated soya bean, Glycine soja was genotyped using tGBS® . Based on 99 085 high-quality SNPs, landraces were classified into three sub-populations which exhibit geographical genetic differentiation. Clustering was inferred from STRUCTURE, principal component analyses and neighbour-joining tree analyses. Using phenotypic data collected at two locations separated by 10 degrees of latitude, 17 trait-associated SNPs (TASs) for flowering time were identified, including a stable locus Chr12:5914898 and previously undetected candidate QTL/genes for flowering time in the vicinity of the previously cloned flowering genes, E1 and E2. Using passport data associated with the collection sites of the landraces, 27 SNPs associated with adaptation to three bioclimatic variables (temperature, daylength, and precipitation) were identified. A series of candidate flowering genes were detected within linkage disequilibrium (LD) blocks surrounding 12 bioclimatic TASs. Nine of these TASs exhibit significant differences in flowering time between alleles within one or more of the three individual sub-populations. Signals of selection during domestication and/or subsequent landrace diversification and adaptation were detected at 38 of the 44 flowering and bioclimatic TASs. Hence, this study lays the groundwork to begin breeding for novel environments predicted to arise following global climate change.

Adaptive Growth: Shaping Auxin-Mediated Root System Architecture

Root system architecture (RSA), governed by the phytohormone auxin, endows plants with an adaptive advantage in particular environments. Using geographically representative arabidopsis (Arabidopsis thaliana) accessions as a resource for GWA mapping, Waidmann et al. and Ogura et al. recently identified two novel components involved in modulating auxin-mediated RSA and conferring plant fitness in particular habitats.

Effects of two centuries of global environmental variation on phenology and physiology of Arabidopsis thaliana

Intraspecific trait variation is caused by genetic and plastic responses to environment. This intraspecific diversity is captured in immense natural history collections, giving us a window into trait variation across continents and through centuries of environmental shifts. Here we tested if hypotheses based on life history and the leaf economics spectrum explain intraspecific trait changes across global spatiotemporal environmental gradients. We measured phenotypes on a 216-year time series of Arabidopsis thaliana accessions from across its native range and applied spatially varying coefficient models to quantify region-specific trends in trait coordination and trait responses to climate gradients. All traits exhibited significant change across space or through time. For example, δ¹⁵ N decreased over time across much of the range and leaf C:N increased, consistent with predictions based on anthropogenic changes in land use and atmosphere. Plants were collected later in the growing season in more recent years in many regions, possibly because populations shifted toward more spring germination and summer flowering as opposed to fall germination and spring flowering. When climate variables were considered, collection dates were earlier in warmer years, while summer rainfall had opposing associations with collection date depending on regions. There was only a modest correlation among traits, indicating a lack of a single life history/physiology axis. Nevertheless, leaf C:N was low for summer- versus spring-collected plants, consistent with a life history-physiology axis from slow-growing winter annuals to fast-growing spring/summer annuals. Regional heterogeneity in phenotype trends indicates complex responses to spatiotemporal environmental gradients potentially due to geographic genetic variation and climate interactions with other aspects of environment. Our study demonstrates how natural history collections can be used to broadly characterize trait responses to environment, revealing heterogeneity in response to anthropogenic change.

Adaptation to Plant Communities across the Genome of Arabidopsis thaliana

Despite the importance of plant-plant interactions on plant community dynamics and crop yield, our understanding of the adaptive genetics underlying these interactions is still limited and deserves to be investigated in the context of complex and diffuse interactions occurring in plant assemblages. Here, based on 145 natural populations of Arabidopsis thaliana located in south-west of France and characterized for plant communities, we conducted a Genome-Environment Association analysis to finely map adaptive genomic regions of A. thaliana associated with plant community descriptors. To control for correlated abiotic environment effects, we also characterized the populations for a set of biologically meaningful climate and soil variables. A nonnegligible fraction of top single nucleotide polymorphisms was associated with both plant community descriptors and abiotic variables, highlighting the importance of considering the actual abiotic drivers of plant communities to disentangle genetic variants for biotic adaptation from genetic variants for abiotic adaptation. The adaptive loci associated with species abundance were highly dependent on the identity of the neighboring species suggesting a high degree of biotic specialization of A. thaliana to members of its plant interaction network. Moreover, the identification of adaptive loci associated with α-diversity and composition of plant communities supports the ability of A. thaliana to interact simultaneously with multiple plant neighbors, which in turn can help to understand the role of community-wide selection. Altogether, our study highlights that dissecting the genetic basis underlying plant-plant interactions at a regional scale while controlling for abiotic confounding factors can help understanding the adaptive mechanisms modulating natural plant assemblages.

TRANSPARENT TESTA GLABRA 1 participates in flowering time regulation in Arabidopsis thaliana

Pleiotropic regulatory factors mediate concerted responses of the plant’s trait network to endogenous and exogenous cues. TRANSPARENT TESTA GLABRA 1 (TTG1) is such a factor that has been predominantly described as a regulator of early developmental traits. Although its closest homologs LIGHT-REGULATED WD1 (LWD1) and LWD2 affect photoperiodic flowering, a role of TTG1 in flowering time regulation has not been reported. Here we reveal that TTG1 is a regulator of flowering time in Arabidopsis thaliana and changes transcript levels of different targets within the flowering time regulatory pathway. TTG1 mutants flower early and TTG1 overexpression lines flower late at long-day conditions. Consistently, TTG1 can suppress the transcript levels of the floral integrators FLOWERING LOCUS T and SUPPRESSOR OF OVEREXPRESSION OF CO1 and can act as an activator of circadian clock components. Moreover, TTG1 might form feedback loops at the protein level. The TTG1 protein interacts with PSEUDO RESPONSE REGULATOR (PRR)s and basic HELIX-LOOP-HELIX 92 (bHLH92) in yeast. In planta, the respective pairs exhibit interesting patterns of localization including a recruitment of TTG1 by PRR5 to subnuclear foci. This mechanism proposes additional layers of regulation by TTG1 and might aid to specify the function of bHLH92. Within another branch of the pathway, TTG1 can elevate FLOWERING LOCUS C (FLC) transcript levels. FLC mediates signals from the vernalization, ambient temperature and autonomous pathway and the circadian clock is pivotal for the plant to synchronize with diurnal cycles of environmental stimuli like light and temperature. Our results suggest an unexpected positioning of TTG1 upstream of FLC and upstream of the circadian clock. In this light, this points to an adaptive value of the role of TTG1 in respect to flowering time regulation.

Genomics of Long- and Short-Term Adaptation in Maize and Teosintes

Maize is an excellent model for the study of plant adaptation. Indeed, post domestication maize quickly adapted to a host of new environments across the globe. And work over the last decade has begun to highlight the role of the wild relatives of maize-the teosintes Zea mays ssp. parviglumis and ssp. mexicana-as excellent models for dissecting long-term local adaptation.Although human-driven selection associated with maize domestication has been extensively studied, the genetic basis of natural variation is still poorly understood. Here we review studies on the genetic basis of adaptation and plasticity in maize and its wild relatives. We highlight a range of different processes that contribute to adaptation and discuss evidence from natural, cultivated, and experimental populations. From an applied perspective, understanding the genetic bases of adaptation and the contribution of plasticity will provide us with new tools to both better understand and mitigate the effect of climate changes on natural and cultivated populations.

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

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

Arabidopsis halleri: a perennial model system for studying population differentiation and local adaptation

Local adaptation is assumed to occur when populations differ in a phenotypic trait or a set of traits, and such variation has a genetic basis. Here, we introduce Arabidopsis halleri and its life history as a perennial model system to study population differentiation and local adaptation. Studies on altitudinal adaptation have been conducted in two regions: Mt. Ibuki in Japan and the European Alps. Several studies have demonstrated altitudinal adaptation in ultraviolet-B (UV-B) tolerance, leaf water repellency against spring frost and anti-herbivore defences. Studies on population differentiation in A. halleri have also focused on metal hyperaccumulation and tolerance to heavy metal contamination. In these study systems, genome scans to identify candidate genes under selection have been applied. Lastly, we briefly discuss how RNA-Seq can broaden phenotypic space and serve as a link to underlying mechanisms. In conclusion, A. halleri provides us with opportunities to study population differentiation and local adaptation, and relate these to the genetic systems underlying target functional traits.

Natural Root Cellular Variation in Responses to Osmotic Stress in Arabidopsis thaliana Accessions

Arabidopsis naturally occurring populations have allowed for the identification of considerable genetic variation remodeled by adaptation to different environments and stress conditions. Water is a key resource that limits plant growth, and its availability is initially sensed by root tissues. The root’s ability to adjust its physiology and morphology under water deficit makes this organ a useful model to understand how plants respond to water stress. Here, we used hyperosmotic shock stress treatments in different Arabidopsis accessions to analyze the root cell morphological responses. We found that osmotic stress conditions reduced root growth and root apical meristem (RAM) size, promoting premature cell differentiation without affecting the stem cell niche morphology. This phenotype was accompanied by a cluster of small epidermal and cortex cells with radial expansion and root hairs at the transition to the elongation zone. We also found this radial expansion with root hairs when plants are grown under hypoosmotic conditions. Finally, root growth was less affected by osmotic stress in the Sg-2 accession followed by Ws, Cvi-0, and Col-0; however, after a strong osmotic stress, Sg-2 and Cvi-0 were the most resilience accessions. The sensitivity differences among these accessions were not explained by stress-related gene expression. This work provides new cellular insights on the Arabidopsis root phenotypic variability and plasticity to osmotic stress.

Distributed phenomics with the unPAK project reveals the effects of mutations

Determining how genes are associated with traits in plants and other organisms is a major challenge in modern biology. The unPAK project - undergraduates phenotyping Arabidopsis knockouts - has generated phenotype data for thousands of non-lethal insertion mutation lines within a single Arabidopsis thaliana genomic background. The focal phenotypes examined by unPAK are complex macroscopic fitness-related traits, which have ecological, evolutionary and agricultural importance. These phenotypes are placed in the context of the wild-type and also natural accessions (phytometers), and standardized for environmental differences between assays. Data from the unPAK project are used to describe broad patterns in the phenotypic consequences of insertion mutation, and to identify individual mutant lines with distinct phenotypes as candidates for further study. Inclusion of undergraduate researchers is at the core of unPAK activities, and an important broader impact of the project is providing students an opportunity to obtain research experience.

Natural selection on the Arabidopsis thaliana genome in present and future climates

Through the lens of evolution, climate change is an agent of natural selection that forces populations to change and adapt, or face extinction. However, current assessments of the risk of biodiversity associated with climate change¹ do not typically take into account how natural selection influences populations differently depending on their genetic makeup². Here we make use of the extensive genome information that is available for Arabidopsis thaliana and measure how manipulation of the amount of rainfall affected the fitness of 517 natural Arabidopsis lines that were grown in Spain and Germany. This allowed us to directly infer selection along the genome³. Natural selection was particularly strong in the hot-dry location in Spain, where 63% of lines were killed and where natural selection substantially changed the frequency of approximately 5% of all genome-wide variants. A significant portion of this climate-driven natural selection of variants was predictable from signatures of local adaptation (R² = 29-52%), as genetic variants that were found in geographical areas with climates more similar to the experimental sites were positively selected. Field-validated predictions across the species range indicated that Mediterranean and western Siberian populations-at the edges of the environmental limits of this species-currently experience the strongest climate-driven selection. With more frequent droughts and rising temperatures in Europe⁴, we forecast an increase in directional natural selection moving northwards from the southern end of Europe, putting many native A. thaliana populations at evolutionary risk.

A Missense Mutation of the HSPB7 Gene Associated with Heat Tolerance in Chinese Indicine Cattle

Simple Summary

A missense mutation (NC_037329.1: g.136054902 C > G: p. Ala69Gly) was identified in the heat shock protein family B (small) member 7 (HSPB7) gene in indicine cattle, which might be a candidate mutation associated with the heat tolerance. Here, Polymerase Chain Reaction and DNA sequencing methods were used to detect this mutation in 774 individuals belonging to 32 Chinese indigenous cattle breeds. The distribution of alleles of NC_037329.1: g.136054902 C > G displays significant geographical difference across native Chinese cattle breeds and cattle carrying allele G distributed in regions with higher mean annual temperature, relative humidity, and temperature humidity index. Our results demonstrate that the mutation of the HSPB7 gene in Chinese indicine cattle might be a candidate gene associated with the heat tolerance.

Abstract

The small heat shock proteins (HSPB) are expressed in response to heat stress, and the heat shock protein family B (small) member 7 (HSPB7) gene has been reported to play an important role in heat tolerance pathways. Only a missense mutation (NC_037329.1: g.136054902 C > G: p.Ala69Gly) was identified in the HSPB7 gene in indicine cattle, which might be a candidate mutation associated with the heat tolerance. Here, we explore the allele frequency of this mutation in 774 individuals belonging to 32 Chinese indigenous cattle breeds using polymerase chain reaction (PCR) and DNA sequencing methods. The distribution of alleles of NC_037329.1: g.136054902 C > G displays significant geographical difference across native Chinese cattle breeds that the allele C was dominant in northern cattle groups, while allele G was dominant in southern indicine cattle groups. Additionally, the association analysis indicated that the G allele was significantly associated with mean annual temperature (T), relative humidity (RH), and temperature humidity index (THI) (p < 0.01), suggesting that cattle carrying allele G were distributed in regions with higher T, RH, and THI. Our results demonstrate that the mutation of the HSPB7 gene in Chinese indicine cattle might be a candidate gene associated with the heat tolerance.

Population climatic history predicts phenotypic responses in novel environments for Arabidopsis thaliana in North America

PREMISE

Determining how species perform in novel climatic environments is essential for understanding (1) responses to climate change and (2) evolutionary consequences of biological invasions. For the vast majority of species, the number of population characteristics that will predict performance and patterns of natural selection in novel locations in the wild remains limited.

METHODS

We evaluated phenological, vegetative, architectural, and fitness-related traits in experimental gardens in contrasting climates (Ontario, Canada, and South Carolina, USA) in the North American non-native distribution of Arabidopsis thaliana. We assessed the effects of climatic distance, geographic distance, and genetic features of history on performance and patterns of natural selection in the novel garden settings.

RESULTS

We found that plants had greater survivorship, flowered earlier, were larger, and produced more fruit in the south, and that genotype-by-environment interactions were significant between gardens. However, our analyses revealed similar patterns of natural selection between gardens in distinct climate zones. After accounting for genetic ancestry, we also detected that population climatic distance best predicted performance within gardens.

CONCLUSIONS

These data suggest that colonization success in novel, non-native environments is determined by a combination of climate and genetic history. When performance at novel sites was assessed with seed sources from geographically and genetically disparate, established non-native populations, proximity to the garden alone was insufficient to predict performance. Our study highlights the need to evaluate seed sources from diverse origins to describe comprehensively phenotypic responses to novel environments, particularly for taxa in which many source populations may contribute to colonization.

The response profile to chronic radiation exposure based on the transcriptome analysis of Scots pine from Chernobyl affected zone

Radioactive contamination of the natural areas is one of the most long-lasting anthropogenic impacts on the environment. Scots pine (Pinus sylvestris L.) is a promising organism for radiation-related research because of its high radiosensitivity, but the genome size of Pinacea species has imposed obstacles for high-throughput studies so far. In this work, we conducted the analysis of the de novo assembled transcriptome of Scots pine populations growing in the Chernobyl-affected zone, which is still today contaminated with radionuclides because of the accident at the nuclear power plant in 1986. The transcriptome profiles indicate a clear pattern of adaptive stress response, which seems to be dose-dependent. The transcriptional response indicates a continuous modulation of the cellular redox system, enhanced expression of chaperones and histones, along with the control of ions balance. Interestingly, the activity of transposable element families is inversely correlated to the exposure levels to radiation. These adaptive responses, which are triggered by radiation doses 30 times lower than the one accepted as a safe for biota species by international regulations, suggest that the environmental management in radiation protection should be reviewed.

Habitat preference differentiates the Holocene range dynamics but not barrier effects on two sympatric, congeneric trees (Tristaniopsis, Myrtaceae)

Niche partitioning can lead to differences in the range dynamics of plant species through its impacts on habitat availability, dispersal, or selection for traits that affect colonization and persistence. We investigated whether niche partitioning into upland and riparian habitats differentiates the range dynamics of two closely related and sympatric eastern Australian trees: the mountain water gum (Tristaniopsis collina) and the water gum (T. laurina). Using genomic data from SNP genotyping of 480 samples, we assessed the impact of biogeographic barriers and tested for signals of range expansion. Circuit theory was used to model isolation-by-resistance across three palaeo-environment scenarios: the Last Glacial Maximum, the Holocene Climate Optimum and present-day (1950-2014). Both trees showed similar genetic structure across historically dry barriers, despite evidence of significant environmental niche differentiation and different post-glacial habitat shifts. Tristaniopsis collina exhibits the signature of serial founder effects consistent with recent or rapid range expansion, whilst T. laurina has genetic patterns consistent with long-term persistence in geographically isolated populations despite occupying a broader bioclimatic niche. We found the minor influence of isolation-by-resistance on both species, though other unknown factors appear to shape genetic variation. We postulate that specialized recruitment traits (adapted to flood-disturbance regimes) rather than habitat availability limited post-glacial range expansion in T. laurina. Our findings indicate that niche breadth does not always facilitate range expansion through colonization and migration across barriers, though it can promote long-term persistence in situ.

Rapid Divergence Followed by Adaptation to Contrasting Ecological Niches of Two Closely Related Columbine Species Aquilegia japonica and A. oxysepala

Elucidating the mechanisms underlying the genetic divergence between closely related species is crucial to understanding the origin and evolution of biodiversity. The genus Aquilegia L. has undergone rapid adaptive radiation, generating about 70 well-recognized species that are specialized to distinct habitats and pollinators. In this study, to address the underlying evolutionary mechanisms that drive the genetic divergence, we analyzed the whole genomes of two ecologically isolated Aquilegia species, A. oxysepala and A. japonica as well as their putative hybrid. Our comparative genomic analyses reveal that while the two species diverged only recently and experienced recurrent gene flow, a high level of genetic divergence is observed in their nuclear genomes. In particular, candidate genomic regions that show signature of selection differ dramatically between the two species. Given that the splitting time of the two species is broadly matched with the decrease in effective population sizes, we propose that allopatric isolation together with natural selection have preceded the interspecific gene flow in the process of speciation. The observed high genetic divergence is likely an outcome of combined effects of natural selection, genetic drift and divergent sorting of ancestral polymorphisms. Our study provides a genome-wide view of how genetic divergence has evolved between closely related species.

Genome-wide variation in nucleotides and retrotransposons in alpine populations of Arabis alpina (Brassicaceae)

Advances in high-throughput sequencing have promoted the collection of reference genomes and genome-wide diversity. However, the assessment of genomic variation among populations has hitherto mainly been surveyed through single-nucleotide polymorphisms (SNPs) and largely ignored the often major fraction of genomes represented by transposable elements (TEs). Despite accumulating evidence supporting the evolutionary significance of TEs, comprehensive surveys remain scarce. Here, we sequenced the full genomes of 304 individuals of Arabis alpina sampled from four nearby natural populations to genotype SNPs as well as polymorphic long terminal repeat retrotransposons (polymorphic TEs; i.e., presence/absence of TE insertions at specific loci). We identified 291,396 SNPs and 20,548 polymorphic TEs, comparing their contributions to genomic diversity and divergence across populations. Few SNPs were shared among populations and overall showed high population-specific variation, whereas most polymorphic TEs segregated among populations. The genomic context of these two classes of variants further highlighted candidate adaptive loci having a putative impact on functional genes. In particular, 4.96% of the SNPs were identified as nonsynonymous or affecting start/stop codons. In contrast, 43% of the polymorphic TEs were present next to Arabis genes enriched in functional categories related to the regulation of reproduction and responses to biotic as well as abiotic stresses. This unprecedented data set, mapping variation gained from SNPs and complementary polymorphic TEs within and among populations, will serve as a rich resource for addressing microevolutionary processes shaping genome variation.

Using natural variation to understand the evolutionary pressures on plant photosynthesis

Photosynthesis is the gateway of the Sun's energy into the biosphere and the source of the ozone layer; thus it is both provider and protector of life as we know it. Despite its pivotal role we know surprisingly little about the genetic basis of variation in photosynthesis and the selective pressures giving rise to or maintaining this variation. In this review, I will briefly summarise our current knowledge of intraspecific and interspecific variation in photosynthesis to understand the main selective constraints on photosynthesis and what this means for the future of nature and agriculture in a changing world.

Natural variation of BSK3 tunes brassinosteroid signaling to regulate root foraging under low nitrogen

Developmental plasticity of root system architecture is crucial for plant performance in nutrient-poor soils. Roots of plants grown under mild nitrogen (N) deficiency show a foraging response characterized by increased root length but mechanisms underlying this developmental plasticity are still elusive. By employing natural variation in Arabidopsis accessions, we show that the brassinosteroid (BR) signaling kinase BSK3 modulates root elongation under mild N deficiency. In particular, a proline to leucine substitution in the predicted kinase domain of BSK3 enhances BR sensitivity and signaling to increase the extent of root elongation. We further show that low N specifically upregulates transcript levels of the BR co-receptor BAK1 to activate BR signaling and stimulate root elongation. Altogether, our results uncover a role of BR signaling in root elongation under low N. The BSK3 alleles identified here provide targets for improving root growth of crops growing under limited N conditions.

Plant roots elongate under mild nitrogen deficiency as part of a foraging response that facilitates nutrient uptake. Here the authors show that natural variation in this response among Arabidopsis accessions depends on the brassinosteroid (BR) signaling kinase BSK3, which can enhance BR sensitivity and root growth.

WhoGEM: an admixture-based prediction machine accurately predicts quantitative functional traits in plants

The explosive growth of genomic data provides an opportunity to make increased use of sequence variations for phenotype prediction. We have developed a prediction machine for quantitative phenotypes (WhoGEM) that overcomes some of the bottlenecks limiting the current methods. We demonstrated its performance by predicting quantitative disease resistance and quantitative functional traits in the wild model plant species, Medicago truncatula, using geographical locations as covariates for admixture analysis. The method's prediction reliability equals or outperforms all existing algorithms for quantitative phenotype prediction. WhoGEM analysis produces evidence that variation in genome admixture proportions explains most of the phenotypic variation for quantitative phenotypes.

Identification of Arabis alpina genomic regions associated with climatic variables along an elevation gradient through whole genome scan

We performed a pooled whole-genome sequencing on samples of the alpine plant Arabis alpina, harvested in ten populations along an elevation gradient in the French Alps. A large dataset of genetic variations was produced as single nucleotide polymorphisms (SNPs). A combined genome scan approach enabled detecting genomic regions associated with a synthetic environmental variable characterizing the climate at each sampling location. Positive loci detected by two methods were retained and belong to 19 regions in the Arabis alpina genome. The most significant region harbors an ortholog of the AtNAC062 gene, encoding a membrane-bound transcription factor described as linking the cold response and pathogen resistance that may confer protection to plants under extended snow coverage at high elevations. Other genes involved in the stress response or in flowering regulation were also detected. Altogether, our results indicated that Arabis alpina represent a suitable model for studying genomic adaptation in alpine perennial plants.

Natural allelic variation of GVS1 confers diversity in the regulation of leaf senescence in Arabidopsis

Leaf senescence affects plant fitness. Plants that evolve in different environments are expected to acquire distinct regulations of leaf senescence. However, the adaptive and evolutionary roles of leaf senescence are largely unknown. We investigated leaf senescence in 259 natural accessions of Arabidopsis by quantitatively assaying dark-induced senescence responses using a high-throughput chlorophyll fluorescence imaging system. A meta-analysis of our data with phenotypic and climatic information demonstrated biological and environmental links with leaf senescence. We further performed genome-wide association mapping to identify the genetic loci underlying the diversity of leaf senescence responses. We uncovered a new locus, Genetic Variants in leaf Senescence (GVS1), with high similarity to reductase, where a single nonsynonymous nucleotide substitution at GVS1 mediates the diversity of the senescence trait. Loss-of-function mutations of GVS1 in Columbia-0 delayed leaf senescence and increased sensitivity to oxidative stress, suggesting that this GVS1 variant promotes optimal responses to developmental and environmental signals. Intriguingly, gvs1 loss-of-function mutants display allele- and accession-dependent phenotypes, revealing the functional diversity of GVS1 alleles not only in leaf senescence, but also oxidative stress. Our discovery of GVS1 as the genetic basis of natural variation in senescence programs reinforces its adaptive potential in modulating life histories across diverse environments.

Genome-wide association study reveals new loci involved in Arabidopsis thaliana and Turnip mosaic virus (TuMV) interactions in the field

The genetic architecture of plant response to viruses has often been studied in model nonnatural pathosystems under controlled conditions. There is an urgent need to elucidate the genetic architecture of the response to viruses in a natural setting. A field experiment was performed in each of two years. In total, 317 Arabidopsis thaliana accessions were inoculated with its natural Turnip mosaic virus (TuMV). The accessions were phenotyped for viral accumulation, frequency of infected plants, stem length and symptoms. Genome-wide association mapping was performed. Arabidopsis thaliana exhibits extensive natural variation in its response to TuMV in the field. The underlying genetic architecture reveals a more quantitative picture than in controlled conditions. Ten genomic regions were consistently identified across the two years. RTM3 (Restricted TEV Movement 3) is a major candidate for the response to TuMV in the field. New candidate genes include Dead box helicase 1, a Tim Barrel domain protein and the eukaryotic translation initiation factor eIF3b. To our knowledge, this study is the first to report the genetic architecture of quantitative response of A. thaliana to a naturally occurring virus in a field environment, thereby highlighting relevant candidate genes involved in plant virus interactions in nature.

Soil microbes alter plant fitness under competition and drought

Plants exist across varying biotic and abiotic environments, including variation in the composition of soil microbial communities. The ecological effects of soil microbes on plant communities are well known, whereas less is known about their importance for plant evolutionary processes. In particular, the net effects of soil microbes on plant fitness may vary across environmental contexts and among plant genotypes, setting the stage for microbially mediated plant evolution. Here, we assess the effects of soil microbes on plant fitness and natural selection on flowering time in different environments. We performed two experiments in which we grew Arabidopsis thaliana genotypes replicated in either live or sterilized soil microbial treatments, and across varying levels of either competition (isolation, intraspecific competition or interspecific competition) or watering (well-watered or drought). We found large effects of competition and watering on plant fitness as well as the expression and natural selection of flowering time. Soil microbes increased average plant fitness under interspecific competition and drought and shaped the response of individual plant genotypes to drought. Finally, plant tolerance to either competition or drought was uncorrelated between soil microbial treatments suggesting that the plant traits favoured under environmental stress may depend on the presence of soil microbes. In summary, our experiments demonstrate that soil microbes can have large effects on plant fitness, which depend on both the environment and individual plant genotype. Future work in natural systems is needed for a complete understanding of the evolutionary importance of interactions between plants and soil microorganisms.

Linking genes with ecological strategies in Arabidopsis thaliana

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

Cis-regulated alternative splicing divergence and its potential contribution to environmental responses in Arabidopsis

Alternative splicing (AS) plays key roles in plant development and the responses of plants to environmental changes. However, the mechanisms underlying AS divergence (differential expression of transcript isoforms resulting from AS) in plant accessions and its contribution to responses to environmental stimuli remain unclear. In this study, we investigated genome-wide variation of AS in Arabidopsis thaliana accessions Col-0, Bur-0, C24, Kro-0 and Ler-1, as well as their F₁ hybrids, and characterized the regulatory mechanisms for AS divergence by RNA sequencing. We found that most of the divergent AS events in Arabidopsis accessions were cis-regulated by sequence variation, including those in core splice site and splicing motifs. Many genes that differed in AS between Col-0 and Bur-0 were involved in stimulus responses. Further genome-wide association analyses of 22 environmental variables showed that single nucleotide polymorphisms influencing known splice site strength were also associated with environmental stress responses. These results demonstrate that cis-variation in genomic sequences among Arabidopsis accessions was the dominant contributor to AS divergence, and it may contribute to differences in environmental responses among Arabidopsis accessions.

Phenotypic and genome-wide association with the local environment of Arabidopsis

The environment imposes critical selective forces on all living organisms, and the sessile nature of plants makes them particularly useful for investigating the relationship between genetic variation and environmental adaptation. In the model plant Arabidopsis thaliana, extensive information on phenotypic and genotypic variation is available, but comparable information on environmental variation within the native range of the species is lacking. Here, we compile 204 geoclimatic variables to describe the local environments of Arabidopsis accessions with known collection sites encompassing a wide geo-environmental range, and fully sequenced genomes from the 1001 Genomes Project. We identify candidate adaptive genetic variation associated with these environmental variables, and validate this approach through comparison with previous experimental studies, and by targeted confirmation of a role of the heterotrimeric G-protein γ subunit, AGG3, in cold tolerance, as newly predicted from our environmental genome wide association study (GWAS). To facilitate identification of adaptive variation, we created Arabidopsis CLIMtools : interactive web-based databases of the environment × genome associations and correlations between the local environments and 131 phenotypes compiled from previous experimental GWASs. Our study presents an extensive analysis of the local environments, landscape genomics and phenotypic variation of Arabidopsis, and illustrates how 'in silico GWAS' approaches can inform and complement experimental phenomics studies.

Adapting legume crops to climate change using genomic approaches

Our agricultural system and hence food security is threatened by combination of events, such as increasing population, the impacts of climate change, and the need to a more sustainable development. Evolutionary adaptation may help some species to overcome environmental changes through new selection pressures driven by climate change. However, success of evolutionary adaptation is dependent on various factors, one of which is the extent of genetic variation available within species. Genomic approaches provide an exceptional opportunity to identify genetic variation that can be employed in crop improvement programs. In this review, we illustrate some of the routinely used genomics-based methods as well as recent breakthroughs, which facilitate assessment of genetic variation and discovery of adaptive genes in legumes. Although additional information is needed, the current utility of selection tools indicate a robust ability to utilize existing variation among legumes to address the challenges of climate uncertainty.

Variation in ICE1 Methylation Primarily Determines Phenotypic Variation in Freezing Tolerance in Arabidopsis thaliana

Cold stress is a major abiotic factor plants face during their life cycle. Although plants often exhibit phenotypic variation in cold tolerance, the underlying mechanism remains poorly understood. In the present study, the 50% lethal temperature (LT50) values of 37 Arabidopsis thaliana accessions at latitudes from 15° to 58° ranged from -13.2°C to -4.9°C and were closely correlated with the cold climates of the collection sites. According to a methylation analysis of all C-repeat (CRT)-binding factor (CBF) pathway genes, the coding and promoter regions of AtICE1, a regulator of CBF genes, exhibited the greatest variability in methylation levels among the accessions and included 5-122 methylated cytosine residues. In contrast, unmethylated or only slightly methylated genes in the CBF pathway showed little variation among the accessions. According to a gene expression analysis of four selected A. thaliana populations with distinct methylation patterns, except for the down-regulated gene AtCBF2, the expression levels of all members of the CBF pathway were negatively correlated with AtICE1 gene methylation levels. Treatment of the four A. thaliana populations with the DNA methylation inhibitory reagent 5-azacytidine resulted in a 30.0-78.3% enhancement of freezing tolerance and decreases in LT50 values of approximately 1.9-3.6°C. Similar effects were observed in drm2 mutants, including 30.0-48.3% increases in freezing tolerance and decreases in LT50 values of approximately 0.7-3.4°C. Thus, the AtICE1 methylation-regulated transcription of CBF pathway genes is responsible for the phenotypic variation in the freezing tolerance observed in A. thaliana.

Evolutionary and ecological functional genomics, from lab to the wild

Plant phenotypes are the result of both genetic and environmental forces that act to modulate trait expression. Over the last few years, numerous approaches in functional genomics and systems biology have led to a greater understanding of plant phenotypic variation and plant responses to the environment. These approaches, and the questions that they can address, have been loosely termed evolutionary and ecological functional genomics (EEFG), and have been providing key insights on how plants adapt and evolve. In particular, by bringing these studies from the laboratory to the field, EEFG studies allow us to gain greater knowledge of how plants function in their natural contexts.

Phenotypic Trait Variation as a Response to Altitude-Related Constraints in Arabidopsis Populations

Natural variations help in identifying genetic mechanisms of morphologically and developmentally complex traits. Mountainous habitats provide an altitudinal gradient where one species encounters different abiotic conditions. We report the study of 341 individuals of Arabidopsis thaliana derived from 30 natural populations not belonging to the 1001 genomes, collected at increasing altitudes, between 200 and 1800 m in the Pyrenees. Class III peroxidases and ribosomal RNA sequences were used as markers to determine the putative genetic relationships among these populations along their altitudinal gradient. Using Bayesian-based statistics and phylogenetic analyses, these Pyrenean populations appear with significant divergence from the other regional accessions from 1001 genome (i.e., from north Spain or south France). Individuals of these populations exhibited varying phenotypic changes, when grown at sub-optimal temperature (22 vs. 15°C). These phenotypic variations under controlled conditions reflected intraspecific morphological variations. This study could bring new information regarding the west European population structure of A. thaliana and its phenotypic variations at different temperatures. The integrative analysis combining genetic, phenotypic variation and environmental datasets is used to analyze the acclimation of population in response to temperature changes. Regarding their geographical proximity and environmental diversity, these populations represent a tool of choice for studying plant response to temperature variation.

HIGHLIGHTS

-Studying the natural diversity of A. thaliana in the Pyrenees mountains helps to understand European population structure and to evaluate the phenotypic trait variation in response to climate change.

On the Causes of Rapid Diversification in the Páramos: Isolation by Ecology and Genomic Divergence in Espeletia

How diversity arises and what is the relative role of allopatric and ecological divergence are among the most persistent questions in evolution and ecology. Here, we assessed whether ecological divergence has enhanced the diversification of the Neotropical alpine plant complex Espeletia, also known as frailejones. This genus has one of the highest diversification rates ever reported and is distributed in the world's fastest evolving biodiversity hotspot, the Páramo (Neotropical alpine grasslands at elevations of c. 2800-4700 m). Our goal was to determine whether ecology plays a role in divergence within the Espeletia complex by quantifying genome-wide patterns of ecological divergence. We characterized 162 samples of the three most common and contrasting ecotypes (distinct morphotypes occupying particular habitats) co-occurring in six localities in the northern Andes using Genotyping by Sequencing. Contrasting ecotypes were caulescent cloud forest populations, caulescent populations from wind-sheltered and well-irrigated depressions and acaulescent populations from wind-exposed drier slopes. We found high polymorphism with a total of 1,273 single nucleotide polymorphisms (SNPs) that defined the relationships among nine genetic clusters. We quantified allelic associations of these markers with localities and habitats using 18 different general and mixed-effects statistical models that accounted for phylogenetic distance. Despite that these models always yielded more SNPs associated with the localities, markers associated with the habitat types were recovered too. We found strong evidence for isolation-by-distance (IBD) across populations despite rampant gene flow, as expected for plant groups with limited seed dispersal. Contrasts between populations of different habitat types showed that an isolation-by-environment (IBE) trend emerged and masked the IBD signal. Maximum likelihood estimation of the number of migrants per generation (Nem) among ecotypes confirmed the IBE pattern. This result illustrates the importance of mountains' environmental variation at a local scale in generating rapid morphological radiations and maintaining multiple adaptations in a fast-evolving ecosystem like the Páramo.

On the Road to Breeding 4.0: Unraveling the Good, the Bad, and the Boring of Crop Quantitative Genomics

Understanding the quantitative genetics of crops has been and will continue to be central to maintaining and improving global food security. We outline four stages that plant breeding either has already achieved or will probably soon achieve. Top-of-the-line breeding programs are currently in Breeding 3.0, where inexpensive, genome-wide data coupled with powerful algorithms allow us to start breeding on predicted instead of measured phenotypes. We focus on three major questions that must be answered to move from current Breeding 3.0 practices to Breeding 4.0: ( a) How do we adapt crops to better fit agricultural environments? ( b) What is the nature of the diversity upon which breeding can act? ( c) How do we deal with deleterious variants? Answering these questions and then translating them to actual gains for farmers will be a significant part of achieving global food security in the twenty-first century.

Plastic and genetic responses of a common sedge to warming have contrasting effects on carbon cycle processes

Climate warming affects plant physiology through genetic adaptation and phenotypic plasticity, but little is known about how these mechanisms influence ecosystem processes. We used three elevation gradients and a reciprocal transplant experiment to show that temperature causes genetic change in the sedge Eriophorum vaginatum. We demonstrate that plants originating from warmer climate produce fewer secondary compounds, grow faster and accelerate carbon dioxide (CO₂) release to the atmosphere. However, warmer climate also caused plasticity in E. vaginatum, inhibiting nitrogen metabolism, photosynthesis and growth and slowing CO₂ release into the atmosphere. Genetic differentiation and plasticity in E. vaginatum thus had opposing effects on CO₂ fluxes, suggesting that warming over many generations may buffer, or reverse, the short‐term influence of this species over carbon cycle processes. Our findings demonstrate the capacity for plant evolution to impact ecosystem processes, and reveal a further mechanism through which plants will shape ecosystem responses to climate change.

McSwan: A joint site frequency spectrum method to detect and date selective sweeps across multiple population genomes

Inferring the mode and tempo of natural selection helps further our understanding of adaptation to past environmental changes. Here, we introduce McSwan, a method to detect and date past and recent natural selection events in the case of a hard sweep. The method is based on the comparison of site frequency spectra obtained under various demographic models that include selection. McSwan demonstrated high power (high sensitivity and specificity) in capturing hard selective sweep events without requiring haplotype phasing. It performed slightly better than SweeD when the recent effective population size was low and the genomic region was small. We then applied our method to a European (CEU) and an African (LWK) human re-sequencing data set. Most hard sweeps were detected in the CEU population (96%). Moreover, hard sweeps in the African population were estimated to have occurred further back in time (mode: 43,625 years BP) compared to those of Europeans (mode: 24,850 years BP). Most of the estimated ages of hard sweeps in Europeans were associated with the Last Glacial Maximum and were enriched in immunity-associated genes.

RAD-Seq Reveals Patterns of Additive Polygenic Variation Caused by Spatially-Varying Selection in the American Eel (Anguilla rostrata)

The American Eel (Anguilla rostrata) has an exceptional life cycle characterized by panmictic reproduction at the species scale, random dispersal, and selection in a highly heterogeneous habitat extending from subtropical to subarctic latitudes. The genetic consequences of spatially-varying selection in this species have been investigated for decades, revealing subtle clines in allele frequency at a few loci that contrast with complete panmixia on the vast majority of the genome. Because reproduction homogenizes allele frequencies every generation, sampling size, and genomic coverage are critical to reach sufficient power to detect selected loci in this context. Here, we used a total of 710 individuals from 12 sites and 12,098 high-quality single nucleotide polymorphisms to re-evaluate the extent to which local selection affects the spatial distribution of genetic diversity in this species. We used environmental association methods to identify markers under spatially-varying selection, which indicated that selection affects ∼1.5% of the genome. We then evaluated the extent to which candidate markers collectively vary with environmental factors using additive polygenic scores. We found significant correlations between polygenic scores and latitude, longitude and temperature which are consistent with polygenic selection acting against maladapted genotypes in different habitats occupied by eels throughout their range of distribution. Gene functions associated with outlier markers were significantly enriched for the insulin signaling pathway, indicating that the trade-offs inherent to occupying such a large distribution range involve the regulation of metabolism. Overall, this study highlights the potential of the additive polygenic scores approach in detecting selective effects in a complex environment.

Population genomics and climate adaptation of a C4 perennial grass, Panicum hallii (Poaceae)

Background

Understanding how and why genetic variation is partitioned across geographic space is of fundamental importance to understanding the nature of biological species. How geographical isolation and local adaptation contribute to the formation of ecotypically differentiated groups of plants is just beginning to be understood through population genomic studies. We used whole genome sequencing combined with association study of climate to discover the drivers of differentiation in the perennial C4 grass Panicum hallii.

Results

Sequencing of 89 natural accessions of P.hallii revealed complex population structure across the species range. Major population genomic separation was found between subspecies P.hallii var. hallii and var. filipes as well as between at least four major unrecognized subgroups within var. hallii. At least 139 genomic SNPs were significantly associated with temperature or precipitation across the range and these SNPs were enriched for non-synonymous substitutions. SNPs associated with temperature and aridity were more often found in or near genes than expected by chance and enriched for putative involvement in dormancy processes, seed maturation, response to hyperosmosis and salinity, abscisic acid metabolism, hormone metabolism, and drought recovery.

Conclusions

Both geography and climate adaptation contribute significantly to patterns of genome-wide variation in P.hallii. Population subgroups within P.hallii may represent early stages in the formation of ecotypes. Climate associated loci identified here represent promising targets for future research in this and other perennial grasses.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5179-7) contains supplementary material, which is available to authorized users.

Genomic signatures of local adaptation to the degree of environmental predictability in rotifers

Environmental fluctuations are ubiquitous and thus essential for the study of adaptation. Despite this, genome evolution in response to environmental fluctuations —and more specifically to the degree of environmental predictability– is still unknown. Saline lakes in the Mediterranean region are remarkably diverse in their ecological conditions, which can lead to divergent local adaptation patterns in the inhabiting aquatic organisms. The facultatively sexual rotifer Brachionus plicatilis shows diverging local adaptation in its life-history traits in relation to estimated environmental predictability in its habitats. Here, we used an integrative approach —combining environmental, phenotypic and genomic data for the same populations– to understand the genomic basis of this diverging adaptation. Firstly, a novel draft genome for B. plicatilis was assembled. Then, genome-wide polymorphisms were studied using genotyping by sequencing on 270 clones from nine populations in eastern Spain. As a result, 4,543 high-quality SNPs were identified and genotyped. More than 90 SNPs were found to be putatively under selection with signatures of diversifying and balancing selection. Over 140 SNPs were correlated with environmental or phenotypic variables revealing signatures of local adaptation, including environmental predictability. Putative functions were associated to most of these SNPs, since they were located within annotated genes. Our results reveal associations between genomic variation and the degree of environmental predictability, providing genomic evidence of adaptation to local conditions in natural rotifer populations.