Convergently evolved metabolites are new to me but not to my attackers

This article is a Commentary on Younkin et al. (2024), 242: 2719–2733.

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.

Specificity and breadth of plant specialized metabolite-microbe interactions

Plant specialized metabolites shape plant interactions with the environment including plant-microbe interactions. While we often group compounds into generic classes, it is the precise structure of a compound that creates a specific role in plant-microbe or-pathogen interactions. Critically, the structure guides definitive targets in individual interactions, yet single compounds are not limited to singular mechanistic targets allowing them to influence interactions across broad ranges of attackers, from bacteria to fungi to animals. Further, the direction of the effect can be altered by counter evolution within the interacting organism leading to single compounds being both beneficial and detrimental. Thus, the benefit of a single compound to a host needs to be assessed by measuring the net benefit across all interactions while in each specific interaction. Factoring this complexity for single compounds in plant-microbe interactions with the massive expansion in our identification of specialized metabolite pathways means that we need systematic studies to classify the full breadth of activities. Only with this full biological knowledge we can develop mechanistic, ecological, and evolutionary models to understand how plant specialized metabolites fully influence plant-microbe and plant-biotic interactions more broadly.

Genetic and molecular landscapes of the generalist phytopathogen Botrytis cinerea

Botrytis cinerea Pers. Fr. (teleomorph: Botryotinia fuckeliana) is a necrotrophic fungal pathogen that attacks a wide range of plants. This updated pathogen profile explores the extensive genetic diversity of B. cinerea, highlights the progress in genome sequencing, and provides current knowledge of genetic and molecular mechanisms employed by the fungus to attack its hosts. In addition, we also discuss recent innovative strategies to combat B. cinerea.

TAXONOMY

Kingdom: Fungi, phylum: Ascomycota, subphylum: Pezizomycotina, class: Leotiomycetes, order: Helotiales, family: Sclerotiniaceae, genus: Botrytis, species: cinerea.

HOST RANGE

B. cinerea infects almost all of the plant groups (angiosperms, gymnosperms, pteridophytes, and bryophytes). To date, 1606 plant species have been identified as hosts of B. cinerea.

GENETIC DIVERSITY

This polyphagous necrotroph has extensive genetic diversity at all population levels shaped by climate, geography, and plant host variation.

PATHOGENICITY

Genetic architecture of virulence and host specificity is polygenic using multiple weapons to target hosts, including secretory proteins, complex signal transduction pathways, metabolites, and mobile small RNA.

DISEASE CONTROL STRATEGIES

Efforts to control B. cinerea, being a high-diversity generalist pathogen, are complicated. However, integrated disease management strategies that combine cultural practices, chemical and biological controls, and the use of appropriate crop varieties will lessen yield losses. Recently, studies conducted worldwide have explored the potential of small RNA as an efficient and environmentally friendly approach for combating grey mould. However, additional research is necessary, especially on risk assessment and regulatory frameworks, to fully harness the potential of this technology.

Large-scale identification of novel transcriptional regulators of the aliphatic glucosinolate pathway in Arabidopsis

Aliphatic glucosinolates are a large group of plant secondary metabolites characteristic of Brassicaceae, including the model plant Arabidopsis. The diverse and complex degradation products of aliphatic glucosinolates contribute to plant responses to herbivory, pathogen attack, and environmental stresses. Most of the biosynthesis genes in the aliphatic glucosinolate pathway have been cloned in Arabidopsis, and the research focus has recently shifted to the regulatory mechanisms controlling aliphatic glucosinolate accumulation. Up till now, more than 40 transcriptional regulators have been identified as regulating the aliphatic glucosinolate pathway, but many more novel regulators likely remain to be discovered based on research evidence over the past decade. In the current study, we took a systemic approach to functionally test 155 candidate transcription factors in Arabidopsis identified by yeast one-hybrid assay, and successfully validated at least 30 novel regulators that could significantly influence the accumulation of aliphatic glucosinolates in our experimental set-up. We also showed that the regulators of the aliphatic glucosinolate pathway have balanced positive and negative effects, and glucosinolate metabolism and plant development can be coordinated. Our work is the largest scale effort so far to validate transcriptional regulators of a plant secondary metabolism pathway, and provides new insights into how the highly diverse plant secondary metabolism is regulated at the transcriptional level.

Making watercress (Nasturtium officinale) cropping sustainable: genomic insights into enhanced phosphorus use efficiency in an aquatic crop

Watercress (Nasturtium officinale) is a nutrient-dense salad crop with high antioxidant capacity and glucosinolate concentration and with the potential to contribute to nutrient security as a locally grown outdoor aquatic crop in northern temperate climates. However, phosphate-based fertilizers used to support plant growth contribute to the eutrophication of aquatic habitats, often pristine chalk streams, downstream of farms, increasing pressure to minimize fertilizer use and develop a more phosphorus-use efficient (PUE) crop. Here, we grew genetically distinct watercress lines selected from a bi-parental mapping population on a commercial watercress farm either without additional phosphorus (P-) or under a commercial phosphate-based fertilizer regime (P+), to decipher effects on morphology, nutritional profile, and the transcriptome. Watercress plants sustained shoot yield in P- conditions, through enhanced root biomass, but with shorter stems and smaller leaves. Glucosinolate concentration was not affected by P- conditions, but both antioxidant capacity and the concentration of sugars and starch in shoot tissue were enhanced. We identified two watercress breeding lines, with contrasting strategies for enhanced PUE: line 60, with highly plastic root systems and increased root growth in P-, and line 102, maintaining high yield irrespective of P supply, but less plastic. RNA-seq analysis revealed a suite of genes involved in cell membrane remodeling, root development, suberization, and phosphate transport as potential future breeding targets for enhanced PUE. We identified watercress gene targets for enhanced PUE for future biotechnological and breeding approaches enabling less fertilizer inputs and reduced environmental damage from watercress cultivation.

Epigenomic divergence correlates with sequence polymorphism in Arabidopsis paralogs

Processes affecting rates of sequence polymorphism are fundamental to the evolution of gene duplicates. The relationship between gene activity and sequence polymorphism can influence the likelihood that functionally redundant gene copies are co-maintained in stable evolutionary equilibria vs other outcomes such as neofunctionalization. Here, we investigate genic variation in epigenome-associated polymorphism rates in Arabidopsis thaliana and consider whether these affect the evolution of gene duplicates. We compared the frequency of sequence polymorphism and patterns of genetic differentiation between genes classified by exon methylation patterns: unmethylated (unM), gene-body methylated (gbM), and transposon-like methylated (teM) states, which reflect divergence in gene expression. We found that the frequency of polymorphism was higher in teM (transcriptionally repressed, tissue-specific) genes and lower in gbM (active, constitutively expressed) genes. Comparisons of gene duplicates were largely consistent with genome-wide patterns - gene copies that exhibit teM accumulate more variation, evolve faster, and are in chromatin states associated with reduced DNA repair. This relationship between expression, the epigenome, and polymorphism may lead to the breakdown of equilibrium states that would otherwise maintain genetic redundancies. Epigenome-mediated polymorphism rate variation may facilitate the evolution of novel gene functions in duplicate paralogs maintained over evolutionary time.

