Molecular Characterization and Phylogenetic Analysis of Two Novel Regio-specific Flavonoid Prenyltransferases from Morus alba and Cudrania tricuspidata*

Background: Plant flavonoid prenyltransferases (FPTs) transfer prenyl moiety to flavonoid cores and have previously been identified only in Leguminosae. Results: The newly identified moraceous FPTs, MaIDT, and CtIDT, are distantly related to leguminous FPTs and feature catalytic regioselectivity and promiscuity. Conclusion: MaIDT and CtIDT evolved independently from leguminous FPTs. Significance: These findings are valuable for identifying additional evolutionarily different non-Leguminosae FPTs. Prenylated flavonoids are attractive specialized metabolites with a wide range of biological activities and are distributed in several plant families. The prenylation catalyzed by prenyltransferases represents a Friedel-Crafts alkylation of the flavonoid skeleton in the biosynthesis of natural prenylated flavonoids and contributes to the structural diversity and biological activities of these compounds. To date, all identified plant flavonoid prenyltransferases (FPTs) have been identified in Leguminosae. In the present study two new FPTs, Morus alba isoliquiritigenin 3′-dimethylallyltransferase (MaIDT) and Cudrania tricuspidata isoliquiritigenin 3′-dimethylallyltransferase (CtIDT), were identified from moraceous plants M. alba and C. tricuspidata, respectively. MaIDT and CtIDT shared low levels of homology with the leguminous FPTs. MaIDT and CtIDT are predicted to be membrane-bound proteins with predicted transit peptides, seven transmembrane regions, and conserved functional domains that are similar to other homogentisate prenyltransferases. Recombinant MaIDT and CtIDT were able to regioselectively introduce dimethylallyl diphosphate into the A ring of three flavonoids with different skeleton types (chalcones, isoflavones, and flavones). Phylogenetic analysis revealed that MaIDT and CtIDT are distantly related to their homologs in Leguminosae, which suggests that FPTs in Moraceae and Leguminosae might have evolved independently. MaIDT and CtIDT represent the first two non-Leguminosae FPTs to be identified in plants and could thus lead to the identification of additional evolutionarily varied FPTs in other non-Leguminosae plants and could elucidate the biosyntheses of prenylated flavonoids in various plants. Furthermore, MaIDT and CtIDT might be used for regiospecific prenylation of flavonoids to produce bioactive compounds for potential therapeutic applications due to their high efficiency and catalytic promiscuity.

