Adaptation to climate across the Arabidopsis thaliana genome

Modularity of genes involved in local adaptation to climate despite physical linkage

BACKGROUND

Linkage among genes experiencing different selection pressures can make natural selection less efficient. Theory predicts that when local adaptation is driven by complex and non-covarying stresses, increased linkage is favored for alleles with similar pleiotropic effects, with increased recombination favored among alleles with contrasting pleiotropic effects. Here, we introduce a framework to test these predictions with a co-association network analysis, which clusters loci based on differing associations. We use this framework to study the genetic architecture of local adaptation to climate in lodgepole pine, Pinus contorta, based on associations with environments.

RESULTS

We identify many clusters of candidate genes and SNPs associated with distinct environments, including aspects of aridity and freezing, and discover low recombination rates among some candidate genes in different clusters. Only a few genes contain SNPs with effects on more than one distinct aspect of climate. There is limited correspondence between co-association networks and gene regulatory networks. We further show how associations with environmental principal components can lead to misinterpretation. Finally, simulations illustrate both benefits and caveats of co-association networks.

CONCLUSIONS

Our results support the prediction that different selection pressures favor the evolution of distinct groups of genes, each associating with a different aspect of climate. But our results went against the prediction that loci experiencing different sources of selection would have high recombination among them. These results give new insight into evolutionary debates about the extent of modularity, pleiotropy, and linkage in the evolution of genetic architectures.

Comparative Genomics Approaches Accurately Predict Deleterious Variants in Plants

Recent advances in genome resequencing have led to increased interest in prediction of the functional consequences of genetic variants. Variants at phylogenetically conserved sites are of particular interest, because they are more likely than variants at phylogenetically variable sites to have deleterious effects on fitness and contribute to phenotypic variation. Numerous comparative genomic approaches have been developed to predict deleterious variants, but the approaches are nearly always assessed based on their ability to identify known disease-causing mutations in humans. Determining the accuracy of deleterious variant predictions in nonhuman species is important to understanding evolution, domestication, and potentially to improving crop quality and yield. To examine our ability to predict deleterious variants in plants we generated a curated database of 2,910 Arabidopsis thaliana mutants with known phenotypes. We evaluated seven approaches and found that while all performed well, their relative ranking differed from prior benchmarks in humans. We conclude that deleterious mutations can be reliably predicted in A. thaliana and likely other plant species, but that the relative performance of various approaches does not necessarily translate from one species to another.

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

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

Geopolitical species revisited: genomic and morphological data indicate that the roundtail chub Gila robusta species complex (Teleostei, Cyprinidae) is a single species

The Gila robusta species complex in the Lower Colorado River Basin has a complicated taxonomic history. Recent authors have separated this group into three nominal taxa, G. robusta, G. intermedia, and G. nigra, however aside from location, no reliable method of distinguishing individuals of these species currently exists. To assess relationships within this group, we examined morphology of type specimens and fresh material, and used RADseq methods to assess phylogenetic relationship among these nominal species. Maximum likelihood and Bayesian inference tree building methods reveal high concordance between tree topologies based on the mitochondrial and nuclear datasets. Coalescent SNAPP analysis resolved a similar tree topology. Neither morphological nor molecular data reveal diagnostic differences between these species as currently defined. As such, G. intermedia and G. nigra should be considered synonyms of the senior G. robusta. We hypothesize that climate driven wet and dry cycles have led to periodic isolation of population subunits and subsequent local divergence followed by reestablished connectivity and mixing. Management plans should therefore focus on retaining genetic variability and viability of geographic populations to preserve adaptability to changing climate conditions.

The genomic basis of environmental adaptation in house mice

House mice (Mus musculus) arrived in the Americas only recently in association with European colonization (~400-600 generations), but have spread rapidly and show evidence of local adaptation. Here, we take advantage of this genetic model system to investigate the genomic basis of environmental adaptation in house mice. First, we documented clinal patterns of phenotypic variation in 50 wild-caught mice from a latitudinal transect in Eastern North America. Next, we found that progeny of mice from different latitudes, raised in a common laboratory environment, displayed differences in a number of complex traits related to fitness. Consistent with Bergmann's rule, mice from higher latitudes were larger and fatter than mice from lower latitudes. They also built bigger nests and differed in aspects of blood chemistry related to metabolism. Then, combining exomic, genomic, and transcriptomic data, we identified specific candidate genes underlying adaptive variation. In particular, we defined a short list of genes with cis-eQTL that were identified as candidates in exomic and genomic analyses, all of which have known ties to phenotypes that vary among the studied populations. Thus, wild mice and the newly developed strains represent a valuable resource for future study of the links between genetic variation, phenotypic variation, and climate.

Assortment of Flowering Time and Immunity Alleles in Natural Arabidopsis thaliana Populations Suggests Immunity and Vegetative Lifespan Strategies Coevolve

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

Genome-wide signatures of flowering adaptation to climate temperature: Regional analyses in a highly diverse native range of Arabidopsis thaliana

Current global change is fueling an interest to understand the genetic and molecular mechanisms of plant adaptation to climate. In particular, altered flowering time is a common strategy for escape from unfavourable climate temperature. In order to determine the genomic bases underlying flowering time adaptation to this climatic factor, we have systematically analysed a collection of 174 highly diverse Arabidopsis thaliana accessions from the Iberian Peninsula. Analyses of 1.88 million single nucleotide polymorphisms provide evidence for a spatially heterogeneous contribution of demographic and adaptive processes to geographic patterns of genetic variation. Mountains appear to be allele dispersal barriers, whereas the relationship between flowering time and temperature depended on the precise temperature range. Environmental genome-wide associations supported an overall genome adaptation to temperature, with 9.4% of the genes showing significant associations. Furthermore, phenotypic genome-wide associations provided a catalogue of candidate genes underlying flowering time variation. Finally, comparison of environmental and phenotypic genome-wide associations identified known (Twin Sister of FT, FRIGIDA-like 1, and Casein Kinase II Beta chain 1) and new (Epithiospecifer Modifier 1 and Voltage-Dependent Anion Channel 5) genes as candidates for adaptation to climate temperature by altered flowering time. Thus, this regional collection provides an excellent resource to address the spatial complexity of climate adaptation in annual plants.

