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

Low frequency of the wild-type freezing-tolerance LsCBF7 allele among lettuce population suggests a negative selection during domestication and breeding

KEY MESSAGE

Sustainable winter production in lettuce requires freezing tolerant varieties. This study identified a wild-type allele of LsCBF7 that could contribute to freezing tolerance improvement in lettuce. Lettuce is one of the most consumed vegetables globally. While ideally grown in 13-21 °C, its cultivation extends into winter in milder climates. However, occasional freezing temperatures can significantly reduce yields. Therefore, the development of freezing-tolerant lettuce varieties has become a long-term goal of lettuce breeding programs. Despite its significance, our understanding of freezing tolerance in lettuce remains limited. Plants have evolved a coping mechanism against freezing, known as cold acclimation, whereby they can increase freezing tolerance when pre-exposed to low nonfreezing temperatures. The CBF pathway is well-known for its central role in cold acclimation. Previously, we identified 14 CBF genes in lettuce and discovered that one of them, LsCBF7, had a loss-of-function mutation. In this study, we uncovered that accessions from colder regions carried the wild-type allele of LsCBF7 and this allele likely contributed to increased freezing tolerance, with 14% of the lettuce population carrying this allele. Interestingly, in wild lettuce (L. serriola) that is considered a progenitor of cultivated lettuce, this wild-type allele was much more common, with a frequency of 90%. This finding suggests that this wild-type allele may have undergone negative selection during the domestication or breeding of lettuce. Our data strongly indicate that this allele could be linked to early bolting, an undesirable trait in lettuce, which may have driven the negative selection. While this wild-type allele shows promise for improving freezing tolerance in lettuce, it is crucial to decouple it from the early bolting trait to fully harness its potential in lettuce breeding.

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

AIM

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

LOCATION

Global.

TAXON

Arabidopsis thaliana ("Arabidopsis").

METHODS

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

RESULTS

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

MAIN CONCLUSIONS

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

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.

Functional Evolution of Pseudofabraea citricarpa as an Adaptation to Temperature Change

Citrus target spot, caused by Pseudofabraea citricarpa, was formerly considered a cold-tolerant fungal disease. However, it has now spread from high-latitude regions to warmer low-latitude regions. Here, we conducted physiological observations on two different strains of the fungus collected from distinct regions, and evaluated their pathogenicity. Interestingly, the CQWZ collected from a low-latitude orchard, exhibited higher temperature tolerance and pathogenicity when compared to the SXCG collected from a high-latitude orchard. To further understand the evolution of temperature tolerance and virulence in these pathogens during the spread process, as well as the mechanisms underlying these differences, we performed genomic comparative analysis. The genome size of CQWZ was determined to be 44,004,669 bp, while the genome size of SXCG was determined to be 45,377,339 bp. Through genomic collinearity analysis, we identified two breakpoints and rearrangements during the evolutionary process of these two strains. Moreover, gene annotation results revealed that the CQWZ possessed 376 annotated genes in the "Xenobiotics biodegradation and metabolism" pathway, which is 79 genes more than the SXCG. The main factor contributing to this difference was the presence of salicylate hydroxylase. We also observed variations in the oxidative stress pathways and core pathogenic genes. The CQWZ exhibited the presence of a heat shock protein (HSP SSB), a catalase (CAT2), and 13 core pathogenic genes, including a LysM effector, in comparison to the SXCG. Furthermore, there were significant disparities in the gene clusters responsible for the production of seven metabolites, such as Fumonisin and Brefeldin. Finally, we identified the regulatory relationship, with the HOG pathway at its core, that potentially contributes to the differences in thermotolerance and virulence. As the global climate continues to warm, crop pathogens are increasingly expanding to new territories. Our findings will enhance understanding of the evolution mechanisms of pathogens under climate change.

Integrative analysis of physiology and genomics provides insights into freeze tolerance adaptations of Acacia koa along an elevational cline