Prenylated flavonoids are an attractive class of plant-specialized metabolites that are primarily distributed in Leguminosae (Fabaceae), Moraceae, Umbelliferae, Guttiferae, Euphorbiaceae, Celastraceae, Compositae, Paulowniaceae, and Zingiberaceae (1)(2)(3)(4)(5). These compounds are hybrid molecules that contain structurally divergent flavonoid skeletons and extensively modified prenyl moieties with different chain lengths, and these compounds exhibit a wide range of biological activities, e.g. antioxidative, antitumor, antibacterial, antiviral, and estrogenic activities (1,3). Notably, the substitution of the flavonoid ring system with prenyl group(s) contributes strongly to these biological activities because the prenyl group(s) increases the lipophilicity of the flavonoids and confers to these molecules a strong affinity to biological membranes that results in greater bioavailability (6).
Due to their diverse chemical structures and impressive biological activities, prenylated flavonoids have been studied in many fields including plant physiology, natural product chemistry, and synthetic chemistry. The prenylation reactions catalyzed by flavonoid prenyltransferases (FPTs) 2 have received particular attention and significantly contribute to the structural diversity and further modification of the prenylated flavonoids (7). For example, prenylation is a key step in the biosynthesis of diverse Diels-Alder-type adducts in mulberry trees (8,9). Consequently, intensive biochemical studies of FPTs * The study was supported by National Natural Science Foundation of China have been performed for more than two decades (10,11). However, the genes coding for these prenyltransferases (PTs) have not been isolated from plants until the recent identification of the first flavonoid-specific PT gene SfN8DT-1 from Sophora flavescens, which is responsible for the specific prenylation of the flavanone naringenin at C-8 (12). Subsequently, further progress was achieved in the molecular biological and biochemical investigations of plant FPTs, and additional relevant genes have been cloned and functionally characterized (Table  1). Moreover, due to their substrate promiscuity in vitro, some soluble types of aromatic PTs from microbes have been demonstrated to be capable of catalyzing the prenylation of flavonoids; e.g. 7-DMATS, a member of DMATS superfamily from Aspergillus fumigatus (18,19), and some PTs of the CloQ/ NphB group from the actinomycetes Streptomyces (20,21).
However, to date all of the functionally characterized plant FPTs have been cloned from Leguminosae ( Table 1). The genes, phylogeny (in relation to the reported FPTs), and substrate spectra (specificity or promiscuity) of the non-Leguminosae FPTs remain enigmatic. These questions prompted us to search for such enzymes.
Moraceous plants are rich sources of isoprenoid-substituted phenolic compounds and their Diels-Alder-type adducts (22,23). In our previous work cell suspension cultures of two closely related moraceous trees, Morus alba (mulberry) and Cudrania tricuspidata, were employed to prenylate different types of flavonoids, and the same bioconversions were observed when microsomes of the cell cultures were used (24). Further chemical investigations confirmed that the prenylated flavonoids formed the main chemical compositions of both of the two plant cell suspension cultures and that Diels-Alder-type adducts were present in large amounts in the cell suspension cultures of M. alba (25,26). These observations indicated that FPTs existed in the cell suspension cultures of each of the two plants. Therefore, in our ongoing work on plant FPTs (14,15), cell suspension cultures of M. alba and C. tricuspidata were selected as materials from which to isolate non-Leguminosae FPTs.
In this report we describe the characterization of two new FPTs, MaIDT and CtIDT, from the cell suspension cultures of M. alba and C. tricuspidata, respectively. As the first identified non-Leguminosae FPTs, these FPTs share fairly low identities with their homologs in Leguminosae and feature catalytic regioselectivity and promiscuity. These findings not only imply the diversity of FPTs in the plant kingdom but may also illuminate the discoveries of other FPTs in non-Leguminosae plants.

