Combinatorial Control of Plant Specialized Metabolism: Mechanisms, Functions, and Consequences.

Plants constantly perceive internal and external cues, many of which they need to address to safeguard their proper development and survival. They respond to these cues by selective activation of specific metabolic pathways involving a plethora of molecular players that act and interact in complex networks. In this review, we illustrate and discuss the complexity in the combinatorial control of plant specialized metabolism. We hereby go beyond the intuitive concept of combinatorial control as exerted by modular-acting complexes of transcription factors that govern expression of specialized metabolism genes. To extend this discussion, we also consider all known hierarchical levels of regulation of plant specialized metabolism and their interfaces by referring to reported regulatory concepts from the plant field. Finally, we speculate on possible yet-to-be-discovered regulatory principles of plant specialized metabolism that are inspired by knowledge from other kingdoms of life and areas of biological research.

The ALOG family members OsG1L1 and OsG1L2 regulate inflorescence branching in rice.

The architecture of the rice inflorescence is an important determinant of crop yield. The length of the inflorescence and the number of branches are among the key factors determining the number of spikelets, and thus grains, that a plant will develop. In particular, the timing of the identity transition from indeterminate branch meristem to determinate spikelet meristem governs the complexity of the inflorescence. In this context, the ALOG gene TAWAWA1 (TAW1) has been shown to delay the transition to determinate spikelet development in Oryza sativa (rice). Recently, by combining precise laser microdissection of inflorescence meristems with RNA-seq, we observed that two ALOG genes, OsG1-like 1 (OsG1L1) and OsG1L2, have expression profiles similar to that of TAW1. Here, we report that osg1l1 and osg1l2 loss-of-function CRISPR mutants have similar phenotypes to the phenotype of the previously published taw1 mutant, suggesting that these genes might act on related pathways during inflorescence development. Transcriptome analysis of the osg1l2 mutant suggested interactions of OsG1L2 with other known inflorescence architecture regulators and the data sets were used for the construction of a gene regulatory network (GRN), proposing interactions among genes potentially involved in controlling inflorescence development in rice. In this GRN, we selected the homeodomain-leucine zipper transcription factor encoding the gene OsHOX14 for further characterization. The spatiotemporal expression profiling and phenotypical analysis of CRISPR loss-of-function mutants of OsHOX14 suggests that the proposed GRN indeed serves as a valuable resource for the identification of new proteins involved in rice inflorescence development.

The saponin bomb: a nucleolar-localized ß-glucosidase hydrolyzes triterpene saponins in Medicago truncatula.

Plants often protect themselves from their own bioactive defense metabolites by storing them in less active forms. Consequently, plants also need systems allowing correct spatiotemporal reactivation of such metabolites, for instance under pathogen or herbivore attack. Via co-expression analysis with public transcriptomes, we determined that the model legume Medicago truncatula has evolved a two-component system composed of a ß-glucosidase, denominated G1, and triterpene saponins, which are physically separated from each other in intact cells. G1 expression is root-specific, stress-inducible, and coregulated with that of the genes encoding the triterpene saponin biosynthetic enzymes. However, the G1 protein is stored in the nucleolus and is released and united with its typically vacuolar-stored substrates only upon tissue damage, partly mediated by the surfactant action of the saponins themselves. Subsequently, enzymatic removal of carbohydrate groups from the saponins creates a pool of metabolites with an increased broad-spectrum antimicrobial activity. The evolution of this defense system benefited from both the intrinsic condensation abilities of the enzyme and the bioactivity properties of its substrates. We dub this two-component system the saponin bomb, in analogy with the mustard oil and cyanide bombs, commonly used to describe the renowned ß-glucosidase-dependent defense systems for glucosinolates and cyanogenic glucosides.

A Seed-Specific Regulator of Triterpene Saponin Biosynthesis in Medicago truncatula.

