- Research article
- Open Access
Metabolome and transcriptome analyses reveal chlorophyll and anthocyanin metabolism pathway associated with cucumber fruit skin color
BMC Plant Biology volume 20, Article number: 386 (2020)
Fruit skin color play important role in commercial value of cucumber, which is mainly determined by the content and composition of chlorophyll and anthocyanins. Therefore, understanding the related genes and metabolomics involved in composition of fruit skin color is essential for cucumber quality and commodity value.
The results showed that chlorophyll a, chlorophyll b and carotenoid content in fruit skin were higher in Lv (dark green skin) than Bai (light green skin) on fruit skin. Cytological observation showed more chloroplast existed in fruit skin cells of Lv. A total of 162 significantly different metabolites were found between the fruit skin of the two genotypes by metabolome analysis, including 40 flavones, 9 flavanones, 8 flavonols, 6 anthocyanins, and other compounds. Crucial anthocyanins and flavonols for fruit skin color, were detected significantly decreased in fruit skin of Bai compared with Lv. By RNA-seq assay, 4516 differentially expressed genes (DEGs) were identified between two cultivars. Further analyses suggested that low expression level of chlorophyll biosynthetic genes, such as chlM, por and NOL caused less chlorophylls or chloroplast in fruit skin of Bai. Meanwhile, a predicted regulatory network of anthocyanin biosynthesis was established to illustrate involving many DEGs, especially 4CL, CHS and UFGT.
This study uncovered significant differences between two cucumber genotypes with different fruit color using metabolome and RNA-seq analysis. We lay a foundation to understand molecular regulation mechanism on formation of cucumber skin color, by exploring valuable genes, which is helpful for cucumber breeding and improvement on fruit skin color.
Fruit skin color is an essential trait with commercial values, mainly determined by content and composition of anthocyanins and chlorophyll [1, 2]. Chlorophyll provides green pigmentation and comprises with chlorophyll a and chlorophyll b molecules. Chlorophyll metabolism can be classified into three major steps: chlorophyll synthesis, chlorophyll cycle and chlorophyll degradation. A series of important enzymes were involved in chlorophyll metabolism, such as glutamyl-tRNA reductase (HemA), porphobilinogen synthase (HemB), magnesium chelatase subunit H (chlH), magnesium-protoporphyrin O-methyltransferase (chlM), protochlorophyllide reductase (por), chlorophyll b reductase (NOL) [3, 4]. Most fruit skin was caused by chlorophyll metabolism, which exhibit green color during the fruit early development, whereas the predominant colorations of yellow, orange and red show in the post stage [5,6,7,8].
Anthocyanins, the most prominent pigment influencing fruit color, were catalyzed by complex enzymes from phenylpropanoid and flavonoid biosynthetic pathways. A wide range of constructive genes were involved in the anthocyanin biosynthesis, such as phenylalanine ammonia lyase (PAL), 4-coumarate: coenzyme a ligase (4CL), chalcone synthase (CHS) and anthocyanidin synthase (ANS) [9,10,11]. Among them, PAL is an essential factor during the anthocyanin synthesis . Flavonoid secondary metabolites are synthesized by a branched pathway of flavonols and anthocyanins synthesis. Previous study reported that various flavonoids exert crucial roles in protecting against UV-light and phytopathogens, development of male fertility, and transport of auxin . Enzymes involved in anthocyanins and flavonoid synthesis are multi-enzyme complex , and pigments tend to accumulate in vacuole (anthocyanins and proanthocyanidins) or cell wall (phlobaphenes) .
Cucumber fruit skin color has great effect on commodity sale and varietal improvement. Previous studies concerning cucumber fruit skin color mainly focus on inheritance and gene primary mapping, such as white fruit skin gene (w), dark green fruit skin gene (DG), green fruit skin gene (dg), yellow green fruit skin gene (yg), and dull fruit skin light green fruit skin gene [16, 17]. The w was rapidly mapped to a 33.0-kb region by two SNP-based markers, ASPCR39262 and ASPCR39229 . However, the molecular mechanism and pigment metabolism of fruit skin color in cucumber is unclear.
The combination of different omics helps us deeply understand several crucial genes involved in plant growth, development, and responses to different stresses [19, 20]. For instance, combined transcriptomic and metabolomics profiling offered some cues in explaining plant phenotype [21,22,23]. Through comparative transcriptomic analysis, reports showed that several novel genes functions were involved in the flavonoid  and other biochemical pathways . In addition, metabolome efficiently analyzed genes roles involved in metabolic pathway and provided essential information on genes exploring . The comparative omics has been successfully applied in fruits to clarify the relationship between different secondary metabolites and expressed genes . However, until now, reports on regulation mechanism of cucumber fruit skin color by transcriptomic and metabolomics analysis still lack.
The aim of our study was to excavate the genes involved in development of cucumber fruit skin color using conjoint analysis. Two high-inbred cucumber genotypes, ‘Lv’ with dark green skin and ‘Bai’ with light green skin from South China type cucumber variety were applied. Comparison results showed that much more content of anthocyanins, flavone, and flavonols in the fruit skin of Lv compared with Bai. In addition, we detected that the key structural genes, transcription factors and other regulators during chlorophyll and anthocyanins biosynthetic pathways. We offered crucial information on fruit skin color and its complex effect on cucumber fruit quality.
Phenotype analysis of lv and Bai
Obvious differences were found between Lv and Bai in the young fruit skin color, the fruit skin color of Lv is dark green but Bai is light green (Fig. 1a, b). The content of chlorophyll a and chlorophyll b were 0.99 mg/g and 0.90 mg/g in Bai, respectively, which were significantly lower than the Lv (Fig. 1c). The result of carotenoid is consists with chlorophyll a and chlorophyll b, the carotenoid content was higher in Lv than Bai (Fig. 1c). These results indicating more pigments accumulated in Lv fruit skin.
The above results indicated that more pigments accumulated in Lv fruit skin, which prompted us to further determine whether difference of chloroplasts in Lv and Bai cell. Through transmission electron microscopy (TEM) assay, we found that less chloroplast existed in Bai cells than Lv (Fig. 2a-c), and the number of thylakoid in a chloroplast of Bai (Fig. 2d-f) was less than Lv, these result was consistent with quantitative analysis of chlorophyll a and chlorophyll b.
The paraffin section assay was carried out to observe arrangement of skin epidermal cells. The results showed that epidermal cells in Lv were more closely arranged than Bai (Fig. 3a, b). The single cell area and single cell perimeter of Bai were both lager than Lv (Fig. 3c, d). In addition, the surface cells on the Lv fruit skin were smaller than Bai in a same field of view by scanning electron microscope (SEM) assay (Fig. 3e, f, S1).
In order to excavate metabolites during the process of cucumber fruit development (Fig. 1), a metabolome program was performed in this study. Combing detection of total ions current (TIC) and multiple reactions monitoring (MRM) profiles, we finally identified 162 significant metabolites (135 up- regulated and 27 down-regulated) between Lv and Bai samples (Fig. 4a), including: 40 flavones, 9 flavanones, 7 flavonols, 6 anthocyanins, and other compounds (Table S1). The representative metabolites, especially anthocyanins, flavones, and flavonols were listed in Table 1.
Functional analysis of metabolites
Six rosinidin O-hexoside, cyanidin O-acetylhexoside, malvidin 3-O-glucoside, malvidin 3, 5- diglucoside, peonidin O-hexoside, and peonidin were identified and all these anthocyanins were significantly decreased in Bai fruit skin compared with Lv. In Bai, peonidin and cyanidin O-malonylhexoside were decreased with 0.00035- and 0.16-fold increments in contrast to Lv, indicating that lower content of anthocyanin partly caused slight hue of Bai (Table 1). Most flavonols were found with 0.006- to 0.16- fold augment in Bai except fustin, while content of fustin was prominently increased 981.85-fold in Bai compared with Lv. Flavones were detected to be the maximum number of metabolites among metabolites with the significant content changes between two cucumber genotypes. Among these, chrysoeriol O-hexosyl-O-rutinoside, and tricetin O-malonylhexoside, luteolin O-sinapoylhexoside demonstrated significantly higher content in Lv, while only tricin O-glucuronic acid was 3.62-fold increase in Bai (Table 1). In addition, KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis demonstrated that different metabolites were mostly enriched in flavonoid biosynthesis and tryptophan metabolism, indicating flavonoid influenced fruit skin color development to some extent (Fig. 4b).
Identification of differently expressed genes (DEGs) by transcriptome
Total RNA from cucumber fruit skin were used for construction of cDNA libraries. After removing adaptor-containing raw reads and low-quality reads, the total number of clean reads was about 24 million for Lv and Bai (Table S2). These clean reads were subsequently mapped to cucumber 9930 genome (Huang et al., 2009). Approximately 90% clean reads were mapped to the reference cucumber genome, with more than 98% uniquely mapped (Table S2). The correlation coefficients in gene expression level from three biological replicates of each line were more than 0.84 (Fig. S2A), and principal component analysis (PCA) showed that biological replications clustered together (Fig. S2 B). The correlation coefficients and PCA suggested that expression patterns have similarity between replicate samples (Fig. S2). In total, 4516 DEGs with 2417 up-regulated and 2099 down-regulated genes were identified in Lv vs Bai.(Fig. 5a; Table S3). Combing transcriptome analysis, 205 DEGs belonged to 44 families encoding transcription factors (TFs), including 87 and 118 DEGs expressed down-regulation and up-regulation in Bai compared with Lv, respectively (Fig. S3). The AP2/ERF, bHLH, MYB, NAC and WRKY families were the top five TF in DEGs (Fig. S3). A total of 15 genes were selected to confirm RNA-seq data by using qRT-PCR, including 9 and 6 genes were selected from down-regulatin and up-regulation, respectively. The qRT-PCR results were consistent with RNA-seq data (Fig. S4). In addition, Csa3G904140 was detected different expressed in the Lv and Bai, and Csa3G904140 is control immature fruit color of cultivated cucumber .
Functional analysis of DEGs
In order to understand the role of DEGs in the formation of fruit skin color, three categories were classified including biological process, molecular function, and cellular components using GO (gene ontology) standardized classification system, and total of 67 GO were significantly enriched. In biological processes category, 46 GO terms were significantly enriched in DEGs, such as thylakoid membrane organization, photosynthesis and chlorophyll biosynthetic process. In molecular function category, two GO categories, including pigment binding and chlorophyll binding were found to be enriched. In cellular component category, 19 GO terms, such as photosystem I, photosystem II, plastoglobule, chloroplast envelope, chloroplast, microtubule, chloroplast stroma and chloroplast thylakoid, were identified to enrich in DEGs (Table S4). Then, we used KEGG pathway database to examine the DEGs-associated pathways. The top 20 pathway enrichment of annotated DEGs across the comparisons of Lv and Bai was shown in Fig. 5b. Related genes of carbon mechanism, amino sugar and nucleotide sugar metabolism, photosynthesis, porphyrin chlorophyll metabolism and phenoylpropanoid biosynthesis were intensively enriched (Fig. 5b).
The GO and KEGG analysis results indicated that DEGs involved in chlorophyll metabolism-related pathway, these results are consist with chlorophyll a and chlorophyll b difference between Lv and Bai. Therefore, we further studied DEGs participate in chlorophyll metabolism in detail and established a predicted chlorophyll biosynthetic pathway (Fig. 6). Fourteen DEGs were identified in chlorophyll biosynthetic pathway. Interestingly, most these DEGs were down-regulated expression in Bai compared to Lv, except one DEG (Csa7G068600).
Regulatory network of predicted flavonoid, and anthocyanidin biosynthetic pathways
In order to better understand the relationship between metabolites and genes in predicted flavonoid biosynthesis between Lv and Bai, the metabolites and gene were combined to establish a predicted network (Fig. 7). The 13 metabolites were significantly expressed difference between Lv and Bai, including 10 down-regulated metabolites (Naringenin chalcone, Naringenin, Tricetin, Kaemferide and six Anthocyanins (Cyanidin O-acetylhexoside, Peonidin, Malvidin 3-O-glucoside, Malvidin 3, 5- diglucoside, Peonidin O-hexoside, and Rosinidin O-hexoside)) and three up-regulated metabolites ((+)-Catechin,(−)-Epicatechin, Gallocatechin). Seven PAL genes (Csa4G008250, Csa4G008760, Csa6G445240, Csa6G445760, Csa6G445770, Csa6G445780 and Csa6G446280) and one F3H (Csa6G108510) gene were up-regulated in the Bai compared to Lv. In addition, two structural genes 4CL (Csa2G433350 and Csa3G638510) were showed − 2.44- and − 1.70-fold decrement, and CHS (Csa3G600020) was − 1.14-fold down-regulation, this could largely explain the high accumulation of Naringenin chalcone and Naringenin in the Lv. Simultaneously, three UFGT genes (Csa3G172390, Csa6G109730 and Csa6G109750) showed − 5.30-, − 5.86- and − 1.65fold down-regulation in the Bai, which also supports six anthocyanins significantly down-regulated in Bai compared with Lv.
Combining omics analysis of diverse genetic resources provides crucial information in understanding molecular basis of plant traits such as fig fruit color , Lilium “Tiny Padhye” bicolor development , peanut resistance on salt stress . The cucumber shows a large variation in fruit skin colour, such as dark green, yellow, light green and milk white, these colours are characteristic of species or specific genotypes. In particular, the dark green and light green skin color cucumber cultivars have generated great interest in customer. In the study, we characterized two different cucumber on fruit skin color (Lv and Bai) using RNA-seq and metabolome. Lv exerted dark green with much more chlorophyll content and more closely arranged epidermal cells. Through analysis of different metabolites, flavones, flavanones, flavonols, and anthocyanins were mostly responsible for skin color differences. In addition, combining transcript level by RNA-seq, we found that several DEGs related to chlorophyll synthesis, anthocyanins synthesis and TFs were possibly involved in the color development.
Regulatory network of DEGs associated with chlorophyll synthesis pathway for skin color in lv and Bai
Chlorophyll is an important pigment for determined the skin color of many fruits. Chlorophyll synthesis has been well studied and important related genes for chlorophyll synthesis have been found in leave and fruits [8, 28]. Gang et al.  found that BpGLK1 the function for decreased chlorophyll content and defective chloroplast development by physiological and ultrastructural analysis. In addition, many key genes of coding enzymes were involved in chlorophyll synthesis pathway, such as HemA, HemB, chlH, chlM, por, NOL [3, 4]. For example, HemA, which is initiated enzyme for chlorophyll synthesis in plastid, catalyzes biosynthesis of 5-aminolevulinic acid from glutamyl-tRNA . The ChlH catalyzes protoporphyrin IX to form Mg-protoporphyrin IX. The magnesium protoporphyrin IX monomethyl ester formation was catalyzed magnesium protoporphyrin IX in chlorophyll synthesis pathway by ChlM . The por is an important enzyme that catalyzes protochlorophyllide to generate chlorophyllide, and this step is a critical intermediate step in converting chlorophyll . Here, 14 DEGs were identified in chlorophyll synthesis pathway. The expression of DEGs in synthesis of chlorophylls synthesis pathway, including one HemA, one HemB, one HemC, two HemE, one HemF, one chlH, one chlM, one chlE, one por, one CLH, two NOL, were down-regulated in Bai compared to Lv. These down-regulated expressions of many key genes involved in chlorophyll synthesis pathway may lead to inhibition of chlorophyll a and chlorophyll b synthesis. These results were consistent with higher accumulation of chlorophyll and more chloroplast in Lv than Bai.
Analysis of anthocyanins and flavonols synthesis for fruit skin color
Metabolites are the final products of cell biological regulation process  and metabolomic analysis enables us investigate the relationship between biological processes and plant characteristic  . The content of anthocyanins and flavonoids has crucial effect on fruit color and taste [22, 35]. The metabolome data combining with transcriptome profiling were discovered genes involved in anthocyanins and flavonols synthesis, thus searching for useful information to illustrate phenomenon of different color in cucumber fruit. Anthocyanins are the final products of the flavonoid biosynthetic pathways, our search showed many DEGs are differently expressed between Lv and Bai in this pathway, such as upstream 4CL, CHS, F3H and UFGT. Previous studies showed 4CL genes play an essential role at the divergence point flavonols aynthesis . The CHS has been found responsible for the anthocyanin biosynthesis during petal coloration in Malus crabapple . Our study identified two 4CL (Csa2G433350 and Csa3G638510) and CHS (Csa3G600020) genes were down-regulated in Bai compared with Lv, and two metabolites (Naringenin chalcone and Naringenin) also down-regulated in Bai. It indicated that CHS was significantly repressed in Bai, and lead to down-regulation of two important metabolites in anthocyanin synthesis. In addition, we detected six types of anthocyanins have differently expressed between Bai and Lv. In anthocyanins biosynthesis, the glycosyl is a crucial progress, which catalyzed by UFGT in Arabidopsis . The UFGT expression was associated with anthocyanin accumulation in different plant [39, 40]. Our results showed that three UFGT expressions are suppressed in Bai, it maybe explain six types of anthocyanins down-regulation in Bai compared to Lv. Other searcher found that the Cyanidin-3-O-rhamnoglucoside, one type of anthocyanins is main anthocyanin and played an important role in skin of figs [41, 42]., while cyanidin-3-O-rhamnoglucoside was not detected in our data, indicating it might be not main anthocyanin in cucumber fruit skin.
Analysis of TFs involved in biosynthesis of anthocyanin in lv fruit skin
Anthocyanins and flavonoid synthesis are regulated by several structural genes and TFs such as MYB, bHLH and WDR proteins. The bHLH proteins can interact with R2R3-MYBs from various subgroups, and form ternary complexes with WDR. The MBW (MYB–bHLH–WDR) complexes participated in flavonols, anthocyanins, and proanthocyanidins (PAs) biosynthesis pathway [43,44,45]. Among these, MYB as major determinant element for anthocyanin accumulation regulation, could activate some pivotal anthocyanin biosynthetic genes by interacting with bHLH respectively [46, 47]. Ectopic overexpression of pear PyMYB10 in Arabidopsis contributed to its pigmentation in immature seeds, indicating PyMYB10 as positive factor in regulating anthocyanin accumulation . Overexpression of peach PpMYB10.1 in tobacco could increase the expression of UFGT, leading to higher anthocyanin accumulation and deeper red flowers in transgenic tobacco . Similarly, MYB could regulate anthocyanin biosynthesis by regulating the expression of UFGT in grape  and apple . In our research, 16 MYB TFs were detected by transcriptome, and expression levels of eight MYBs were up-regulated in fruit skin of Lv compared with Bai, indicating MYBs in Lv contributed the expression of related genes involved in anthocyanin synthesis.
The bHLH played an important role in anthocyanin synthesis by forming a complex with MYBs . Overexpression of SlPRE2, an atypical bHLH, accelerated seedling morphogenesis and produced yellowing ripen fruits with reduced chlorophyll and carotenoid in tomato fruit . Overexpressing Arabidopsis GLABRA3 (bHLH) exhibited higher anthocyanin accumulation than control sample in tomato fruit . In this study, 11 bHLHs were up-regulated in Lv fruit skin, while seven bHLHs were significantly down-regulated compared with Bai, suggesting bHLHs function as different roles in biosynthesis of anthocyanin.
Overall, the regulation mechanism of fruit skin color on cucumber was firstly carried out by metabolome and RNA-Seq. The content of chlorophyll a, chlorophyll b and carotenoid were higher in Lv than Bai, and cytological observation showed more chloroplast existed in Lv. Crucial anthocyanins and flavonols responsible for fruit skin color development showed significantly different expression between two cucumber genotypes by metabolome. Several genes, especially por and NOL, CHS and UFGT, which play important roles in chlorophyll synthesis and anthocyanins biosynthesis pathway, respectively, were differently expressed between Bai and Lv fruit skin. Taken together, these different metabolites and genes identified in our study provide an important metabolic and functional role for chlorophyll synthesis and anthocyanins biosynthesis pathway in cucumber skin color.
Plant materials and growth conditions
Two cucumber high inbred lines (Lv and Bai) were used in this study, and were inbred line selected by our research group after multi-generation self-crossing. Lv and Bai were both South China type variety with contrasting differences in fruit skin color. Seeds were germinated on culture dish in a dark environment. Then, the seedlings were grown in a culture room under 14 h/10 h with 28 °C/18 °C in day/night. When plants were grown to two true leaf stages, and were transferred to the open field in Baiyun Area, Guangzhou City, China.
Analysis of chlorophyll and carotenoid content in fruit skin between Bai and lv
Chlorophyll and carotenoid content of fruit skin from Lv and Bai were measured on the basis of the procedure described by Xie et al. (2019) . Approximately, 0.2 g fruit skin were placed in 5 ml solution (9:1 = acetone: 0.1 M NH4OH).The samples were centrifuged at 3000 r for 20 min, and supernatants were collected. The same process was repeated thrice and the supernatants were collected using hexane. Finally, the mixed supernatant was measured by spectrophotometer at the absorption wavelengths of 663 nm and 645 nm (Beckman Coulter DU-800, CITY, USA). The measurements were performed with biological replicates.
Scanning and transmission electron microscopy
After cucumber fruit skin was air-dried, the epidermis cells were visualized under a HITACHI SU8020 variable pressure SEM (Hitachi, Japan). For TEM assay, fruit skin were cut into small pieces, and were collected for fixation, and the process was performed as according to Wang et al. (2019) .
Metabolite profiling was performed using a widely targeted metabolome method by Wuhan Metware Biotechnology Co., Ltd. (Wuhan, China) (http://www.metware.cn/). Freezing-dried fruit skin was crushed into powder using a mixer mill (MM 400, Retsch). The fruit skin (1 cm wide and 0.2 cm thick along the fruit lengthwise) were sampled 10–15 days after female flowers open, and three replicates each of Lv and Bai. A total of 100 mg powder was extracted overnight at 4 °C with 1.0 ml 70% aqueous methanol, then centrifuged at 10, 000 g for 10 min. After that, these extracts were absorbed, filtrated, and analyzed by an LC-ESI-MS/MS system. Analytical conditions were based on the procedures as described in Wang et al. (2017) . Quantification of metabolites was carried out using a MRM method . Metabolites with significant differences in content were set with thresholds of variable importance in projection (VIP) ≥1 and fold change ≥2 or ≤ 0.5 .
The fruit skin (1 cm wide and 0.2 cm thick along the fruit skin lengthwise in the middle part) were sampled 10–15 days after female flowers open. A total of twelve samples (three replicates each of Lv and Bai) were prepared for RNA extraction based on the instruction of TRIZOL reagent (TaKaRa, Japan). RNA was purified and concentrated using an RNeasy MinElute clean up kit (Qiagen, Germany) after RNA extraction. Then, about 2.5 μg RNA from each sample was prepared for constructing sequencing libraries and the library quality was detected by Agilent Bioanalyzer 2100 system. The library preparations were sequenced on Illumina Hiseq2500 platform and 125/150 bp paired-end reads were generated. Index of the reference genome was built using Bowtie v2.2.3 and paired-end clean reads were aligned to the reference genome using TopHat v2.0.12 .
Gene expression level was analysis by FPKM (fragments per kilobase per million reads) method . The FPKM of genes were calculated by Cuffquant and cuffnorm (v2.2.1) (v2.2.1) . DESeq2 was used to identify DEGs according to the two criteria(fold change ≥2 or ≤ 0.5and q ≤ 0.01). WEGO software and KEGG database were employed to GO enrichment and bigological pathway enrichment, respectively [59, 60].
Quantitative real-time PCR (qRT-PCR) validation
The qRT-PCR reaction was performed on ABI PRISM 7900HT machine (Applied Biosystems, USA) by using the SYBR Premix Ex Taq Kit (TaKaRa, Japan), and qRT-PCR reaction process was performed according to Wang et al. (2019) . All primers used in qRT-PCR were listed in Table S5.
Availability of data and materials
All the sequencing data were submitted to NCBI Sequence Read Archive database under accession number PRJNA647135.
Quantitative real-Time PCR
Differentially expressed genes
Transmission electron microscopy
Scanning electron microscope
Total ions current
Multiple reactions monitoring
Kyoto Encyclopedia of Genes and Genomes
Principal component analysis
Coproporphyrinogen III oxidase
Magnesium chelatase subunit H
Magnesium-protoporphyrin IX monomethyl ester
Divinyl chlorophyllide a 8-vinyl-reductase
Chlorophyllide a oxygenase
Chlorophyll/bacteriochlorophyll a synthase
Chlorophyll(ide) b reductase
7-hydroxymethyl chlorophyll a reductase
4-coumarate: CoA ligase
UDP glucose-flavonoid 3-O-glcosyl-transferase
Rosianskey Y, Dahan Y, Yadav S, Freiman ZE, Milo-Cochavi S, Kerem Z, Eyal Y, Flaishman MA. Chlorophyll metabolism in pollinated vs. parthenocarpic fig fruits throughout development and ripening. Planta. 2016;244(2):491–504.
Kayesh E, Shangguan L, Korir NK, Sun X, Bilkish N, Zhang Y, Han J, Song C, Cheng Z, Fang J. Fruit skin color and the role of anthocyanin. Acta Physiol Plant. 2013;35(10):2879–90.
Tanaka R, Tanaka A. Tetrapyrrole biosynthesis in higher plants. Annu Rev Plant Biol. 2007;58:321–46.
Tanaka A, Tanaka R. Chlorophyll metabolism. Curr Opin Plant Biol. 2006;9(3):248–55.
Lee SB, Kim JE, Kim HT, Lee G, Kim B, Lee JM. Genetic mapping of the c1 locus by GBS-based BSA-seq revealed Pseudo-Response Regulator 2 as a candidate gene controlling pepper fruit color. Theor Appl Genet. 2020;133(6):1897–910.
Xie J, Yao S, Ming J, Deng L, Zeng K. Variations in chlorophyll and carotenoid contents and expression of genes involved in pigment metabolism response to oleocellosis in citrus fruits. Food Chem. 2019;272:49–57.
Meng L, Fan Z, Zhang Q, Wang C, Gao Y, Deng Y, Zhu B, Zhu H, Chen J, Shan W, et al. BEL1-LIKE HOMEODOMAIN 11 regulates chloroplast development and chlorophyll synthesis in tomato fruit. Plant J. 2018;94(6):1126–40.
Lai B, Hu B, Qin Y, Zhao J, Wang H, Hu G. Transcriptomic analysis of Litchi chinensis pericarp during maturation with a focus on chlorophyll degradation and flavonoid biosynthesis. BMC Genomics. 2015;16(1):225.
Falcone FM, Rius SP, Casati P. Flavonoids: biosynthesis, biological functions, and biotechnological applications. Front Plant Sci. 2012;3:222.
Martens S, Preuss A, Matern U. Multifunctional flavonoid dioxygenases: flavonol and anthocyanin biosynthesis in Arabidopsis thaliana L. Phytochemistry. 2010;71(10):1040–9.
Pelletier MK, Murrell JR, Shirley BW. Characterization of flavonol synthase and leucoanthocyanidin dioxygenase genes in Arabidopsis (further evidence for differential regulation of “early” and “late” genes). Plant Physiol. 1997;113(4):1437–45.
Boss PK, Davies C, Robinson SP. Analysis of the expression of anthocyanin pathway genes in developing vitis vinifera L. cv shiraz grape berries and the implications for pathway regulation. Plant Physiol. 1996;111(4):1059–66.
Winkel-Shirley B. Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. Plant Physiol. 2001;126(2):485–93.
Saslowsky D, Winkel-Shirley B. Localization of flavonoid enzymes in Arabidopsis roots. Plant J. 2001;27(1):37–48.
Grotewold E. The challenges of moving chemicals within and out of cells: insights into the transport of plant natural products. Planta. 2004;219(5):906–9.
Pierce LK, Wehner TC. Review of genes and linkage groups in cucumber. Hortscience. 1990;25(6):605–15.
Call A, Wehner T. Gene list 2010 for cucumber. Cucurbit Genet Cooperative Rep. 2010;33-34:69–103.
Liu H, Meng H, Pan Y, Liang X, Jiao J, Li Y, Chen S, Cheng Z. Fine genetic mapping of the white immature fruit color gene w to a 33.0-kb region in cucumber (Cucumis sativus L.). Theor Appl Genet. 2015;128(12):2375–85.
Moreno-Risueno MA, Busch W, Benfey PN. Omics meet networks - using systems approaches to infer regulatory networks in plants. Curr Opin Plant Biol. 2010;13(2):126–31.
Fridman E, Pichersky E. Metabolomics, genomics, proteomics, and the identification of enzymes and their substrates and products. Curr Opin Plant Biol. 2005;8(3):242–8.
Li Y, Fang J, Qi X, Lin M, Zhong Y, Sun L, Cui W. Combined analysis of the fruit metabolome and transcriptome reveals candidate genes involved in flavonoid biosynthesis in Actinidia arguta. Int J Mol Sci. 2018;19(5):1471.
Wang Z, Cui Y, Vainstein A, Chen S, Ma H. Regulation of fig (Ficus carica L.) fruit color: Metabolomic and transcriptomic analyses of the flavonoid biosynthetic pathway. Front. Plant Sci. 2017;8:1990.
Xu L, Yang P, Feng Y, Xu H, Cao Y, Tang Y, Yuan S, Liu X, Ming J. Spatiotemporal transcriptome analysis provides insights into bicolor tepal development in lilium “tiny Padhye”. Front Plant Sci. 2017;8:398.
Yonekura-Sakakibara K, Tohge T, Matsuda F, Nakabayashi R, Takayama H, Niida R, Watanabe-Takahashi A, Inoue E, Saito K. Comprehensive flavonol profiling and transcriptome coexpression analysis leading to decoding gene-metabolite correlations in Arabidopsis. Plant Cell. 2008;20(8):2160.
Saito K, Hirai MY, Yonekura-Sakakibara K. Decoding genes with coexpression networks and metabolomics – ‘majority report by precogs’. Trends Plant Sci. 2008;13(1):36–43.
Tang HY, Dong X, Wang JK, Xia JH, Xie F, Zhang Y, Yao X, Xu YJ, Wang ZJ. Fine mapping and candidate gene prediction for white immature fruit skin in cucumber (Cucumis sativus L.). Int J Mol Sci. 2018;19(5):1493.
Cui F, Sui N, Duan G, Liu Y, Han Y, Liu S, Wan S, Li G. Identification of metabolites and transcripts involved in salt stress and recovery in peanut. Front Plant Sci. 2018;9:217.
Wen C, Lin S, Chu F. Transcriptome analysis of a subtropical deciduous tree: autumn leaf senescence gene expression profile of Formosan gum. Plant Cell Physiol. 2015;56(1):163–74.
Gang H, Li R, Zhao Y, Liu G, Chen S, Jiang J. Loss of GLK1 transcription factor function reveals new insights in chlorophyll biosynthesis and chloroplast development. J Exp Bot. 2019;70(12):3125–38.
Eckhardt U, Grimm B, Hortensteiner S. Recent advances in chlorophyll biosynthesis and breakdown in higher plants. Plant Mol Biol. 2004;56(1):1–14.
Wang Z, Hong X, Hu K, Wang Y, Wang X, Du S, Li Y, Hu D, Cheng K, An B, et al. Impaired magnesium protoporphyrin IX methyltransferase (ChlM) impedes chlorophyll synthesis and plant growth in rice. Front Plant Sci. 2017;8:1694.
Kwon C, Kim S, Song G, Kim D, Paek N. Two NADPH: protochlorophyllide oxidoreductase (POR) isoforms play distinct roles in environmental adaptation in rice. Rice. 2017;10(1):1.
Fiehn O. Metabolomics--the link between genotypes and phenotypes. Plant Mol Biol. 2002;48(1–2):155–71.
Yang D, Zhang J, Li M, Shi L. Metabolomics analysis reveals the salt-tolerant mechanism in Glycine soja. J Plant Growth Regul. 2017;36(2):460–71.
Wu SB, Dastmalchi K, Long C, Kennelly EJ. Metabolite profiling of Jaboticaba (Myrciaria cauliflora) and other dark-colored fruit juices. J Agric Food Chem. 2012;60(30):7513–25.
Sun H, Li Y, Feng S, Zou W, Guo K, Fan C, Si S, Peng L. Analysis of five rice 4-coumarate:coenzyme a ligase enzyme activity and stress response for potential roles in lignin and flavonoid biosynthesis in rice. Biochem Bioph Res Co. 2013;430(3):1151–6.
Tai D, Tian J, Zhang J, Song T, Yao Y. A Malus crabapple chalcone synthase gene, McCHS, regulates red petal color and flavonoid biosynthesis. PLoS One. 2014;9(10):e110570.
Saito K, Yonekura-Sakakibara K, Nakabayashi R, Higashi Y, Yamazaki M, Tohge T, Fernie AR. The flavonoid biosynthetic pathway in Arabidopsis: structural and genetic diversity. Plant Physiol Biochem. 2013;72:21–34.
Wang H, Wang C, Fan W, Yang J, Appelhagen I, Wu Y, Zhang P. A novel glycosyltransferase catalyses the transfer of glucose to glucosylated anthocyanins in purple sweet potato. J Exp Bot. 2018;69(22):5444–59.
Song C, Zhao S, Hong X, Liu J, Schulenburg K, Schwab W. A UDP-glucosyltransferase functions in both acylphloroglucinol glucoside and anthocyanin biosynthesis in strawberry (Fragaria x ananassa). Plant J. 2016;85(6):730–42.
Ercisli S, Tosun M, Karlidag H, Dzubur A, Hadziabulic S, Aliman Y. Color and antioxidant characteristics of some fresh fig (Ficus carica L.) genotypes from northeastern Turkey. Plant Foods Hum Nutr. 2012;67(3):271–6.
Solomon A, Golubowicz S, Yablowicz Z, Grossman S, Bergman M, Gottlieb HE, Altman A, Kerem Z, Flaishman MA. Antioxidant activities and anthocyanin content of fresh fruits of common fig (Ficus carica L.). J Agric Food Chem. 2006;54(20):7717–23.
Xu W, Dubos C, Lepiniec L. Transcriptional control of flavonoid biosynthesis by MYB-bHLH-WDR complexes. Trends Plant Sci. 2015;20(3):176–85.
Xu W, Grain D, Bobet S, Le Gourrierec J, Thevenin J, Kelemen Z, Lepiniec L, Dubos C. Complexity and robustness of the flavonoid transcriptional regulatory network revealed by comprehensive analyses of MYB-bHLH-WDR complexes and their targets in Arabidopsis seed. New Phytol. 2014;202(1):132–44.
Koes R, Verweij W, Quattrocchio F. Flavonoids: a colorful model for the regulation and evolution of biochemical pathways. Trends Plant Sci. 2005;10(5):236–42.
Petroni K, Tonelli C. Recent advances on the regulation of anthocyanin synthesis in reproductive organs. Plant Sci. 2011;181(3):219–29.
Allan AC, Hellens RP, Laing WA. MYB transcription factors that colour on fruit. Trends Plant Sci. 2008;13(3):99–102.
Feng S, Wang Y, Yang S, Xu Y, Chen X. Anthocyanin biosynthesis in pears is regulated by a R2R3-MYB transcription factor PyMYB10. Planta. 2010;232(1):245–55.
Tuan PA, Bai S, Yaegaki H, Tamura T, Hihara S, Moriguchi T, Oda K. The crucial role of PpMYB10.1 in anthocyanin accumulation in peach and relationships between its allelic type and skin color phenotype. BMC Plant Biol. 2015;15:280.
Koyama K, Numata M, Nakajima I, Goto-Yamamoto N, Matsumura H, Tanaka N. Functional characterization of a new grapevine MYB transcription factor and regulation of proanthocyanidin biosynthesis in grapes. J Exp Bot. 2014;65(15):4433–49.
Chagne D, Lin-Wang K, Espley RV, Volz RK, How NM, Rouse S, Brendolise C, Carlisle CM, Kumar S, De Silva N, et al. An ancient duplication of apple MYB transcription factors is responsible for novel red fruit-flesh phenotypes. Plant Physiol. 2013;161(1):225–39.
Zhu Z, Chen G, Guo X, Yin W, Yu X, Hu J, Hu Z. Overexpression of SlPRE2, an atypical bHLH transcription factor, affects plant morphology and fruit pigment accumulation in tomato. Sci Rep. 2017;7(1):5786.
Nukumizu Y, Wada T, Tominaga-Wada R. Tomato (Solanum lycopersicum) homologs of TRIPTYCHON (SlTRY) and GLABRA3 (SlGL3) are involved in anthocyanin accumulation. Plant Signal Behav. 2013;8(7):e24575.
Wang M, Jiang B, Liu W, Lin Y, Liang Z, He X, Peng Q. Transcriptome analyses provide novel Insights into heat stress responses in Chieh-Qua (Benincasa hispida Cogn. var. Chieh-Qua How). Int J Mol Sci. 2019;20(4):883.
Yun Y, Liang F, Deng B, Lai G, Vicente Gonçalves CM, Lu H, Yan J, Huang X, Yi L, Liang Y. Informative metabolites identification by variable importance analysis based on random variable combination. Metabolomics. 2015;11(6):1539–51.
Fehlmann T, Reinheimer S, Geng C, Su X, Drmanac S, Alexeev A, Zhang C, Backes C, Ludwig N, Hart M, et al. cPAS-based sequencing on the BGISEQ-500 to explore small non-coding RNAs. Clin Epigenetics. 2016;8:123.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.
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:W293–7.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:D480–4.
The authors thank other lab members for assistance.
This work was supported by the National Key Research and Development Program (2018YFD0100700), the Natural Science Foundation of Guangdong Province (2018A030310196), Special Fund for Agro-scientific Research in the Public Interest (201503110–07), the Presidential foundation of Guangdong Academy of Agricultural Sciences (201813), Special fund for scientific innovation strategy-construction of high level Academy of Agriculture Science (R2019YJ-YB3004), Science and Technology Program of Guangdong (2020B020220001), Discipline team building projects of Guangdong Academy of Agricultural Sciences in the 13th Five-Year Period, Special funds for agricultural industry development of Guangzhou (1710023), the Key laboratory open fund project of Vegetable Institution of GAAS (201701). The funders had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Single cell on average between Lv and Bai.
The correlation analysis and principal component analysis (PCA) in Lv and Bai fruit.
The number of DEGs belonging to different transcription factor families detected in Lv and Bai.
Relative expression of genes related to transcriptional factors. Data is presented as the mean ± standard deviation (n = 9).
Differentially metabolites in the skin of Lv and Bai fruit.
Overview of reads from Lv and Bai fruit skin by RNA-seq.
Differentially expressed genes between Lv and Bai fruit skin.
GO enrichment analysis for DEG in Lv and Bai fruit skin.
List of primer of qRT-PCR.
About this article
Cite this article
Wang, M., Chen, L., Liang, Z. et al. Metabolome and transcriptome analyses reveal chlorophyll and anthocyanin metabolism pathway associated with cucumber fruit skin color. BMC Plant Biol 20, 386 (2020). https://doi.org/10.1186/s12870-020-02597-9
- Cucumis sativus L.