QTL mapping of freezing tolerance: links to fitness and adaptive trade-offs

Conditionally Deleterious Mutation Load Accumulates in Genomic Islands of Local Adaptation but Can Be Purged with Sufficient Genotypic Redundancy

AbstractLocal adaptation frequently evolves in patches or environments that are connected via migration. In these cases, genomic regions that are linked to a locally adapted locus experience reduced effective migration rates. Via individual-based simulations of a two-patch system, we show that this reduced effective migration results in the accumulation of conditionally deleterious mutations, but not universally deleterious mutations, adjacent to adaptive loci. When there is redundancy in the genetic basis of local adaptation (i.e., genotypic redundancy), turnover of locally adapted polymorphisms allows conditionally deleterious mutation load to be purged. The amount of mutational load that accumulates adjacent to locally adapted loci is dependent on redundancy, recombination rate, migration rate, population size, strength of selection, and the phenotypic effect size of adaptive alleles. Our results highlight the need to be cautious when interpreting patterns of local adaptation at the level of phenotype or fitness, as the genetic basis of local adaptation can be transient, and evolution may confer a degree of maladaptation to nonlocal environments.

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.

Strategies to identify and dissect trade-offs in plants

Trade-offs between traits arise and reflect constraints imposed by the environment and physicochemical laws. Trade-off situations are expected to be highly relevant for sessile plants, which have to respond to changes in the environment to ensure survival. Despite increasing interest in determining the genetic and molecular basis of plant trade-offs, there are still gaps and differences with respect to how trade-offs are defined, how they are measured, and how their genetic architecture is dissected. The first step to fill these gaps is to establish what is meant by trade-offs. In this review we provide a classification of the existing definitions of trade-offs according to: (1) the measures used for their quantification, (2) the dependence of trade-offs on environment, and (3) experimental designed used (i.e. a single individual across different environments or a population of individuals in single or multiple environments). We then compare the approaches for quantification of trade-offs based on phenotypic, between-individual, and genetic correlations, and stress the need for developing further quantification indices particularly for trade-offs between multiple traits. Lastly, we highlight the genetic mechanisms underpinning trade-offs and experimental designs that facilitate their discovery in plants, with focus on usage of natural variability. This review also offers a perspective for future research aimed at identification of plant trade-offs, dissection of their genetic architecture, and development of strategies to overcome trade-offs, with applications in crop breeding.

Predicting gene expression responses to environment in Arabidopsis thaliana using natural variation in DNA sequence

The evolution of gene expression responses are a critical component of adaptation to variable environments. Predicting how DNA sequence influences expression is challenging because the genotype to phenotype map is not well resolved for cis regulatory elements, transcription factor binding, regulatory interactions, and epigenetic features, not to mention how these factors respond to environment. We tested if flexible machine learning models could learn some of the underlying cis- regulatory genotype to phenotype map. We tested this approach using cold-responsive transcriptome profiles in 5 diverse Arabidopsis thaliana accessions. We first tested for evidence that cis regulation plays a role in environmental response, finding 14 and 15 motifs that were significantly enriched within the up- and down-stream regions of cold-responsive differentially regulated genes (DEGs). We next applied convolutional neural networks (CNNs), which learn de novo cis- regulatory motifs in DNA sequences to predict expression response to environment. We found that CNNs predicted differential expression with moderate accuracy, with evidence that predictions were hindered by biological complexity of regulation and the large potential regulatory code. Overall, DEGs between specific environments can be predicted based on variation in cis- regulatory sequences, although more information needs to be incorporated and better models may be required.

Dominance reversals: the resolution of genetic conflict and maintenance of genetic variation

Beneficial reversals of dominance reduce the costs of genetic trade-offs and can enable selection to maintain genetic variation for fitness. Beneficial dominance reversals are characterized by the beneficial allele for a given context (e.g. habitat, developmental stage, trait or sex) being dominant in that context but recessive where deleterious. This context dependence at least partially mitigates the fitness consequence of heterozygotes carrying one non-beneficial allele for their context and can result in balancing selection that maintains alternative alleles. Dominance reversals are theoretically plausible and are supported by mounting empirical evidence. Here, we highlight the importance of beneficial dominance reversals as a mechanism for the mitigation of genetic conflict and review the theory and empirical evidence for them. We identify some areas in need of further research and development and outline three methods that could facilitate the identification of antagonistic genetic variation (dominance ordination, allele-specific expression and allele-specific ATAC-Seq (assay for transposase-accessible chromatin with sequencing)). There is ample scope for the development of new empirical methods as well as reanalysis of existing data through the lens of dominance reversals. A greater focus on this topic will expand our understanding of the mechanisms that resolve genetic conflict and whether they maintain genetic variation.

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.

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.

Foliar Phenotypic Plasticity Reflects Adaptation to Environmental Variability