Is specialized metabolite regulation specialized?

Recent technical and theoretical advances have generated an explosion in the identification of specialized metabolite pathways. In comparison, our understanding of how these pathways are regulated is relatively lagging. This and the relatively young age of specialized metabolite pathways has partly contributed to a default and common paradigm whereby specialized metabolite regulation is theorized as relatively simple with a few key transcription factors and the compounds are non-regulatory end-products. In contrast, studies into model specialized metabolites, such as glucosinolates, are beginning to identify a new understanding whereby specialized metabolites are highly integrated into the plants' core metabolic, physiological, and developmental pathways. This model includes a greatly extended compendium of transcription factors controlling the pathway, key transcription factors that co-evolve with the pathway and simultaneously control core metabolic and developmental components, and finally the compounds themselves evolve regulatory connections to integrate into the plants signaling machinery. In this review, these concepts are illustrated using studies in the glucosinolate pathway within the Brassicales. This suggests that the broader community needs to reconsider how they do or do not integrate specialized metabolism into the regulatory network of their study species.

Polygenic pathogen networks influence transcriptional plasticity in the Arabidopsis-Botrytis pathosystem

Bidirectional flow of information shapes the outcome of the host-pathogen interactions and depends on the genetics of each organism. Recent work has begun to use co-transcriptomic studies to shed light on this bidirectional flow, but it is unclear how plastic the co-transcriptome is in response to genetic variation in both the host and pathogen. To study co-transcriptome plasticity, we conducted transcriptomics using natural genetic variation in the pathogen, Botrytis cinerea, and large-effect genetic variation abolishing defense signaling pathways within the host, Arabidopsis thaliana. We show that genetic variation in the pathogen has a greater influence on the co-transcriptome than mutations that abolish defense signaling pathways in the host. Genome-wide association mapping using the pathogens' genetic variation and both organisms' transcriptomes allowed an assessment of how the pathogen modulates plasticity in response to the host. This showed that the differences in both organism's responses were linked to trans-expression quantitative trait loci (eQTL) hotspots within the pathogen's genome. These hotspots control gene sets in either the host or pathogen and show differential allele sensitivity to the host's genetic variation rather than qualitative host specificity. Interestingly, nearly all the trans-eQTL hotspots were unique to the host or pathogen transcriptomes. In this system of differential plasticity, the pathogen mediates the shift in the co-transcriptome more than the host.

Causes of Mutation Rate Variability in Plant Genomes

Mutation is the source of all heritable diversity, the essential material of evolution and breeding. While mutation rates are often regarded as constant, variability in mutation rates has been observed at nearly every level-varying across mutation types, genome locations, gene functions, epigenomic contexts, environmental conditions, genotypes, and species. This mutation rate variation arises from differential rates of DNA damage, repair, and transposable element activation and insertion that together produce what is measured by DNA mutation rates. We review historical and recent investigations into the causes and consequences of mutation rate variability in plants by focusing on the mechanisms shaping this variation. Emerging mechanistic models point to the evolvability of mutation rate variation across genomes via mechanisms that target DNA repair, shaping the diversification of plants at phenotypic and genomic scales.

The rhizosphere microbiome and host plant glucosinolates exhibit feedback cycles in Brassica rapa

The rhizosphere microbiome influences many aspects of plant fitness, including production of secondary compounds and defence against insect herbivores. Plants also modulate the composition of the microbial community in the rhizosphere via secretion of root exudates. We tested both the effect of the rhizosphere microbiome on plant traits, and host plant effects on rhizosphere microbes using recombinant inbred lines (RILs) of Brassica rapa that differ in production of glucosinolates (GLS), secondary metabolites that contribute to defence against insect herbivores. First, we investigated the effect of genetic variation in GLS production on the composition of the rhizosphere microbiome. Using a Bayesian Dirichlet-multinomial regression model (DMBVS), we identified both negative and positive associations between bacteria from six genera and the concentration of five GLS compounds produced in plant roots. Additionally, we tested the effects of microbial inoculation (an intact vs. disrupted soil microbiome) on GLS production and insect damage in these RILs. We found a significant microbial treatment × genotype interaction, in which total GLS was higher in the intact relative to the disrupted microbiome treatment in some RILs. However, despite differences in GLS production between microbial treatments, we observed no difference in insect damage between treatments. Together, these results provide evidence for a full feedback cycle of plant-microbe interactions mediated by GLS; that is, GLS compounds produced by the host plant "feed-down" to influence rhizosphere microbial community and rhizosphere microbes "feed-up" to influence GLS production.

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.

Integrated omics reveal novel functions and underlying mechanisms of the receptor kinase FERONIA in Arabidopsis thaliana

The receptor kinase FERONIA (FER) is a versatile regulator of plant growth and development, biotic and abiotic stress responses, and reproduction. To gain new insights into the molecular interplay of these processes and to identify new FER functions, we carried out quantitative transcriptome, proteome, and phosphoproteome profiling of Arabidopsis (Arabidopsis thaliana) wild-type and fer-4 loss-of-function mutant plants. Gene ontology terms for phytohormone signaling, abiotic stress, and biotic stress were significantly enriched among differentially expressed transcripts, differentially abundant proteins, and/or misphosphorylated proteins, in agreement with the known roles for FER in these processes. Analysis of multiomics data and subsequent experimental evidence revealed previously unknown functions for FER in endoplasmic reticulum (ER) body formation and glucosinolate biosynthesis. FER functions through the transcription factor NAI1 to mediate ER body formation. FER also negatively regulates indole glucosinolate biosynthesis, partially through NAI1. Furthermore, we found that a group of abscisic acid (ABA)-induced transcription factors is hypophosphorylated in the fer-4 mutant and demonstrated that FER acts through the transcription factor ABA INSENSITIVE5 (ABI5) to negatively regulate the ABA response during cotyledon greening. Our integrated omics study, therefore, reveals novel functions for FER and provides new insights into the underlying mechanisms of FER function.

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.

Aphid Species and Feeding Location on Canola Influences the Impact of Glucosinolates on a Native Lady Beetle Predator