Natural selection for plant species in heterogeneous environments creates genetic variation for traits such as cold tolerance. While physiological or molecular analyses have been used to evaluate stress tolerance adaptations, combining these approaches may provide deeper insight. Acacia koa (koa) occurs from sea level to 2300 m in Hawai'i, USA. At high elevations, natural koa populations have declined due to deforestation, and freeze tolerance is a limiting factor for tree regeneration. We used physiology and molecular analyses to evaluate cold tolerance of koa populations from low (300-750 m), middle (750-1500 m), and high elevations (1500-2100 m). Half of the seedlings were cold acclimated by exposure to progressively lowered air temperatures for eight weeks (from 25.6/22.2°C to 8/4°C, day/night). Using the whole plant physiology-freezing test and koa C-repeat Binding Factor CBF genes, our results indicated that koa can be cold-acclimated when exposed to low, non-freezing temperatures. Seedlings from high elevations had consistently higher expression of Koa CBF genes associated with cold tolerance, helping to explain variation in cold-hardy phenotypes. Evaluation of the genetic background of 22 koa families across the elevations with low coverage RNA sequencing indicated that high elevation koa had relatively low values of heterozygosity, suggesting that adaptation is more likely to arise in the middle and low elevation sources. This physiology and molecular data for cold tolerance of koa across the elevation gradient of the Hawaiian Islands provides insights into natural selection processes and may help to support guidelines for conservation and seed transfer in forest restoration efforts.

AP2/EREBP Pathway Plays an Important Role in Chaling Wild Rice Tolerance to Cold Stress

Cold stress is the main factor limiting rice production and distribution. Chaling wild rice can survive in cold winters. AP2/EREBP is a known transcription factor family associated with abiotic stress. We identified the members of the AP2/EREBP transcription factor family in rice, maize, and Arabidopsis, and conducted collinearity analysis and gene family analysis. We used Affymetrix array technology to analyze the expression of AP2/EREBP family genes in Chaling wild rice and cultivated rice cultivar Pei'ai64S, which is sensitive to cold. According to the GeneChip results, the expression levels of AP2/EREBP genes in Chaling wild rice were different from those in Pei'ai64S; and the increase rate of 36 AP2/EREBP genes in Chaling wild rice was higher than that in Pei'ai64S. Meanwhile, the MYC elements in cultivated rice and Chaling wild rice for the Os01g49830, Os03g08470, and Os03g64260 genes had different promoter sequences, resulting in the high expression of these genes in Chaling wild rice under low-temperature conditions. Furthermore, we analyzed the upstream and downstream genes of the AP2/EREBP transcription factor family and studied the conservation of these genes. We found that the upstream transcription factors were more conserved, indicating that these upstream transcription factors may be more important in regulating cold stress. Meanwhile, we found the expression of AP2/EREBP pathway genes was significantly increased in recombinant inbred lines from Nipponbare crossing with Chaling wild rice, These results suggest that the AP2/EREBP signaling pathway plays an important role in Chaling wild rice tolerance to cold stress.

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.

Potential and limits of (mal)adaptive mutation rate plasticity in plants

Genetic mutations provide the heritable material for plant adaptation to their environments. At the same time, the environment can affect the mutation rate across plant genomes. However, the extent to which environmental plasticity in mutation rates can facilitate or hinder adaptation remains a longstanding and unresolved question. Emerging discoveries of mechanisms affecting mutation rate variability provide opportunities to consider this question in a new light. Links between chromatin states, transposable elements, and DNA repair suggest cases of adaptive mutation rate plasticity could occur. Yet, numerous evolutionary and biological forces are expected to limit the impact of any such mutation rate plasticity on adaptive evolution. Persistent uncertainty about the significance of mutation rate plasticity on adaptation motivates new experimental and theoretical research relevant to understanding plant responses in changing environments.

Genotype-environment associations to reveal the molecular basis of environmental adaptation

A fundamental goal in plant biology is to identify and understand the variation underlying plants' adaptation to their environment. Climate change has given new urgency to this goal, as society aims to accelerate adaptation of ecologically important plant species, endangered plant species, and crops to hotter, less predictable climates. In the pre-genomic era, identifying adaptive alleles was painstaking work, leveraging genetics, molecular biology, physiology, and ecology. Now, the rise of genomics and new computational approaches may facilitate this research. Genotype-environment associations (GEAs) use statistical associations between allele frequency and environment of origin to test the hypothesis that allelic variation at a given gene is adapted to local environments. Researchers may scan the genome for GEAs to generate hypotheses on adaptive genetic variants (environmental genome-wide association studies). Despite the rapid adoption of these methods, many important questions remain about the interpretation of GEA findings, which arise from fundamental unanswered questions on the genetic architecture of adaptation and limitations inherent to association-based analyses. We outline strategies to ground GEAs in the underlying hypotheses of genetic architecture and better test GEA-generated hypotheses using genetics and ecophysiology. We provide recommendations for new users who seek to learn about the molecular basis of adaptation. When combined with a rigorous hypothesis testing framework, GEAs may facilitate our understanding of the molecular basis of climate adaptation for plant improvement.