Arabidopsis thaliana ecotypes adapted to native habitats with different daylengths, temperatures, and precipitation were grown experimentally under seven combinations of light intensity and leaf temperature to assess their acclimatory phenotypic plasticity in foliar structure and function. There were no differences among ecotypes when plants developed under moderate conditions of 400 µmol photons m-2 s-1 and 25 °C. However, in response to more extreme light or temperature regimes, ecotypes that evolved in habitats with pronounced differences in either the magnitude of changes in daylength or temperature or in precipitation level exhibited pronounced adjustments in photosynthesis and transpiration, as well as anatomical traits supporting these functions. Specifically, when grown under extremes of light intensity (100 versus 1000 µmol photons m-2 s-1) or temperature (8 °C versus 35 °C), ecotypes from sites with the greatest range of daylengths and temperature over the growing season exhibited the greatest differences in functional and structural features related to photosynthesis (light- and CO2-saturated capacity of oxygen evolution, leaf dry mass per area or thickness, phloem cells per minor vein, and water-use efficiency of CO2 uptake). On the other hand, the ecotype from the habitat with the lowest precipitation showed the greatest plasticity in features related to water transport and loss (vein density, ratio of water to sugar conduits in foliar minor veins, and transpiration rate). Despite these differences, common structure-function relationships existed across all ecotypes and growth conditions, with significant positive, linear correlations (i) between photosynthetic capacity (ranging from 10 to 110 µmol O2 m-2 s-1) and leaf dry mass per area (from 10 to 75 g m-2), leaf thickness (from 170 to 500 µm), and carbohydrate-export infrastructure (from 6 to 14 sieve elements per minor vein, from 2.5 to 8 µm2 cross-sectional area per sieve element, and from 16 to 82 µm2 cross-sectional area of sieve elements per minor vein); (ii) between transpiration rate (from 1 to 17 mmol H2O m-2 s-1) and water-transport infrastructure (from 3.5 to 8 tracheary elements per minor vein, from 13.5 to 28 µm2 cross-sectional area per tracheary element, and from 55 to 200 µm2 cross-sectional area of tracheary elements per minor vein); (iii) between the ratio of transpirational water loss to CO2 fixation (from 0.2 to 0.7 mol H2O to mmol-1 CO2) and the ratio of water to sugar conduits in minor veins (from 0.4 to 1.1 tracheary to sieve elements, from 4 to 6 µm2 cross-sectional area of tracheary to sieve elements, and from 2 to 6 µm2 cross-sectional area of tracheary elements to sieve elements per minor vein); (iv) between sugar conduits and sugar-loading cells; and (v) between water conducting and sugar conducting cells. Additionally, the proportion of water conduits to sugar conduits was greater for all ecotypes grown experimentally under warm-to-hot versus cold temperature. Thus, developmental acclimation to the growth environment included ecotype-dependent foliar structural and functional adjustments resulting in multiple common structural and functional relationships.

Homozygosity mapping identified loci and candidate genes responsible for freezing tolerance in Camelina sativa

Homozygosity mapping is an effective tool for detecting genomic regions responsible for a given trait when the phenotype is controlled by a limited number of dominant or co-dominant loci. Freezing tolerance is a major attribute in agricultural crops such as camelina. Previous studies indicated that freezing tolerance differences between a tolerant (Joelle) and susceptible (CO46) variety of camelina were controlled by a small number of dominant or co-dominant genes. We performed whole genome homozygosity mapping to identify markers and candidate genes responsible for freezing tolerance difference between these two genotypes. A total of 28 F3 RILs were sequenced to ∼30× coverage, and parental lines were sequenced to >30-40× coverage with Pacific Biosciences high fidelity technology and 60× coverage using Illumina whole genome sequencing. Overall, about 126k homozygous single nucleotide polymorphism markers were identified that differentiate both parents. Moreover, 617 markers were also homozygous in F3 families fixed for freezing tolerance/susceptibility. All these markers mapped to two contigs forming a contiguous stretch of chromosome 11. The homozygosity mapping detected 9 homozygous blocks among the selected markers and 22 candidate genes with strong similarity to regions in or near the homozygous blocks. Two such genes were differentially expressed during cold acclimation in camelina. The largest block contained a cold-regulated plant thionin and a putative rotamase cyclophilin 2 gene previously associated with freezing resistance in arabidopsis (Arabidopsis thaliana). The second largest block contains several cysteine-rich RLK genes and a cold-regulated receptor serine/threonine kinase gene. We hypothesize that one or more of these genes may be primarily responsible for freezing tolerance differences in camelina varieties.

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.

FER and LecRK show haplotype-dependent cold-responsiveness and mediate freezing tolerance in Lotus japonicus

Many plant species have succeeded in colonizing a wide range of diverse climates through local adaptation, but the underlying molecular genetics remain obscure. We previously found that winter survival was a direct target of selection during colonization of Japan by the perennial legume Lotus japonicus and identified associated candidate genes. Here, we show that two of these, FERONIA-receptor like kinase (LjFER) and a S-receptor-like kinase gene (LjLecRK), are required for non-acclimated freezing tolerance and show haplotype-dependent cold-responsive expression. Our work suggests that recruiting a conserved growth regulator gene, FER, and a receptor-like kinase gene, LecRK, into the set of cold-responsive genes has contributed to freezing tolerance and local climate adaptation in L. japonicus, offering functional genetic insight into perennial herb evolution.