Aphids that attack canola (Brassica napus L.) exhibit feeding preferences for different parts of canola plants, which may be associated with brassica-specific glucosinolates. However, this idea remains untested. Furthermore, canola aphid species employ different strategies for tolerating glucosinolates. While the green peach aphid, Myzus persicae (Sulzer) (Hemiptera: Aphididae), excretes glucosinolates, the cabbage aphid Brevicoryne brassicae (L.) (Hemiptera: Aphididae) sequesters them. Given the different detoxification mechanisms, we predicted that both aphid species and aphid feeding location would affect prey suitability for larvae of the predator, Hippodamia convergens (Guérin-Méneville) (Coleoptera: Coccinellidae). We hypothesized that aphids, specifically glucosinolate-sequestering cabbage aphid, reared on reproductive structures that harbor higher glucosinolates concentrations would have greater negative effects on predators than those reared on vegetative structures which have lower levels of glucosinolates, and that the impact of aphid feeding location would vary depending on the prey detoxification mechanism. To test these predictions, we conducted experiments to compare 1) glucosinolates profiles between B. brassicae and M. persicae reared on reproductive and vegetative canola structures, 2) aphid population growth on each structure, and 3) their subsequent impact on fitness traits of H. convergens. Results indicate that the population growth of both aphids was greater on reproductive structures, with B. brassicae having the highest population growth. B. brassicae reared on reproductive structures had the highest concentrations of glucosinolates, and the greatest adverse effects on H. convergens. These findings suggest that both aphid-prey species and feeding location on canola could influence populations of this predator and, thus, its potential for biological control of canola aphids.

Plant Responses Underlying Timely Specialized Metabolites Induction of Brassica Crops

A large subset of plant stress-signaling pathways, including those related with chemical defense production, exhibit diurnal or circadian oscillations. However the extent to which diurnal or circadian time influences the stress mediated accumulation of plant specialized metabolites remains largely unknown. Because plant responses to physical stress (e.g., wounding) is considered a common component of mounting a response against a broad range of environmental stresses, including herbivory, we have utilized mechanical wounding as the stress stimulus to determine the direct contribution of time of day on the induced defenses of Brassica crops. We analyzed glucosinolates (GSLs) from leaves of broccoli (Brassica oleracea) and turnip greens (Brassica rapa) following exposure to mechanical wounding at dawn (ZT0), mid-day (ZT4), and dusk (ZT8). Several GSLs differentially accumulated and their changes depended upon the time of day at wounding was performed. This response varied considerably between species. In a parallel experiment, we investigated whether diurnal activation of Brassica phytochemicals in response to wounding might prime plants against herbivore attack. Results showed that maximal response of plant chemical defense against larvae of the generalist pest Mamestra brassicae occurred at ZT0 in broccoli and ZT8 in turnip greens. Metabolome analysis for global trends of time dependent compounds showed that sulfur-containing phytochemicals, GSL hydrolysis products, auxin-signaling components, and other metabolites activators of plant disease resistance (nicotinamide and pipecolate) had important contributions to the responses of M. brassicae feeding behavior in broccoli at morning. Overall, the findings in this study highlight a significant role for time of day in the wound stress responsive metabolome, which can in turn affect plant-herbivore interactions.

A genome-scale TF-DNA interaction network of transcriptional regulation of Arabidopsis primary and specialized metabolism

Plant metabolism is more complex relative to individual microbes. In single-celled microbes, transcriptional regulation by single transcription factors (TFs) is sufficient to shift primary metabolism. Corresponding genome-level transcriptional regulatory maps of metabolism reveal the underlying design principles responsible for these shifts as a model in which master regulators largely coordinate specific metabolic pathways. Plant primary and specialized metabolism occur within innumerable cell types, and their reactions shift depending on internal and external cues. Given the importance of plants and their metabolites in providing humanity with food, fiber, and medicine, we set out to develop a genome-scale transcriptional regulatory map of Arabidopsis metabolic genes. A comprehensive set of protein-DNA interactions between Arabidopsis thaliana TFs and gene promoters in primary and specialized metabolic pathways were mapped. To demonstrate the utility of this resource, we identified and functionally validated regulators of the tricarboxylic acid (TCA) cycle. The resulting network suggests that plant metabolic design principles are distinct from those of microbes. Instead, metabolism appears to be transcriptionally coordinated via developmental- and stress-conditional processes that can coordinate across primary and specialized metabolism. These data represent the most comprehensive resource of interactions between TFs and metabolic genes in plants.

Quantitative interactions: the disease outcome of Botrytis cinerea across the plant kingdom

Botrytis cinerea is a fungal pathogen that causes necrotic disease on more than a thousand known hosts widely spread across the plant kingdom. How B. cinerea interacts with such extensive host diversity remains largely unknown. To address this question, we generated an infectivity matrix of 98 strains of B. cinerea on 90 genotypes representing eight host plants. This experimental infectivity matrix revealed that the disease outcome is largely explained by variations in either the host resistance or pathogen virulence. However, the specific interactions between host and pathogen account for 16% of the disease outcome. Furthermore, the disease outcomes cluster among genotypes of a species but are independent of the relatedness between hosts. When analyzing the host specificity and virulence of B. cinerea, generalist strains are predominant. In this fungal necrotroph, specialization may happen by a loss in virulence on most hosts rather than an increase of virulence on a specific host. To uncover the genetic architecture of Botrytis host specificity and virulence, a genome-wide association study (GWAS) was performed and revealed up to 1492 genes of interest. The genetic architecture of these traits is widespread across the B. cinerea genome. The complexity of the disease outcome might be explained by hundreds of functionally diverse genes putatively involved in adjusting the infection to diverse hosts.

Genome size evolution is associated with climate seasonality and glucosinolates, but not life history, soil nutrients or range size, across a clade of mustards

BACKGROUND AND AIMS

We investigate patterns of evolution of genome size across a morphologically and ecologically diverse clade of Brassicaceae, in relation to ecological and life history traits. While numerous hypotheses have been put forward regarding autecological and environmental factors that could favour small vs. large genomes, a challenge in understanding genome size evolution in plants is that many hypothesized selective agents are intercorrelated.

METHODS

We contribute genome size estimates for 47 species of Streptanthus Nutt. and close relatives, and take advantage of many data collections for this group to assemble data on climate, life history, soil affinity and composition, geographic range and plant secondary chemistry to identify simultaneous correlates of variation in genome size in an evolutionary framework. We assess models of evolution across clades and use phylogenetically informed analyses as well as model selection and information criteria approaches to identify variables that can best explain genome size variation in this clade.

KEY RESULTS

We find differences in genome size and heterogeneity in its rate of evolution across subclades of Streptanthus and close relatives. We show that clade-wide genome size is positively associated with climate seasonality and glucosinolate compounds. Model selection and information criteria approaches identify a best model that includes temperature seasonality and fraction of aliphatic glucosinolates, suggesting a possible role for genome size in climatic adaptation or a role for biotic interactions in shaping the evolution of genome size. We find no evidence supporting hypotheses of life history, range size or soil nutrients as forces shaping genome size in this system.

CONCLUSIONS

Our findings suggest climate seasonality and biotic interactions as potential forces shaping the evolution of genome size and highlight the importance of evaluating multiple factors in the context of phylogeny to understand the effect of possible selective agents on genome size.

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.