EXPERIMENTAL PROCEDURES
General Experimental Procedures-1 H and 13 C NMR spectra were recorded on Varian NMR System 600 spectrometers (Varian Inc., Palo Alto, CA). The HPLC-UV/ESI-MS analyses of the enzymatic products were performed as previously described (15) with the exception that the solvent system consisted of a linear gradient from 50% to 100% (v/v) methanol in water with 0.1% formic acid over a period of 30 min followed by an isocratic elution with 100% methanol for 10 min and that the UV detector was set at 350 or 265 nm. For the isolation of the enzymatic products, an approach similar to that previously described was employed (15) with the exception that the semi-preparative reverse-phase HPLC was performed using a linear gradient from 50% to 100% (v/v) methanol in water over a period of 20 min.
Isolation of cDNAs Homologous to AtVTE2-1 from M. alba-Ten-day-old M. alba cell cultures were treated with 0.1 mM methyl jasmonate for 20 h, and total RNA was subsequently extracted with an E.Z.N.A. TM Plant RNA kit (Omega Bio-tek Inc., Doraville, GA) and reverse-transcribed (RT) using a SMARTer TM RACE cDNA Amplification kit (Clontech Inc., Mountain View, CA). The RT products were subjected to rapid amplification of cDNA ends (RACE) according to the manufacturer's protocol. The 3Ј ends of MaIDT was obtained using the primer expressed sequence tag (EST)-3ЈGSP (5Ј-TGATGCCG-ATATTGACAGGATAAATAAGCC-3Ј) that was specific for an EST sequence of the mulberry root (GenBank TM accession  GT734921), which is homologous to AtVTE2-1 (Gen-Bank TM accession number AY089963). The 5Ј ends of MaIDT were re-isolated by RT-PCR using the primer MaIDT-5ЈGSP (5Ј-AGAGGAGGAGCAGAATAGAAGTGGGCCA-3Ј). Its full-length clone was acquired by RT-PCR using the gene-specific primer pairs MaIDT-Fw (5Ј-GAATTCATGGAGCTCTC-AATCTCTCACTCT-3Ј) and MaIDT-Rv (5Ј-GCGGCCGCT-TATATGAAAGGAAATATGACAAACTCCA-3Ј). In these sequences the restriction sites for subcloning are underlined. The PCR products were cloned into pEASY-Blunt Simple vector (TransGen Biotech Co., Ltd., Beijing, China) for sequencing and then subcloned into pESC-HIS vector (Stratagene Inc., La Jolla, CA), which resulted in the expression construct pESC-HIS-MaIDT.
Heterologous Expression in Yeast-The expression vectors pESC-HIS-MaIDT and pESC-HIS-CtIDT were introduced into the yeast strain YPH499. The microsomal fractions of YPH499 expressing pESC-HIS-MaIDT and pESC-HIS-CtIDT were prepared as previously described (14,15) and resuspended in 100 mM Tris-HCl (pH 9.0). The total protein concentration was determined by the Bradford method (30).
Prenyltransferase Activity-For the quantitative determination of the PT activity, the basic reaction mixture (100 l) contained 100 mM Tris-HCl (pH 9.0), 10 mM MgCl 2 , 200 M prenyl acceptor, 400 M DMAPP, and 100 g of protein of the recombinant yeast microsome. The reaction mixtures were incubated at 30°C, and the reactions were terminated by the addition of 200 l of methanol. The protein was removed by centrifugation at 14,000 ϫ g for 20 min. The enzymatic products were analyzed by HPLC-UV/ESI-MS under the conditions described above. For the quantitative measurements of the enzyme activity, three parallel assays were carried out routinely, and the values were obtained in duplicate for each assay.
Biochemical Properties of the Recombinant Enzymes-The assays for the determinations of the kinetic parameters (100 l) of the prenyl acceptors contained 400 M DMAPP, 100 g of protein of the recombinant yeast microsome, and prenyl acceptors at final concentrations of 5, 10, 20, 40, 80, 160, and 400 M. For the determinations of the kinetic parameters of DMAPP, isoliquiritigenin (1) at 400 M and DMAPP at final concentrations of 5, 10, 20, 40, 80, 160, and 400 M were used. The incu-bation time was 30 min. The apparent K m values were calculated from Lineweaver-Burk plots using the Hyper32 software.
Preparative Synthesis of the Enzymatic Products for Structural Elucidation-The assays for the isolation of the enzymatic products (10 -15 ml) contained 100 mM Tris-HCl (pH 9.0), 10 mM MgCl 2 , 1 mM DMAPP, 500 M prenyl acceptors, and 30 mg of protein of the recombinant yeast microsome. The reaction mixtures were incubated at 30°C for 16 h and subsequently extracted with ethyl acetate (20 ml ϫ 5). After evaporation of the solvent, the residues were dissolved in methanol and purified by reverse-phase semi-preparative HPLC under the conditions described above. The isolated products were subjected to MS and 1 H and 13 C NMR spectroscopic analyses, which yielded the following results.  Fig. S10) (38).
Computer-assisted Sequence Analysis-The sequence identities were obtained via alignment of the amino acid sequences performed with the program BLAST2 SEQUENCES (www. ncbi.nlm.nih.gov). The PT proteins of the plants were aligned using ClustalW (39). From the alignments, a consensus phylogenetic tree was generated by the neighbor-joining method using MEGA6 (40). Bootstrap values are indicated in percentages (only those Ͼ70% are presented) on the nodes. The bootstrap values were obtained from 1000 bootstrap replicates. The scale bar corresponds to 0.1 estimated amino acid changes per site. The exon-intron structure diagram was generated using the online Gene Structure Display Server (GSDS) and the exon position and gene length method.
Nucleotide Sequence Accession Number-The nucleotide sequences of MaIDT and CtIDT have been deposited in the GenBank TM database under the accession numbers KM262659 and KM262660, respectively.

