https://orcid.org/0000-0001-8663-3041
Figures
Abstract
Pennisetum setaceum ‘Rubrum’ is an ornamental grass plant that produces purple leaves in high-light environments and light purple or green leaves in low-light environments, the latter of which greatly reduces its aesthetic appeal. Therefore, we aimed to identify the key genes associated with leaf coloration and elucidate the molecular mechanisms involved in the color changes in P. setaceum ‘Rubrum’ leaves. We performed transcriptome sequencing of P. setaceum ‘Rubrum’ leaves before and after shading. A total of 19,043 differentially expressed genes were identified, and the numbers of upregulated and downregulated genes at T1 stage, when compared with their expression at the T0 stage, were 10,761 and 8,642, respectively. The possible pathways that determine P. setaceum ‘Rubrum’ leaf color included flavonoid biosynthesis, flavone and flavonol biosynthesis, and carotenoid biosynthesis. There were 31 differentially expressed genes related to chlorophyll metabolism, of which 21 were related to chlorophyll biosynthesis and 10 to chlorophyll degradation, as well as three transcription factors that may be involved in the regulation of chlorophyll degradation. There were 31 key enzyme genes involved in anthocyanin synthesis and accumulation in P. setaceum ‘Rubrum’ leaves, with four transcription factors that may be involved in the regulation of anthocyanin metabolism. The transcriptome data were verified and confirmed reliable by real-time fluorescence quantitative PCR analysis. These findings provide a genetic basis for improving leaf color in P. setaceum ‘Rubrum.’
Citation: Zhu T, Wang X, Xu Z, Xu J, Li R, Liu N, et al. (2020) Screening of key genes responsible for Pennisetum setaceum ‘Rubrum’ leaf color using transcriptome sequencing. PLoS ONE 15(11): e0242618. https://doi.org/10.1371/journal.pone.0242618
Editor: Haitao Shi, Hainan University, CHINA
Received: August 9, 2020; Accepted: November 5, 2020; Published: November 23, 2020
Copyright: © 2020 Zhu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: Sequencing data are uploaded to the NCBI repository under BioProject PRJNA669318 (http://www.ncbi.nlm.nih.gov/bioproject/669318). All other relevant data are within the paper and its Supporting Information files.
Funding: This research was financially supported by the Forest Park Engineering Center open project of State Forestry Administration, China (Grant No. PTJH1500210) and the Fujian Agriculture and Forestry University Landscape Architecture School discipline construction project (Grant No.YSYL-bdpy-6). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
In recent years, ornamental grass with colorful leaves and novel morphologies has become popular for landscaping and gardening. Owing to its strong ecological adaptability and unique artistic use, ornamental grass adds freshness to existing urban plant design. Ornamental grass with colored leaves, especially deep reds, purples, and other warm tones can create a mysterious or cozy atmosphere [1]. P. setaceum ‘Rubrum’ is a cultivated grass with an elegant shape and light and delicate panicles. Its natural and simple appearance brings unique beauty and draws attention to the landscaping design. Plant color, as an essential factor of sensory quality, has a noticeable effect on ornamental value and commercial value [2]. Thus, the purplish-red leaves of P. setaceum ‘Rubrum’ contribute to its value; however, in our previous study, we show that when P. setaceum ‘Rubrum’ is grown in the shade, its leaf color turns green and the pigment content changes [3].
Chlorophyll and anthocyanins are the primary pigments that determine plant leaf color. Chlorophyll gives leaves their green color, while anthocyanin is the primary pigment responsible for red, purple, and blue colors [4]. Chloroplasts produce various pigments, including carotenoids and chlorophyll, which are the main pigment components in regular green leaves [4]. At present, most research has been mainly focused on the biosynthesis and degradation of chlorophyll. In the model plant Arabidopsis thaliana, 27 genes encode the 15 enzymes involved in chlorophyll biosynthesis from glutamyl-tRNA(Glu) to chlorophyll b [5]. The chlorophyll biosynthesis pathway starts with glutamyl-tRNA, progresses through a series of reactions catalyzed by glutamyl-tRNA reductase (HemA), Mg-chelatase, light-dependent NADPH: rotochlorophyllide oxidoreductase (POR), chlorophyllide a oxygenase (CAO), and other enzymes before it ends with chlorophyll a and chlorophyll b as products [5, 6]. The anthocyanin biosynthesis pathway is a branch of the flavonoid biosynthesis pathway in many plants [7], and most of the genes in this pathway have been identified and cloned [8]. The main enzymes related to anthocyanin biosynthesis are phenylalanine ammonia-lyase (PAL), chalcone isomerase (CHI), dihydroflavonol-4-reductase (DFR), and UDPG-flavonoid 3-O-glucosyltransferase (UFGT). Among these, PAL is a key enzyme in the metabolic pathway of phenylpropionic acid in plants and plays a significant role in the biosynthesis of lignin and flavonoids [9]. CHI catalyzes the stereospecific and intramolecular isomerization of naringenin chalcone into corresponding (2S)-flavanones, and improves the conversion rate of the isomerization reaction [10]. DFR can reduce dihydrokaempferol to produce leucopelargonidin [11, 12], and UFGT, which is the last enzyme in anthocyanin biosynthesis, converts unstable anthocyanins into stable anthocyanins [13]. Transcription factors (TFs) involved in the anthocyanin biosynthetic pathway mainly belong to the three protein families, R2R3-MYB, basic helix loop helix (bHLH), and WD40 repeat (WDR) [14, 15]. Anthocyanin biosynthesis is mainly regulated by the interaction between two or three of these proteins [14–16].
Among the external environmental factors, light intensity is the most crucial factor affecting leaf color in plants [17]. Under different light intensities and light durations, the contents of photosynthetic pigments and anthocyanins in plant leaves change, thus affecting the color of leaves. There have been reports on the effects of light on the color of a model plant, fruits, and several woody plants. For example, in the model plant Petunia, anthocyanin synthesis was repressed under shade conditions [18]. In fruit studies, light significantly increased the accumulation of flavonoids and the expression of biosynthetic genes in grape skin, while dark treatment decreased the accumulation of flavonoids in grape skin [19, 20]. Studies on woody plants of crabapple showed that prolonged sunlight exposure was beneficial to the accumulation of anthocyanins in the leaves of crabapple trees, resulting in red coloration [21]. The effect of light on anthocyanin accumulation at the molecular level is concentrated in model plants, while the mechanism of anthocyanin accumulation and red leaf coloration in most woody and herbaceous plants remains elusive [21]. Beckwith et al. used high-performance liquid chromatography and nuclear magnetic resonance to study the effects of different light intensities and various light sources on anthocyanin accumulation in the leaves of P. setaceum ‘Rubrum’ [2]; nonetheless, the molecular mechanisms of these correlations are rarely explored.
In our previous studies, we found that shading followed by light treatment significantly increased the values of the leaf color parameters L* (lightness) and b* (red-greenness), while a* values (blue-yellowness) decreased [3]. The physiological indexes chlorophyll a and chlorophyll b content significantly increased, whereas anthocyanin, flavonoid, soluble sugar content, and PAL and CHI activities significantly decreased [3]. We also found that the activities of PPO and POD did not change significantly [3]. This work indicated that light conditions could alter both chlorophyll and anthocyanin biosynthesis, although the molecular mechanism of leaf color changes still remains unclear. Therefore, the aims of the present study were to elucidate the molecular mechanisms that regulate leaf color parameters and physiological indicators. We performed transcriptome sequencing of leaves before and after shading and screened for key genes that control leaf coloration. This study will facilitate the horticulture and ornamental quality of P. setaceum ‘Rubrum’ and may also provide valuable information for further analysis of the molecular mechanisms involved in leaf coloration in herbaceous plant species.
Materials and methods
Plant materials and treatments
Seedlings of P. setaceum ‘Rubrum’ were transplanted from the ornamental grass base in Minqing County (Fuzhou City, Fujian Province) to the experimental base of Fujian Agriculture and Forestry University (Fujian Province, China) in mid-March 2018 and cultivated in open flower pots. The experiment was conducted in late June which is in the summer season for the region. The plants were shaded with a sunshade net that had a 10% transmittance, which was measured with an illuminance meter (transmittance = indoor illuminance/open—air illuminance). After the new leaves turned completely green, the sunshade net was removed and the plants were exposed to light until the leaf color returned to a purple color. On the basis of the leaf color change in the pre-experiment (Fig 1) [3], we selected T0, T1, T2, T3, and T4 time points to determine the leaf color parameters and physiological indexes: T0, untreated state; T1, stage when new leaves became completely green after 12 days of shading; T2, one day after the sunshade nets were removed; T3, three days after the sunshade nets were removed; and T4, stage at which the leaf color completely reverted to purple after 6 days in full light. The samples collected at stages T0 and T1 were used for transcriptome sequencing. To improve the reliability of the experimental results, the experiment was conducted in triplicate, with each experiment consisting of six pots. Leaf samples collected from the six pots were pooled together. The third or fourth leaf from the top of the plant was sampled at 9:00 every morning. The collected samples were immediately frozen in liquid nitrogen and stored at -80°C for further use.
The three pictures in each row are three different plants randomly selected in each stage.
Experimental method
There were six mixed samples in the T0 and T1 groups. Total RNA was extracted using an RNAprep Pure Plant RNA Purification kit (Tiangen Biotech, Beijing, China) following the manufacturer’s instructions. The integrity of RNA was verified by RNase free agarose gel electrophoresis and the concentration was measured using an Agilent Biological Analyzer 2100 system (Agilent Technologies, Palo Alto, CA, USA). The Illumina Hi-Seq 4000 platform was used for sequencing, and the read length was 150 bp. Low-quality reads or reads containing adapter sequences were removed before bioinformatics analysis. The quality of transcriptome sequencing was confirmed by sequencing saturation and unigene coverage before conducting the next sequencing. Transcriptional group assembly was completed using Trinity software (http://trinityrnaseq.sf.net) [22]. First, the samples from the T0 and T1 stages were used to verify their reproducibility. A comprehensive analysis was performed using Principal Component Analysis and Pearson correlation coefficient (r2) as evaluation indices. Gene expression was calculated by reads per kilobase of transcript per million mapped reads (RPKM). Differentially expressed genes (DEGs) were screened according to the general filtering criteria of edgeR (|1og2Fold Change| > 1 and false discovery rate [FDR] < 0.05) [23]. The gene ontology analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs in the T0 and T1 stages were based on the methods described by Young et al. [24] and Wu et al. [25], respectively, with a q ≤ 0.05 threshold value. Transcription factors were predicted based on PlantTFDB (http://planttfdb.cbi.pku.edu.cn) [26]. Ten key DEGs related to leaf color and two internal reference genes were selected for further analysis, and the reliability of the transcriptome results was verified by real-time quantitative PCR (qPCR) analysis [27]. Sequencing data are uploaded to the NCBI repository under BioProject PRJNA669318 (http://www.ncbi.nlm.nih.gov/bioproject/669318).
Results
Transcriptomic analysis of purple and green P. setaceum ‘Rubrum’ leaves
The average Q20 and Q30 scores of RNA from samples collected in the T0 and T1 stages and sequenced on the Illumina Hi-Seq platform after filtration were greater than 98.35% and 94.82%, respectively, and had a GC content of approximately 56.20% (Table 1). The sequencing saturation (S1 Fig) and gene coverage (S2 Fig) results of each sample confirmed the good quality of this transcriptome sequencing, and the next step could proceed. A total of 93,750 unigenes were obtained using Trinity software, with N50 of 1362 nt and the average length of 916 nt; the overall assembly effect was good (Table 2). A total of 51,735 genes were annotated using four databases (Fig 2), namely the NCBI nonredundant protein database (Nr), Clusters of EuKaryotic Orthologous Groups of proteins database (KOG), KEGG, and Swiss-Prot database. The annotation success rate was the highest in the Nr database (55.00%) and the lowest in the KEGG database (19.23%).
Numbers indicate annotated unigenes.
Identification of differentially expressed genes in purple and green P. setaceum ‘Rubrum’ leaves
Principal component and r2 analyses were used to evaluate sample repeatability and to test reliability. The first and second principal component explained 96.4% of the overall variance (Fig 3A), and the repeatability test showed a correlation greater than 0.98 between each sample (Fig 3B). The sample selection was reasonable, and the sample repeatability and test reliability were high. A volcano plot of 19,403 DEGs is presented in Fig 3C. There were 10,761 upregulated (red) and 8642 downregulated (blue) unigenes in T1 relative to their expression in T0.
(A) Principal components analysis. (B) r2 analysis to test repeatability. (C) Statistics of differentially expressed genes. Blue represents downregulated unigenes, and red represents upregulated unigenes.
KEGG enrichment analysis identifies biosynthetic pathways responsible for leaf color
KEGG pathway analysis revealed that 9839 unigenes were associated with 131 pathways. A total of 386 unigenes were associated with metabolic pathways related to leaf color (Table 3). The main metabolic pathways related to leaf color were porphyrin and chlorophyll metabolism, photosynthesis, flavonoid biosynthesis, and carotenoid biosynthesis. The genes related to photosynthesis (27 genes) and anthocyanin (2 genes) were the least represented in this analysis.
A total of 1976 DEGs annotated in the KEGG database were significantly enriched in 34 pathways. Biosynthesis of secondary metabolites was the most significantly enriched pathway, followed by metabolic pathways and photosynthesis (Fig 4A). Among the pathways related to plant pigmentation, 33 DEGs were enriched in flavonoid biosynthesis pathway, four DEGs were in the flavone and flavonol biosynthesis pathway, 34 DEGs were in the carotenoid biosynthesis pathway, and two DEGs were in the anthocyanin biosynthesis pathway. These results suggested that leaf coloration in P. setaceum ‘Rubrum’ is associated with changes in these metabolic pathways.
(A) All DEGs enriched in KEGG pathways; (B) Pathway enrichment of upregulated DEGs; (C) Pathway enrichment of downregulated DEGs.
Most of the DEGs that were upregulated in the T1 stage (Fig 4B) were enriched in the following pathways: photosynthesis, porphyrin and chlorophyll metabolism, photosynthesis-antenna proteins, and carbon fixation in photosynthetic organisms. Among the downregulated genes (Fig 4C), DEGs were significantly enriched in flavonoid biosynthesis, carotenoid biosynthesis, phenylalanine metabolism, and phenylpropanoid biosynthesis.
Genes related to chlorophyll metabolism.
A previous study found that chlorophyll, anthocyanins, flavonoids, and PAL and CHI activities are closely related to leaf color [3]. Therefore, in the present study, we focused on KEGG database annotation and prediction of TFs to investigate the regulatory relationship of these structural genes and TFs on leaf color in P. setaceum ‘Rubrum.’
Based on the pathway analysis of the DEGs, 31 DEGs encode enzymes that regulate chlorophyll metabolism (Fig 5A and 5B). In order to analyze the pathways modulated by chlorophyll-related DEGs, we analyzed chlorophyll metabolism pathway and heat map of P. setaceum ‘Rubrum’. Twenty-one unigenes were related to chlorophyll synthesis and 10 to chlorophyll degradation (Fig 5C). Among the DEGs involved in chlorophyll biosynthesis (Fig 5A), six genes were predicted to belong to the HemA enzyme family, which functions in the first stage of chlorophyll biosynthesis. In the second stage of chlorophyll biosynthesis, two porphobilinogen deaminase (HemC) genes were predicted to encode porphobilinogen (PBG). In the last stage of chlorophyll biosynthesis, three genes were annotated to the magnesium chelatase subunit ChlD, two genes to magnesium chelatase subunit ChlH, and one gene to magnesium-chelatase subunit ChlI of the Mg-chelatase enzyme. Five genes were annotated to POR, one gene was annotated to chlorophyll synthase (ChlG), and one gene was annotated to CAO. These DEGs are the key regulatory enzymes in chlorophyll biosynthesis and were upregulated at the T1 stage.
(A) Chlorophyll biosynthesis pathway [5, 6], (B) chlorophyll degradation pathway [28–30], and (C) heat map of pathway-related differentially expressed genes. The red box indicates upregulated genes, the blue box indicates downregulated genes, and the orange box indicates genes that are both upregulated and downregulated in T1 stage.
A gene annotated to NON-YELLOWCOLORING1 (NYC1), which is related to chlorophyll degradation (Fig 5B), was downregulated, whereas a gene annotated to NYC1-LIKE (NOL) was upregulated in the T1 stage. There were three downregulated genes in the chlorophyllase (CLH) family. A total of five genes that were related to pheophorbide a oxygenase (PaO) were downregulated in the T1 stage.
Based on PlantTFDB [26], the annotation results for the TFs predicted five genes that regulate chlorophyll biosynthesis (Table 4). These five TFs were all annotated in the phytochrome interacting factors (PIFs) family. According to the gene annotation, as long as the genes coding helicase/helix were downregulated, the chlorophyll content in the T1 stage was higher than that in T0, indicating that the helicase encoded by these genes may lead to chlorophyll degradation.
Genes related to anthocyanin metabolism involved in P. setaceum ‘Rubrum’ leaf color.
The biosynthesis of the leaf pigment anthocyanin is a very complex process that includes several metabolic pathways, including phenylalanine (Fig 6A), flavonoid (Fig 6B), and anthocyanin (Fig 6C) metabolic pathways [31, 32]. Based on the pathway analysis of the DEGs, a total of 31 key enzyme-encoding genes may be involved in the synthesis and accumulation of anthocyanins in P. setaceum ‘Rubrum’ leaves. Of the seven genes annotated to the PAL family, only two genes were upregulated (Fig 6D). Nine genes annotated to chalcone synthase (CHS), one gene annotated to CHI, one gene annotated to flavonol synthase (FLS), and one gene annotated to DFR were downregulated. Within the anthocyanidin reductase (ANR) family, which has nine annotated genes, only two genes were upregulated in the T1 stage. One gene annotated to leucoanthocyanidin dioxygenase (LDOX) was downregulated (Fig 6D). In the anthocyanin metabolism pathway, two downregulated genes were annotated to the UFGT family downstream of anthocyanin biosynthesis (Fig 6D).
The first (A), second (B), and third (C) stages of anthocyanin biosynthesis [31, 32]. The blue box indicates downregulated genes, and the orange box indicates both upregulated and downregulated genes.
We predicted the genes encoding three protein families of anthocyanin synthesis, namely R2R3-MYB, bHLH and WDR, and four genes were predicted to be involved in P. setaceum ‘Rubrum’ leaf color formation (Table 5). A total of 34 genes were differentially expressed in T0 and T1 stages, two of which were R2R3-MYB TFs and related to anthocyanin. Among them, the expression of MYB1 was upregulated, resulting in the inhibition of anthocyanin synthesis. The annotation of the C1 gene indicated it encodes an anthocyanin regulatory protein whose downregulation promotes anthocyanin synthesis. A total of 33 differentially expressed genes were predicted in the bHLH TF family, and two genes that were related to anthocyanin were annotated as R-S: one was upregulated, while the other one was downregulated. From a total of 43 genes from the WD40 family of TFs, four were differentially expressed and downregulated in T1, but none of these genes were related to anthocyanins. These results suggested that the leaf color in P. setaceum ‘Rubrum’ may be jointly regulated by a complex formed by bHLH and R2R3-MYB.
Real-time quantitative PCR validation of genes identified in the P. setaceum ‘Rubrum’ transcriptome
To further verify the reliability of the transcriptome sequencing, we randomly selected 10 DEGs closely related to chlorophyll metabolism and anthocyanin metabolism for verification. MDH and ACP were selected as internal reference genes based on the results of the gene expression study on Pennisetum glaucum [33]. The dissolution curves of the 12 genes were unimodal. The results of qPCR detection of these genes in samples from T0 and T1 were consistent with the upregulated and downregulated expression in the transcriptome dataset (Fig 7), confirming the reliability of the transcriptomic sequencing.
The log2 Ratio(T1/T0) represents the relative gene expression.
Discussion
KEGG pathways for chlorophyll and anthocyanin biosynthesis are associated with leaf color changes in response to light condition
Transcriptome sequencing is a highly efficient and affordable method for the fast generation of nearly all cellular transcripts under different physiological conditions. It has been widely used in studies using A. thaliana, maize, and rice, as well as ornamental plants such as Lagerstroemia [34] and Lilium ‘Tiny Padhye’ [35].
The KEGG database is used to analyze the metabolic pathways and functions of gene products to elucidate the complex biological behaviors of genes. The KEGG enrichment analysis conducted in the present study was consistent with our previous results [3]. The red-greenness a* values of the color index decreased significantly in the T1 stage, when the leaf color reached a blue-yellowness grade. Physiological experiments showed increased chlorophyll content at T1 stage [3], suggesting that photosynthesis was enhanced, which was verified with the results that showed the upregulated DEGs were mainly enriched in KEGG pathways related to chlorophyll and photosynthesis. Similarly, the downregulated DEGs were enriched in the pathways related to anthocyanin biosynthesis in the T1 stage, indicating a decrease in anthocyanin content, which was consistent with the previous leaf color parameters and physiological experimental results [3].
Expression of genes related to chlorophyll metabolism regulate leaf color changes in response to light
Chlorophyll is the pigment responsible for the green color of plant leaves, and its biosynthesis is well studied [5, 6]. HemA is the first enzyme of the chlorophyll biosynthesis pathway and catalyzes the transformation of L-glutamyl-tRNA to δ-glutamyl RNA [6]. HemC genes encode porphobilinogen deaminase, which functions in the second stage of chlorophyll biosynthesis. In the last stage of chlorophyll biosynthesis, Mg-chelatase is the critical enzyme at the branching point in the chlorophyll biosynthesis pathway. The three subunits of the Mg-chelatase enzyme are magnesium chelatase subunit ChlD, magnesium chelatase subunit ChlH, and magnesium-chelatase subunit ChlI [36]. After the shading treatment, the P. setaceum ‘Rubrum’ leaves turned green. We found that the expression of six HemA, two HemC, three ChlD, two ChlH, one ChlI, five POR, one ChlG, and one CAO gene was upregulated in the T1 stage, which resulted in an increase in the chlorophyll content of the leaves, promoting their green coloration. This finding was consistent with the results of previous physiological experiments in which chlorophyll content in the T1 stage was significantly higher than that in the T0 stage, consistent with the results of leaf color parameters [3]. These results all suggest the expression of genes related to chlorophyll synthesis is sensitive to light intensity and affects the dynamic color change observed in P. setaceum ‘Rubrum’ leaves.
The above genes that play a key role in chlorophyll synthesis have been studied in other plants. A study on yellow leaves of Lagerstroemia indica found that downregulation of genes such as HemA and HemD hindered normal chlorophyll biosynthesis [37]. Sakuraba et al. studied faded green leaf mutants in rice and found that POR is an essential gene for chlorophyll biosynthesis in rice leaves [38]. The literature shows that HemA, HemD, and POR are critical genes in chlorophyll biosynthesis. Combined with the results of previous studies, the results presented herein confirm the reliability of our conjecture about the key genes in the P. setaceum ‘Rubrum’ chlorophyll synthesis process.
In the chlorophyll degradation process, NYC1 and NOL catalyze the degradation of chlorophyll b to chlorophyll a [39], while the CLH enzyme catalyzes the removal of the phytol from chlorophyll a to produce chlorophyllide [40]. PaO is another enzyme involved in the chlorophyll degradation pathway. Our analysis identified 10 DEGs that encode enzymes related to chlorophyll degradation, namely one NOL, one NYC1, three CHL, and five PaO genes. These genes were significantly downregulated in the T1 stage and the degradation rate of chlorophyll was reduced. Consequently, the chlorophyll content in leaves from the T1 stage was greater than that in stage T0. These results were consistent with previous leaf color parameters and physiological observations [3], further suggesting these genes are related to chlorophyll degradation and leaf color changes in P. setaceum ‘Rubrum’ in response to changes in light conditions.
Previous studies have explored the degradation of chlorophyll in different plants. Sato et al. studied the chlorophyll degradation pathway during rice senescence and found that nol and nyc1 rice mutants showed a stay-green phenotype and chlorophyll b degradation was severely inhibited, which indicated that NOL and NYC1 promoted chlorophyll degradation [39]. The results of present study show that at the T1 stage, one NYC1 was downregulated to promote chlorophyll degradation, while one NYC1-LIKE (NOL) gene was upregulated to inhibit chlorophyll degradation, indicating that the regulation of NOL on leaf color may be more complex in P. setaceum ‘Rubrum’. Takamiya et al. showed that both CLH and PaO are essential enzymes for chlorophyll degradation in plants [30], while Schenk et al. found that the chlorophyllases genes AtCLH1 and AtCLH2 are not required for chlorophyll decomposition during leaf senescence in A. thaliana [41]. The results of the present study suggest that CLH and PaO expression is related to chlorophyll degradation and that their gene products are likely necessary for changes in P. setaceum ‘Rubrum’ leaf color.
Among the TFs related to chlorophyll biosynthesis, five regulatory genes that are members of the PIFs family were predicted to be related to chlorophyll degradation. In the PIFs family, we predicted that the helicases/helix genes would be downregulated as chlorophyll content increased, suggesting that these genes promote chlorophyll degradation. Huq et al. reported a negative correlation between PIF1 and chlorophyll biosynthesis in plants [42], and Shin et al. found that PIF3 reduced the rate of chlorophyll biosynthesis and photosynthetic efficiency by inhibiting the expression of GUN5 and other essential genes for chlorophyll biosynthesis in the dark [43]. We found that under shading conditions, the downregulation of PIFs promotes the synthesis of chlorophyll. Previous research and our results show that the PIFs family can regulate the biosynthesis of chlorophyll.
Genes differentially regulated by light are related to anthocyanin metabolism
PAL is an essential enzyme for the first stage of anthocyanin biosynthesis, where it catalyzes L-phenylalanine to form trans-cinnamic acid [31]. CHS and CHI are two important enzymes in the flavonoid biosynthesis pathway. Flavonols are a side product of naringenin chalcone produced in a reaction catalyzed by FLS [32]. The DFR enzyme can use dihydroflavonols, dihydromyricetin, dihydrokaempferol, and/or dihydroquercetin as substrates to synthesize leucoanthocyanidins [11]. After the shading treatment, the leaf color of P. setaceum ‘Rubrum’ changed from purple to green, which correlated with the significant downregulation of anthocyanin biosynthesis genes. These included seven PAL, nine CHS, one FLS, one CHI, one DFR, nine ANR, one LDOX, and two UFGT genes. The transcriptome sequencing data for anthocyanin biosynthesis pathway members showed that most of these genes were downregulated in the T1 stage, which was consistent with results from physiological experiments on flavonoid content, anthocyanin content, PAL enzyme activity, and CHI enzyme activity previously reported [3]. When the genes encoding key enzymes in the anthocyanin biosynthesis pathway were downregulated, the leaves in the T1 stage lacked purple or red coloration, which was also consistent with previous leaf color parameter results that confirmed the purplish-red leaves in the T0 stage and green leaves in the T1 stage [3].
Similar conclusions were reported in studies on flower and pericarp color. For example, Sun et al. found that the expression pattern of FhCHS1 in Freesia hybrida flowers was significantly related to the accumulation patterns of anthocyanins during flower development [44]. Furthermore, the function of FhCHS1 was verified by ectopic expression in Petunia hybrida, and it was found that CHS1 was responsible for the flower color transition from the original white to pink [44]. Christopher et al. reported higher anthocyanin content in orange than in yellow Clivia miniata, and observed the same trend for the expression of CHS and DFR [45]. Zhao et al. found that the transcription of FLS was downregulated throughout apple fruit development under reduced UV [46]. Zhao et al. showed that the red peel color of Litchi chinensis fruits is attributed to UFGT and its role in anthocyanin biosynthesis [47]. Therefore, we believe that these genes play an important role in the purple leaf phenotype of P. setaceum ‘Rubrum.’
The MYB-bHLH-WD40 (MBW) transcription complexes mainly regulate anthocyanin biosynthesis [13, 14]. However, the results of our study showed that during the process of leaf color change from purple to green in P. setaceum ‘Rubrum,’ the DEGs related to anthocyanins were not annotated with WD40, whereas the expression of anthocyanin genes in MYB and bHLH TFs changed significantly. Therefore, we hypothesized that the purple leaf color of P. setaceum ‘Rubrum’ is determined by the complex of bHLH and MYB family TFs. In our study, C1 of the R2R3-MYB gene family and one R-S gene from the bHLH gene family were downregulated, which was positively correlated with anthocyanin synthesis. Shin et al. expressed maize C1 and R-S genes in rice endosperm, and the resulting kernels produced several different kinds of flavonoids [48], which was similar to our results. Therefore, we hypothesized that C1 and R-S can activate structural genes in the flavonoid biosynthesis pathway individually or synergistically. However, our sequencing data showed that only one R-S gene in the bHLH gene family was upregulated, suggesting its role in the inhibition of anthocyanin synthesis. This may be because the R2R3-MYB gene family plays a key role in the regulation of anthocyanin synthesis mediated by the MBW complex [49]. Hence, the upregulated expression of one R-S gene has a limited effect on anthocyanin synthesis, although the specific mechanisms of these gene products warrant further study. The expression of MYB1 from the MYB family of TFs was upregulated after the leaf color turned green, indicating that MYB1 may negatively regulate anthocyanin synthesis in P. setaceum ‘Rubrum.’ Most of the R2R3-MYB genes in Arabidopsis positively regulate anthocyanin synthesis, while a few genes are negative regulators of this pathway [50]. For example, Schwinn found that MYB1 positively regulates anthocyanin biosynthesis in onion [51], which is contradictory to the results of the present study. Thus, there may exist different mechanisms in different species by which MYB1 regulates anthocyanin biosynthesis.
Conclusions
In this study, Illumina RNA-Seq was used to conduct transcriptome sequencing analysis, which was used to assemble 93,750 unigenes. Genetic analysis based on the Illumina sequencing was used to identify DEGs that revealed many candidate genes related to leaf color, including genes involved in the metabolic pathway of chlorophyll and anthocyanins. This transcriptome analysis provides important information for the continued study of leaf color-related genes and offers not only a basis for future functional genomics research in P. setaceum ‘Rubrum,’ but also resources for genomics research in other P. setaceum cultivars. The study is of considerable significance for horticulture practices and can help to improve the ornamental quality of P. setaceum ‘Rubrum.’
Supporting information
S1 Fig. Saturation of six samples of Pennisetum setaceum ‘Rubrum’.
https://doi.org/10.1371/journal.pone.0242618.s001
(TIF)
S2 Fig. Statistics of gene coverage for six samples of Pennisetum setaceum ‘Rubrum’.
https://doi.org/10.1371/journal.pone.0242618.s002
(TIF)
Acknowledgments
Thanks to all the authors for their contributions to plant material care, data analysis and helpful suggestions in manuscript preparation. We would like to thank Editage (www.editage.cn) for English language editing and GENE DENOVO (www.genedenovo.com) for Illumina sequencing.
References
- 1. Klinkenborg V. Grasses: Versatile partners for uncommon garden design. Ny Times Bk Rev. 2002:13-+. WOS:000175885600051.
- 2. Beckwith AG, Zhang YJ, Seeram NP, Cameron AC, Nair MG. Relationship of light quantity and anthocyanin production in Pennisetum setaceum Cvs. Rubrum and Red Riding Hood. J Agr Food Chem. 2004;52(3):456–61. WOS:000188830300009. pmid:14759132
- 3. Zhu T ZL, Xu ZM, et al. Effects of Shading and Regaining Natural Light Treatment on Leaf Color and Physiological Changes of Pennisetum setaceum ‘Rubrum’[J]. ACTA AGRESTIA SINICA. 2019;27(01):266–271.
- 4. Yang YX, Chen XX, Xu B, Li YX, Ma YH, Wang GG. Phenotype and transcriptome analysis reveals chloroplast development and pigment biosynthesis together influenced the leaf color formation in mutants of Anthurium andraeanum 'Sonate'. Front Plant Sci. 2015;6. WOS:000351007000001.
- 5. Nagata N, Tanaka R, Satoh S, Tanaka A. Identification of a vinyl reductase gene for chlorophyll synthesis in Arabidopsis thaliana and implications for the evolution of Prochlorococcus species. Plant Cell. 2005;17(1):233–40. WOS:000226317000018. pmid:15632054
- 6. Beale SI. Green genes gleaned. Trends Plant Sci. 2005;10(7):309–12. WOS:000230906100003. pmid:15951223
- 7. Jaakola L. New insights into the regulation of anthocyanin biosynthesis in fruits. Trends Plant Sci. 2013;18(9):477–83. WOS:000325058400003. pmid:23870661
- 8. Tanaka Y, Ohmiya A. Seeing is believing: engineering anthocyanin and carotenoid biosynthetic pathways. Curr Opin Biotech. 2008;19(2):190–7. WOS:000255982100017. pmid:18406131
- 9. Lister CE, Lancaster JE, Walker JRL. Phenylalanine ammonia-lyase (PAL) activity and its relationship to anthocyanin and flavonoid levels in New Zealand-grown apple cultivars. J Am Soc Hortic Sci. 1996;121(2):281–5. WOS:A1996TX78200020.
- 10. Cheng AX, Zhang XB, Han XJ, Zhang YY, Gao S, Liu CJ, et al. Identification of chalcone isomerase in the basal land plants reveals an ancient evolution of enzymatic cyclization activity for synthesis of flavonoids. New Phytol. 2018;217(2):909–24. WOS:000419324000038. pmid:29083033
- 11. Murray JR, Smith AG, Hackett WP. Differential Dihydroflavonol Reductase Transcription and Anthocyanin Pigmentation in the Juvenile and Mature Phases of Ivy (Hedera-Helix L). Planta. 1994;194(1):102–9. WOS:A1994NR40200014.
- 12. Xie S, Zhao T, Zhang ZW, Meng JF. Reduction of Dihydrokaempferol by Vitis vinfera Dihydroflavonol 4-Reductase to Produce Orange Pelargonidin-Type Anthocyanins. J Agr Food Chem. 2018;66(13):3524–32. WOS:000429508400031. pmid:29554804
- 13. Ju ZG, Liu CL, Yuan YB. Activities of Chalcone Synthase and Udpgal—Flavonoid-3-O-Glycosyltransferase in Relation to Anthocyanin Synthesis in Apple. Sci Hortic-Amsterdam. 1995;63(3–4):175–85. WOS:A1995TD53800005.
- 14. Koes R, Verweij W, Quattrocchio F. Flavonoids: a colorful model for the regulation and evolution of biochemical pathways. Trends Plant Sci. 2005;10(5):236–42. WOS:000229465500007. pmid:15882656
- 15. Xu WJ, Dubos C, Lepiniec L. Transcriptional control of flavonoid biosynthesis by MYB-bHLH-WDR complexes. Trends Plant Sci. 2015;20(3):176–85. WOS:000350933100007. pmid:25577424
- 16. Li CH, Qiu J, Yang GS, Huang SR, Yin JM. Isolation and characterization of a R2R3-MYB transcription factor gene related to anthocyanin biosynthesis in the spathes of Anthurium andraeanum (Hort.). Plant Cell Rep. 2016;35(10):2151–65. WOS:000383708500012. pmid:27424029
- 17. Albert NW, Lewis DH, Zhang H, Irving LJ, Jameson PE, Davies KM. Light-induced vegetative anthocyanin pigmentation in Petunia. J Exp Bot. 2009;60(7):2191–202. WOS:000266115300027. pmid:19380423
- 18. Albert NW, Lewis DH, Zhang H, Schwinn KE, Jameson PE, Davies KM. Members of an R2R3-MYB transcription factor family in Petunia are developmentally and environmentally regulated to control complex floral and vegetative pigmentation patterning. Plant J. 2011;65(5):771–84. WOS:000287701200009. pmid:21235651
- 19. Azuma A, Yakushiji H, Koshita Y, Kobayashi S. Flavonoid biosynthesis-related genes in grape skin are differentially regulated by temperature and light conditions. Planta. 2012;236(4):1067–80. WOS:000309229300010. pmid:22569920
- 20. Matus JT, Loyola R, Vega A, Pena-Neira A, Bordeu E, Arce-Johnson P, et al. Post-veraison sunlight exposure induces MYB-mediated transcriptional regulation of anthocyanin and flavonol synthesis in berry skins of Vitis vinifera. J Exp Bot. 2009;60(3):853–67. WOS:000264189900015. pmid:19129169
- 21. Lu YF, Zhang ML, Meng XN, Wan HH, Zhang J, Tian J, et al. Photoperiod and shading regulate coloration and anthocyanin accumulation in the leaves of malus crabapples. Plant Cell Tiss Org. 2015;121(3):619–32. WOS:000354098700009.
- 22. Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8(8):1494–512. WOS:000322228600003. pmid:23845962
- 23. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40. WOS:000273116100025. pmid:19910308
- 24. Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2). WOS:000276434300013. pmid:20132535
- 25. Wu QJ, Chen ZD, Sun WJ, Deng TT, Chen MJ. De novo Sequencing of the Leaf Transcriptome Reveals Complex Light-Responsive Regulatory Networks in Camellia sinensis cv. Baijiguan. Front Plant Sci. 2016;7. WOS:000372503600001. pmid:26870046
- 26. Guo AY, Chen X, Gao G, Zhang H, Zhu QH, Liu XC, et al. PlantTFDB: a comprehensive plant transcription factor database. Nucleic Acids Res. 2008;36:D966–D9. WOS:000252545400169. pmid:17933783
- 27. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE Guidelines: Minimum Information for Publication of Quantitative Real-Time PCR Experiments. Clin Chem. 2009;55(4):611–22. WOS:000264774200006. pmid:19246619
- 28. Hortensteiner S. Chlorophyll degradation during senescence. Annu Rev Plant Biol. 2006;57:55–77. Epub 2006/05/04. pmid:16669755.
- 29. Guyer L, Salinger K, Krugel U, Hortensteiner S. Catalytic and structural properties of pheophytinase, the phytol esterase involved in chlorophyll breakdown. J Exp Bot. 2018;69(4):879–89. WOS:000425379300014. pmid:29036670
- 30. Takamiya K, Tsuchiya T, Ohta H. Degradation pathway(s) of chlorophyll: what has gene cloning revealed? Trends Plant Sci. 2000;5(10):426–31. WOS:000089920700009. pmid:11044719
- 31. Achnine L, Blancaflor EB, Rasmussen S, Dixon RA. Colocalization of L-phenylalanine ammonia-lyase and cinnarnate 4-hydroxylase for metabolic channeling in phenylpropanoid biosynthesis. The Plant cell. 2004;16(11):3098–109. WOS:000225228500019. pmid:15472080
- 32. Liu Y, Tikunov Y, Schouten RE, Marcelis LFM, Visser RGF, Bovy A. Anthocyanin Biosynthesis and Degradation Mechanisms in Solanaceous Vegetables: A Review. Front Chem. 2018;6. WOS:000427037300001. pmid:29594099
- 33. Reddy PS, Reddy DS, Sharma KK, Bhatnagar-Mathur P, Vadez V. Cloning and validation of reference genes for normalization of gene expression studies in pearl millet [Pennisetum glaucum (L.) R. Br.] by quantitative real-time PCR. Plant Gene. 2015;1:35–42.
- 34. Li SZ, Zheng TC, Zhuo XK, Li ZJ, Li LL, Li P, et al. Transcriptome profiles reveal that gibberellin-related genes regulate weeping traits in crape myrtle. Hortic Res-England. 2020;7(1). WOS:000524501200016. pmid:32257240
- 35. Xu LF, Yang PP, Feng YY, Xu H, Cao YW, Tang YC, et al. Spatiotemporal Transcriptome Analysis Provides Insights into Bicolor Tepal Development in Lilium "Tiny Padhye". Front Plant Sci. 2017;8. WOS:000397279500001. pmid:28392796
- 36. Papenbrock J, Grafe S, Kruse E, Hanel F, Grimm B. Mg-chelatase of tobacco: identification of a Chl D cDNA sequence encoding a third subunit, analysis of the interaction of the three subunits with the yeast two-hybrid system, and reconstitution of the enzyme activity by co-expression of recombinant CHL D, CHL H and CHL I. Plant J. 1997;12(5):981–90. WOS:000070997800002. pmid:9418040
- 37. Li Y, Zhang ZY, Wang P, Wang SA, Ma LL, Li LF, et al. Comprehensive transcriptome analysis discovers novel candidate genes related to leaf color in a Lagerstroemia indica yellow leaf mutant. Genes Genom. 2015;37(10):851–63. WOS:000362020200006.
- 38. Sakuraba Y, Rahman ML, Cho SH, Kim YS, Koh HJ, Yoo SC, et al. The rice faded green leaf locus encodes protochlorophyllide oxidoreductase B and is essential for chlorophyll synthesis under high light conditions. Plant J. 2013;74(1):122–33. Epub 2013/01/08. pmid:23289852.
- 39. Sato Y, Morita R, Katsuma S, Nishimura M, Tanaka A, Kusaba M. Two short-chain dehydrogenase/reductases, NON-YELLOW COLORING 1 and NYC1-LIKE, are required for chlorophyll b and light-harvesting complex II degradation during senescence in rice. Plant J. 2009;57(1):120–31. WOS:000261962700010. pmid:18778405
- 40. Shemer TA, Harpaz-Saad S, Belausov E, Lovat N, Krokhin O, Spicer V, et al. Citrus chlorophyllase dynamics at ethylene-induced fruit color-break: A study of chlorophyllase expression, posttranslational processing kinetics, and in situ intracellular localization. Plant Physiol. 2008;148(1):108–18. WOS:000258947600011. pmid:18633118
- 41. Schenk N, Schelbert S, Kanwischer M, Goldschmidt EE, Dormann P, Hortensteiner S. The chlorophyllases AtCLH1 and AtCLH2 are not essential for senescence-related chlorophyll breakdown in Arabidopsis thaliana. Febs Lett. 2007;581(28):5517–25. WOS:000253487900031. pmid:17996203
- 42. Huq E, Al-Sady B, Hudson M, Kim C, Apel K, Quail PH. Phytochrome-interacting factor 1 is a critical bHLH regulator of chlorophyll biosynthesis. Science. 2004;305(5692):1937–41. Epub 2004/09/28. pmid:15448264.
- 43. Shin J, Kim K, Kang H, Zulfugarov IS, Bae G, Lee CH, et al. Phytochromes promote seedling light responses by inhibiting four negatively-acting phytochrome-interacting factors. P Natl Acad Sci USA. 2009;106(18):7660–5. WOS:000265783600071. pmid:19380720
- 44. Sun W, Meng XY, Liang LJ, Jiang WS, Huang YF, He J, et al. Molecular and Biochemical Analysis of Chalcone Synthase from Freesia hybrid in Flavonoid Biosynthetic Pathway. Plos One. 2015;10(3). WOS:000350688100077. pmid:25742495
- 45. Viljoen CD, Snyman MC, Spies JJ. Identification and expression analysis of chalcone synthase and dihydroflavonol 4-reductase in Clivia miniata. S Afr J Bot. 2013;87:18–21. WOS:000321163500004.
- 46. Henry-Kirk RA, Plunkett B, Hall M, McGhie T, Allan AC, Wargent JJ, et al. Solar UV light regulates flavonoid metabolism in apple (Malus x domestica). Plant Cell and Environment. 2018;41(3):675–88. WOS:000424841000016.
- 47. Zhao ZC, Hu GB, Hu FC, Wang HC, Yang ZY, Lai B. The UDP glucose: flavonoid-3-O-glucosyltransferase (UFGT) gene regulates anthocyanin biosynthesis in litchi (Litchi chinesis Sonn.) during fruit coloration. Mol Biol Rep. 2012;39(6):6409–15. WOS:000303351900002. pmid:22447536
- 48. Shin YM, Park HJ, Yim SD, Baek NI, Lee CH, An GH, et al. Transgenic rice lines expressing maize C1 and R-S regulatory genes produce various flavonoids in the endosperm. Plant Biotechnol J. 2006;4(3):303–15. WOS:000236584800004. pmid:17147636
- 49. Zhi JJ, Liu XX, Li DJ, Huang YGM, Yan SS, Cao BH, et al. CRISPR/Cas9-mediated SlAN2 mutants reveal various regulatory models of anthocyanin biosynthesis in tomato plant. Plant Cell Rep. 2020;39(6):799–809. WOS:000521916100001. pmid:32221665
- 50. Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010;15(10):573–81. Epub 2010/08/03. pmid:20674465.
- 51. Schwinn KE, Ngo H, Kenel F, Brummell DA, Albert NW, McCallum JA, et al. The Onion (Allium cepa L.) R2R3-MYB Gene MYB1 Regulates Anthocyanin Biosynthesis. Front Plant Sci. 2016;7. WOS:000443589700001. pmid:28018399