Next Article in Journal
Synthesis and Biological Activity of Novel (Z)- and (E)-Verbenone Oxime Esters
Next Article in Special Issue
Headspace Solid-Phase Microextraction and Ultrasonic Extraction with the Solvent Sequences in Chemical Profiling of Allium ursinum L. Honey
Previous Article in Journal
Diallylthiosulfinate (Allicin), a Volatile Antimicrobial from Garlic (Allium sativum), Kills Human Lung Pathogenic Bacteria, Including MDR Strains, as a Vapor
Previous Article in Special Issue
Gypmacrophin A, a Rare Pentacyclic Sesterterpenoid, Together with Three Depsides, Functioned as New Chemical Evidence for Gypsoplaca macrophylla (Zahlbr.) Timdal Identification
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 22
Issue 10
10.3390/molecules22101706
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization, Function, and Transcriptional Profiling Analysis of 3-Hydroxy-3-methylglutaryl-CoA Synthase Gene (GbHMGS1) towards Stresses and Exogenous Hormone Treatments in Ginkgo biloba

by
Xiangxiang Meng
,
Qiling Song
,
Jiabao Ye
,
Lanlan Wang
and
Feng Xu
*
College of Horticulture and Gardening, Yangtze University, Jingzhou 434025, Hubei, China
*
Author to whom correspondence should be addressed.
Molecules 2017, 22(10), 1706; https://doi.org/10.3390/molecules22101706
Submission received: 26 August 2017 / Accepted: 8 October 2017 / Published: 12 October 2017
(This article belongs to the Special Issue Diversity of Terpenoids)
Browse Figures
Versions Notes

Abstract

:
3-Hydroxy-3-methylglutaryl-CoA synthase (HMGS) is one of the rate-limiting enzymes in the mevalonate pathway as it catalyzes the condensation of acetoacetyl-CoA to form 3-hydroxy-3-methylglutaryl-CoA. In this study, A HMGS gene (designated as GbHMGS1) was cloned from Ginkgo biloba for the first time. GbHMGS1 contained a 1422-bp open-reading frame encoding 474 amino acids. Comparative and bioinformatics analysis revealed that GbHMGS1 was extensively homologous to HMGSs from other plant species. Phylogenetic analysis indicated that the GbHMGS1 belonged to the plant HMGS superfamily, sharing a common evolutionary ancestor with other HMGSs, and had a further relationship with other gymnosperm species. The yeast complement assay of GbHMGS1 in HMGS-deficient Saccharomyces cerevisiae strain YSC6274 demonstrated that GbHMGS1 gene encodes a functional HMGS enzyme. The recombinant protein of GbHMGS1 was successfully expressed in E. coli. The in vitro enzyme activity assay showed that the kcat and Km values of GbHMGS1 were 195.4 min−1 and 689 μM, respectively. GbHMGS1 was constitutively expressed in all tested tissues, including the roots, stems, leaves, female flowers, male flowers and fruits. The transcript accumulation for GbHMGS1 was highest in the leaves. Expression profiling analyses revealed that GbHMGS1 expression was induced by abiotic stresses (ultraviolet B and cold) and hormone treatments (salicylic acid, methyl jasmonate, and ethephon) in G. biloba, indicating that GbHMGS1 gene was involved in the response to environmental stresses and plant hormones.

1. Introduction

Ginkgo biloba L., commonly called the maidenhair tree, is the oldest relict plant species in the world and is considered a “living fossil” [1]. Many bioactive constituents in the leaf extract of G. biloba are beneficial to human health; these components include flavonoids and terpene trilactones (TTLs including ginkgolides and bilobalide), which are proven to reduce blood pressure, dilate the peripheral blood vessels, discourage clot formation, promote blood circulation and cerebral metabolism, reinforce the walls of the capillaries, and protect nerve cells from harm during oxygen deprivation [2]. Flavonoids are the preferred natural medications for the treatment of dementia and can be obtained from many other plants, whereas TTLs are unique components of G. biloba [3] that play a key role in the active ingredients. Ginkgolide B has a specific antagonistic effect on platelet-activating factor [4], whereas bilobalide has a protective function against neuronal damage and can be used to treat demyelinating brain, spinal cord, and neurological diseases [5,6]. However, the dry leaves of G. biloba have very low content of TTLs [7]. The unique chemical structure of TTLs renders them difficult to synthesize through chemical and biological methods, therefore increasing TTL content in G. biloba is of important significance in the ginkgo industry. Many studies have shown that the key to increasing the content of TTLs is to increase the expression level of the key genes in their metabolic pathways [8,9,10].
Terpenoids possess multiple ecophysiological functions, such as regulating plant growth and development, resisting various stresses, and providing direct and indirect plant defenses [11]. Considering their diverse functions, TTLs have been extensively applied as important raw materials in medicine, food, cosmetics, spices, and health care [12]. The terpenoids of plants, including TTLs, are synthesized mainly by two pathways: the classical cytosolic mevalonate (MVA) pathway, which facilitates the synthesis of sesquiterpenoids, triterpenoids, and steroids; and the plastidial methyl-erythritol-4-phosphate (MEP) pathway, which is MVA-independent and normally facilitates the formation of monoterpenoids, diterpenoids, and tetraterpenoids [13]. As a sesquiterpene, bilobalide is thought to be derived from farnesyl diphosphate (FPP) in the MVA pathway, but more authors suggested that bilobalide is most likely a product of partially degraded ginkgolide [14], whereas ginkgolide biosynthesis is thought to be achieved by the MEP pathway. However, a cross-talk seemed to occur between the two biosynthetic pathways because the universal five carbon skeleton–isopentenyl diphosphate (IPP) of terpenoid biosynthesis is derived from both MVA and MEP pathways [15]. As the second catalyzing enzyme, 3-hydroxy-3-methylglutaryl coenzyme-A synthase (HMGS) catalyzes the irreversible conversion of acetoacetyl-CoA into 3-hydroxy-3-methylglutaryl-CoA (HMG-CoA), which is the committed step in the MVA pathway of terpenoid biosynthesis in plants. 3-Hydroxy-3-methylglutaryl-CoA reductase (HMGR) catalyzes HMG-CoA to MVA and further converts it into isopentenyl pyrophosphate (IPP) and its isomer dimethylallyl pyrophosphate (DMAPP), which are the common precursors of many types of terpenoid biosynthesis. Afterward, IPP and DMAPP as the substrates are catalyzed by their corresponding enzymes to synthesize the following intermediates: geranyl diphosphate, FPP, and geranylgeranyl diphosphate, which are further catalyzed by the corresponding terpene synthase to produce all kinds of terpenoids [16,17] (Figure 1). Accumulated evidence suggests that HMGS plays an important role in the plant defense mechanism by participating in the MVA pathway, and the synthesis of plant terpenoids also positively correlates with the activity of HMGS. For example, a previous study reported that a positive correlation exists between HMGS activity and the dry rubber content of the latex from rubber trees [18]. Ren et al. [19] also stated that the ganoderic acid content was approximately 15.1%–24.2% higher than that of the control group, whereas the GlHMGS gene was overexpressed in Ganoderma lucidum. Another transgenic study also indicated that the overexpression of BjHMGS in Arabidopsis thaliana increases sterol concentration and enhances Botrytis cinerea resistance [20]. Therefore, we speculated that the GbHMGS1 gene may play an important regulatory role in the resistance mechanism to abiotic stresses as well as in the biosynthesis of sesquiterpenoids (e.g., farnesene and bisabolene) in G. biloba [21].
Each gene in the MEP and MVA pathways that contributes to the biosynthesis of isoprenoid should be identified and characterized to determine the overall biosynthetic pathway of terpenoids in G. biloba. To date, many genes involved in the MEP pathways are isolated from G. biloba, such as 1-deoxy-d-xylulose 5-phosphate synthase (DXS) [22], 1-deoxy-d-xylulose 5-phosphate reducto-isomerase [23], 2-C-methyl-d-erythritol 4-phosphate cytidyltransferase [24], 4-(cytidine 5’-diphospho)-2-C-methyl-d-erythritol kinase (CMK) [8], 2-C-methyl-d-erythritiol 2,4-cyclo-diphosphate synthase [25], 1-hydroxy-2-methyl-2-(E)-butenyl-4-diphosphate synthase [26], and 1-hydroxy-2-methyl-2-(E)-butenyl-4-diphosphate reductase [27]. However, the presence of IPP cross-talk further complicates the origin of TTL building block in G. biloba [28]. Therefore, the MVA pathway may also have a minimal contribution to the biosynthesis of TTLs. Up to now, some of genes encoding enzymes involved in MVA pathway have been cloned and characterized from G. biloba, including HMGR [29], mevalonate diphosphate decarboxylase (MVD) [10], mevalonate kinase (MVK) and acetyl-CoA C-acetyltransferase (AACT) [30]. However, the HMGS genes in G. biloba are poorly studied. In this study, we cloned a HMGS gene (GbHMGS1) from G. biloba for the first time. The GbHMGS1 gene was characterized by prokaryotic expression, and yeast complement assay were performed to confirm the function of GbHMGS1. Furthermore, the expression profiling of GbHMGS1 in different tissues and under abiotic stresses (cold and ultraviolet B (UV-B)) and hormonal treatments (salicylic acid (SA), methyl jasmonate (MJ), and ethephon (Eth)) were investigated by qRT-PCR analysis. These findings will be useful supplements for unveiling the overall biosynthetic pathway as well as for the large-scale production of isoprenoid with pharmacological properties of this living fossil.

