Natural genetic variation in Arabidopsis thaliana defense metabolism genes modulates field fitness

Divergent Flowering Time Responses to Increasing Temperatures Are Associated With Transcriptome Plasticity and Epigenetic Modification Differences at FLC Promoter Region of Arabidopsis thaliana

Understanding the genetic, and transcriptomic changes that drive the phenotypic plasticity of fitness traits is a central question in evolutionary biology. In this study, we utilised 152 natural Swedish Arabidopsis thaliana accessions with re-sequenced genomes, transcriptomes and methylomes and measured flowering times (FTs) under two temperature conditions (10°C and 16°C) to address this question. We revealed that the northern accessions exhibited advanced flowering in response to decreased temperature, whereas the southern accessions delayed their flowering, indicating a divergent flowering response. This contrast in flowering responses was associated with the isothermality of their native ranges, which potentially enables the northern accessions to complete their life cycle more rapidly in years with shorter growth seasons. At the transcriptome level, we observed extensive rewiring of gene co-expression networks, with the expression of 25 core genes being associated with the mean FT and its plastic variation. Notably, variations in FLC expression sensitivity between northern and southern accessions were found to be associated with the divergence FT response. Further analysis suggests that FLC expression sensitivity is associated with differences in CG, CHG and CHH methylation at the promoter region. Overall, our study revealed the association between transcriptome plasticity and flowering time plasticity among different accessions, providing evidence for its relevance in ecological adaptation. These findings offer deeper insights into the genetics of rapid responses to environmental changes and ecological adaptation.

6Pgdh polymorphism in wild bulb mite populations: prevalence, environmental correlates and life history trade-offs

Genetic polymorphism in key metabolic genes plays a pivotal role in shaping phenotypes and adapting to varying environments. Polymorphism in the metabolic gene 6-phosphogluconate dehydrogenase (6Pgdh) in bulb mites, Rhizoglyphus robini is characterized by two alleles, S and F, that differ by a single amino acid substitution and correlate with male reproductive fitness. The S-bearing males demonstrate a reproductive advantage. Although the S allele rapidly fixes in laboratory settings, the persistence of polymorphic populations in the wild is noteworthy. This study examines the prevalence and stability of 6Pgdh polymorphism in natural populations across Poland, investigating potential environmental influences and seasonal variations. We found widespread 6Pgdh polymorphism in natural populations, with allele frequencies varying across locations and sampling dates but without clear geographical or seasonal clines. This widespread polymorphism and spatio-temporal variability may be attributed to population demography and gene flow between local populations. We found some correlation between soil properties, particularly cation content (Na, K, Ca, and Mg) and 6Pgdh allele frequencies, showcasing the connection between mite physiology and soil characteristics and highlighting the presence of environment-dependent balancing selection. We conducted experimental fitness assays to determine whether the allele providing the advantage in male-male competition has antagonistic effects on life-history traits and if these effects are temperature-dependent. We found that temperature does not differentially influence development time or juvenile survival in different 6Pgdh genotypes. This study reveals the relationship between genetic variation, environmental factors, and reproductive fitness in natural bulb mite populations, shedding light on the dynamic mechanisms governing 6Pgdh polymorphism.

Genome-wide association identifies a BAHD acyltransferase activity that assembles an ester of glucuronosylglycerol and phenylacetic acid

Genome-wide association studies (GWAS) are an effective approach to identify new specialized metabolites and the genes involved in their biosynthesis and regulation. In this study, GWAS of Arabidopsis thaliana soluble leaf and stem metabolites identified alleles of an uncharacterized BAHD-family acyltransferase (AT5G57840) associated with natural variation in three structurally related metabolites. These metabolites were esters of glucuronosylglycerol, with one metabolite containing phenylacetic acid as the acyl component of the ester. Knockout and overexpression of AT5G57840 in Arabidopsis and heterologous overexpression in Nicotiana benthamiana and Escherichia coli demonstrated that it is capable of utilizing phenylacetyl-CoA as an acyl donor and glucuronosylglycerol as an acyl acceptor. We, thus, named the protein Glucuronosylglycerol Ester Synthase (GGES). Additionally, phenylacetyl glucuronosylglycerol increased in Arabidopsis CYP79A2 mutants that overproduce phenylacetic acid and was lost in knockout mutants of UDP-sulfoquinovosyl: diacylglycerol sulfoquinovosyl transferase, an enzyme required for glucuronosylglycerol biosynthesis and associated with glycerolipid metabolism under phosphate-starvation stress. GGES is a member of a well-supported clade of BAHD family acyltransferases that arose by duplication and neofunctionalized during the evolution of the Brassicales within a larger clade that includes HCT as well as enzymes that synthesize other plant-specialized metabolites. Together, this work extends our understanding of the catalytic diversity of BAHD acyltransferases and uncovers a pathway that involves contributions from both phenylalanine and lipid metabolism.

Drought induces moderate, diverse changes in the odour of grassland species

Plants react to drought stress with numerous changes including altered emissions of volatile organic compounds (VOC) from leaves, which provide protection against oxidative tissue damage and mediate numerous biotic interactions. Despite the share of grasslands in the terrestrial biosphere, their importance as carbon sinks and their contribution to global biodiversity, little is known about the influence of drought on VOC profiles of grassland species. Using coupled gas chromatography-mass spectrometry, we analysed the odorants emitted by 22 European grassland species exposed to an eight-week-lasting drought treatment (DT; 30% water holding capacity, WHC). We focused on the odorants emitted during the light phase from whole plant shoots in their vegetative stage. Emission rates were standardised to the dry weight of each shoot. Well-watered (WW) plants (70% WHC) served as control. Drought-induced significant changes included an increase in total emission rates of plant VOC in six and a decrease in three species. Diverging effects on the number of emitted VOC (chemical richness) or on the Shannon diversity of the VOC profiles were detected in 13 species. Biosynthetic pathways-targeted analyses revealed 13 species showing drought-induced higher emission rates of VOC from one, two, three, or four major biosynthetic pathways (lipoxygenase, shikimate, mevalonate and methylerythritol phosphate pathway), while six species exhibited reduced emission rates from one or two of these pathways. Similarity trees of odorant profiles and their drought-induced changes based on a biosynthetically informed distance metric did not match species phylogeny. However, a phylogenetic signal was detected for the amount of terpenoids released by the studied species under WW and DT conditions. A comparative analysis of emission rates of single compounds released by WW and DT plants revealed significant VOC profile dissimilarities in four species only. The moderate drought-induced changes in the odorant emissions of grassland species are discussed with respect to their impact on trophic interactions across the food web. (294 words).

Complementing model species with model clades

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant tree of life continues to improve. The intersection of these 2 research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade." These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

Tomato root specialized metabolites evolved through gene duplication and regulatory divergence within a biosynthetic gene cluster

