Orthology-based analysis helps map evolutionary diversification and predict substrate class use of BAHD acyltransferases

Applications of ancestral sequence reconstruction for understanding the evolution of plant specialized metabolism

Studies of enzymes in modern-day plants have documented the diversity of metabolic activities retained by species today but only provide limited insight into how those properties evolved. Ancestral sequence reconstruction (ASR) is an approach that provides statistical estimates of ancient plant enzyme sequences which can then be resurrected to test hypotheses about the evolution of catalytic activities and pathway assembly. Here, I review the insights that have been obtained using ASR to study plant metabolism and highlight important methodological aspects. Overall, studies of resurrected plant enzymes show that (i) exaptation is widespread such that even low or undetectable levels of ancestral activity with a substrate can later become the apparent primary activity of descendant enzymes, (ii) intramolecular epistasis may or may not limit evolutionary paths towards catalytic or substrate preference switches, and (iii) ancient pathway flux often differs from modern-day metabolic networks. These and other insights gained from ASR would not have been possible using only modern-day sequences. Future ASR studies characterizing entire ancestral metabolic networks as well as those that link ancient structures with enzymatic properties should continue to provide novel insights into how the chemical diversity of plants evolved. This article is part of the theme issue 'The evolution of plant metabolism'.

Phylogenomic and synteny analysis of BAHD and SCP/SCPL gene families reveal their evolutionary histories in plant specialized metabolism

Plant chemical diversity is largely owing to a number of enzymes which catalyse reactions involved in the assembly, and in the subsequent chemical modifications, of the core structures of major classes of plant specialized metabolites. One such reaction is acylation. With this in mind, to study the deep evolutionary history of BAHD and the serine-carboxypeptidase-like (SCPL) acyltransferase genes, we assembled phylogenomic synteny networks based on a large-scale inference analysis of orthologues across whole-genome sequences of 126 species spanning Stramenopiles and Archaeplastida, including Arabidopsis thaliana, tomato (Solanum lycopersicum) and maize (Zea mays). As such, this study combined the study of genomic location with changes in gene sequences. Our analyses revealed that serine-carboxypeptidase (SCP)/serine-carboxypeptidase-like (SCPL) genes had a deeper evolutionary origin than BAHD genes, which expanded massively on the transition to land and with the development of the vascular system. The two gene families additionally display quite distinct patterns of copy number variation across phylogenies as well as differences in cross-phylogenetic syntenic network components. In unlocking the above observations, our analyses demonstrate the possibilities afforded by modern phylogenomic (syntenic) networks, but also highlight their current limitations, as demonstrated by the inability of phylogenetic methods to separate authentic SCPL acyltransferases from standard SCP peptide hydrolases.This article is part of the theme issue 'The evolution of plant metabolism'.

Unlocking saponin biosynthesis in soapwort

Soapwort (Saponaria officinalis) is a flowering plant from the Caryophyllaceae family with a long history of human use as a traditional source of soap. Its detergent properties are because of the production of polar compounds (saponins), of which the oleanane-based triterpenoid saponins, saponariosides A and B, are the major components. Soapwort saponins have anticancer properties and are also of interest as endosomal escape enhancers for targeted tumor therapies. Intriguingly, these saponins share common structural features with the vaccine adjuvant QS-21 and, thus, represent a potential alternative supply of saponin adjuvant precursors. Here, we sequence the S. officinalis genome and, through genome mining and combinatorial expression, identify 14 enzymes that complete the biosynthetic pathway to saponarioside B. These enzymes include a noncanonical cytosolic GH1 (glycoside hydrolase family 1) transglycosidase required for the addition of D-quinovose. Our results open avenues for accessing and engineering natural and new-to-nature pharmaceuticals, drug delivery agents and potential immunostimulants.

The impact of the GLOSSY2 and GLOSSY2-LIKE BAHD-proteins in affecting the product profile of the maize fatty acid elongase

