Overexpression of arogenate dehydratase reveals an upstream point of metabolic control in phenylalanine biosynthesis.

Out of the three aromatic amino acids, the highest flux in plants is directed towards phenylalanine, which is utilized to synthesize proteins and thousands of phenolic metabolites contributing to plant fitness. Phenylalanine is produced predominantly in plastids via the shikimate pathway and subsequent arogenate pathway, both of which are subject to complex transcriptional and post-transcriptional regulation. Previously, it was shown that allosteric feedback inhibition of arogenate dehydratase (ADT), which catalyzes the final step of the arogenate pathway, restricts flux through phenylalanine biosynthesis. Here, we show that in petunia (Petunia hybrida) flowers, which typically produce high phenylalanine levels, ADT regulation is relaxed, but not eliminated. Moderate expression of a feedback-insensitive ADT increased flux towards phenylalanine, while high overexpression paradoxically reduced phenylalanine formation. This reduction could be partially, but not fully, recovered by bypassing other known metabolic flux control points in the aromatic amino acid network. Using comparative transcriptomics, reverse genetics, and metabolic flux analysis, we discovered that transcriptional regulation of the d-ribulose-5-phosphate 3-epimerase gene in the pentose phosphate pathway controls flux into the shikimate pathway. Taken together, our findings reveal that regulation within and upstream of the shikimate pathway shares control over phenylalanine biosynthesis in the plant cell.

Seedling growth and fall armyworm feeding preference influenced by dhurrin production in sorghum.

Cyanogenic glucosides (CGs) play a key role in host-plant defense to insect feeding; however, the metabolic tradeoffs between synthesis of CGs and plant growth are not well understood. In this study, genetic mutants coupled with nondestructive phenotyping techniques were used to study the impact of the CG dhurrin on fall armyworm [Spodoptera frugiperda (J.E. Smith)] (FAW) feeding and plant growth in sorghum [Sorghum bicolor (L.) Moench]. A genetic mutation in CYP79A1 gene that disrupts dhurrin biosynthesis was used to develop sets of near-isogenic lines (NILs) with contrasting dhurrin contents in the Tx623 bmr6 genetic background. The NILs were evaluated for differences in plant growth and FAW feeding damage in replicated greenhouse and field trials. Greenhouse studies showed that dhurrin-free Tx623 bmr6 cyp79a1 plants grew more quickly than wild-type plants but were more susceptible to insect feeding based on changes in green plant area (GPA), total leaf area, and total dry weight over time. The NILs exhibited similar patterns of growth in field trials with significant differences in leaf area and dry weight of dhurrin-free plants between the infested and non-infested treatments. Taken together, these studies reveal a significant metabolic tradeoff between CG biosynthesis and plant growth in sorghum seedlings. Disruption of dhurrin biosynthesis produces plants with higher growth rates than wild-type plants but these plants have greater susceptibility to FAW feeding.

SDG8-Mediated Histone Methylation and RNA Processing Function in the Response to Nitrate Signaling.

Chromatin modification has gained increased attention for its role in the regulation of plant responses to environmental changes, but the specific mechanisms and molecular players remain elusive. Here, we show that the Arabidopsis (Arabidopsis thaliana) histone methyltransferase SET DOMAIN GROUP8 (SDG8) mediates genome-wide changes in H3K36 methylation at specific genomic loci functionally relevant to nitrate treatments. Moreover, we show that the specific H3K36 methyltransferase encoded by SDG8 is required for canonical RNA processing, and that RNA isoform switching is more prominent in the sdg8-5 deletion mutant than in the wild type. To demonstrate that SDG8-mediated regulation of RNA isoform expression is functionally relevant, we examined a putative regulatory gene, CONSTANS, CO-like, and TOC1 101 (CCT101), whose nitrogen-responsive isoform-specific RNA expression is mediated by SDG8. We show by functional expression in shoot cells that the different RNA isoforms of CCT101 encode distinct regulatory proteins with different effects on genome-wide transcription. We conclude that SDG8 is involved in plant responses to environmental nitrogen supply, affecting multiple gene regulatory processes including genome-wide histone modification, transcriptional regulation, and RNA processing, and thereby mediating developmental and metabolic processes related to nitrogen use.

Convergent evolution of plant specialized 1,4-naphthoquinones: metabolism, trafficking, and resistance to their allelopathic effects.