2. Results

2.1. Identification and Characterization of GbHMGS1

The conserved sequence composed of 542 bp was obtained by PCR using primers designed from the homologous regions of the reported HMGS proteins. Based on this sequence, a putative HMGS gene was obtained from G. biloba by RACE method. The full-length cDNA sequence of GbHMGS1 was 1945 bp and contained a 1422-bp ORF encoding a 474 amino acid protein. A 98 bp 5′-UTR was upstream of the start codon, and the coding region was followed by a 425 bp 3′-UTR downstream from the stop codon (Supplementary Material Figure S1). A BLASTn search on NCBI showed that the cDNA sequence of GbHMGS1 was highly similar to those of HMGS genes from other plants, indicating that the cloned gene was a member of the HMGS family. Therefore, the gene was designated as GbHMGS1 (GenBank accession number: MF281039).

2.2. Bioinformatic Analysis of GbHMGS1 Protein

ExPASy online (http://web.expasy.org/compute_pi/) analysis results showed that the pI and theoretical molecular weight of the GbHMGS1 protein were 5.23 and 52.67 kDa, respectively, which are similar to those of previously reported plant HMGSs [31,32,33,34]. A BLASTp search in GenBank database (https://blast.ncbi. nlm.nih.gov) showed that deduced GbHMGS1 protein exhibited 74%–85% similarity to HMGS proteins of many other plants (Supplementary Material Table S2); specifically, 85%, 83%, 79%, 78%, 77%, 75% and 74% similarity to the counterparts of Taxus media, Pinus sylvestris, Narcissus tazetta, Sorghum bicolor, Theobroma cacao, Hevea brasiliensis, and Camptotheca acuminata, respectively, which further verified that GbHMGS1 belonged to plant HMGS family.
The multiple alignments of the deduced GbHMGS1 protein sequences with other HMGSs using software Vector NTI11.5 revealed that GbHMGS1 was highly homologous to HMGSs from other plants. Further sequencing analysis showed that GbHMGS1 protein contained the conserved motif “NxD/NE/VEGI/VDx(2)NACF/YxG” (Figure 2, marked with a red box) and five conserved sites (amino acids Glu86, Cys120, Ser251, Gly328, and Ser362) (Figure 2, marked with “asterisk”), which are crucial for HMGS activity [35]. These results indicated that GbHMGS1 catalyzes a similar enzymatic reaction to other plant HMGSs.
A comparative modeling of the 3D structure of GbHMGS1 (Figure 3) was generated based on the highest query coverage of the template Brassica juncea HMGS (2f82.1.A) with SWISS-MODEL [36] to further elucidate the GbHMGS1 protein. Similar to those of other HMGSs, the predicted 3D model of GbHMGS1 consisted of two structural regions called the lower and upper regions [35]. In addition, the interface of the upper and lower regions defined the acetoacetyl-CoA-binding sites [37], and all the conserved motifs and active sites were localized in a five-layered core structure of the upper region.
A phylogenetic tree was constructed using the amino acid sequences of GbHMGS1 and HMGS protein sequences from other species to investigate the sequence diversity and evolutionary relationships among HMGS proteins. As shown in Figure 4, the phylogenetic tree was divided into four branches of Monocotyledons, Dicotyledons, Gymnospermae, and Algae. GbHMGS1 belonged to the Gymnospermae cluster and had the closest relationship with HMGSs from T. media and P. sylvestris, and they all belong to the Gymnospermae, which is in accordance with the species classification of plants. These relationships reflect the evolutionary conservation and evolutionary diversity of plant HMGS.

2.3. Prokaryotic Expression and In Vitro Enzyme Activity Analysis of GbHMGS1

To express GbHMGS1 in E. coli, we constructed the pET32a-GbHMGS1 expression vector with the T7 promoter and a His-tag by ligating the coding region of GbHMGS1 into pET-32a(+). The expression construct was evaluated for in-frame fusion through restriction enzyme digestion and DNA sequencing. Undergoing IPGE induction, GbHMGS1 was expressed as a major protein product in the insoluble, soluble and total cellular protein (Supplementary Material Figure S2). The molecular weight of the expressed recombinant protein with the His-tag was estimated to be about 52.67 kDa, coinciding with the molecular weight of GbHMGS1 protein predicted through bioinformatics. Purification of the recombinant GbHMGS1 protein was achieved by Nickel-CL agarose affinity chromatography and used for enzyme activity assay. The rates of catalysis by purified GbHMGS1 were measured as a function of acetoacetyl-CoA concentrations to evaluate the kinetic parameters of GbHMGS1. The initial velocity V0 plotted against substrate concentration yielded a rectangular hyperbola with Vmax = 37.1 ± 5.4 nM·min−1 (Table 1). Sirinupong et al. [38] reported Km = 530 ± 50 μM and a specific activity of 41.3 μmol·min−1·mg−1 using acetoacetyl-CoA as substrate for a recombinant HMGS from H. brasiliensis. Similar to HbHMGS, the purified GbHMGS1 had a higher Km of 689 μΜ and a lower specific activity of 15.5 μmol·min−1·mg−1, indicating that GbHMGS1 possesses HMGS activity.

2.4. Functional Complementation of GbHMGS1 in Saccharomyces cerevisiae

To determine the function of GbHMGS1, we successfully constructed the expression vector pYES2-GbHMGS1. The GbHMGS1 cDNA driven by a galactose-dependent promoter was expressed in the HMGS-deficient haploid yeast strain YSC6274, and the mevalonate derived from MVA pathway is essential for yeast survival. Wild type S. cerevisiae strain YSC1020 could grow on either induction medium YPG with G418 (Figure 5A) or non-induction medium YPD with G418 (Figure 5B). By contrast, the transformed yeast YSC6274 strain harboring pYES2-GbHMGS1 could grow on YPG expression medium (Figure 5C) but not on the non-induction medium YPD with G418 (Figure 5D). The results showed that the expression of GbHMGS1 can rescue the functional defect of the hmgs knockout yeast by encoding the HMGS enzyme to catalyze the synthesis of mevalonic acid. This data further confirmed that GbHMGS1 have HMGS activity.

2.5. Differential Expression Profiles of GbHMGS1 in Various Tissues of G. biloba

As shown in Figure 6, the GbHMGS1 gene was constitutively expressed in all the tested tissues of G. biloba. In addition, GbHMGS1 was highly expressed in the mature leaves, followed by young fruits, male flowers, and roots, but was lowly expressed in the stems and hardly in the female flowers. These results revealed that GbHMGS1 had a preferential expression pattern in mature leaves.

2.6. GbHMGS1 Transcript Level and TTL Content Changes in G. biloba under the Induction of UV-B, Cold, MJ, Eth, SA Elicitor

A time-course expression analysis of the GbHMGS1 gene and TTL content under cold conditions, UV-B, MJ, Eth, and SA treatments was conducted on the ginkgo seedlings at 4–5 leaf stages to explore the roles of the GbHMGS1 genes in response to different abiotic stresses and hormonal treatments. As shown in Figure 7, all elicitor inductions significantly increased GbHMGS1 expression, and the TTL content was significantly enhanced after the corresponding various abiotic stresses and hormone treatments. Under UV-B, the transcript level of GbHMGS1 increased gradually until 16 h post-treatment and decreased at 24 h, however, the level was still significantly higher than that of the control. Afterward, the level gradually increased again, GbHMGS1 transcripts peaked at 48 h time point with about 2-fold increase compared with that of the control (Figure 7C). The TTL content was continuously induced and increased by 20% compared with that of the control after treatment for 48 h (Figure 7D). When exposed to the cold treatment, the expression level of the GbHMGS1 gene was slightly down-regulated during short-term cold treatment and dropped 3.4% lower than that of the control 2 days after treatment. However, when the treatment time was extended, the expression level was up-regulated and significantly increased to about 3.98-fold of the control for up to 8 days of post-treatment (Figure 7A). The TTL content did not change significantly in the initial period of treatment but still on an upward trend, and peaked at 8 days (11.16% compared with that of the control) (Figure 7B).
In response to the exogenous hormone treatments, Figure 7F shows that TTL content all continuously increased with the increasing treatment time and reached maximum at 8 days after MJ, SA, and Eth treatments, with about 20%, 14.14% and 15% increase than that of the control, respectively. When exposed to the Eth treatment (Figure 7E), the transcript level of GbHMGS1 was up-regulated within 6 days with the highest level of 458.79%, and then decreased over time. Figure 7E showed that the exogenous application of MJ caused a rapid increase in GbHMGS1 transcripts, a significant increase in GbHMGS1 expression at 4 days of post-treatment by 189.77%, a slight decrease thereafter, and a gradual increase again until the 8th day. The GbHMGS1 expression was effectively induced by SA (Figure 7E), the transcript level of GbHMGS1 continuously increased with the increasing treatment time and reached the highest level (181.61% as compared with that of the control) at 8 days of post-treatment.

3. Discussion

As a major rate-limiting enzyme in the terpenoid biosynthesis, HMGS converts acetoacetyl-CoA into HMGS-CoA in the MVA pathway. Furthermore, the HMGS gene family is related to the development and defense mechanism of plants [39]. Given its importance in MVA biosynthesis and stress resistance, HMGS received increasing attention and has been widely studied in angiosperms at the chemical, enzymological, and genetic levels [40,41,42]. However, few reports were found on HMGS enzymes or genes in gymnosperm species at the molecular level. The characterization and functional analysis of GbHMGS1 will be helpful to further understand its role in terpenoid biosynthesis in G. biloba.

3.1. GbHMGS1 is a Member of HMGS Gene Family

In the present study, we successfully cloned a HMGS gene from G. biloba and analyzed its function through in vitro activity assay and yeast complement assay. GbHMGS1 protein contained 474 amino acids and weighed 52.67 kDa; this finding coincides with the result of the SDS-PAGE profile in prokaryotic expression and verifies that the plant HMGS protein is generally composed of 460–500 amino acid residues and has a relative molecular mass of 50–60 kDa [31,32,33,34]. Multiple alignments showed that the deduced GbHMGS1 was highly identical to other plant HMGSs and has the highest homology with T. media because both plants belong to Gymnosperms, which was speculated to have similar HMGS functions. GbHMGS1 contained the conserved substrate binding motifs of HMGSs “NxD/NE/VEGI/VDx(2)NACF/YxG” and five active sites (amino acids Glu86, Cys120, Ser251, Gly328 and Ser362). The conserved motif “NxD/NE/VEGI/VDx(2)NACF/YxG” forms the core of the enzyme structure, which is localized at the entrance of the active site, and controls the catalysis of substrates by HMGS. Mutation of this motif reduces the catalytic activity of the enzyme or leads to the formation of abnormal products, which is conserved across species [35]. In addition, three conserved amino acid residues, namely, Cys120, His250 and Asn329 (residue numbering corresponds to GbHMGS1) were known to be active residues essential for the catalytic activity in HMGS [43]. The polypeptide chain of GbHMGS1 also contains three domains, namely, N-terminus, catalytic region, and C-terminus. The N-terminus contains a conservative signal peptide sequence and shows a high degree of similarity across most plants. The C-terminus contains an important catalytic cysteine residue that acts as a nucleophile in the first step of reaction to accelerate the formation of HMG-CoA [37]. HMGS proteins also have highly similar 3D structures consisting of the upper and lower regions. These findings have deepened our understanding of the correlation between the structure and function of GbHMGS1, which also provides important basic data for further studying the metabolic process of terpenoids in G. biloba. In vitro enzyme activity assay verified that GbHMGS1 possesses HMGS activity. Complement assay revealed that the expression of GbHMGS1 supplies the basic material for yeast survival, thereby confirming the catalytic function of GbHMGS1. Taken together, the results mentioned above suggested that GbHMGS1-encoded HMGS enzyme is a member of the plant HMGS protein family and has a similar catalytic activity of catalyzing the conversion of acetoacetyl-CoA into HMG-CoA in the MVA pathway in G. biloba.

3.2. Expression Pattern of GbHMGS1 in Different Tissues of G. biloba

The most unique components of G. biloba extracts are TTLs. Therefore, researchers have focused their studies on the TTLs of G. biloba. Strong evidence suggested that TTLs are synthesized in the roots and subsequently translocated to aerial parts through the phloem; the leaves only act as a storage location for ginkgolides and bilobalide [9,44,45]. Despite this progress, very little is known concerning the biosynthesis and identification of sesquiterpenes in G. biloba. Recently, two genes encoding (E,E)-farnesol and α-bisabolene synthases involved in sesquiterpene were cloned and functionally characterized from G. biloba. In future studies, whether the spatial expression profile of GbHMGS1 is positively correlated with other sesquiterpenes (e.g., farnesene and bisabolene) [21] content in different tissues of G. biloba will be investigated in light of the fact that HMGS played an important role in sesquiterpene biosynthesis.
Up to now, tissue-dependent differences in HMGS expression level have also been studied for many medicinal plants. The expression pattern of HMGS in plant tissues greatly varies across different plants. For example, HMGS is constitutively expressed, strongest in the leaves, moderate in the stems, and weakest in the roots of S. miltiorrhiza and C. acuminata [31,32]. In addition, TwHMGS exhibited the highest expression level in the stem, followed by the leaf, and the expression level in the root was the lowest [33]. A common feature is that the constitutive high expression level of HMGS was found in above-ground tissues, including the stems, leaves, and flowers but lowest in the roots. Similarly, the constitutive high-level expression of GbHMGS1 in ginkgo was found in the leaves and fruits because the aboveground organs function in terpenoid production and protection from solar radiation. On the contrary, the basic-level expression of GbHMGS1 was found in the roots, which might be because the underground organs function in terpenoid production to prevent microbial infection or nutritional stress. One unexpected observation was that GbHMGS1 had the lowest expression in female flowers. Alex et al. [46] studied the developmental expression pattern of HMGS and found that HMGS expression is the highest at early stages in the flower, seed and the seedling of B. juncea. Multiple homologous HMGS genes may participate in terpenoid biosynthesis in different tissues or during different developmental stages in one type of tissue. For instance, HMGS1 and HMGS2 from H. brasiliensis function as critical genes in the synthetic pathway of latex but are expressed in tissue-specific manner [38]. Our data revealed that GbHMGS1 was a tissue-specifically expressed gene, suggesting that GbHMGS1 may also participate in other unknown functions.

3.3. Transcript Level of GbHMGS1 and TTL Content Changes under Different Elicitor Treatments in Ginkgo Seedlings

Plant HMGS genes participate in the resistance of plants to adverse environmental stresses, and many plants exhibit a level of tolerance when exposed to extreme temperature for a period of time. A low temperature could up-regulates many genes, including GbPAL [47], GbANS [48], and GbFLS [49] that are involved in the flavonoid biosynthesis in G. biloba. However, few studies have investigated the relationship between cold stress and terpenoid biosynthesis. In the present study, GbHMGS1 was markedly up-regulated by cold treatment. The up-regulated expression pattern of GbHMGS1 by cold treatment confirmed its role in the response to low temperature stress.
Plants encounter a variety of abiotic stresses from the environment in their life. Among these stresses, UV-B stress is a critical environmental factor that threatens plant development and growth because the levels of UV-B on the surface of the earth increase continuously due to the decreased ozone layer. Plants have developed various metabolites, including UV-absorbing phenolic compounds and terpenoids, to protect against the deleterious effects of UV-B stress. Ultraviolet irradiation has been found to induce the expression of genes (e.g., CAD and CHI) [50,51]. However, few information concerning HMGS enzymes under UV stress is available. In the present study, UV-B treatment induced a steady up-regulation of GbHMGS1 transcription during 48 h treatment, except for a decrease at 24 h in G. biloba. Therefore, delving into the regulation mechanism of GbHMGS1 gene products may provide new insights into the processes regulated by UV-B.
MJ and SA are important phytohormones and signal molecules that regulate the response of plants against various stresses, including ozone exposure, UV radiation, and herbivore and pathogen attacks [52]. MJ and SA are elicitors that trigger the pathway of secondary metabolism in plants [53,54]. SA is also a key signal molecule in the systemic acquired resistance (SAR) of plants [55]. The exogenous application of SA has been conducted to enhance the biosynthesis of GA, GB, and bilobalide in G. biloba cell cultures [56]. Meanwhile, a set of MEP pathway genes, including GbIDS2 [57], GbDXS [23] and GbCMK2 [8], that are involved in ginkgolide biosynthesis in G. biloba were induced, they were all positively responsive to SA treatment. Furthermore, SA could up-regulate HMGS expression in G. lucidum [20], S. miltiorrhiza [31], C. acuminata [32], T. wilfordii [33], and C. nobile [34]. Similarly, the expression of the putative GbHMGS1 was also strongly induced by SA in G. biloba (Figure 7E). Thus, GbHMGS1 may be a potentially active member in SAR and participate in the synthesis of terpenoids, especially sesquiterpene for defending against local pathogen attack. Moreover, the SA-induced increase in TTL contents in the present study might be attributed to an integrated effect of a cluster of genes related to TTLs biosynthesis.
A number of reports have been published on the relationship between MJ and terpene metabolism [53,58]. Exogenously applied MJ can induce secondary defense metabolites along with the corresponding synthases in a wide range of plant species. For instance, MJ treatment enhanced the production of the mono- and diterpenoids and their respective terpene synthase activities in Picea abies [59,60]. In S. miltiorrhiza, MJ treatment enhanced the production of diterpenoid tanshinone and its terpene synthase activities [54,61]. In the present study, the accumulation of TTLs and the GbHMGS1 transcription were up-regulated by MJ (Figure 7E,F). Nevertheless, MJ is not a specific inducer of the GbHMGS1 gene because it also induces other genes. In addition, a previous study has shown that ginsenoside accumulation is far more raised by combining the elicitation of Eth and MJ than by single MJ elicitation, without causing any apparent growth inhibition [62]. This finding provides a new insight into the high accumulation of TTLs by applying diverse synergistic elicitors at the optimal concentration ratio.
Ethylene is a powerful natural regulating substance in plant metabolism and regulates many processes, including seed germination, root growth, fruit ripening, and many kinds of stress phenomena, during plant growth and development. Ethylene also modulates the expression of terpenoid biosynthesis-related genes, including GbLPS, GbMVD, GbHMGR, GbDXS, and GbGGPPS genes in G. biloba [63], by acting and interacting with other recognized plant hormones in trace amounts. The present study is the first to report that Eth increased TTL contents and up-regulated GbHMGS1 expression in G. biloba (Figure 7E,F). Addicott et al. [64] found that temperature influences the tissue response to ethylene, and applying Eth at low temperatures significantly lowers the regulatory effect on the tissue. Additionally, Bae et al. [62] found that Eth at 50 μM enhanced both root growth and ginsenoside accumulation in ginseng adventitious root cultures. Therefore, we speculate that the transcript level of GbHMGS1 may be improved by applying Eth at a proper high temperature, or the selection of the appropriate Eth concentration may further increase the TTLs content.
Considering the upregulated expression pattern of GbHMGS1 under induction of all the tested elicitors (Cold, UV-B, MJ, SA, and Eth), we suggest that GbHMGS1 is an elicitor-responsive gene, which participate in the stress resistance of ginkgo plant to adverse stresses and the signaling pathway involved in the regulation of exogenous hormones. However, the induction effects of various elicitors were diverse, which might be associated with many factors, such as the species of the elicitor, treatment concentration and time, the growth state of tissues, and different action mechanisms.
Accumulation of TTLs was enhanced by abiotic stresses and hormone treatments, may be associated with the upregulated expression of the genes in the MEP pathway, in light of the fact that the MVA pathway minimally contributes to TTL biosynthesis. For instance, Ethylene modulates the expression of terpenoid biosynthesis-related genes, including GbDXS [23], and the transcripts of GbCMK2 [8], GbDXS2 [23], and GbIDS2 [59] are also positively responsive to SA and MeJA treatments. Therefore, further work will be carried out to clarify the relationship between the TTLs biosynthesis and all gene expression involved in the MEP pathway in G. biloba.

4. Materials and Methods

4.1. Plant Materials and Multiple Stress Treatments

The plant materials harvested from 27 year-old grafts of G. biloba cultivar “Jiafoshou” grown in the Botanical Garden of Yangtze University (Jingzhou, China) were used for gene cloning and tissue expression from roots, stems, leaves, young fruits, male flowers, and female flowers. All materials were frozen immediately in liquid nitrogen and stored at −80 °C prior to total RNA extraction.
One year-old seedlings of G. biloba cultivar “Jiafoshou” grown in a greenhouse were used as exogenous stress test materials. The growth conditions were set as follows: temperature, 23 ± 0.5 °C; humidity, 60–70%; light, 100 μmol·m−2·s−1; lighting time, 12 h day/12 h night; and watering with 100 mL of liquid fertilizer (containing 0.36 g N, 0.12 g P, and 0.14 g K) weekly. Ginkgo seedlings at 4–5 leaf stages were subjected to treatments with cold, UV-B, MJ, Eth, and SA, and three biological replicates were conducted for each sampling time point under various stress conditions. For UV-B treatment, the seedlings were exposed to UV-B (1500 μJ/m2) irradiation in a closed chamber, whereas the control was placed in a dark-closed chamber, leaves of both the treated and control seedlings were sampled at 0, 8, 16, 24 and 48 h. For cold stress, the seedlings were placed in a 4 °C illuminating incubator (400 μmol−2·s−1) without changing the other growth conditions. The control was treated at 25 °C light growth chamber. For exogenous hormone treatments, 100 μM MJ, 10 mM Eth, and 10 mM SA containing 0.01% (v/v) Tween 20 were sprayed onto the leaves of ginkgo seedlings while all other aspects of the growth conditions were similar to those of the controls, which are the plants treated with 0.01% (v/v) Tween 20 as control. Leaves from the treated and control plants were collected on the 0, 2nd, 4th, 6th and 8th days after the cold and exogenous hormone treatments. All leaves of test ginkgo seedlings were frozen immediately in liquid nitrogen after sampling, followed by storage at −80 °C for the analysis of the gene transcription level and TTL contents.

4.2. Cloning of the Full-length cDNA of GbHMGS1

Total RNA was extracted from the frozen leaves of G. biloba using the CTAB method described by Cai et al. [65]. The quality and concentration of the RNA were all determined by agarose gel electrophoresis and spectrophotometer analysis. The extracted RNA was reverse transcribed into 5′-RACE-Ready cDNA and 3′-RACE-Ready cDNA using PowerScriptTM Reverse Transcriptase (Clontech, Palo Alto, CA, USA) according to the Kit User Manual. A pair of degenerate primers HMGS-FP1 and HMGS-RP2 was designed based on the conserved amino acid and nucleotide sequences of the plant HMGS genes to obtain the internal conserved fragment by using the one-step RT-PCR kit (Dalian TaKaRa, Dalian, China) under the following conditions: 94 °C for 3 min; 32 cycles of amplification at 94 °C for 30 s, 54 °C for 30 s, and 72 °C for 2 min; and extended at 72 °C for 10 min. The amplified PCR products were purified using AxyPrepTM DNA Gel Extraction Ki, ligated into the pMD18-T vector (Takara, Dalian, China), and introduced into Escherichia coli TOP10. Positive clones were selected and sent to Shanghai Sangon Biotechnology Company (Shanghai, China) for sequencing. Subsequent BLAST results confirmed that the amplified product was a partial fragment of the HMGS gene. Based on the sequence of the cloned HMGS fragment, the specific primer pairs (HMGS-5GSP1, HMGS-3GSP1) and the nested primer pairs (HMGS-5GSP2, HMGS-3GSP2) were designed to amplify the 5′ and 3′ ends of the cDNA of the GbHMGS1 gene using the SMARTTM RACE cDNA Amplification Kit (Clontech, Palo Alto, CA, USA). 5′-RACE-PCR and 3′-RACE-PCR were performed according to the manufacturer's instructions. The PCR products were purified and cloned into the pMD18-T vector for sequencing. After comparing and aligning the sequences of the 5′ RACE and 3′ RACE, a pair of primers HMGSES1 and HMGSET2 was designed based on the stitching sequence of 5′ RACE, 3′ RACE, and the internal region products for amplifying the full-length cDNA sequence of GbHMGS1 using 3′-Ready cDNA as the template. PCR was performed under the following conditions: 94 °C for 3 min, followed by 32 cycles of amplification (94 °C for 30 s, 56 °C for 30 s, and 72 °C for 2 min), and an extension for 10 min at 72 °C. After sequencing, the full-length cDNA of GbHMGS1 was analyzed for molecular characterization. The sequences of all primers in the present study were shown in in Table S1 of the Supplementary Material.

4.3. Bioinformatics and Molecular Evolution Analysis

GbHMGS1 gene sequences were translated into amino acid sequences using DNAMAN software. The open-reading frame (ORF) of the GbHMGS1 gene was predicted using Vector NTI 11.5. Protein homology searches were performed using the NCBI BLAST server (http://www.ncbi.nlm.nih.gov/BLAST/). Multiple sequence alignment was performed with the software Vector NTI11.5 program. Compute pI/Mw tool (http://web.expasy.org/protparam/) was used to predict the molecular weight and isoelectric point (pI). The comparative modeling of the 3D structure of GbHMGS1 was generated with SWISS-MODEL. The 3D structural analyses were performed with Weblab Viewerlite. A phylogenetic tree was constructed using HMGS protein sequences with neighbor-joining method using Clustal X2.0 and MEGA5.0 with 1000 bootstrap replicates.

4.4. Prokaryotic Expression and In Vitro Enzyme Activity Analysis of GbHMGS1

A pair of primers, HMGSYS2 and HMGSYT3, incorporated with restriction enzyme sites (Eco RI and Hind III, respectively) and protective bases were designed and synthesized to amplify the coding region of GbHMGS1 by RT-PCR and clone the coding region into an expression vector pET-32a(+), yielding His-tagged GbHMGS1 (pET32a-GbHMGS1). After confirmation by sequencing, the resulting recombinant plasmid was introduced into the E. coli strain BL21 (DE3) using heat shock method. A single colony of E. coli BL21 (DE3) cells harboring the expression plasmid pET32a-GbHMGS1 was inoculated in Luria–Bertani (LB) medium at 37 °C containing kanamycin (50 mg/L) and grown with 150 rpm shaking at 37 °C until the OD600 reached 0.6. For induction, the cells were induced by isopropyl-β-d-1-thiogalactopyranoside at a final concentration of 0.5 mM and further cultured at 30 °C for 3 h. The cells were lysed by sonication for 10 s at 4 °C and centrifuged at 25,000× g for 15 min. The supernates (soluble fraction) and pellets were examined by SDS-PAGE analysis, followed by Coomassie Brilliant Blue G250 staining. The recombinant GbHMGS1 protein from the induced cells was purified using Nickel-CL agarose affinity chromatography (Bangalore Genei, Bangalore, India) and used for in vitro enzyme activity determination.
The protein content was estimated by the Lowry assay [66] using bovine serum albumin as the standard. HMGS activity was determined by a modified method as described by Suvachittanont and Wititsuwannakul [18]. For the kinetics of acetoacetyl-CoA condensation by GbHMGS1, various concentrations of the substrate (20, 40, 60, 80 and 100 μM acetoacetyl-CoA) were added to 500 μL of the reaction buffer. A Michaelis–Menten (V0 = Vmax[S]/(Km + [S], R2 = 0.99) curve was drawn based on the data. The specific activity, Km, and Vmax of the recombinant GbHMGS1 protein were calculated using DPS ver. 9.5. (DPS Software Inc., Hangzhou, China).

4.5. Functional Complementation of GbHMGS1 in Yeast

The coding regions of the GbHMGS1 were amplified via PCR with the following two pairs of primers: HMGS-YS2 and HMGS-YT3. The forward primers contained the Bam HI restriction site, and the reverse primers contained the Xba I restriction site. The amplified products and the pYES2 vector (Invitrogen, Carlsbad, CA, USA) were digested with Bam HI and Xba I, and then ligated and transformed into E. coli DH5α. The constructed pYES2-GbHMGS1 plasmid was extracted from positive clones and subsequently transformed into YSC6274 with the Frozen-EZ Yeast Transformation II Kit (Zymo Research, Orange, CA, USA). The transformants were spotted on SC (-Ura) medium (6.7% yeast nitrogen base without amino acid; 2% galactose). Positive clones were further confirmed through PCR. The transformed diploid cells were induced to sporulate and form haploid cells containing pYES2-GbHMGS1. The diploid S. cerevisiae strain YSC1021 and the haploid S. cerevisiae strain YSC6274 lacking the HMGS allele were purchased from the Open Biosystems Yeast Knock-out Strain Collection (Open Biosystems, Huntsville, AL, USA). The diploid S. cerevisiae strain YSC1021 and haploid strain YSC6274 cells were separately grown on YPD (1% yeast extract, 2% bacto-peptone, and 2% glucose) + G418 (Geneticin, Thermo Fisher Scientific Inc., Waltham, MA, USA) and YPG (1% yeast extract, 2% bacto-peptone, and 2% galactose) + G418 medium to observe their growth status.

4.6. Quantitative RT-PCR Analysis of GbHMGS1 Expression Levels

The transcription levels of GbHMGS1 in different ginkgo tissues and in the treated leaf samples from the ginkgo seedlings collected at different time points after stress treatments were detected by qRT-PCR. Total RNA from each sample was extracted using the CTAB method. The PrimeScriptTM RT reagent Kit with gDNA Eraser was used with 500 ng of each total RNA to synthesize single-strand cDNA (Dalian TaKaRa). qRT-PCR was conducted using a Bio-Rad Mini OpticonTM Real-time PCR Mini Cycler (BioRad, Hercules, CA, USA) with SYBR® Premix Ex TaqTM II Kit (Dalian TaKaRa) according to the method of Xu et al. [67]. The PCR program was performed as follows: 95 °C for 5 min, followed by 40 cycles of 95 °C for 10 s and 60 °C for 30 s, with the melting curves stage (95 °C, 15 s; 60 °C, 1 min; and 95 °C, 15 s). G. biloba 18S gene (GenBank accession No. D16448) was used as the internal control for the normalization of all reactions. The primers for GbHMGS1 (GbHMGSRTS and GbHMGSRTA) and Gb18S (Gb18SF and Gb18SR) for qRT-PCR were designed using Primer Premier 5.0 (Premier Biosoft International, Palo Alto, CA, USA). The efficiency of these primers was investigated by applying primer melting curve analysis and gel electrophoresis. The two genes (GbHMGS1 and 18S) were analyzed simultaneously for each plant sample. GbHMGS1 expression was relatively quantified by calibrating its transcription level against 18S. For the post-abiotic treatment analysis of GbHMGS1 expression at different time points, the transcription levels of GbHMGS1 in the treated ginkgo seedlings were compared with those of the controls to obtain the induction level of GbHMGS1. Each sample was run in triplicate and repeated thrice. 2−ΔΔCt was calculated to analyze the relative expression level of genes [68]. qRT-PCR data were technically replicated with error bars, representing mean ± SE (n = 3).

4.7. Determination of TTL Content in Ginkgo under Induction of Abiotic Stresses

Ginkgo seedlings were collected and subjected to freeze-drying. Ginkgolide A (GA), ginkgolide B (GB), ginkgolide C (GC), and bilobalide (BB) were extracted and quantified using gas chromatography (GC-14C, Shimadzu, Kyoto, Japan) with a wide bore capillary column [69]. The content of TTLs was calculated as the sum of the contents of GA, GB, GC, and BB and was expressed as dry weight (DW) percentages. All tests were conducted in triplicate, and data were presented as mean ± SD (n = 3).

4.8. Statistical Analysis

Data were analyzed with one-way ANOVA using SPSS 11.0 for Windows (SPSS Inc., Chicago, IL, USA). Comparisons between multiple treatment groups were performed with Tukey’s honestly significant difference test. p < 0.05 was considered to be statistically significant.

5. Conclusions

This study first isolated and characterized GbHMGS1-encoding HMGS enzyme in the MVA pathway from G. biloba. Multiple alignments showed that the deduced GbHMGS1 was highly identical to other HMGSs and contained all the conserved motifs and conserved sites of HMGSs. Prokaryotic expression and functional complementation assays revealed that GbHMGS1 encoded a functional HMGS and could catalyze the conversion of acetoacetyl-CoA to HMG-CoA involved in MVA pathway. The transcript level of GbHMGS1 preferentially expressed in leaves of G. biloba. The GbHMGS1 expression under different elicitor treatments (cold, UV-B, MJ, Eth, and SA) in ginkgo seedlings further verified the crucial roles of the GbHMGS1 in stress response and sesquiterpenoid biosynthesis. Therefore, the present work on GbHMGS1 is a useful supplement to elucidate the overall molecular regulatory mechanism underlying terpenoid especially sesquiterpene biosynthesis in G. biloba.

Supplementary Materials

Supplementary data associated with this article can be found in the online version.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (31370680).

Author Contributions

X.M. performed the experiments and analyzed the data. X.M. and F.X. drafted the manuscript. Q.S. and L.W. performed the qRT-PCR analysis. J.Y. contributed the functional complementation. F.X. designed the experiments. All authors read and approved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Carrier, D.J.; van Beek, T.A.; van der Heijden, R.; Verpoorte, R. Distribution of ginkgolides and terpenoid biosynthetic activity in Ginkgo biloba. Phytochemistry 1998, 48, 89–92. [Google Scholar] [CrossRef]
  2. Mohanta, T.K.; Tamboli, Y.T.P.K. Phytochemical and medicinal importance of Ginkgo biloba L. Nat. Prod. Res. 2004, 28, 746–752. [Google Scholar] [CrossRef] [PubMed]
  3. Liao, H.J.; Zheng, Y.F.; Li, H.Y.; Peng, G.P. Two new ginkgolides from the leaves of Ginkgo biloba. Planta Med. 2011, 77, 1818–1821. [Google Scholar] [CrossRef] [PubMed]
  4. Nakanishi, K. Terpene trilactones from Gingko biloba: From ancient times to the 21st century. Bioorgan. Med. Chem. 2005, 13, 4987–5000. [Google Scholar] [CrossRef] [PubMed]
  5. Diamond, B.J.; Bailey, M.R. Ginkgo biloba: Indications, mechanisms, and safety. Psychiatr. Clin. North. Am. 2013, 36, 73–83. [Google Scholar] [CrossRef] [PubMed]
  6. Serranogarcia, N.; Pedrazachaverri, J.; Maressamano, J.J.; Orozcoibarra, M.; Cruzsalgado, A.; Jiménezanguiano, A.; Sotelo, J.; Trejosolís, C. Antiapoptotic effects of EGb 761. Evid-based Compl. Alt. 2013, 2013, 495703. [Google Scholar]
  7. Beek, T.A.V.; Montoro, P.; Xie, P.S.; Van Beek, T.A. Chemical analysis and quality control of Ginkgo biloba leaves, extracts, and phytopharmaceuticals. J. Chromatogr. A 2009, 1216, 2002–2032. [Google Scholar] [CrossRef] [PubMed]
  8. Kim, S.M.; Kim, Y.B.; Kuzuyama, T.; Kim, S. Two copies of 4-(cytidine 50-diphospho)-2-C-methyl-d-erythritol kinase (CMK) gene in Ginkgo biloba: Molecular cloning and functional characterization. Planta 2008, 228, 941–950. [Google Scholar] [CrossRef] [PubMed]
  9. Kim, J.H.; Lee, K.I.; Chang, Y.J.; Kim, S.U. Developmental pattern of Ginkgo biloba levopimaradiene synthase (GbLPS) as probed by promoter analysis in Arabidopsis thaliana. Plant Cell Rep. 2012, 31, 1119–1127. [Google Scholar] [CrossRef] [PubMed]
  10. Liao, Y.L.; Xu, F.; Huang, X.H.; Zhang, W.W.; Cheng, H.; Wang, X.H.; Cheng, S.Y.; Shen, Y.B. Characterization and transcriptional profiling of Ginkgo biloba mevalonate diphosphate decarboxylase gene (GbMVD) promoter towards light and exogenous hormone treatments. Plant Mol. Biol. Rep. 2016, 34, 566–581. [Google Scholar] [CrossRef]
  11. Tholl, D. Terpene synthases and the regulation, diversity and biological roles of terpene metabolism. Curr. Opin. Plant Biol. 2006, 9, 297–304. [Google Scholar] [CrossRef] [PubMed]
  12. Chadwick, M.; Trewin, H.; Gawthrop, F.; Wagstaff, C. Sesquiterpenoids lactones: Benefits to plants and people. Int. J. Mol. Sci. 2013, 14, 12780–12805. [Google Scholar] [CrossRef] [PubMed]
  13. Shen, G.A.; Pang, Y.Z.; Wu, W.S.; Liao, Z.H.; Zhao, L.X.; Sun, X.F.; Tang, K.X. Cloning and characterization of a root-specific expressing gene encoding 3-hydroxy-3-methylglutaryl coenzyme a reductase from Ginkgo biloba. Mol. Biol. Rep. 2006, 33, 117–127. [Google Scholar] [CrossRef] [PubMed]
  14. Schwarz, M.; Arigoni, D. Ginkgolide biosynthesis. Compr. Nat. Prod. Chem. 1999, 2, 367–399. [Google Scholar]
  15. Hemmerlin, A.; Hoeffler, J.F.; Meyer, O.; Tritsch, D.; Kagan, I.A.; Grosdemange-Billiard, C.; Rohmer, M.; Bach, T.J. Cross-talk between the cytosolic mevalonate and the plastidial alonate and the plastidial methylerythritol phosphate pathways in tobacco bright yellow-2 cells. J. Biol. Chem. 2003, 278, 26666–26676. [Google Scholar] [CrossRef] [PubMed]
  16. Patra, B.; Schluttenhofer, C.; Wu, Y.; Pattanaik, S.; Yuan, L. Transcriptional regulation of secondary metabolite biosynthesis in plants. BBA-Gene Regul. Mech. 2013, 1829, 1236–1247. [Google Scholar] [CrossRef] [PubMed]
  17. Vranová, E.; Coman, D.; Gruissem, W. Network analysis of the MVA and MEP pathways for isoprenoid synthesis. Annu. Rev. Plant Biol. 2013, 64, 665–700. [Google Scholar] [CrossRef] [PubMed]
  18. Suvachittanont, W.; Wititsuwannakul, R. 3-Hydroxy 3-methylglutaryl-CoA synthase in Hevea brasiliensis. Phytochemistry 1995, 40, 757–761. [Google Scholar] [CrossRef]
  19. Ren, A.; Ouyang, X.; Shi, L.; Jiang, A.L.; Mu, D.S.; Li, M.J.; Han, Q.; Zhao, M.W. Molecular characterization and expression analysis of GlHMGS, a gene encoding hydroxymethylglutaryl-CoA synthase from Ganoderma lucidum (Ling-zhi) in ganoderic acid biosynthesis pathway. World J. Microb. Biot. 2013, 29, 523–531. [Google Scholar] [CrossRef] [PubMed]
  20. Wang, H.; Nagegowda, D.A.; Rawat, R.; Bouvier-Nave, P.; Guo, D.; Bach, T.J.; Chye, M.L. Overexpression of Brassica juncea wildtype and mutant HMG-CoA synthase 1 in Arabidopsis up-regulates genes in sterol biosynthesis and enhances sterol production and stress tolerance. Plant Biotechnol. J. 2012, 10, 31–42. [Google Scholar] [CrossRef] [PubMed]
  21. Parveen, I.; Wang, M.; Zhao, J.P.; Chittiboyina, A.G.; Tabanca, N.; Ali, A.; Baerson, S.R.; Techen, N.; Chappell, J.; Khan, I.A.; et al. Investigating sesquiterpene biosynthesis in Ginkgo biloba: Molecular cloning and functional characterization of (E,E)-farnesol and α-bisabolene synthases. Plant Mol. Biol. 2015, 89, 451–462. [Google Scholar] [CrossRef] [PubMed]
  22. Gong, Y.F.; Liao, Z.H.; Guo, B.H.; Sun, X.F.; Tang, K.X. Molecular cloning and expression profile analysis of Ginkgo biloba DXS gene encoding 1-deoxy-d-xylulose 5-phosphate synthase, the first committed enzyme of the 2-C-methyl-d-erythritol 4-phosphate pathway. Planta Med. 2006, 72, 329–335. [Google Scholar] [CrossRef] [PubMed]
  23. Kim, S.M.; Kuzuyama, T.; Chang, Y.J.; Song, K.S.; Kim, S.U. Identification of class 2 1-deoxy-d-xylulose 5-phosphate synthase and 1-deoxy-d-xylulose 5-phosphate reductoisomerase genes from Ginkgo biloba and their transcription in embryo culture with respect to ginkgolide biosynthesis. Planta Med. 2006, 72, 234–240. [Google Scholar] [CrossRef] [PubMed]
  24. Kim, S.M.; Kuzuyama, T.; Chang, Y.J.; Kwon, H.J.; Kim, S.U. Cloning and functional characterization of 2-C-methyl-d-erythritol 4-phosphate cytidyltransferase (GbMECT) gene from Ginkgo biloba. Phytochemistry 2006, 67, 1435–1441. [Google Scholar] [CrossRef] [PubMed]
  25. Kim, S.M.; Zuyama, T.; Chang, Y.J.; Kim, S.U. Cloning and characterization of 2-C-methyl-d-erythritol 2,4-cyclodiphosphate synthase (MECS) gene from Ginkgo biloba. Plant Cell Rep. 2006, 25, 829–835. [Google Scholar] [CrossRef] [PubMed]
  26. Kim, S.M.; Kim, S.U. Characterization of 1-hydroxy-2-methyl-2-(E)-butenyl-4-diphosphate synthase (HDS) gene from Ginkgo biloba. Mol. Biol. Rep. 2010, 37, 973–979. [Google Scholar] [CrossRef]
  27. Kim, S.M.; Kuzuyama, T.; Kobayashi, A.; Sando, T.; Chang, Y.J.; Kim, S.U. 1-Hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase (IDS) is encoded by multicopy genes in gymnosperms Ginkgo biloba and Pinus taeda. Planta 2008, 227, 287–298. [Google Scholar] [CrossRef] [PubMed]
  28. Lichtenthaler, H.K.; Rohmer, M.; Schwender, J. Two independent biochemical pathways for isopentenyl diphosphate and isoprenoid biosynthesis in higher plants. Physiol. Plantarum 1997, 101, 643–652. [Google Scholar] [CrossRef]
  29. Liao, Y.L.; Xu, F.; Huang, X.H.; Zhang, W.W.; Cheng, H.; Li, L.L.; Cheng, S.Y.; Shen, Y.B. Promoter analysis and transcriptional profiling of Ginkgo biloba 3-hydroxy-3-methylglutaryl coenzyme a reductase (GbHMGR) gene in abiotic stress responses. Not. Bot. Horti Agrobot. 2015, 43, 25–34. [Google Scholar] [CrossRef]
  30. Chen, Q.W.; Yan, J.P.; Meng, X.X.; Xu, F.; Zhang, W.W.; Liao, Y.L.; Q, J.W. Molecular cloning, characterization, and functional analysis of acetyl-CoA C-acetyltransferase and mevalonate kinase genes involved in terpene trilactone biosynthesis from Ginkgo biloba. Molecules 2017, 22, 74. [Google Scholar] [CrossRef] [PubMed]
  31. Zhang, L.; Yan, X.; Wang, J.; Li, S.; Liao, P.; Kai, G. Molecular cloning and expression analysis of a new putative gene encoding 3-hydroxy-3-methylglutaryl-CoA synthase from Salvia miltiorrhiza. Acta Physiol. Plant 2011, 33, 953–961. [Google Scholar] [CrossRef]
  32. Kai, G.; Li, S.; Wang, W.; Lu, Y.; Wang, J.; Liao, P.; Cui, L. Molecular cloning and expression analysis of a gene encoding 3-hydroxy-3-methylglutaryl-CoA synthase from Camptotheca acuminata. Russ. J. Plant Physiol. 2013, 60, 131–138. [Google Scholar] [CrossRef]
  33. Liu, Y.J.; Zhao, Y.J.; Zhang, M.; Su, P.; Wang, X.J.; Zhang, X.N.; Gao, W.; Huang, L.Q. Cloning and characterisation of the gene encoding 3-hydroxy-3-methylglutaryl-CoA synthase in Tripterygium wilfordii. Molecules 2014, 19, 19696–19707. [Google Scholar] [CrossRef] [PubMed]
  34. Cheng, S.; Wang, X.; Xu, F.; Chen, Q.; Tao, T.; Lei, J.; Zhang, W.; Liao, Y.; Chang, J.; Li, X. Cloning, expression profiling and functional analysis of CnHMGS, a gene encoding 3-hydroxy-3-methylglutaryl coenzyme A synthase from Chamaemelum nobile. Molecules 2016, 21, 316. [Google Scholar] [CrossRef] [PubMed]
  35. Campobasso, N.; Patel, M.; Wilding, I.E.; Kallender, H.; Rosenberg, M.; Gwynn, M.N. Staphylococcus aureus 3-hydroxy-3-methylglutaryl-CoA synthase: Crystal structure and mechanism. J. Biol. Chem. 2004, 279, 44883–44888. [Google Scholar] [CrossRef] [PubMed]
  36. Schwede, T.; Kopp, J.; Guex, N.; Peitsch, M.C. SWISS-MODEL: An automated protein homology-modeling server. Nucleic Acids Res. 2003, 31, 3381–3385. [Google Scholar] [CrossRef] [PubMed]
  37. Theisen, M.J.; Misra, I.; Saadat, D.; Campobasso, N.; Miziorko, H.M.; Harrison, D.H. 3-hydroxy-3-methylglutaryl-CoA synthase intermediate complex observed in “real-time”. Proc. Natl. Acad. Sci. USA 2004, 101, 16442–16447. [Google Scholar] [CrossRef] [PubMed]
  38. Sirinupong, N.; Suwanmanee, P.; Doolittle, R.F.; Suvachitanont, W. Molecular cloning of a new cDNA and expression of 3-hydroxy-3-methylglutaryl-CoA synthase gene from Hevea brasiliensis. Planta 2005, 221, 502–512. [Google Scholar] [CrossRef] [PubMed]
  39. Chang, J.; Ning, Y.; Xu, F.; Cheng, S.; Li, X. Research advance of 3-hydroxy-3-methylglutaryl-coenzyme a synthase in plant isoprenoid biosynthesis. J. Anim. Plant Sci. 2015, 25, 1441–1450. [Google Scholar]
  40. Holton, T.A.; Brugliera, F.; Tanaka, Y. Cloning and expression of flavonol synthase from Petunia hybrida. Plant J. 1993, 4, 1003–1010. [Google Scholar] [CrossRef] [PubMed]
  41. Lin, G.; Lian, Y.; Ryu, J.; Sung, M.; Park, J.; Park, H.; Park, B.; Shin, J.; Lee, M.; Cheon, C. Expression and purification of His-tagged flavonol synthase of Camellia sinensis from Escherichia coli. Protein Expres. Purif. 2007, 55, 287–292. [Google Scholar] [CrossRef] [PubMed]
  42. Owens, D.K.; Alerding, A.B.; Crosby, K.C.; Bandara, A.B.; Westwood, J.H.; Winkel, B.S.J. Functional analysis of a predicted flavonol synthase gene family in Arabidopsis. Plant Physiol. 2008, 147, 1046–1061. [Google Scholar] [CrossRef] [PubMed]
  43. Misra, I.; Wang, C.Z.; Miziorko, H.M. The influence of conserved aromatic residues in 3-hydroxy-3-methylglutaryl-CoA synthase. J. Biol. Chem. 2003, 278, 26443–26449. [Google Scholar] [CrossRef] [PubMed]
  44. Cartayrade, A.; Neau, E.; Sohier, C.; Balz, J.P.; Carde, J.P.; Walter, J. Ginkgolide and bilobalide biosynthesis in Ginkgo biloba. 1. Sites of synthesis, translocation and accumulation of ginkgolides and bilobalide. Plant Physiol. Biochem. 1997, 35, 859–868. [Google Scholar]
  45. Lu, X.; Yang, H.; Liu, X.; Shen, O.; Wang, N.; Qi, L.W.; Li, P. Combining metabolic profiling and gene expression analysis to reveal the biosynthesis site and transport of ginkgolides in Ginkgo biloba L. Front. Plant Sci. 2017, 8, 872. [Google Scholar] [CrossRef] [PubMed]
  46. Alex, D.; Bach, T.J.; Chye, M.L. Expression of Brassica juncea 3-hydroxy-3-methylglutaryl CoA synthase is developmentally regulated and stress-responsive. Plant J. 2000, 22, 415–426. [Google Scholar] [CrossRef] [PubMed]
  47. Xu, F.; Cai, R.; Cheng, S.; Du, H.; Wang, Y.; Cheng, S. Molecular cloning, characterization and expression of phenylalanine ammonialyase gene from Ginkgo biloba. Afr. J. Biotechnol. 2008, 7, 721–729. [Google Scholar]
  48. Xu, F.; Cheng, H.; Cai, R.; Li, L.L.; Chang, J.; Zhu, J.; Zhang, F.X.; Chen, L.J.; Wang, Y.; Cheng, S.H.; et al. Molecular cloning and function analysis of an anthocyanidin synthase gene from Ginkgo biloba, and its expression in abiotic stress responses. Mol. Cells 2008, 26, 536–547. [Google Scholar] [PubMed]
  49. Xu, F.; Li, L.; Zhang, W.; Cheng, H.; Sun, N.; Cheng, S.; Wang, Y. Isolation, characterization, and function analysis of a flavonol synthase gene from Ginkgo biloba. Mol. Biol. Rep. 2012, 39, 2285–2296. [Google Scholar] [CrossRef] [PubMed]
  50. Cheng, H.; Li, L.; Xu, F.; Cheng, S.; Cao, F.; Wang, Y.; Yuan, H.; Jiang, D.; Wu, C. Expression patterns of a cinnamyl alcohol dehydrogenase gene involved in lignin biosynthesis and environmental stress in Ginkgo biloba. Mol. Biol. Rep. 2013, 40, 707–721. [Google Scholar] [CrossRef] [PubMed]
  51. Cheng, H.; Li, L.; Cheng, S.; Cao, F.; Wang, Y.; Yuan, H. Molecular cloning and function assay of a chalcone isomerase gene (GbCHI) from Ginkgo biloba. Plant Cell Rep. 2011, 30, 49–62. [Google Scholar] [CrossRef] [PubMed]
  52. Robert-Seilaniantz, A.; Grant, M.; Jones, J.D.G. Hormone crosstalk in plant disease and defense: More than just jasmonate-salicylate antagonism. Annu. Rev. Phytopathol. 2011, 49, 317–343. [Google Scholar] [CrossRef] [PubMed]
  53. Pauwels, L.; Morreel, K.; Witte, E.D.; Lammertyn, F.; Montagu, M.V.; Boerjan, W.; Inzée, D.; Goossens, A. Mapping methyl jasmonate-mediated transcriptional reprogramming of metabolism and cell cycle progression in cultured Arabidopsis cells. Proc. Natl. Acad. Sci. USA 2008, 105, 1380–1385. [Google Scholar] [CrossRef] [PubMed]
  54. Dong, J.; Wan, G.; Liang, Z. Accumulation of salicylic acid-induced phenolic compounds and raised activities of secondary metabolic and antioxidative enzymes in Salvia miltiorrhiza cell culture. J. Biotechnol. 2010, 148, 99–104. [Google Scholar] [CrossRef] [PubMed]
  55. Gaffney, T.; Friedrich, L.; Vernooij, B.; Negrotto, D.; Nye, D.; Uknes, S.; Ward, E.; Kessmann, H.; Ryals, J. Requirement of salicylic acid for the induction of systemic acquired resistance. Science 1993, 261, 754–756. [Google Scholar] [CrossRef] [PubMed]
  56. Kang, S.M.; Min, J.Y.; Kim, Y.D. Effects of methyl jasmonate and salicylic acid on the production of bilobalide and ginkgolides in cell cultures of Ginkgo biloba. In Vitro Cell. Dev. Biol. Plant 2006, 42, 44–49. [Google Scholar] [CrossRef]
  57. Kang, M.K.; Nargis, S.; Kim, S.M.; Kim, S.U. Distinct expression patterns of two Ginkgo biloba, 1-hydroxy-2-methyl-2-(E)-butenyl-4-diphosphate reductase/isopentenyl diphospahte synthase (HDR/IDS) promoters in Arabidopsis model. Plant Physiol. Biochem. 2013, 62, 47–53. [Google Scholar] [CrossRef] [PubMed]
  58. Wang, Q.J.; Zheng, L.P.; Zhao, P.F.; Zhao, Y.L.; Wang, J.W. Cloning and characterization of an elicitor-responsive gene encoding 3-hydroxy-3-methylglutaryl coenzyme A reductase involved in 20-hydroxyecdysone production in cell cultures of Cyanotis arachnoidea. Plant Physiol. Biochem. 2014, 84, 1–9. [Google Scholar] [CrossRef] [PubMed]
  59. Martin, D.; Tholl, D.; Gershenzon, J.; Bohlmann, J. Methyl jasmonate induces traumatic resin ducts, terpenoid resin biosynthesis, and terpenoid accumulation in developing xylem of Norway spruce stems. Plant Physiol. 2002, 129, 1003–1018. [Google Scholar] [CrossRef] [PubMed]
  60. Keeling, C.I.; Bohlmann, J. Diterpene resin acids in conifers. Phytochemistry 2006, 67, 2415–2423. [Google Scholar] [CrossRef] [PubMed]
  61. Dai, Z.; Cui, G.; Zhou, S.F.; Zhang, X.; Huang, L. Cloning and characterization of a novel 3-hydroxy-3-methylglutaryl coenzyme A reductase gene from Salvia miltiorrhiza involved in diterpenoid tanshinone accumulation. J. Plant Physiol. 2011, 168, 148–157. [Google Scholar] [CrossRef] [PubMed]
  62. Bae, K.H.; Choi, Y.; Shin, C.G.; Kim, Y.Y.; Kim, Y.S. Enhanced ginsenoside productivity by combination of ethephon and methyl jasmoante in ginseng (Panax ginseng C.A. Meyer) adventitious root cultures. Biotechnol. Lett. 2006, 28, 1163. [Google Scholar] [CrossRef] [PubMed]
  63. Huang, X.H. Molecular Cloning and Function Analysis of Promoters of the Terpene Lactones Synthase Key Gene from Ginkgo biloba, and Their Expressions in Abiotic Stress Responses. M.D. Thesis, Yangtze University, Jingzhou, China, 2013. [Google Scholar]
  64. Woolf, A.B.; Clemens, J.; Plummer, J.A. Leaf maturity and temperature affect the selective removal of floral buds from Camellia with ethephon. Laryngoscope 1995, 123, 2176–2179. [Google Scholar]
  65. Cai, R.; Xu, F.; Chen, L.; Cheng, S. Modification of total RNA isolation method from different Ginkgo biloba organs. Biotechnology 2007, 17, 38–41. [Google Scholar]
  66. Lowry, O.H.; Rosebrough, N.J.; Farr, A.L.; Randall, R.J. Protein measurement with the Folin phenol reagent. J. Biol. Chem. 1951, 193, 265. [Google Scholar] [PubMed]
  67. Xu, F.; Ning, Y.J.; Zhang, W.W.; Liao, Y.L.; Li, L.L.; Chengm, H.; Cheng, S.Y. An R2R3-MYB transcription factor as a negative regulator of the flavonoid biosynthesis pathway in Ginkgo biloba. Funct. Integr. Genom. 2014, 14, 177–189. [Google Scholar] [CrossRef] [PubMed]
  68. Schmittgen, T.D.; Livak, K.J. Analyzing real-time PCR data by the comparative CT method. Nat. Protoc. 2008, 3, 1101–1108. [Google Scholar] [CrossRef] [PubMed]
  69. Liao, Y.; Xu, F.; Zhu, J.; Wang, Y.; Cheng, S. Separation and determination of terpene trilactiones by gas chromatography with wide bore capillary column. Acta Agric. Boreali-occidentalis Sin. 2008, 17, 146–149. [Google Scholar]
Sample Availability: Samples of Ginkgo biloba and cDNAs of GbHMGS1 are available from the authors.
Molecules 22 01706 g001 550
Figure 1. MVA and MEP pathways for terpenoid biosynthesis in Ginkgo biloba.
Figure 1. MVA and MEP pathways for terpenoid biosynthesis in Ginkgo biloba.
Molecules 22 01706 g001
Molecules 22 01706 g002 550
Figure 2. The multiple alignments of GbHMGS1 amino acid sequence with other HMGS proteins. The conservative motif, “NxD/NE/VEGI/VDx(2)NACF/YxG” is boxed. The active sites are indicated with asterisk. The species, protein names and GenBank accession numbers are as follows: G. biloba (GbHMGS1, MF281039); T. media (TmHMGS, AAT73206.1); P. sylvestris (PsHMGS, CAA65250.1); H. brasiliensis (HbHMGS, AAK73854.1); C. acuminate (CaHMGS, ACD87446.1); A. thaliana (AtHMGS, CAA58763.1); T. cacao (TcHMGS, XP_007040101.2); N. tazetta (NtHMGS, AHF81872.1). Dark blue: identity = 100%; red: 75% ≤ identity < 100%; light blue: 50% ≤ identity < 75%.
Figure 2. The multiple alignments of GbHMGS1 amino acid sequence with other HMGS proteins. The conservative motif, “NxD/NE/VEGI/VDx(2)NACF/YxG” is boxed. The active sites are indicated with asterisk. The species, protein names and GenBank accession numbers are as follows: G. biloba (GbHMGS1, MF281039); T. media (TmHMGS, AAT73206.1); P. sylvestris (PsHMGS, CAA65250.1); H. brasiliensis (HbHMGS, AAK73854.1); C. acuminate (CaHMGS, ACD87446.1); A. thaliana (AtHMGS, CAA58763.1); T. cacao (TcHMGS, XP_007040101.2); N. tazetta (NtHMGS, AHF81872.1). Dark blue: identity = 100%; red: 75% ≤ identity < 100%; light blue: 50% ≤ identity < 75%.
Molecules 22 01706 g002
Molecules 22 01706 g003 550
Figure 3. The three-dimensional model of GbHMGS1. The helices are indicated by helices in red and conservative area are indicated by helices in yellow. The sheets are indicated by patches in blue and yellow. Turns and loops are indicated by lines.
Figure 3. The three-dimensional model of GbHMGS1. The helices are indicated by helices in red and conservative area are indicated by helices in yellow. The sheets are indicated by patches in blue and yellow. Turns and loops are indicated by lines.
Molecules 22 01706 g003
Molecules 22 01706 g004 550
Figure 4. Phylogenetic tree of HMGSs from different species using the neighbor-joining method. The reliability of the tree is measured by bootstrap analysis with 1000 replications. The bars represent evolutionary distance.
Figure 4. Phylogenetic tree of HMGSs from different species using the neighbor-joining method. The reliability of the tree is measured by bootstrap analysis with 1000 replications. The bars represent evolutionary distance.
Molecules 22 01706 g004
Molecules 22 01706 g005 550
Figure 5. Functional complementation for the growth of the yeast strain YSC6274 complemented with GbHMGS1. (A) Diploid YSC1021 strain on YPG + G418 medium grew within 2 days; (B) Diploid YSC1021 strain on YPD + G418 medium grew within 2 days; (C) Haploid YSC6274 strain containing pYES2-GbHMGS1 on YPG + G418 medium grew within 2 days; (D) Haploid YSC6274 strain containing pYES2-GbHMGS1 on YPD + G418 medium failed to grow.
Figure 5. Functional complementation for the growth of the yeast strain YSC6274 complemented with GbHMGS1. (A) Diploid YSC1021 strain on YPG + G418 medium grew within 2 days; (B) Diploid YSC1021 strain on YPD + G418 medium grew within 2 days; (C) Haploid YSC6274 strain containing pYES2-GbHMGS1 on YPG + G418 medium grew within 2 days; (D) Haploid YSC6274 strain containing pYES2-GbHMGS1 on YPD + G418 medium failed to grow.
Molecules 22 01706 g005
Molecules 22 01706 g006 550
Figure 6. Spatial expression analysis of GbHMGS1 in various tissues of G. biloba. R: root; S: stem; L: leaf; Ff: female flower; Mf: male flower; Yf: young fruit. The gene expression level of GbHMGS1 in root was set to 1. Data are mean ± SD from triplicate experiments (n = 3). Means with different letters are significantly different at p < 0.05.
Figure 6. Spatial expression analysis of GbHMGS1 in various tissues of G. biloba. R: root; S: stem; L: leaf; Ff: female flower; Mf: male flower; Yf: young fruit. The gene expression level of GbHMGS1 in root was set to 1. Data are mean ± SD from triplicate experiments (n = 3). Means with different letters are significantly different at p < 0.05.
Molecules 22 01706 g006
Molecules 22 01706 g007 550
Figure 7. GbHMGS1 transcript level and total terpene trilactone (TTL) content changes in G. biloba under the induction of Cold, UV-B, SA, MJ, Eth. (A) GbHMGS1 transcript level changes by Cold; (B) TTL content changes by Cold; (C) GbHMGS1 transcript level changes by UV-B; (D) TTL content changes by UV-B; (E) GbHMGS1 transcript level changes by SA, MJ and Eth; (F) TTL content changes by SA, MJ and Eth. The expression levels were normalized to Gb18S gene. Data were analyzed as expression ratios relative to the level of control (CK). Data are mean ± SD from triplicate experiments (n = 3). Means with different letters from each time of post-treatment are significantly different at p < 0.05.
Figure 7. GbHMGS1 transcript level and total terpene trilactone (TTL) content changes in G. biloba under the induction of Cold, UV-B, SA, MJ, Eth. (A) GbHMGS1 transcript level changes by Cold; (B) TTL content changes by Cold; (C) GbHMGS1 transcript level changes by UV-B; (D) TTL content changes by UV-B; (E) GbHMGS1 transcript level changes by SA, MJ and Eth; (F) TTL content changes by SA, MJ and Eth. The expression levels were normalized to Gb18S gene. Data were analyzed as expression ratios relative to the level of control (CK). Data are mean ± SD from triplicate experiments (n = 3). Means with different letters from each time of post-treatment are significantly different at p < 0.05.
Molecules 22 01706 g007
Table
Table 1. The enzyme activity determination of recombinant GbHMGS1.
Table 1. The enzyme activity determination of recombinant GbHMGS1.
Enzymekcat (min−1)Km (μM)kcat/Km (μM−1·min−1)Specific Activity (μmol·min−1·mg−1)Vmax (nM·min−1)
No HMGS in pET32NDNDND3.4 ± 0.5ND
GbHMGS1195.4 ± 13.1689 ± 32.60.28 ± 0.0115.5 ± 4.837.1 ± 5.4
HvHMGSND530 ± 50ND41.3 ± 2.9130 ± 30
ND, not determined.

Share and Cite

MDPI and ACS Style

Meng, X.; Song, Q.; Ye, J.; Wang, L.; Xu, F. Characterization, Function, and Transcriptional Profiling Analysis of 3-Hydroxy-3-methylglutaryl-CoA Synthase Gene (GbHMGS1) towards Stresses and Exogenous Hormone Treatments in Ginkgo biloba. Molecules 2017, 22, 1706. https://doi.org/10.3390/molecules22101706

AMA Style

Meng X, Song Q, Ye J, Wang L, Xu F. Characterization, Function, and Transcriptional Profiling Analysis of 3-Hydroxy-3-methylglutaryl-CoA Synthase Gene (GbHMGS1) towards Stresses and Exogenous Hormone Treatments in Ginkgo biloba. Molecules. 2017; 22(10):1706. https://doi.org/10.3390/molecules22101706

Chicago/Turabian Style

Meng, Xiangxiang, Qiling Song, Jiabao Ye, Lanlan Wang, and Feng Xu. 2017. "Characterization, Function, and Transcriptional Profiling Analysis of 3-Hydroxy-3-methylglutaryl-CoA Synthase Gene (GbHMGS1) towards Stresses and Exogenous Hormone Treatments in Ginkgo biloba" Molecules 22, no. 10: 1706. https://doi.org/10.3390/molecules22101706

APA Style

Meng, X., Song, Q., Ye, J., Wang, L., & Xu, F. (2017). Characterization, Function, and Transcriptional Profiling Analysis of 3-Hydroxy-3-methylglutaryl-CoA Synthase Gene (GbHMGS1) towards Stresses and Exogenous Hormone Treatments in Ginkgo biloba. Molecules, 22(10), 1706. https://doi.org/10.3390/molecules22101706

Article Metrics

Back to TopTop