Patterns of pan-genome occupancy and gene coexpression under water-deficit in Brachypodium distachyon

Natural populations are characterized by abundant genetic diversity driven by a range of different types of mutation. The tractability of sequencing complete genomes has allowed new insights into the variable composition of genomes, summarized as a species pan-genome. These analyses demonstrate that many genes are absent from the first reference genomes, whose analysis dominated the initial years of the genomic era. Our field now turns towards understanding the functional consequence of these highly variable genomes. Here, we analysed weighted gene coexpression networks from leaf transcriptome data for drought response in the purple false brome Brachypodium distachyon and the differential expression of genes putatively involved in adaptation to this stressor. We specifically asked whether genes with variable "occupancy" in the pan-genome - genes which are either present in all studied genotypes or missing in some genotypes - show different distributions among coexpression modules. Coexpression analysis united genes expressed in drought-stressed plants into nine modules covering 72 hub genes (87 hub isoforms), and genes expressed under controlled water conditions into 13 modules, covering 190 hub genes (251 hub isoforms). We find that low occupancy pan-genes are under-represented among several modules, while other modules are over-enriched for low-occupancy pan-genes. We also provide new insight into the regulation of drought response in B. distachyon, specifically identifying one module with an apparent role in primary metabolism that is strongly responsive to drought. Our work shows the power of integrating pan-genomic analysis with transcriptomic data using factorial experiments to understand the functional genomics of environmental response.

Distinct Cold Acclimation of Productivity Traits in Arabidopsis thaliana Ecotypes

Improvement of crop climate resilience will require an understanding of whole-plant adaptation to specific local environments. This review places features of plant form and function related to photosynthetic productivity, as well as associated gene-expression patterns, into the context of the adaptation of Arabidopsis thaliana ecotypes to local environments with different climates in Sweden and Italy. The growth of plants under common cool conditions resulted in a proportionally greater emphasis on the maintenance of photosynthetic activity in the Swedish ecotype. This is compared to a greater emphasis on downregulation of light-harvesting antenna size and upregulation of a host of antioxidant enzymes in the Italian ecotype under these conditions. This differential response is discussed in the context of the climatic patterns of the ecotypes’ native habitats with substantial opportunity for photosynthetic productivity under mild temperatures in Italy but not in Sweden. The Swedish ecotype’s response is likened to pushing forward at full speed with productivity under low temperature versus the Italian ecotype’s response of staying safe from harm (maintaining redox homeostasis) while letting productivity decline when temperatures are transiently cold. It is concluded that either strategy can offer directions for the development of climate-resilient crops for specific locations of cultivation.

Genotype-dependent contribution of CBF transcription factors to long-term acclimation to high light and cool temperature

When grown under cool temperature, winter annuals upregulate photosynthetic capacity as well as freezing tolerance. Here, the role of three cold-induced C-repeat-binding factor (CBF1-3) transcription factors in photosynthetic upregulation and freezing tolerance was examined in two Arabidopsis thaliana ecotypes originating from Italy (IT) or Sweden (SW), and their corresponding CBF1-3-deficient mutant lines it:cbf123 and sw:cbf123. Photosynthetic, morphological and freezing-tolerance phenotypes, as well as gene expression profiles, were characterized in plants grown from the seedling stage under different combinations of light level and temperature. Under high light and cool (HLC) growth temperature, a greater role of CBF1-3 in IT versus SW was evident from both phenotypic and transcriptomic data, especially with respect to photosynthetic upregulation and freezing tolerance of whole plants. Overall, features of SW were consistent with a different approach to HLC acclimation than seen in IT, and an ability of SW to reach the new homeostasis through the involvement of transcriptional controls other than CBF1-3. These results provide tools and direction for further mechanistic analysis of the transcriptional control of approaches to cold acclimation suitable for either persistence through brief cold spells or for maximisation of productivity in environments with continuous low temperatures.

Temperature-induced dynamics of plant carbohydrate metabolism