Plant 1,4-naphthoquinones encompass a class of specialized metabolites known to mediate numerous plant-biotic interactions. This class of compounds also presents a remarkable case of convergent evolution. The 1,4-naphthoquinones are synthesized by species belonging to nearly 20 disparate orders spread throughout vascular plants, and their production occurs via one of four known biochemically distinct pathways. Recent developments from large-scale biology and genetic studies corroborate the existence of multiple pathways to synthesize plant 1,4-naphthoquinones and indicate that extraordinary events of metabolic innovation and links to respiratory and photosynthetic quinone metabolism probably contributed to their independent evolution. Moreover, because many 1,4-naphthoquinones are excreted into the rhizosphere and they are highly reactive in biological systems, plants that synthesize these compounds also needed to independently evolve strategies to deploy them and to resist their effects. In this review, we highlight new progress made in understanding specialized 1,4-naphthoquinone biosynthesis and trafficking with a focus on how these discoveries have shed light on the convergent evolution and diversification of this class of compounds in plants. We also discuss how emerging themes in metabolism-based herbicide resistance may provide clues to mechanisms plants employ to tolerate allelopathic 1,4-naphthoquinones.

Evaluation of a Stable Isotope-Based Direct Quantification Method for Dicamba Analysis from Air and Water Using Single-Quadrupole LC-MS.

Dicamba is a moderately volatile herbicide used for post-emergent control of broadleaf weeds in corn, soybean, and a number of other crops. With increased use of dicamba due to the release of dicamba-resistant cotton and soybean varieties, growing controversy over the effects of spray drift and volatilization on non-target crops has increased the need for quantifying dicamba collected from water and air sampling. Therefore, this study was designed to evaluate stable isotope-based direct quantification of dicamba from air and water samples using single-quadrupole liquid chromatography-mass spectrometry (LC-MS). The sample preparation protocols developed in this study utilize a simple solid-phase extraction (SPE) protocol for water samples and a single-step concentration protocol for air samples. The LC-MS detection method achieves sensitive detection of dicamba based on selected ion monitoring (SIM) of precursor and fragment ions and relies on the use of an isotopically labeled internal standard (IS) (D3-dicamba), which allows for calculating recoveries and quantification using a relative response factor (RRF). Analyte recoveries of 106-128% from water and 88-124% from air were attained, with limits of detection (LODs) of 0.1 ng mL-1 and 1 ng mL-1, respectively. The LC-MS detection method does not require sample pretreatment such as ion-pairing or derivatization to achieve sensitivity. Moreover, this study reveals matrix effects associated with sorbent resin used in air sample collection and demonstrates how the use of an isotopically labeled IS with RRF-based analysis can account for ion suppression. The LC-MS method is easily transferrable and offers a robust alternative to methods relying on more expensive tandem LC-MS/MS-based options.

A peroxisomal thioesterase plays auxiliary roles in plant ß-oxidative benzoic acid metabolism.

Peroxisomal ß-oxidative degradation of compounds is a common metabolic process in eukaryotes. Reported benzoyl-coenzyme A (BA-CoA) thioesterase activity in peroxisomes from petunia flowers suggests that, like mammals and fungi, plants contain auxiliary enzymes mediating ß-oxidation. Here we report the identification of Petunia hybrida thioesterase 1 (PhTE1), which catalyzes the hydrolysis of aromatic acyl-CoAs to their corresponding acids in peroxisomes. PhTE1 expression is spatially, developmentally and temporally regulated and exhibits a similar pattern to known benzenoid metabolic genes. PhTE1 activity is inhibited by free coenzyme A (CoA), indicating that PhTE1 is regulated by the peroxisomal CoA pool. PhTE1 downregulation in petunia flowers led to accumulation of BA-CoA with increased production of benzylbenzoate and phenylethylbenzoate, two compounds which rely on the presence of BA-CoA precursor in the cytoplasm, suggesting that acyl-CoAs can be exported from peroxisomes. Furthermore, PhTE1 downregulation resulted in increased pools of cytoplasmic phenylpropanoid pathway intermediates, volatile phenylpropenes, lignin and anthocyanins. These results indicate that PhTE1 influences (i) intraperoxisomal acyl-CoA/CoA levels needed to carry out ß-oxidation, (ii) efflux of ß-oxidative products, acyl-CoAs and free acids, from peroxisomes, and (iii) flux distribution within the benzenoid/phenylpropanoid metabolic network. Thus, this demonstrates that plant thioesterases play multiple auxiliary roles in peroxisomal ß-oxidative metabolism.

