Evolution and Function of the Populus SABATH Family Reveal That a Single Amino Acid Change Results in a Substrate Switch

Molecular basis of one-step methyl anthranilate biosynthesis in grapes, sweet orange, and maize

Plants synthesize an array of volatile compounds, many of which serve ecological roles in attracting pollinators, deterring herbivores, and communicating with their surroundings. Methyl anthranilate (MeAA) is an anti-herbivory defensive volatile responsible for grape aroma that is emitted by several agriculturally relevant plants, including citrus, grapes, and maize. Unlike maize, which uses a one-step anthranilate methyltransferase (AAMT), grapes have been thought to use a two-step pathway for MeAA biosynthesis. By mining available transcriptomics data, we identified two AAMTs in Vitis vinifera (wine grape), as well as one ortholog in "Concord" grape. Many angiosperms methylate the plant hormone salicylic acid (SA) to produce methyl salicylate, which acts as a plant-to-plant communication molecule. Because the Citrus sinensis (sweet orange) SA methyltransferase can methylate both anthranilate (AA) and SA, we used this enzyme to examine the molecular basis of AA activity by introducing rational mutations, which identified several active site residues that increase activity with AA. Reversing this approach, we introduced mutations that imparted activity with SA in the maize AAMT, which uncovered different active site residues from those in the citrus enzyme. Sequence and phylogenetic analysis revealed that one of the Vitis AAMTs shares an ancestor with jasmonic acid methyltransferases, similar to the AAMT from strawberry (Frageria sp.). Collectively, these data demonstrate the molecular mechanisms underpinning AA activity across methyltransferases and identify one-step enzymes by which grapes synthesize MeAA.

The Mutational Road not Taken: Using Ancestral Sequence Resurrection to Evaluate the Evolution of Plant Enzyme Substrate Preferences

We investigated the flowering plant salicylic acid methyl transferase (SAMT) enzyme lineage to understand the evolution of substrate preference change. Previous studies indicated that a single amino acid replacement to the SAMT active site (H150M) was sufficient to change ancestral enzyme substrate preference from benzoic acid to the structurally similar substrate, salicylic acid (SA). Yet, subsequent studies have shown that the H150M function-changing replacement did not likely occur during the historical episode of enzymatic divergence studied. Therefore, we reinvestigated the origin of SA methylation preference here and additionally assessed the extent to which epistasis may act to limit mutational paths. We found that the SAMT lineage of enzymes acquired preference to methylate SA from an ancestor that preferred to methylate benzoic acid as previously reported. In contrast, we found that a different amino acid replacement, Y267Q, was sufficient to change substrate preference with others providing small positive-magnitude epistatic improvements. We show that the kinetic basis for the ancestral enzymatic change in substate preference by Y267Q appears to be due to both a reduced specificity constant, kcat/KM, for benzoic acid and an improvement in KM for SA. Therefore, this lineage of enzymes appears to have had multiple mutational paths available to achieve the same evolutionary divergence. While the reasons remain unclear for why one path was taken, and the other was not, the mutational distance between ancestral and descendant codons may be a factor.

Identification and Characterization of Jasmonic Acid Methyltransferase Involved in the Formation of Floral Methyl Jasmonate in Hedychium coronarium

Hedychium coronarium is a popular ornamental flower in tropical and subtropical areas due to its elegant appearance and inviting fragrance. Methyl jasmonate (MeJA) is one of the volatile compounds in the blooming flowers of H. coronarium. However, the molecular mechanism underlying floral MeJA formation is still unclear in H. coronarium. In this study, a total of 12 SABATH family genes were identified in the genome of H. coronarium, and their encoded proteins range from 366 to 387 amino acids. Phylogenetic analysis revealed seven clades in the SABATH family and a JMT ortholog clade, including two HcSABATH members. Combined with expression profiling of HcSABATH members, HcJMT1 was identified as the top candidate gene for floral MeJA biosynthesis. In vitro enzyme assays showed that HcJMT1 can catalyze the production of MeJA from jasmonic acid. Gene expression analysis indicated that HcJMT1 exhibited the highest expression in the labella and lateral petals, the major sites of MeJA emission. During flower development, the two MeJA isomers, major isomers in the products of the HcJMT1 protein, were released after anthesis, in which stage HcJMT1 displayed high expression. Our results indicated that HcJMT1 is involved in the formation of floral MeJA in H. coronarium.

Genomic and structural basis for evolution of tropane alkaloid biosynthesis