Fine mapping identifies NAD-ME1 as a candidate underlying a major locus controlling temporal variation in primary and specialized metabolism in Arabidopsis

Plant metabolism is modulated by a complex interplay between internal signals and external cues. A major goal of all quantitative metabolomic studies is to clone the underlying genes to understand the mechanistic basis of this variation. Using fine-scale genetic mapping, in this work we report the identification and initial characterization of NAD-DEPENDENT MALIC ENZYME 1 (NAD-ME1) as the candidate gene underlying the pleiotropic network Met.II.15 quantitative trait locus controlling variation in plant metabolism and circadian clock outputs in the Bay × Sha Arabidopsis population. Transcript abundance and promoter analysis in NAD-ME1^Bay-0 and NAD-ME1^Sha alleles confirmed allele-specific expression that appears to be due a polymorphism disrupting a putative circadian cis-element binding site. Analysis of transfer DNA insertion lines and heterogeneous inbred families showed that transcript variation of the NAD-ME1 gene led to temporal shifts of tricarboxylic acid cycle intermediates, glucosinolate (GSL) accumulation, and altered regulation of several GSL biosynthesis pathway genes. Untargeted metabolomic analyses revealed complex regulatory networks of NAD-ME1 dependent upon the daytime. The mutant led to shifts in plant primary metabolites, cell wall components, isoprenoids, fatty acids, and plant immunity phytochemicals, among others. Our findings suggest that NAD-ME1 may act as a key gene to coordinate plant primary and secondary metabolism in a time-dependent manner.

The Arabidopsis thaliana nucleotide sugar transporter GONST2 is a functional homolog of GONST1

Glycosylinositolphosphorylceramides (GIPCs) are the predominant lipid in the outer leaflet of the plasma membrane. Characterized GIPC glycosylation mutants have severe or lethal plant phenotypes. However, the function of the glycosylation is unclear. Previously, we characterized Arabidopsis thaliana GONST1 and showed that it was a nucleotide sugar transporter which provides GDP-mannose for GIPC glycosylation. gonst1 has a severe growth phenotype, as well as a constitutive defense response. Here, we characterize a mutant in GONST1's closest homolog, GONST2. The gonst2-1 allele has a minor change to GIPC headgroup glycosylation. Like other reported GIPC glycosylation mutants, gonst1-1gonst2-1 has reduced cellulose, a cell wall polymer that is synthesized at the plasma membrane. The gonst2-1 allele has increased resistance to a biotrophic pathogen Golovinomyces orontii but not the necrotrophic pathogen Botrytis cinerea. Expression of GONST2 under the GONST1 promoter can rescue the gonst1 phenotype, indicating that GONST2 has a similar function to GONST1 in providing GDP-D-Man for GIPC mannosylation.

Red-light is an environmental effector for mutualism between begomovirus and its vector whitefly

Environments such as light condition influence the spread of infectious diseases by affecting insect vector behavior. However, whether and how light affects the host defense which further affects insect preference and performance, remains unclear, nor has been demonstrated how pathogens co-adapt light condition to facilitate vector transmission. We previously showed that begomoviral βC1 inhibits MYC2-mediated jasmonate signaling to establish plant-dependent mutualism with its insect vector. Here we show red-light as an environmental catalyzer to promote mutualism of whitefly-begomovirus by stabilizing βC1, which interacts with PHYTOCHROME-INTERACTING FACTORS (PIFs) transcription factors. PIFs positively control plant defenses against whitefly by directly binding to the promoter of terpene synthase genes and promoting their transcription. Moreover, PIFs interact with MYC2 to integrate light and jasmonate signaling and regulate the transcription of terpene synthase genes. However, begomovirus encoded βC1 inhibits PIFs' and MYC2' transcriptional activity via disturbing their dimerization, thereby impairing plant defenses against whitefly-transmitted begomoviruses. Our results thus describe how a viral pathogen hijacks host external and internal signaling to enhance the mutualistic relationship with its insect vector.

Plant Secondary Metabolites as Defenses, Regulators, and Primary Metabolites: The Blurred Functional Trichotomy

The plant kingdom produces hundreds of thousands of low molecular weight organic compounds. Based on the assumed functions of these compounds, the research community has classified them into three overarching groups: primary metabolites, which are directly required for plant growth; secondary (or specialized) metabolites, which mediate plant-environment interactions; and hormones, which regulate organismal processes and metabolism. For decades, this functional trichotomy of plant metabolism has shaped theory and experimentation in plant biology. However, exact biochemical boundaries between these different metabolite classes were never fully established. A new wave of genetic and chemical studies now further blurs these boundaries by demonstrating that secondary metabolites are multifunctional; they can function as potent regulators of plant growth and defense as well as primary metabolites sensu lato. Several adaptive scenarios may have favored this functional diversity for secondary metabolites, including signaling robustness and cost-effective storage and recycling. Secondary metabolite multifunctionality can provide new explanations for ontogenetic patterns of defense production and can refine our understanding of plant-herbivore interactions, in particular by accounting for the discovery that adapted herbivores misuse plant secondary metabolites for multiple purposes, some of which mirror their functions in plants. In conclusion, recent work unveils the limits of our current functional classification system for plant metabolites. Viewing secondary metabolites as integrated components of metabolic networks that are dynamically shaped by environmental selection pressures and transcend multiple trophic levels can improve our understanding of plant metabolism and plant-environment interactions.

flasher, a novel mutation in a glucosinolate modifying enzyme, conditions changes in plant architecture and hormone homeostasis

Meristem function is underpinned by numerous genes that affect hormone levels, ultimately controlling phyllotaxy, the transition to flowering and general growth properties. Class I KNOX genes are major contributors to this process, promoting cytokinin biosynthesis but repressing gibberellin production to condition a replication competent state. We identified a suppressor mutant of the KNOX1 mutant brevipedicellus (bp) that we termed flasher (fsh), which promotes stem and pedicel elongation, suppresses early senescence, and negatively affects reproductive development. Map-based cloning and complementation tests revealed that fsh is due to an E40K change in the flavin monooxygenase GS-OX5, a gene encoding a glucosinolate (GSL) modifying enzyme. In vitro enzymatic assays revealed that fsh poorly converts substrate to product, yet the levels of several GSLs are higher in the suppressor line, implicating FSH in feedback control of GSL flux. FSH is expressed predominantly in the vasculature in patterns that do not significantly overlap those of BP, implying a non-cell autonomous mode of meristem control via one or more GSL metabolites. Hormone analyses revealed that cytokinin levels are low in bp, but fsh restores cytokinin levels to near normal by activating cytokinin biosynthesis genes. In addition, jasmonate levels in the fsh suppressor are significantly lower than in bp, which is likely due to elevated expression of JA inactivating genes. These observations suggest the involvement of the GSL pathway in generating one or more negative effectors of growth that influence inflorescence architecture and fecundity by altering the balance of hormonal regulators.