Plants produce a vast array of defense compounds to protect themselves from pathogen attack or herbivore predation. Saponins are a specific class of defense compounds comprising bioactive glycosides with a steroidal or triterpenoid aglycone backbone. The model legume Medicago truncatula synthesizes two types of saponins, hemolytic saponins and nonhemolytic soyasaponins, which accumulate as specific blends in different plant organs. Here, we report the identification of the seed-specific transcription factor TRITERPENE SAPONIN ACTIVATION REGULATOR3 (TSAR3), which controls hemolytic saponin biosynthesis in developing M. truncatula seeds. Analysis of genes that are coexpressed with TSAR3 in transcriptome data sets from developing M. truncatula seeds led to the identification of CYP88A13, a cytochrome P450 that catalyzes the C-16α hydroxylation of medicagenic acid toward zanhic acid, the final oxidation step of the hemolytic saponin biosynthesis branch in M. truncatula In addition, two uridine diphosphate glycosyltransferases, UGT73F18 and UGT73F19, which glucosylate hemolytic sapogenins at the C-3 position, were identified. The genes encoding the identified biosynthetic enzymes are present in clusters of duplicated genes in the M. truncatula genome. This appears to be a common theme among saponin biosynthesis genes, especially glycosyltransferases, and may be the driving force of the metabolic evolution of saponins.

Transcriptome analysis reveals rice MADS13 as an important repressor of the carpel development pathway in ovules.

In angiosperms, floral homeotic genes encoding MADS-domain transcription factors regulate the development of floral organs. Specifically, members of the SEPALLATA (SEP) and AGAMOUS (AG) subfamilies form higher-order protein complexes to control floral meristem determinacy and to specify the identity of female reproductive organs. In rice, the AG subfamily gene OsMADS13 is intimately involved in the determination of ovule identity, since knock-out mutant plants develop carpel-like structures in place of ovules, resulting in female sterility. Little is known about the regulatory pathways at the base of rice gynoecium development. To investigate molecular mechanisms acting downstream of OsMADS13, we obtained transcriptomes of immature inflorescences from wild-type and Osmads13 mutant plants. Among a total of 476 differentially expressed genes (DEGs), a substantial overlap with DEGs from the SEP-family Osmads1 mutant was found, suggesting that OsMADS1 and OsMADS13 may act on a common set of target genes. Expression studies and preliminary analyses of two up-regulated genes encoding Zinc-finger transcription factors indicated that our dataset represents a valuable resource for the identification of both OsMADS13 target genes and novel players in rice ovule development. Taken together, our study suggests that OsMADS13 is an important repressor of the carpel pathway during ovule development.

The MYB transcription factor Emission of Methyl Anthranilate 1 stimulates emission of methyl anthranilate from Medicago truncatula hairy roots.

Plants respond to herbivore or pathogen attacks by activating specific defense programs that include the production of bioactive specialized metabolites to eliminate or deter the attackers. Volatiles play an important role in the interaction of a plant with its environment. Through transcript profiling of jasmonate-elicited Medicago truncatula cells, we identified Emission of Methyl Anthranilate (EMA) 1, a MYB transcription factor that is involved in the emission of the volatile compound methyl anthranilate when expressed in M. truncatula hairy roots, giving them a fruity scent. RNA sequencing (RNA-Seq) analysis of the fragrant roots revealed the upregulation of a methyltransferase that was subsequently characterized to catalyze the O-methylation of anthranilic acid and was hence named M. truncatula anthranilic acid methyl transferase (MtAAMT) 1. Given that direct activation of the MtAAMT1 promoter by EMA1 could not be unambiguously demonstrated, we further probed the RNA-Seq data and identified the repressor protein M. truncatula plant AT-rich sequence and zinc-binding (MtPLATZ) 1. Emission of Methyl Anthranilate 1 binds a tandem repeat of the ACCTAAC motif in the MtPLATZ1 promoter to transactivate gene expression. Overexpression of MtPLATZ1 in transgenic M. truncatula hairy roots led to transcriptional silencing of EMA1, indicating that MtPLATZ1 may be part of a negative feedback loop to control the expression of EMA1. Finally, application of exogenous methyl anthranilate boosted EMA1 and MtAAMT1 expression dramatically, thus also revealing a positive amplification loop. Such positive and negative feedback loops seem to be the norm rather than the exception in the regulation of plant specialized metabolism.

Clustered regularly interspaced short palindromic repeats tools for plant metabolic engineering: achievements and perspectives.

The plant kingdom represents the biggest source of feedstock, food, and added-value compounds. Engineering plant metabolic pathways to increase the phytochemical production or improve the nutraceutical value of crops is challenging because of the intricate interaction networks that link multiple genes, enzymatic steps, and metabolites, even when pathways are fully elucidated. The development of clustered regularly interspaced short palindromic repeats - CRISPR-associated (CRISPR-Cas) technologies has helped to overcome limitations in metabolic engineering, providing efficient and versatile tools for multigene editing. CRISPR approaches in plants were shown to have a remarkable efficiency in genome editing of different species to improve agronomic and metabolic traits. Here, we give an overview of the different achievements and perspectives of CRISPR technology in plant metabolic engineering.