The Origin and Biosynthesis of the Benzenoid Moiety of Ubiquinone (Coenzyme Q) in Arabidopsis.

It is not known how plants make the benzenoid ring of ubiquinone, a vital respiratory cofactor. Here, we demonstrate that Arabidopsis thaliana uses for that purpose two separate biosynthetic branches stemming from phenylalanine and tyrosine. Gene network modeling and characterization of T-DNA mutants indicated that acyl-activating enzyme encoded by At4g19010 contributes to the biosynthesis of ubiquinone specifically from phenylalanine. CoA ligase assays verified that At4g19010 prefers para-coumarate, ferulate, and caffeate as substrates. Feeding experiments demonstrated that the at4g19010 knockout cannot use para-coumarate for ubiquinone biosynthesis and that the supply of 4-hydroxybenzoate, the side-chain shortened version of para-coumarate, can bypass this blockage. Furthermore, a trans-cinnamate 4-hydroxylase mutant, which is impaired in the conversion of trans-cinnamate into para-coumarate, displayed similar defects in ubiquinone biosynthesis to that of the at4g19010 knockout. Green fluorescent protein fusion experiments demonstrated that At4g19010 occurs in peroxisomes, resulting in an elaborate biosynthetic architecture where phenylpropanoid intermediates have to be transported from the cytosol to peroxisomes and then to mitochondria where ubiquinone is assembled. Collectively, these results demonstrate that At4g19010 activates the propyl side chain of para-coumarate for its subsequent ß-oxidative shortening. Evidence is shown that the peroxisomal ABCD transporter (PXA1) plays a critical role in this branch.

Contribution of CoA ligases to benzenoid biosynthesis in petunia flowers.

Biosynthesis of benzoic acid from Phe requires shortening of the side chain by two carbons, which can occur via the ß-oxidative or nonoxidative pathways. The first step in the ß-oxidative pathway is cinnamoyl-CoA formation, likely catalyzed by a member of the 4-coumarate:CoA ligase (4CL) family that converts a range of trans-cinnamic acid derivatives into the corresponding CoA thioesters. Using a functional genomics approach, we identified two potential CoA-ligases from petunia (Petunia hybrida) petal-specific cDNA libraries. The cognate proteins share only 25% amino acid identity and are highly expressed in petunia corollas. Biochemical characterization of the recombinant proteins revealed that one of these proteins (Ph-4CL1) has broad substrate specificity and represents a bona fide 4CL, whereas the other is a cinnamate:CoA ligase (Ph-CNL). RNA interference suppression of Ph-4CL1 did not affect the petunia benzenoid scent profile, whereas downregulation of Ph-CNL resulted in a decrease in emission of benzylbenzoate, phenylethylbenzoate, and methylbenzoate. Green fluorescent protein localization studies revealed that the Ph-4CL1 protein is localized in the cytosol, whereas Ph-CNL is in peroxisomes. Our results indicate that subcellular compartmentalization of enzymes affects their involvement in the benzenoid network and provide evidence that cinnamoyl-CoA formation by Ph-CNL in the peroxisomes is the committed step in the ß-oxidative pathway.

Completion of the core ß-oxidative pathway of benzoic acid biosynthesis in plants.