Carbohydrates are direct products of photosynthetic CO₂ assimilation. Within a changing temperature regime, both photosynthesis and carbohydrate metabolism need tight regulation to prevent irreversible damage of plant tissue and to sustain energy metabolism, growth and development. Due to climate change, plants are and will be exposed to both long-term and short-term temperature changes with increasing amplitude. Particularly sudden fluctuations, which might comprise a large temperature amplitude from low to high temperature, pose a challenge for plants from the cellular to the ecosystem level. A detailed understanding of fundamental regulatory processes, which link photosynthesis and carbohydrate metabolism under such fluctuating environmental conditions, is essential for an estimate of climate change consequences. Further, understanding these processes is important for biotechnological application, breeding and engineering. Environmental light and temperature regimes are sensed by a molecular network that comprises photoreceptors and molecular components of the circadian clock. Photosynthetic efficiency and plant productivity then critically depend on enzymatic regulation and regulatory circuits connecting plant cells with their environment and re-stabilising photosynthetic efficiency and carbohydrate metabolism after temperature-induced deflection. This review summarises and integrates current knowledge about re-stabilisation of photosynthesis and carbohydrate metabolism after perturbation by changing temperature (heat and cold).

Nothing in Evolution Makes Sense Except in the Light of Biology

A key question in biology is the predictability of the evolutionary process. If we can correctly predict the outcome of evolution, we may be better equipped to anticipate and manage species' adaptation to climate change, habitat loss, invasive species, or emerging infectious diseases, as well as improve our basic understanding of the history of life on Earth. In the present article, we ask the questions when, why, and if the outcome of future evolution is predictable. We first define predictable and then discuss two conflicting views: that evolution is inherently unpredictable and that evolution is predictable given the ability to collect the right data. We identify factors that generate unpredictability, the data that might be required to make predictions at some level of precision or at a specific timescale, and the intellectual and translational value of understanding when prediction is or is not possible.

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.

Competition × drought interactions change phenotypic plasticity and the direction of selection on Arabidopsis traits

Populations often exhibit genetic diversity in traits involved in responses to abiotic stressors, but what maintains this diversity is unclear. Arabidopsis thaliana exhibits high within-population variation in drought response. One hypothesis is that competition, varying at small scales, promotes diversity in resource use strategies. However, little is known about natural variation in competition effects on Arabidopsis physiology. We imposed drought and competition treatments on diverse genotypes. We measured resource economics traits, physiology, and fitness to characterize plasticity and selection in response to treatments. Plastic responses to competition differed depending on moisture availability. We observed genotype-drought-competition interactions for relative fitness: competition had little effect on relative fitness under well-watered conditions, whereas competition caused rank changes in fitness under drought. Early flowering was always selected. Higher δ¹³ C was selected only in the harshest treatment (drought and competition). Competitive context significantly changed the direction of selection on aboveground biomass and inflorescence height in well-watered environments. Our results highlight how local biotic conditions modify abiotic selection, in some cases promoting diversity in abiotic stress response. The ability of populations to adapt to environmental change may thus depend on small-scale biotic heterogeneity.

Indications for a Central Role of Hexokinase Activity in Natural Variation of Heat Acclimation in Arabidopsis thaliana

Diurnal and seasonal changes of abiotic environmental factors shape plant performance and distribution. Changes of growth temperature and light intensity may vary significantly on a diurnal, but also on a weekly or seasonal scale. Hence, acclimation to a changing temperature and light regime is essential for plant survival and propagation. In the present study, we analyzed photosynthetic CO₂ assimilation and metabolic regulation of the central carbohydrate metabolism in two natural accessions of Arabidopsis thaliana that originate from north western Russia and south Italy during exposure to heat and a combination of heat and high light. Our findings indicate that it is hardly possible to predict photosynthetic capacities under combined stress from single stress experiments. Further, capacities of hexose phosphorylation were found to be significantly lower in the Italian than in the Russian accession, which could explain an inverted sucrose-to-hexose ratio. Together with the finding of significantly stronger accumulation of anthocyanins under heat/high light, these observations indicate a central role of hexokinase activity in the stabilization of photosynthesis and carbohydrate metabolism during environmental changes.

The cold response regulator CBF1 promotes Arabidopsis hypocotyl growth at ambient temperatures

