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

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

Climate change challenges, plant science solutions.

Climate change is a defining challenge of the 21st century, and this decade is a critical time for action to mitigate the worst effects on human populations and ecosystems. Plant science can play an important role in developing crops with enhanced resilience to harsh conditions (e.g. heat, drought, salt stress, flooding, disease outbreaks) and engineering efficient carbon-capturing and carbon-sequestering plants. Here, we present examples of research being conducted in these areas and discuss challenges and open questions as a call to action for the plant science community.

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

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

The population genomics of adaptive loss of function.

Discoveries of adaptive gene knockouts and widespread losses of complete genes have in recent years led to a major rethink of the early view that loss-of-function alleles are almost always deleterious. Today, surveys of population genomic diversity are revealing extensive loss-of-function and gene content variation, yet the adaptive significance of much of this variation remains unknown. Here we examine the evolutionary dynamics of adaptive loss of function through the lens of population genomics and consider the challenges and opportunities of studying adaptive loss-of-function alleles using population genetics models. We discuss how the theoretically expected existence of allelic heterogeneity, defined as multiple functionally analogous mutations at the same locus, has proven consistent with empirical evidence and why this impedes both the detection of selection and causal relationships with phenotypes. We then review technical progress towards new functionally explicit population genomic tools and genotype-phenotype methods to overcome these limitations. More broadly, we discuss how the challenges of studying adaptive loss of function highlight the value of classifying genomic variation in a way consistent with the functional concept of an allele from classical population genetics.

Combining population genomics and fitness QTLs to identify the genetics of local adaptation in Arabidopsis thaliana.

Evidence for adaptation to different climates in the model species Arabidopsis thaliana is seen in reciprocal transplant experiments, but the genetic basis of this adaptation remains poorly understood. Field-based quantitative trait locus (QTL) studies provide direct but low-resolution evidence for the genetic basis of local adaptation. Using high-resolution population genomic approaches, we examine local adaptation along previously identified genetic trade-off (GT) and conditionally neutral (CN) QTLs for fitness between locally adapted Italian and Swedish A. thaliana populations [Ågren J, et al. (2013) Proc Natl Acad Sci USA 110:21077-21082]. We find that genomic regions enriched in high FST SNPs colocalize with GT QTL peaks. Many of these high FST regions also colocalize with regions enriched for SNPs significantly correlated to climate in Eurasia and evidence of recent selective sweeps in Sweden. Examining unfolded site frequency spectra across genes containing high FST SNPs suggests GTs may be due to more recent adaptation in Sweden than Italy. Finally, we collapse a list of thousands of genes spanning GT QTLs to 42 genes that likely underlie the observed GTs and explore potential biological processes driving these trade-offs, from protein phosphorylation, to seed dormancy and longevity. Our analyses link population genomic analyses and field-based QTL studies of local adaptation, and emphasize that GTs play an important role in the process of local adaptation.

Mutation Accumulation in an Asexual Relative of Arabidopsis.