Tremendous plant metabolic diversity arises from phylogenetically restricted specialized metabolic pathways. Specialized metabolites are synthesized in dedicated cells or tissues, with pathway genes sometimes colocalizing in biosynthetic gene clusters (BGCs). However, the mechanisms by which spatial expression patterns arise and the role of BGCs in pathway evolution remain underappreciated. In this study, we investigated the mechanisms driving acylsugar evolution in the Solanaceae. Previously thought to be restricted to glandular trichomes, acylsugars were recently found in cultivated tomato roots. We demonstrated that acylsugars in cultivated tomato roots and trichomes have different sugar cores, identified root-enriched paralogs of trichome acylsugar pathway genes, and characterized a key paralog required for root acylsugar biosynthesis, SlASAT1-LIKE (SlASAT1-L), which is nested within a previously reported trichome acylsugar BGC. Last, we provided evidence that ASAT1-L arose through duplication of its paralog, ASAT1, and was trichome-expressed before acquiring root-specific expression in the Solanum genus. Our results illuminate the genomic context and molecular mechanisms underpinning metabolic diversity in plants.

Genome-wide association studies identify loci controlling specialized seed metabolites in Arabidopsis

Plants synthesize specialized metabolites to facilitate environmental and ecological interactions. During evolution, plants diversified in their potential to synthesize these metabolites. Quantitative differences in metabolite levels of natural Arabidopsis (Arabidopsis thaliana) accessions can be employed to unravel the genetic basis for metabolic traits using genome-wide association studies (GWAS). Here, we performed metabolic GWAS on seeds of a panel of 315 A. thaliana natural accessions, including the reference genotypes C24 and Col-0, for polar and semi-polar seed metabolites using untargeted ultra-performance liquid chromatography-mass spectrometry. As a complementary approach, we performed quantitative trait locus (QTL) mapping of near-isogenic introgression lines between C24 and Col-0 for specific seed specialized metabolites. Besides common QTL between seeds and leaves, GWAS revealed seed-specific QTL for specialized metabolites, indicating differences in the genetic architecture of seeds and leaves. In seeds, aliphatic methylsulfinylalkyl and methylthioalkyl glucosinolates associated with the ALKENYL HYDROXYALKYL PRODUCING loci (GS-ALK and GS-OHP) on chromosome 4 containing alkenyl hydroxyalkyl producing 2 (AOP2) and 3 (AOP3) or with the GS-ELONG locus on chromosome 5 containing methylthioalkyl malate synthase (MAM1) and MAM3. We detected two unknown sulfur-containing compounds that were also mapped to these loci. In GWAS, some of the annotated flavonoids (kaempferol 3-O-rhamnoside-7-O-rhamnoside, quercetin 3-O-rhamnoside-7-O-rhamnoside) were mapped to transparent testa 7 (AT5G07990), encoding a cytochrome P450 75B1 monooxygenase. Three additional mass signals corresponding to quercetin-containing flavonols were mapped to UGT78D2 (AT5G17050). The association of the loci and associating metabolic features were functionally verified in knockdown mutant lines. By performing GWAS and QTL mapping, we were able to leverage variation of natural populations and parental lines to study seed specialized metabolism. The GWAS data set generated here is a high-quality resource that can be investigated in further studies.

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.

A persistent major mutation in canonical jasmonate signaling is embedded in an herbivory-elicited gene network

When insect herbivores attack plants, elicitors from oral secretions and regurgitants (OS) enter wounds during feeding, eliciting defense responses. These generally require plant jasmonate (JA) signaling, specifically, a jasmonoyl-L-isoleucine (JA-Ile) burst, for their activation and are well studied in the native tobacco Nicotiana attenuata. We used intraspecific diversity captured in a 26-parent MAGIC population planted in nature and an updated genome assembly to impute natural variation in the OS-elicited JA-Ile burst linked to a mutation in the JA-Ile biosynthetic gene NaJAR4. Experiments revealed that NaJAR4 variants were associated with higher fitness in the absence of herbivores but compromised foliar defenses, with two NaJAR homologues (4 and 6) complementing each other spatially and temporally. From decade-long seed collections of natural populations, we uncovered enzymatically inactive variants occurring at variable frequencies, consistent with a balancing selection regime maintaining variants. Integrative analyses of OS-induced transcriptomes and metabolomes of natural accessions revealed that NaJAR4 is embedded in a nonlinear complex gene coexpression network orchestrating responses to OS, which we tested by silencing four hub genes in two connected coexpressed networks and examining their OS-elicited metabolic responses. Lines silenced in two hub genes (NaGLR and NaFB67) co-occurring in the NaJAR4/6 module showed responses proportional to JA-Ile accumulations; two from an adjacent module (NaERF and NaFB61) had constitutively expressed defenses with high resistance. We infer that mutations with large fitness consequences can persist in natural populations due to compensatory responses from gene networks, which allow for diversification in conserved signaling pathways and are generally consistent with predictions of an omnigene model.

Phenotypic Variation in Phytochemical Defense of Trembling Aspen in Western North America: Genetics, Development, and Geography

Trembling aspen (Populus tremuloides) is arguably the most important deciduous tree species in the Intermountain West of North America. There, as elsewhere in its range, aspen exhibits remarkable genetic variation in observable traits such as morphology and phenology. In contrast to Great Lakes populations, however, relatively little is known about phytochemical variation in western aspen. This survey of phytochemistry in western aspen was undertaken to assess how chemical expression varies among genotypes, cytotypes (diploid vs. triploid), and populations, and in response to development and mammalian browsing. We measured levels of foliar nitrogen, salicinoid phenolic glycosides (SPGs) and condensed tannins (CTs), as those constituents influence organismal interactions and ecosystem processes. Results revealed striking genotypic variation and considerable population variation, but minimal cytotype variation, in phytochemistry of western aspen. Levels of SPGs and nitrogen declined, whereas levels of CTs increased, with tree age. Browsed ramets had much higher levels of SPGs, and lower levels of CTs, than unbrowsed ramets of the same genotype. We then evaluated how composite chemical profiles of western aspen differ from those of Great Lakes aspen (assessed in earlier research). Interestingly, mature western aspen trees maintain much higher levels of SPGs, and lower levels of CTs, than Great Lakes aspen. Phenotypic variation in chemical composition of aspen - a foundation species - in the Intermountain West likely has important consequences for organismal interactions and forest ecosystem dynamics. Moreover, those consequences likely play out over spatial and temporal scales somewhat differently than have been documented for Great Lakes aspen.

Allyl isothiocyanate and 6-(methylsulfinyl) hexyl isothiocyanate contents vary among wild and cultivated wasabi (Eutrema japonium)