Light and temperature are two core environmental factors that coordinately regulate plant growth and survival throughout their entire life cycle. However, the mechanisms integrating light and temperature signaling pathways in plants remain poorly understood. Here, we report that CBF1, an AP2/ERF-family transcription factor essential for plant cold acclimation, promotes hypocotyl growth under ambient temperatures in Arabidopsis. We show that CBF1 increases the protein abundance of PIF4 and PIF5, two phytochrome-interacting bHLH-family transcription factors that play pivotal roles in modulating plant growth and development, by directly binding to their promoters to induce their gene expression, and by inhibiting their interaction with phyB in the light. Moreover, our data demonstrate that CBF1 promotes PIF4/PIF5 protein accumulation and hypocotyl growth at both 22°C and 17°C, but not at 4°C, with a more prominent role at 17°C than at 22°C. Together, our study reveals that CBF1 integrates light and temperature control of hypocotyl growth by promoting PIF4 and PIF5 protein abundance in the light, thus providing insights into the integration mechanisms of light and temperature signaling pathways in plants.

Adaptation to climate change through genetic accommodation and assimilation of plastic phenotypes

Theory suggests that evolutionary changes in phenotypic plasticity could either hinder or facilitate evolutionary rescue in a changing climate. Nevertheless, the actual role of evolving plasticity in the responses of natural populations to climate change remains unresolved. Direct observations of evolutionary change in nature are rare, making it difficult to assess the relative contributions of changes in trait means versus changes in plasticity to climate change responses. To address this gap, this review explores several proxies that can be used to understand evolving plasticity in the context of climate change, including space for time substitutions, experimental evolution and tests for genomic divergence at environmentally responsive loci. Comparisons among populations indicate a prominent role for divergence in environmentally responsive traits in local adaptation to climatic gradients. Moreover, genomic comparisons among such populations have identified pervasive divergence in the regulatory regions of environmentally responsive loci. Taken together, these lines of evidence suggest that divergence in plasticity plays a prominent role in adaptation to climatic gradients over space, indicating that evolving plasticity is also likely to play a key role in adaptive responses to climate change through time. This suggests that genetic variation in plastic responses to the environment (G × E) might be an important predictor of species' vulnerabilities to climate-driven decline or extinction. This article is part of the theme issue 'The role of plasticity in phenotypic adaptation to rapid environmental change'.

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

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

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

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

Genetic and physiological mechanisms of freezing tolerance in locally adapted populations of a winter annual

PREMISE

Despite myriad examples of local adaptation, the phenotypes and genetic variants underlying such adaptive differentiation are seldom known. Recent work on freezing tolerance and local adaptation in ecotypes of Arabidopsis thaliana from Italy and Sweden provides an essential foundation for uncovering the genotype-phenotype-fitness map for an adaptive response to a key environmental stress.

METHODS

We examined the consequences of a naturally occurring loss-of-function (LOF) mutation in an Italian allele of the gene that encodes the transcription factor CBF2, which underlies a major freezing-tolerance locus. We used four lines with a Swedish genetic background, each containing a LOF CBF2 allele. Two lines had introgression segments containing the Italian CBF2 allele, and two contained deletions created using CRISPR-Cas9. We used a growth chamber experiment to quantify freezing tolerance and gene expression before and after cold acclimation.

RESULTS

Freezing tolerance was lower in the Italian (11%) compared to the Swedish (72%) ecotype, and all four experimental CBF2 LOF lines had reduced freezing tolerance compared to the Swedish ecotype. Differential expression analyses identified 10 genes for which all CBF2 LOF lines, and the IT ecotype had similar patterns of reduced cold responsive expression compared to the SW ecotype.

CONCLUSIONS

We identified 10 genes that are at least partially regulated by CBF2 that may contribute to the differences in cold-acclimated freezing tolerance between the Italian and Swedish ecotypes. These results provide novel insight into the molecular and physiological mechanisms connecting a naturally occurring sequence polymorphism to an adaptive response to freezing conditions.

Drought adaptation in Arabidopsis thaliana by extensive genetic loss-of-function