The maize glossy2 and glossy2-like genes are homologs, which encode proteins that belong to the BAHD family of acyltransferases. In planta genetic studies have demonstrated that these genes may be involved in the elongation of very long chain fatty acids (VLCFAs) that are precursors of the cuticular wax fraction of the plant cuticle. VLCFAs are synthesized by a fatty acyl-CoA elongase complex (FAE) that consists of four component enzymes. Previously, we functionally identified the maize FAE component enzymes by their ability to complement haploid Saccharomyces cerevisiae strains that carry lethal deletion alleles for each FAE component enzyme. In this study we used these complemented haploid strains and wild-type diploid strains to evaluate whether the co-expression of either GLOSSY2 or GLOSSY2-LIKE with individual maize FAE component enzymes affects the VLCFA product-profile of the FAE system. Wild-type diploid strains produced VLCFAs of up to 28-carbon chain length. Co-expression of GLOSSY2 or GLOSSY2-LIKE with a combination of maize 3-ketoacyl-CoA synthases stimulated the synthesis of longer VLCFAs, up to 30-carbon chain lengths. However, such results could not be recapitulated when these co-expression experiments were conducted in the yeast haploid mutant strains that lacked individual components of the endogenous FAE system. Specifically, lethal yeast mutant strains that are genetically complemented by the expression of maize FAE-component enzymes produce VLCFAs that range between 20- and 26-carbon chain lengths. However, expressing either GLOSSY2 or GLOSSY2-LIKE in these complemented strains does not enable the synthesis of longer chain VLCFAs. These results indicate that the apparent stimulatory role of GLOSSY2 or GLOSSY2-LIKE to enable the synthesis of longer chain VLCFAs in diploid yeast cells may be associated with mixing plant enzyme components with the endogenous FAE complex.

Functional Characterization and Catalytic Activity Improvement of Borneol Acetyltransferase from Wurfbainia longiligularis

In plant secondary metabolite biosynthesis, acylation is a diverse physiological process, with BAHD acyltransferases playing an essential role. Borneol acetyltransferase (BAT) is an alcohol acetyltransferase, which catalyzes borneol and acetyl-CoA to synthesize bornyl acetate (BA). However, the enzymes involved in the biosynthesis of BA have so far only been characterized in Wurfbainia villosa, the studies on the WvBATs have only been conducted in vitro, and the catalytic activity was relatively low. In this research, three genes (WlBAT1, WlBAT2, and WlBAT3) have been identified to encode BATs that are capable of acetylating borneol to synthesize BA in vitro. We also determined that WlBAT1 has the highest catalytic efficiency for borneol-type substrates, including (+)-borneol, (-)-borneol, and isoborneol. Furthermore, we found that BATs could catalyze a wide range of substrate types in vitro, but in vivo, they exclusively catalyzed borneol-type substrates. Through molecular simulations and site-directed mutagenesis, it was revealed that residues D32, N36, H168, N297, N355, and H384 are crucial for the catalytic activity of WlBAT1, while the R382I-D385R double mutant of WlBAT1 exhibited an increasing acylation efficiency for borneol-type substrates in vitro and in vivo. These findings offer key genetic elements for the metabolic engineering of plants and synthetic biology to produce BA.

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.

Seeing double: two different homospermidine oxidases are involved in pyrrolizidine alkaloid biosynthesis in different organs of comfrey (Symphytum officinale)

Pyrrolizidine alkaloids (PAs) are toxic specialized metabolites produced in several plant species and frequently contaminate herbal teas or livestock feed. In comfrey (Symphytum officinale, Boraginaceae), they are produced in two different organs of the plant, the root and young leaves. In this study, we demonstrate that homospermidine oxidase (HSO), a copper-containing amine oxidase (CuAO) responsible for catalyzing the formation of the distinctive pyrrolizidine ring in PAs, is encoded by two individual genes. Specifically, SoCuAO1 is expressed in young leaves, while SoCuAO5 is expressed in roots. CRISPR/Cas9-mediated knockout of socuao5 resulted in hairy roots (HRs) unable to produce PAs, supporting its function as HSO in roots. Plants regenerated from socuao5 knockout HRs remained completely PA-free until the plants began to develop inflorescences, indicating the presence of another HSO that is expressed only during flower development. Stable expression of SoCuAO1 in socuao5 knockout HRs rescued the ability to produce PAs. In vitro assays of both enzymes transiently expressed in Nicotiana benthamiana confirmed their HSO activity and revealed the ability of HSO to control the stereospecific cyclization of the pyrrolizidine backbone. The observation that the first specific step of PA biosynthesis catalyzed by homospermidine synthase requires only one gene copy, while two independent paralogs are recruited for the subsequent homospermidine oxidation in different tissues of the plant, suggests a complex regulation of the pathway. This adds a new level of complexity to PA biosynthesis, a system already characterized by species-specific, tight spatio-temporal regulation, and independent evolutionary origins in multiple plant lineages.