Cloning of MaIDT from the Cultured Cells of M. alba-
According to previous studies (12)(13)(14)(15)(16)(17), plant FPTs from leguminous plants are closely related to each other and cluster together in the phylogenetic tree. Therefore, we initially hypothesized that the FPTs of M. alba would cluster together with the leguminous FPTs and used degenerated primer pairs based on these FPTs to clone the candidate gene(s). Similar cloning strategies proved to be successful in the isolation of SfFPT and GuA6DT from S. flavescens and Glycyrrhiza uralensis in our previous work (14,15). However, after RT-PCR and RACE, we failed to achieve active PTs. A similar cloning strategy also failed to lead to the identification of active FPT from Epimedium acuminatum (41). Hence, we switched to the Morus EST library for the cloning of these candidate genes. Utilizing the internal sequence of an EST clone that is homologous to the homogentisate PT AtVTE2-1, we succeeded in amplifying the full-length cDNA of MaIDT, and PT activity was observed after the heterologous expression of MaIDT in yeast.
The encoded polypeptide of MaIDT (402 amino acids) had seven putative transmembrane ␣-helices as predicted by the TMHMM 2.0 program. The polypeptide possessed putative transit peptide sequences for targeting plastids as predicted by Chlorop 1.1, TargetP 1.1, and iPSORT. A multiple alignment of the encoded polypeptide and the reported leguminous FPTs is shown in Fig. 1 are not yet known. Moreover, the polypeptide shared identities of only 24 -30% at the amino acid level with the reported leguminous FPTs, although the above data fulfilled the three criteria for selecting FPT candidates that have been adopted in a previous report (12).
Enzymatic Activity of MaIDT and Identification of the Enzymatic Products-In previous studies mulberry trees and related plants have been found to be rich sources of prenylated flavonoids (22,23), and prenylation activity has been confirmed in microsomal fractions of mulberry tree cell cultures (42,43). Accordingly, various flavonoids, including chalcones, flavones, flavanones, and isoflavones have been employed as acceptors to determine enzymatic activities. Of these, isoliquiritigenin (1) is accepted to have the highest efficiency. HPLC-UV/ESI-MS analysis revealed that a product at 24.4 min with a 68-atomic mass unit higher shift in its molecular peak was detected in the incubation mixture with the microsomes from the yeast cells transformed with MaIDT but not with an empty vector (Fig. 2). Product formation was strictly dependent on the presence of MaIDT, DMAPP, Mg 2ϩ , and 1. A linear dependence of the product formation on the amount of microsome protein was found up to 200 g per 100-l assay and on the reaction time was up to 30 min at 30°C.
To clarify the prenylation properties, the enzymatic product (1a) was subsequently prepared and subjected to NMR analysis. In the 1 H NMR spectrum, the presence of the prenyl moiety was suggested by the appropriate signals for two methyls (␦ H 1.64 and ␦ H 1.78, s, 3H each), a methylene (␦ H 3.37, d, J ϭ 7.2 Hz, 2H), and a methine group (␦ H 5.28, t, J ϭ 7.2 Hz, 1H). The location of the prenyl moiety was indicated by the disappearance of the H-3Ј signal (1H at ␦ H 6.35) and the change in the H-5Ј peak shape, which became an ortho-coupled doublet (␦ H 6.53, d, J ϭ 8.4 Hz) (supplemental Fig. S1). This observation is supported by a downfield chemical shift of C-3Ј from ␦ C 103.7 in 1 to ␦ C 114.5 in 1a (supplemental Fig. S2). The NMR data of 1a were in good agreement with those of 3Ј-dimethylallyl isoliquiritigenin (31,32). Therefore, the enzymatic product of MaIDT was unambiguously identified as 3Ј-dimethylallyl isoliquiritigenin (1a), which indicates that this enzyme regiospecifically transferred a dimethylallyl moiety to 1 at the C-3Ј position; thus, this recombinant enzyme was designated as MaIDT (M. alba isoliquiritigenin 3Ј-dimethylallyltransferase).
Cloning and Identification of CtIDT from Cultured Cells of C. tricuspidata-To increase understanding of the diversity of the FPTs in moraceous plants, a cell suspension culture of C. tricuspidata was chosen for further study. Due to the absence of genetic information about C. tricuspidata, homologous PCR was used to clone candidate gene(s) as previously described (14,15). Similarly, degenerated primer pairs based on the conserved regions of the leguminous FPTs failed to produce  DECEMBER 26, 2014 • VOLUME 289 • NUMBER 52 an active FPT. Next, we designed degenerated primer pairs based on the conserved amino acid regions of MaIDT and its closely related genes that are obtained via a BLASTN search of GenBank TM and eventually cloned a full-length cDNA possessing FPT activity (CtIDT). The polypeptide encoded by CtIDT (398 amino acids) fulfilled the three criteria for selecting FPTs (12). This polypeptide shared sequence identities of 66% with MaIDT and only 24 -30% with the reported leguminous FPTs. Enzymatic assays containing various flavonoids and HPLC-UV/ ESI-MS analyses revealed that CtIDT exhibited a function similar to that of MaIDT (Fig. 3, A and B); thus, this recombinant enzyme was designated as CtIDT (C. tricuspidata isoliquiritigenin 3Ј-dimethylallyltransferase).