Interdisciplinary syntheses are needed to scale up discovery of the environmental drivers and molecular basis of adaptation in nature. Here we integrated novel approaches using whole genome sequences, satellite remote sensing, and transgenic experiments to study natural loss-of-function alleles associated with drought histories in wild Arabidopsis thaliana. The genes we identified exhibit population genetic signatures of parallel molecular evolution, selection for loss-of-function, and shared associations with flowering time phenotypes in directions consistent with longstanding adaptive hypotheses seven times more often than expected by chance. We then confirmed predicted phenotypes experimentally in transgenic knockout lines. These findings reveal the importance of drought timing to explain the evolution of alternative drought tolerance strategies and further challenge popular assumptions about the adaptive value of genetic loss-of-function in nature. These results also motivate improved species-wide sequencing efforts to better identify loss-of-function variants and inspire new opportunities for engineering climate resilience in crops.

Water shortages caused by droughts lead to crop losses that affect billions of people around the world each year. By discovering how wild plants adapt to drought, it may be possible to identify traits and genes that help to improve the growth of crop plants when water is scarce. It has been suggested that plants have adapted to droughts by flowering at times of the year when droughts are less likely to occur. For example, if droughts are more likely to happen in spring, the plants may delay flowering until the summer.

Arabidopsis thaliana is a small plant that is found across Eurasia, Africa and North America, including in areas that are prone to drought at different times of the year. Individual plants of the same species may carry different versions of the same gene (known as alleles). Some of these alleles may not work properly and are referred to as loss-of-function alleles. Monroe et al. investigated whether A. thaliana plants carry any loss-of-function alleles that are associated with droughts happening in the spring or summer, and whether they are linked to when those plants will flower.

Monroe et al. analyzed satellite images collected over the last 30 years to measure when droughts have occurred. Next, they searched genome sequences of Arabidopsis thaliana for alleles that might help the plants to adapt to droughts in the spring or summer. Combining the two approaches revealed that loss-of-function alleles associated with spring droughts were strongly predicted to be associated with the plants flowering later in the year. Similarly, loss-of-function alleles associated with summer droughts were predicted to be associated with the plants flowering earlier in the year.

These findings support the idea that plants can adapt to drought by changing when they produce flowers, and suggest that loss-of-function alleles play a major role in this process. New techniques for editing genes mean it is easier than ever to generate new loss-of-function alleles in specific genes. Therefore, the results presented by Monroe et al. may help researchers to develop new varieties of crop plants that are better adapted to droughts.

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.

Optimization of Photosynthetic Productivity in Contrasting Environments by Regulons Controlling Plant Form and Function

We review the role of a family of transcription factors and their regulons in maintaining high photosynthetic performance across a range of challenging environments with a focus on extreme temperatures and water availability. Specifically, these transcription factors include CBFs (C-repeat binding factors) and DREBs (dehydration-responsive element-binding), with CBF/DREB1 primarily orchestrating cold adaptation and other DREBs serving in heat, drought, and salinity adaptation. The central role of these modulators in plant performance under challenging environments is based on (i) interweaving of these regulators with other key signaling networks (plant hormones and redox signals) as well as (ii) their function in integrating responses across the whole plant, from light-harvesting and sugar-production in the leaf to foliar sugar export and water import and on to the plant's sugar-consuming sinks (growth, storage, and reproduction). The example of Arabidopsisthaliana ecotypes from geographic origins with contrasting climates is used to describe the links between natural genetic variation in CBF transcription factors and the differential acclimation of plant anatomical and functional features needed to support superior photosynthetic performance in contrasting environments. Emphasis is placed on considering different temperature environments (hot versus cold) and light environments (limiting versus high light), on trade-offs between adaptations to contrasting environments, and on plant lines minimizing such trade-offs.

Identification of genes from the ICE-CBF-COR pathway under cold stress in Aegilops-Triticum composite group and the evolution analysis with those from Triticeae

Adverse environmental conditions limit various aspects of plant growth, productivity, and ecological distribution. To get more insights into the signaling pathways under low temperature, we identified 10 C-repeat binding factors (CBFs), 9 inducer of CBF expression (ICEs) and 10 cold-responsive (CORs) genes from Aegilops-Triticum composite group under cold stress. Conserved amino acids analysis revealed that all CBF, ICE, COR contained specific and typical functional domains. Phylogenetic analysis of CBF proteins from Triticeae showed that these CBF homologs were divided into 11 groups. CBFs from Triticum were found in every group, which shows that these CBFs generated prior to the divergence of the subfamilies of Triticeae. The evolutionary relationship among the ICE and COR proteins in Poaceae were divided into four groups with high multispecies specificity, respectively. Moreover, expression analysis revealed that mRNA accumulation was altered by cold treatment and the genes of three types involved in the ICE-CBF-COR signaling pathway were induced by cold stress. Together, the results make CBF, ICE, COR genes family in Triticeae more abundant, and provide a starting point for future studies on transcriptional regulatory network for improvement of chilling tolerance in crop.

