Maize specialized metabolome networks reveal organ-preferential mixed glycosides

Functional characterization of a maize UDP-glucosyltransferase involved in the biosynthesis of flavonoid 7-O-glucosides and di-O-glucosides

Flavonoids are glycosylated in the final steps of their synthesis by UDP-dependent glycosyltransferase enzymes. We present the functional characterization of the first maize flavonoid O-glucosyltransferase enzyme from family 85, which exhibits properties not previously described. ZmUGT85W2 catalyzes the O-glucosylation of flavonols, flavones and flavanones, to form primarily 7-O-glucosides, but also flavonol O-glucoside positional isomers, flavones and flavonol di-O-glucosides. ZmUGT85W2 exhibited a differential kinetic behavior depending on the flavonoid acceptor, showing hyperbolic dependence for flavonols and sigmoidal response for flavanones and flavones. Structural and molecular docking analyses predicted conserved residues interacting with the sugar donor, with close contact with the 7-hydroxyl of the flavonoid acceptors, consistent with enzymatic activity results. In addition, ZmUGT85W2 is induced by UV-B radiation, and its expression is controlled by the B and PL1 transcription factors. Consistently, higher levels of flavone and flavonol O-glycosides are accumulated in leaves of plants exposed to solar UV-B compared to control plants, suggesting that ZmUGT85W2 is involved in the biosynthesis of these metabolites in maize leaves, contributing to UV-B tolerance. The activity of ZmUGT85W2, along with its elevated expression in silks and pericarps expressing the R2R3-MYB transcription factor P1, highlights its critical role in the accumulation of flavonoid O-glucosides in these tissues. Together, our findings reveal a key step in maize flavonoid O-glycosides biosynthesis, with the observed positive cooperative behaviors suggesting that ZmUGT85W2 plays a crucial role in finely regulating metabolic flux towards these compounds in planta.

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

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

Exploration of the Diversity of Vicine and Convicine Derivatives in Faba Bean (Vicia faba L.) Cultivars: Insights from LC-MS/MS Spectra

While numerous Fabaceae seeds are a good nutritional source of high-quality protein, the use of some species is hampered by toxic effects caused by exposure to metabolites that accumulate in the seeds. One such species is the faba or broad bean (Vicia faba L.), which accumulates vicine and convicine. These two glycoalkaloids cause favism, the breakdown of red blood cells in persons with a glucose-6-phosphate dehydrogenase deficiency. Because this is the most common enzyme deficiency worldwide, faba bean breeding efforts have focused on developing cultivars with low levels of these alkaloids. Consequently, quantification methods have been developed; however, they quantify vicine and convicine only and not the derivatives of these compounds that potentially generate the same bio-active molecules. Based on the recognition of previously unknown (con)vicine-containing compounds, we screened the fragmentation spectra of LC-MS/MS data from five faba bean cultivars using the characteristic fragments generated by (con)vicine. This resulted in the recognition of more than a hundred derivatives, of which 89 were tentatively identified. (Con)vicine was mainly derivatized through the addition of sugars, hydroxycinnamic acids, and dicarboxylic acids, with a group of compounds composed of two (con)vicine residues linked by dicarboxyl fatty acids. In general, the abundance profiles of the different derivatives in the five cultivars mimicked that of vicine and convicine, but some showed a derivative-specific profile. The description of the (con)vicine diversity will impact the interpretation of future studies on the biosynthesis of (con)vicine, and the content in potentially bio-active alkaloids in faba beans may be higher than that represented by the quantification of vicine and convicine alone.

Pour some sugar on me: The diverse functions of phenylpropanoid glycosylation