Wasabi (Japanese horseradish, Eutrema japonicum) is the only cultivated species in the genus Eutrema with functional components that provide a strong pungent flavor. To evaluate genetic resources for wasabi breeding, we surveyed variations in the two most abundant isothiocyanate (ITC) components in wasabi, allyl isothiocyanate (AITC) and 6-methylsulfinyl (hexyl) isothiocyanate (6-MSITC, hexaraphane). We also examined the phylogenetic relationships among 36 accessions of wild and cultivated wasabi in Japan using chloroplast DNA analysis. Our results showed that (i) the 6-MSITC content in currently cultivated wasabi accessions was significantly higher than in escaped cultivars, whereas the AITC content was not significantly different. (ii) Additionally, the 6-MSITC content in cultivated wasabi was significantly lower in the spring than during other seasons. This result suggested that the 6-MSITC content responds to environmental conditions. (iii) The phylogenetic position and the 6-MSITC content of accessions from Rebun, Hokkaido Prefecture had different profiles compared with those from southern Honshu, Japan, indicating heterogeneity of the Rebun populations from other Japanese wasabi accessions. (iv) The total content of AITC and 6-MSITC in cultivated wasabi was significantly higher than that of wild wasabi. In conclusion, old cultivars or landraces of wasabi, "zairai", are the most suitable candidates for immediate use as genetic resources.

Plant-herbivore interactions: Experimental demonstration of genetic variability in plant-plant signalling

Plant-herbivore interactions mediated by plant-plant signalling have been documented in different species but its within-species variability has hardly been quantified. Here, we tested if herbivore foraging activity on plants was influenced by a prior contact with a damaged plant and if the effect of such plant-plant signalling was variable across 113 natural genotypes of Arabidopsis thaliana. We filmed the activity of the generalist herbivore Cornu aspersum during 1 h on two plants differing only in a prior contact with a damaged plant or not. We recorded each snails' first choice, and measured its first duration on a plant, the proportion of time spent on both plants and leaf consumption. Overall, plant-plant signalling modified the foraging activity of herbivores in A. thaliana. On average, snails spent more time and consumed more of plants that experienced a prior contact with a damaged plant. However, the effects of plant-plant signalling on snail behaviour was variable: depending on genotype identity, plant-plant signalling made undamaged plants more repellant or attractive to snails. Genome-wide associations revealed that genes related to stress coping ability and jasmonate pathway were associated to this variation. Together, our findings highlight the adaptive significance of plant-plant signalling for plant-herbivore interactions.

Genetic variation underlying differential ammonium and nitrate responses in Arabidopsis thaliana

Nitrogen is an essential element required for plant growth and productivity. Understanding the mechanisms and natural genetic variation underlying nitrogen use in plants will facilitate the engineering of plant nitrogen use to maximize crop productivity while minimizing environmental costs. To understand the scope of natural variation that may influence nitrogen use, we grew 1,135 Arabidopsis thaliana natural genotypes on two nitrogen sources, nitrate and ammonium, and measured both developmental and defense metabolite traits. By using different environments and focusing on multiple traits, we identified a wide array of different nitrogen responses. These responses are associated with numerous genes, most of which were not previously associated with nitrogen responses. Only a small portion of these genes appear to be shared between environments or traits, while most are predominantly specific to a developmental or defense trait under a specific nitrogen source. Finally, by using a large population, we were able to identify unique nitrogen responses, such as preferring ammonium or nitrate, which appear to be generated by combinations of loci rather than a few large-effect loci. This suggests that it may be possible to obtain novel phenotypes in complex nitrogen responses by manipulating sets of genes with small effects rather than solely focusing on large-effect single gene manipulations.

High-energy-level metabolism and transport occur at the transition from closed to open flowers

During the maturation phase of flower development, the onset of anthesis visibly marks the transition from buds to open flowers, during which petals stretch out, nectar secretion commences, and pollination occurs. Analysis of the metabolic changes occurring during this developmental transition has primarily focused on specific classes of metabolites, such as pigments and scent emission, and far less on the whole network of primary and secondary metabolites. To investigate the metabolic changes occurring at anthesis, we performed multi-platform metabolomics alongside RNA sequencing in individual florets harvested from the main inflorescence of Arabidopsis (Arabidopsis thaliana) ecotype Col-0. To trace metabolic fluxes at the level of the whole inflorescence and individual florets, we further integrated these studies with radiolabeled experiments. These extensive analyses revealed high-energy-level metabolism and transport of carbohydrates and amino acids, supporting intense metabolic rearrangements occurring at the time of this floral transition. These comprehensive data are discussed in the context of our current understanding of the metabolic shifts underlying flower opening. We envision that this analysis will facilitate the introgression of floral metabolic traits promoting pollination in crop species for which a comprehensive knowledge of flower metabolism is still limited.

Genome-wide association mapping within a local Arabidopsis thaliana population more fully reveals the genetic architecture for defensive metabolite diversity

A paradoxical finding from genome-wide association studies (GWAS) in plants is that variation in metabolite profiles typically maps to a small number of loci, despite the complexity of underlying biosynthetic pathways. This discrepancy may partially arise from limitations presented by geographically diverse mapping panels. Properties of metabolic pathways that impede GWAS by diluting the additive effect of a causal variant, such as allelic and genetic heterogeneity and epistasis, would be expected to increase in severity with the geographical range of the mapping panel. We hypothesized that a population from a single locality would reveal an expanded set of associated loci. We tested this in a French Arabidopsis thaliana population (less than 1 km transect) by profiling and conducting GWAS for glucosinolates, a suite of defensive metabolites that have been studied in depth through functional and genetic mapping approaches. For two distinct classes of glucosinolates, we discovered more associations at biosynthetic loci than the previous GWAS with continental-scale mapping panels. Candidate genes underlying novel associations were supported by concordance between their observed effects in the TOU-A population and previous functional genetic and biochemical characterization. Local populations complement geographically diverse mapping panels to reveal a more complete genetic architecture for metabolic traits. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

High-throughput measurement of plant fitness traits with an object detection method using Faster R-CNN

Revealing the contributions of genes to plant phenotype is frequently challenging because loss-of-function effects may be subtle or masked by varying degrees of genetic redundancy. Such effects can potentially be detected by measuring plant fitness, which reflects the cumulative effects of genetic changes over the lifetime of a plant. However, fitness is challenging to measure accurately, particularly in species with high fecundity and relatively small propagule sizes such as Arabidopsis thaliana. An image segmentation-based method using the software ImageJ and an object detection-based method using the Faster Region-based Convolutional Neural Network (R-CNN) algorithm were used for measuring two Arabidopsis fitness traits: seed and fruit counts. The segmentation-based method was error-prone (correlation between true and predicted seed counts, r²  = 0.849) because seeds touching each other were undercounted. By contrast, the object detection-based algorithm yielded near perfect seed counts (r²  = 0.9996) and highly accurate fruit counts (r²  = 0.980). Comparing seed counts for wild-type and 12 mutant lines revealed fitness effects for three genes; fruit counts revealed the same effects for two genes. Our study provides analysis pipelines and models to facilitate the investigation of Arabidopsis fitness traits and demonstrates the importance of examining fitness traits when studying gene functions.