The tropane alkaloids (TAs) cocaine and hyoscyamine have been used medicinally for thousands of years. To understand the evolutionary origins and trajectories of serial biosynthetic enzymes of TAs and especially the characteristic tropane skeletons, we generated the chromosome-level genome assemblies of cocaine-producing Erythroxylum novogranatense (Erythroxylaceae, rosids clade) and hyoscyamine-producing Anisodus acutangulus (Solanaceae, asterids clade). Comparative genomic and phylogenetic analysis suggested that the lack of spermidine synthase/N-methyltransferase (EnSPMT1) in ancestral asterids species contributed to the divergence of polyamine (spermidine or putrescine) methylation in cocaine and hyoscyamine biosynthesis. Molecular docking analysis and key site mutation experiments suggested that ecgonone synthases CYP81AN15 and CYP82M3 adopt different active-site architectures to biosynthesize the same product ecgonone from the same substrate in Erythroxylaceae and Solanaceae. Further synteny analysis showed different evolutionary origins and trajectories of CYP81AN15 and CYP82M3, particularly the emergence of CYP81AN15 through the neofunctionalization of ancient tandem duplication genes. The combination of structural biology and comparative genomic analysis revealed that ecgonone methyltransferase, which is responsible for the biosynthesis of characteristic 2-substituted carboxymethyl group in cocaine, evolved from the tandem copies of salicylic acid methyltransferase by the mutations of critical E216 and S153 residues. Overall, we provided strong evidence for the independent origins of serial TA biosynthetic enzymes on the genomic and structural level, underlying the chemotypic convergence of TAs in phylogenetically distant species.

Discovery of specific catalytic activity toward IAA/FA by LaSABATHs based on genome-wide phylogenetic and enzymatic analysis of SABATH gene family from Larix kaempferi

The SABATH methyltransferases catalyze methylation of small-molecule metabolites, which participate in plant growth, development and defense response. Given lack of genome-wide studies on gymnosperms SABATH family, the formation and functional differentiation mechanism of the Larix kaempferi SABATH gene family was systematically and exhaustively explored by analyzing gene sequence characteristics, phylogenetic relationship, expression pattern, and enzyme activities. Phylogenetic analysis showed that 247 SABATH genes from 14 land plants were divided into 4 clades, and lineage-specific gene duplication events were important factors that contributed to the evolution of the SABATH gene family in gymnosperms and angiosperms. Substrate specificity analysis of 18 Larix SABATH proteins showed that LaSABATHs could catalyze O-methylation of indole-3-acetic acid (IAA) and farnesic acid (FA), N-methylation of theobromine, and S-methylation of thiobenzoic acid. Furthermore, only LaSABATH2 and LaSABATH29 could catalyze O-methylation of FA, and only LaSABATH30 could catalyze O-methylation of IAA. Homology modeling and molecular docking studies showed the hydrogen bond formed between the His188 of LaSABATH30 and IAA and the noticeable hydrophobic IAA-binding pocket may be helpful for IAA methylation. In this study, identification of proteins with significant specific catalytic activity toward FA and IAA provided high-quality candidate genes for forest genetics and breeding.

A collaborative classroom investigation of the evolution of SABATH methyltransferase substrate preference shifts over 120 million years of flowering plant history

Next-generation sequencing has resulted in an explosion of available data, much of which remains unstudied in terms of biochemical function; yet, experimental characterization of these sequences has the potential to provide unprecedented insight into the evolution of enzyme activity. One way to make inroads into the experimental study of the voluminous data available is to engage students by integrating teaching and research in a college classroom such that eventually hundreds or thousands of enzymes may be characterized. In this study, we capitalize on this potential to focus on SABATH methyltransferase enzymes that have been shown to methylate the important plant hormone, salicylic acid (SA), to form methyl salicylate. We analyze data from 76 enzymes of flowering plant species in 23 orders and 41 families to investigate how widely conserved substrate preference is for SA methyltransferase orthologs. We find a high degree of conservation of substrate preference for SA over the structurally similar metabolite, benzoic acid, with recent switches that appear to be associated with gene duplication and at least three cases of functional compensation by paralogous enzymes. The presence of Met in active site position 150 is a useful predictor of SA methylation preference in SABATH methyltransferases but enzymes with other residues in the homologous position show the same substrate preference. Although our dense and systematic sampling of SABATH enzymes across angiosperms has revealed novel insights, this is merely the “tip of the iceberg” since thousands of sequences remain uncharacterized in this enzyme family alone.

Identification and analysis of key genes involved in methyl salicylate biosynthesis in different birch species