Local adaptation: Causal agents of selection and adaptive trait divergence

Divergent selection across the landscape can favor the evolution of local adaptation in populations experiencing contrasting conditions. Local adaptation is widely observed in a diversity of taxa, yet we have a surprisingly limited understanding of the mechanisms that give rise to it. For instance, few have experimentally confirmed the biotic and abiotic variables that promote local adaptation, and fewer yet have identified the phenotypic targets of selection that mediate local adaptation. Here, we highlight critical gaps in our understanding of the process of local adaptation and discuss insights emerging from in-depth investigations of the agents of selection that drive local adaptation, the phenotypes they target, and the genetic basis of these phenotypes. We review historical and contemporary methods for assessing local adaptation, explore whether local adaptation manifests differently across life history, and evaluate constraints on local adaptation.

Climate drives rhizosphere microbiome variation and divergent selection between geographically distant Arabidopsis populations

Disentangling the contribution of climatic and edaphic factors to microbiome variation and local adaptation in plants requires an experimental approach to uncouple their effects and test for causality. We used microbial inocula, soil matrices and plant genotypes derived from two natural Arabidopsis thaliana populations in northern and southern Europe in an experiment conducted in climatic chambers mimicking seasonal changes in temperature, day length and light intensity of the home sites of the two genotypes. The southern A. thaliana genotype outperformed the northern genotype in the southern climate chamber, whereas the opposite was true in the northern climate chamber. Recipient soil matrix, but not microbial composition, affected plant fitness, and effects did not differ between genotypes. Differences between chambers significantly affected rhizosphere microbiome assembly, although these effects were small in comparison with the shifts induced by physicochemical differences between soil matrices. The results suggest that differences in seasonal changes in temperature, day length and light intensity between northern and southern Europe have strongly influenced adaptive differentiation between the two A. thaliana populations, whereas effects of differences in soil factors have been weak. By contrast, below-ground differences in soil characteristics were more important than differences in climate for rhizosphere microbiome differentiation.

Locally adaptive temperature response of vegetative growth in Arabidopsis thaliana

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

The influence of genetic architecture on responses to selection under drought in rice

Accurately predicting responses to selection is a major goal in biology and important for successful crop breeding in changing environments. However, evolutionary responses to selection can be constrained by such factors as genetic and cross‐environment correlations, linkage, and pleiotropy, and our understanding of the extent and impact of such constraints is still developing. Here, we conducted a field experiment to investigate potential constraints to selection for drought resistance in rice (Oryza sativa) using phenotypic selection analysis and quantitative genetics. We found that traits related to drought response were heritable, and some were under selection, including selection for earlier flowering, which could allow drought escape. However, patterns of selection generally were not opposite under wet and dry conditions, and we did not find individual or closely linked genes that influenced multiple traits, indicating a lack of evidence that antagonistic pleiotropy, linkage, or cross‐environment correlations would constrain selection for drought resistance. In most cases, genetic correlations had little influence on responses to selection, with direct and indirect selection largely congruent. The exception to this was seed mass under drought, which was predicted to evolve in the opposite direction of direct selection due to correlations. Because of this indirect effect on selection on seed mass, selection for drought resistance was not accompanied by a decrease in seed mass, and yield increased with fecundity. Furthermore, breeding lines with high fitness and yield under drought also had high fitness and yield under wet conditions, indicating that there was no evidence for a yield penalty on drought resistance. We found multiple genes in which expression influenced both water use efficiency (WUE) and days to first flowering, supporting a genetic basis for the trade‐off between drought escape and avoidance strategies. Together, these results can provide helpful guidance for understanding and managing evolutionary constraints and breeding stress‐resistant crops.

Seed dormancy varies widely among Arabidopsis thaliana populations both between and within Fennoscandia and Italy

The timing of germination is a key life‐history trait in plants, which is strongly affected by the strength of seed dormancy. Continental‐wide genetic variation in seed dormancy has been related to differences in climate and the timing of conditions suitable for seedling establishment. However, for predictions of adaptive potential and consequences of climatic change, information is needed regarding the extent to which seed dormancy varies within climatic regions and the factors driving such variation.

We quantified dormancy of seeds produced by 17 Italian and 28 Fennoscandian populations of Arabidopsis thaliana when grown in the greenhouse and at two field sites in Italy and Sweden. To identify possible drivers of among‐population variation in seed dormancy, we examined the relationship between seed dormancy and climate at the site of population origin, and between seed dormancy and flowering time.

Seed dormancy was on average stronger in the Italian compared to the Fennoscandian populations, but also varied widely within both regions. Estimates of seed dormancy in the three maternal environments were positively correlated. Among Fennoscandian populations, seed dormancy tended to increase with increasing summer temperature and decreasing precipitation at the site of population origin. In the smaller sample of Italian populations, no significant association was detected between mean seed dormancy and climate at the site of origin. The correlation between population mean seed dormancy and flowering time was weak and not statistically significant within regions.

The correlation between seed dormancy and climatic factors in Fennoscandia suggests that at least some of the among‐population variation is adaptive and that climate change will affect selection on this trait. The weak correlation between population mean seed dormancy and flowering time indicates that the two traits can evolve independently.

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.