Moraceous Flavonoid Prenyltransferases
Substrate Specificities of MaIDT and CtIDT-To investigate the substrate specificities of the moraceous FPTs, we tested the activities of MaIDT and CtIDT via HPLC-UV/ESI-MS analysis utilizing 18 potential substrates that included various types of flavonoids and structurally simple phenols (Fig. 3, A and B). In addition to isoliquiritigenin (1), MaIDT and CtIDT were also able to catalyze the prenylation of three other chalcones (2Ј,4Јdihydroxychalcone (2), 2,4,2Ј,4Ј-tetrahydroxychalcone (3), and butein (4)) and two isoflavonoids (genistein (7), 2Ј-hydroxygenistein (8)), and the flavone type substrate apigenin (11) was only accepted by MaIDT. The enzymatic products were prepared from larger scale assays and subjected to MS and NMR spectroscopic analyses with the exception of apigenin (11) (supplemental Fig. S3-S10), the prenylated product of which was identified as 6-dimethylallylapigenin (11a) in a comparison with an authentic standard using a HPLC-UV/ESI-MS analysis (data not shown). The results clearly showed that although the three types of flavonoids could be accepted by MaIDT and CtIDT, one dimethylallyl moiety was regioselectively introduced into the A ring of the skeletons at the same position (C-3Ј of the chalcones and C-6 of the isoflavones and flavones).
The relative reactivities of MaIDT and CtIDT toward different flavonoids also provided insight into the "substrate structure-enzyme selectivity relationships." MaIDT preferred 2Ј,4Ј-dihydroxychalcones to isoflavonoids and flavone. The existence of hydroxyls at C-2Ј,4Ј (C-5,7 for isoflavones) of the A ring was crucial for prenylation because the lack of any one of these hydroxyl groups led to no enzymatic product generation (5, 6, and 9). Regarding the B ring, the presence of a hydroxyl group at C-4 position increased the prenylation efficiency, as the conversion rate of 1 by MaIDT was higher than that of 2 (Fig. 3, A and B). In contrast, the hydroxyl group at the C-2/3 (C-2Ј/3Ј for isoflavones) position of B ring reduced the prenylation efficiency as revealed by comparisons of the relative yields of 3 and 4 to 1 and of 8 and 10 to 7 (Fig. 3, A and B).
The prenylation catalyzed by FPTs represents a Friedel-Crafts alkylation of flavonoid skeleton. A carbocation is generated from DMAPP and carries out an electrophilic attack at flavonoid ring, leading to the regio-specific C-C bond formation. Hence, the electron density around attacked carbon at flavonoid ring might play an important role in the course of prenylation. For example, the insufficient electron density around C-6 of daidzein (9) without C-5 hydroxyl group compared with genistein (7) might lead to the non-occurrence of prenylation.
Regarding CtIDT, its substrate specificity was slightly different from that of MaIDT; first, CtIDT accepted isoflavones with higher efficiency than did MaIDT, and this finding is partially in accordance with the fact that 7a has been identified in large amounts in C. tricuspidata cell suspension cultures (26); second, the hydroxyl group at C-2 of chalcones decreased the activity of CtIDT as revealed by comparisons of the relative yields of 3 to 1, whereas the corresponding hydroxyl group at the C-2Ј position of the isoflavone increased the activity as revealed by the comparison of the relative yields of 8 to 7.
Other prenyl donors (including GPP, farnesyl diphosphate, geranylgeranyl diphosphate, and phytyl diphosphate) were also tested in the enzymatic assay. The results revealed that GPP could also act as prenyl donor when 1 is used as the acceptor (Fig. 3C).