A keystone gene underlies the persistence of an experimental food web

Genes encode information that determines an organism's fitness. Yet we know little about whether genes of one species influence the persistence of interacting species in an ecological community. Here, we experimentally tested the effect of three plant defense genes on the persistence of an insect food web and found that a single allele at a single gene promoted coexistence by increasing plant growth rate, which in turn increased the intrinsic growth rates of species across multiple trophic levels. Our discovery of a "keystone gene" illustrates the need to bridge between biological scales, from genes to ecosystems, to understand community persistence.

Whole-Genome Duplication and Host Genotype Affect Rhizosphere Microbial Communities

The composition of microbial communities found in association with plants is influenced by host phenotype and genotype. However, the ways in which specific genetic architectures of host plants shape microbiomes are unknown. Genome duplication events are common in the evolutionary history of plants and influence many important plant traits, and thus, they may affect associated microbial communities. Using experimentally induced whole-genome duplication (WGD), we tested the effect of WGD on rhizosphere bacterial communities in Arabidopsis thaliana. We performed 16S rRNA amplicon sequencing to characterize differences between microbiomes associated with specific host genetic backgrounds (Columbia versus Landsberg) and ploidy levels (diploid versus tetraploid). We modeled relative abundances of bacterial taxa using a hierarchical Bayesian approach. We found that host genetic background and ploidy level affected rhizosphere community composition. We then tested to what extent microbiomes derived from a specific genetic background or ploidy level affected plant performance by inoculating sterile seedlings with microbial communities harvested from a prior generation. We found a negative effect of the tetraploid Columbia microbiome on growth of all four plant genetic backgrounds. These findings suggest an interplay between host genetic background and ploidy level and bacterial community assembly with potential ramifications for host fitness. Given the prevalence of ploidy-level variation in both wild and managed plant populations, the effects on microbiomes of this aspect of host genetic architecture could be a widespread driver of differences in plant microbiomes. IMPORTANCE Plants influence the composition of their associated microbial communities, yet the underlying host-associated genetic determinants are typically unknown. Genome duplication events are common in the evolutionary history of plants and affect many plant traits. Using Arabidopsis thaliana, we characterized how whole-genome duplication affected the composition of rhizosphere bacterial communities and how bacterial communities associated with two host plant genetic backgrounds and ploidy levels affected subsequent plant growth. We observed an interaction between ploidy level and genetic background that affected both bacterial community composition and function. This research reveals how genome duplication, a widespread genetic feature of both wild and crop plant species, influences bacterial assemblages and affects plant growth.

Contrasting gene-level signatures of selection with reproductive fitness

Selection leaves signatures in the DNA sequence of genes, with many test statistics devised to detect its action. While these statistics are frequently used to support hypotheses about the adaptive significance of particular genes, the effect these genes have on reproductive fitness is rarely quantified experimentally. Consequently, it is unclear how gene-level signatures of selection are associated with empirical estimates of gene effect on fitness. Eukaryotic datasets that permit this comparison are very limited. Using the model plant Arabidopsis thaliana, for which these resources are available, we calculated seven gene-level substitution and polymorphism-based statistics commonly used to infer selection (dN/dS, NI, DOS, Tajima's D, Fu and Li's D*, Fay and Wu's H, and Zeng's E) and, using knockout lines, compared these to gene-level estimates of effect on fitness. We found that consistent with expectations, essential genes were more likely to be classified as negatively selected. By contrast, using 379 Arabidopsis genes for which data was available, we found no evidence that genes predicted to be positively selected had a significantly different effect on fitness than genes evolving more neutrally. We discuss these results in the context of the analytic challenges posed by Arabidopsis, one of the only systems in which this study could be conducted, and advocate for examination in additional systems. These results are relevant to the evaluation of genome-wide studies across species where experimental fitness data is unavailable, as well as highlighting an increasing need for the latter.

Genetic trade-offs underlie divergent life history strategies for local adaptation in white clover

Local adaptation is common in plants, yet characterization of its underlying genetic basis is rare in herbaceous perennials. Moreover, while many plant species exhibit intraspecific chemical defence polymorphisms, their importance for local adaptation remains poorly understood. We examined the genetic architecture of local adaptation in a perennial, obligately-outcrossing herbaceous legume, white clover (Trifolium repens). This widespread species displays a well-studied chemical defence polymorphism for cyanogenesis (HCN release following tissue damage) and has evolved climate-associated cyanogenesis clines throughout its range. Two biparental F₂  mapping populations, derived from three parents collected in environments spanning the U.S. latitudinal species range (Duluth, MN, St. Louis, MO and Gainesville, FL), were grown in triplicate for two years in reciprocal common garden experiments in the parental environments (6,012 total plants). Vegetative growth and reproductive fitness traits displayed trade-offs across reciprocal environments, indicating local adaptation. Genetic mapping of fitness traits revealed a genetic architecture characterized by allelic trade-offs between environments, with 100% and 80% of fitness QTL in the two mapping populations showing significant QTL×E interactions, consistent with antagonistic pleiotropy. Across the genome there were three hotspots of QTL colocalization. Unexpectedly, we found little evidence that the cyanogenesis polymorphism contributes to local adaptation. Instead, divergent life history strategies in reciprocal environments were major fitness determinants: selection favoured early investment in flowering at the cost of multiyear survival in the southernmost site versus delayed flowering and multiyear persistence in the northern environments. Our findings demonstrate that multilocus genetic trade-offs contribute to contrasting life history characteristics that allow for local adaptation in this outcrossing herbaceous perennial.

Plant Transcription Factors Involved in Drought and Associated Stresses

Transcription factors (TFs) play a significant role in signal transduction networks spanning the perception of a stress signal and the expression of corresponding stress-responsive genes. TFs are multi-functional proteins that may simultaneously control numerous pathways during stresses in plants-this makes them powerful tools for the manipulation of regulatory and stress-responsive pathways. In recent years, the structure-function relationships of numerous plant TFs involved in drought and associated stresses have been defined, which prompted devising practical strategies for engineering plants with enhanced stress tolerance. Vast data have emerged on purposely basic leucine zipper (bZIP), WRKY, homeodomain-leucine zipper (HD-Zip), myeloblastoma (MYB), drought-response elements binding proteins/C-repeat binding factor (DREB/CBF), shine (SHN), and wax production-like (WXPL) TFs that reflect the understanding of their 3D structure and how the structure relates to function. Consequently, this information is useful in the tailored design of variant TFs that enhances our understanding of their functional states, such as oligomerization, post-translational modification patterns, protein-protein interactions, and their abilities to recognize downstream target DNA sequences. Here, we report on the progress of TFs based on their interaction pathway participation in stress-responsive networks, and pinpoint strategies and applications for crops and the impact of these strategies for improving plant stress tolerance.