Discovering a mitochondrion-localized BAHD acyltransferase involved in calystegine biosynthesis and engineering the production of 3β-tigloyloxytropane

Solanaceous plants produce tropane alkaloids (TAs) via esterification of 3α- and 3β-tropanol. Although littorine synthase is revealed to be responsible for 3α-tropanol esterification that leads to hyoscyamine biosynthesis, the genes associated with 3β-tropanol esterification are unknown. Here, we report that a BAHD acyltransferase from Atropa belladonna, 3β-tigloyloxytropane synthase (TS), catalyzes 3β-tropanol and tigloyl-CoA to form 3β-tigloyloxytropane, the key intermediate in calystegine biosynthesis and a potential drug for treating neurodegenerative disease. Unlike other cytosolic-localized BAHD acyltransferases, TS is localized to mitochondria. The catalytic mechanism of TS is revealed through molecular docking and site-directed mutagenesis. Subsequently, 3β-tigloyloxytropane is synthesized in tobacco. A bacterial CoA ligase (PcICS) is found to synthesize tigloyl-CoA, an acyl donor for 3β-tigloyloxytropane biosynthesis. By expressing TS mutant and PcICS, engineered Escherichia coli synthesizes 3β-tigloyloxytropane from tiglic acid and 3β-tropanol. This study helps to characterize the enzymology and chemodiversity of TAs and provides an approach for producing 3β-tigloyloxytropane.

Genome-wide identification of tea plant (Camellia sinensis) BAHD acyltransferases reveals their role in response to herbivorous pests

BACKGROUND

BAHD acyltransferases are among the largest metabolic protein domain families in the genomes of terrestrial plants and play important roles in plant growth and development, aroma formation, and biotic and abiotic stress responses. Little is known about the BAHDs in the tea plant, a cash crop rich in secondary metabolites.

RESULTS

In this study, 112 BAHD genes (CsBAHD01-CsBAHD112) were identified from the tea plant genome, with 85% (98/112) unevenly distributed across the 15 chromosomes. The number of BAHD gene family members has significantly expanded from wild tea plants to the assamica type to the sinensis type. Phylogenetic analysis showed that they could be classified into seven subgroups. Promoter cis-acting element analysis revealed that they contain a large number of light, phytohormones, and stress-responsive elements. Many members displayed tissue-specific expression patterns. CsBAHD05 was expressed at more than 500-fold higher levels in purple tea leaves than in green tea leaves. The genes exhibiting the most significant response to MeJA treatment and feeding by herbivorous pests were primarily concentrated in subgroups 5 and 6. The expression of 23 members of these two subgroups at different time points after feeding by tea green leafhoppers and tea geometrids was examined via qPCR, and the results revealed that the expression of CsBAHD93, CsBAHD94 and CsBAHD95 was significantly induced after the tea plants were subjected to feeding by both pricking and chewing pests. Moreover, based on the transcriptome data for tea plants being fed on by these two pests, a transcriptional regulatory network of different transcription factor genes coexpressed with these 23 members was constructed.

CONCLUSIONS

Our study provides new insights into the role of BAHDs in the defense response of tea plants, and will facilitate in-depth studies of the molecular function of BAHDs in resistance to herbivorous pests.

Complete biosynthesis of the potent vaccine adjuvant QS-21

QS-21 is a potent vaccine adjuvant currently sourced by extraction from the Chilean soapbark tree. It is a key component of human vaccines for shingles, malaria, coronavirus disease 2019 and others under development. The structure of QS-21 consists of a glycosylated triterpene scaffold coupled to a complex glycosylated 18-carbon acyl chain that is critical for immunostimulant activity. We previously identified the early pathway steps needed to make the triterpene glycoside scaffold; however, the biosynthetic route to the acyl chain, which is needed for stimulation of T cell proliferation, was unknown. Here, we report the biogenic origin of the acyl chain, characterize the series of enzymes required for its synthesis and addition and reconstitute the entire 20-step pathway in tobacco, thereby demonstrating the production of QS-21 in a heterologous expression system. This advance opens up unprecedented opportunities for bioengineering of vaccine adjuvants, investigating structure-activity relationships and understanding the mechanisms by which these compounds promote the human immune response.

Improved production of the antidiabetic metabolite montbretin A in Nicotiana benthamiana: discovery, characterization, and use of Crocosmia shikimate shunt genes