FRS7 and FRS12 recruit NINJA to regulate expression of glucosinolate biosynthesis genes

The sessile lifestyle of plants requires accurate physiology adjustments to be able to thrive in a changing environment. Plants integrate environmental timing signals to control developmental and stress responses. Here, we identified Far1 Related Sequence (FRS) 7 and FRS12, two transcriptional repressors that accumulate in short-day conditions, as regulators of Arabidopsis glucosinolate (GSL) biosynthesis. Loss of function of FRS7 and FRS12 results in plants with increased amplitudes of diurnal expression of GSL pathway genes. Protein interaction analyses revealed that FRS7 and FRS12 recruit the NOVEL INTERACTOR OF JAZ (NINJA) to assemble a transcriptional repressor complex. Genetic and molecular evidence demonstrated that FRS7, FRS12 and NINJA jointly regulate the expression of GSL biosynthetic genes, and thus constitute a molecular mechanism that modulates specialized metabolite accumulation.

Diverse Allyl Glucosinolate Catabolites Independently Influence Root Growth and Development

Glucosinolates (GSLs) are sulfur-containing defense metabolites produced in the Brassicales, including the model plant Arabidopsis (Arabidopsis thaliana). Previous work suggests that specific GSLs may function as signals to provide direct feedback regulation within the plant to calibrate defense and growth. These GSLs include allyl-GSL, a defense metabolite that is one of the most widespread GSLs in Brassicaceae and has also been associated with growth inhibition. Here we show that at least three separate potential catabolic products of allyl-GSL or closely related compounds affect growth and development by altering different mechanisms influencing plant development. Two of the catabolites, raphanusamic acid and 3-butenoic acid, differentially affect processes downstream of the auxin signaling cascade. Another catabolite, acrylic acid, affects meristem development by influencing the progression of the cell cycle. These independent signaling events propagated by the different catabolites enable the plant to execute a specific response that is optimal to any given environment.

mGWAS Uncovers Gln-Glucosinolate Seed-Specific Interaction and its Role in Metabolic Homeostasis

Gln is a key player in plant metabolism. It is one of the major free amino acids that is transported into the developing seed and is central for nitrogen metabolism. However, Gln natural variation and its regulation and interaction with other metabolic processes in seeds remain poorly understood. To investigate the latter, we performed a metabolic genome-wide association study (mGWAS) of Gln-related traits measured from the dry seeds of the Arabidopsis (Arabidopsis thaliana) diversity panel using all potential ratios between Gln and the other members of the Glu family as traits. This semicombinatorial approach yielded multiple candidate genes that, upon further analysis, revealed an unexpected association between the aliphatic glucosinolates (GLS) and the Gln-related traits. This finding was confirmed by an independent quantitative trait loci mapping and statistical analysis of the relationships between the Gln-related traits and the presence of specific GLS in seeds. Moreover, an analysis of Arabidopsis mutants lacking GLS showed an extensive seed-specific impact on Gln levels and composition that manifested early in seed development. The elimination of GLS in seeds was associated with a large effect on seed nitrogen and sulfur homeostasis, which conceivably led to the Gln response. This finding indicates that both Gln and GLS play key roles in shaping the seed metabolic homeostasis. It also implies that select secondary metabolites might have key functions in primary seed metabolism. Finally, our study shows that an mGWAS performed on dry seeds can uncover key metabolic interactions that occur early in seed development.

Pathogen Genetic Control of Transcriptome Variation in the Arabidopsis thaliana - Botrytis cinerea Pathosystem

In plant-pathogen relations, disease symptoms arise from the interaction of the host and pathogen genomes. Host-pathogen functional gene interactions are well described, whereas little is known about how the pathogen genetic variation modulates both organisms' transcriptomes. To model and generate hypotheses on a generalist pathogen control of gene expression regulation, we used the Arabidopsis thaliana-Botrytis cinerea pathosystem and the genetic diversity of a collection of 96 B. cinerea isolates. We performed expression-based genome-wide association (eGWA) for each of 23,947 measurable transcripts in Arabidopsis (host), and 9267 measurable transcripts in B. cinerea (pathogen). Unlike other eGWA studies, we detected a relative absence of locally acting expression quantitative trait loci (cis-eQTL), partly caused by structural variants and allelic heterogeneity hindering their identification. This study identified several distantly acting trans-eQTL linked to eQTL hotspots dispersed across Botrytis genome that altered only Botrytis transcripts, only Arabidopsis transcripts, or transcripts from both species. Gene membership in the trans-eQTL hotspots suggests links between gene expression regulation and both known and novel virulence mechanisms in this pathosystem. Genes annotated to these hotspots provide potential targets for blocking manipulation of the host response by this ubiquitous generalist necrotrophic pathogen.

Using networks to identify and interpret natural variation

Studies on natural variation and network biology inherently work to summarize vast amounts of information and data. The combination of these two areas of study while creating datasets of immense complexity is critical to their mutual progress. Networks are necessary as a way to work to reduce the dimensionality inherent in natural variation with 100 s to 1000 s of genotypes. Correspondingly natural variation is essential for testing how networks may or may not be shared across individuals or species. Advances in this area of cross-fertilization including using networks directly as phenotypes and the use of networks to help in prioritizing candidate gene validation efforts. Interesting new observations on frequent presence-absence variation in gene content and adaptation is beginning to highlight the potential for natural variation in network presence-absence. This review attempts to delve into these new insights.

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.

Auxin-sensitive Aux/IAA proteins mediate drought tolerance in Arabidopsis by regulating glucosinolate levels

A detailed understanding of abiotic stress tolerance in plants is essential to provide food security in the face of increasingly harsh climatic conditions. Glucosinolates (GLSs) are secondary metabolites found in the Brassicaceae that protect plants from herbivory and pathogen attack. Here we report that in Arabidopsis, aliphatic GLS levels are regulated by the auxin-sensitive Aux/IAA repressors IAA5, IAA6, and IAA19. These proteins act in a transcriptional cascade that maintains expression of GLS levels when plants are exposed to drought conditions. Loss of IAA5/6/19 results in reduced GLS levels and decreased drought tolerance. Further, we show that this phenotype is associated with a defect in stomatal regulation. Application of GLS to the iaa5,6,19 mutants restores stomatal regulation and normal drought tolerance. GLS action is dependent on the receptor kinase GHR1, suggesting that GLS may signal via reactive oxygen species. These results provide a novel connection between auxin signaling, GLS levels and drought response.

Auxin-sensitive Aux/IAA proteins mediate drought tolerance in Arabidopsis by regulating glucosinolate levels