The phenylpropanoid metabolism is the source of a vast array of specialized metabolites that play diverse functions in plant growth and development and contribute to all aspects of plant interactions with their surrounding environment. These compounds protect plants from damaging ultraviolet radiation and reactive oxygen species, provide mechanical support for the plants to stand upright, and mediate plant-plant and plant-microorganism communications. The enormous metabolic diversity of phenylpropanoids is further expanded by chemical modifications known as "decorative reactions", including hydroxylation, methylation, glycosylation, and acylation. Among these modifications, glycosylation is the major driving force of phenylpropanoid structural diversification, also contributing to the expansion of their properties. Phenylpropanoid glycosylation is catalyzed by regioselective uridine diphosphate (UDP)-dependent glycosyltransferases (UGTs), whereas glycosyl hydrolases known as β-glucosidases are the major players in deglycosylation. In this article, we review how the glycosylation process affects key physicochemical properties of phenylpropanoids, such as molecular stability and solubility, as well as metabolite compartmentalization/storage and biological activity/toxicity. We also summarize the recent knowledge on the functional implications of glycosylation of different classes of phenylpropanoid compounds. A balance of glycosylation/deglycosylation might represent an essential molecular mechanism to regulate phenylpropanoid homeostasis, allowing plants to dynamically respond to diverse environmental signals.

A comprehensive metabolomics and lipidomics atlas for the legumes common bean, chickpea, lentil and lupin

Legumes represent an important component of human and livestock diets; they are rich in macro- and micronutrients such as proteins, dietary fibers and polyunsaturated fatty acids. Whilst several health-promoting and anti-nutritional properties have been associated with grain content, in-depth metabolomics characterization of major legume species remains elusive. In this article, we used both gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS) to assess the metabolic diversity in the five legume species commonly grown in Europe, including common bean (Phaseolus vulgaris), chickpea (Cicer arietinum), lentil (Lens culinaris), white lupin (Lupinus albus) and pearl lupin (Lupinus mutabilis), at the tissue level. We were able to detect and quantify over 3400 metabolites covering major nutritional and anti-nutritional compounds. Specifically, the metabolomics atlas includes 224 derivatized metabolites, 2283 specialized metabolites and 923 lipids. The data generated here will serve the community as a basis for future integration to metabolomics-assisted crop breeding and facilitate metabolite-based genome-wide association studies to dissect the genetic and biochemical bases of metabolism in legume species.

QT-GWAS: A novel method for unveiling biosynthetic loci affecting qualitative metabolic traits

Although the plant kingdom provides an enormous diversity of metabolites with potentially beneficial applications for humankind, a large fraction of these metabolites and their biosynthetic pathways remain unknown. Resolving metabolite structures and their biosynthetic pathways is key to gaining biological understanding and to allow metabolic engineering. In order to retrieve novel biosynthetic genes involved in specialized metabolism, we developed a novel untargeted method designated as qualitative trait GWAS (QT-GWAS) that subjects qualitative metabolic traits to a genome-wide association study, while the conventional metabolite GWAS (mGWAS) mainly considers the quantitative variation of metabolites. As a proof of the validity of QT-GWAS, 23 and 15 of the retrieved associations identified in Arabidopsis thaliana by QT-GWAS and mGWAS, respectively, were supported by previous research. Furthermore, seven gene-metabolite associations retrieved by QT-GWAS were confirmed in this study through reverse genetics combined with metabolomics and/or in vitro enzyme assays. As such, we established that CYTOCHROME P450 706A5 (CYP706A5) is involved in the biosynthesis of chroman derivatives, UDP-GLYCOSYLTRANSFERASE 76C3 (UGT76C3) is able to hexosylate guanine in vitro and in planta, and SULFOTRANSFERASE 202B1 (SULT202B1) catalyzes the sulfation of neolignans in vitro. Collectively, our study demonstrates that the untargeted QT-GWAS method can retrieve valid gene-metabolite associations at the level of enzyme-encoding genes, even new associations that cannot be found by the conventional mGWAS, providing a new approach for dissecting qualitative metabolic traits.

Characterization of Corn Silk Extract Using HPLC/HRMS/MS Analyses and Bioinformatic Data Processing