Cold stress and freezing tolerance negatively affect the fitness of Arabidopsis thaliana accessions under field and controlled conditions

Higher acclimated freezing tolerance improved winter survival, but reduced reproductive fitness of Arabidopsis thaliana accessions under field and controlled conditions. Low temperature is one of the most important abiotic factors influencing plant fitness and geographical distribution. In addition, cold stress is known to influence crop yield and is therefore of great economic importance. Increased freezing tolerance can be acquired by the process of cold acclimation, but this may be associated with a fitness cost. To assess the influence of cold stress on the fitness of plants, long-term field trials over 5 years were performed with six natural accessions of Arabidopsis thaliana ranging from very tolerant to very sensitive to freezing. Fitness parameters, as seed yield and 1000 seed mass, were measured and correlation analyses with temperature and freezing tolerance data performed. The results were compared with fitness parameters from controlled chamber experiments over 3 years with application of cold priming and triggering conditions. Winter survival and seed yield per plant were positively correlated with temperature in field experiments. In addition, winter survival and 1000 seed mass were correlated with the cold-acclimated freezing tolerance of the selected Arabidopsis accessions. The results provide strong evidence for a trade-off between higher freezing tolerance and reproductive fitness in A. thaliana, which might have ecological impacts in the context of global warming.

Identifying loci under selection via explicit demographic models

Adaptive genetic variation is a function of both selective and neutral forces. To accurately identify adaptive loci, it is thus critical to account for demographic history. Theory suggests that signatures of selection can be inferred using the coalescent, following the premise that genealogies of selected loci deviate from neutral expectations. Here, we build on this theory to develop an analytical framework to identify loci under selection via explicit demographic models (LSD). Under this framework, signatures of selection are inferred through deviations in demographic parameters, rather than through summary statistics directly, and demographic history is accounted for explicitly. Leveraging the property of demographic models to incorporate directionality, we show that LSD can provide information on the environment in which selection acts on a population. This can prove useful in elucidating the selective processes underlying local adaptation, by characterizing genetic trade-offs and extending the concepts of antagonistic pleiotropy and conditional neutrality from ecological theory to practical application in genomic data. We implement LSD via approximate Bayesian computation and demonstrate, via simulations, that LSD (a) has high power to identify selected loci across a large range of demographic-selection regimes, (b) outperforms commonly applied genome-scan methods under complex demographies and (c) accurately infers the directionality of selection for identified candidates. Using the same simulations, we further characterize the behaviour of isolation-with-migration models conducive to the study of local adaptation under regimes of selection. Finally, we demonstrate an application of LSD by detecting loci and characterizing genetic trade-offs underlying flower colour in Antirrhinum majus.

Mapping freezing tolerance QTL in alfalfa: based on indoor phenotyping

BACKGROUND

Winter freezing temperature impacts alfalfa (Medicago sativa L.) persistence and seasonal yield and can lead to the death of the plant. Understanding the genetic mechanisms of alfalfa freezing tolerance (FT) using high-throughput phenotyping and genotyping is crucial to select suitable germplasm and develop winter-hardy cultivars. Several clones of an alfalfa F1 mapping population (3010 x CW 1010) were tested for FT using a cold chamber. The population was genotyped with SNP markers identified using genotyping-by-sequencing (GBS) and the quantitative trait loci (QTL) associated with FT were mapped on the parent-specific linkage maps. The ultimate goal is to develop non-dormant and winter-hardy alfalfa cultivars that can produce extended growth in the areas where winters are often mild.

RESULTS

Alfalfa FT screening method optimized in this experiment comprises three major steps: clone preparation, acclimation, and freezing test. Twenty clones of each genotype were tested, where 10 samples were treated with freezing temperature, and 10 were used as controls. A moderate positive correlation (r ~ 0.36, P < 0.01) was observed between indoor FT and field-based winter hardiness (WH), suggesting that the indoor FT test is a useful indirect selection method for winter hardiness of alfalfa germplasm. We detected a total of 20 QTL associated with four traits; nine for visual rating-based FT, five for percentage survival (PS), four for treated to control regrowth ratio (RR), and two for treated to control biomass ratio (BR). Some QTL positions overlapped with WH QTL reported previously, suggesting a genetic relationship between FT and WH. Some favorable QTL from the winter-hardy parent (3010) were from the potential genic region for a cold tolerance gene CBF. The BLAST alignment of a CBF sequence of M. truncatula, a close relative of alfalfa, against the alfalfa reference showed that the gene's ortholog resides around 75 Mb on chromosome 6.

CONCLUSIONS

The indoor freezing tolerance selection method reported is useful for alfalfa breeders to accelerate breeding cycles through indirect selection. The QTL and associated markers add to the genomic resources for the research community and can be used in marker-assisted selection (MAS) for alfalfa cold tolerance improvement.

Ecological factors influence balancing selection on leaf chemical profiles of a wildflower