Biochemical Properties of the Recombinant MaIDT and
CtIDT-Previous studies have reported that the plant FPT activity is dependent on divalent cations and strongly affected by pH and temperature (10). Further investigations using 1 as an acceptor and DMAPP as a prenyl donor in this study revealed that the highest activities occurred at ϳ30°C and that the activity decreased rapidly at temperatures above 40°C (Fig.  4A). The analysis of the enzyme activity within the pH range of 6.0 to 11.0 revealed that the optimal pH value was between 8.0 and 9.0 (Fig. 4B). MaIDT activity was observed to decrease in the order of Mg 2ϩ Ͼ Ba 2ϩ Ͼ Ca 2ϩ Ͼ Mn 2ϩ Ͼ Fe 2ϩ Ͼ Ni 2ϩ , and Co 2ϩ and Cu 2ϩ did not lead to the production of 1a; CtIDT activity was found to decrease in the order of Mg 2ϩ Ͼ Mn 2ϩ Ͼ Ca 2ϩ Ͼ Fe 2ϩ Ͼ Ba 2ϩ , whereas Ni 2ϩ , Co 2ϩ , and Cu 2ϩ did not lead to the production of 1a (Fig. 4C). No product was detected in the buffer after the addition of EDTA. Phylogenetic and Gene Structure Analysis of the Flavonoid Dimethylallyltransferases and Related Enzymes-A neighborjoining phylogenetic tree was constructed to analyze the evolutionary relationships of MaIDT and CtIDT with the other plant aromatic PTs (Fig. 5A). As the first identified FPTs from non-Leguminosae plant species, MaIDT and CtIDT will have particularly significant effects on our understanding of the molecular evolution of plant FPTs. The two moraceous FPTs MaIDT and CtIDT appeared to cluster together to form a new branch with a common origin in the homogentisate phytyltransferases for plastoquinone formation (AtVTE2-2 and GmVTE2-2). This branch is clearly divergent from that of the leguminous FPTs, which were recruited from the vitamin E biosynthetic pathway. Considering the similar substrate spectra and sequence divergence of the leguminous and moraceous FPTs, it can be concluded that they have evolved independently via recruitment from the vitamin E and plastoquinone biosynthetic pathways, respectively. The sequence distance might reflect the taxonomical distance between Moraceae from Leguminosae plants (44). This situation would represent a somewhat subtler example of convergent evolution. Moreover, this type of convergent evolution has been termed "repeated evolution," which indicates situations in which "a similar function arose not from a completely unrelated sequence but from a homologous, although not orthologous gene" (45,46). Consequently, the search for new FPTs should not focus solely on the degrees of sequence identity with known enzymes. Furthermore, in evo-lutionarily distant plants, it is possible that as-yet unidentified FPTs might share low homology with the known FPTs.
Based on a BLASTP search of the Morus Genome Database (MorusDB) (47,48), 1 hit with MorusDB accession number Morus014065 was found and shared an identity of 97% with MaIDT, and its genome structure was compared with those of the leguminous FPT GmG4DT and homogentisate PTs involved in vitamin E and plastoquinone biosyntheses, which are available in GenBank TM . The exon-intron structure of MaIDT was quite different from those of the leguminous FPT GmG4DT and the related homogentisate PTs involved in vitamin E biosynthesis, whereas the exon-intron organizations of MaIDT and the homogentisate PTs for plastoquinone biosynthesis were very similar (Fig. 5B). These observations confirm the results of the phylogenetic analysis, which suggested that the leguminous and moraceous FPTs were derived from homogentisate PTs involved in vitamin E and plastoquinone biosynthesis, respectively, and have evolved independently.