Phenotypic plasticity in response to temperature fluctuations is genetically variable, and relates to climatic variability of origin, in Arabidopsis thaliana

Under current climate change, increasing mean temperatures are not only causing hotter summers, but temperature variability is increasing as well. Phenotypic plasticity can help plants to overcome negative effects of temperature variability and allow them to rapidly adjust traits to adverse conditions. Moreover, genetic variation in such plasticity could provide potential for adaptive evolution in response to changing climate variability. Here, we conducted an experiment with 11 Arabidopsis thaliana genotypes to investigate intraspecific variation in plant responses to two aspects of variable temperature stress: timing and frequency. We found that the timing but not frequency of temperature stress affected the phenology, growth, reproduction and allocation strategy of plants, and that genotypes differed substantially in their responses. Moreover, trait plasticity was positively related to precipitation variability of origin, suggesting an adaptive role of plasticity. Our results indicate that the developmental stage of a plant during heat stress is a key determinant of its response, and that plasticity to temperature variability is an evolving and possibly adaptive trait in natural populations of A. thaliana. More generally, our study demonstrates the usefulness of studying plant responses to climatic variability per se, given that climatic variability is predicted to increase in the future.

Measuring the dynamic photosynthome

Background

Photosynthesis underpins plant productivity and yet is notoriously sensitive to small changes in environmental conditions, meaning that quantitation in nature across different time scales is not straightforward. The 'dynamic' changes in photosynthesis (i.e. the kinetics of the various reactions of photosynthesis in response to environmental shifts) are now known to be important in driving crop yield.

Scope

It is known that photosynthesis does not respond in a timely manner, and even a small temporal 'mismatch' between a change in the environment and the appropriate response of photosynthesis toward optimality can result in a fall in productivity. Yet the most commonly measured parameters are still made at steady state or a temporary steady state (including those for crop breeding purposes), meaning that new photosynthetic traits remain undiscovered.

Conclusions

There is a great need to understand photosynthesis dynamics from a mechanistic and biological viewpoint especially when applied to the field of 'phenomics' which typically uses large genetically diverse populations of plants. Despite huge advances in measurement technology in recent years, it is still unclear whether we possess the capability of capturing and describing the physiologically relevant dynamic features of field photosynthesis in sufficient detail. Such traits are highly complex, hence we dub this the 'photosynthome'. This review sets out the state of play and describes some approaches that could be made to address this challenge with reference to the relevant biological processes involved.

Influence of intergenotypic competition on multigenerational persistence of abiotic stress resistance transgenes in populations of Arabidopsis thaliana

Reducing crop losses due to abiotic stresses is a major target of agricultural biotechnology that will increase with climate change and global population growth. Concerns, however, have been raised about potential ecological impacts if transgenes become established in wild populations and cause increased competitiveness of weedy or invasive species. Potential risks will be a function of transgene movement, population sizes, and fitness effects on the recipient population. While key components influencing gene flow have been extensively investigated, there have been few studies on factors subsequent to transgene movement that can influence persistence and competitiveness. Here, we performed multiyear, multigenerational, assessment to examine fitness effects and persistence of three mechanistically different abiotic stress tolerance genes: C-repeat binding factor 3/drought responsive element binding factor 1a (CBF3/DREB1a); Salt overly sensitive 1 (SOS1); and Mannose-6-phosphate reductase (M6PR). Transgenic Arabidopsis thaliana overexpressing these genes were grown in pure populations and in competition with wild-type (WT) parents for six generations spanning a range of field environment conditions. Growth, development, biomass, seed production, and transgene frequency were measured at each generation. Seed planted for each generation was obtained from the previous generation as would occur during establishment of a new genotype in the environment. The three transgenes exhibited different fitness effects and followed different establishment trajectories. In comparison with pure populations, CBF3 lines exhibited reduced dry weight, seed yield, and viable seed yield, relative to WT background. In contrast, overexpression of SOS1 and M6PR did not significantly impact productivity measures in pure populations. In competition with WT, negative fitness effects were magnified. Transgene frequencies were significantly reduced for CBF3 and SOS1 while frequencies of M6PR appeared to be subject to genetic drift. These studies demonstrate the importance of fitness effects and intergenotype competition in influencing persistence of transgenes conferring complex traits.

Two-way mixed-effects methods for joint association analysis using both host and pathogen genomes

Infectious diseases are often affected by specific pairings of hosts and pathogens and therefore by both of their genomes. The integration of a pair of genomes into genome-wide association mapping can provide an exquisitely detailed view of the genetic landscape of complex traits. We present a statistical method, ATOMM (Analysis with a Two-Organism Mixed Model), that maps a trait of interest to a pair of genomes simultaneously; this method makes use of whole-genome sequence data for both host and pathogen organisms. ATOMM uses a two-way mixed-effect model to test for genetic associations and cross-species genetic interactions while accounting for sample structure including interactions between the genetic backgrounds of the two organisms. We demonstrate the applicability of ATOMM to a joint association study of quantitative disease resistance (QDR) in the Arabidopsis thaliana-Xanthomonas arboricola pathosystem. Our method uncovers a clear host-strain specificity in QDR and provides a powerful approach to identify genetic variants on both genomes that contribute to phenotypic variation.

Genomic signals of selection predict climate-driven population declines in a migratory bird