Positive selective sweeps of epigenetic mutations regulating specialized metabolites in plants

DNA methylation is an important factor regulating gene expression in organisms. However, whether DNA methylation plays a key role in adaptive evolution is unknown. Here, we show evidence of naturally selected DNA methylation in Arabidopsis thaliana. In comparison with single nucleotide polymorphisms, three types of methylation—methylated CGs (mCGs), mCHGs, and mCHHs—contributed highly to variable gene expression levels among an A. thaliana population. Such variably expressed genes largely affect a large variation of specialized metabolic quantities. Among the three types of methylations, only mCGs located in promoter regions of genes associated with specialized metabolites show a selective sweep signature in the A. thaliana population. Thus, naturally selected mCGs appear to be key mutations that cause the expressional diversity associated with specialized metabolites during plant evolution.

Genetic variation, environment and demography intersect to shape Arabidopsis defense metabolite variation across Europe

Plants produce diverse metabolites to cope with the challenges presented by complex and ever-changing environments. These challenges drive the diversification of specialized metabolites within and between plant species. However, we are just beginning to understand how frequently new alleles arise controlling specialized metabolite diversity and how the geographic distribution of these alleles may be structured by ecological and demographic pressures. Here, we measure the variation in specialized metabolites across a population of 797 natural Arabidopsis thaliana accessions. We show that a combination of geography, environmental parameters, demography and different genetic processes all combine to influence the specific chemotypes and their distribution. This showed that causal loci in specialized metabolism contain frequent independently generated alleles with patterns suggesting potential within-species convergence. This provides a new perspective about the complexity of the selective forces and mechanisms that shape the generation and distribution of allelic variation that may influence local adaptation.

Since plants cannot move, they have evolved chemical defenses to help them respond to changes in their surroundings. For example, where animals run from predators, plants may produce toxins to put predators off. This approach is why plants are such a rich source of drugs, poisons, dyes and other useful substances. The chemicals plants produce are known as specialized metabolites, and they can change a lot between, and even within, plant species. The variety of specialized metabolites is a result of genetic changes and evolution over millions of years.

Evolution is a slow process, yet plants are able to rapidly develop new specialized metabolites to protect them from new threats. Even different populations of the same species produce many distinct metabolites that help them survive in their surroundings. However, the factors that lead plants to produce new metabolites are not well understood, and it is not known how this affects genetic variation.

To gain a better understanding of this process, Katz et al. studied 797 European variants of a common weed species called Arabidopsis thaliana, which is widely studied. The investigation found that many factors affect the range of specialized metabolites in each variant. These included local geography and environment, as well as genetics and population history (demography). Katz et al. revealed a pattern of relationships between the variants that could mirror their evolutionary history as the species spread and adapted to new locations.

These results highlight the complex network of factors that affect plant evolution. Rapid diversification is key to plant survival in new and changing environments and has resulted in a wide range of specialized metabolites. As such they are of interest both for studying plant evolution and for understanding their ecology. Expanding similar work to more populations and other species will broaden the scope of our ability to understand how plants adapt to their surroundings.

Arsenic Stress-Related F-Box (ASRF) gene regulates arsenic stress tolerance in Arabidopsis thaliana

Arsenic (As), a non-biodegradable contaminant, is extremely toxic to plants and animals in its inorganic form. As negatively affects plant growth and development, primarily by inducing oxidative stress through redox imbalance. Here we characterized the Arabidopsis F-box protein gene AT2G16220 (Arsenic Stress-Related F-box (ASRF)) that we identified in the genome-wide association study. The asrf mutant seedlings showed high sensitivity to arsenate (As^V) stress. As^V significantly affected asrf seedling growth when germinated on or exposed to As^V-supplemented growth regimes. As^V stress significantly induced production of reactive oxygen species and proline accumulation in asrf, so the asrf maintained high proline content, possibly for cellular protection and redox homeostasis. Heterozygous seedlings (Col-0 x asrf, F1 progeny) were relatively less affected by As^V stress than asrf mutant but showed slightly reduced growth compared with the Col-0 wild type, which suggests that the homozygous ASRF locus is important for As^V stress resistance. Transcriptome analysis involving the mutant and wild type revealed altered phosphate homeostasis in asrf seedlings, which implies that ASRF is required for maintaining phosphate and cellular- homeostasis under excess As^V. Our findings confirm the roles of ASRF in As stress tolerance in plants, for a novel way to mitigate arsenic stress.

Metabolomic Variation Aligns with Two Geographically Distinct Subpopulations of Brachypodium Distachyon before and after Drought Stress

Brachypodium distachyon (Brachypodium) is a non-domesticated model grass that has been used to assess population level genomic variation. We have previously established a collection of 55 Brachypodium accessions that were sampled to reflect five different climatic regions of Turkey; designated 1a, 1c, 2, 3 and 4. Genomic and methylomic variation differentiated the collection into two subpopulations designated as coastal and central (respectively from regions 1a, 1c and the other from 2, 3 and 4) which were linked to environmental variables such as relative precipitation. Here, we assessed how far genomic variation would be reflected in the metabolomes and if this could be linked to an adaptive trait. Metabolites were extracted from eight-week-old seedlings from each accession and assessed using flow infusion high-resolution mass spectrometry (FIE-HRMS). Principal Component Analysis (PCA) of the derived metabolomes differentiated between samples from coastal and central subpopulations. The major sources of variation between seedling from the coastal and central subpopulations were identified. The central subpopulation was typified by significant increases in alanine, aspartate and glutamate metabolism and the tricarboxylic acid (TCA) cycle. Coastal subpopulation exhibited elevated levels of the auxin, indolacetic acid and rhamnose. The metabolomes of the seedling were also determined following the imposition of drought stress for seven days. The central subpopulation exhibited a metabolomic shift in response to drought, but no significant changes were seen in the coastal one. The drought responses in the central subpopulation were typified by changes in amino acids, increasing the glutamine that could be functioning as a stress signal. There were also changes in sugars that were likely to be an osmotic counter to drought, and changes in bioenergetic metabolism. These data indicate that genomic variation in our Turkish Brachypodium collection is largely reflected as distinctive metabolomes (“metabolotypes”) through which drought tolerance might be mediated.

Regulation of glucosinolate biosynthesis