A detailed understanding of abiotic stress tolerance in plants is essential to provide food security in the face of increasingly harsh climatic conditions. Glucosinolates (GLSs) are secondary metabolites found in the Brassicaceae that protect plants from herbivory and pathogen attack. Here we report that in Arabidopsis, aliphatic GLS levels are regulated by the auxin-sensitive Aux/IAA repressors IAA5, IAA6, and IAA19. These proteins act in a transcriptional cascade that maintains expression of GLS levels when plants are exposed to drought conditions. Loss of IAA5/6/19 results in reduced GLS levels and decreased drought tolerance. Further, we show that this phenotype is associated with a defect in stomatal regulation. Application of GLS to the iaa5,6,19 mutants restores stomatal regulation and normal drought tolerance. GLS action is dependent on the receptor kinase GHR1, suggesting that GLS may signal via reactive oxygen species. These results provide a novel connection between auxin signaling, GLS levels and drought response.

Brassicaceae produce glucosinolates to protect against herbivory and pathogens. Here the authors show that auxin-sensitive Aux/IAA repressor proteins regulate aliphatic glucosinolate levels in Arabidopsis and this promotes stomatal closure via reactive oxygen species during drought stress.

Plant Networks as Traits and Hypotheses: Moving Beyond Description

Biology relies on the central thesis that the genes in an organism encode molecular mechanisms that combine with stimuli and raw materials from the environment to create a final phenotypic expression representative of the genomic programming. While conceptually simple, the genotype-to-phenotype linkage in a eukaryotic organism relies on the interactions of thousands of genes and an environment with a potentially unknowable level of complexity. Modern biology has moved to the use of networks in systems biology to try to simplify this complexity to decode how an organism's genome works. Previously, biological networks were basic ways to organize, simplify, and analyze data. However, recent advances are allowing networks to move beyond description and become phenotypes or hypotheses in their own right. This review discusses these efforts, like mapping responses across biological scales, including relationships among cellular entities, and the direct use of networks as traits or hypotheses.

Viruses mobilize plant immunity to deter nonvector insect herbivores

A parasite-infected host may promote performance of associated insect vectors; but possible parasite effects on nonvector insects have been largely unexplored. Here, we show that Begomovirus, the largest genus of plant viruses and transmitted exclusively by whitefly, reprogram plant immunity to promote the fitness of the vector and suppress performance of nonvector insects (i.e., cotton bollworm and aphid). Infected plants accumulated begomoviral βC1 proteins in the phloem where they were bound to the plant transcription factor WRKY20. This viral hijacking of WRKY20 spatiotemporally redeployed plant chemical immunity within the leaf and had the asymmetrical benefiting effects on the begomoviruses and its whitefly vectors while negatively affecting two nonvector competitors. This type of interaction between a parasite and two types of herbivores, i.e., vectors and nonvectors, occurs widely in various natural and agricultural ecosystems; thus, our results have broad implications for the ecological significance of parasite-vector-host tripartite interactions.

Plant-necrotroph co-transcriptome networks illuminate a metabolic battlefield

A central goal of studying host-pathogen interaction is to understand how host and pathogen manipulate each other to promote their own fitness in a pathosystem. Co-transcriptomic approaches can simultaneously analyze dual transcriptomes during infection and provide a systematic map of the cross-kingdom communication between two species. Here we used the Arabidopsis-B. cinerea pathosystem to test how plant host and fungal pathogen interact at the transcriptomic level. We assessed the impact of genetic diversity in pathogen and host by utilization of a collection of 96 isolates infection on Arabidopsis wild-type and two mutants with jasmonate or salicylic acid compromised immunities. We identified ten B. cinereagene co-expression networks (GCNs) that encode known or novel virulence mechanisms. Construction of a dual interaction network by combining four host- and ten pathogen-GCNs revealed potential connections between the fungal and plant GCNs. These co-transcriptome data shed lights on the potential mechanisms underlying host-pathogen interaction.

Interactions of Tomato and Botrytis cinerea Genetic Diversity: Parsing the Contributions of Host Differentiation, Domestication, and Pathogen Variation

Although the impacts of crop domestication on specialist pathogens are well known, less is known about the interaction of crop variation and generalist pathogens. To study how genetic variation within a crop affects plant resistance to generalist pathogens, we infected a collection of wild and domesticated tomato accessions with a genetically diverse population of the generalist pathogen Botrytis cinerea We quantified variation in lesion size of 97 B. cinerea genotypes (isolates) on six domesticated tomato genotypes (Solanum lycopersicum) and six wild tomato genotypes (Solanum pimpinellifolium). Lesion size was significantly affected by large effects of the host and pathogen's genotype, with a much smaller contribution of domestication. This pathogen collection also enables genome-wide association mapping of B. cinerea Genome-wide association mapping of the pathogen showed that virulence is highly polygenic and involves a diversity of mechanisms. Breeding against this pathogen would likely require the use of diverse isolates to capture all possible mechanisms. Critically, we identified a subset of B. cinerea genes where allelic variation was linked to altered virulence against wild versus domesticated tomato, as well as loci that could handle both groups. This generalist pathogen already has a large collection of allelic variation that must be considered when designing a breeding program.

Digital Imaging Combined with Genome-Wide Association Mapping Links Loci to Plant-Pathogen Interaction Traits

Plant resistance to generalist pathogens with broad host ranges, such as Botrytis cinerea (Botrytis), is typically quantitative and highly polygenic. Recent studies have begun to elucidate the molecular genetic basis of plant-pathogen interactions using commonly measured traits, including lesion size and/or pathogen biomass. However, with the advent of digital imaging and high-throughput phenomics, there are a large number of additional traits available to study quantitative resistance. In this study, we used high-throughput digital imaging analysis to investigate previously poorly characterized visual traits of plant-pathogen interactions related to disease resistance using the Arabidopsis (Arabidopsis thaliana)/Botrytis pathosystem. From a large collection of visual lesion trait measurements, we focused on color, shape, and size to test how these aspects of the Arabidopsis/Botrytis interaction are genetically related. Through genome-wide association mapping in Arabidopsis, we show that lesion color and shape are genetically separable traits associated with plant disease resistance. Moreover, by employing defined mutants in 23 candidate genes identified from the genome-wide association mapping, we demonstrate links between loci and each of the different plant-pathogen interaction traits. These results expand our understanding of the functional mechanisms driving plant disease resistance.

Transcriptional regulation of nitrogen-associated metabolism and growth

Nitrogen is an essential macronutrient for plant growth and basic metabolic processes. The application of nitrogen-containing fertilizer increases yield, which has been a substantial factor in the green revolution¹. Ecologically, however, excessive application of fertilizer has disastrous effects such as eutrophication². A better understanding of how plants regulate nitrogen metabolism is critical to increase plant yield and reduce fertilizer overuse. Here we present a transcriptional regulatory network and twenty-one transcription factors that regulate the architecture of root and shoot systems in response to changes in nitrogen availability. Genetic perturbation of a subset of these transcription factors revealed coordinate transcriptional regulation of enzymes involved in nitrogen metabolism. Transcriptional regulators in the network are transcriptionally modified by feedback via genetic perturbation of nitrogen metabolism. The network, genes and gene-regulatory modules identified here will prove critical to increasing agricultural productivity.