The ongoing loss of biodiversity caused by rapid climatic shifts requires accurate models for predicting species' responses. Despite evidence that evolutionary adaptation could mitigate climate change impacts, evolution is rarely integrated into predictive models. Integrating population genomics and environmental data, we identified genomic variation associated with climate across the breeding range of the migratory songbird, yellow warbler (Setophaga petechia). Populations requiring the greatest shifts in allele frequencies to keep pace with future climate change have experienced the largest population declines, suggesting that failure to adapt may have already negatively affected populations. Broadly, our study suggests that the integration of genomic adaptation can increase the accuracy of future species distribution models and ultimately guide more effective mitigation efforts.

Linking genome signatures of selection and adaptation in non-model plants: exploring potential and limitations in the angiosperm Amborella

Selective sweeps may be caused by environmental conditions that select for a gene function or trait at one locus, causing reduced variation at neighboring sites due to linkage, with specific non-selected variants being swept along with the selected variant. For many species, genomic and environmental data are available to test hypotheses that environmental conditions are correlated with selected regions. Most genomic studies relating selection to environment use model organisms or crop species; typically, these studies have genomic data from large numbers of individuals and extensive environmental data. Here, we review studies associating selective sweeps with environment and consider the impediments to successful application of these methods to non-model species. We present an initial investigation into linking genomic regions of selection to environmental conditions in the narrowly distributed, non-model plant Amborella trichopoda (Amborellaceae), the sister species to all other living flowering plants and one of over 2500 plant species endemic to New Caledonia.

Adaptive diversification of growth allometry in the plant Arabidopsis thaliana

Seed plants vary tremendously in size and morphology; however, variation and covariation in plant traits may be governed, at least in part, by universal biophysical laws and biological constants. Metabolic scaling theory (MST) posits that whole-organismal metabolism and growth rate are under stabilizing selection that minimizes the scaling of hydrodynamic resistance and maximizes the scaling of resource uptake. This constrains variation in physiological traits and in the rate of biomass accumulation, so that they can be expressed as mathematical functions of plant size with near-constant allometric scaling exponents across species. However, the observed variation in scaling exponents calls into question the evolutionary drivers and the universality of allometric equations. We have measured growth scaling and fitness traits of 451 Arabidopsis thaliana accessions with sequenced genomes. Variation among accessions around the scaling exponent predicted by MST was correlated with relative growth rate, seed production, and stress resistance. Genomic analyses indicate that growth allometry is affected by many genes associated with local climate and abiotic stress response. The gene with the strongest effect, PUB4, has molecular signatures of balancing selection, suggesting that intraspecific variation in growth scaling is maintained by opposing selection on the trade-off between seed production and abiotic stress resistance. Our findings suggest that variation in allometry contributes to local adaptation to contrasting environments. Our results help reconcile past debates on the origin of allometric scaling in biology and begin to link adaptive variation in allometric scaling to specific genes.

Navigating the Interface Between Landscape Genetics and Landscape Genomics

As next-generation sequencing data become increasingly available for non-model organisms, a shift has occurred in the focus of studies of the geographic distribution of genetic variation. Whereas landscape genetics studies primarily focus on testing the effects of landscape variables on gene flow and genetic population structure, landscape genomics studies focus on detecting candidate genes under selection that indicate possible local adaptation. Navigating the transition between landscape genomics and landscape genetics can be challenging. The number of molecular markers analyzed has shifted from what used to be a few dozen loci to thousands of loci and even full genomes. Although genome scale data can be separated into sets of neutral loci for analyses of gene flow and population structure and putative loci under selection for inference of local adaptation, there are inherent differences in the questions that are addressed in the two study frameworks. We discuss these differences and their implications for study design, marker choice and downstream analysis methods. Similar to the rapid proliferation of analysis methods in the early development of landscape genetics, new analytical methods for detection of selection in landscape genomics studies are burgeoning. We focus on genome scan methods for detection of selection, and in particular, outlier differentiation methods and genetic-environment association tests because they are the most widely used. Use of genome scan methods requires an understanding of the potential mismatches between the biology of a species and assumptions inherent in analytical methods used, which can lead to high false positive rates of detected loci under selection. Key to choosing appropriate genome scan methods is an understanding of the underlying demographic structure of study populations, and such data can be obtained using neutral loci from the generated genome-wide data or prior knowledge of a species' phylogeographic history. To this end, we summarize recent simulation studies that test the power and accuracy of genome scan methods under a variety of demographic scenarios and sampling designs. We conclude with a discussion of additional considerations for future method development, and a summary of methods that show promise for landscape genomics studies but are not yet widely used.

Genetic diversity of Avena ventricosa populations along an ecogeographical transect in Cyprus is correlated to environmental variables

Avena ventricosa Balansa ex Coss. is considered the C-genome donor of the cultivated hexaploid oat and is a ‘priority’ species for conservation, since it has limited geographic distribution and the only recorded populations in Europe are present in Cyprus. The current study attempts to characterize the genetic structure and fragmentation of the species via the application of genotypic markers. It was revealed that the genetic variety was mainly allocated among the populations collected, since clustering obtained was according to the geographic origin of the samples and the habitat. Species distribution modeling showed that the most important climatic variable defining A. ventricosa distribution is the mean diurnal temperature. Furthermore, significant association of the genetic structure to environmental variables was detected; overall, a negative association to precipitation was confirmed, while significant correlations of genetic structure and the temperature at the time of anthesis and germination were established. The safeguarding of this valuable genetic resource is discussed.

AtHMA4 Drives Natural Variation in Leaf Zn Concentration of Arabidopsis thaliana