In addition to having different biological activities of interest, corn silks play a role in the defense of plants. While benzoxamines and flavonoids have already been identified as molecules of plant defense and growth mechanisms, knowledge on the phytochemical composition of corn silk is lacking. Such knowledge would make it possible to better select the most effective varieties to improve resistance or bioactive properties. In this article, an approach was implemented to map a corn silk extract in two complementary ways. The first one involved working with UHPLC/HRMS data and Kendrick and van Krevelen plots to highlight a homologous series of compounds, such as lipids from 17 to 23 carbons, monoglycosylated flavonoids from 21 to 24 carbons, diglycosylated flavonoids of 26 to 28 carbons and organic acids of 14 to 19 carbons. The second way was to analyze the sample in UHPLC/HRMS² and to plot mass spectral similarity networks with the GNPS platform and Cytoscape software to refine identification. By combining the information obtained, we were able to propose an identification for 104 detected molecules, including 7 nitrogenous, 28 lipidic and 67 phenolic compounds, leading to the first detailed phytochemical analysis of corn silk extract.

Energy-Resolved Mass Spectrometry as a Tool for Identification of Lignin Depolymerization Products

Lignin is the largest source of bio-based aromatic compounds in nature, and its valorization is essential to the sustainability of lignocellulosic biorefining. Characterizing lignin-derived compounds remains challenging due to the heterogeneity of this biopolymer. Tandem mass spectrometry is a promising tool for lignin structural analytics, as fragmentation patterns of model compounds can be extrapolated to identify characteristic moieties in complex samples. This work extended previous resonance excitation-type collision-induced dissociation (CID) methods that identified lignin oligomers containing β-O-4, β-5, and β-β bonds, to also identify characteristics of 5-5, β-1, and 4-O-5 dimers, enabled by quadrupole time-of-flight (QTOF) CID with energy-resolved mass spectrometry (ERMS). Overall, QTOF-ERMS offers in-depth structural information and could ultimately contribute to tools for high-throughput lignin dimer identification.

Glycoside-specific metabolomics combined with precursor isotopic labeling for characterizing plant glycosyltransferases

Glycosylation by uridine diphosphate-dependent glycosyltransferases (UGTs) in plants contributes to the complexity and diversity of secondary metabolites. UGTs are generally promiscuous in their use of acceptors, making it challenging to reveal the function of UGTs in vivo. Here, we described an approach that combined glycoside-specific metabolomics and precursor isotopic labeling analysis to characterize UGTs in Arabidopsis. We revisited the UGT72E cluster, which has been reported to catalyze the glycosylation of monolignols. Glycoside-specific metabolomics analysis reduced the number of differentially accumulated metabolites in the ugt72e1e2e3 mutant by at least 90% compared with that from traditional untargeted metabolomics analysis. In addition to the two previously reported monolignol glycosides, a total of 62 glycosides showed reduced accumulation in the ugt72e1e2e3 mutant, 22 of which were phenylalanine-derived glycosides, including 5-OH coniferyl alcohol-derived and lignan-derived glycosides, as confirmed by isotopic tracing of [¹³C₆]-phenylalanine precursor. Our method revealed that UGT72Es could use coumarins as substrates, and genetic evidence showed that UGT72Es endowed plants with enhanced tolerance to low iron availability under alkaline conditions. Using the newly developed method, the function of UGT78D2 was also evaluated. These case studies suggest that this method can substantially contribute to the characterization of UGTs and efficiently investigate glycosylation processes, the complexity of which have been highly underestimated.

Targeted and Untargeted Metabolomic Analyses Reveal Organ Specificity of Specialized Metabolites in the Model Grass Brachypodium distachyon

Brachypodium distachyon, because of its fully sequenced genome, is frequently used as a model grass species. However, its metabolome, which constitutes an indispensable element of complex biological systems, remains poorly characterized. In this study, we conducted comprehensive, liquid chromatography-mass spectrometry (LC-MS)-based metabolomic examination of roots, leaves and spikes of Brachypodium Bd21 and Bd3-1 lines. Our pathway enrichment analysis emphasised the accumulation of specialized metabolites representing the flavonoid biosynthetic pathway in parallel with processes related to nucleotide, sugar and amino acid metabolism. Similarities in metabolite profiles between both lines were relatively high in roots and leaves while spikes showed higher metabolic variance within both accessions. In roots, differences between Bd21 and Bd3-1 lines were manifested primarily in diterpenoid metabolism, while differences within spikes and leaves concerned nucleotide metabolism and nitrogen management. Additionally, sulphate-containing metabolites differentiated Bd21 and Bd3-1 lines in spikes. Structural analysis based on MS fragmentation spectra enabled identification of 93 specialized metabolites. Among them phenylpropanoids and flavonoids derivatives were mainly determined. As compared with closely related barley and wheat species, metabolic profile of Brachypodium is characterized with presence of threonate derivatives of hydroxycinnamic acids.