Despite the importance of benzoic acid (BA) as a precursor for a wide array of primary and secondary metabolites, its biosynthesis in plants has not been fully elucidated. BA formation from phenylalanine requires shortening of the C(3) side chain by two carbon units, which can occur by a non-ß-oxidative route and/or a ß-oxidative pathway analogous to the catabolism of fatty acids. Enzymes responsible for the first and last reactions of the core BA ß-oxidative pathway (cinnamic acid → cinnamoyl-CoA → 3-hydroxy-3-phenylpropanoyl-CoA → 3-oxo-3-phenylpropanoyl-CoA → BA-CoA) have previously been characterized in petunia, a plant with flowers rich in phenylpropanoid/benzenoid volatile compounds. Using a functional genomics approach, we have identified a petunia gene encoding cinnamoyl-CoA hydratase-dehydrogenase (PhCHD), a bifunctional peroxisomal enzyme responsible for two consecutively occurring unexplored intermediate steps in the core BA ß-oxidative pathway. PhCHD spatially, developmentally, and temporally coexpresses with known genes in the BA ß-oxidative pathway, and correlates with emission of benzenoid volatiles. Kinetic analysis of recombinant PhCHD revealed it most efficiently converts cinnamoyl-CoA to 3-oxo-3-phenylpropanoyl-CoA, thus forming the substrate for the final step in the pathway. Down-regulation of PhCHD expression in petunia flowers resulted in reduced CHD enzyme activity, as well as decreased formation of BA-CoA, BA and their derived volatiles. Moreover, transgenic lines accumulated the PhCHD substrate cinnamoyl-CoA and the upstream pathway intermediate cinnamic acid. Discovery of PhCHD completes the elucidation of the core BA ß-oxidative route in plants, and together with the previously characterized CoA-ligase and thiolase enzymes, provides evidence that the whole pathway occurs in peroxisomes.

MutS HOMOLOG1 is a nucleoid protein that alters mitochondrial and plastid properties and plant response to high light.

Mitochondrial-plastid interdependence within the plant cell is presumed to be essential, but measurable demonstration of this intimate interaction is difficult. At the level of cellular metabolism, several biosynthetic pathways involve both mitochondrial- and plastid-localized steps. However, at an environmental response level, it is not clear how the two organelles intersect in programmed cellular responses. Here, we provide evidence, using genetic perturbation of the MutS Homolog1 (MSH1) nuclear gene in five plant species, that MSH1 functions within the mitochondrion and plastid to influence organellar genome behavior and plant growth patterns. The mitochondrial form of the protein participates in DNA recombination surveillance, with disruption of the gene resulting in enhanced mitochondrial genome recombination at numerous repeated sequences. The plastid-localized form of the protein interacts with the plastid genome and influences genome stability and plastid development, with its disruption leading to variegation of the plant. These developmental changes include altered patterns of nuclear gene expression. Consistency of plastid and mitochondrial response across both monocot and dicot species indicate that the dual-functioning nature of MSH1 is well conserved. Variegated tissues show changes in redox status together with enhanced plant survival and reproduction under photooxidative light conditions, evidence that the plastid changes triggered in this study comprise an adaptive response to naturally occurring light stress.

Transcriptional repression of GTL1 under water-deficit stress promotes anthocyanin biosynthesis to enhance drought tolerance.

The transcription factor GT2-LIKE 1 (GTL1) has been implicated in orchestrating a transcriptional network of diverse physiological, biochemical, and developmental processes. In response to water-limiting conditions, GTL1 is a negative regulator of stomatal development, but its potential rolein other water-deficit responses is unknown. We hypothesized that GTL1 regulates transcriptome changes associated with drought tolerance over leaf developmental stages. To test the hypothesis, gene expression was profiled by RNA-seq analysis in emerging and expanding leaves of wild-type and a drought-tolerant gtl1-4 knockout mutant under well-watered and water-deficit conditions. Our comparative analysis of genotype-treatment combinations within leaf developmental age identified 459 and 1073 differentially expressed genes in emerging and expanding leaves, respectively, as water-deficit responsive GTL1-regulated genes. Transcriptional profiling identified a potential role of GTL1 in two important pathways previously linked to drought tolerance: flavonoid and polyamine biosynthesis. In expanding leaves, negative regulation of GTL1 under water-deficit conditions promotes biosynthesis of flavonoids and anthocyanins that may contribute to drought tolerance. Quantification of polyamines did not support a role for GTL1 in these drought-responsive pathways, but this is likely due to the complex nature of polyamine synthesis and turnover. Our global transcriptome analysis suggests that transcriptional repression of GTL1 by water deficit allows plants to activate diverse pathways that collectively contribute to drought tolerance.

Phylloquinone (vitamin K(1) ) biosynthesis in plants: two peroxisomal thioesterases of Lactobacillales origin hydrolyze 1,4-dihydroxy-2-naphthoyl-CoA.