Zinc (Zn) is an essential element for plant growth and development, and Zn derived from crop plants in the diet is also important for human health. Here, we report that genetic variation in Heavy Metal-ATPase 4 (HMA4) controls natural variation in leaf Zn content. Investigation of the natural variation in leaf Zn content in a world-wide collection of 349 Arabidopsis thaliana wild collected accessions identified two accessions, Van-0 and Fab-2, which accumulate significantly lower Zn when compared with Col-0. Both quantitative trait loci (QTL) analysis and bulked segregant analysis (BSA) identified HMA4 as a strong candidate accounting for this variation in leaf Zn concentration. Genetic complementation experiments confirmed this hypothesis. Sequence analysis revealed that a 1-bp deletion in the third exon of HMA4 from Fab-2 is responsible for the lose of function of HMA4 driving the low Zn observed in Fab-2. Unlike in Fab-2 polymorphisms in the promoter region were found to be responsible for the weak function of HMA4 in Van-0. This is supported by both an expression analysis of HMA4 in Van-0 and through a series of T-DNA insertion mutants which generate truncated HMA4 promoters in the Col-0 background. In addition, we also observed that Fab-2, Van-0 and the hma4-2 null mutant in the Col-0 background show enhanced resistance to a combination of high Zn and high Cd in the growth medium, raising the possibility that variation at HMA4 may play a role in environmental adaptation.

Genotyping by Sequencing and Genome-Environment Associations in Wild Common Bean Predict Widespread Divergent Adaptation to Drought

Drought will reduce global crop production by >10% in 2050 substantially worsening global malnutrition. Breeding for resistance to drought will require accessing crop genetic diversity found in the wild accessions from the driest high stress ecosystems. Genome-environment associations (GEA) in crop wild relatives reveal natural adaptation, and therefore can be used to identify adaptive variation. We explored this approach in the food crop Phaseolus vulgaris L., characterizing 86 geo-referenced wild accessions using genotyping by sequencing (GBS) to discover single nucleotide polymorphisms (SNPs). The wild beans represented Mesoamerica, Guatemala, Colombia, Ecuador/Northern Peru and Andean groupings. We found high polymorphism with a total of 22,845 SNPs across the 86 accessions that confirmed genetic relationships for the groups. As a second objective, we quantified allelic associations with a bioclimatic-based drought index using 10 different statistical models that accounted for population structure. Based on the optimum model, 115 SNPs in 90 regions, widespread in all 11 common bean chromosomes, were associated with the bioclimatic-based drought index. A gene coding for an ankyrin repeat-containing protein and a phototropic-responsive NPH3 gene were identified as potential candidates. Genomic windows of 1 Mb containing associated SNPs had more positive Tajima's D scores than windows without associated markers. This indicates that adaptation to drought, as estimated by bioclimatic variables, has been under natural divergent selection, suggesting that drought tolerance may be favorable under dry conditions but harmful in humid conditions. Our work exemplifies that genomic signatures of adaptation are useful for germplasm characterization, potentially enhancing future marker-assisted selection and crop improvement.

New insights into the regulation of leaf senescence in Arabidopsis

Plants undergo developmental changes throughout their life history. Senescence, the final stage in the life history of a leaf, is an important and unique developmental process whereby plants relocate nutrients from leaves to other developing organs, such as seeds, stems, or roots. Recent attempts to answer fundamental questions about leaf senescence have employed a combination of new ideas and advanced technologies. As senescence is an integral part of a plant's life history that is linked to earlier developmental stages, age-associated leaf senescence may be analysed from a life history perspective. The successful utilization of multi-omics approaches has resolved the complicated process of leaf senescence, replacing a component-based view with a network-based molecular mechanism that acts in a spatial-temporal manner. Senescence and death are critical for fitness and are thus evolved characters. Recent efforts have begun to focus on understanding the evolutionary basis of the developmental process that incorporates age information and environmental signals into a plant's survival strategy. This review describes recent insights into the regulatory mechanisms of leaf senescence in terms of systems-level spatiotemporal changes, presenting them from the perspectives of life history strategy and evolution.

Landscape genomics and pathway analysis to understand genetic adaptation of South African indigenous goat populations

In Africa, extensively raised livestock populations in most smallholder farming communities are exposed to harsh and heterogeneous climatic conditions and disease pathogens that they adapt to in order to survive. Majority of these livestock species, including goats, are of non-descript and uncharacterized breeds and their response to natural selection presented by heterogeneous environments is still unresolved. This study investigated genetic diversity and its association with environmental and geographic conditions in 194 South African indigenous goats from different geographic locations genotyped on the Illumina goat SNP50K panel. Population structure analysis revealed a homogeneous genetic cluster of the Tankwa goats, restricted to the Northern Cape province. Overall, the Boer, Kalahari Red, and Savanna showed a wide geographic spread of shared genetic components, whereas the village ecotypes revealed a longitudinal distribution. The relative importance of environmental factors on genetic variation of goat populations was assessed using redundancy analysis (RDA). Climatic and geographic variables explained 22% of the total variation while climatic variables alone accounted for 17% of the diversity. Geographic variables solitarily explained 1% of the total variation. The first axis (Model I) of the RDA analysis revealed 329 outlier SNPs. Landscape genomic approaches of spatial analysis method (SAM) identified a total of 843 (1.75%) SNPs, while latent factor mixed models (LFMM) identified 714 (1.48%) SNPs significantly associated with environmental variables. Significant markers were within genes involved in biological functions potentially important for environmental adaptation. Overall, the study suggested environmental factors to have some effect in shaping the genetic variation of South African indigenous goat populations. Loci observed to be significant and under selection may be responsible for the adaption of the goat populations to local production systems.

The genetics underlying natural variation of plant-plant interactions, a beloved but forgotten member of the family of biotic interactions