Network-Guided Discovery of Extensive Epistasis between Transcription Factors Involved in Aliphatic Glucosinolate Biosynthesis

Plants use diverse mechanisms influenced by vast regulatory networks of indefinite scale to adapt to their environment. These regulatory networks have an unknown potential for epistasis between genes within and across networks. To test for epistasis within an adaptive trait genetic network, we generated and tested 47 Arabidopsis thaliana double mutant combinations for 20 transcription factors, which all influence the accumulation of aliphatic glucosinolates, the defense metabolites that control fitness. The epistatic combinations were used to test if there is more or less epistasis depending on gene membership within the same or different phenotypic subnetworks. Extensive epistasis was observed between the transcription factors, regardless of subnetwork membership. Metabolite accumulation displayed antagonistic epistasis, suggesting the presence of a buffering mechanism. Epistasis affecting enzymatic estimated activity was highly conditional on the tissue and environment and shifted between both antagonistic and synergistic forms. Transcriptional analysis showed that epistasis shifts depend on how the trait is measured. Because the 47 combinations described here represent a small sampling of the potential epistatic combinations in this genetic network, there is potential for significantly more epistasis. Additionally, the main effect of the individual gene was not predictive of the epistatic effects, suggesting that there is a need for further studies.

An evolutionarily young defense metabolite influences the root growth of plants via the ancient TOR signaling pathway

To optimize fitness a plant should monitor its metabolism to appropriately control growth and defense. Primary metabolism can be measured by the universally conserved TOR (Target of Rapamycin) pathway to balance growth and development with the available energy and nutrients. Recent work suggests that plants may measure defense metabolites to potentially provide a strategy ensuring fast reallocation of resources to coordinate plant growth and defense. There is little understanding of mechanisms enabling defense metabolite signaling. To identify mechanisms of defense metabolite signaling, we used glucosinolates, an important class of plant defense metabolites. We report novel signaling properties specific to one distinct glucosinolate, 3-hydroxypropylglucosinolate across plants and fungi. This defense metabolite, or derived compounds, reversibly inhibits root growth and development. 3-hydroxypropylglucosinolate signaling functions via genes in the ancient TOR pathway. If this event is not unique, this raises the possibility that other evolutionarily new plant metabolites may link to ancient signaling pathways.

Plants, like all organisms, must invest their resources carefully. Growing new roots or shoots may allow a plant to better exploit its environment. But a plant should never leave itself vulnerable to disease. As such, there must be a balance between allocating resources to growth or to defense.

Brassicas like cabbage, Brussels sprouts and wasabi use unique compounds called glucosinolates to protect themselves against pests and disease-causing microbes. These same compounds give these vegetables their distinctive flavors, and they are the source of many of the health benefits linked to eating these vegetables. Yet it was not known if glucosinolates could also affect a plant’s growth and development.

Malinovsky et al. tested a number of purified glucosinolates with the model plant Arabidopsis thaliana, and found that one (called 3-hydroxypropylglucosinolate) caused the plants to grow with stunted roots. When 10 other species of plant were grown with this glucosinolate, almost all had shorter-than-normal roots. The effect was not limited to plants; baker’s yeast also grew less when its liquid media contained the plant-derived compound.

The reason glucosinolates can protect plants against insect pests, provide us with health benefits, and widely inhibit growth is most likely because they have evolved to interact with proteins that are found in many different organisms.Indeed, through experiments with mutant Arabidopsis plants, Malinovsky et al. revealed that their glucosinolate influences the TOR complex. This complex of proteins works in an ancient and widespread signaling pathway that balances growth and development with the available energy and nutrients in organisms ranging from humans to yeast to plants.

The TOR complex plays such a vital role in living cells that problems with this complex have been linked to diseases such as cancer and heart disease. Importantly, the chemical structure of this glucosinolate is unlike other compounds that have already been tested against the TOR complex. As such, it is possible that this glucosinolate might lead to new drugs for a range of human diseases. Further, as this compound affects plant growth, it could also act as a starting point for new herbicides.

Together these findings show how studying molecules made in model organisms and understanding how they function can lead to the identification of new compounds and targets with an unexpectedly wide range of potential uses.

Initiation of ER Body Formation and Indole Glucosinolate Metabolism by the Plastidial Retrograde Signaling Metabolite, MEcPP

Plants have evolved tightly regulated signaling networks to respond and adapt to environmental perturbations, but the nature of the signaling hub(s) involved have remained an enigma. We have previously established that methylerythritol cyclodiphosphate (MEcPP), a precursor of plastidial isoprenoids and a stress-specific retrograde signaling metabolite, enables cellular readjustments for high-order adaptive functions. Here, we specifically show that MEcPP promotes two Brassicaceae-specific traits, namely endoplasmic reticulum (ER) body formation and induction of indole glucosinolate (IGs) metabolism selectively, via transcriptional regulation of key regulators NAI1 for ER body formation and MYB51/122 for IGs biosynthesis). The specificity of MEcPP is further confirmed by the lack of induction of wound-inducible ER body genes as well as IGs by other altered methylerythritol phosphate pathway enzymes. Genetic analyses revealed MEcPP-mediated COI1-dependent induction of these traits. Moreover, MEcPP signaling integrates the biosynthesis and hydrolysis of IGs through induction of nitrile-specifier protein1 and reduction of the suppressor, ESM1, and production of simple nitriles as the bioactive end product. The findings position the plastidial metabolite, MEcPP, as the initiation hub, transducing signals to adjust the activity of hard-wired gene circuitry to expand phytochemical diversity and alter the associated subcellular structure required for functionality of the secondary metabolites, thereby tailoring plant stress responses.

Plastic Transcriptomes Stabilize Immunity to Pathogen Diversity: The Jasmonic Acid and Salicylic Acid Networks within the Arabidopsis/Botrytis Pathosystem

To respond to pathogen attack, selection and associated evolution has led to the creation of plant immune system that are a highly effective and inducible defense system. Central to this system are the plant defense hormones jasmonic acid (JA) and salicylic acid (SA) and crosstalk between the two, which may play an important role in defense responses to specific pathogens or even genotypes. Here, we used the Arabidopsis thaliana-Botrytis cinerea pathosystem to test how the host's defense system functions against genetic variation in a pathogen. We measured defense-related phenotypes and transcriptomic responses in Arabidopsis wild-type Col-0 and JA- and SA-signaling mutants, coi1-1 and npr1-1, individually challenged with 96 diverse B. cinerea isolates. Those data showed genetic variation in the pathogen influences on all components within the plant defense system at the transcriptional level. We identified four gene coexpression networks and two vectors of defense variation triggered by genetic variation in B. cinerea This showed that the JA and SA signaling pathways functioned to constrain/canalize the range of virulence in the pathogen population, but the underlying transcriptomic response was highly plastic. These data showed that plants utilize major defense hormone pathways to buffer disease resistance, but not the metabolic or transcriptional responses to genetic variation within a pathogen.

