1. Introduction
Flower color, a vital trait of ornamental plants, is mainly determined by flavonoid [
1,
2]. Flavonoids are a large group of plant natural pigments which comprise of chalcones, aurones, flavones, isoflavones, flavandiols, flavonols, proanthocyanins, anthocyanins and so on [
3]. In addition to tissue pigmentation, flavonoids also fulfill significant roles in multiple physiological processes, such as UV protection, auxin transport, phytopathogens and herbivores defense, signaling between plants and microbes as well as pollen development [
4,
5]. Importantly, flavonoids have a wide variety of health benefits to human, including delay the aging of immune organs and nervous system, eyesight improvement, prevention of the cancer, Alzheimer’s diseasea and cardiovascular disease [
6,
7]. Extensive studies about flavonoid biosynthetic pathway have been finished through using
Petunia hybrida (petunia) [
8],
Antirrhinum majus (snapdragon) [
9] and
Arabidopsis thaliana (Arabidopsis) [
10] as models, and a great deal enzymes that taken part in flavonoid biosynthesis have been characterized [
11].
Chalcone synthase (CHS), a polyketide synthase, is crucial in biosynthetic pathway of flavonoid and serves as the gatekeeper to regulate its biosynthesis [
12]. CHS catalyzes a three-step condensative reaction from trimolecular malonyl-CoA and monomolecular
p-coumaroyl-CoA to produce naringenin chalcone [
13]. The first step, coumaroyl moiety is loaded to the site of active cysteine (Cys164). Next, condensation reactions occur based on the decarboxylation of malonyl-CoA to produce nucleophile for chain elongation. Finally, the reaction generate thioester-linked tetraketide which cyclizes into a hydroxylated aromatic ring to yield chalcone [
14]. Due to the important roles on starting flavonoid biosynthesis, CHSs have been intensively researched in many higher plants, including petunia [
8],
Arabidopsis [
10], snapdragon [
9],
Oncidium orchid [
15], apple [
16], grape [
17],
Gerbera hybrida [
18],
Dianthus chinensis [
19] and so on. In
Antirrhinum majus, the first
CHS mutant called white nivea was described, which was generated by deleting its single
CHS gene [
20]. For petunia, there are more than 8 copies of
CHS genes in its genome, but only
CHS-A and
CHS-J transcript in floral tissues and conspicuously down-regulated in the white parts of ‘Red Star’ [
21,
22]. In addition, the spatial suppression of
CHSA also caused the natural bicolor floral phenotype as well as inability to generate functional pollen tubes [
23]. In contrast to petunia, the maize
CHS mutants can still initiate pollen tube growth, but is incapable of sustaining its growth exceed 12 h, ultimately leading to male sterile [
24]. The
Arabidopsis tt4 mutant which lacks the brown tannins in seed and anthocyanin in the cotyledons and hypocotyls are demonstrated causing by a mutation in the
CHS gene, and other phenotypic effects such as change in pollen fluorescence and root morphology were also reported later [
25,
26]. In
Gerbera hybrida, three like genes of
CHS are expressed during its flower development, whereas only
CHS1 transcript corresponds with the synthesis of anthocyanins and flavanols [
27]. As for Asiatic hybrid lily, also three
CHS genes (
CHSA-CHSC) are active in colored tepals, but transcript patterns of them are diverse [
28,
29]. Therefore, all the mentioned results indicate that CHS plays indispensable roles during plant development and is crucial for petal color formation in some plants.
Rhododendron delavayi (
R. delavayi) belonging to Ericaceae family is one of the significant ornamental plant species. Considering its vivid flowers and high horticultural values,
R. delavayi has been becoming increasingly popular in the world [
30]. However, the first rate-limiting enzyme, chalcone synthase, which is vital for petal coloration has not been isolated and analyzed from
R. delavayi. Here a
CHS gene (named as
RdCHS1) was isolated from petals of
R. delavayi and its functional roles on flavonoid biosynthesis had been comprehensively demonstrated. Expression analysis of
RdCHS1 was conducted showing that it might participate in the biosynthesis of flavonoid (not only anthocyanin) in all detected tissues. Subsequently, catalytic property of RdCHS1 was confirmed and its function in planta was verified by transferring into
Arabidopsis tt4 mutant. At the same time,
RdCHS1 was also ectopicly expressed in tobacco, the data displayed that it could strengthen the pink color of the corolla from pale pink to dark pink. In this paper, we comprehensively identify the function of chalcone synthase in
R. delavayi, and prove
RdCHS1is a
CHS gene working on flavonoid biosynthesis in
R. delavayi.
3. Discussion
The CHS enzyme family is crucial for plant growth and development. They are ubiquitous in different plant species and encoded by multiple genes. For example, 3, 4, 12 and 14
CHS genes have been cloned and identified from
Lilium spp. [
28,
29],
Dahlia variabilis [
37],
Zea mays [
38] and
Petunia hybrida [
39], several of them are the true CHSs while others involve in various metabolic pathway [
40].
Arabidopsis has four
CHS genes, of which one is the true CHS and takes part in flavonoid biosynthesis [
41]. Of RdCHSs, Protein sequence alignment and phylogenetic analysis showed that RdCHS1 was closely related to the bona fide CHSs indicating its ability to produce naringenin chalcone (
Figure 1, 2).
CHSs of plant share highly similarity in their amino acid sequence. RdCHS1 amino acid are 82.6-91.4% identical with
Arabidopsis CHS and
Vitis vinifera CHS (
Figure 1). An increasing number of PKSs (except for CHS), such as stilbene synthase (STS), acridone synthase (ACS) and 2-pyrone synthase (2PS) have been demonstrated possess the identical catalytic mechanism with CHS but are different in the intramolecular cyclization patterns and predilection for starter substrate [
42,
43]. The above enzymes are also quite similar at amino acid level with CHS. Thus, plentiful CHS sequences in public databases, identified only through their sequence similarity, may encode other related enzymes in fact [
44]. Alternatively, the kinds of amino acid in protein sequence could help to define the CHS. The Phe residues (Phe215 and Phe265) are two gatekeepers in CHSs that facilitate substrate loading and appropriate folding during cyclization process [
14,
42]. Phe265, a critical residue for the substrate selectivity, is conserved in the sequence of CHSs but changes in other plant PKSs [
43,
45]. In OsCHS9, Phe265 is substituted with Gly which displays undetectable CHS activity, instead encodes GUS catalyzing the production of bisdemethoxycurcumin [
46]. RdCHS1 contains both Phe residues in amino acid sequence implying it was functional CHSs in
R. delavayi and exhibits CHS activity. While in RdCHS2 and RdCHS3, Phe265 is replaced by Val indicating their function as other plant PKSs (
Figure 1), thus RdCHS1 was selected for further analysis.
Gene expression patterns are correlated with their functions, differential analysis of gene expression can provide key information for the study of gene features, regulations as well as origin [
47,
48]. The expression patterns of
RdCHS1 are spatially regulated. Its highest mRNA level was detected in leaves followed by flower tissues, and least in roots. Meanwhile,
RdCHS1 expression in flowers is also developmentally regulated, but not associated with accumulation of total anthocyanins (
Figure 3). These are different from the results in
Petunia hybrida and
Gerbera hybrida, in which CHS expression is excessive in flower and coupled with anthocyanin pigmentation [
21,
27]. Together with the high expression in leaves,
RdCHS1 expression profiling in this research suggest that it might not be the dominant CHS enzyme in
R. delavayi petals.
In vitro enzymatic assays showed RdCHS1 was an authentic CHS similar to CHS enzymes in
Gerbera hybrida (GCHS1 and GCHS3) that could convert malonyl-CoA and
p-coumaroyl-CoA molecules into naringenin chalcone (
Figure 4) [
18]. Although RdCHS1 performed typical CHS function, it does not seem to be a major functional CHS in petal pigmentation according to its highest expression in leaf and lower catalytic efficiency towards malonyl-CoA and
p-coumaroyl-CoA. Recent studies have reported that CHS gene in
Physcomitrella patens can accept dihydro-p-coumaroyl-CoA and cinnamoyl-CoA to produce relevant chalcones [
49]. Similarly,
CHS gene from
Scutellaria baicalensis also has the ability to convert isovaleryl-CoA, phenylacetyl-CoA, isobutyryl-CoA and benzoyl-CoA into a variety of products containing the aromatic polyketide which is unnatural [
50]. Therefore, CHS is generally a promiscuous enzyme in the light of substrate specificity which suggests its functional diversification during the process of evolution. In view of above discoveries, further researches are needed to explore the catalytic properties of RdCHS1 so as to determine its function in plants and also lay the foundation for functional divergence study of
CHS gene family in
R. delavayi.
A few studies have observed the influence of reducing CHS activity in various plants. Such as in
Arabidopsis,
tt4 was the CHS gene mutant, which displayed a deficiency in the synthesis of anthocyanin, and also exhibited an absence on the pigment of seed coat [
51]. Thus,
tt4 mutant is a suitable model for verifying whether
RdCHS1 is taken part in anthocyanin and proanthocyanidin biosynthesis. As expected,
RdCHS1 completely restored the purple coloration of
tt4 in hypocotyls and cotyledons and the pigment in seed coats, which confirmed the function of RdCHS1 as CHS
in vivo (
Table S2,
Figure 5). These results are similar to the case of maize, its
C2 gene encoding CHS was also overexpressed in
tt4 mutant and exhibited alike pigmentation phenotype as well as accumulation patterns of flavonoid [
52]. Meanwhile, complementation of
Arabidopsis flavonoid mutants were also conducted by other enzymes involved in flavonoid biosynthesis such as chalcone isomerase from
Ophiorrhiza japonica and dihydroflavonol 4-reductase from
Dryopteris erythrosora, all these findings demonstrated that the function of enzymes participated in flavonoid biosynthesis are exchangeable among different plant species [
53,
54].
Comparing to wild-type tobacco, overexpression of
RdCHS1 gene resulted in dark-pink flowers and increased cyanidin-type anthocyanins in transgenic lines (
Table S4,
Figure 7), and this similar phenomenon had also been looked at in
Solanum lycopersicum [
55]. Interestingly, the enhanced biosynthesis of anthocyanins in
RdCHS1 transgenic tobacco is due to the increased expression of endogenous structure genes (
NtCHS,
NtCHI,
NtF3H,
NtF3’H,
NtANS and
NtUFGT) and regulatory factors (
NtAN1a,
NtAN1b,
NtAN2) (
Figure 8). Up-regulatuon of these endogenous gene expressions may be due to the positive feedback regulation of flavonoid biosynthetic genes through pathway intermediates, or that RdCHS1 protein may directly interact with transcriptional regulation proteins to increase anthocyanin accumulation [
56].
NtAN1a/
NtAN1b (two bHLH transcription factors) and
NtAN2 (a R2R3-Myb transcription factor) had been demonstrated regulating anthocyanin synthesis in tobacco, thus their massive expressions promote the total metabolic flux in transgenic tobacco [
57,
58].Unexpectedly, the transcript level of
NtF3’5’H remained constant, and this may be because the host tobacco is incapable to produce delphinidin-type anthocyanins on account of dihydromyricetin deficiency [
59]. Overall, the above results prove that
RdCHS1 is an anthocyanin-related gene and involved in flower color formation in
R. delavayi.
Figure 1.
Amino acid sequences alignment of CHS protein in R. delavayi with proteins from other species. The yellow box represents three conserved catalytic residues in CHS. The blue frame amino acid determines the specificity of CHS substrate. The green rectangular box indicates the highly conserved domains of CHS. The black triangles represent important residues in binding to coumarinyl coenzyme A and residues specific to the cyclic reaction of CHS.
Figure 1.
Amino acid sequences alignment of CHS protein in R. delavayi with proteins from other species. The yellow box represents three conserved catalytic residues in CHS. The blue frame amino acid determines the specificity of CHS substrate. The green rectangular box indicates the highly conserved domains of CHS. The black triangles represent important residues in binding to coumarinyl coenzyme A and residues specific to the cyclic reaction of CHS.
Figure 2.
Phylogenetic analyses of RdCHSs. Plant species and GenBank accession numbers are as follows: FhCHS (Freesia hybrida, AEO45114.1), IgCHS (Iris germanica, BAE53636.1), VvCHS (Vitis vinifera, BAA31259.1), AtCHS (Arabidopsis thaliana, AAA32771.1), ZmCHS (Zea mays, CAA42763.1), OsCHS (Oryza sativa, BAA19186.2), GbCHS (Ginkgo biloba, AAT68477.1), MsCHS2 (Medicago sativa, P30074.1), PnCHS (Psilotum rudum, BAA87922), EaCHS (Equisetum arvense, Q9MBB1.1), AmQNS (Aegle marmelos, AGE44110), RgACS (Ruta graveolens, CAC14058.1), RpBAS (Rheuam palmatum, AAK82824.1), VvSTS (Vitis vinifera, ABV82966.1), PsSTS (Pinus sylvestris, CAA43165), Gh2PS (Gerbera hybrida, P48391.2), RpALS (Rheum palmatum, AAS87170), CsOLS (Cannabis sativa, B1Q2B6), HIVPS (Huamulus lupulus, ACD69659.1), HaBPS (Hypericum androsaemum, AAL79808.1), MdBIS1 (Malus domestica, NP001315967), BfBBS (Bromheadia finiaysoniana, CAA10514.1), WtPKS1 (Wachendorfia thyrsiflora, AAW50921), RdORS (Rhododendron dauricum, BAV83003), AtPKSA (Arabidopsis thaliana, O23674), and AtPKSB (Arabidopsis thaliana,Q8LDM2). QNS, OLS, ALS, BIS, and ORS stand for quinolone synthase, olivetol synthase, aloesone synthase, 3, 5-dihydroxybiphenol synthase, and orcinol synthase, respectively.
Figure 2.
Phylogenetic analyses of RdCHSs. Plant species and GenBank accession numbers are as follows: FhCHS (Freesia hybrida, AEO45114.1), IgCHS (Iris germanica, BAE53636.1), VvCHS (Vitis vinifera, BAA31259.1), AtCHS (Arabidopsis thaliana, AAA32771.1), ZmCHS (Zea mays, CAA42763.1), OsCHS (Oryza sativa, BAA19186.2), GbCHS (Ginkgo biloba, AAT68477.1), MsCHS2 (Medicago sativa, P30074.1), PnCHS (Psilotum rudum, BAA87922), EaCHS (Equisetum arvense, Q9MBB1.1), AmQNS (Aegle marmelos, AGE44110), RgACS (Ruta graveolens, CAC14058.1), RpBAS (Rheuam palmatum, AAK82824.1), VvSTS (Vitis vinifera, ABV82966.1), PsSTS (Pinus sylvestris, CAA43165), Gh2PS (Gerbera hybrida, P48391.2), RpALS (Rheum palmatum, AAS87170), CsOLS (Cannabis sativa, B1Q2B6), HIVPS (Huamulus lupulus, ACD69659.1), HaBPS (Hypericum androsaemum, AAL79808.1), MdBIS1 (Malus domestica, NP001315967), BfBBS (Bromheadia finiaysoniana, CAA10514.1), WtPKS1 (Wachendorfia thyrsiflora, AAW50921), RdORS (Rhododendron dauricum, BAV83003), AtPKSA (Arabidopsis thaliana, O23674), and AtPKSB (Arabidopsis thaliana,Q8LDM2). QNS, OLS, ALS, BIS, and ORS stand for quinolone synthase, olivetol synthase, aloesone synthase, 3, 5-dihydroxybiphenol synthase, and orcinol synthase, respectively.
Figure 3.
Expression profiles of RdCHS1 in R. delavayi. (A) Relative expression levels of RdCHS1 gene in different tissues; Pe, petals; To, toruses; Sc, scapes; Pi, pistils; St, stamens; Ro, roots; Le, leaves. (B) Relative expression levels of RdCHS1 at five flower developmental stages; S1, flower buds about 1 cm; S2, flower buds about 1.5 cm; S3, flower buds about 2 cm; S4, freshly opened flowers; S5, blooming flowers. Results represent means ± SE from three biological replicates.
Figure 3.
Expression profiles of RdCHS1 in R. delavayi. (A) Relative expression levels of RdCHS1 gene in different tissues; Pe, petals; To, toruses; Sc, scapes; Pi, pistils; St, stamens; Ro, roots; Le, leaves. (B) Relative expression levels of RdCHS1 at five flower developmental stages; S1, flower buds about 1 cm; S2, flower buds about 1.5 cm; S3, flower buds about 2 cm; S4, freshly opened flowers; S5, blooming flowers. Results represent means ± SE from three biological replicates.
Figure 4.
Biochemical assays of recombinant RdCHS1. (A) Expression of RdCHS1 in E. coli. (1) Maker (2) Total soluble protein from E. coli expressing pET-32a (+) vector (3) Total soluble protein from E. coli expressing RdCHS1 prior to induction by IPTG (4) Purified RdCHS1. (B) Standard of naringin chalcone. (C) The control (empty pET-32a vector). (D) HPLC profiles of the reaction products of RdCHS1.
Figure 4.
Biochemical assays of recombinant RdCHS1. (A) Expression of RdCHS1 in E. coli. (1) Maker (2) Total soluble protein from E. coli expressing pET-32a (+) vector (3) Total soluble protein from E. coli expressing RdCHS1 prior to induction by IPTG (4) Purified RdCHS1. (B) Standard of naringin chalcone. (C) The control (empty pET-32a vector). (D) HPLC profiles of the reaction products of RdCHS1.
Figure 5.
Complementation of RdCHS1 function in Arabidopsis tt4 mutant. (A) Phenotype of the wild type, mutant (tt4) and T2 transgenic lines. (B) Expressional analysis of RdCHS1 in wild-type, mutant and transgenic lines. (C) Total contents of anthocyanins in Arabidopsis seedlings. Results correspond to means from three biological replicates. Asterisks indicate significant differences between means of wild-type and transgenic plants calculated by Tukey HSD test (***P<0.001).
Figure 5.
Complementation of RdCHS1 function in Arabidopsis tt4 mutant. (A) Phenotype of the wild type, mutant (tt4) and T2 transgenic lines. (B) Expressional analysis of RdCHS1 in wild-type, mutant and transgenic lines. (C) Total contents of anthocyanins in Arabidopsis seedlings. Results correspond to means from three biological replicates. Asterisks indicate significant differences between means of wild-type and transgenic plants calculated by Tukey HSD test (***P<0.001).
Figure 6.
HPLC analyses of anthocyanins in Arabidopsis seedlings. HPLC chromatograms of the samples from seedlings of wild-type (A), mutant (B), RdCHS1-1 (C) and RdCHS1-3 (D) .
Figure 6.
HPLC analyses of anthocyanins in Arabidopsis seedlings. HPLC chromatograms of the samples from seedlings of wild-type (A), mutant (B), RdCHS1-1 (C) and RdCHS1-3 (D) .
Figure 7.
Effect of RdCHS1 on anthocyanin accumulation in transgenic tobacco flowers. (A) Tobacco flowers of wild-type and transgenic lines. (B) Expression profiles of RdCHS1 in flowers of transgenic tobacco. (C) Quantitation of anthocyanin accumulation levels in transgenic tobacco flowers with HPLC. (D) HPLC chromatograms of the samples from flowers of wild-type and transgenic tobacco. Results correspond to means from three biological replicates. Asterisks indicate significant differences between means of wild-type and transgenic plants calculated by Tukey HSD test (**P<0.01; *P<0.05).
Figure 7.
Effect of RdCHS1 on anthocyanin accumulation in transgenic tobacco flowers. (A) Tobacco flowers of wild-type and transgenic lines. (B) Expression profiles of RdCHS1 in flowers of transgenic tobacco. (C) Quantitation of anthocyanin accumulation levels in transgenic tobacco flowers with HPLC. (D) HPLC chromatograms of the samples from flowers of wild-type and transgenic tobacco. Results correspond to means from three biological replicates. Asterisks indicate significant differences between means of wild-type and transgenic plants calculated by Tukey HSD test (**P<0.01; *P<0.05).
Figure 8.
Expression analyses of endogenous anthocyanin biosynthetic genes in corollas of transgenic tobacco. (A) Expression profiles of structure genes in corollas of transgenic tobacco. (B) Expression profiles of regulatory genes in corollas of transgenic tobacco. Results represent means ± SE from three biological replicates. Different letters above the bars indicate significant difference between the samples judged by Tukey HSD test (P<0.01).
Figure 8.
Expression analyses of endogenous anthocyanin biosynthetic genes in corollas of transgenic tobacco. (A) Expression profiles of structure genes in corollas of transgenic tobacco. (B) Expression profiles of regulatory genes in corollas of transgenic tobacco. Results represent means ± SE from three biological replicates. Different letters above the bars indicate significant difference between the samples judged by Tukey HSD test (P<0.01).