Gene Regulatory Networks Mediating Cold Acclimation: The CBF Pathway

Under low nonfreezing temperature conditions, plants from temperate climates undergo physiological and biochemical adjustments that increase their tolerance to freezing temperatures. This response, termed cold acclimation, is largely regulated by changes in gene expression. Molecular and genetic studies have identified a small family of transcription factors, called C-repeat binding factors (CBFs), as key regulators of the transcriptomic rearrangement that leads to cold acclimation. The function of these proteins is tightly controlled, and an inadequate supply of CBF activity may be detrimental to the plant. Accumulated evidence has revealed an extremely intricate network of positive and negative regulators of cold acclimation that coalesce at the level of CBF promoters constituting a central hub where multiple internal and external signals are integrated. Moreover, CBF expression is also controlled at posttranscriptional and posttranslational levels further refining CBF regulation. Recently, natural variation studies in Arabidopsis have demonstrated that mutations resulting in changes in CBF expression have an adaptive value for wild populations. Intriguingly, CBF genes are also present in plant species that do not cold acclimate, which suggest that they may also have additional functions. For instance, CBFs are required for some cold-related abiotic stress responses. In addition, their involvement in plant development deserves further study. Although more studies are necessary to fully harness CBF biotechnological potential, these transcription factors are meant to be key for a rational design of crops with enhanced tolerance to abiotic stress.

Transgenerational effects of mild heat in Arabidopsis thaliana show strong genotype specificity that is explained by climate at origin

Transgenerational environmental effects can trigger strong phenotypic variation. However, it is unclear how cues from different preceding generations interact. Also, little is known about the genetic variation for these life history traits. Here, we present the effects of grandparental and parental mild heat, and their combination, on four traits of the third-generation phenotype of 14 Arabidopsis thaliana genotypes. We tested for correlations of these effects with climate and constructed a conceptual model to identify the environmental conditions that favour the parental effect on flowering time. We observed strong evidence for genotype-specific transgenerational effects. On average, A. thaliana accustomed to mild heat produced more seeds after two generations. Parental effects overruled grandparental effects in all traits except reproductive biomass. Flowering was generally accelerated by all transgenerational effects. Notably, the parental effect triggered earliest flowering in genotypes adapted to dry summers. Accordingly, this parental effect was favoured in the model when early summer heat terminated the growing season and environments were correlated across generations. Our results suggest that A. thaliana can partly accustom to mild heat over two generations and genotype-specific parental effects show non-random evolutionary divergence across populations that may support climate change adaptation in the Mediterranean.

Genome-wide analysis of UDP-glycosyltransferase super family in Brassica rapa and Brassica oleracea reveals its evolutionary history and functional characterization

Background

Glycosyltransferases comprise a highly divergent and polyphyletic multigene family that is involved in widespread modification of plant secondary metabolites in a process called glycosylation. According to conserved domains identified in their amino acid sequences, these glycosyltransferases can be classified into a single UDP-glycosyltransferase (UGT) 1 superfamily.

Results

We performed genome-wide comparative analysis of UGT genes to trace evolutionary history in algae, bryophytes, pteridophytes, and angiosperms; then, we further investigated the expansion mechanisms and function characterization of UGT gene families in Brassica rapa and Brassica oleracea. Using Hidden Markov Model search, we identified 3, 21, 140, 200, 115, 147, and 147 UGTs in Chlamydomonas reinhardtii, Physcomitrella patens, Selaginella moellendorffii, Oryza sativa, Arabidopsis thaliana, B. rapa, and B. oleracea, respectively. Phylogenetic analysis revealed that UGT80 gene family is an ancient gene family, which is shared by all plants and UGT74 gene family is shared by ferns and angiosperms, but the remaining UGT gene families were shared by angiosperms. In dicot lineage, UGTs among three species were classified into three subgroups containing 3, 6, and 12 UGT gene families. Analysis of chromosomal distribution indicates that 98.6 and 71.4% of UGTs were located on B. rapa and B. oleracea pseudo-molecules, respectively. Expansion mechanism analyses uncovered that whole genome duplication event exerted larger influence than tandem duplication on expansion of UGT gene families in B. rapa, and B. oleracea. Analysis of selection forces of UGT orthologous gene pairs in B. rapa, and B. oleracea compared to A. thaliana suggested that orthologous genes in B. rapa, and B. oleracea have undergone negative selection, but there were no significant differences between A. thaliana –B. rapa and A. thaliana –B. oleracea lineages. Our comparisons of expression profiling illustrated that UGTs in B. rapa performed more discrete expression patterns than these in B. oleracea indicating stronger function divergence. Combing with phylogeny and expression analysis, the UGTs in B. rapa and B. oleracea experienced parallel evolution after they diverged from a common ancestor.