Structural and functional diversity in plant specialized metabolism signals and products: The case of oxylipins and triterpenes.

Metabolic enzymes tend to evolve towards catalytic efficacy, precision and speed. This seems particularly true for ancient and conserved enzymes involved in fundamental cellular processes that are present virtually in every cell and organism and converting and producing relatively limited metabolite numbers. Nevertheless, sessile organisms like plants have an astonishing repertoire of specific (specialized) metabolites that, by numbers and chemical complexity, by far exceed primary metabolites. Most theories agree that early gene duplication, subsequent positive selection and diversifying evolution have allowed relaxed selection of duplicated metabolic genes, thus facilitating the accumulation of mutations that could broaden substrate/product specificity and lower activation barriers and kinetics. Here, we use oxylipins, oxygenated fatty acids of plastidial origin to which the phytohormone jasmonate belongs, and triterpenes, a large group of specialized metabolites whose biosynthesis is often elicited by jasmonates, to showcase the structural and functional diversity of chemical signals and products in plant metabolism.

Identification and characterization of CYP71 subclade cytochrome P450 enzymes involved in the biosynthesis of bitterness compounds in Cichorium intybus.

Industrial chicory (Cichorium intybus var. sativum) and witloof (C. intybus var. foliosum) are crops with an important economic value, mainly cultivated for inulin production and as a leafy vegetable, respectively. Both crops are rich in nutritionally relevant specialized metabolites with beneficial effects for human health. However, their bitter taste, caused by the sesquiterpene lactones (SLs) produced in leaves and taproot, limits wider applications in the food industry. Changing the bitterness would thus create new opportunities with a great economic impact. Known genes encoding enzymes involved in the SL biosynthetic pathway are GERMACRENE A SYNTHASE (GAS), GERMACRENE A OXIDASE (GAO), COSTUNOLIDE SYNTHASE (COS) and KAUNIOLIDE SYNTHASE (KLS). In this study, we integrated genome and transcriptome mining to further unravel SL biosynthesis. We found that C. intybus SL biosynthesis is controlled by the phytohormone methyl jasmonate (MeJA). Gene family annotation and MeJA inducibility enabled the pinpointing of candidate genes related with the SL biosynthetic pathway. We specifically focused on members of subclade CYP71 of the cytochrome P450 family. We verified the biochemical activity of 14 C. intybus CYP71 enzymes transiently produced in Nicotiana benthamiana and identified several functional paralogs for each of the GAO, COS and KLS genes, pointing to redundancy in and robustness of the SL biosynthetic pathway. Gene functionality was further analyzed using CRISPR/Cas9 genome editing in C. intybus. Metabolite profiling of mutant C. intybus lines demonstrated a successful reduction in SL metabolite production. Together, this study increases our insights into the C. intybus SL biosynthetic pathway and paves the way for the engineering of C. intybus bitterness.

Jasmonate: A hormone of primary importance for plant metabolism.

Over the years, jasmonates (JAs) have become recognized as one of the main plant hormones that regulate stress responses by activating defense programs and the production of specialized metabolites. High JA levels have been associated with reduced plant growth, supposedly as a result of the reallocation of carbon sources from primary growth to the biosynthesis of defense compounds. Recent advances suggest however that tight regulatory networks integrate several sensing pathways to steer plant metabolism, and thereby drive the trade-off between growth and defense. In this review, we discuss how JA influences primary metabolism and how it is connected to light-regulated processes, nutrient sensing and energy metabolism. Finally, we speculate that JA, in a conceptual parallelism with adrenaline for humans, overall boosts cellular processes to keep up with an increased metabolic demand during harsh times.

Interference between ER stress-related bZIP-type and jasmonate-inducible bHLH-type transcription factors in the regulation of triterpene saponin biosynthesis in Medicago truncatula.