Balancing selection is frequently invoked as a mechanism that maintains variation within and across populations. However, there are few examples of balancing selection operating on loci underpinning complex traits, which frequently display high levels of variation. We investigated mechanisms that may maintain variation in a focal polymorphism-leaf chemical profiles of a perennial wildflower (Boechera stricta, Brassicaceae)-explicitly interrogating multiple ecological and genetic processes including spatial variation in selection, antagonistic pleiotropy and frequency-dependent selection. A suite of common garden and greenhouse experiments showed that the alleles underlying variation in chemical profile have contrasting fitness effects across environments, implicating two ecological drivers of selection on chemical profile: herbivory and drought. Phenotype-environment associations and molecular genetic analyses revealed additional evidence of past selection by these drivers. Together, these data are consistent with balancing selection on chemical profile, probably caused by pleiotropic effects of secondary chemical biosynthesis genes on herbivore defence and drought response.

Life-history trade-offs and the genetic basis of fitness in Arabidopsis thaliana

Resources allocated to survival cannot be used to increase fecundity, but the extent to which this trade-off constrains adaptation depends on overall resource status. Adaptation to local environmental conditions may therefore entail the evolution of traits that increase the amount of resources available to individuals (their resource status or 'condition'). We examined the relative contribution of trade-offs and increased condition to adaptive evolution in a recombinant inbred line population of Arabidopsis thaliana planted at the native sites of the parental ecotypes in Italy and Sweden in 2 years. We estimated genetic correlations among fitness components based on genotypic means and explored their causes with QTL mapping. The local ecotype produced more seeds per fruit than did the non-local ecotype, reflected in stronger adaptive differentiation than was previously shown based on survival and fruit number only. Genetic correlations between survival and overall fecundity, and between number of fruits and number of seeds per fruit, were positive, and there was little evidence of a trade-off between seed size and number. Quantitative trait loci for these traits tended to map to the same regions of the genome and showed positive pleiotropic effects. The results indicate that adaptive differentiation between the two focal populations largely reflects the evolution of increased ability to acquire resources in the local environment, rather than shifts in the relative allocation to different life-history traits. Differentiation both in phenology and in tolerance to cold is likely to contribute to the advantage of the local genotype at the two sites.

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

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

Genome-Wide Association Studies and Transcriptome Changes during Acclimation and Deacclimation in Divergent Brassica napus Varieties

Information concerning genes and signals regulating cold acclimation processes in plants is abundant; however, less is known about genes and signals regulating the deacclimation process. A population of primarily winter B. napus varieties was used to conduct a genome-wide association study and to compare the transcriptomes from two winter B. napus varieties showing time-dependent differences in response to cold acclimation and deacclimation treatments. These studies helped to identify loci, candidate genes, and signaling processes impacting deacclimation in B. napus. GWAS identified polymorphisms at five different loci associated with freezing tolerance following deacclimation. Local linkage decay rates near these polymorphisms identified 38 possible candidate genes. Several of these genes have been reported as differentially regulated by cold stress in arabidopsis (Arabidopsis thaliana), including a calcium-binding EF-hand family protein (encoded by BnaCnng10250D) that was also differentially expressed during deacclimation in this study. Thousands of other genes differentially expressed during the acclimation and deacclimation treatments implicated processes involving oxidative stress, photosynthesis, light-regulated diurnal responses, and growth regulation. Generally, responses observed during acclimation were reversed within one week of deacclimation. The primary differences between the two winter B. napus varieties with differential deacclimation responses involved protection from oxidative stress and the ability to maintain photosynthesis.

Strong stabilizing selection on timing of germination in a Mediterranean population of Arabidopsis thaliana

PREMISE

Timing of germination can strongly influence plant fitness by affecting seedling survival and by having cascading effects on later life-history traits. In seasonal environments, the period favorable for seedling establishment and growth is limited, and timing of germination is likely to be under stabilizing selection because of conflicting selection through survival and fecundity. Moreover, optimal germination time may vary among genotypes because of inherent differences in later life-history traits.

METHODS

To examine how germination time affects survival, fecundity, and the relative fitness of two genotypes differing in time to first flower, we conducted a field experiment in an Italian population of the winter annual Arabidopsis thaliana, in which seedling establishment occurs mainly in November. We transplanted seedlings of the local genotype and of a Swedish genotype monthly from August to December and monitored survival and fecundity.

RESULTS

Only seedlings transplanted in November and December survived until reproduction, and fitness of the November cohort was 35 times higher than that of the December cohort, indicating strong stabilizing selection on timing of germination. There was no evidence of conflicting selection: seedling survival, adult survival, and fecundity were all highest in the November cohort. Moreover, the relative fitness of the two genotypes did not differ significantly between cohorts.

CONCLUSIONS

The very narrow window of opportunity for seedling establishment was related to rapid seasonal changes in soil moisture and temperature, suggesting that rate of seasonal change is an important aspect to consider for understanding spatiotemporal variation in selection on phenological traits.

Genetic dynamics of earlier maturity group emergence in south-to-north extension of Northeast China soybeans

KEY MESSAGE