Flexible and digestible wood caused by viral-induced alteration of cell wall composition

Woody plant material represents a vast renewable resource that has the potential to produce biofuels and other bio-based products with favorable net CO₂ emissions.^1,2 Its potential has been demonstrated in a recent study that generated novel structural materials from flexible moldable wood.³ Apple rubbery wood (ARW) disease is the result of a viral infection that causes woody stems to exhibit increased flexibility.⁴ Although ARW disease is associated with the presence of an RNA virus⁵ known as apple rubbery wood virus (ARWV), how the unique symptoms develop is unknown. We demonstrate that the symptoms of ARWV infections arise from reduced lignification within the secondary cell wall of xylem fibers and result in increased wood digestibility. In contrast, the mid-lamellae region and xylem ray cells are largely unaffected by the infection. Gene expression and proteomic data from symptomatic xylem clearly show the downregulation of phenylalanine ammonia lyase (PAL), the enzyme catalyzing the first committed step in the phenylpropanoid pathway leading to lignin biosynthesis. A large increase in soluble phenolics in symptomatic xylem, including the lignin precursor phenylalanine, is also consistent with PAL downregulation. ARWV infection results in the accumulation of many host-derived virus-activated small interfering RNAs (vasiRNAs). PAL-derived vasiRNAs are among the most abundant vasiRNAs in symptomatic xylem and are likely the cause of reduced PAL activity. Apparently, the mechanism used by the virus to alter lignin exhibits similarities to the RNAi strategy used to alter lignin in genetically modified trees to generate comparable improvements in wood properties.6, 7, 8

Video abstract

•

Apple rubbery wood (ARW) symptoms are caused by decreased lignin in woody tissue

•

RNA-seq, proteomics, and metabolomics suggest phenylalanine levels decrease

•

Virus-activated small interfering RNAs (vasiRNAs) are generated in response to ARWV infection

•

VasiRNAs cause siRNA-based downregulation of phenylalanine ammonia

Hydroxycinnamic acid-modified xylan side chains and their cross-linking products in rice cell walls are reduced in the Xylosyl arabinosyl substitution of xylan 1 mutant

The intricate architecture of cell walls and the complex cross-linking of their components hinders some industrial and agricultural applications of plant biomass. Xylan is a key structural element of grass cell walls, closely interacting with other cell wall components such as cellulose and lignin. The main branching points of grass xylan, 3-linked l-arabinosyl substitutions, can be modified by ferulic acid (a hydroxycinnamic acid), which cross-links xylan to other xylan chains and lignin. XAX1 (Xylosyl arabinosyl substitution of xylan 1), a rice (Oryza sativa) member of the glycosyltransferase family GT61, has been described to add xylosyl residues to arabinosyl substitutions modified by ferulic acid. In this study, we characterize hydroxycinnamic acid-decorated arabinosyl substitutions present on rice xylan and their cross-linking, in order to decipher the role of XAX1 in xylan synthesis. Our results show a general reduction of hydroxycinnamic acid-modified 3-linked arabinosyl substitutions in xax1 mutant rice regardless of their modification with a xylosyl residue. Moreover, structures resembling the direct cross-link between xylan and lignin (ferulated arabinosyl substitutions bound to lignin monomers and dimers), together with diferulates known to cross-link xylan, are strongly reduced in xax1. Interestingly, apart from feruloyl and p-coumaroyl modifications on arabinose, putative caffeoyl and oxalyl modifications were characterized, which were also reduced in xax1. Our results suggest an alternative function of XAX1 in the transfer of hydroxycinnamic acid-modified arabinosyl substitutions to xylan, rather than xylosyl transfer to arabinosyl substitutions. Ultimately, XAX1 plays a fundamental role in cross-linking, providing a potential target for the improvement of use of grass biomass.