Asexual populations experience weaker responses to natural selection, which causes deleterious mutations to accumulate over time. Additionally, stochastic loss of individuals free of deleterious mutations can lead to an irreversible increase in mutational load in asexuals (the "click" in Muller's Ratchet). Here we report on the genomic divergence and distribution of mutations across eight sympatric pairs of sexual and apomictic (asexual) Boechera (Brassicaceae) genotypes. We show that apomicts harbor a greater number of derived mutations than sympatric sexual genotypes. Furthermore, in phylogenetically constrained sites that are subject to contemporary purifying selection, the ancestral, conserved allele is more likely to be retained in sexuals than apomicts. These results indicate that apomictic lineages accumulate mutations at otherwise conserved sites more often than sexuals, and support the conclusion that deleterious mutation accumulation can be a powerful force in the evolution of asexual higher plants.

Drought regimens predict life history strategies in Heliophila.

Explaining variation in life history strategies is an enduring goal of evolutionary biology and ecology. Early theory predicted that for plants, annual and perennial life histories reflect adaptations to environments that experience alternative drought regimens. Nevertheless, empirical support for this hypothesis from comparative analyses remains lacking. Here, we test classic life history theory in Heliophila L. (Brassicaceae), a diverse genus of flowering plants native to Africa, controlling for phylogeny and integrating 34 yr of satellite-based drought detection with 2192 herbaria occurrence records. We find that the common ancestor of these species was likely to be an annual, and that perenniality and annuality have repeatedly evolved, an estimated seven and five times, respectively. By comparing historical drought regimens, we show that annuals rather than perennial species occur in environments where droughts are significantly more frequent. We also find evidence that annual plants adapt to predictable drought regimens by escaping drought-prone seasons as seeds. These results yield compelling support for longstanding theoretical predictions by revealing the importance of drought frequency and predictability to explain plant life history. More broadly, this work highlights scalable approaches, integrating herbaria records and remote sensing to address outstanding questions in evolutionary ecology.

Deployment of Lidar from a Ground Platform: Customizing a Low-Cost, Information-Rich and User-Friendly Application for Field Phenomics Research.

Using sensors and electronic systems for characterization of plant traits provides valuable digital inputs to support complex analytical modeling in genetics research. In field applications, frequent sensor deployment enables the study of the dynamics of these traits and their interaction with the environment. This study focused on implementing lidar (light detection and ranging) technology to generate 2D displacement data at high spatial resolution and extract plant architectural parameters, namely canopy height and cover, in a diverse population of 252 maize (Zea mays L.) genotypes. A prime objective was to develop the mechanical and electrical subcomponents for field deployment from a ground vehicle. Data reduction approaches were implemented for efficient same-day post-processing to generate by-plot statistics. The lidar system was successfully deployed six times in a span of 42 days. Lidar data accuracy was validated through independent measurements in a subset of 75 experimental units. Manual and lidar-derived canopy height measurements were compared resulting in root mean square error (RMSE) = 0.068 m and r2 = 0.81. Subsequent genome-wide association study (GWAS) analyses for quantitative trait locus (QTL) identification and comparisons of genetic correlations and heritabilities for manual and lidar-based traits showed statistically significant associations. Low-cost, field-ready lidar of computational simplicity make possible timely phenotyping of diverse populations in multiple environments.

Exploiting Differential Gene Expression and Epistasis to Discover Candidate Genes for Drought-Associated QTLs in Arabidopsis thaliana.

Soil water availability represents one of the most important selective agents for plants in nature and the single greatest abiotic determinant of agricultural productivity, yet the genetic bases of drought acclimation responses remain poorly understood. Here, we developed a systems-genetic approach to characterize quantitative trait loci (QTLs), physiological traits and genes that affect responses to soil moisture deficit in the TSUxKAS mapping population of Arabidopsis thaliana. To determine the effects of candidate genes underlying QTLs, we analyzed gene expression as a covariate within the QTL model in an effort to mechanistically link markers, RNA expression, and the phenotype. This strategy produced ranked lists of candidate genes for several drought-associated traits, including water use efficiency, growth, abscisic acid concentration (ABA), and proline concentration. As a proof of concept, we recovered known causal loci for several QTLs. For other traits, including ABA, we identified novel loci not previously associated with drought. Furthermore, we documented natural variation at two key steps in proline metabolism and demonstrated that the mitochondrial genome differentially affects genomic QTLs to influence proline accumulation. These findings demonstrate that linking genome, transcriptome, and phenotype data holds great promise to extend the utility of genetic mapping, even when QTL effects are modest or complex.

Genetic Costs of Domestication and Improvement.

The "cost of domestication" hypothesis posits that the process of domesticating wild species can result in an increase in the number, frequency, and/or proportion of deleterious genetic variants that are fixed or segregating in the genomes of domesticated species. This cost may limit the efficacy of selection and thus reduce genetic gains in breeding programs for these species. Understanding when and how deleterious mutations accumulate can also provide insight into fundamental questions about the interplay of demography and selection. Here we describe the evolutionary processes that may contribute to deleterious variation accrued during domestication and improvement, and review the available evidence for "the cost of domestication" in animal and plant genomes. We identify gaps and explore opportunities in this emerging field, and finally offer suggestions for researchers and breeders interested in understanding or avoiding the consequences of an increased number or frequency of deleterious variants in domesticated species.

Screening for Natural Variation in Water Use Efficiency Traits in a Diversity Set of Brassica napus L. Identifies Candidate Variants in Photosynthetic Assimilation.

Seed yield and quality of crop species are significantly reduced by water deficit. Stable isotope screening (δ13C) of a diversity set of 147 accessions of Brassica napus grown in the field identified several accessions with extremes in water use efficiency (WUE). We next conducted an investigation of the physiological characteristics of selected natural variants with high and low WUE to understand how these characteristics translate to differences in WUE. We identified an interesting Spring accession, G302 (Mozart), that exhibited the highest WUE in the field and high CO2 assimilation rates coupled with an increased electron transport capacity (Jmax) under the imposed conditions. Differences in stomatal density and stomatal index did not translate to differences in stomatal conductance in the investigated accessions. Stomatal conductance response to exogenous ABA was analyzed in selected high and low WUE accessions. Spring lines showed little variation in response to exogenous ABA, while one Semi-Winter line (SW047) showed a significantly more rapid response to exogenous ABA, that corresponded to the high WUE indicated by δ13C measurements. This research illustrates the importance of examining natural variation at a physiological level for investigation of the underlying mechanisms influencing the diversity of carbon isotope discrimination values in the field and identifies natural variants in B. napus with improved WUE and potential relevant traits.

Variation in MPK12 affects water use efficiency in Arabidopsis and reveals a pleiotropic link between guard cell size and ABA response.

Plant water relations are critical for determining the distribution, persistence, and fitness of plant species. Studying the genetic basis of ecologically relevant traits, however, can be complicated by their complex genetic, physiological, and developmental basis and their interaction with the environment. Water use efficiency (WUE), the ratio of photosynthetic carbon assimilation to stomatal conductance to water, is a dynamic trait with tremendous ecological and agricultural importance whose genetic control is poorly understood. In the present study, we use a quantitative trait locus-mapping approach to locate, fine-map, clone, confirm, and characterize an allelic substitution that drives differences in WUE among natural accessions of Arabidopsis thaliana. We show that a single amino acid substitution in an abscisic acid-responsive kinase, AtMPK12, causes reduction in WUE, and we confirm its functional role using transgenics. We further demonstrate that natural alleles at AtMPK12 differ in their response to cellular and environmental cues, with the allele from the Cape Verde Islands (CVI) being less responsive to hormonal inhibition of stomatal opening and more responsive to short-term changes in vapor pressure deficit. We also show that the CVI allele results in constitutively larger stomata. Together, these differences cause higher stomatal conductance and lower WUE compared with the common allele. These physiological changes resulted in reduced whole-plant transpiration efficiency and reduced fitness under water-limited compared with well-watered conditions. Our work demonstrates how detailed analysis of naturally segregating functional variation can uncover the molecular and physiological basis of a key trait associated with plant performance in ecological and agricultural settings.

Molecular and systems approaches towards drought-tolerant canola crops.

1169 I. 1170 II. 1170 III. 1172 IV. 1176 V. 1181 VI. 1182 1183 References 1183 SUMMARY: Modern agriculture is facing multiple challenges including the necessity for a substantial increase in production to meet the needs of a burgeoning human population. Water shortage is a deleterious consequence of both population growth and climate change and is one of the most severe factors limiting global crop productivity. Brassica species, particularly canola varieties, are cultivated worldwide for edible oil, animal feed, and biodiesel, and suffer dramatic yield loss upon drought stress. The recent release of the Brassica napus genome supplies essential genetic information to facilitate identification of drought-related genes and provides new information for agricultural improvement in this species. Here we summarize current knowledge regarding drought responses of canola, including physiological and -omics effects of drought. We further discuss knowledge gained through translational biology based on discoveries in the closely related reference species Arabidopsis thaliana and through genetic strategies such as genome-wide association studies and analysis of natural variation. Knowledge of drought tolerance/resistance responses in canola together with research outcomes arising from new technologies and methodologies will inform novel strategies for improvement of drought tolerance and yield in this and other important crop species.

Adaptation to warmer climates by parallel functional evolution of CBF genes in Arabidopsis thaliana.

The evolutionary processes and genetics underlying local adaptation at a specieswide level are largely unknown. Recent work has indicated that a frameshift mutation in a member of a family of transcription factors, C-repeat binding factors or CBFs, underlies local adaptation and freezing tolerance divergence between two European populations of Arabidopsis thaliana. To ask whether the specieswide evolution of CBF genes in Arabidopsis is consistent with local adaptation, we surveyed CBF variation from 477 wild accessions collected across the species' range. We found that CBF sequence variation is strongly associated with winter temperature variables. Looking specifically at the minimum temperature experienced during the coldest month, we found that Arabidopsis from warmer climates exhibit a significant excess of nonsynonymous polymorphisms in CBF genes and revealed a CBF haplotype network whose structure points to multiple independent transitions to warmer climates. We also identified a number of newly described mutations of significant functional effect in CBF genes, similar to the frameshift mutation previously indicated to be locally adaptive in Italy, and find that they are significantly associated with warm winters. Lastly, we uncover relationships between climate and the position of significant functional effect mutations between and within CBF paralogs, suggesting variation in adaptive function of different mutations. Cumulatively, these findings support the hypothesis that disruption of CBF gene function is adaptive in warmer climates, and illustrate how parallel evolution in a transcription factor can underlie adaptation to climate.

Does low stomatal conductance or photosynthetic capacity enhance growth at elevated CO2 in Arabidopsis?

The objective of this study was to determine if low stomatal conductance (g) increases growth, nitrate (NO3 (-)) assimilation, and nitrogen (N) utilization at elevated CO2 concentration. Four Arabidopsis (Arabidopsis thaliana) near isogenic lines (NILs) differing in g were grown at ambient and elevated CO2 concentration under low and high NO3 (-) supply as the sole source of N. Although g varied by 32% among NILs at elevated CO2, leaf intercellular CO2 concentration varied by only 4% and genotype had no effect on shoot NO3 (-) concentration in any treatment. Low-g NILs showed the greatest CO2 growth increase under N limitation but had the lowest CO2 growth enhancement under N-sufficient conditions. NILs with the highest and lowest g had similar rates of shoot NO3 (-) assimilation following N deprivation at elevated CO2 concentration. After 5 d of N deprivation, the lowest g NIL had 27% lower maximum carboxylation rate and 23% lower photosynthetic electron transport compared with the highest g NIL. These results suggest that increased growth of low-g NILs under N limitation most likely resulted from more conservative N investment in photosynthetic biochemistry rather than from low g.

Expression quantitative trait locus mapping across water availability environments reveals contrasting associations with genomic features in Arabidopsis.

The regulation of gene expression is crucial for an organism's development and response to stress, and an understanding of the evolution of gene expression is of fundamental importance to basic and applied biology. To improve this understanding, we conducted expression quantitative trait locus (eQTL) mapping in the Tsu-1 (Tsushima, Japan) × Kas-1 (Kashmir, India) recombinant inbred line population of Arabidopsis thaliana across soil drying treatments. We then used genome resequencing data to evaluate whether genomic features (promoter polymorphism, recombination rate, gene length, and gene density) are associated with genes responding to the environment (E) or with genes with genetic variation (G) in gene expression in the form of eQTLs. We identified thousands of genes that responded to soil drying and hundreds of main-effect eQTLs. However, we identified very few statistically significant eQTLs that interacted with the soil drying treatment (GxE eQTL). Analysis of genome resequencing data revealed associations of several genomic features with G and E genes. In general, E genes had lower promoter diversity and local recombination rates. By contrast, genes with eQTLs (G) had significantly greater promoter diversity and were located in genomic regions with higher recombination. These results suggest that genomic architecture may play an important a role in the evolution of gene expression.

Genetic mapping of adaptation reveals fitness tradeoffs in Arabidopsis thaliana.

Organisms inhabiting different environments are often locally adapted, and yet despite a considerable body of theory, the genetic basis of local adaptation is poorly understood. Unanswered questions include the number and effect sizes of adaptive loci, whether locally favored loci reduce fitness elsewhere (i.e., fitness tradeoffs), and whether a lack of genetic variation limits adaptation. To address these questions, we mapped quantitative trait loci (QTL) for total fitness in 398 recombinant inbred lines derived from a cross between locally adapted populations of the highly selfing plant Arabidopsis thaliana from Sweden and Italy and grown for 3 consecutive years at the parental sites (>40,000 plants monitored). We show that local adaptation is controlled by relatively few genomic regions of small to modest effect. A third of the 15 fitness QTL we detected showed evidence of tradeoffs, which contrasts with the minimal evidence for fitness tradeoffs found in previous studies. This difference may reflect the power of our multiyear study to distinguish conditionally neutral QTL from those that reflect fitness tradeoffs. In Sweden, but not in Italy, the local genotype underlying fitness QTL was often maladaptive, suggesting that adaptation there is constrained by a lack of adaptive genetic variation, attributable perhaps to genetic bottlenecks during postglacial colonization of Scandinavia or to recent changes in selection regime caused by climate change. Our results suggest that adaptation to markedly different environments can be achieved through changes in relatively few genomic regions, that fitness tradeoffs are common, and that lack of genetic variation can limit adaptation.

Natural variation in abiotic stress responsive gene expression and local adaptation to climate in Arabidopsis thaliana.

Gene expression varies widely in natural populations, yet the proximate and ultimate causes of this variation are poorly known. Understanding how variation in gene expression affects abiotic stress tolerance, fitness, and adaptation is central to the field of evolutionary genetics. We tested the hypothesis that genes with natural genetic variation in their expression responses to abiotic stress are likely to be involved in local adaptation to climate in Arabidopsis thaliana. Specifically, we compared genes with consistent expression responses to environmental stress (expression stress responsive, "eSR") to genes with genetically variable responses to abiotic stress (expression genotype-by-environment interaction, "eGEI"). We found that on average genes that exhibited eGEI in response to drought or cold had greater polymorphism in promoter regions and stronger associations with climate than those of eSR genes or genomic controls. We also found that transcription factor binding sites known to respond to environmental stressors, especially abscisic acid responsive elements, showed significantly higher polymorphism in drought eGEI genes in comparison to eSR genes. By contrast, eSR genes tended to exhibit relatively greater pairwise haplotype sharing, lower promoter diversity, and fewer nonsynonymous polymorphisms, suggesting purifying selection or selective sweeps. Our results indicate that cis-regulatory evolution and genetic variation in stress responsive gene expression may be important mechanisms of local adaptation to climatic selective gradients.

2b-RAD: a simple and flexible method for genome-wide genotyping.

We describe 2b-RAD, a streamlined restriction site-associated DNA (RAD) genotyping method based on sequencing the uniform fragments produced by type IIB restriction endonucleases. Well-studied accessions of Arabidopsis thaliana were genotyped to validate the method's accuracy and to demonstrate fine-tuning of marker density as needed. The simplicity of the 2b-RAD protocol makes it particularly suitable for high-throughput genotyping as required for linkage mapping and profiling genetic variation in natural populations.

Combining quantitative trait loci analysis with physiological models to predict genotype-specific transpiration rates.

Transpiration is controlled by evaporative demand and stomatal conductance (gs ), and there can be substantial genetic variation in gs . A key parameter in empirical models of transpiration is minimum stomatal conductance (g0 ), a trait that can be measured and has a large effect on gs and transpiration. In Arabidopsis thaliana, g0 exhibits both environmental and genetic variation, and quantitative trait loci (QTL) have been mapped. We used this information to create a genetically parameterized empirical model to predict transpiration of genotypes. For the parental lines, this worked well. However, in a recombinant inbred population, the predictions proved less accurate. When based only upon their genotype at a single g0 QTL, genotypes were less distinct than our model predicted. Follow-up experiments indicated that both genotype by environment interaction and a polygenic inheritance complicate the application of genetic effects into physiological models. The use of ecophysiological or 'crop' models for predicting transpiration of novel genetic lines will benefit from incorporating further knowledge of the genetic control and degree of independence of core traits/parameters underlying gs variation.