This population genetic study is characterized with direct comparisons of days to flowering QTL-allele matrices between newly evolved and originally old maturity groups of soybeans to explore its evolutionary dynamics using the RTM-GWAS procedure. The Northeast China (NEC) soybeans are the major germplasm source of modern soybean production in Americas (> 80% of the world total). NEC is a relatively new soybean area in China, expanded after its nomadic status in the seventeenth century. At nine sites of four ecoregions in NEC, 361 varieties were tested for their days to flowering (DTF), a geography-sensitive trait as an indicator for maturity groups (MGs). The DTF reduced obviously along with soybeans extended to higher latitudes, ranging in 41-83 days and MG 000-III. Using the RTM-GWAS (restricted two-stage multi-locus model genome-wide association study) procedure, 81 QTLs with 342 alleles were identified, accounting for 77.85% genetic contribution (R2 = 0.01-7.74%/locus), and other 20.75% (98.60-77.85%, h2 = 98.60%) genetic variation was due to a collective of unmapped QTLs. With soybeans northward, breeding effort made the original MG I-III evolved to MG 0-00-000. In direct comparisons of QTL-allele matrices among MGs, the genetic dynamics are identified with local exotic introduction/migration (58.48%) as the first and selection against/exclusion of positive alleles causing new recombination (40.64%) as the second, while only a few allele emergence/mutation happened (0.88%, limited in MG 0, not in MG 00-000). In new MG emergence, 24 QTLs with 19 candidate genes are the major sources. A genetic potential of further DTF shortening (13-21 days) is predicted for NEC population. The QTL detection in individual ecoregions showed various ecoregion-specific QTLs-alleles/genes after co-localization treatment (removing the random environment shifting ones).

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.

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.

Introduction: Plant Cold Acclimation and Winter Survival

This introductory chapter provides a brief overview of plant freezing tolerance, cold acclimation, including subzero acclimation, and the subsequent deacclimation when plants return to warm conditions favoring growth and development. We describe the basic concepts and approaches that are currently followed to investigate these phenomena. We highlight the multidisciplinary nature of these investigations and the necessity to use methodologies from different branches of science, such as ecology, genetics, physiology, cell biology, biochemistry, and biophysics to gain a complete understanding of the complex adaptive mechanisms ultimately underlying plant winter survival.

Less photoprotection can be good in some genetic and environmental contexts

Antioxidant systems modulate oxidant-based signaling networks and excessive removal of oxidants can prevent beneficial acclimation responses. Evidence from mutant, transgenic, and locally adapted natural plant systems is used to interpret differences in the capacity for antioxidation and formulate hypotheses for future inquiry. We focus on the first line of chloroplast antioxidant defense, pre-emptive thermal dissipation of excess absorbed light (monitored as nonphotochemical fluorescence quenching, NPQ) as well as on tocopherol-based antioxidation. Findings from NPQ-deficient and tocopherol-deficient mutants that exhibited enhanced biomass production and/or enhanced foliar water-transport capacity are reviewed and discussed in the context of the impact of lower levels of antioxidation on plant performance in hot/dry conditions, under cool temperature, and in the presence of biotic stress. The complexity of cellular redox-signaling networks is related to the complexity of environmental and endogenous inputs as well as to the need for intensified training and collaboration in the study of plant-environment interactions across biological sub-disciplines.

Freezing resistance, safety margins, and survival vary among big sagebrush populations across the western United States

PREMISE

Physiological responses to temperature extremes are considered strong drivers of species' demographic responses to climate variability. Plants are typically classified as either avoiders or tolerators in their freezing-resistance mechanism, but a gradient of physiological-threshold freezing responses may exist among individuals of a species. Moreover, adaptive significance of physiological freezing responses is poorly characterized, particularly under warming conditions that relax selection on cold hardiness.

METHODS

Freezing responses were measured in winter and again for new foliage in spring for 14 populations of Artemisia tridentata collected throughout its range and planted in a warm common garden. The relationships of the freezing responses to survival were evaluated in the warm garden and in two colder gardens.

RESULTS

Winter and spring freezing resistance were not correlated and appeared to be under differing selection regimes, as evident in correlations with different population climate of origin variables. All populations resisted considerably lower temperatures in winter than in spring, with populations from more continental climates showing narrower freezing safety margins (difference in temperatures at which ice-nucleation occurs and 50% reduction in chlorophyll fluorescence occurs) in spring. Populations with greater winter freezing resistance had lower survivorship in the warmest garden, while populations with greater spring freezing resistance had lower survivorship in a colder garden.

CONCLUSIONS

These survivorship patterns relative to physiological thresholds suggest excess freezing resistance may incur a survival cost that likely relates to a trade-off between carbon gain and freezing resistance during critical periods of moisture availability. This cost has implications for seed moved from cooler to warmer environments and for plants growing in warming environments.

Unpacking Conditional Neutrality: Genomic Signatures of Selection on Conditionally Beneficial and Conditionally Deleterious Mutations

