Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Transcriptome profiling of genes related to light-induced anthocyanin biosynthesis in eggplant (Solanum melongena L.) before purple color becomes evident

BMC Genomics volume 19, Article number: 201 (2018) Cite this article

A Correction to this article was published on 04 May 2018

This article has been updated

Abstract

Background

The anthocyanins are highly enriched in eggplants (Solanum melongena L.) with purple peel. However, our previous study showed that anthocyanins biosynthesis in eggplant cultivar ‘Lanshan Hexian’ was completely regulated by light and color becomes evident at most 2 days after exposure to light. In the present investigation, transcriptome study was made to explore the underlying molecular mechanisms of light-induced anthocyanin biosynthesis in eggplant (Solanum melongena L.) before color becomes evident.

Results

RNA-Seq was performed for four time points (0, 0.5, 4 and 8 h after bags removal) where concerted changes happened. A total of 32,630 genes or transcripts were obtained by transcriptome sequencing, from which 1956 differentially expressed genes (DEGs) were found. Gene Ontology analysis showed that the 1956 DEGs covered a wide range of cellular components, molecular functions and biological processes. All the DEGs were further divided into 26 clusters based on their distinct expression patterns. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis found out 24 structural anthocyanin biosynthesis genes which distributing in seven clusters. In addition, 102 transcription factors, which exhibited highly dynamic changes in response to light, were found in the seven clusters. Three photoreceptors, UV Resistance Locus 8 (UVR8), Cryptochrome 3 (CRY3) and UVR3, were identified as DEGs. The light signal transduction elements, COP1 and two SPAs, might be responsible for anthocyanin biosynthesis regulation.

Conclusion

Based on the transcriptome data, the anthocyanin biosynthesis structural genes, transcription factors, photoreceptors and light signal transduction elements were quickly screened which may act as the key regulatory factors in anthocyanin biosynthesis pathway. By comparing the transcriptome data with our previous studies, 869 genes were confirmed to participate in the light-induced anthocyanin biosynthesis. These results expand our knowledge of light-induced anthocyanin biosynthesis in plants, which allowing for fruit coloration to be improved under low-light conditions in future.

Background

Anthocyanins are natural pigments and widely found in plants. They not only confer plants vivid colors to attract pollinators and seed dispersers, but also enhance plants resistance to stress [1]. Anthocyanins are synthesized from phenylalanine, and then under the continuous catalytic by phenylalanine ammonia lyase (PAL), 4-coumarate-CoA ligase (4CL), chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxyl enzyme (F3H), dihydroflavonol reductase (DFR) and anthocyanin synthase (ANS) or leucoanthocyanidin dioxygenase (LDOX). This pathway creates unstable anthocyanins [2] that are stabilized by glycosylation (GT), methylation (MT) or acylation (AT) before transportation to the cell vacuole for storage [3]. Anthocyanins are mainly of six types: pelargonidin, cyanidin, delphinidin, peonidin, petunidin and malvidin [4]. The structural anthocyanin biosynthesis genes of many plant species have been identified and are mostly polygenes, such as PAL, CHI, CHS, and DFR [5,6,7].