Despite the importance of plant-plant interactions on crop yield and plant community dynamics, our understanding of the genetic and molecular bases underlying natural variation of plant-plant interactions is largely limited in comparison with other types of biotic interactions. By listing 63 quantitative trait loci (QTL) mapping and global gene expression studies based on plants directly challenged by other plants, we explored whether the genetic architecture and the function of the candidate genes underlying natural plant-plant interactions depend on the type of interactions between two plants (competition versus commensalism versus reciprocal helping versus asymmetry). The 16 transcriptomic studies are unevenly distributed between competitive interactions (n = 12) and asymmetric interactions (n = 4, all focusing on response to parasitic plants). By contrast, 17 and 30 QTL studies were identified for competitive interactions and asymmetric interactions (either weed suppressive ability or response to parasitic plants), respectively. Surprisingly, no studies have been carried out on the identification of genetic and molecular bases underlying natural variation in positive interactions. The candidate genes underlying natural plant-plant interactions can be classified into seven categories of plant function that have been identified in artificial environments simulating plant-plant interactions either frequently (photosynthesis, hormones), only recently (cell wall modification and degradation, defense pathways against pathogens) or rarely (ABC transporters, histone modification and meristem identity/life history traits). Finally, we introduce several avenues that need to be explored in the future to obtain a thorough understanding of the genetic and molecular bases underlying plant-plant interactions within the context of realistic community complexity.

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

By following the evolution of populations that are initially genetically homogeneous, much can be learned about core biological principles. For example, it allows for detailed studies of the rate of emergence of de novo mutations and their change in frequency due to drift and selection. Unfortunately, in multicellular organisms with generation times of months or years, it is difficult to set up and carry out such experiments over many generations. An alternative is provided by "natural evolution experiments" that started from colonizations or invasions of new habitats by selfing lineages. With limited or missing gene flow from other lineages, new mutations and their effects can be easily detected. North America has been colonized in historic times by the plant Arabidopsis thaliana, and although multiple intercrossing lineages are found today, many of the individuals belong to a single lineage, HPG1. To determine in this lineage the rate of substitutions-the subset of mutations that survived natural selection and drift-, we have sequenced genomes from plants collected between 1863 and 2006. We identified 73 modern and 27 herbarium specimens that belonged to HPG1. Using the estimated substitution rate, we infer that the last common HPG1 ancestor lived in the early 17th century, when it was most likely introduced by chance from Europe. Mutations in coding regions are depleted in frequency compared to those in other portions of the genome, consistent with purifying selection. Nevertheless, a handful of mutations is found at high frequency in present-day populations. We link these to detectable phenotypic variance in traits of known ecological importance, life history and growth, which could reflect their adaptive value. Our work showcases how, by applying genomics methods to a combination of modern and historic samples from colonizing lineages, we can directly study new mutations and their potential evolutionary relevance.

A Genomic Map of Climate Adaptation in Arabidopsis thaliana at a Micro-Geographic Scale

Understanding the genetic bases underlying climate adaptation is a key element to predict the potential of species to face climate warming. Although substantial climate variation is observed at a micro-geographic scale, most genomic maps of climate adaptation have been established at broader geographical scales. Here, by using a Pool-Seq approach combined with a Bayesian hierarchical model that control for confounding by population structure, we performed a genome-environment association (GEA) analysis to investigate the genetic basis of adaptation to six climate variables in 168 natural populations of Arabidopsis thaliana distributed in south-west of France. Climate variation among the 168 populations represented up to 24% of climate variation among 521 European locations where A. thaliana inhabits. We identified neat and strong peaks of association, with most of the associated SNPs being significantly enriched in likely functional variants and/or in the extreme tail of genetic differentiation among populations. Furthermore, genes involved in transcriptional mechanisms appear predominant in plant functions associated with local climate adaptation. Globally, our results suggest that climate adaptation is an important driver of genomic variation in A. thaliana at a small spatial scale and mainly involves genome-wide changes in fundamental mechanisms of gene regulation. The identification of climate-adaptive genetic loci at a micro-geographic scale also highlights the importance to include within-species genetic diversity in ecological niche models for projecting potential species distributional shifts over short geographic distances.

Eco-Metabolomics and Metabolic Modeling: Making the Leap From Model Systems in the Lab to Native Populations in the Field

Experimental high-throughput analysis of molecular networks is a central approach to characterize the adaptation of plant metabolism to the environment. However, recent studies have demonstrated that it is hardly possible to predict in situ metabolic phenotypes from experiments under controlled conditions, such as growth chambers or greenhouses. This is particularly due to the high molecular variance of in situ samples induced by environmental fluctuations. An approach of functional metabolome interpretation of field samples would be desirable in order to be able to identify and trace back the impact of environmental changes on plant metabolism. To test the applicability of metabolomics studies for a characterization of plant populations in the field, we have identified and analyzed in situ samples of nearby grown natural populations of Arabidopsis thaliana in Austria. A. thaliana is the primary molecular biological model system in plant biology with one of the best functionally annotated genomes representing a reference system for all other plant genome projects. The genomes of these novel natural populations were sequenced and phylogenetically compared to a comprehensive genome database of A. thaliana ecotypes. Experimental results on primary and secondary metabolite profiling and genotypic variation were functionally integrated by a data mining strategy, which combines statistical output of metabolomics data with genome-derived biochemical pathway reconstruction and metabolic modeling. Correlations of biochemical model predictions and population-specific genetic variation indicated varying strategies of metabolic regulation on a population level which enabled the direct comparison, differentiation, and prediction of metabolic adaptation of the same species to different habitats. These differences were most pronounced at organic and amino acid metabolism as well as at the interface of primary and secondary metabolism and allowed for the direct classification of population-specific metabolic phenotypes within geographically contiguous sampling sites.

High-frequency temperature variability mirrors fixed differences in thermal limits of the massive coral Porites lobata (Dana, 1846)