Glucosinolates are secondary defense metabolites produced by plants of the order Brassicales, which includes the model species Arabidopsis and many crop species. In the past 13 years, the regulation of glucosinolate synthesis in plants has been intensively studied, with recent research revealing complex molecular mechanisms that connect glucosinolate production with responses to other central pathways. In this review, we discuss how the regulation of glucosinolate biosynthesis is ecologically relevant for plants, how it is controlled by transcription factors, and how this transcriptional machinery interacts with hormonal, environmental, and epigenetic mechanisms. We present the central players in glucosinolate regulation, MYB and basic helix-loop-helix transcription factors, as well as the plant hormone jasmonate, which together with other hormones and environmental signals allow the coordinated and rapid regulation of glucosinolate genes. Furthermore, we highlight the regulatory connections between glucosinolates, auxin, and sulfur metabolism and discuss emerging insights and open questions on the regulation of glucosinolate biosynthesis.

Making sense of Integrated Pest Management (IPM) in the light of evolution

Integrated Pest Management (IPM) is a holistic approach to combat pests (including herbivores, pathogens, and weeds) using a combination of preventive and curative actions, and only applying synthetic pesticides when there is an urgent need. Just as the recent recognition that an evolutionary perspective is useful in medicine to understand and predict interactions between hosts, diseases, and medical treatments, we argue that it is crucial to integrate an evolutionary framework in IPM to develop efficient and reliable crop protection strategies that do not lead to resistance development in herbivores, pathogens, and weeds. Such a framework would not only delay resistance evolution in pests, but also optimize each element of the management and increase the synergies between them. Here, we outline key areas within IPM that would especially benefit from a thorough evolutionary understanding. In addition, we discuss the difficulties and advantages of enhancing communication among research communities rooted in different biological disciplines and between researchers and society. Furthermore, we present suggestions that could advance implementation of evolutionary principles in IPM and thus contribute to the development of sustainable agriculture that is resilient to current and emerging pests.

Independent evolution of ancestral and novel defenses in a genus of toxic plants (Erysimum, Brassicaceae)

Phytochemical diversity is thought to result from coevolutionary cycles as specialization in herbivores imposes diversifying selection on plant chemical defenses. Plants in the speciose genus Erysimum (Brassicaceae) produce both ancestral glucosinolates and evolutionarily novel cardenolides as defenses. Here we test macroevolutionary hypotheses on co-expression, co-regulation, and diversification of these potentially redundant defenses across this genus. We sequenced and assembled the genome of E. cheiranthoides and foliar transcriptomes of 47 additional Erysimum species to construct a phylogeny from 9868 orthologous genes, revealing several geographic clades but also high levels of gene discordance. Concentrations, inducibility, and diversity of the two defenses varied independently among species, with no evidence for trade-offs. Closely related, geographically co-occurring species shared similar cardenolide traits, but not glucosinolate traits, likely as a result of specific selective pressures acting on each defense. Ancestral and novel chemical defenses in Erysimum thus appear to provide complementary rather than redundant functions.

Epistatic Transcription Factor Networks Differentially Modulate Arabidopsis Growth and Defense

Plants integrate internal and external signals to finely coordinate growth and defense for maximal fitness within a complex environment. A common model suggests that growth and defense show a trade-offs relationship driven by energy costs. However, recent studies suggest that the coordination of growth and defense likely involves more conditional and intricate connections than implied by the trade-off model. To explore how a transcription factor (TF) network may coordinate growth and defense, we used a high-throughput phenotyping approach to measure growth and flowering in a set of single and pairwise mutants previously linked to the aliphatic glucosinolate (GLS) defense pathway. Supporting a link between growth and defense, 17 of the 20 tested defense-associated TFs significantly influenced plant growth and/or flowering time. The TFs' effects were conditional upon the environment and age of the plant, and more critically varied across the growth and defense phenotypes for a given genotype. In support of the coordination model of growth and defense, the TF mutant's effects on short-chain aliphatic GLS and growth did not display a simple correlation. We propose that large TF networks integrate internal and external signals and separately modulate growth and the accumulation of the defensive aliphatic GLS.

Natural variation in Arabidopsis shoot branching plasticity in response to nitrate supply affects fitness

The capacity of organisms to tune their development in response to environmental cues is pervasive in nature. This phenotypic plasticity is particularly striking in plants, enabled by their modular and continuous development. A good example is the activation of lateral shoot branches in Arabidopsis, which develop from axillary meristems at the base of leaves. The activity and elongation of lateral shoots depends on the integration of many signals both external (e.g. light, nutrient supply) and internal (e.g. the phytohormones auxin, strigolactone and cytokinin). Here, we characterise natural variation in plasticity of shoot branching in response to nitrate supply using two diverse panels of Arabidopsis lines. We find extensive variation in nitrate sensitivity across these lines, suggesting a genetic basis for variation in branching plasticity. High plasticity is associated with extreme branching phenotypes such that lines with the most branches on high nitrate have the fewest under nitrate deficient conditions. Conversely, low plasticity is associated with a constitutively moderate level of branching. Furthermore, variation in plasticity is associated with alternative life histories with the low plasticity lines flowering significantly earlier than high plasticity lines. In Arabidopsis, branching is highly correlated with fruit yield, and thus low plasticity lines produce more fruit than high plasticity lines under nitrate deficient conditions, whereas highly plastic lines produce more fruit under high nitrate conditions. Low and high plasticity, associated with early and late flowering respectively, can therefore be interpreted alternative escape vs mitigate strategies to low N environments. The genetic architecture of these traits appears to be highly complex, with only a small proportion of the estimated genetic variance detected in association mapping.

Transcriptional Variation in Glucosinolate Biosynthetic Genes and Inducible Responses to Aphid Herbivory on Field-Grown Arabidopsis thaliana

Recently, increasing attempts have been made to understand how plant genes function in natura. In this context, transcriptional profiles represent plant physiological status in response to environmental stimuli. Herein, we combined high-throughput RNA-Seq with insect survey data on 19 accessions of Arabidopsis thaliana grown at a field site in Switzerland. We found that genes with the gene ontology (GO) annotations of "glucosinolate biosynthetic process" and "response to insects" were most significantly enriched, and the expression of these genes was highly variable among plant accessions. Nearly half of the total expression variation in the glucosinolate biosynthetic genes (AOPs, ESM1, ESP, and TGG1) was explained by among-accession variation. Of these genes, the expression level of AOP3 differed among Col-0 accession individuals depending on the abundance of the mustard aphid (Lipaphis erysimi). We also found that the expression of the major cis-jasmone activated gene CYP81D11 was positively correlated with the number of flea beetles (Phyllotreta striolata and Phyllotreta atra). Combined with the field RNA-Seq data, bioassays confirmed that AOP3 was up-regulated in response to attack by mustard aphids. The combined results from RNA-Seq and our ecological survey illustrate the feasibility of using field transcriptomics to detect an inducible defense, providing a first step towards an in natura understanding of biotic interactions involving phenotypic plasticity.

RNA-Seq analysis of soft rush (Juncus effusus): transcriptome sequencing, de novo assembly, annotation, and polymorphism identification

BACKGROUND