These structural genes are simultaneously regulated by an activated transcriptional MBW complex, consisting of the R2R3-MYB and bHLH transcription factors, and the WD40 repeat protein [8, 9]. The MYB family is one of the largest gene families in plant. It is broadly divided into four sectors according to the number of conserved domains [10]. The R2R3-MYB transcription factors play a central role in anthocyanin biosynthesis through their influence on the transcriptional level of the structural anthocyanin biosynthesis genes [11, 12]. Most are positive regulators [e.g., AtMYB75, PpMYB9 and NnMYB5 [13,14,15]. However, others are inhibitory in the context of anthocyanin synthesis [e.g., EsMYBF1, PtMYB182 and VvMYB4 [16,17,18]. BHLHs belong to a multi-gene family that can be divided into 26 subgroups, and flavonoid related bHLHs were grouped into subgroup IIIf [19]. In Arabidopsis, some bHLHs, including TT8, GL3 and EGL3, have been identified as participants in the biosynthesis of different flavonoids [20,21,22]. In addition, the bHLHs which were associated with anthocyanin biosynthesis have been identified from a variety of plant species, such as CmbHLH2 [23], LcbHLH3 [24], MtbHLH (MtTT8) [25], LebHLH (AH) [26] and VvMYC1 [27]. Some other transcription factors, including Elongated Hypocotyl 5 (HY5), light-regulated zinc finger protein (LZF), CONSTANS-like (COL) and squamosa promoter binding protein-like 9 (SPL9), have also been shown to be associated with anthocyanin biosynthesis [28,29,30,31].

Light is a key environmental factor that affects the synthesis of anthocyanins [5, 32,33,34,35,36]. This process is facilitated by special photoreceptors. They transmit a signal to the downstream transcription factors that regulate plant growth and metabolic processes, including anthocyanin biosynthesis. The photoreceptors in higher plants have been identified to date including phytochromes (PHYs), phototropins (PHOTs), cryptochromes (CRYs) and UV-resistance locus 8 (UVR8). The light receptors are sensitive to specific wavelengths: PHYs detect red/far red, CRYs detect blue [37], PHOTs detect UV-A, and UVR8 detects the ultraviolet wavelengths [38]. The significance of the photoreceptors in anthocyanin biosynthesis has been confirmed by experiments with Arabidopsis. For example, phyA mutants could not induce the expression of CHS (an enzyme in the anthocyanin biosynthesis pathway) and the synthesis of anthocyanins under far red light [39]. The expression level of CHS and the content of flavonoids [40] were significantly reduced in the uvr8 mutants. In addition, CRY1 induced anthocyanins synthesis under blue light through the activation of anthocyanin biosynthesis enzymes [41].

Downstream of these photoreceptors, Constitutive Photomorphogenic 1 (COP1) acts as a central repressor in the anthocyanin biosynthesis pathway, including PAP1, PAP2, HY5, Hypocotyl Homolog (HYH), Long Hypocotyl in Far-Red 1 (HFR1) and Long After Far-Red Light 1 (LAF1) [42,43,44,45,46,47]. COP1 is considered to be a ‘molecular switch’ in metabolic processes which are stimulated by light [48]. Under dark conditions, COP1 was localized in the nucleus. COP1 promotes the degradation of the photomorphogenesis-promoting transcription factors [49]. Under light, the concentration of COP1 in the nucleus decreased rapidly [50] and the light-activated photoreceptors inhibited COP1 activities, thereby causing accumulation of the photomorphogenesis-promoting transcription factors. Maier et al. [46] reported that large amounts of anthocyanins were synthesized in cop1 mutants under dark conditions. In addition, COP1 formed a complex with phytochrome A-105 (SPA) in the nucleus to inhibit photomorphogenesis [46, 51, 52]. In Arabidopsis, spa mutants exhibited similar features with COP1 mutants under dark conditions [53, 54].

Eggplant (Solanum melongena L.) is a commercially significant crop cultivated and consumed in many countries [55]. The varieties of eggplant with dark purple peel, which is caused by a high anthocyanin concentration, are more attractive to consumers than the paler types. The anthocyanin concentration in the purple eggplant cultivars is high in comparison to that of other deeply colored fruits and vegetables, e.g., 2.34× that of grapes, and 7.08× that of red onions [56].

The observation that anthocyanin biosynthesis is stimulated by light is paralleled in other species, e.g. lithci, Chinese bayberry, grape berry, apple, and red pear [6, 36, 47, 57]. Previous researches focused on exploring the regulatory mechanisms associated with the accumulation of anthocyanins at the molecular level [6, 7]. Recently, the effect of light on the whole process of anthocyanins biosynthesis was preliminary studied in our previous study [35]. However, the underlying molecular mechanisms of light-induced anthocyanins biosynthesis before color becomes evident are not well known.

Our previous study found out that anthocyanin biosynthesis in the peel of eggplant (cv Lanshan Hexian) is regulated by light and color becomes evident at most 2 days after exposure to light. [5, 35]. In this study, the expression patterns of a variety of anthocyanin biosynthesis structural genes and regulator genes under light over an 8 h period were analyzed by qRT-PCR. In the present study, RNA-seq was used to capture the majority of transcriptional changes in anthocyanin biosynthesis related genes. This analysis was attempted to elucidate the molecular mechanisms of light-induced anthocyanin biosynthesis before color becomes evident. This study will facilitate genetic engineering protocols which would improve the anthocyanin content of eggplant or other plants under low-light condition.

Methods

Plant materials and treatments

The eggplants were grown in the horticultural farm of Shanghai Jiao Tong University, Shanghai China. After full bloom, the sepals were covered with paper bags which were completely light-impermeable. The bags were removed from the fully-grown fruits after 24 days at 8:00 am. The peels of the eggplants were collected immediately after bags removal, and then after a further 0.5, 4 and 8 h. Three eggplant fruits were mixed as a sample and three samples were set as biological replicates. The pooled samples were immediately frozen in liquid nitrogen, and stored at − 80 °C prior to further processing.

RNA extraction, library construction and transcriptome sequencing

Total RNA was extracted from pooled samples using RNAiso Plus (TaKaRa, Otsu, Shiga, Japan) according to the manufacturer’s instructions. After the total RNA extraction and DNase I treatment, the concentration and quality of each sample was examined using a NANODROP2000 (Thermo Scientific, Wilmington, DE) and Agilent 2100 Bioanalyzer. The mRNA was isolated by magnetic beads with Oligo (dT) and then synthesized to cDNA. Short fragments were purified and resolved with EB buffer for end reparation and single nucleotide A (adenine) addition. After that, the short fragments were connected with adapters. During the quality control (QC) steps, Agilent 2100 Bioanaylzer and ABI StepOnePlus Real-Time PCR System were used in quantification and qualification of the sample libraries. Finally, the libraries were sequenced using Illumina HiSeq™ 2000. Four time points (0, 0.5, 4 and 8 h) were selected for RNA-Seq analysis to obtain a general overview transcriptome of the eggplant peel samples that had been subjected to different light exposures. Three biological replicates were used in the RNA-Seq experiments involving all the four time points, and it means that a total of 12 samples were harvested to construct libraries. Pearson’s correlation was exploited to measure the relevance between each two samples.

Bioinformatics analysis

The raw reads produced by Illumina HiSeq™ 2000 were subjected to QC which can determine whether a resequencing step is needed or not. After read evaluated and filtered of low QC reads, reads were aligned to the eggplant genome sequences (http://eggplant.kazusa.or.jp/) [58]. And the alignment data was utilized to calculate distribution of reads on reference genes and mapping ratio. After alignment result passing QC, gene and isoform expression analysis was carried out. The differentially expressed genes (DEGs) were selected by a threshold of False Discovery Rate (FDR) ≤ 0.001 and an absolute Log2Ratio value ≥1 among the three biological replicates based on the analysis method of the poisson distribution. And the DEGs between groups were screened according to the following criteria: Fold change ≥2 and diverge probability ≥0.8 using Noiseq package method [59]. Further, the sequences of DEGs were compared with the NCBI non-redundant (Nr) database [60], Gene Ontology (GO) [61] and Kyoto Encyclopedia of Genes and Genomes (KEGG) [62] databases to identify and annotate the obtained DEGs using Blast software [63]. In addition, the Log2 (folds of mean RPKM values to the 0 h time point) were used to generate cluster diagrams by MultiExperiment Viewer software with a colour scale (MeV v4.4.1, http://www.tm4.org/) using the hierarchical clustering method [64].

RNA-seq data validation

qRT-PCR analysis was performed to validate the accuracy of the gene expression obtained from the assembled transcriptome via RNA-Seq. The correlation between the RNA-seq data and qRT-PCR result was analyzed by using SigmaStat 3.0 (SPSS, Chicago, IL). Total RNA were extracted from the peel samples by the MiniBEST Universal RNA Extraction Kit (TaKaRa) and examined using a NANODROP2000 (Thermo Scientific, Wilmington, DE) and characterized on 1.2% agarose gel electrophoresis. And the eligible RNA was selected based on a threshold of absorbance ratio 260/280 ≥ 1.8 and 260/230 ≥ 1.8. 500 ng RNA was synthesized into cDNA with the PrimeScript™ RT Master Mix (Perfect Real Time) (Takara). According to the manufacturer’s instructions of SYBR® Premix Ex Taq™ II (Tli RNaseH Plus) (Takara), qRT-PCR was performed on CFX Real Time PCR Detection System (BioRAD) using the following program: 95 °C for 30 s, followed by 40 cycles of 95 °C for 5 s and 60 °C for 30 s. The Actin gene (Sme2.5_00072.1_g00003.1) from eggplant was amplified in parallel as an internal reference gene [65]. The relative expression levels of the amplified products were analyzed using the 2-ΔΔCT method [66]. All qRT-PCR analyses were performed in three biological and three technical replications. The primers employed for this study are listed in Additional file 1.

Results

Anthocyanin biosynthesis structural genes were induced by light within an 8 h exposure period

The transcriptions of CRY1, CRY2, HY5 and MYB1 were previously shown to be triggered by exposure to light over an 8 h period [5]. Here, the expression levels of six anthocyanin biosynthesis structural genes (CHS, CHI, F3H, F3’5’H, DFR and ANS which were homologous with the genes reported in Arabidopsis) were analyzed by qRT-PCR (Fig. 1). The expression patterns of the six structural genes were similar to those of CRY1, CRY2, HY5 and MYB1 [5]. The qRT-PCR analysis showed that the expression levels of most of anthocyanin biosynthesis structural genes were up-regulated slightly at 0.5 h, peaked at 4 h and trended back at 8 h. And similar expression patterns were found in CRY1 and HY5 [5]. Thus the four time points were chosen to analyze the transcriptome of ‘Lanshan Hexian’ after bags removal using RNA-seq expecting to capture majority of transcriptional changes in anthocyanin biosynthesis related genes.

Fig. 1
figure 1

The expression patterns of 6 structural genes and correlation analysis of the RNA-seq data and qRT-PCR result. Correlation analysis showed correlation between RNA-seq data and qRT-PCR was good (Pearson R > 0.99). The left y-axis indicates the corresponding expression data of RNA-seq (black lines). The right y-axis shows the relative gene expression levels detected by qRT-PCR (red lines). The x-axis shows the length of time (h) the eggplant peel was exposed to light

Full size image

Transcriptional response of light by RNA-seq

An average of 49 million clean reads was obtained from each sample, of which 79~ 81% could be mapped to the eggplant genome (http://eggplant.kazusa.or.jp/): 46~ 50% could be mapped to eggplant genes (Additional file 2). Finally, the reads were assembled into a total of 32,629 genes or transcripts (Additional file 3). To test the samples chosen was reliable, the correlation value between each two samples was calculated based on FPKM (Fig. 2 and Additional file 4). A false discovery rate ≤ 0.001, an absolute Log2Ratio value ≥1 and diverge probability ≥0.8 were the thresholds applied to assess the significance of differences in gene expression. There were 1956 DEGs detected by means of time point to time point comparisons, including 843 (0.5 h vs 0 h), 948 (4 h vs 0 h), 776 (8 h vs 0 h), 865 (4 h vs 0.5 h), 941 (8 h vs 0.5 h) and 223 (8 h vs 4 h) (Fig. 3 and Additional file 5). All the known anthocyanin biosynthesis related genes were involved in the DEGs. The RNA-Seq and qRT-PCR results for these genes were highly correlated (r = 0.97~ 1.000) (Fig. 1). In addition, 10 selected genes, which had different expression patterns after bag removal, were chosen to validate the RNA-Seq results by qRT-PCR (r = 0.91~ 1.000) (Additional file 6). These results provided further validation of the reliability of the transcriptome data.

Fig. 2
figure 2

Heat map representation of the correlations coefficient between each two samples. The color represents the correlation coefficient (the darker the color, the higher the correlation, the lighter the color, the lower the correlation)

Full size image
Fig. 3
figure 3

Overview of differentially expressed genes (DEGs)

Full size image

GO analysis of the DEGs

The GO consortium provides a standard vocabulary that is dynamic and controlled, and that can be applied to all eukaryotes even though knowledge of the roles of genes and cell proteins is accumulating and expanding [67]. We therefore used GO assignments to classify the functions of the DEGs in eggplant peel. Of the total of 1956 DEGs, 873 have been annotated according to the GO database and are classified as ‘biological process’, ‘cellular component’ and ‘molecular function’. The results covered 34 important functional groups including 17 for biological process, 9 for cellular components, and 8 for molecular function (Fig. 4). The two largest subcategories in the ‘biological process’ category were ‘metabolic process’ including 500 DEGs and ‘cellular process’ including 344 DEGs. In the ‘cellular component’ category, the two largest subcategories were ‘cell’ and ‘cell part’ and both included 367 DEGs. For the ‘molecular function’ category, ‘catalytic activity’ including 514 DEGs and ‘binding activity’ including 429 DEGs were the most abundant.

Fig. 4
figure 4

GO categories assigned to DEGs. The left y-axis shows the percentages of proteins identified, and the right y-axis shows the protein number. The proteins were categorized according to the annotation of GO, and the number of each category is displayed based on biological process, cellular components, and molecular functions

Full size image

Cluster analysis of expression patterns in the DEGs

The 1956 DEGs were subjected to complete-linkage hierarchical clustering using a Euclidean distance metric by MeV v4.4.1 software and divided into 26 clusters (Additional file 7). We also present an overview of the KEGG pathway enrichment to provide a global description of the enriched biological pathways in each cluster of similarly regulated transcripts. The structural genes in anthocyanin biosynthesis pathway were divided into three pathways, including ‘phenylpropanoid biosynthesis’, ‘flavonoid biosynthesis,’ and ‘anthocyanin biosynthesis’. These pathways were enriched in clusters 7, 9, 12, 15, 18, 19 and 24 (Fig. 5). Therefore, genes in these clusters were the focus of further attention in this study. In addition, ‘metabolic pathways’ and ‘biosynthesis of secondary metabolites’ were highly enriched in six out of seven clusters.

Fig. 5
figure 5

Cluster analysis of DEGs with significant expression profile changes and KEGG pathway enrichment analysis. All the DEGs were subjected to complete-linkage hierarchical clustering using a Euclidean distance metric and divided into 26 clusters. The x-axis showed the time point to time point comparison. The y-axis shows the relative log2(ratio) of each comparison. Enriched KEGG pathways are listed to the right of each cluster

Full size image

As shown in Fig. 5, the expression level of genes in cluster 7 declined during this period. The expression level of genes in cluster 9 were down-regulated at 0.5 h but up-regulated to peak level in the samples taken after 4 h exposure to light, and maintained a high expression level. The genes in cluster 12, 15 and 19 were up-regulated slightly at 0.5 h, peaked after 4 h exposure to light, and maintained a high expression level at 8 h. Cluster 18 and 24 contained genes that were up-regulated to peak level in the samples taken after 0.5 h, but then the expression level of genes decreased after 4 and 8 h (Fig. 5). From the above, we found the genes in cluster 12, 15 and 19 have similar expression pattern but with different relative Log2ratio value at the same time point after eggplants exposure to light. Therefore, the genes in cluster 12, 15 and 19 could be analyzed together.

The structural anthocyanin biosynthesis genes

Twenty-four structural anthocyanin biosynthesis genes were identified. They participated in each step of anthocyanin metabolic pathway. 2 PAL, 4CL, CHS, 2 CHIs, F3H, 2 F3’Hs, F3’5’H, DFR, ANS, 2 3GT and 5GT presented in cluster 12, 15 and 19. The expression levels of these genes were increased and peaked after 4 h. In cluster 18, another three PAL genes were identified which expression level with peaked expression level at 0.5 h time point. The expression level of the exception, Sme2.5_03336.1_g00008.1 (PAL), another 4CL and 5GT were specifically induced after 0.5 h light exposure but then their activity declined. Another two CHS genes were down-regulated after 0.5 h, but then up-regulated were divided into cluster 9. The transcriptional level of another F3’H gene was down-regulated over the period after exposure to light (Fig. 6a).

Fig. 6
figure 6

Heat map representation of the genes expression patterns. a The structural anthocyanin biosynthesis genes, b All the transcription factors response to light, c The photoreceptors and genes related to light signal transduction. The color scale represents the Log2 (fold-change to the 0 h time point). Genes that are upregulated appear in yellow, and those that are downregulated appear in blue, with the relative Log2Ratio reflected by the intensity of the color

Full size image

Transcription factors

One hundred and thirty-eight transcription factors that exhibited highly dynamic changes in response to light were found, 102 of which were present in the 7 clusters. Among the 102 transcription factors, MYB (18 genes) constituted the largest group, followed by AP2-EREBP (16 genes), WRKY (15 genes), GRAS (6 genes), bHLH (6 genes), C2C2-CO-like (5 genes), C2C2-Dof (3 genes), HSF (3 genes), SPL (3 genes), zf-HD (2 genes), MADS (2 genes), mTERF (2 genes), TCP (2 genes), NAC (2 genes), C3H (2 genes), Tify (2 genes), C2C2-YABBY (2 genes), G2-like, Sigma70-like, C2C2-GATA, C2H2, Trihelix, HRT, HY5, bZIP, TAZ, TGA1and BES1 (Fig. 6b). Strikingly, HY5, MYB113, TT8 and WRKY44 (TTG2), which were known as anthocyanin biosynthesis regulated genes [25, 31, 46, 68, 69], were divided into clusters 12, 15 and 19 (Fig. 6b). In the three clusters, another 8 MYBs, 3 C2C2-CO-likes, 2 SPLs, bZIP, WRKY and bHLH were found. Previous studies have reported that genes which belong to the same gene families with these transcription factors participated in anthocyanin biosynthesis [24, 28, 31, 69,70,71,72,73]. Thus, we suggested that these transcription factors might be involved in light-induced anthocyanin biosynthesis.

Plant light signal perception and transduction

Under light condition, photoreceptors act as central promotors to substitute ubiquitination and degradation of the positive regulators in anthocyanin biosynthesis pathway caused by COP1, such as HY5 [42] and MYB1 [46, 47]. Here, UVR3, CRY3 and UVR8 were identified in clusters 9, 15 and 24, respectively (Fig. 6c). The expression levels of CRY3 and UVR3 were changed slightly at 0.5 h time point, increased to the highest level at 4 h time point but then declined slightly at 8 h time point. By contrast, the expression level of UVR8 was increased to a peak after 0.5 h exposure to light, declined towards the control level after 4 h and continued to decrease for the next 4 h. The COP1/SPA complex act as a molecular switch of light-induced anthocyanin biosynthesis located downstream of the photoreceptors. In this study, COP1 and two SPAs were found in clusters 12 and 15 (Fig. 6c). KEGG pathway analysis showed that the three genes were enriched in ‘ubiquitin mediated proteolysis’.

Discussion

In our previous study, anthocyanins couldn’t accumulate in eggplant peel under shaded conditions. After eggplants were exposed to light, anthocyanins rapidly synthesized within a short time [5, 35]. These findings indicated that the coloration of eggplants peel completely depend on light, as observed in previous studies on lithci, Chinese bayberry, grape berry, apple, and red pear [6, 36, 47, 57]. The existing studies focused on exploring the molecular regulatory mechanisms associated with the accumulation of anthocyanins, including the paper we recently published [6, 7, 35]. However, whether and how does light participate in regulating anthocyanin biosynthesis before color becomes evident have not been studied so far. Here, RNA-seq was used to explore the relationship between light and the expression of candidate genes, which might be involved in anthocyanin biosynthesis over an 8 h period in eggplant.

Jiang et al. has reported that the transcriptional levels of CRY1, CRY2, HY5 and MYB1 were significantly induced by exposure to light over an 8 h period [5]. In addition, qRT-PCR analysis showed that the expression levels of six anthocyanin biosynthesis structural genes were also affected in our pre-preparation work. Then another 18 structural anthocyanin biosynthesis genes were identified as DEGs by RNA-seq analysis and divided into 7 clusters. Most of them were up-regulated after exposure to light for 4 h. These results suggested that the light response of these genes is dependent upon some upstream transcription factors. For instance, HY5 is necessary for the rapid transcription of F3H and CHS during the dark-to-light transition, which eventually allows the accumulation of anthocyanins [74]. In addition, the expression patterns of these anthocyanin biosynthesis structural genes corresponded to the light intensity fluctuation which was parabolic in the day. This result might reveal the direct reason why eggplant peel shows poor color production under weak light conditions.

In this study, 102 transcription factors which were significantly regulated by light were involved in the 7 clusters. Among these transcription factors, the well-known anthocyanin biosynthesis regulating genes including HY5 [31], MYB113 [31], TT8 [9] and WRKY44 (TTG2) [69] were found. As is well known, HY5 [42, 51] and MYB [46] were repressed by the activity of COP1 in darkness, but they triggered expression of light-inducible genes, such as the gene encoding CHS after light exposure [75]. In addition, MYB1 [76, 77], MYB2 [12], MYB5 [70], COL11 (C2C2-CO-like) [28], MADS [78], MYC2 (bHLH) [79] and SPL [29] were reported to participate in anthocyanin biosynthesis in various plant species. In addition, the connections between other transcription factors and anthocyanin synthesis have attracted much attention. For example, the family genes of WRKY [69], C2C2-CO-like [28, 80] and NAC [81] were also reported to be involved in anthocyanin biosynthesis. In the present study, 15 WRKYs, 6 GRASs, 5 C2C2-CO-likes, 3 SPLs and NAC were screened in the seven clusters that included anthocyanin biosynthesis structural genes. These genes might be involved in the regulation of anthocyanin biosynthesis. Also, these genes point to the direction for further studies of anthocyanin biosynthesis in eggplants.

COP1 is an E3 ubiquitin ligase that represses photomorphogenesis by mediating light-dependent degradation of various transcription factors involved in light signaling. The subcellular localization of COP1 is light dependent [82]. In darkness, COP1 is situated in the nucleus where it interacts with photomorphogenesis-promoting transcription factors. When plants are exposed to light, the concentration of COP1 in the nucleus declines rapidly [50]. SmCOP1, which is annotated as Sme2.5_00128.1_g00013.1, is functionally complemented in phenotypes of corresponding Arabidopsis mutants where it interacts with SmHY5 and SmMYB1 (SmMYB113) [5]. In this study, another COP1 annotated as Sme2.5_00499.1_g00010.1, and two SPA genes were located in clusters 12 and 15. Their expression levels kept control level at 0.5 h time point, increased to the highest level at 4 h time point but then declined a little at 8 h time point. Strikingly, the expression level of Sme2.5_00128.1_g00013.1 showed little change after eggplants were exposed to light when compared to control level. Stacey et al. [82] reported that light is the control of COP1 partitioning between nucleus and cytoplasm. Maier and Hoecker [83] found that mutations in the COP1/SPA complex, which is the most important repressor of light signaling, cause a hyperaccumulation of anthocyanins in normal and high intensity light. Taking into account the previous information from Arabidopsis, we concluded that the expression level of COP1 may be not be regulated by light and the function of the COP1 protein was to combine the photoreceptors but not the positive regulators, in the context of anthocyanin biosynthesis. These results suggest that Sme2.5_00499.1_g00010.1 may participate in light-induced anthocyanin biosynthesis in a different way.

Three photoreceptors in eggplant, CRY3, UVR3, and UVR8, were identified as DEGs in this study. Jiang et al. [5] suggested that a blue-light-triggered CRY1/CRY2-COP1 interaction allows HY5 and MYB1 to combine with the downstream anthocyanin biosynthesis genes (CHS and DFR) in eggplant. These results from this later experiment suggest that the expression of CRY3 was stimulated by light, but CRY1 and CRY2 were not thus stimulated. Also, the expression pattern of CRY3 was similar to those of CRY1 and CRY2 [5]. The anthocyanin content and the expression level of UVR8 were also elevated under UV-B in peach, strawberry, lettuce and radish sprouts [84,85,86,87]. In this study, UVR8 was up-regulated by light after 0.5 h exposure and was down-regulated after 4 and 8 h. It is possible that unknown genes negatively regulate the expression of UVR8 under sunlight. In addition, UVR3 encodes a photolyase that is specific for 6,4 photoproducts in Arabidopsis thaliana but which is not activated by light [88]. In eggplant peel the UVR3 mRNA level was increased within 4 h after bags removal suggesting that UVR3 responds to light differently in eggplant.

Comparing the transcriptome data with our previous studies [35], we found that the 1956 DEGs could be divided into two parts: 869 genes could be identified as DEGs in both the two studies (Additional file 8) and the other 1087 genes were only identified as DEGs in this study. Then the biological processes and pathways involved in the 869 DEGs were identified by GO and KEGG pathway enrichment analysis (Additional files 9 and 10). Except for the ‘anatomical structure formation’ in the ‘biological process’ category was not found, the functional groups involved in the 869 DEGs were same with which involved in all the DEGs (Additional file 9 and Fig. 4). KEGG pathway enrichment analysis showed that ‘phenylpropanoid biosynthesis’, ‘flavonoid biosynthesis,’ and ‘anthocyanin biosynthesis’ pathways were enriched in the 869 DEGs (Additional file 10). Further analysis found that 22 structural genes, 37 transcription factors and 3 photoreceptors were included in the 869 DEGs. Strikingly, 3 AP2-EREBPs, 3 WRKYs, 2 bHLHs, 2 MYBs, 2 C2C2-CO-likes, C2C2-YABBY, GRAS, HSF, MADS, mTERF, SPL, TAZ and TGA1 were found in the 37 transcription factors except for HY5, MYB113, TT8 and WRKY44. These results suggested that these 869 DEGs participated in the light-induced anthocyanin biosynthesis. As for the identification of the other 1087 DEGs, this might be consequence of the light condition change, the circadian clock or other factors.

Conclusion

In this study, RNA-seq was applied to elucidate the underlying molecular mechanism of light-induced anthocyanin biosynthesis in eggplant peel before color becomes evident. A total of 1956 DEGs were identified and divided into 26 clusters based on the distinct expression patterns. By comparing the transcriptome data with our previous study [35], the structural genes, photoreceptors, light signal transduction elements and transcription factors involved in the anthocyanin biosynthesis pathway could be found in the intersection of the two transcriptome data. These results suggested that light induce anthocyanin biosynthesis from the moment of exposure to light to the saturation of anthocyanin. This study will expand our understanding of the mechanism of light-induced anthocyanin biosynthesis in plant, and point to new directions in the study of anthocyanin biosynthesis.

Change history

  • 04 May 2018

    After publication of the original article [1] it was noted that in Additional file 1: Table S1, and Fig. 1, specific primer sequences were incorrect, and taken from Sme2.5_02154.1_g00001.1 rather than Sme2.5_13923.1_g00001.1.

Abbreviations

4CL:

4-coumarate-CoA ligase

ANS:

Anthocyanin synthase

AT:

Acylation

CHI:

Chalcone isomerase

CHS:

Chalcone synthase

COL:

CONSTANS-like

COP1:

Constitutive Photomorphogenic 1

CRY3:

Cryptochrome 3

DEGs:

Differentially expressed genes

DFR:

Dihydroflavonol reductase

F3H:

Flavanone 3-hydroxyl enzyme

FDR:

False discovery rate

GO:

Gene ontology

GT:

Glycosylation

HFR:

Hypocotyl in Far-Red 1

HY5:

Hypocotyl 5

HYH:

Hypocotyl homolog

KEGG:

Kyoto encyclopedia of genes and genomes

LAF1:

Long After Far-Red Light 1

LDOX:

Leucoanthocyanidin dioxygenase

LZF:

Light-regulated zinc finger protein

MT:

Methylation

Nr:

NCBI non-redundant database

PAL:

Phenylalanine ammonia lyase

PHOTs:

Phototropins

PHYs:

Phytochromes

QC:

Quality control

SPL9:

Squamosa promoter binding protein-like 9

SRA:

Sequence read archive

UVR8:

UV Resistance Locus 8

References

  1. Stintzing FC, Carle R. Functional properties of anthocyanins and betalains in plants, food, and in human nutrition. Trends Food Sci Technol. 2004;15:19–38.

    Article  CAS  Google Scholar 

  2. Sakuta M. Diversity in plant red pigments: anthocyanins and betacyanins. Plant Biotechnol Rep. 2014;8:37–48.

    Article  Google Scholar 

  3. Winkel-Shirley B. Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. Plant Physiol. 2001;126(February 2015):485–93.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  4. Nishihara M, Nakatsuka T, Yamamura S. Flavonoid components and flower color change in transgenic tobacco plants by suppression of chalcone isomerase gene. FEBS Lett. 2005;579(27):6074–8.

    Article  PubMed  CAS  Google Scholar 

  5. Jiang M, Ren L, Lian H, Liu Y, Chen H. Novel insight into the mechanism underlying light-controlled anthocyanin accumulation in eggplant (Solanum melongena L.). Plant Sci. 2016;249:46–58.

    Article  PubMed  CAS  Google Scholar 

  6. Zhang HN, Li WC, Wang HC, Shi SY, Shu B, Liu LQ, Wei YZ, Xie JH. Transcriptome profiling of light-regulated anthocyanin biosynthesis in the pericarp of Litchi. Front Plant Sci. 2016;7:963.

    PubMed  PubMed Central  Google Scholar 

  7. Wang Z, Du H, Zhai R, Song L, Ma F, Xu L. Transcriptome analysis reveals candidate genes related to color fading of ‘Red Bartlett’ (Pyrus communis L.). Front Plant Sci. 2017;8:455.

    PubMed  PubMed Central  Google Scholar 

  8. Ramsay NA, Glover BJ. MYB-bHLH-WD40 protein complex and the evolution of cellular diversity. Trends Plant Sci. 2005;10:63–70.

    Article  PubMed  CAS  Google Scholar 

  9. Baudry A, Heim MA, Dubreucq B, Caboche M, Weisshaar B, Lepiniec L. TT2, TT8, and TTG1 synergistically specify the expression of BANYULS and proanthocyanidin biosynthesis in Arabidopsis thaliana. Plant J. 2004;39(3):366–80.

    Article  PubMed  CAS  Google Scholar 

  10. 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.

    Article  PubMed  CAS  Google Scholar 

  11. Liu JY, Osbourn A, Ma PD. MYB transcription factors as regulators of Phenylpropanoid metabolism in plants. Mol Plant. 2015;8(5):689–708.

    Article  PubMed  CAS  Google Scholar 

  12. Chiu LW, Li L. Characterization of the regulatory network of BoMYB2 in controlling anthocyanin biosynthesis in purple cauliflower. Planta. 2012;236(4):1153–64.

    Article  PubMed  CAS  Google Scholar 

  13. Sun S-S, Gugger PF, Wang Q-F, Chen J-M. Identification of a R2R3-MYBgene regulating anthocyanin biosynthesis and relationships between its variation and flower color difference in lotus ( Nelumbo Adans.). PeerJ. 2016;4:e2369.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Li S, Wang W, Gao J, Yin K, Wang R, Wang C, Petersen M, Mundy J, Qiu J-L. MYB75 phosphorylation by MPK4 is required for light-induced anthocyanin accumulation in Arabidopsis. Plant Cell. 2016;28(11):2866–83.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Zhou H, Peng Q, Zhao J, Owiti A, Ren F, Liao L, Wang L, Deng X, Jiang Q, Han Y. Multiple R2R3-MYB transcription factors involved in the regulation of anthocyanin accumulation in peach flower. Front Plant Sci. 2016;7(October):1–11.

    Google Scholar 

  16. Pérez-Díaz JR, Pérez-Díaz J, Madrid-Espinoza J, González-Villanueva E, Moreno Y, Ruiz-Lara S. New member of the R2R3-MYB transcription factors family in grapevine suppresses the anthocyanin accumulation in the flowers of transgenic tobacco. Plant Mol Biol. 2016;90(1–2):63–76.

    Article  PubMed  CAS  Google Scholar 

  17. Yoshida K, Ma D, Constabel CP. The MYB182 protein down-regulates Proanthocyanidin and anthocyanin biosynthesis in poplar by repressing both structural and regulatory flavonoid genes. Plant Physiol. 2015;167(3):693–710.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Huang W, Khaldun ABM, Chen J, Zhang C, Lv H, Yuan L, Wang Y. A R2R3-MYB transcription factor regulates the Flavonol biosynthetic pathway in a traditional Chinese medicinal plant, Epimedium sagittatum. Front Plant Sci. 2016;7(July):1–12.

    Google Scholar 

  19. Pires N, Dolan L. Origin and diversification of basic-helix-loop-helix proteins in plants. Mol Biol Evol. 2010;27(4):862–74.

    Article  PubMed  Google Scholar 

  20. Zhang F, Gonzalez A, Zhao M, Payne CT, Lloyd A. A network of redundant bHLH proteins functions in all TTG1-dependent pathways of Arabidopsis. Development. 2003;130(20):4859–69.

    Article  PubMed  CAS  Google Scholar 

  21. Payne CT, Zhang F, Lloyd AM. GL3 encodes a bHLH protein that regulates trichome development in arabidopsis through interaction with GL1 and TTG1. Genetics. 2000;156(3):1349–62.

    PubMed  PubMed Central  CAS  Google Scholar 

  22. Shirley BW, Kubasek WL, Storz G, Bruggemann E, Koornneef M, Ausubel FM, Goodman HM. Analysis of Arabidopsis mutants deficient in flavonoid biosynthesis. Plant J. 1995;8(5):659–71.

    Article  PubMed  CAS  Google Scholar 

  23. Xiang LL, Liu XF, Li X, Yin XR, Grierson D, Li F, Chen KS. A novel bHLH transcription factor involved in regulating anthocyanin biosynthesis in chrysanthemums (chrysanthemum morifolium ramat). PLoS One. 2015;10(11):1–17.

    Article  CAS  Google Scholar 

  24. Lai B, Du L-N, Liu R, Hu B, Su W-N, Qin Y-H, Zhao J-T, Wang H-C, Hu G-B. Two LcbHLH transcription factors interacting with LcMYB1 in regulating late structural genes of anthocyanin biosynthesis in Nicotiana and Litchi chinensis during anthocyanin accumulation. Front Plant Sci. 2016;7(166):1–15.

    Google Scholar 

  25. Li P, Chen B, Zhang G, Chen L, Dong Q, Wen J, Mysore KS, Zhao J. Regulation of anthocyanin and proanthocyanidin biosynthesis by Medicago truncatula bHLH transcription factor MtTT8. New Phytol. 2016;210(3):905–21.

    Article  PubMed  CAS  Google Scholar 

  26. Qiu Z, Wang X, Gao J, Guo Y, Huang Z, Du Y. The tomato Hoffman's anthocyaninless gene encodes a bHLH transcription factor involved in anthocyanin biosynthesis that is developmentally regulated and induced by low temperatures. PLoS One. 2016;11(3):1–22.

    Google Scholar 

  27. Hichri I, Heppel SC, Jrm P, Léon C, Czemmel S, Delrot S, Lauvergeat V, Bogs J. The basic helix-loop-helix transcription factor MYC1 is involved in the regulation of the flavonoid biosynthesis pathway in grapevine. Mol Plant. 2010;3(3):509–23.

    Article  PubMed  CAS  Google Scholar 

  28. Bai S, Saito T, Honda C, Hatsuyama Y, Ito A, Moriguchi T. An apple B-box protein, MdCOL11, is involved in UV-B- and temperature-induced anthocyanin biosynthesis. Planta. 2014;240(5):1051–62.

    Article  PubMed  CAS  Google Scholar 

  29. Gou J-Y, Felippes FF, Liu C-J, Weigel D, Wang J-W. Negative regulation of anthocyanin biosynthesis in Arabidopsis by a miR156-targeted SPL transcription factor. Plant Cell. 2011;23(4):1512–22.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Chang CSJ, Li YH, Chen LT, Chen WC, Hsieh WP, Shin J, Jane WN, Chou SJ, Choi G, Hu JM, et al. LZF1, a HY5-regulated transcriptional factor, functions in Arabidopsis de-etiolation. Plant J. 2008;54(2):205–19.

    Article  PubMed  CAS  Google Scholar 

  31. Shin DH, Choi M, Kim K, Bang G, Cho M, Choi SB, Choi G, Park YI. HY5 regulates anthocyanin biosynthesis by inducing the transcriptional activation of the MYB75/PAP1 transcription factor in Arabidopsis. FEBS Lett. 2013;587(10):1543–7.

    Article  PubMed  CAS  Google Scholar 

  32. Takos AM, Jaffe FW, Jacob SR, Bogs J, Robinson SP, Walker AR. Light-induced expression of a MYB gene regulates anthocyanin biosynthesis in red apples. Plant Physiol. 2006;142(3):1216–32.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Feng FJ, Li MJ, Ma FW, Cheng LL. The effects of bagging and debagging on external fruit quality, metabolites, and the expression of anthocyanin biosynthetic genes in ‘Jonagold’ apple (Malus domestica Borkh.). Sci Hortic. 2014;165:123–31.

    Article  CAS  Google Scholar 

  34. Zoratti L, Karppinen K, Escobar AL, Haggman H, Jaakola L. Light-controlled flavonoid biosynthesis in fruits. Front Plant Sci. 2014;5:534.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Li J, Ren L, Gao Z, Jiang M, Liu Y, Zhou L, He Y, Chen H. Combined transcriptomic and proteomic analysis constructs a new model for light-induced anthocyanin biosynthesis in eggplant (Solanum melongena L.). Plant Cell Environ. 2017;40(12):3069–87.

    Article  PubMed  CAS  Google Scholar 

  36. 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.

    Article  PubMed  CAS  Google Scholar 

  37. Lin C. Plant blue-light receptors. Trends Plant Sci. 2000;5:337–42.

    Article  PubMed  CAS  Google Scholar 

  38. Rizzini L, Favory JJ, Cloix C, Faggionato D, O'Hara A, Kaiserli E, Baumeister R, Schafer E, Nagy F, Jenkins GI, et al. Perception of UV-B by the Arabidopsis UVR8 protein. Science. 2011;332(6025):103–6.

    Article  PubMed  CAS  Google Scholar 

  39. Kunkel T, Neuhaus G, Batschauer A, Chua NH, Schäfer E. Functional analysis of yeast-derived phytochrome A and B phycocyanobilin adducts. Plant J. 1996;10:625–36.

    Article  PubMed  CAS  Google Scholar 

  40. Kliebenstein DJ, Lim JE, Landry LG, Last RL. Arabidopsis UVR8 regulates ultraviolet-B signal transduction and tolerance and contains sequence similarity to human regulator of chromatin condensation 1. Plant Physiol. 2002;130(1):234–43.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  41. Ahmad M, Lin C, Cashmore AR. Mutations throughout an Arabidopsis blue-light photoreceptor impair blue-light-responsive anthocyanin accumulation and inhibition of hypocotyl elongation. Plant J. 1995;8(5):653–8.

    Article  PubMed  CAS  Google Scholar 

  42. Osterlund MT, Hardtke CS, Wei N, Deng XW. Targeted destabilization of HY5 during light-regulated development of Arabidopsis. Nature. 2000;405(6785):462–6.

    Article  PubMed  CAS  Google Scholar 

  43. Shalitin D, Yang H, Mockler TC, Maymon M, Guo H, Whitelam GC, Lin C. Regulation of Arabidopsis cryptochrome 2 by blue-light-dependent phosphorylation. Nature. 2002;417(6890):763–7.

    Article  PubMed  CAS  Google Scholar 

  44. Seo HS, Yang JY, Ishikawa M, Bolle C, Ballesteros ML, Chua NH. LAF1 ubiquitination by COP1 controls photomorphogenesis and is stimulated by SPA1. Nature. 2003;423(6943):995–9.

    Article  PubMed  CAS  Google Scholar 

  45. Yang JP, Lin RC, James S, Hoecker U, Liu BL, Xu L, Deng XW, Wang HY. Light regulates COP1-mediated degradation of HFR1, a transcription factor essential for light signaling in arabidopsis. Plant Cell. 2005;17(3):804–21.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  46. Maier A, Schrader A, Kokkelink L, Falke C, Welter B, Iniesto E, Rubio V, Uhrig JF, Hulskamp M, Hoecker U. Light and the E3 ubiquitin ligase COP1/SPA control the protein stability of the MYB transcription factors PAP1 and PAP2 involved in anthocyanin accumulation in Arabidopsis. Plant J. 2013;74(4):638–51.

    Article  PubMed  CAS  Google Scholar 

  47. Li YY, Mao K, Zhao C, Zhao XY, Zhang HL, Shu HR, Hao YJ. MdCOP1 ubiquitin E3 ligases interact with MdMYB1 to regulate light-induced anthocyanin biosynthesis and red fruit coloration in apple. Plant Physiol. 2012;160(2):1011–22.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  48. Yin R, Skvortsova MY, Loubéry S, Ulm R. COP1 is required for UV-B–induced nuclear accumulation of the UVR8 photoreceptor. Proc Natl Acad Sci. 2016;113(30):201607074.

    Article  CAS  Google Scholar 

  49. Datta S, Johansson H, Hettiarachchi C, Irigoyen ML, Desai M, Rubio V, Holm M. LZF1/SALT TOLERANCE HOMOLOG3, an Arabidopsis B-box protein involved in light-dependent development and gene expression, undergoes COP1-mediated ubiquitination. Plant Cell. 2008;20(9):2324–38.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. Hong SH, Kim HJ, Ryu JS, Choi H, Jeong S, Shin J, Choi G, Nam HG. CRY1 inhibits COP1-mediated degradation of BIT1, a MYB transcription factor, to activate blue light-dependent gene expression in Arabidopsis. Plant J. 2008;55(3):361–71.

    Article  PubMed  CAS  Google Scholar 

  51. Saijo Y, Sullivan JA, Wang H, Yang J, Shen Y, Rubio V, Ma L, Hoecker U, Deng XW. The COP1-SPA1 interaction defines a critical step in phytochrome A-mediated regulation of HY5 activity. Genes Dev. 2003;17(21):2642–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  52. Liu B, Zuo Z, Liu H, Liu X, Lin C. Arabidopsis cryptochrome 1 interacts with SPA1 to suppress COP1 activity in response to blue light. Genes Dev. 2011;25(10):1029–34.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Hoecker U, Xu Y, Quail PH. SPA1: a new genetic locus involved in phytochrome A-specific signal transduction. Plant Cell. 1998;10(January):19–33.

    PubMed  PubMed Central  CAS  Google Scholar 

  54. Laubinger S, Fittinghoff K, Hoecker U. The SPA quartet: a family of WD-repeat proteins with a central role in suppression of photomorphogenesis in arabidopsis. Plant Cell. 2004;16(9):2293–306.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  55. Cericola F, Portis E, Lanteri S, Toppino L, Barchi L, Acciarri N, Pulcini L, Sala T, Rotino GL. Linkage disequilibrium and genome-wide association analysis for anthocyanin pigmentation and fruit color in eggplant. BMC Genomics. 2014;15:896.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Wu X, Beecher GR, Holden JM, Haytowitz DB, Gebhardt SE, Prior RL. Concentrations of anthocyanins in common foods in the United States and estimation of normal consumption. J Agric Food Chem. 2006;54(11):4069–75.

    Article  PubMed  CAS  Google Scholar 

  57. Niu SS, Xu CJ, Zhang WS, Zhang B, Li X, Lin-Wang K, Ferguson IB, Allan AC, Chen KS. Coordinated regulation of anthocyanin biosynthesis in Chinese bayberry (Myrica rubra) fruit by a R2R3 MYB transcription factor. Planta. 2010;231(4):887–99.

    Article  PubMed  CAS  Google Scholar 

  58. Hirakawa H, Shirasawa K, Miyatake K, Nunome T, Negoro S, Ohyama A, Yamaguchi H, Sato S, Isobe S, Tabata S, et al. Draft genome sequence of eggplant (Solanum melongena L.): the representative solanum species indigenous to the old world. DNA Res. 2014;21(6):649–60.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  59. Tarazona S, Furio-Tari P, Turra D, di Pietro A, Nueda MJ, Ferrer A, Conesa A. Data quality aware analysis of differential expression in RNA-seq with NOISeq R/Bioc package. Nucleic Acids Res. 2015;43(21):e140.

    PubMed  PubMed Central  Google Scholar 

  60. Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35(Database):D61–5.

    Article  PubMed  CAS  Google Scholar 

  61. Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L, et al. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006;34(Web Server):W293–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  62. Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004;32(Database issue):D277–80.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  63. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  64. Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, et al. TM4: a free, open-source system for microarray data management and analysis. BioTechniques. 2003;34(2):374–8.

    Article  PubMed  CAS  Google Scholar 

  65. Liu H, Fu D, Zhu B, Yan H, Shen X, Zuo J, Zhu Y, Luo Y. Virus-induced gene silencing in eggplant (Solanum melongena). J Integr Plant Biol. 2012;54(6):422–9.

    Article  PubMed  CAS  Google Scholar 

  66. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001;25(4):402–8.

    Article  PubMed  CAS  Google Scholar 

  67. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat Genet. 2000;25(1):25–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  68. Nesi N, Debeaujon I, Jond C, Pelletier G, Caboche M, Lepiniec L. The TT8 gene encodes a basic helix-loop-helix domain protein required for expression of DFR and BAN genes in Arabidopsis siliques. Plant Cell. 2000;12(10):1863–78.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  69. Gonzalez A, Brown M, Hatlestad G, Akhavan N, Smith T, Hembd A, Moore J, Montes D, Mosley T, Resendez J, et al. TTG2 controls the developmental regulation of seed coat tannins in Arabidopsis by regulating vacuolar transport steps in the proanthocyanidin pathway. Dev Biol. 2016;419(1):54–63.

    Article  PubMed  CAS  Google Scholar 

  70. Deluc L, Barrieu F, Marchive C, Lauvergeat V, Decendit A, Richard T, Carde JP, Merillon JM, Hamdi S. Characterization of a grapevine R2R3-MYB transcription factor that regulates the phenylpropanoid pathway. Plant Physiol. 2006;140(2):499–511.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  71. Cui LG, Shan JX, Shi M, Gao JP, Lin HX. The miR156-SPL9-DFR pathway coordinates the relationship between development and abiotic stress tolerance in plants (vol 80, pg 1108, 2014). Plant J. 2015;82(5):901.

    Article  CAS  Google Scholar 

  72. Xu WJ, Dubos C, Lepiniec L. Transcriptional control of flavonoid biosynthesis by MYB-bHLH-WDR complexes. Trends Plant Sci. 2015;20(3):176–85.

    Article  PubMed  CAS  Google Scholar 

  73. Schwinn KE, Boase MR, Bradley JM, Lewis DH, Deroles SC, Martin CR, Davies KM. MYB and bHLH transcription factor transgenes increase anthocyanin pigmentation in petunia and lisianthus plants, and the petunia phenotypes are strongly enhanced under field conditions. Front Plant Sci. 2014;5(November):1–11.

    Google Scholar 

  74. Lee J, He K, Stolc V, Lee H, Figueroa P, Gao Y, Tongprasit W, Zhao H, Lee I, Deng XW. Analysis of transcription factor HY5 genomic binding sites revealed its hierarchical role in light regulation of development. Plant Cell. 2007;19(3):731–49.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  75. Chattopadhyay S, Ang LH, Puente P, Deng XW, Wei N. Arabidopsis bZIP protein HY5 directly interacts with light-responsive promoters in mediating light control of gene expression. Plant Cell. 1998;10(5):673–83.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  76. Lai B, Li XJ, Hu B, Qin YH, Huang XM, Wang HC, Hu GB. LcMYB1 is a key determinant of differential anthocyanin accumulation among genotypes, tissues, developmental phases and ABA and light stimuli in Litchi chinensis. PLoS One. 2014;9(1):e86293.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  77. Aharoni A, de Vos CH, Wein M, Sun Z, Greco R, Kroon A, Mol JN, O'Connell AP. The strawberry FaMYB1 transcription factor suppresses anthocyanin and flavonol accumulation in transgenic tobacco. Plant J. 2001;28(3):319–32.

    Article  PubMed  CAS  Google Scholar 

  78. Jaakola L, Poole M, Jones MO, Kamarainen-Karppinen T, Koskimaki JJ, Hohtola A, Haggman H, Fraser PD, Manning K, King GJ, et al. A SQUAMOSA MADS box gene involved in the regulation of anthocyanin accumulation in bilberry fruits. Plant Physiol. 2010;153(4):1619–29.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  79. An JP, Li HH, Song LQ, Su L, Liu X, You CX, Wang XF, Hao YJ. The molecular cloning and functional characterization of MdMYC2, a bHLH transcription factor in apple. Plant Physiol Biochem. 2016;108:24–31.

    Article  PubMed  CAS  Google Scholar 

  80. Chang CS, Maloof JN, Wu SH. COP1-mediated degradation of BBX22/LZF1 optimizes seedling development in Arabidopsis. Plant Physiol. 2011;156(1):228–39.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  81. Zhou H, Lin-Wang K, Wang H, Gu C, Dare AP, Espley RV, He H, Allan AC, Han Y. Molecular genetics of blood-fleshed peach reveals activation of anthocyanin biosynthesis by NAC transcription factors. Plant J. 2015;82(1):105–21.

    Article  PubMed  CAS  Google Scholar 

  82. Stacey MG, Hicks SN, von Arnim AG. Discrete domains mediate the light-responsive nuclear and cytoplasmic localization of Arabidopsis COP1. Plant Cell. 1999;11(3):349–64.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  83. Maier A, Hoecker U. COP1/SPA ubiquitin ligase complexes repress anthocyanin accumulation under low light and high light conditions. Plant Signal Behav. 2015;10(1):e970440.

    Article  PubMed  CAS  Google Scholar 

  84. Sgherri C, Scattino C, Pinzino C, Tonutti P, Ranieri AM. Ultraviolet-B radiation applied to detached peach fruit: a study of free radical generation by EPR spin trapping. Plant Physiol Biochem. 2015;96:124–31.

    Article  PubMed  CAS  Google Scholar 

  85. Josuttis M, Dietrich H, Treutter D, Will F, Linnemannstons L, Kruger E. Solar UVB response of bioactives in strawberry (Fragaria x ananassa Duch. L.): a comparison of protected and open-field cultivation. J Agric Food Chem. 2010;58(24):12692–702.

    Article  PubMed  CAS  Google Scholar 

  86. Tsormpatsidis E, Henbest R, Davis FJ, Battey N, Hadley P, Wagstaffe A. UV irradiance as a major influence on growth, development and secondary products of commercial importance in Lollo Rosso lettuce ‘Revolution’grown under polyethylene films. Environ Exp Bot. 2008;63(1):232–9.

    Article  CAS  Google Scholar 

  87. Wu Q, Su N, Zhang X, Liu Y, Cui J, Liang Y. Hydrogen peroxide, nitric oxide and UV RESISTANCE LOCUS8 interact to mediate UV-B-induced anthocyanin biosynthesis in radish sprouts. Sci Rep. 2016;6:29164.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  88. Castells E, Molinier J, Drevensek S, Genschik P, Barneche F, Bowler C. det1-1-induced UV-C hyposensitivity through UVR3 and PHR1 photolyase gene over-expression. Plant J. 2010;63(3):392–404.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

The project was supported by the National Natural Science Foundation of China (31471870). The authors declare that they have no competing interests.

Funding

The financial of this research was supported by National Natural Science Foundation of China (31471870). Funding body had no role in the design of the study and collection, analysis, interpretation of data and in writing the manuscript.

Availability of data and materials

Full data accompanying these experiments has been deposited in sequence read archive (SRA)(https://www.ncbi.nlm.nih.gov/sra), the accessions for our submission are: The raw data for the DGE analysis were also deposited in the NCBI SRA under accession numbers SRR5650714 (including 0 h-1, 0 h-2 and 0 h-3), SRR5651526 (including 0.5 h-1, 0.5 h-2 and 0.5 h-3), SRR5658205 (including 4 h-1, 4 h-2 and 4 h-3), and SRR5658226 (including 8 h-1, 8 h-2 and 8 h-3). The login account was provided as following: username: junge; password: 15150531037hyj.

Author information

Authors and Affiliations

  1. School of Agriculture and Biology, Shanghai Jiao Tong University, 800 Dongchuan Road, Minhang District, Shanghai, 200240, China

    Jing Li, Yong-Jun He, Lu Zhou, Yang Liu, Mingmin Jiang, Li Ren & Huoying Chen

  2. Institute for Agri-Food Standards and Testing Technology, Shanghai Academy of Agricultural Sciences, 1000 Jinqi Road, Fengxian District, Shanghai, 201403, China

    Li Ren

Authors
  1. Jing Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yong-Jun He
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Lu Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Yang Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Mingmin Jiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Li Ren
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Huoying Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

H-YC, YL, JL, M-MJ and LR planned and designed the research. JL, Y-JH, LZ and LR performed experiments. JL analyzed data. JL and Y-JH wrote the manuscript. All authors read and approved the manuscript.

Corresponding author

Correspondence to Huoying Chen.

Ethics declarations

Ethics approval

Not applicable, as this study did not involve human or animal subjects, and the seeds of eggplant cultivar ‘Lanshan Hexian’ were stored in School of Agriculture and Biology, Shanghai Jiao Tong University.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. List of primers sequences used in this study. (DOCX 33 kb)

Additional file 2:

Table S2. Summary statistics of sequencing and assembly. (DOCX 28 kb)

Additional file 3:

Table S3. All the assembled genes in this project. Gene_id: gene ID number; FPKM: FPKM value of this gene; the differentially expressed genes (DEGs) were annotated by overstriking. (XLS 5248 kb)

Additional file 4:

Table S4. Correlation value between each two samples. (XLS 25 kb)

Additional file 5;

Table S5. All the identified DEGs in this study. Gene_id: gene ID number; TAIR10: gene ID in Arabidosis corresponding to the degs in eggplant; E-value: The lowest E-value of Blast-nr for each genes; FPKM: FPKM value of this gene; Means-0 h:mean expression (FPKM) of 0 h group; log2Ratio(s2/s1): Log2(folds of mean expression in two groups); Probability:Probability of difference; Blast nr: NR annotation. (XLS 1400 kb)

Additional file 6:

Figure S1. Correlation analysis of the RNA-seq data and qRT-PCR result. Correlation analysis showed correlation between RNA-seq data and qRT-PCR was good (Pearson R > 0.9). The left y-axis indicates the corresponding expression data of RNA-seq (black lines). The right y-axis shows the relative gene expression levels detected by qRT-PCR (red lines). The x-axis shows the length of time (h) the eggplant peel was exposed to light. (DOCX 1112 kb)

Additional file 7:

Figure S2. Cluster analysis of DEGs with significant expression profile changes. All the DEGs were subjected to complete-linkage hierarchical clustering using a Euclidean distance metric and divided into 26 clusters. The x-axis showed the time point to time point comparison. The y-axis shows the relative log2(ratio) of each comparison. (DOCX 1568 kb)

Additional file 8:

Table S6 The 869 DEGs identified in both this study and the previous study [35]. (XLS 636 kb)

Additional file 9:

Figure S3. GO categories assigned to the 869 DEGs. The left y-axis shows the percentages of proteins identified, and the right y-axis shows the protein number. The proteins were categorized according to the annotation of GO, and the number of each category is displayed based on biological process, cellular components, and molecular functions. (DOCX 798 kb)

Additional file 10:

Table S7 KEGG pathway enrichment analysis of the 869 DEGs. (DOCX 40 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, J., He, YJ., Zhou, L. et al. Transcriptome profiling of genes related to light-induced anthocyanin biosynthesis in eggplant (Solanum melongena L.) before purple color becomes evident. BMC Genomics 19, 201 (2018). https://doi.org/10.1186/s12864-018-4587-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-018-4587-z

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us