Spatial heterogeneity in environmental characteristics can drive adaptive differentiation when contrasting environments exert divergent selection pressures. This environmental and genetic heterogeneity can substantially influence population and community resilience to disturbance events. Here, we investigated corals from the highly variable back reef habitats of Ofu Island in American Samoa that thrive in thermal conditions known to elicit widespread bleaching and mortality elsewhere. To investigate the relative importance of acclimation vs. site of origin in shaping previously observed differences in coral tolerance limits in Ofu, specimens of the common Indo-Pacific coral Porites lobata (Dana, 1846) from locations with differing levels of thermal variability were acclimated to low and high thermal variation in controlled common garden aquaria. Overall, there were minimal effects of the acclimation exposure. Corals native to the site with the highest level of daily variability grew fastest, regardless of acclimation treatment. When exposed to lethal thermal stress, corals native to both variable sites contained elevated levels of heat shock proteins and maintained photosynthetic performance for 1-2 days longer than corals from the stable environment. Despite being separated by<5 km, there was significant genetic differentiation among coral colonies (FST=0.206, p<0.0001; nuclear ribosomal DNA), while Symbiodiniaceae were all Cladocopium sp. (ITS-type C15). Our study demonstrates consistent signatures of adaptation in growth and stress resistance in corals from naturally thermally variable habitats, suggesting that differences in the amount of thermal variability may be an important contributor to adaptive differentiation in reef building corals.

Temperature fine-tunes Mediterranean Arabidopsis thaliana life-cycle phenology geographically

To understand how adaptive evolution in life-cycle phenology operates in plants, we need to unravel the effects of geographic variation in putative agents of natural selection on life-cycle phenology by considering all key developmental transitions and their co-variation patterns. We address this goal by quantifying the temperature-driven and geographically varying relationship between seed dormancy and flowering time in the annual Arabidopsis thaliana across the Iberian Peninsula. We used data on genetic variation in two major life-cycle traits, seed dormancy (DSDS50) and flowering time (FT), in a collection of 300 A. thaliana accessions from the Iberian Peninsula. The geographically varying relationship between life-cycle traits and minimum temperature, a major driver of variation in DSDS50 and FT, was explored with geographically weighted regressions (GWR). The environmentally varying correlation between DSDS50 and FT was analysed by means of sliding window analysis across a minimum temperature gradient. Maximum local adjustments between minimum temperature and life-cycle traits were obtained in the southwest Iberian Peninsula, an area with the highest minimum temperatures. In contrast, in off-southwest locations, the effects of minimum temperature on DSDS50 were rather constant across the region, whereas those of minimum temperature on FT were more variable, with peaks of strong local adjustments of GWR models in central and northwest Spain. Sliding window analysis identified a minimum temperature turning point in the relationship between DSDS50 and FT around a minimum temperature of 7.2 °C. Above this minimum temperature turning point, the variation in the FT/DSDS50 ratio became rapidly constrained and the negative correlation between FT and DSDS50 did not increase any further with increasing minimum temperatures. The southwest Iberian Peninsula emerges as an area where variation in life-cycle phenology appears to be restricted by the duration and severity of the hot summer drought. The temperature-driven varying relationship between DSDS50 and FT detected environmental boundaries for the co-evolution between FT and DSDS50 in A. thaliana. In the context of global warming, we conclude that A. thaliana phenology from the southwest Iberian Peninsula, determined by early flowering and deep seed dormancy, might become the most common life-cycle phenotype for this annual plant in the region.

Signatures of natural selection in abiotic stress-responsive genes of Solanum chilense

Environmental conditions are strong selective forces, which may influence adaptation and speciation. The wild tomato species Solanum chilense, native to South America, is exposed to a range of abiotic stress factors. To identify signatures of natural selection and local adaptation, we analysed 16 genes involved in the abiotic stress response and compared the results to a set of reference genes in 23 populations across the entire species range. The abiotic stress-responsive genes are characterized by elevated nonsynonymous nucleotide diversity and divergence. We detected signatures of positive selection in several abiotic stress-responsive genes on both the population and species levels. Local adaptation to abiotic stresses is particularly apparent at the boundary of the species distribution in populations from coastal low-altitude and mountainous high-altitude regions.

Adaptation of Arabidopsis thaliana to the Yangtze River basin

Background

Organisms need to adapt to keep pace with a changing environment. Examining recent range expansion aids our understanding of how organisms evolve to overcome environmental constraints. However, how organisms adapt to climate changes is a crucial biological question that is still largely unanswered. The plant Arabidopsis thaliana is an excellent system to study this fundamental question. Its origin is in the Iberian Peninsula and North Africa, but it has spread to the Far East, including the most south-eastern edge of its native habitats, the Yangtze River basin, where the climate is very different.

Results

We sequenced 118 A. thaliana strains from the region surrounding the Yangtze River basin. We found that the Yangtze River basin population is a unique population and diverged about 61,409 years ago, with gene flows occurring at two different time points, followed by a population dispersion into the Yangtze River basin in the last few thousands of years. Positive selection analyses revealed that biological regulation processes, such as flowering time, immune and defense response processes could be correlated with the adaptation event. In particular, we found that the flowering time gene SVP has contributed to A. thaliana adaptation to the Yangtze River basin based on genetic mapping.

Conclusions

A. thaliana adapted to the Yangtze River basin habitat by promoting the onset of flowering, a finding that sheds light on how a species can adapt to locales with very different climates.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1378-9) contains supplementary material, which is available to authorized users.

Genomic basis and evolutionary potential for extreme drought adaptation in Arabidopsis thaliana