It is common to look for signatures of local adaptation in genomes by identifying loci with extreme levels of allele frequency divergence among populations. This approach to finding genes associated with local adaptation often assumes antagonistic pleiotropy, wherein alternative alleles are strongly favored in alternative environments. Conditional neutrality has been proposed as an alternative to antagonistic pleiotropy, but conditionally neutral polymorphisms are transient, and it is unclear how much outlier signal would be maintained under different forms of conditional neutrality. Here, we use individual-based simulations and a simple analytical heuristic to show that a pattern that mimics local adaptation at the phenotypic level, where each genotype has the highest fitness in its home environment, can be produced by the accumulation of mutations that are neutral in their home environment and deleterious in nonlocal environments. Because conditionally deleterious mutations likely arise at a rate many times higher than conditionally beneficial mutations, they can have a significant cumulative effect on fitness even when individual effect sizes are small. We show that conditionally deleterious mutations driving nonlocal maladaptation may be undetectable by even the most powerful genome scans, as differences in allele frequency between populations are typically small. We also explore the evolutionary effects of conditionally beneficial mutations and find that they can maintain significant signals of local adaptation, and they would be more readily detectable than conditionally deleterious mutations using conventional genome scan approaches. We discuss implications for interpreting outcomes of transplant experiments and genome scans that are used to study the genetic basis of local adaptation.

No Cost of Complexity in Bacteriophages Adapting to a Complex Environment

A long-standing prediction in evolutionary biology is that organisms experience a so-called "cost of complexity" manifested as a decreasing rate of adaptation in populations as organisms or selective environments become increasingly complex. This theory assumes the ubiquity of antagonistic pleiotropy, or trade-offs in fitness, for mutations affecting multiple traits or phenotypes. A particular manifestation of antagonism thought to be at play in adaptive dynamics involves the relationship between viral growth rate and capsid stability, an interaction that may impede the adaptation of viral pathogens to novel hosts and environments. Here, we present a comparison of the genetics of adaptation for populations of bacteriophages undergoing complete adaptive walks under both simple and complex selective conditions, with complexity being determined by the number of traits under directional selection. We found no evidence for a long-term cost of complexity in viruses experiencing complex selection, with on average at least as great a rate of adaptation under more complex conditions, and rampant evidence for synergistic, rather than antagonistic, pleiotropy. The lack of evident trade-offs between multiple phenotypes implies that emerging pathogens may be able to improve their growth in many different hosts or environments simultaneously, and to do so at a faster rate than previously anticipated.

CBF-dependent and CBF-independent regulatory pathways contribute to the differences in freezing tolerance and cold-regulated gene expression of two Arabidopsis ecotypes locally adapted to sites in Sweden and Italy

Arabidopsis thaliana (Arabidopsis) increases in freezing tolerance in response to low nonfreezing temperatures, a phenomenon known as cold acclimation. The CBF regulatory pathway, which contributes to cold acclimation, includes three genes—CBF1, CBF2 and CBF3—encoding closely-related transcription factors that regulate the expression of more than 100 genes—the CBF regulon—that impart freezing tolerance. Here we compare the CBF pathways of two Arabidopsis ecotypes collected from sites in Sweden (SW) and Italy (IT). Previous studies showed that the SW ecotype was more freezing tolerant than the IT ecotype and that the IT ecotype had a nonfunctional CBF2 gene. Here we present results establishing that the difference in CBF2 alleles contributes to the difference in freezing tolerance between the two ecotypes. However, other differences in the CBF pathway as well as CBF-independent pathways contribute the large majority of the difference in freezing tolerance between the two ecotypes. The results also provided evidence that most cold-induced CBF regulon genes in both the SW and IT ecotypes are coregulated by CBF-independent pathways. Additional analysis comparing our results with those published by others examining the Col-0 accession resulted in the identification of 44 CBF regulon genes that were conserved among the three accessions suggesting that they likely have important functions in life at low temperature. The comparison further supported the conclusion that the CBF pathway can account for a large portion of the increase in freezing tolerance that occurs with cold acclimation in a given accession, but that CBF-independent pathways can also make a major contribution.

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.

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

Local adaptation can occur due to individual genetic variants that increase the fitness of individuals in their home environments but decrease fitness in other environments [genetic trade-offs (GTs)] or genetic variants that increase fitness in one environment but have no effect in other environments [conditional neutrality (CN)]. Here, we show that GT quantitative trait loci (QTLs) for fitness between Italian and Swedish Arabidopsis thaliana exhibit strong population genomic signatures of local adaptation, including elevated levels of allele frequency differentiation, correlations to climatic variables, and recent sweeps. Highly divergent genes between Italy and Sweden populations show evidence of more recent selection in Sweden than Italy, and the biological annotations of these genes suggest interesting mechanisms underlying local adaptation.

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 F_ST SNPs colocalize with GT QTL peaks. Many of these high F_ST 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 F_ST 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.

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.

Genetic basis of photosynthetic responses to cold in two locally adapted populations of Arabidopsis thaliana