The plant-specialized metabolite montbretin A (MbA) is being developed as a new treatment option for type-2 diabetes, which is among the ten leading causes of premature death and disability worldwide. MbA is a complex acylated flavonoid glycoside produced in small amounts in below-ground organs of the perennial plant Montbretia (Crocosmia × crocosmiiflora). The lack of a scalable production system limits the development and potential application of MbA as a pharmaceutical or nutraceutical. Previous efforts to reconstruct montbretin biosynthesis in Nicotiana benthamiana (Nb) resulted in low yields of MbA and higher levels of montbretin B (MbB) and montbretin C (MbC). MbA, MbB, and MbC are nearly identical metabolites differing only in their acyl moieties, derived from caffeoyl-CoA, coumaroyl-CoA, and feruloyl-CoA, respectively. In contrast to MbA, MbB and MbC are not pharmaceutically active. To utilize the montbretia caffeoyl-CoA biosynthesis for improved MbA engineering in Nb, we cloned and characterized enzymes of the shikimate shunt of the general phenylpropanoid pathway, specifically hydroxycinnamoyl-CoA: shikimate hydroxycinnamoyl transferase (CcHCT), p-coumaroylshikimate 3'-hydroxylase (CcC3'H), and caffeoylshikimate esterase (CcCSE). Gene expression patterns suggest that CcCSE enables the predominant formation of MbA, relative to MbB and MbC, in montbretia. This observation is supported by results from in vitro characterization of CcCSE and reconstruction of the shikimate shunt in yeast. Using CcHCT together with montbretin biosynthetic genes in multigene constructs resulted in a 30-fold increase of MbA in Nb. This work advances our understanding of the phenylpropanoid pathway and features a critical step towards improved MbA production in bioengineered Nb.

Sequence-Structure Analysis Unlocking the Potential Functional Application of the Local 3D Motifs of Plant-Derived Diterpene Synthases

Plant-derived diterpene synthases (PdiTPSs) play a critical role in the formation of structurally and functionally diverse diterpenoids. However, the specificity or functional-related features of PdiTPSs are not well understood. For a more profound insight, we collected, constructed, and curated 199 functionally characterized PdiTPSs and their corresponding 3D structures. The complex correlations among their sequences, domains, structures, and corresponding products were comprehensively analyzed. Ultimately, our focus narrowed to the geometric arrangement of local structures. We found that local structural alignment can rapidly localize product-specific residues that have been validated by mutagenesis experiments. Based on the 3D motifs derived from the residues around the substrate, we successfully searched diterpene synthases (diTPSs) from the predicted terpene synthases and newly characterized PdiTPSs, suggesting that the identified 3D motifs can serve as distinctive signatures in diTPSs (I and II class). Local structural analysis revealed the PdiTPSs with more conserved amino acid residues show features unique to class I and class II, whereas those with fewer conserved amino acid residues typically exhibit product diversity and specificity. These results provide an attractive method for discovering novel or functionally equivalent enzymes and probing the product specificity in cases where enzyme characterization is limited.

Family characteristics, phylogenetic reconstruction, and potential applications of the plant BAHD acyltransferase family

The BAHD acyltransferase family is a class of proteins in plants that can acylate a variety of primary and specialized secondary metabolites. The typically acylated products have greatly improved stability, lipid solubility, and bioavailability and thus show significant differences in their physicochemical properties and pharmacological activities. Here, we review the protein structure, catalytic mechanism, and phylogenetic reconstruction of plant BAHD acyltransferases to describe their family characteristics, acylation reactions, and the processes of potential functional differentiation. Moreover, the potential applications of the BAHD family in human activities are discussed from the perspectives of improving the quality of economic plants, enhancing the efficacy of medicinal plants, improving plant biomass for use in biofuel, and promoting stress resistance of land plants. This review provides a reference for the research and production of plant BAHD acyltransferases.

A BAHD-type acyltransferase concludes the biosynthetic pathway of non-bitter glycoalkaloids in ripe tomato fruit

Tomato is the highest value fruit and vegetable crop worldwide, yet produces α-tomatine, a renowned toxic and bitter-tasting anti-nutritional steroidal glycoalkaloid (SGA) involved in plant defense. A suite of modifications during tomato fruit maturation and ripening converts α-tomatine to the non-bitter and less toxic Esculeoside A. This important metabolic shift prevents bitterness and toxicity in ripe tomato fruit. While the enzymes catalyzing glycosylation and hydroxylation reactions in the Esculeoside A pathway have been resolved, the proposed acetylating step remains, to date, elusive. Here, we discovered that GAME36 (GLYCOALKALOID METABOLISM36), a BAHD-type acyltransferase catalyzes SGA-acetylation in cultivated and wild tomatoes. This finding completes the elucidation of the core Esculeoside A biosynthetic pathway in ripe tomato, allowing reconstitution of Esculeoside A production in heterologous microbial and plant hosts. The involvement of GAME36 in bitter SGA detoxification pathway points to a key role in the evolution of sweet-tasting tomato as well as in the domestication and breeding of modern cultivated tomato fruit.