Epistasis × environment interactions among Arabidopsis thaliana glucosinolate genes impact complex traits and fitness in the field

Despite the growing number of studies showing that genotype × environment and epistatic interactions control fitness, the influences of epistasis × environment interactions on adaptive trait evolution remain largely uncharacterized. Across three field trials, we quantified aliphatic glucosinolate (GSL) defense chemistry, leaf damage, and relative fitness using mutant lines of Arabidopsis thaliana varying at pairs of causal aliphatic GSL defense genes to test the impact of epistatic and epistasis × environment interactions on adaptive trait variation. We found that aliphatic GSL accumulation was primarily influenced by additive and epistatic genetic variation, leaf damage was primarily influenced by environmental variation and relative fitness was primarily influenced by epistasis and epistasis × environment interactions. Epistasis × environment interactions accounted for up to 48% of the relative fitness variation in the field. At a single field site, the impact of epistasis on relative fitness varied significantly over 2 yr, showing that epistasis × environment interactions within a location can be temporally dynamic. These results suggest that the environmental dependency of epistasis can profoundly influence the response to selection, shaping the adaptive trajectories of natural populations in complex ways, and deserves further consideration in future evolutionary studies.

A novel Filamentous Flower mutant suppresses brevipedicellus developmental defects and modulates glucosinolate and auxin levels

BREVIPEDICELLUS (BP) encodes a class-I KNOTTED1-like homeobox (KNOX) transcription factor that plays a critical role in conditioning a replication competent state in the apical meristem, and it also governs growth and cellular differentiation in internodes and pedicels. To search for factors that modify BP signaling, we conducted a suppressor screen on bp er (erecta) plants and identified a mutant that ameliorates many of the pleiotropic defects of the parent line. Map based cloning and complementation studies revealed that the defect lies in the FILAMENTOUS FLOWER (FIL) gene, a member of the YABBY family of transcriptional regulators that contribute to meristem organization and function, phyllotaxy, leaf and floral organ growth and polarity, and are also known to repress KNOX gene expression. Genetic and cytological analyses of the fil-10 suppressor line indicate that the role of FIL in promoting growth is independent of its previously characterized influences on meristem identity and lateral organ polarity, and likely occurs non-cell-autonomously from superior floral organs. Transcription profiling of inflorescences revealed that FIL downregulates numerous transcription factors which in turn may subordinately regulate inflorescence architecture. In addition, FIL, directly or indirectly, activates over a dozen genes involved in glucosinolate production in part by activating MYB28, a known activator of many aliphatic glucosinolate biosynthesis genes. In the bp er fil-10 suppressor mutant background, enhanced expression of CYP71A13, AMIDASE1 (AMI) and NITRILASE genes suggest that auxin levels can be modulated by shunting glucosinolate metabolites into the IAA biosynthetic pathway, and increased IAA levels in the bp er fil-10 suppressor accompany enhanced internode and pedicel elongation. We propose that FIL acts to oppose KNOX1 gene function through a complex regulatory network that involves changes in secondary metabolites and auxin.

A Global Coexpression Network Approach for Connecting Genes to Specialized Metabolic Pathways in Plants

Plants produce diverse specialized metabolites (SMs), but the genes responsible for their production and regulation remain largely unknown, hindering efforts to tap plant pharmacopeia. Given that genes comprising SM pathways exhibit environmentally dependent coregulation, we hypothesized that genes within a SM pathway would form tight associations (modules) with each other in coexpression networks, facilitating their identification. To evaluate this hypothesis, we used 10 global coexpression data sets, each a meta-analysis of hundreds to thousands of experiments, across eight plant species to identify hundreds of coexpressed gene modules per data set. In support of our hypothesis, 15.3 to 52.6% of modules contained two or more known SM biosynthetic genes, and module genes were enriched in SM functions. Moreover, modules recovered many experimentally validated SM pathways, including all six known to form biosynthetic gene clusters (BGCs). In contrast, bioinformatically predicted BGCs (i.e., those lacking an associated metabolite) were no more coexpressed than the null distribution for neighboring genes. These results suggest that most predicted plant BGCs are not genuine SM pathways and argue that BGCs are not a hallmark of plant specialized metabolism. We submit that global gene coexpression is a rich, largely untapped resource for discovering the genetic basis and architecture of plant natural products.

An integrated RNAseq-1H NMR metabolomics approach to understand soybean primary metabolism regulation in response to Rhizoctonia foliar blight disease

BACKGROUND

Rhizoctonia solani AG1-IA is a devastating phytopathogen causing Rhizoctonia foliar blight (RFB) of soybean worldwide with yield losses reaching 60%. Plant defense mechanisms are complex and information from different metabolic pathways is required to thoroughly understand plant defense regulation and function. Combining information from different "omics" levels such as transcriptomics, metabolomics, and proteomics is required to gain insights into plant metabolism and its regulation. As such, we studied fluctuations in soybean metabolism in response to R. solani infection at early and late disease stages using an integrated transcriptomics-metabolomics approach, focusing on the regulation of soybean primary metabolism and oxidative stress tolerance.

RESULTS

Transcriptomics (RNAseq) and metabolomics (1H NMR) data were analyzed individually and by integration using bidirectional orthogonal projections to latent structures (O2PLS) to reveal possible links between the metabolome and transcriptome during early and late infection stages. O2PLS analysis detected 516 significant transcripts, double that reported in the univariate analysis, and more significant metabolites than detected in partial least squares discriminant analysis. Strong separation of treatments based on integration of the metabolomes and transcriptomes of the analyzed soybean leaves was revealed, similar trends as those seen in analyses done on individual datasets, validating the integration method being applied. Strong fluctuations of soybean primary metabolism occurred in glycolysis, the TCA cycle, photosynthesis and photosynthates in response to R. solani infection. Data were validated using quantitative real-time PCR on a set of specific markers as well as randomly selected genes. Significant increases in transcript and metabolite levels involved in redox reactions and ROS signaling, such as peroxidases, thiamine, tocopherol, proline, L-alanine and GABA were also recorded. Levels of ethanol increased 24 h post-infection in soybean leaves, and alcohol dehydrogenase (ADH) loss-of-function mutants of Arabidopsis thaliana had higher necrosis than wild type plants.

CONCLUSIONS

As a proof-of-concept, this study offers novel insights into the biological correlations and identification of candidate genes and metabolites that can be used in soybean breeding for resistance to R. solani AG1-IA infection. Additionally, these findings imply that alcohol and its associated gene product ADH may have important roles in plant resistance to R. solani AG1-IA causing foliar blight.