It is not known how plants cleave the thioester bond of 1,4-dihydroxy-2-naphthoyl-CoA (DHNA-CoA), a necessary step to form the naphthoquinone ring of phylloquinone (vitamin K(1) ). In fact, only recently has the hydrolysis of DHNA-CoA been demonstrated to be enzyme driven in vivo, and the cognate thioesterase characterized in the cyanobacterium Synechocystis. With a few exceptions in certain prokaryotic (Sorangium and Opitutus) and eukaryotic (Cyanidium, Cyanidioschyzon and Paulinella) organisms, orthologs of DHNA-CoA thioesterase are missing outside of the cyanobacterial lineage. In this study, genomic approaches and functional complementation experiments identified two Arabidopsis genes encoding functional DHNA-CoA thioesterases. The deduced plant proteins display low percentages of identity with cyanobacterial DHNA-CoA thioesterases, and do not even share the same catalytic motif. GFP-fusion experiments demonstrated that the Arabidopsis proteins are targeted to peroxisomes, and subcellular fractionations of Arabidopsis leaves confirmed that DHNA-CoA thioesterase activity occurs in this organelle. In vitro assays with various aromatic and aliphatic acyl-CoA thioester substrates showed that the recombinant Arabidopsis enzymes preferentially hydrolyze DHNA-CoA. Cognate T-DNA knock-down lines display reduced DHNA-CoA thioesterase activity and phylloquinone content, establishing in vivo evidence that the Arabidopsis enzymes are involved in phylloquinone biosynthesis. Extraordinarily, structure-based phylogenies coupled to comparative genomics demonstrate that plant DHNA-CoA thioesterases originate from a horizontal gene transfer with a bacterial species of the Lactobacillales order.

Functional convergence of structurally distinct thioesterases from cyanobacteria and plants involved in phylloquinone biosynthesis.

The synthesis of phylloquinone (vitamin K1) in photosynthetic organisms requires a thioesterase that hydrolyzes 1,4-dihydroxy-2-naphthoyl-CoA (DHNA-CoA) to release 1,4-dihydroxy-2-naphthoate (DHNA). Cyanobacteria and plants contain distantly related hotdog-fold thioesterases that catalyze this reaction, although the structural basis of these convergent enzymatic activities is unknown. To investigate this, the crystal structures of hotdog-fold DHNA-CoA thioesterases from the cyanobacterium Synechocystis (Slr0204) and the flowering plant Arabidopsis thaliana (AtDHNAT1) were determined. These enzymes form distinct homotetramers and use different active sites to catalyze hydrolysis of DHNA-CoA, similar to the 4-hydroxybenzoyl-CoA (4-HBA-CoA) thioesterases from Pseudomonas and Arthrobacter. Like the 4-HBA-CoA thioesterases, the DHNA-CoA thioesterases contain either an active-site aspartate (Slr0204) or glutamate (AtDHNAT1) that are predicted to be catalytically important. Computational modeling of the substrate-bound forms of both enzymes indicates the residues that are likely to be involved in substrate binding and catalysis. Both enzymes are selective for DHNA-CoA as a substrate, but this selectivity is achieved using divergent predicted binding strategies. The Slr0204 binding pocket is predominantly hydrophobic and closely conforms to DHNA, while that of AtDHNAT1 is more polar and solvent-exposed. Considered in light of the related 4-HBA-CoA thioesterases, these structures indicate that hotdog-fold thioesterases using either an active-site aspartate or glutamate diverged into distinct clades prior to the evolution of strong substrate specificity in these enzymes.

Strategies to study the metabolic origins of specialized plant metabolites: The specialized 1,4-naphthoquinones.

One of the hallmarks of specialized plant metabolites is that they are produced using precursors from central metabolism. Therefore, in addition to identifying and characterizing the pathway genes and enzymes involved in synthesizing a specialized compound, it is critical to study its metabolic origins. Identifying what primary metabolic pathways supply precursors to specialized metabolism and how primary metabolism has diversified to sustain fluxes to specialized metabolite pathways is imperative to optimizing synthetic biology strategies for producing high-value plant natural products in crops and microbial systems. Improved understanding of the metabolic origins of specialized plant metabolites has also revealed instances of recurrent evolution of the same compound, or nearly identical compounds, with similar ecological functions, thereby expanding knowledge about the factors driving the chemical diversity in the plant kingdom. In this chapter, we describe detailed methods for performing tracer studies, chemical inhibitor experiments, and reverse genetics. We use examples from investigations of the metabolic origins of specialized plant 1,4-naphthoquinones (1,4-NQs). The plant 1,4-NQs provide an excellent case study for illustrating the importance of investigating the metabolic origins of specialized metabolites. Over half a century of research by many groups has revealed that the pathways to synthesize plant 1,4-NQs are the result of multiple events of convergent evolution across several disparate plant lineages and that plant 1,4-NQ pathways are supported by extraordinary events of metabolic innovation and by various primary metabolic sources.