Dangerous sugars: Structural diversity and functional significance of acylsugar-like defense compounds in flowering plants

Acylsugars constitute a diverse class of secondary metabolites found in many flowering plant families. Comprising sugar cores and acyl groups connected by ester and/or ether linkages, acylsugar structures vary considerably at all taxonomic levels - from populations of the same species to across species of the same family and across flowering plants, with some species producing hundreds of acylsugars in a single organ. Acylsugars have been most well-studied in the Solanaceae family, but structurally analogous compounds have also been reported in the Convolvulaceae, Martyniaceae, Geraniaceae, Rubiaceae, Rosaceae and Caryophyllaceae families. Focusing on Solanaceae and Convolvulaceae acylsugars, this review highlights their structural diversity, the potential biosynthetic mechanisms that produce this diversity, and its functional significance. Finally, we also discuss the possibility that some of this diversity is merely "noise", arising out of enzyme promiscuity and/or non-adaptive evolutionary mechanisms.

Phylogenomic analyses across land plants reveals motifs and coexpression patterns useful for functional prediction in the BAHD acyltransferase family

The BAHD acyltransferase family is one of the largest enzyme families in flowering plants, containing dozens to hundreds of genes in individual genomes. Highly prevalent in angiosperm genomes, members of this family contribute to several pathways in primary and specialized metabolism. In this study, we performed a phylogenomic analysis of the family using 52 genomes across the plant kingdom to gain deeper insights into its functional evolution and enable function prediction. We found that BAHD expansion in land plants was associated with significant changes in various gene features. Using pre-defined BAHD clades, we identified clade expansions in different plant groups. In some groups, these expansions coincided with the prominence of metabolite classes such as anthocyanins (flowering plants) and hydroxycinnamic acid amides (monocots). Clade-wise motif-enrichment analysis revealed that some clades have novel motifs fixed on either the acceptor or the donor side, potentially reflecting historical routes of functional evolution. Co-expression analysis in rice and Arabidopsis further identified BAHDs with similar expression patterns, however, most co-expressed BAHDs belonged to different clades. Comparing BAHD paralogs, we found that gene expression diverges rapidly after duplication, suggesting that sub/neo-functionalization of duplicate genes occurs quickly via expression diversification. Analyzing co-expression patterns in Arabidopsis in conjunction with orthology-based substrate class predictions and metabolic pathway models led to the recovery of metabolic processes of most of the already-characterized BAHDs as well as definition of novel functional predictions for some uncharacterized BAHDs. Overall, this study provides new insights into the evolution of BAHD acyltransferases and sets up a foundation for their functional characterization.

Near-real time determination of BAHD acyl-coenzyme A transferase reaction rates and kinetic parameters using Ellman's reagent

BAHD acyl-coenzyme A (CoA) acyltransferases play key roles in a large number of biosynthetic reactions involved in plant specialized metabolism. One approach to measure reaction rates for these enzymes is to quantify the amide or ester reaction products following chromatographic separation of reaction components, an approach that can be labor intensive and time consuming, and complicated by a lack of pure standards. We previously developed and validated an alternative approach using 5,5'-dithio-bis-(2-nitrobenzoic acid) (DTNB, Ellman's reagent) to spectrophotometrically monitor reaction progress by the release of free CoA in the reaction. This approach allows near-real time measurement of reaction rates, permitting reaction conditions (buffer, reactant, and enzyme concentrations, etc.) to be changed "on the fly." The ease and rapidity of data collection allows a high density of data points to be collected for determination of kinetic parameters. Here we provide a detailed procedure for using DTNB to measure BAHD acyl-CoA acyltransferase reaction rates, and as an example, use it to determine kinetic parameters for red clover hydroxycinnamoyl-CoA:shikimate hydroxycinnamoyltransferase, a BAHD acyl-CoA hydroxycinnamoyltransferase not previously characterized with respect to kinetic parameters. This approach may be more generally applicable to transferases using CoA donors.