Prenylation as a Typical Modification of Aromatic Compounds in Moraceous
Plants-A large variety of biological activities has been reported for prenylated aromatic compounds, and these natural products are attractive resources in the food and pharmaceutical industries. Aromatic PTs catalyze transfer reactions of prenyl moieties from prenyl diphosphates to diverse aromatic acceptors. Theses PTs are involved in both general and specialized metabolism, play important roles in living organisms, and strongly contribute to the structural diversity and biological activities of aromatic natural products (7,49).
In moraceous plants, many prenylated aromatic compounds, such as Diels-Alder-type adducts, flavonoids, 2-arylbenzofurans, xanthones, and stilbenes, have been identified, and their biological/pharmacological activities have been extensively studied (9,22,23). Of these compounds optically active Diels-Alder-type adducts, which are characteristic mulberry constituents with promising biological activities, are attractive compounds from a biosynthetic point of view, and prenylation is a key biosynthetic step in their formation (8,9).
Evolution of FPTs-Plant aromatic PTs are divided into two distinct subgroups: one is the 4-hydroxybenzoate PT, which is involved in the biosyntheses of ubiquinone (benzoquinone) and shikonin (naphthoquinone) derivatives; the second subgroup is the homogentisate prenyltransferase (HGPT) family, which is involved in vitamin E and plastoquinone biosyntheses. Thus far, all the identified plant FPTs belong to HGPT family (11)(12)(13)(14)(15).
Based on the newly identified moraceous FPTs, the phylogenetic tree of the HGPT gene family was rebuilt. A clear phylogenetic separation of the HGPT genes based on primary metabolism (i.e. those involved in tocopherol, tocotrienol, and plastoquinone biosyntheses) versus specialized metabolism can be observed. The formers phylogenetic group is clustered in a pattern of functional clades because the molecular evolution of this group has been severely constrained by natural selection. However, the latter group exhibits a remarkable flexibility of the HGPT family to evolve enzymes with new substrate specificities. Moreover, several functionally divergent enzymes, such as FPTs, phloroglucinol PTs, and coumarin PTs, have emerged from the primary metabolic HGPTs (50 -53). Interestingly, repeated evolution has been observed to exist in the FPTs. The moraceous FPTs MaIDT and CtIDT are grouped in a branch that is clearly different from the branch in which the leguminous FPTs cluster, which implies that FPTs evolved independently in the two plant lineages. Repeated evolution might represent the majority of convergent evolution cases in plantspecialized metabolism, and many such cases have been identified (46); e.g. in each plant lineage the respective linalool synthase is more similar to other terpene synthases, which are responsible for the synthesis of other monoterpenes within that lineage than it is to the linalool synthases from different lineages (54,55). Furthermore, the phylogeny of FPTs seems to exhibit a pattern of species-specific expansion in Leguminosae and Moraceae. Apparently, the identified moraceous FPTs and leguminous FPTs arose individually by repeated gene duplications and divergence that occurred after the split of the two plant families. However, our understanding of the plant FPTs is currently based on just over 10 genes so far; thus, different FPTs need to be identified to clarify the molecular evolution of plant FPTs.
In some studies the abovementioned species-specific patterns of repeated gene duplications and sequence divergence are referred to as "blooms" (56). In practice, the approach of selecting candidate genes from species-specific blooms has proven successful for the discovery of new genes (56,57). As FPTs have now been successfully identified in moraceous plants, this approach is likely worthwhile for attempts to clone other FPTs from moraceous plants. Given the great variety of prenylated flavonoids that have been isolated from moraceous plants, a large family of FPTs can be imagined. Moreover, prenylated flavonoids are active constituents of many other non-Leguminosae medicinal plants. As the first successful identification of FPTs from non-Leguminosae plants, this study might inspire the identification and characterization of additional and evolutionarily varied FPTs that are involved in the biosyntheses of pharmacologically active prenylated flavonoids in plants other than those of the Leguminosae and Moraceae families.
In summary, we have characterized the first two non-Leguminosae plant FPTs MaIDT and CtIDT from M. alba and C. tricuspidata, respectively. These FPTs feature catalytic regioselectivity and promiscuity. Multiple alignments, phylogenetic analyses and gene structure analyses revealed that these novel FPTs evolved independently from their counterparts in Leguminosae. These findings provide useful guidance for identifying additional and evolutionarily varied FPTs in non-Leguminosae plants. Furthermore, based on their substrate promiscuities and high catalytic efficiencies, MaIDT and CtIDT might be useful as an environmentally friendly and efficient biological catalysts of the regiospecific prenylation of flavonoids to produce bioactive compounds for potential therapeutic applications.