A reference genome for the long-term kleptoplast-retaining sea slug Elysia crispata morphotype clarki.

Several species of sacoglossan sea slugs possess the incredible ability to sequester chloroplasts from the algae they consume. These "photosynthetic animals" incorporate stolen chloroplasts, called kleptoplasts, into the epithelial cells of tubules that extend from their digestive tracts throughout their bodies. The mechanism by which these slugs maintain functioning kleptoplasts in the absence of an algal nuclear genome is unknown. Here, we report a draft genome of the sacoglossan slug Elysia crispata morphotype clarki, a morphotype native to the Florida Keys that can retain photosynthetically active kleptoplasts for several months without feeding. We used a combination of Oxford Nanopore Technologies long reads and Illumina short reads to produce a 786-Mb assembly (N50 = 0.459 Mb) containing 68,514 predicted protein-coding genes. A phylogenetic analysis found no evidence of horizontal acquisition of genes from algae. We performed gene family and gene expression analyses to identify E. crispata genes unique to kleptoplast-containing slugs that were more highly expressed in fed versus unfed developmental life stages. Consistent with analyses in other kleptoplastic slugs, our investigation suggests that genes encoding lectin carbohydrate-binding proteins and those involved in regulation of reactive oxygen species and immunity may play a role in kleptoplast retention. Lastly, we identified four polyketide synthase genes that could potentially encode proteins producing UV- and oxidation-blocking compounds in slug cell membranes. The genome of E. crispata is a quality resource that provides potential targets for functional analyses and enables further investigation into the evolution and mechanisms of kleptoplasty in animals.

A dedicated thioesterase of the Hotdog-fold family is required for the biosynthesis of the naphthoquinone ring of vitamin K1.

Phylloquinone (vitamin K(1)) is a bipartite molecule that consists of a naphthoquinone ring attached to a phytyl side chain. The coupling of these 2 moieties depends on the hydrolysis of the CoA thioester of 1,4-dihydroxy-2-naphthoate (DHNA), which forms the naphthalenoid backbone. It is not known whether such a hydrolysis is enzymatic or chemical. In this study, comparative genomic analyses identified orthologous genes of unknown function that in most species of cyanobacteria cluster with predicted phylloquinone biosynthetic genes. The encoded approximately 16-kDa proteins display homology with some Hotdog domain-containing CoA thioesterases that are involved in the catabolism of 4-hydroxybenzoyl-CoA and gentisyl-CoA (2,5-dihydroxybenzoyl-CoA) in certain soil-dwelling bacteria. The Synechocystis ortholog, encoded by gene slr0204, was expressed as a recombinant protein and was found to form DHNA as reaction product. Unlike its homologs in the Hotdog domain family, Slr0204 showed strict substrate specificity. The Synechocystis slr0204 knockout was devoid of DHNA-CoA thioesterease activity and accumulated DHNA-CoA. As a result, knockout cells contained 13-fold less phylloquinone than their wild-type counterparts and displayed the typical photosensitivity to high light associated to phylloquinone deficiency in cyanobacteria.

Allelopathy as an evolutionary game.