Conclusion

We first traced the evolutionary history of UGT gene families in plants and revealed its evolutionary and functional characterization of UGTs in B. rapa, and B. oleracea. This study provides novel insights into the evolutionary history and functional divergence of important traits or phenotype-related gene families in plants.

Electronic supplementary material

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

Topological features of a gene co-expression network predict patterns of natural diversity in environmental response

Molecular interactions affect the evolution of complex traits. For instance, adaptation may be constrained by pleiotropic or epistatic effects, both of which can be reflected in the structure of molecular interaction networks. To date, empirical studies investigating the role of molecular interactions in phenotypic evolution have been idiosyncratic, offering no clear patterns. Here, we investigated the network topology of genes putatively involved in local adaptation to two abiotic stressors-drought and cold-in Arabidopsis thaliana Our findings suggest that the gene-interaction topologies for both cold and drought stress response are non-random, with genes that show genetic variation in drought expression response (eGxE) being significantly more peripheral and cold response genes being significantly more central than genes which do not show GxE. We suggest that the observed topologies reflect different constraints on the genetic pathways involved in environmental response. The approach presented here may inform predictive models linking genetic variation in molecular signalling networks with phenotypic variation, specifically traits involved in environmental response.

Habitat Temperature and Precipitation of Arabidopsis thaliana Ecotypes Determine the Response of Foliar Vasculature, Photosynthesis, and Transpiration to Growth Temperature

Acclimatory adjustments of foliar vascular architecture, photosynthetic capacity, and transpiration rate in Arabidopsis thaliana ecotypes (Italian, Polish [Col-0], Swedish) were characterized in the context of habitat of origin. Temperatures of the habitat of origin decreased linearly with increasing habitat latitude, but habitat precipitation was greatest in Italy, lowest in Poland, and intermediate in Sweden. Plants of the three ecotypes raised under three different growth temperature regimes (low, moderate, and high) exhibited highest photosynthetic capacities, greatest leaf thickness, highest chlorophyll a/b ratio and levels of β-carotene, and greatest levels of wall ingrowths in phloem transfer cells, and, in the Col-0 and Swedish ecotypes, of phloem per minor vein in plants grown at the low temperature. In contrast, vein density and minor vein tracheary to sieve element ratio increased with increasing growth temperature - most strongly in Col-0 and least strongly in the Italian ecotype - and transpirational water loss correlated with vein density and number of tracheary elements per minor vein. Plotting of these vascular features as functions of climatic conditions in the habitat of origin suggested that temperatures during the evolutionary history of the ecotypes determined acclimatory responses of the foliar phloem and photosynthesis to temperature in this winter annual that upregulates photosynthesis in response to lower temperature, whereas the precipitation experienced during the evolutionary history of the ecotypes determined adjustment of foliar vein density, xylem, and transpiration to temperature. In particular, whereas photosynthetic capacity, leaf thickness, and foliar minor vein phloem features increased linearly with increasing latitude and decreasing temperature of the habitats of origin in response to experimental growth at low temperature, transpiration rate, foliar vein density, and minor vein tracheary element numbers and cross-sectional areas increased linearly with decreasing precipitation level in the habitats of origin in response to experimental growth at high temperature. This represents a situation where temperature acclimation of the apparent capacity for water flux through the xylem and transpiration rate in a winter annual responded differently from that of photosynthetic capacity, in contrast to previous reports of strong relationships between hydraulic conductance and photosynthesis in other studies.