As Earth is currently experiencing dramatic climate change, it is of critical interest to understand how species will respond to it. The chance of a species withstanding climate change is likely to depend on the diversity within the species and, particularly, whether there are sub-populations that are already adapted to extreme environments. However, most predictive studies ignore that species comprise genetically diverse individuals. We have identified genetic variants in Arabidopsis thaliana that are associated with survival of an extreme drought event-a major consequence of global warming. Subsequently, we determined how these variants are distributed across the native range of the species. Genetic alleles conferring higher drought survival showed signatures of polygenic adaptation and were more frequently found in Mediterranean and Scandinavian regions. Using geo-environmental models, we predicted that Central European, but not Mediterranean, populations might lag behind in adaptation by the end of the twenty-first century. Further analyses showed that a population decline could nevertheless be compensated by natural selection acting efficiently over standing variation or by migration of adapted individuals from populations at the margins of the species' distribution. These findings highlight the importance of within-species genetic heterogeneity in facilitating an evolutionary response to a changing climate.

Dissecting the Genetic Basis of Local Adaptation in Soybean

Soybean (Glycine max) is the most widely grown oilseed in the world and is an important source of protein for both humans and livestock. Soybean is widely adapted to both temperate and tropical regions, but a changing climate demands a better understanding of adaptation to specific environmental conditions. Here, we explore genetic variation in a collection of 3,012 georeferenced, locally adapted landraces from a broad geographical range to help elucidate the genetic basis of local adaptation. We used geographic origin, environmental data and dense genome-wide SNP data to perform an environmental association analysis and discover loci displaying steep gradients in allele frequency across geographical distance and between landrace and modern cultivars. Our combined application of methods in environmental association mapping and detection of selection targets provide a better understanding of how geography and selection may have shaped genetic variation among soybean landraces. Moreover, we identified several important candidate genes related to drought and heat stress, and revealed important genomic regions possibly involved in the geographic divergence of soybean.

Union of phylogeography and landscape genetics

Phylogeography and landscape genetics have arisen within the past 30 y. Phylogeography is said to be the bridge between population genetics and systematics, and landscape genetics the bridge between landscape ecology and population genetics. Both fields can be considered as simply the amalgamation of classic biogeography with genetics and genomics; however, they differ in the temporal, spatial, and organismal scales addressed and the methodology used. I begin by briefly summarizing the history and purview of each field and suggest that, even though landscape genetics is a younger field (coined in 2003) than phylogeography (coined in 1987), early studies by Dobzhansky on the "microgeographic races" of Linanthus parryae in the Mojave Desert of California and Drosophila pseudoobscura across the western United States presaged the fields by over 40 y. Recent advances in theory, models, and methods have allowed researchers to better synthesize ecological and evolutionary processes in their quest to answer some of the most basic questions in biology. I highlight a few of these novel studies and emphasize three major areas ripe for investigation using spatially explicit genomic-scale data: the biogeography of speciation, lineage divergence and species delimitation, and understanding adaptation through time and space. Examples of areas in need of study are highlighted, and I end by advocating a union of phylogeography and landscape genetics under the more general field: biogeography.

Edaphic history over seedling characters predicts integration and plasticity of integration across geologically variable populations of Arabidopsis thaliana

PREMISE OF THE STUDY

Studies on phenotypic plasticity and plasticity of integration have uncovered functionally linked modules of aboveground traits and seedlings of Arabidopsis thaliana, but we lack details about belowground variation in adult plants. Functional modules can be comprised of additional suites of traits that respond to environmental variation. We assessed whether shoot and root responses to nutrient environments in adult A. thaliana were predictable from seedling traits or population-specific geologic soil characteristics at the site of origin.

METHODS

We compared 17 natural accessions from across the native range of A. thaliana using 14-day-old seedlings grown on agar or sand and plants grown to maturity across nutrient treatments in sand. We measured aboveground size, reproduction, timing traits, root length, and root diameter. Edaphic characteristics were obtained from a global-scale dataset and related to field data.

KEY RESULTS

We detected significant among-population variation in root traits of seedlings and adults and in plasticity in aboveground and belowground traits of adult plants. Phenotypic integration of roots and shoots varied by population and environment. Relative integration was greater in roots than in shoots, and integration was predicted by edaphic soil history, particularly organic carbon content, whereas seedling traits did not predict later ontogenetic stages.

CONCLUSIONS

Soil environment of origin has significant effects on phenotypic plasticity in response to nutrients, and on phenotypic integration of root modules and shoot modules. Root traits varied among populations in reproductively mature individuals, indicating potential for adaptive and integrated functional responses of root systems in annuals.

Long-term balancing selection contributes to adaptation in Arabidopsis and its relatives

BACKGROUND

In contrast to positive selection, which reduces genetic variation by fixing beneficial alleles, balancing selection maintains genetic variation within a population or species and plays crucial roles in adaptation in diverse organisms. However, which genes, genome-wide, are under balancing selection and the extent to which these genes are involved in adaptation are largely unknown.

RESULTS

We performed a genome-wide scan for genes under balancing selection across two plant species, Arabidopsis thaliana and its relative Capsella rubella, which diverged about 8 million generations ago. Among hundreds of genes with shared coding-region polymorphisms, we find evidence for long-term balancing selection in five genes: AT1G35220, AT2G16570, AT4G29360, AT5G38460, and AT5G44000. These genes are involved in the response to biotic and abiotic stress and other fundamental biochemical processes. More intriguingly, for these genes, we detected significant ecological diversification between the two haplotype groups, suggesting that balancing selection has been very important for adaptation.

CONCLUSIONS

Our results indicate that beyond the well-known S-locus genes and resistance genes, many loci are under balancing selection. These genes are mostly correlated with resistance to stress or other fundamental functions and likely play a more important role in adaptation to diverse habitats than previously thought.

Towards the Genomic Basis of Local Adaptation in Landraces

