- Open Access
Metabolome and transcriptome profiling provide insights into green apple peel reveals light- and UV-B-responsive pathway in anthocyanins accumulation
BMC Plant Biology volume 21, Article number: 351 (2021)
In nature, green apple are associated with the accumulation of chlorophyll, while red apple varieties are associated with anthocyanins accumulation. Notably, in this study, the green skin color apple variety ‘white winter pearmain’ treated with ultraviolet-B (UV-B) exhibited red skins and marked anthocyanin accumulation, while visible light could not. But there are few reports on the biosynthesis difference of anthocyanins in green apple by visible light and UV-B-treatment. Here, we explored the difference of metabolites and genes expression level in green apple by transcriptomic and metabolic.
The metabolic analysis revealed that there were 152 and 178 significantly changed metabolites in the visible light and UV-B-treated green apple, respectively, compared to the control, and flavone, flavonol, and anthocyanin were the most significantly increased; and transcriptomic analysis showed that 37,110 and 37,709 differentially expressed genes, including 382 and 475 transcription factors (TFs) were detected in light and UV-B-treatment fruit, respectively. Quantitative reverse transcription PCR (qRT-PCR) results confirmed changes in the expression levels of genes encoding metabolites involved in the flavonoid synthesis pathways. The flavonoid metabolic flux in the UV-B treatment increased the accumulation of cyanidin 3-glucoside and cyanidin 3, 5-diglucoside compared to under the light-treatment. Furthermore, we performed qRT-PCR analysis of anthocyanin biosynthesis genes and predicted the gene of MD00G1134400 (a UDP glucose-flavonoid 3–0-glucosyltransferase) may be a candidate gene for anthocyanins accumulation and highly expressed in UV-B-treatment fruit. Expression profiles of several transcription factors of the families MYB, bHLH, NAC were highly correlated with the content of the anthocyanin.
The composition and contents of anthocyanins in green apple in UV-B-treatment very greatly. A series of metabolites and candidate genes were revealed through combined analysis of metabolome and transcriptome. These results provide an important data for dissecting candidate genes and molecular basis governing green apple color formation in response to visible light and UV-B light.
In recent years, red apple have attracted the attention of consumers due to their attractive appearance and their reported nutritional value [1, 2]. Red apple varieties such as ‘Fuji’  and ‘Fortune’  can synthesize high numbers of anthocyanins and are preferred by more consumers. Conversely, green apple cannot synthesize anthocyanins during natural growth, so that they mainly exhibit green or yellow colors, for example, the ‘white winter pearmain’ [5, 6] and ‘mutsu’  apple varieties.
Annually, during apple growing seasons, particularly in the later stages of growth, sunlight plays important role in the development of fruit color [7,8,9]. In China, paper bagging is a tool applied in apple or pear fruit production [10,11,12], can not only control insect and bird pests, but can also enhance red pigmentation in the red fruit varieties in the late stages of fruit growth . In addition, the degree of red pigmentation in one fruit variety could vary significantly in one location across several years based on varying sunlight intensity and temperature in an environment. Early studies suggested that light and ultraviolet-B (UV-B) radiation regulate anthocyanin synthesis by regulating the activity of key enzymes and the gene transcription levels in anthocyanin biosynthesis pathway in red skin color fruit [13,14,15]. Furthermore, lowing the temperature is more beneficial to anthocyanin synthesis compared to high temperature in fruit peel . Some researchers have attempted to elucidate the mechanisms of fruit coloration using high-throughput techniques [7, 9, 13]. However, to date, the effect of visible light and UV-B on coloration of red skin color and green skin color apple fruit, and their molecular foundation is rarely studied.
The red coloration in red apple varieties is due to the synthesis of anthocyanin in exocarp [4, 15]. Anthocyanins are synthesized through flavonoid biosynthesis pathways, and the biosynthesis of anthocyanins involves several key enzymes and the genes expressions of coding these enzymes, it mainly contains phenylalanine ammonia-lyase (PAL), chalcone synthase (CHS), chalcone isomerase (CHI), flavanone-3-hydroxylase (F3H), dihydroflavonol 4-reductase (DFR), anthocyanidin synthase (ANS), and UDP-glucose: flavonoid 3-glucosyltransferase (UFGT) . Regulating the genes encoding these enzymes are regulated by a protein complex formed by an R2R3-MYB, a basic-helix-loop-helix, and a WD40 . The complex induces gene expression by binding to the promoters of structural gene . In numerous species, MYB transcription factors (TFs) are considered essential for the activation of anthocyanin biosynthesis related genes. In apple, MdMYB1 and MdMYB10 act as positive TFs for anthocyanin biosynthesis in fruit skin [19, 20], while MdbHLH3 plays a regulatory role during anthocyanins biosynthesis . Although the research on apple fruit coloring mainly focuses on red varieties, but few people have studied green varieties, particularly, the response of green apple varieties to visible light and UV-B exposure. Therefore, it is critical to investigate the potential mechanisms by which visible light and UV-B induce anthocyanin synthesis in green apple varieties.
In the present study, we reported the anthocyanin accumulation by metabolome and transcriptome analyses in ‘white winter pearmain’ apple fruit peel samples treated with artificial visible light (63 h), artificial UV-B radiation (63 h), and a control (CK, bagged fruit stored in the dark for 63 h) after harvesting. Based on an integrated analysis of the differential metabolites and gene expression levels, we identified the metabolites involved in visible light and UV-B-responsive reactions, and analyzed the differential regulation of structural genes involved in anthocyanins biosynthesis. The result of the present study could elucidate further the mechanisms by which visible light and UV-B regulate anthocyanin biosynthesis in apple. In addition, the metabolomic and transcriptomic data obtained would offer basic data that would facilitate further investigations on the influence of visible light- and UV-B-induced fruit coloration in green apple varieties.
Changes in fruit pigmentation patterns in the peel of ‘white winter pearmain’ fruit
Following transportation to the laboratory and upon visual inspection, the unbagged ‘white winter pearmain’ fruit were green and the bagged fruit were pale yellow. In addition, the fruit removed from the bags and subjected to light treatment were pale yellow after 63 h, while the fruit removed from the bags and subjected to UV-B radiation treatment were red after 63 h (Fig. 1A). The red coloration was accompanied by an increase of anthocyanin concentration. In addition, the contents of carotenoid, chlorophyll a and chlorophyll b in the bagged and light-treatment fruit were very low, while the chlorophyll b and anthocyanin accumulation levels increased significantly following treatment with UV-B radiation (Fig. 1B-E).
Metabolic differences among the light, UV-B, and CK peel treatments
To compare the differences of metabolite compositions among the light, UV-B, and CK peel treatments, the peel sampled were subjected to LC–MS/MS analysis. In the present study, among the 248 metabolites, there were 58 flavonoids and 8 anthocyanins, a major phenylpropane metabolic pathway, were identified and quantified in apple peel (Table S2). Principal components analysis (PCA) revealed that the CK, light, and UV-B treatments were clearly separated in the PC1–PC2 score plots (Fig. 2A). In addition, a Venn diagram of the dataset revealed that 210 metabolites were expressed differentially between the CK and light treatments, between the CK and the UV-B radiation treatment, as well as between the light and the UV-B radiation treatment (Fig. 2B). Among the 248 metabolites identified in the apple peel, 152 and 178 (61.3 and 71.8%) changed substantially in the light and the UV-B treatments, respectively, compared to in the CK treatment. In addition, 127 and 133 metabolites were upregulated, and 25 and 45 metabolites were downregulated in the light and UV-B treatments, respectively, compared to the metabolites in the CK treatment (Fig. 2C-D). Moreover, 57 metabolites were upregulated and 94 metabolites were downregulated in the UV-B treatment compared with the light treatment (Fig. 2E).
Analysis of metabolic pathways and metabolites under the light- and UV-B-treatment
To identify different metabolites in the peel of light, UV-B compared to those of CK, UV-B compared to those of light, a metabolom comparison of the materials mentioned above was carried out. A KEGG pathway enrichment analysis revealed that the significantly changed metabolites were enriched mainly in ‘phenylpropanoid metabolic processes’ including flavonoid biosynthesis, flavanone biosynthesis, flavone biosynthesis, flavonol biosynthesis, and anthocyanin biosynthesis (Fig. 3A-C). Fifty flavonoids and 8 anthocyanins with significant changes in the metabolites of flavonoid and astaxanthin were screened out from the metabolome data for heat map analysis (Fig. 3D, Table S3). The results showed that the compounds of related to flavonoids and anthocyanins increased significantly under visible light and UV-B induction, such as cyanidin 3,5-O-diglucoside, cyanidin 3-O-galactoside, delphinidin 3-O-rutinoside, quercetin and so on.
Differential expression of genes in the fruit peel
To identify the differentially expressed genes (DEGs) among the light, UV-B, and CK treatments, a transcriptomic comparison was performed. PCA score and Venn plots of gene expression among treatments revealed the effectiveness of data (Figure S1). Based on fold change ≥ 2 and FDR < 0.05 criteria, 37,110 and 37,709 DEGs were identified in the light and UV-B treatments compared with CK, respectively, and 3469 and 4212 genes, respectively, were significantly upregulated (Table S4a-b). In addition, 37,590 DEGs were identified in the UV-B compared with light, and 2762 genes were significantly upregulated (Table S4c). GO enrichment top 50 analysis of the 3355, 3428 and 2659 DEGs showed that the three major biological processes were biological process, cellular component, and molecular function. In the molecular function domain, the DEGs were enriched in catalytic activity and binding. In the cellular component domain, most of the DEGs were enriched in four categories, including apoplast, chloroplast, plastid, and thylakoid lumen (Fig. 4A-C). KEGG pathway enrichment analysis revealed that the DEGs were mainly enriched in various metabolic processes including flavonoid biosynthesis and secondary metabolism biosynthesis (Fig. 4D-F).
Changes of flavonoid and anthocyanin biosynthesis
Transcriptomic data revealed that flavonoid and anthocyanin biosynthesis genes were differentially induced by light and UV-B radiation treatments (Fig. 5A). We selected 15 significantly different expressed genes to validate the results of the transcriptomic analysis. qRT-PCR results revealed that 10 out of the 15 genes were expressed at higher levels under UV-B than under the light treatment. Particularly, the level of expression of MD00G1134400, a UFGT homologue, was enhanced under the UV-B treatment, implying the requirement of such structural enzymes in anthocyanin biosynthesis in apple (Fig. 5B).
Integrated analyses of the metabolomic and transcriptomic datasets
To reveal the regulatory mechanisms of enhanced flavonoid and anthocyanin biosynthesis in peel exposed to light and UV-B radiation, 58 metabolites (flavonoid biosynthesis, 50; and anthocyanin metabolism, 8) and 30 DEGs (flavonoid biosynthesis, 27; and anthocyanin metabolism, 3) were selected for further analyses (Fig. 6). Although a vast number of flavonoids were detected in the peel subjected to light and UV-B radiation treatments (Table S3), UV-B radiation induced higher anthocyanin concentrations than the light treatment (Fig. 1A). Particularly, the compounds that demonstrated changes, including cyanidin 3, 5-O-diglucoside (cyanin) and cyanidin 3-O-glucoside (kuromanin), showed markedly higher abundances under the UV-B than under the visible light treatment (Table S3). In addition, we arranged 30 DEGs based on their corresponding positions in the apple flavonoid biosynthesis pathway (Fig. 6). Among them, some MdCHS, MdCHI, MdF3H, MdFLS, MdDFR, MdANS, and MdUFGT homologues exhibited different expression levels.
Differential expression of anthocyanin biosynthetic genes in fruit peel
Anthocyanin in the skins investigated was biosynthesized via the flavonoid pathway (Fig. 7A). To determine the roles of anthocyanin biosynthetic genes under the light and UV-B conditions, the expression patterns of genes involved in their biosynthesis were analyzed. Significant expression changes were observed in both the light and UV-B responsive genes. Pearson’s correlation analysis showed that the genes expression levels in transcriptome verification experimrnt was significantly correlated with the reported transcription levels of anthocyanin biosynthesis related genes (P < 0.05, and P < 0.01) by light and UV-B treatment in apple peel (Table 1). We further performed qRT-PCR analysis on the reported related genes and candidate genes. The results showed that most of the genes, such as MdCHS, MdF3H, MdDFR, MdANS, and Md00G1134400 were more responsive to UV-B induction (Fig. 7B). The results suggested that the contribution of UV-B induced anthocyanin synthesis may be greater than that of visible light in green apple.
Transcription factors (TF) and light signal transduction related genes of light-, UV-B-responsive
Three hundred eighty-two differentially expressed transcription factors belonging to 52 TF families, and 475 differentially expressed transcription factors belonging to 59 TF families were identified in the light-treatment and the UV-B treatment compared with CK, respectively, and 346 differentially expressed transcription factors belonging to 48 TF families were identified in the UV-B treatment compared with light-treatment (Table S5). In addition, we observed the up-regulation and down-regulation of 28 TF families that were both regulated by light and UV-B radiation treatments, including WRKY, NAC, MYB, bHLH, and AP2/ERF families (Fig. 8). Although rather similar numbers of TFs were observed between the light and the UV-B radiation treatments, most of the TFs were upregulated in UV-B. In the gene set, AP2/ERF was the largest family with 47 members, followed by the MYB, NAC, bHLH, and WRKY families, with 46, 45, 33, and 32 members in the UV-B treatment, respectively (Fig. 8B, Table S5). And the numbers of TFs in UV-B treatment were signifficant increased than light treatment (Fig. 8C). The results showed that UV-B may increasing the accumulation of anthocyanin through upregulation of TFs.
In addition, the transcriptional level of MYB, bHLH, and WD40 group genes were generally higher under UV-B treatment, and UV-B treatment showed tha same trend as the content anthocyanin. Here, we identified MdHY5 (MD08G1147100) gene by transcriptom. The transcript level of MdHY5 decreased under UV-B treatment in apple, which was poorly correlated with anthocyanin level. However, the transcript abundance of MdUVR8 (MD12G1149100) was consistent with anthocyanin content, and UV-B induced the expression of MdUVR8 was higher compared to light treatment in apple (Fig. 9).
Coloration of bagged, unbagged, artificial light- and UV-B-treated apple fruit
Both light and UV-B radiation can induce coloration in red apple fruit [22, 23]. In nature, the red apple varieties can synthesize anthocyanins, which would appear red, and are influenced significantly by sunlight. Conversely, the green apple cannot synthesize anthocyanins and exhibits a green color even following exposure to sunlight . Nevertheless, only a few studies have explored the responses of green apple fruit varieties to visible light and UV-B radiation. Previous studies have shown that bagging could facilitate fruit coloration . Therefore, we bagged ‘white winter pearmain’ apple in the course of development. We observed that prolonged exposure to UV-B radiation could induce the coloration of green apple, while prolonged exposure to visible light had no effect (Fig. 1). The results indicated that green apple could synthesize red pigment under UV-B induction similar to red apple.
Metabolites involved in the induction of coloration under visible light and UV-B radiation
Previous studies have reported that sunlight influences physiological change in fruit [11, 25]. In particular, sunlight influences the synthesis of phenols in late stages of fruit maturation , suggesting sunlight plays an important role in fruit development. There is evidence of an increase in anthocyanins and flavonoid compounds in fruit exposed to sunlight and that acclimation to sunlight could have positive effects on the growth of red apple varieties . However, green fruit varieties do not exhibit responses in anthocyanin synthesis following exposure to sunlight . In the present study, visible light and UV-B radiation could increase the accumulation of metabolites in fruit, mainly in the metabolic pathways of amino acids and their derivatives, as well as the metabolic pathways of anthocyanins, flavonoids (flavanones and flavonol), lipids, and organic acids, and their derivatives. In particular, the metabolites involved in the flavonoid metabolic pathway were up-regulated significantly (Fig. 2C, D). KEGG analysis showed that visible light and UV-B radiation induced metabolic and biosynthesis pathways, particularly flavonoid and anthocyanin biosynthesis pathways (Fig. 3A, B). In addition, visible light and UV-B radiation induced the synthesis of different flavonoids, and anthocyanins in 58 metabolites (Fig. 3D). The results indicated that there were significant differences in metabolite synthesis between visible light and UV-B radiation treatments. We further identified the metabolites associated with the anthocyanin biosynthesis pathway in the metabolome data, and observed that both visible light and UV-B radiation could induce the synthesis of delphinidin 3-rutinoside and cyanidin 3-glucoside, but only UV-B treatment could induce the formation of cyanidin 3, 5-diglucoside (Fig. 6, Table S3). Notably, cyanidin 3, 5-diglucoside is the dominant anthocyanin in rose flowers . In the present study, cyanidin 3-glucoside and cyanidin 3, 5-diglucoside were the major anthocyanin components in the fruit treated with UV-B radiation. In addition, the concentrations of cyanidin 3, 5-diglucoside were 12.3-fold the concentrations of cyanidin 3-glucoside. Conversely, light treatment increased cyanidin 3-glucoside concentrations slightly and did not enhance cyanidin 3, 5-diglucoside production, which indicated that the red skin color could be attributed to the higher level of production of cyanidin 3, 5-diglucoside.
Differentially expressed genes and other transcription factors are involved in fruit coloration
Sunlight facilitates the accumulation of different metabolic compounds by regulating plant gene expression . Several studies have demonstrated that light induces anthocyanin synthesis through the expression of key enzyme genes in the anthocyanin biosynthesis pathway [16, 29]. In addition, numerous transcription factors associated with anthocyanin have been studied such as MdMYB1 and MdMYB10 [19, 20, 30, 31]. In the present study, Go and KEGG analyses revealed that the gene expression changes under visible light and UV-B radiation treatments were mainly in the metabolic pathways and the secondary metabolite pathways, particularly the flavonoid synthesis pathway (the results were consistent with the metabolome data) (Fig. 4). Transcriptome data showed that 6 out of the 15 genes involved in flavonoid biosynthesis were induced by both visible light and UV-B radiation, with one genes induced significantly by light, and the other five genes induced significantly by UV-B. In addition, three genes were closely associated with the anthocyanin biosynthesis pathway and were induced by both light and UV-B. GenBank analysis revealed that the three genes belonged to the same gene family, UFGT. At present, MD01G1234400 function (GenBank: AF117267) has been reported in several studies, and it is regulated by UV-B and temperature . Our results showed that three genes, including MD01G1234400, MD07G1306900, and MD00G1134400, were induced by visible light and UV-B radiation. In particular, MD00G1134400 expression was induced more significantly by UV-B. qRT-PCR was used to verify gene expression, and the results were consistent with transcriptomic data (Fig. 5B).
We applied both metabolomic and transcriptomic techniques in the analysis of the effects of visible light and UV-B radiation on metabolites and genes expression in the flavonoid biosynthesis pathways. We identified 30 DEGs, including MdCHS, MdCHI, MdF3H, MdFLS, MdDFR, MdANS, and MdUFGT that responded to both visible light and UV-B radiation (Fig. 6). In addition, visible light and UV-B radiation induced the expression of TF family genes, mainly including MYB, bHLH, and WRKY. Particularly, the number and the transcription level of TFs induced by UV-B radiation was significant (Figs. 8 and 9). The results indicated that light and UV-B play important roles in the synthesis of phenylpropane, particularly in the flavonoid and anthocyanin synthesis pathways, and the effect of UV-B was apparent in the green apple fruit.
Differential regulation of genes in anthocyanin biosynthesis
Anthocyanins are synthesized mainly via the flavonoid biosynthesis pathway . In recent years, many studies have explored the key enzyme genes involved anthocyanin synthesis. Previous studies have reported that the anthocyanin synthesis pathway genes are regulated by sunlight in ‘Fuji’ apple [3, 9]. Ubi et al.  showed that the genes are induced strongly by low temperature. To verify that the key enzyme genes in the anthocyanin biosynthesis pathway are influenced by visible light and UV-B radiation in green apple, we analyzed the expression of associated genes. According to the results, both visible light and UV-B radiation could induce gene expression in the anthocyanin biosynthesis pathway. Particularly, UV-B radiation induced gene expression downstream significantly (Fig. 7). The results of the metabolomic and transcriptomic analyses revealed that UV-B induces the expression of anthocyanins such as cyanidin 3, 5-diglucoside in green apple, and the identification of the candidate genes involved requires further exploration. The correlation analysis confirmed that flavonoid and anthocyanin-associated genes in the transcriptome were closely related to the reported anthocyanin biosynthesis genes (Table 1). Therefore, it can be said that these genes were related to anthocyanin biosynthesis when apple are exposed to UV-B. In this study, green apple fruit were treated using visible light and UV-B radiation, and these results provide a theoretical basis that could facilitate the understanding of the effect of light and UV-B on the green fruit, and offer novel insights that could facilitate the mining of anthocyanin synthesis genes (Fig. 10).
In this study, we analyzed the metabolome and transcriptome of ‘white winter pearmain’ apple fruit peel after harvesting and exposure to visible light and UV-B radiation. Using an integrated analysis of the differential metabolites and genes expression levels, we identified the metabolites involved in responses to visible light and UV-B, and further analyzed the differential expression of genes involved in anthocyanin biosynthesis using qRT-PCR. After data mining from transcriptome, a apple UDP glucose-flavonoid 3–0-glucosyltransferase gene, namely Md00G1134400, was found to be differentially expression of UV-B-responsive involved in anthocyanin biosynthesis in green apple fruit. In conclusion, the UV-B-induced red color of green apple may be caused by the regulation of UV-B-responsive signals, including UVR8 and members of the MBW complex, to regulate the transcriptional expression of structural genes of anthocyanin biosynthesis. So our results provide insights into the synthesis anthocyanin metabolic products and the basic gene-expression networks in green fruits exposed to UV-B radiation. Further work on gene function validation is required to have deep insights of the genetic and molecular mechanisms underlying green apple fruit coloration.
Materials and methods
Plant materials and experiments
The ‘white winter pearmain’ apple fruit were obtained from a commercial orchard in Taian, Shandong Province, China. In total, 150 fruitlets were covered with double layers of yellow-black paper bags 45 days after full blossom (DAFB). The fruit remained bagged before harvest (160 DAFB) and were transported to the laboratory with the remaining unbagged fruit. The fruit of artificial visible light-treatment (bagged fruit exposed to visible light [300–320 μmol m-2 s-1]) were placed approximately 10 cm (apple surface to visible light lamp), the wavelength of light was 400-700 nm, and the fruit of and artificial UV-B treatment (bagged fruit exposed to UV-B radiation [10–12 kJ m-2 s-1]) were placed approximately 25 cm (apple surface to UV-B lamp), the wavelength of UV-B lamp was 313 nm. The above two groups of treatments were performed on 60 fruit for 63 h after bag removal at room temperature (21 °C), respectively, and each treatment had three biological replicates, each replicate had 20 fruit. The control group (30 fruit) was included the bagged fruit maintained in the dark for 63 h. Fruit peel were scraped and immediately collected, quickly frozen in liquid nitrogen, and stored in a freezer at -80 °C. Three apple were used for each treatment, with three replicates. The phenotypes of the unbagged, bagged, bagged-light treatment, and bagged-UV-B treatment fruit were observed and recorded.
Measurement of chlorophyll, carotenoid, and total anthocyanin concentrations
Pigments were extracted from fruit peel samples (0.3 g) using 3 ml 95% ethanol at room temperature (21 °C) in the dark until the peel tissues bleached completely according to previously described  with minor modification. The extracts were centrifuged at 8000 × g for 2 min, and the supernatants were collected and their absorbance was measured using a UV-1780 spectrophotometer (Shimadzu, Japan) at 665 nm, 649 nm, and 470 nm. The formulae for calculating chlorophyll concentrations were as following: chlorophyll a = 13.95*A665 – 6.88*A649; and chlorophyll b = 24.96*A649 – 7.32*A665. The formulae for calculating carotenoid concentrations were as following: carotenoids = (1000*A470 – 2.05*Ca – 114.8*Cb)/245 and the results were expressed as mg g−1.
Anthocyanin was extracted using 0.3 g of the frozen peel tissues in 3 ml 95% methanol:1.5 M hydrochloric acid (85:15) for 24 h at 4 °C in the dark according to previously described  with minor modifications. The extract was centrifuged at 12,000 × g for 5 min, and the supernatant was collected and analyzed using a UV-1780 spectrophotometer at 530, 620, and 650 nm. The anthocyanin concentrations were calculated using the following formula: 46,200*[(A530 -A620)-0.1*(A650-A620)] and the results were expressed as ug g−1.
RNA-Seq analysis, cDNA synthesis, and data analysis
Total RNA was extracted from fruit peel of the CK, and the Light and UV-B treatment samples, and cDNA synthesis was performed using the RNA library prep kit (Tiangen, Beijing, China). At least three biological repeats were collected and mixed. The RNA-seq and data analysis were performed by Metware Biotechnology Co., Ltd. (Wuhan, China) as previously described [35,36,37,38,39,40,41,42,43] based on available genetic databases, including NCBI non-redundant protein database, the SwissProt database, the TrEMBL database, the Pfam database, the KOG database, the KEGG database, and the GO database, based on log2 fold changes > 1, and a false discovery rates (FDR) < 0.05. The experiments were repeated three times.
Multiple reaction monitoring (MRM)
Multiple reaction monitoring (MRM) was performed by Metware Biotechnology Co., Ltd. (Wuhan, China). The freeze-dried apple peel were crushed using an MM 400 mixer mill (Retsch, Haan, Germany) for 1.5 min at 30 Hz. The powder (100 mg) was weighed and extracted overnight at 4 °C. After centrifugation at 10,000 × g for 10 min, the extracts were absorbed and filtered. The sample extracts were analyzed using a liquid chromatography coupled with tandem mass spectrometry (LC–MS/MS) system (HPLC, High Performance Liquid Chromatography, Shim-pack UFLC SHIMADZU CBM30A system, https://www.shimadzu.com.cn/; MS/MS, Tandem mass spectrometry, Applied Biosystems 6500 QTRAP, http://www.appliedbiosystems.com.cn/). The analytical conditions were as follows: HPLC: column, Waters ACQUITY UPLC HSS T3 C18 (1.8 μm, 2.1 mm × 100 mm); solvent system, water (0.04% acetic acid): acetonitrile (0.04% acetic acid); gradient program, 95:5 V/V at 0 min, 5:95 V/V at 11.0 min, 5:95 V/V at 12.0 min, 95:5 V/V at 12.1 min, 95:5 V/V at 15.0 min; flow rate, 0.4 mL/min; temperature, 40 °C; injection volume: 2 μl. The mass spectrum conditions mainly included an electrospray ionization (ESI) temperature of 500 °C, a mass spectrum voltage of 5500 V, and curtain gas pressure of 25 psi, with the collision-activated dissociation parameter set as high. In the triple quadrupole mass spectrometer, each ion pair was scanned according to the optimized declustering potential and collision energy . A specific set of MRM transitions was monitored for each period according to the metabolites eluted within the period. The MRM was performed in triplicate for each treatment. The experiments were repeated three times.
RNA extraction, DNA sequencing, and real-time quantitative reverse transcription (qRT-PCR) analysis
Samples stored in the freezer at -80 °C were used for total RNA extraction using a total RNA isolation system (Tiangen, Beijing, China). First-strand cDNAs were synthesized using a First-strand cDNA Synthesis Kit (Tiangen, Beijing, China). The full-length coding sequences of genes were isolated from ‘white winter pearmain’ cDNA. Real-time quantitative reverse transcription (qRT-PCR) analysis was performed according to MIQE guidelines [45, 46]. qRT-PCR was performed on a Bio-Rad CFX96TM Real-time PCR System using SYBR Real Master Mix (Transgen, Beijing, China) under the following PCR thermal cycling conditions: predenaturation at 95 °C for 30 s; followed by 39 cycles of 95 °C for 5 s, 60 °C for 15 s, and 72 °C for 20 s. MdActin (GenBank Accession Number: AB638619) was used as the housekeeping gene. The sequences of the primers used are listed in Supplementary Table S1. Three biological replicates were performed for each gene, and the standard curve method was applied in statistical analysis.
Least significant differences (p < 0.05) were calculated for mean separation using Data Processing System procedures (Zhejiang University, Zhejiang, China). Standard errors were calculated and figures illustrated using Sigma Plot v12.0 (Systat Software Inc., San Jose, CA, USA). The data pearson correlation coefficients was analyzed using SPSS 16.0 software. For RNA-seq, differential expression analysis between two samples was performed using the DESeq2 package. After the differential expression analysis, the multiple hypothesis test correction of the P-value was performed using the Benjamini–Hochberg procedure based on a log2 fold change > 1 and an FDR < 0.05. The procedures were repeated three times.
Availability of data and materials
We have uploaded the RNA sequencing data to NCBI and the link is https://www.ncbi.nlm.nih.gov/bioproject/PRJNA623676. The data sets supporting the results of this article are include within the article and supplementary table.
Quantitative reverse transcription PCR
UDP glucose-flavonoid 3–0-glucosyltransferase
Days after full blossom
Multiple reaction monitoring
Principal components analysis
Differentially expressed genes
Significantly changed metabolites
Kyoto Encyclopedia of Genes and Genomes
Espley RV, Bovy A, Bava C, Jaeger SR, Tomes S, Norling C, Putterill J. Analysis of genetically modified red-fleshed apples reveals effects on growth and consumer attributes. Plant Biotechnol J. 2013;11(4):408–19.
Sansavini S, Donati F, Costa F, Tartarini S. Advances in apple breeding for enhanced fruit quality and resistance to biotic stresses: new varieties for the European market. J Fruit Ornam Plant Res. 2004;12:13–52.
Jakopic J, Stampar F, Veberic R. The influence of exposure to light on the phenolic content of ‘Fuji’apple. Sci Hortic. 2009;123(2):234–9.
Feng F, Li M, Ma F, Cheng L. Phenylpropanoid metabolites and expression of key genes involved in anthocyanin biosynthesis in the shaded peel of apple fruit in response to sun exposure. Plant Physiol Biochem. 2013;69:54–61.
Ding RR, Che XK, Liu HZ, Du BY, Dong KT, Zhang YH. Effects of 1-MCP and storage temperature on transcription of mevalonate (MVA) enzyme genes of α-farnesene in ‘White Winter Pearmain’apples fruit. Sci Hortic. 2020;259:108841.
Ding RR, Du BY, Zhang YH. Conjugated trienols and programmed cell death are more closely related to superficial scald than reactive oxygen species in apple fruit stored at low temperature. Sci Hortic. 2019;246:597–603.
Bai S, Sun Y, Qian M, Yang F, Ni J, Tao R, Teng Y. Transcriptome analysis of bagging-treated red Chinese sand pear peels reveals light-responsive pathway functions in anthocyanin accumulation. Sci Rep. 2017;7:63.
Rumainum IM, Worarad K, Yamaki Y, Yamane K. Effects of developmental stages, light, and an auxin polar transport inhibitor on the skin and flesh pigmentation of red-fleshed peach fruit. Hortic J. 2016;85(2).
Vimolmangkang S, Zheng D, Han Y, Khan MA, Soria-Guerra RE, Korban SS. Transcriptome analysis of the exocarp of apple fruit identifies light-induced genes involved in red color pigmentation. Gene. 2014;534(1):78–87.
Chen J, Gao W, Lv D, Cai M, Qin S, Ma H, Liu G. The effects of bagging on the fruit skin development process of Red Fuji apple. J Fruit Sci. 2009;26:217–21.
Huang C, Yu B, Teng Y, Su J, Shu Q, Cheng Z, Zeng L. Effects of fruit bagging on coloring and related physiology, and qualities of red Chinese sand pears during fruit maturation. Sci Hortic. 2009;121:149–58.
Xia J, Zhang Z, Qu S. Changes of peel pigments and flesh sugar contents in bagging Malus pumila during development stage. Acta Bot Boreal-Occid Sin. 2010;30:1675–80.
Czemmel S, Höll J, Loyola R, Arce-Johnson P, Alcalde JA, Matus JT, Bogs J. Transcriptome-wide identification of novel UV-B-and light modulated flavonol pathway genes controlled by VviMYBF1. Front Plant Sci. 2017;8:1084.
Honda C, Kotoda N, Wada M, Kondo S, Kobayashi S, Soejima J, Moriguchi T. Anthocyanin biosynthetic genes are coordinately expressed during red coloration in apple skin. Plant Physiol Biochem. 2002;40:955–62.
Ubi BE, Honda C, Bessho H, Kondo S, Wada M, Kobayashi S, Moriguchi T. Expression analysis of anthocyanin biosynthetic genes in apple skin: effect of UV-B and temperature. Plant Sci. 2006;170(3):571–8.
Xie R, Zheng L, He S, Zheng Y, Yi S, Deng L. Anthocyanin biosynthesis in fruit tree crops: genes and their regulation. Afr J Biotech. 2011;10(86):19890–7.
Xu W, Dubos C, Lepiniec L. Transcriptional control of flavonoid biosynthesis by MYB–bHLH–WDR complexes. Trends Plant Sci. 2015;20(3):176–85.
Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010;5(10):573–81.
Espley RV, Hellens RP, Putterill J, Stevenson DE, Kutty-Amma S, Allan AC. Red colouration in apple fruit is due to the activity of the MYB transcription factor, MdMYB10. Plant J. 2007;49:414–27.
Hu DG, Sun CH, Ma QJ, You CX, Cheng L, Hao YJ. MdMYB1 regulates anthocyanin and malate accumulation by directly facilitating their transport into vacuoles in apples. Plant Physiol. 2016;170:1315–30.
Xie XB, Li S, Zhang RF, Zhao J, Chen YC, Zhao Q, Hao YJ. The bHLH transcription factor MdbHLH3 promotes anthocyanin accumulation and fruit colouration in response to low temperature in apples. Plant Cell Environ. 2012;35:1884–97.
Jaakola L. New insights into the regulation of anthocyanin biosynthesis in fruit. Trends Plant Sci. 2013;18(9):477–83.
Takos AM, Jaffé 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:1216–32.
Bakhshi D, Arakawa O. Induction of phenolic compounds biosynthesis with light irradiation in the flesh of red and yellow apples. J Appl Hortic. 2005;8(2):101–4.
Opara LU, Studman CJ, Banks NH. Sunlight affects the incidence of internal ring cracking and other physical attributes of ‘Gala’apples. J Tree Fruit Prod. 1997;2(1):45–52.
Sun S, Xin L, Gao H, Wang J, Li P. Response of phenolic compounds in ‘Golden Delicious’ and ‘Red Delicious’ apples peel to fruit bagging and subsequent sunlight re-exposure. Sci Hortic. 2014;168:161–7.
Mikanagi Y, Sait N, Yokoi M, Tatsuzawa F. Anthocyanins in flowers of genus rosa, sections cinnamomeae (=rosa), chinenses, gallicanae and some modern garden roses. Biochem Syst Ecol. 2000;28(9):887–902.
Stracke R, Ishihara H, Huep G, Barsch A, Mehrtens F, Niehaus K, Weisshaar B. Differential regulation of closely related R2R3-MYB transcription factors controls flavonol accumulation in different parts of the Arabidopsis thaliana seedling. Plant J. 2010;50:660–77.
Henry-Kirk RA, Plunkett B, Hall M, McGhie T, Allan AC, Wargent JJ, Espley RV. Solar UV light regulates flavonoid metabolism in apple (Malus x domestica). Plant Cell Environ. 2018;41:675–88.
Bai S, Tao R, Tang Y, Yin L, Teng Y. Bbx 16, a b-box protein, positively regulates light-nduced anthocyanin accumulation by activating myb 10 in red pear. Plant Biotechnol J. 2019;17(10):1985–97.
Ma C, Jing C, Chang B, Yan J, Liang B, Liu L, Zhao Z. The effect of promoter methylation on MdMYB1 expression determines the level of anthocyanin accumulation in skins of two non-red apple cultivars. BMC Plant Biol. 2018;18(1):108.
Kim SH, Lee JR, Hong ST, Yoo YK, An G, Kim SR. Molecular cloning and analysis of anthocyanin biosynthesis genes preferentially expressed in apple skin. Plant Sci. 2003;165(2):403–13.
Lichtenthaler HK. Chlorophylls and carotenoids: pigments of photosynthetic biomembranes. Method Enzymol. 1987;148C:350–82.
Zhao J, Xie X, Shen X, Wang Y. Effect of sunlight-exposure on antioxidants and antioxidant enzyme activities in ‘d’Anjou’pear in relation to superficial scald development. Food Chem. 2016;210:18–25.
Conesa A, Madrigal P, Tarazona S, Gomez-Cabrero D, Cervera A, McPherson A, Mortazavi A. A survey of best practices for RNA-seq data analysis. Genome Biol. 2016;17(1):13.
Gel B, Serra E. karyoploteR: an R/Bioconductor package to plot customizable genomes displaying arbitrary data. Bioinformatics. 2017;33(19):3088–90.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Yamanishi Y. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36(1):D480–4.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2013;30:923–30.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Tatusov RL, Koonin EV, Lipman DJ. A genomic perspective on protein families. Science. 1997;278(5338):631–7.
Tatusov RL, Galperin MY, Natale DA, Koonin EV. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000;28(1):33–6.
Varet H, Brillet-Guéguen L, Coppée JY, Dillies MA. SARTools: a DESeq2-and edgeR-based R pipeline for comprehensive differential analysis of RNA-Seq data. PLoS One. 2016;11:e0157022.
Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57.
Chen W, Gong L, Guo Z, Wang W, Zhang H, Liu X, Luo J. A novel integrated method for large-scale detection, identification, and quantification of widely targeted metabolites: application in the study of rice metabolomics. Mol Plant. 2013;6(006):1769–80.
Bustin SA, Beaulieu JF, Huggett J, Jaggi R, Kibenge FS, Olsvik PA, Toegel S. MIQE precis: practical implementation of minimum standard guidelines for fluorescence-based quantitative real-time PCR experiments. BMC Mol Biol. 2010;11(1):74.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Vandesompele J. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4):611–22.
This work was supported by the National Natural Science Foundation of China (31370359) and Natural Science Foundation of Shandong Provincial [grant number ZR2020MC028]. Funds were used for the experimental conduction and the open access payment. No funder has participated in data analysis and interpretation or writing of the manuscript.
Ethics approval and consent to participate
This article did not contain any studies with human participants or animals performed by any of the authors.
Consent for publication
This research is not applicable to consent for publication.
The authors declare no competing financial interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Primers for qRT-PCR.
a. Significantly changed metabolites in Light compared with CK. b. Significantly changed metabolites in UV-B compared with CK. c. Significantly changed metabolites in Light compared with UV-B.
Abundance of flavonoids and anthocyanins quantified by LC-MS/MS.
a-1. All differentially changed genes in Light compared with CK. a-2. Significantly differentially expressed genes in Light compared with CK. b-1. All differentially changed genes in UV-B compared with CK. b-2. Significantly differentially expressed genes in UV-B compared with CK. c-1. All differentially changed genes in UV-B compared with Light. c-2. Significantly differentially expressed genes in UV-B compared with Light.
a. Differentially expression transcription factors of Light compared with CK. b. Differentially expression transcription factors of UV-B compared with CK. c. Differentially expression transcription factors of Light compared with UV-B.
Differentially expressed genes (DEGs) under the dark (CK), visible light and UV-B radiation treatments. (A) PCA score plot of genes profiles from the Light, UV-B, and CK. (B) Venn plot of genes profiles from the Light, UV-B, and CK.
About this article
Cite this article
Ding, R., Che, X., Shen, Z. et al. Metabolome and transcriptome profiling provide insights into green apple peel reveals light- and UV-B-responsive pathway in anthocyanins accumulation. BMC Plant Biol 21, 351 (2021). https://doi.org/10.1186/s12870-021-03121-3
- Flavonoids and anthocyanins accumulation
- Visible light
- Metabolomic and transcriptome