Rice metabolic regulatory network spanning the entire life cycle

As one of the most important crops in the world, rice (Oryza sativa) is a model plant for metabolome research. Although many studies have focused on the analysis of specific tissues, the changes in metabolite abundance across the entire life cycle have not yet been determined. In this study, combining both targeted and nontargeted metabolite profiling methods, a total of 825 annotated metabolites were quantified in rice samples from different tissues covering the entire life cycle. The contents of metabolites in different tissues of rice were significantly different, with various metabolites accumulating in the plumule and radicle during seed germination. Combining these data with transcriptome data obtained from the same time period, we constructed the Rice Metabolic Regulation Network. The metabolites and co-expressed genes were further divided into 12 clusters according to their accumulation patterns, with members within each cluster displaying a uniform and clear pattern of abundance across development. Using this dataset, we established a comprehensive metabolic profile of the rice life cycle and used two independent strategies to identify novel transcription factors-namely the use of known regulatory genes as bait to screen for new networks underlying lignin metabolism and the unbiased identification of new glycerophospholipid metabolism regulators on the basis of tissue specificity. This study thus demonstrates how guilt-by-association analysis of metabolome and transcriptome data spanning the entire life cycle in cereal crops provides novel resources and tools to aid in understanding the mechanisms underlying important agronomic traits.

Origin and Function of Structural Diversity in the Plant Specialized Metabolome

The plant specialized metabolome consists of a multitude of structurally and functionally diverse metabolites, variable from species to species. The specialized metabolites play roles in the response to environmental changes and abiotic or biotic stresses, as well as in plant growth and development. At its basis, the specialized metabolism is built of four major pathways, each starting from a few distinct primary metabolism precursors, and leading to distinct basic carbon skeleton core structures: polyketides and fatty acid derivatives, terpenoids, alkaloids, and phenolics. Structural diversity in specialized metabolism, however, expands exponentially with each subsequent modification. We review here the major sources of structural variety and question if a specific role can be attributed to each distinct structure. We focus on the influences that various core structures and modifications have on flavonoid antioxidant activity and on the diversity generated by oxidative coupling reactions. We suggest that many oxidative coupling products, triggered by initial radical scavenging, may not have a function in se, but could potentially be enzymatically recycled to effective antioxidants. We further discuss the wide structural variety created by multiple decorations (glycosylations, acylations, prenylations), the formation of high-molecular weight conjugates and polyesters, and the plasticity of the specialized metabolism. We draw attention to the need for untargeted methods to identify the complex, multiply decorated and conjugated compounds, in order to study the functioning of the plant specialized metabolome.

CRISPR-Cas9 editing of CAFFEOYL SHIKIMATE ESTERASE 1 and 2 shows their importance and partial redundancy in lignification in Populus tremula × P. alba

Lignins are cell wall-located aromatic polymers that provide strength and hydrophobicity to woody tissues. Lignin monomers are synthesized via the phenylpropanoid pathway, wherein CAFFEOYL SHIKIMATE ESTERASE (CSE) converts caffeoyl shikimate into caffeic acid. Here, we explored the role of the two CSE homologs in poplar (Populus tremula × P. alba). Reporter lines showed that the expression conferred by both CSE1 and CSE2 promoters is similar. CRISPR-Cas9-generated cse1 and cse2 single mutants had a wild-type lignin level. Nevertheless, CSE1 and CSE2 are not completely redundant, as both single mutants accumulated caffeoyl shikimate. In contrast, the cse1 cse2 double mutants had a 35% reduction in lignin and associated growth penalty. The reduced-lignin content translated into a fourfold increase in cellulose-to-glucose conversion upon limited saccharification. Phenolic profiling of the double mutants revealed large metabolic shifts, including an accumulation of p-coumaroyl, 5-hydroxyferuloyl, feruloyl and sinapoyl shikimate, in addition to caffeoyl shikimate. This indicates that the CSEs have a broad substrate specificity, which was confirmed by in vitro enzyme kinetics. Taken together, our results suggest an alternative path within the phenylpropanoid pathway at the level of the hydroxycinnamoyl-shikimates, and show that CSE is a promising target to improve plants for the biorefinery.