Juncus effusus L. (family: Juncaceae; order: Poales) is a helophytic rush growing in temperate damp or wet terrestrial habitats and is of almost cosmopolitan distribution. The species has been studied intensively with respect to its interaction with co-occurring plants as well as microbes being involved in major biogeochemical cycles. J. effusus has biotechnological value as component of Constructed Wetlands where the plant has been employed in phytoremediation of contaminated water. Its genome has not been sequenced.

RESULTS

In this study we carried out functional annotation and polymorphism analysis of de novo assembled RNA-Seq data from 18 genotypes using 249 million paired-end Illumina HiSeq reads and 2.8 million 454 Titanium reads. The assembly comprised 158,591 contigs with a mean contig length of 780 bp. The assembly was annotated using the dammit! annotation pipeline, which queries the databases OrthoDB, Pfam-A, Rfam, and runs BUSCO (Benchmarking Single-Copy Ortholog genes). In total, 111,567 contigs (70.3%) were annotated with functional descriptions, assigned gene ontology terms, and conserved protein domains, which resulted in 30,932 non-redundant gene sequences. Results of BUSCO and KEGG pathway analyses were similar for J. effusus as for the well-studied members of the Poales, Oryza sativa and Sorghum bicolor. A total of 566,433 polymorphisms were identified in transcribed regions with an average frequency of 1 polymorphism in every 171 bases.

CONCLUSIONS

The transcriptome assembly was of high quality and genome coverage was sufficient for global analyses. This annotated knowledge resource can be utilized for future gene expression analysis, genomic feature comparisons, genotyping, primer design, and functional genomics in J. effusus.

Plant trichomes and a single gene GLABRA1 contribute to insect community composition on field-grown Arabidopsis thaliana

BACKGROUND

Genetic variation in plants alters insect abundance and community structure in the field; however, little is known about the importance of a single gene among diverse plant genotypes. In this context, Arabidopsis trichomes provide an excellent system to discern the roles of natural variation and a key gene, GLABRA1, in shaping insect communities. In this study, we transplanted two independent glabrous mutants (gl1-1 and gl1-2) and 17 natural accessions of Arabidopsis thaliana to two localities in Switzerland and Japan.

RESULTS

Fifteen insect species inhabited the plant accessions, with the insect community composition significantly attributed to variations among plant accessions. The total abundance of leaf-chewing herbivores was negatively correlated with trichome density at both field sites, while glucosinolates had variable effects on leaf chewers between the sites. Interestingly, there was a parallel tendency for the abundance of leaf chewers to be higher on gl1-1 and gl1-2 than on their different parental accessions, Ler-1 and Col-0, respectively. Furthermore, the loss of function in the GLABRA1 gene significantly decreased the resistance of plants to the two predominant chewers; flea beetles and turnip sawflies.

CONCLUSIONS

Overall, our results indicate that insect community composition significantly varies among A. thaliana accessions across two distant field sites, with GLABRA1 playing a key role in altering the abundance of leaf-chewing herbivores. Given that such a trichome variation is widely observed in Brassicaceae plants, the present study exemplifies the community-wide effect of a single plant gene on crucifer-feeding insects in the field.

Nonlinear phenotypic variation uncovers the emergence of heterosis in Arabidopsis thaliana

Heterosis describes the phenotypic superiority of hybrids over their parents in traits related to agronomic performance and fitness. Understanding and predicting nonadditive inheritance such as heterosis is crucial for evolutionary biology as well as for plant and animal breeding. However, the physiological bases of heterosis remain debated. Moreover, empirical data in various species have shown that diverse genetic and molecular mechanisms are likely to explain heterosis, making it difficult to predict its emergence and amplitude from parental genotypes alone. In this study, we examined a model of physiological dominance initially proposed by Sewall Wright to explain the nonadditive inheritance of traits like metabolic fluxes at the cellular level. We evaluated Wright's model for two fitness-related traits at the whole-plant level, growth rate and fruit number, using 450 hybrids derived from crosses among natural accessions of A. thaliana. We found that allometric relationships between traits constrain phenotypic variation in a nonlinear and similar manner in hybrids and accessions. These allometric relationships behave predictably, explaining up to 75% of heterosis amplitude, while genetic distance among parents at best explains 7%. Thus, our findings are consistent with Wright's model of physiological dominance and suggest that the emergence of heterosis on plant performance is an intrinsic property of nonlinear relationships between traits. Furthermore, our study highlights the potential of a geometric approach of phenotypic relationships for predicting heterosis of major components of crop productivity and yield.

Soil microbes alter plant fitness under competition and drought

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

A high-throughput method of analyzing multiple plant defensive compounds in minimized sample mass

PREMISE OF THE STUDY

Current methods for quantifying herbivore-induced alterations in plant biochemistry are often unusable by researchers due to practical constraints. We present a cost-effective, high-throughput protocol to quantify multiple biochemical responses from small plant tissue samples using spectrophotometric techniques.

METHODS AND RESULTS

Using Solanum lycopersicum and Medicago polymorpha leaves pre- and post-herbivory, we demonstrate that our protocol quantifies common plant defense responses: peroxidase production, polyphenol oxidase production, reactive oxygen species production, total protein production, and trypsin-like protease inhibition activity.

CONCLUSIONS

Current protocols can require 500 mg of tissue, but our assays detect activity in less than 10 mg. Our protocol takes two people 6 h to run any of the assays on 300 samples in triplicate, or all of the assays on 20 samples. Our protocol enables researchers to plan complex experiments that compare local versus systemic plant responses, quantify environmental and genetic variation, and measure population-level variation.

Regulation of Root Angle and Gravitropism

Regulation of plant root angle is critical for obtaining nutrients and water and is an important trait for plant breeding. A plant’s final, long-term root angle is the net result of a complex series of decisions made by a root tip in response to changes in nutrient availability, impediments, the gravity vector and other stimuli. When a root tip is displaced from the gravity vector, the short-term process of gravitropism results in rapid reorientation of the root toward the vertical. Here, we explore both short- and long-term regulation of root growth angle, using natural variation in tomato to identify shared and separate genetic features of the two responses. Mapping of expression quantitative trait loci mapping and leveraging natural variation between and within species including Arabidopsis suggest a role for PURPLE ACID PHOSPHATASE 27 and CELL DIVISION CYCLE 73 in determining root angle.

Extreme copy number variation at a tRNA ligase gene affecting phenology and fitness in yellow monkeyflowers