In plants, most competition is resource competition, where one plant simply preempts the resources away from its neighbors. Interference competition, as the name implies, is a form of direct interference to prevent resource access. Interference competition is common among animals that can physically fight, but in plants, one of the main mechanisms of interference competition is allelopathy. Allelopathic plants release cytotoxic chemicals into the environment which can increase their ability to compete with surrounding organisms for limited resources. The circumstances and conditions favoring the development and maintenance of allelochemicals, however, are not well understood. Particularly, despite the obvious benefits of allelopathy, current data suggest it seems to have only rarely evolved. To gain insight into the cost and benefit of allelopathy, we have developed a 2 × 2 matrix game to model the interaction between plants that produce allelochemicals and plants that do not. Production of an allelochemical introduces novel cost associated with both synthesis and detoxifying a toxic chemical but may also convey a competitive advantage. A plant that does not produce an allelochemical will suffer the cost of encountering one. Our model predicts three cases in which the evolutionarily stable strategies are different. In the first, the nonallelopathic plant is a stronger competitor, and not producing allelochemicals is the evolutionarily stable strategy. In the second, the allelopathic plant is the better competitor, and production of allelochemicals is the more beneficial strategy. In the last case, neither is the evolutionarily stable strategy. Instead, there are alternating stable states, depending on whether the allelopathic or nonallelopathic plant arrived first. The generated model reveals circumstances leading to the evolution of allelochemicals and sheds light on utilizing allelochemicals as part of weed management strategies. In particular, the wide region of alternative stable states in most parameterizations, combined with the fact that the absence of allelopathy is likely the ancestral state, provides an elegant answer to the question of why allelopathy seems to rarely evolve despite its obvious benefits. Allelopathic plants can indeed outcompete nonallelopathic plants, but this benefit is simply not great enough to allow them to go to fixation and spread through the population. Thus, most populations would remain purely nonallelopathic.

Integrative analysis of the shikonin metabolic network identifies new gene connections and reveals evolutionary insight into shikonin biosynthesis.

Plant specialized 1,4-naphthoquinones present a remarkable case of convergent evolution. Species across multiple discrete orders of vascular plants produce diverse 1,4-naphthoquinones via one of several pathways using different metabolic precursors. Evolution of these pathways was preceded by events of metabolic innovation and many appear to share connections with biosynthesis of photosynthetic or respiratory quinones. Here, we sought to shed light on the metabolic connections linking shikonin biosynthesis with its precursor pathways and on the origins of shiknoin metabolic genes. Downregulation of Lithospermum erythrorhizon geranyl diphosphate synthase (LeGPPS), recently shown to have been recruited from a cytoplasmic farnesyl diphosphate synthase (FPPS), resulted in reduced shikonin production and a decrease in expression of mevalonic acid and phenylpropanoid pathway genes. Next, we used LeGPPS and other known shikonin pathway genes to build a coexpression network model for identifying new gene connections to shikonin metabolism. Integrative in silico analyses of network genes revealed candidates for biochemical steps in the shikonin pathway arising from Boraginales-specific gene family expansion. Multiple genes in the shikonin coexpression network were also discovered to have originated from duplication of ubiquinone pathway genes. Taken together, our study provides evidence for transcriptional crosstalk between shikonin biosynthesis and its precursor pathways, identifies several shikonin pathway gene candidates and their evolutionary histories, and establishes additional evolutionary links between shikonin and ubiquinone metabolism. Moreover, we demonstrate that global coexpression analysis using limited transcriptomic data obtained from targeted experiments is effective for identifying gene connections within a defined metabolic network.

A bimodular oxidoreductase mediates the specific reduction of phylloquinone (vitamin K1) in chloroplasts.

Plants and certain species of cyanobacteria are the only organisms capable of synthesizing phylloquinone (vitamin K1 for vertebrates), which they use as an electron carrier during photosynthesis. Recent studies, however, have identified a plastidial pool of non-photoactive phylloquinone that could be involved in additional cellular functions. Here, we characterized an Arabidopsis bimodular enzyme--the At4g35760 gene product--comprising an integral domain homologous to the catalytic subunit of mammalian vitamin K1 epoxide reductase (VKORC1, EC 1.1.4.1) that is fused to a soluble thioredoxin-like moiety. GFP-fusion experiments in tobacco mesophyll cells established that the plant protein is targeted to plastids, and analyses of transcript and protein levels showed that expression is maximal in leaf tissues. The fused and individual VKORC1 domains were separately expressed in yeast, removing their chloroplast targeting pre-sequence and adding a C-terminal consensus signal for retention in the endoplasmic reticulum. The corresponding microsomal preparations were equally effective at mediating the dithiotreitol-dependent reduction of phylloquinone and menaquinone into their respective quinol forms. Strikingly, unlike mammalian VKORC1, the Arabidopsis enzyme did not reduce phylloquinone epoxide, and was resistant to inhibition by warfarin. The isoprenoid benzoquinone conjugates plastoquinone and ubiquinone were not substrates, establishing that the plant enzyme evolved strict specificity for the quinone form of naphthalenoid conjugates. In vitro reconstitution experiments established that the soluble thioredoxin-like domain can function as an electron donor for its integral VKORC1 partner.