Local adaptation is common, but the traits and genes involved are often unknown. Physiological responses to cold probably contribute to local adaptation in wide-ranging species, but the genetic basis underlying natural variation in these traits has rarely been studied. Using a recombinant inbred (495 lines) mapping population from locally adapted populations of Arabidopsis thaliana from Sweden and Italy, we grew plants at low temperature and mapped quantitative trait loci (QTLs) for traits related to photosynthesis: maximal quantum efficiency (Fv/Fm), rapidly reversible photoprotection (NPQfast), and photoinhibition of PSII (NPQslow) using high-throughput, whole-plant measures of chlorophyll fluorescence. In response to cold, the Swedish line had greater values for all traits, and for every trait, large effect QTLs contributed to parental differences. We found one major QTL affecting all traits, as well as unique major QTLs for each trait. Six trait QTLs overlapped with previously published locally adaptive QTLs based on fitness measured in the native environments over 3 years. Our results demonstrate that photosynthetic responses to cold can vary dramatically within a species, and may predominantly be caused by a few QTLs of large effect. Some photosynthesis traits and QTLs probably contribute to local adaptation in this system.

Specificity in Arabidopsis thaliana recruitment of root fungal communities from soil and rhizosphere

Biotic and abiotic conditions in soil pose major constraints on growth and reproductive success of plants. Fungi are important agents in plant soil interactions but the belowground mycobiota associated with plants remains poorly understood. We grew one genotype each from Sweden and Italy of the widely-studied plant model Arabidopsis thaliana. Plants were grown under controlled conditions in organic topsoil local to the Swedish genotype, and harvested after ten weeks. Total DNA was extracted from three belowground compartments: endosphere (sonicated roots), rhizosphere and bulk soil, and fungal communities were characterized from each by amplification and sequencing of the fungal barcode region ITS2. Fungal species diversity was found to decrease from bulk soil to rhizosphere to endosphere. A significant effect of plant genotype on fungal community composition was detected only in the endosphere compartment. Despite A. thaliana being a non-mycorrhizal plant, it hosts a number of known mycorrhiza fungi in its endosphere compartment, which is also colonized by endophytic, pathogenic and saprotrophic fungi. Species in the Archaeorhizomycetes were most abundant in rhizosphere samples suggesting an adaptation to environments with high nutrient turnover for some of these species. We conclude that A. thaliana endosphere fungal communities represent a selected subset of fungi recruited from soil and that plant genotype has small but significant quantitative and qualitative effects on these communities.

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.

The Rapid Evolution of an Ohnolog Contributes to the Ecological Specialization of Incipient Yeast Species

Identifying the molecular changes that lead to ecological specialization during speciation is one of the major goals of molecular evolution. One question that remains to be thoroughly investigated is whether ecological specialization derives strictly from adaptive changes and their associated trade-offs, or from conditionally neutral mutations that accumulate under relaxed selection. We used whole-genome sequencing, genome annotation and computational analyses to identify genes that have rapidly diverged between two incipient species of Saccharomyces paradoxus that occupy different climatic regions along a south-west to north-east gradient. As candidate loci for ecological specialization, we identified genes that show signatures of adaptation and accelerated rates of amino acid substitutions, causing asymmetric evolution between lineages. This set of genes includes a glycyl-tRNA-synthetase, GRS2, which is known to be transcriptionally induced under heat stress in the model and sister species S. cerevisiae. Molecular modelling, expression analysis and fitness assays suggest that the accelerated evolution of this gene in the Northern lineage may be caused by relaxed selection. GRS2 arose during the whole-genome duplication (WGD) that occurred 100 million years ago in the yeast lineage. While its ohnolog GRS1 has been preserved in all post-WGD species, GRS2 has frequently been lost and is evolving rapidly, suggesting that the fate of this ohnolog is still to be resolved. Our results suggest that the asymmetric evolution of GRS2 between the two incipient S. paradoxus species contributes to their restricted climatic distributions and thus that ecological specialization derives at least partly from relaxed selection rather than a molecular trade-off resulting from adaptive evolution.

Acclimation of Swedish and Italian ecotypes of Arabidopsis thaliana to light intensity

This study addressed whether ecotypes of Arabidopsis thaliana from Sweden and Italy exhibited differences in foliar acclimation to high versus low growth light intensity, and compared CO₂ uptake under growth conditions with light- and CO₂-saturated intrinsic photosynthetic capacity and leaf morphological and vascular features. Differential responses between ecotypes occurred mainly at the scale of leaf architecture, with thicker leaves with higher intrinsic photosynthetic capacities and chlorophyll contents per leaf area, but no difference in photosynthetic capacity on a chlorophyll basis, in high light-grown leaves of the Swedish versus the Italian ecotype. Greater intrinsic photosynthetic capacity per leaf area in the Swedish ecotype was accompanied by a greater capacity of vascular infrastructure for sugar and water transport, but this was not associated with greater CO₂ uptake rates under growth conditions. The Swedish ecotype with its thick leaves is thus constructed for high intrinsic photosynthetic and vascular flux capacity even under growth chamber conditions that may not permit full utilization of this potential. Conversely, the Swedish ecotype was less tolerant of low growth light intensity than the Italian ecotype, with smaller rosette areas and lesser aboveground biomass accumulation in low light-grown plants. Foliar vein density and stomatal density were both enhanced by high growth light intensity with no significant difference between ecotypes, and the ratio of water to sugar conduits was also similar between the two ecotypes during light acclimation. These findings add to the understanding of the foliar vasculature's role in plant photosynthetic acclimation and adaptation.