Triterpene saponins (TS) are a structurally diverse group of metabolites that are widely distributed in plants. They primarily serve as defense compounds and their production is often triggered by biotic stresses through signaling cascades that are modulated by phytohormones such as the jasmonates (JA). Two JA-modulated basic helix-loop-helix (bHLH) transcription factors (TFs), triterpene saponin biosynthesis activating regulator 1 (TSAR1) and TSAR2, have previously been identified as direct activators of TS biosynthesis in the model legume Medicago truncatula. Here, we report on the involvement of the core endoplasmic reticulum (ER) stress-related basic leucine zipper (bZIP) TFs bZIP17 and bZIP60 in the regulation of TS biosynthesis. Expression and processing of M. truncatula bZIP17 and bZIP60 proteins were altered in roots with perturbed TS biosynthesis or treated with JA. Accordingly, such roots displayed an altered ER network structure. M. truncatula bZIP17 and bZIP60 proteins were shown to localize in the nucleus and appeared to be capable of interfering with the TSAR-mediated transactivation of TS biosynthesis genes. Furthermore, interference between ER stress-related bZIP and JA-modulated bHLH TFs in the regulation of JA-dependent terpene biosynthetic pathways may be widespread in the plant kingdom, as we demonstrate that it also occurs in the regulation of monoterpene indole alkaloid biosynthesis in the medicinal plant Catharanthus roseus.

Corrigendum: Interference between ER stress-related bZIP-type and jasmonate-inducible bHLH-type transcription factors in the regulation of triterpene saponin biosynthesis in Medicago truncatula.

[This corrects the article DOI: 10.3389/fpls.2022.903793.].

CRISPR-mediated accelerated domestication of African rice landraces.

African Oryza glaberrima and Oryza sativa landraces are considered valuable resources for breeding traits due to their adaptation to local environmental and soil conditions. They often possess superior resistance to endemic pests and tolerance to drought and nutrient deficiencies when compared to the "imported" high production Asian rice varieties. In contrast, "domestication traits" such as seed shattering, lodging, and seed yield are not well established in these African landraces. Therefore, the use of these African varieties for high production agriculture is limited by unpredictable yield and grain quality. We are addressing this shortcoming by developing protocols for genetically transforming African landraces to allow the use of CRISPR-Cas mediated breeding approaches. Here we use as proof of concept the cultivated African landrace Kabre to target selected known "domestication loci" and improve the agronomic potential of Kabre rice. Stable genetic transformation with CRISPR-Cas9-based vectors generated single and simultaneous multiple gene knockouts. Plants with reduced stature to diminish lodging were generated by disrupting the HTD1 gene. Furthermore, three loci shown to control seed size and/or yield (GS3, GW2 and GN1A) were targeted using a multiplex CRISPR-Cas9 construct. This resulted in mutants with significantly improved seed yield. Our study provides an example of how new breeding technologies can accelerate the development of highly productive African landrace rice varieties, an important advancement considering that Africa is a hotspot for worldwide population growth and therefore prone to food shortage.

A Genomic View of Alternative Splicing of Long Non-coding RNAs during Rice Seed Development Reveals Extensive Splicing and lncRNA Gene Families.

Alternative splicing (AS) is a key modulator of development in many eukaryotic organisms. In plants, alternative splice forms of non-coding RNAs (ncRNAs) are known to modulate flowering time in Arabidopsis and fertility in rice. Here we demonstrate that alternative splicing of coding and long non-coding RNAs occurs during rice seed development by comparing AS in immature seeds vs. embryo and endosperm of mature seeds. Based on computational predictions of AS events determined from a Bayesian analysis of junction counts of RNA-seq datasets, differential splicing of protein-coding, and non-coding RNAs was determined. In contrast to roots, leaves, flowers, buds, and reproductive meristems, developing seeds had 5.8-57 times more predicted AS. Primers designed to span introns and exons were used to detect AS events predicted by rMATs in cDNA derived from early (milk) seed, embryo, and endosperm. Comparing milk seed vs. mature embryo and endosperm, AS of MORC7 (a gene implicated in epigenetic gene silencing), was markedly different. Long non-coding RNAs (lncRNAs) also underwent AS during the transition from milk seed to mature embryo and endosperm, with a complex gene structure, and were more extensively processed than predicted by current genome annotation. Exon retention of lncRNAs was enhanced in embryos. Searching all 5,515 lncRNAs in the NCBI genome annotation uncovered gene families based on highly conserved regions shared by groups of 3-35 lncRNAs. The homologies to other lncRNAs, as well as homologies to coding sequences, and the genomic context of lncRNAs provide inroads for functional analysis of multi-exonic lncRNAs that can be extensively processed during seed development.