Species of the perennial woody plant genus Betula dominate subalpine forests and play a significant role in preserving biological diversity. In addition to their conventional benefits, birches synthesize a wide range of secondary metabolites having pharmacological significance. Methyl salicylate (MeSA) is one of these naturally occurring compounds constitutively produced by different birch species. MeSA is therapeutically important in human medicine for muscle injuries and joint pain. However, MeSA is now mainly produced synthetically due to a lack of information relating to MeSA biosynthesis and regulation. In this study, we performed a comprehensive bioinformatics analysis of two candidate genes mediating MeSA biosynthesis, SALICYLIC ACID METHYLTRANSFERASE (SAMT) and SALICYLIC ACID-BINDING PROTEIN 2 (SABP2), of high (B. lenta, B. alleghaniensis, B. medwediewii, and B. grossa) and low (B. pendula, B. utilis, B. alnoides, and B. nana) MeSA-producing birch species. Phylogenetic analyses of SAMT and SABP2 genes and homologous genes from other plant species confirmed their evolutionary relationships. Multiple sequence alignments of the amino acid revealed the occurrence of important residues for substrate specificity in SAMT and SABP2. The analysis of cis elements in different birches indicated a functional multiplicity of SAMT and SABP2 and provided insights into the regulation of both genes. We successfully developed six prominent single nucleotide substitution markers that were validated with 38 additional birch individuals to differentiate high and low MeSA-producing birch species. Relative tissue-specific expression analysis of SAMT in leaf and bark tissue of two high and two low MeSA-synthesizing birches revealed a high expression in the bark of both high MeSA-synthesizing birches. In contrast, SABP2 expression in tissues revealed indifferent levels of expression between species belonging to the two groups. The comparative expression and bioinformatics analyses provided vital information that could be used to apply plant genetic engineering technology in the mass production of organic MeSA.

Genome-wide identification and expression analysis of SABATH methyltransferases in tea plant (Camellia sinensis): insights into their roles in plant defense responses

SABATH methyltransferases convent plant small-molecule metabolites into volatile methyl esters, which play important roles in many biological processes and defense reactions in plants. In this study, a total of 32 SABATH genes were identified in the Camellia sinensis var. sinensis (CSS) genome, which were renamed CsSABATH1 to CsSABATH32. Genome location annotation suggested that tandem duplication was responsible for the expansion of SABATH genes in tea plant. Multiple sequence alignment and phylogenetic analysis showed that the CsSABATHs could be classified into three groups (I, II and III), which were also supported by gene structures and conserved motifs analysis. Group II contained only two CsSABATH proteins, which were closely related to PtIAMT, AtIAMT and OsIAMT. The group III SABATH genes of tea plant exhibited expansion on the CSS genome compared with Camellia sinensis var. assamica (CSA) genome. Based on RNA-seq data, the CsSABATHs exhibited tissue-specific expression patterns, and the members with high expression in buds and young leaves were also obviously upregulated after MeJA treatment. The expression of many transcription factors was significantly correlated with that of different members of the CsSABATH gene family, suggesting a potential regulatory relationship between them. Quantitative real-time PCR (qPCR) expression analysis showed that CsSABATHs could respond to exogenous JA, SA and MeSA treatments in tea plants. RNA-seq data analysis and qPCR validation suggested that CsSABATH8, 11, 16, 25, 29 and 32 might play a special role in plant defense against insect herbivory. These results provide references for evolutionary studies of the plant SABATH family and the exploration of the potential roles of CsSABATHs in tea plant defense responses.

Salicylic acid in Populus tomentosa is a remote signalling molecule induced by Botryosphaeria dothidea infection

The salicylic acid (SA) plays a critical role during the establishment of systemic acquired resistance (SAR) in uninfected plant tissues after localised exposure to a pathogen. Here, we studied SA in Populus tomentosa infected by the plant pathogen Botryosphaeria dothidea. The accumulation of SA and methyl salicylate (MeSA) occurred in chronological order in P. tomentosa. The SA and MeSA contents were greater at infected than uninfected sites. Additionally, a gene expression analysis indicated that SA might be accumulated by phenylalanine ammonialyase (PAL) and converted to MeSA by SA carboxyl methyltransferase (SAMT), while MeSA might convert to SA by SA-binding protein 2 (SABP2). The expressions of SAMT at infected sites and SABP2 at uninfected sites, respectively, were significantly up-regulated. Thus, SA might be converted to MeSA at infected sites and transported as a signalling molecule to uninfected sites, where it is converted to SA for SAR. Moreover, the expressions of pathogenesis-related genes PR-1, PR-2 and PR-5 in P. tomentosa were up-regulated by the B. dothidea infection. Our study determined that variations in SA and MeSA contents occur at infected and uninfected sites in poplar after pathogen infection and contributed to the remote signals for poplar SAR.