Landraces are key elements of agricultural biodiversity that have long been considered a source of useful traits. Their importance goes beyond subsistence agriculture and the essential need to preserve genetic diversity, because landraces are farmer-developed populations that are often adapted to environmental conditions of significance to tackle environmental concerns. It is therefore increasingly important to identify adaptive traits in crop landraces and understand their molecular basis. This knowledge is potentially useful for promoting more sustainable agricultural techniques, reducing the environmental impact of high-input cropping systems, and diminishing the vulnerability of agriculture to global climate change. In this review, we present an overview of the opportunities and limitations offered by landraces’ genomics. We discuss how rapid advances in DNA sequencing techniques, plant phenotyping, and recombinant DNA-based biotechnology encourage both the identification and the validation of the genomic signature of local adaptation in crop landraces. The integration of ‘omics’ sciences, molecular population genetics, and field studies can provide information inaccessible with earlier technological tools. Although empirical knowledge on the genetic and genomic basis of local adaptation is still fragmented, it is predicted that genomic scans for adaptation will unlock an intraspecific molecular diversity that may be different from that of modern varieties.

Genomic models predict successful coral adaptation if future ocean warming rates are reduced

Population genomic surveys suggest that climate-associated genetic variation occurs widely across species, but whether it is sufficient to allow population persistence via evolutionary adaptation has seldom been quantified. To ask whether rapid adaptation in reef-building corals can keep pace with future ocean warming, we measured genetic variation at predicted warm-adapted loci and simulated future evolution and persistence in a high-latitude population of corals from Rarotonga, Cook Islands. Alleles associated with thermal tolerance were present but at low frequencies in this cooler, southerly locality. Simulations based on predicted ocean warming in Rarotonga showed rapid evolution of heat tolerance resulting in population persistence under mild warming scenarios consistent with low CO₂ emission plans, RCP2.6 and RCP4.5. Under more severe scenarios, RCP6.0 and RCP8.5, adaptation was not rapid enough to prevent extinction. Population adaptation was faster for models based on smaller numbers of additive loci that determine thermal tolerance and for higher population growth rates. Finally, accelerated migration via transplantation of thermally tolerant individuals (1 to 5%/year) sped adaptation. These results show that cool-water corals can adapt to warmer oceans but only under mild scenarios resulting from international emissions controls. Incorporation of genomic data into models of species response to climate change offers a promising method for estimating future adaptive processes.

Detecting gene subnetworks under selection in biological pathways

Advances in high throughput sequencing technologies have created a gap between data production and functional data analysis. Indeed, phenotypes result from interactions between numerous genes, but traditional methods treat loci independently, missing important knowledge brought by network-level emerging properties. Therefore, detecting selection acting on multiple genes affecting the evolution of complex traits remains challenging. In this context, gene network analysis provides a powerful framework to study the evolution of adaptive traits and facilitates the interpretation of genome-wide data. We developed a method to analyse gene networks that is suitable to evidence polygenic selection. The general idea is to search biological pathways for subnetworks of genes that directly interact with each other and that present unusual evolutionary features. Subnetwork search is a typical combinatorial optimization problem that we solve using a simulated annealing approach. We have applied our methodology to find signals of adaptation to high-altitude in human populations. We show that this adaptation has a clear polygenic basis and is influenced by many genetic components. Our approach, implemented in the R package signet, improves on gene-level classical tests for selection by identifying both new candidate genes and new biological processes involved in adaptation to altitude.

A genome-wide association study of 23 agronomic traits in Chinese wheat landraces

Uncovering the genetic basis of agronomic traits in wheat landraces is important for ensuring global food security via the development of improved varieties. Here, 723 wheat landraces from 10 Chinese agro-ecological zones were evaluated for 23 agronomic traits in six environments. All accessions could be clustered into five subgroups based on phenotypic data via discriminant function analysis, which was highly consistent with genotypic classification. A genome-wide association study was conducted for these traits using 52 303 DArT-seq markers to identify marker-trait associations and candidate genes. Using both the general linear model and the mixed linear model, 149 significant markers were identified for 21 agronomic traits based on best linear unbiased prediction values. Considering the linkage disequilibrium decay distance in this study, significant markers within 10 cM were combined as a quantitative trait locus (QTL), with a total of 29 QTL identified for 15 traits. Of these, five QTL for heading date, flag leaf width, peduncle length, and thousand kernel weight had been reported previously. Twenty-five candidate genes associated with significant markers were identified. These included the known vernalization genes VRN-B1 and vrn-B3 and the photoperiod response genes Ppd and PRR. Overall, this study should be helpful in elucidating the underlying genetic mechanisms of complex agronomic traits and performing marker-assisted selection in wheat.

The developmental transcriptome landscape of receptive endometrium during embryo implantation in dairy goats

Under natural conditions, some embryos cannot implant successfully because of the dysfunction of receptive endometrium (RE). Thus, it is imperative for us to study the molecular mechanisms involved in the formation of the RE from pre-receptive endometrium (PE). In this study, the endometrium from gestational day 5 (D5, PE) and gestational day 15 (D15, RE) dairy goats were selected to systematically analyze the transcriptome using strand-specific Ribo-Zero RNA-Seq, >120 million high-quality paired-end reads were generated and 47,616 transcripts were identified in the endometrium of dairy goats. A total of 810 mRNAs were differentially expressed genes (DEGs) between the RE and PE meeting the criteria of P-values<0.05. Bioinformatics analysis of the DEGs revealed that a number of biological processes and pathways were potentially involved in the establishment of the RE, notably energy metabolism and amino acid metabolism. Furthermore, we speculated that CXCL14, IGFBP3, and LGALS15 potentially participated in the development of endometrium. What's more, putative SNPs, InDels and AS events were identified and analyzed in the endometrium. In a word, this resulting view of the transcriptome greatly enhances the comprehensive transcript catalog and uncovers the global trends in gene expression during the formation of receptive endometrium in dairy goats.