Copy number variation (CNV) is a major part of the genetic diversity segregating within populations, but remains poorly understood relative to single nucleotide variation. Here, we report on a tRNA ligase gene (Migut.N02091; RLG1a) exhibiting unprecedented, and fitness-relevant, CNV within an annual population of the yellow monkeyflower Mimulus guttatus. RLG1a variation was associated with multiple traits in pooled population sequencing (PoolSeq) scans of phenotypic and phenological cohorts. Resequencing of inbred lines revealed intermediate-frequency three-copy variants of RLG1a (trip+; 5/35 = 14%), and trip+ lines exhibited elevated RLG1a expression under multiple conditions. trip+ carriers, in addition to being over-represented in late-flowering and large-flowered PoolSeq populations, flowered later under stressful conditions in a greenhouse experiment (p < 0.05). In wild population samples, we discovered an additional rare RLG1a variant (high+) that carries 250-300 copies of RLG1a totalling ~5.7 Mb (20-40% of a chromosome). In the progeny of a high+ carrier, Mendelian segregation of diagnostic alleles and qPCR-based copy counts indicate that high+ is a single tandem array unlinked to the single-copy RLG1a locus. In the wild, high+ carriers had highest fitness in two particularly dry and/or hot years (2015 and 2017; both p < 0.01), while single-copy individuals were twice as fecund as either CNV type in a lush year (2016: p < 0.005). Our results demonstrate fluctuating selection on CNVs affecting phenological traits in a wild population, suggest that plant tRNA ligases mediate stress-responsive life-history traits, and introduce a novel system for investigating the molecular mechanisms of gene amplification.

Modulation of ACD6 dependent hyperimmunity by natural alleles of an Arabidopsis thaliana NLR resistance gene

Plants defend themselves against pathogens by activating an array of immune responses. Unfortunately, immunity programs may also cause unintended collateral damage to the plant itself. The quantitative disease resistance gene ACCELERATED CELL DEATH 6 (ACD6) serves to balance growth and pathogen resistance in natural populations of Arabidopsis thaliana. An autoimmune allele, ACD6-Est, which strongly reduces growth under specific laboratory conditions, is found in over 10% of wild strains. There is, however, extensive variation in the strength of the autoimmune phenotype expressed by strains with an ACD6-Est allele, indicative of genetic modifiers. Quantitative genetic analysis suggests that ACD6 activity can be modulated in diverse ways, with different strains often carrying different large-effect modifiers. One modifier is SUPPRESSOR OF NPR1-1, CONSTITUTIVE 1 (SNC1), located in a highly polymorphic cluster of nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes, which are prototypes for qualitative disease resistance genes. Allelic variation at SNC1 correlates with ACD6-Est activity in multiple accessions, and a common structural variant affecting the NL linker sequence can explain differences in SNC1 activity. Taken together, we find that an NLR gene can mask the activity of an ACD6 autoimmune allele in natural A. thaliana populations, thereby linking different arms of the plant immune system.

Plants defend themselves against pathogens by activating immune responses. Unfortunately, these can cause unintended collateral damage to the plant itself. Nevertheless, some wild plants have genetic variants that confer a low threshold for the activation of immunity. While these enable a plant to respond particularly quickly to pathogen attack, such variants might be potentially dangerous. We are investigating one such variant of the immune gene ACCELERATED CELL DEATH 6 (ACD6) in the plant Arabidopsis thaliana. We discovered that there are variants at other genetic loci that can mask the effects of an overly active ACD6 gene. One of these genes, SUPPRESSOR OF NPR1-1, CONSTITUTIVE 1 (SNC1), codes for a known immune receptor. The SNC1 variant that attenuates ACD6 activity is rather common in A. thaliana populations, suggesting that new combinations of the hyperactive ACD6 variant and this antagonistic SNC1 variant will often arise by natural crosses. Similarly, because the two genes are unlinked, outcrossing will often lead to the hyperactive ACD6 variants being unmasked again. We propose that allelic diversity at SNC1 contributes to the maintenance of the hyperactive ACD6 variant in natural A. thaliana populations.

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

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

Localization of the glucosinolate biosynthetic enzymes reveals distinct spatial patterns for the biosynthesis of indole and aliphatic glucosinolates

Glucosinolates constitute the primary defense metabolites in Arabidopsis thaliana (Arabidopsis). Indole and aliphatic glucosinolates, biosynthesized from tryptophan and methionine, respectively, are known to serve distinct biological functions. Although all genes in the biosynthetic pathways are identified, and it is known where glucosinolates are stored, it has remained elusive where glucosinolates are produced at the cellular and tissue level. To understand how the spatial organization of the different glucosinolate biosynthetic pathways contributes to their distinct biological functions, we investigated the localization of enzymes of the pathways under constitutive conditions and, for indole glucosinolates, also under induced conditions, by analyzing the spatial distribution of several fluorophore-tagged enzymes at the whole plant and the cellular level. We show that key steps in the biosynthesis of the different types of glucosinolates are localized in distinct cells in separate as well as overlapping vascular tissues. The presence of glucosinolate biosynthetic enzymes in parenchyma cells of the vasculature may assign new defense-related functions to these cell types. The knowledge gained in this study is an important prerequisite for understanding the orchestration of chemical defenses from site of synthesis to site of storage and potential (re)mobilization upon attack.

Comparison of the Relative Potential for Epigenetic and Genetic Variation To Contribute to Trait Stability

The theoretical ability of epigenetic variation to influence the heritable variation of complex traits is gaining traction in the study of adaptation. This theory posits that epigenetic marks can control adaptive phenotypes but the relative potential of epigenetic variation in comparison to genetic variation in these traits is not presently understood. To compare the potential of epigenetic and genetic variation in adaptive traits, we analyzed the influence of DNA methylation variation on the accumulation of chemical defense compounds glucosinolates from the order Brassicales. Several decades of work on glucosinolates has generated extensive knowledge about their synthesis, regulation, genetic variation and contribution to fitness establishing this pathway as a model pathway for complex adaptive traits. Using high-throughput phenotyping with a randomized block design of ddm1 derived Arabidopsis thaliana epigenetic Recombinant Inbred Lines, we measured the correlation between DNA methylation variation and mean glucosinolate variation and within line stochastic variation. Using this information, we identified epigenetic Quantitative Trait Loci that contained specific Differentially Methylated Regions associated with glucosinolate traits. This showed that variation in DNA methylation correlates both with levels and variance of glucosinolates and flowering time with trait-specific loci. By conducting a meta-analysis comparing the results to different genetically variable populations, we conclude that the influence of DNA methylation variation on these adaptive traits is much lower than the corresponding impact of standing genetic variation. As such, selective pressure on these traits should mainly affect standing genetic variation to lead to adaptation.

Plant Defenses against Herbivory: Closing the Fitness Gap

Many morphological and chemical features of plants are classified as plant defenses against herbivores. By definition, plant defenses should increase a plant's fitness (i.e., its contribution to the gene pool of the next generation) as a function of herbivory. Over the past years, substantial progress has been made in understanding and manipulating the mechanistic basis of many putative plant defense traits. However, most plant defenses are still characterized by proximate variables such as herbivore performance or plant damage rather than actual fitness. Determining fitness benefits as a function of herbivory therefore remains a major knowledge gap that must be filled to understand the ecology and evolution of plant defenses.