- Research article
- Open Access
Study of 'Redhaven' peach and its white-fleshed mutant suggests a key role of CCD4 carotenoid dioxygenase in carotenoid and norisoprenoid volatile metabolism
BMC Plant Biology volume 11, Article number: 24 (2011)
Carotenoids are plant metabolites which are not only essential in photosynthesis but also important quality factors in determining the pigmentation and aroma of flowers and fruits. To investigate the regulation of carotenoid metabolism, as related to norisoprenoids and other volatile compounds in peach (Prunus persica L. Batsch.), and the role of carotenoid dioxygenases in determining differences in flesh color phenotype and volatile composition, the expression patterns of relevant carotenoid genes and metabolites were studied during fruit development along with volatile compound content. Two contrasted cultivars, the yellow-fleshed 'Redhaven' (RH) and its white-fleshed mutant 'Redhaven Bianca' (RHB) were examined.
The two genotypes displayed marked differences in the accumulation of carotenoid pigments in mesocarp tissues. Lower carotenoid levels and higher levels of norisoprenoid volatiles were observed in RHB, which might be explained by differential activity of carotenoid cleavage dioxygenase (CCD) enzymes. In fact, the ccd4 transcript levels were dramatically higher at late ripening stages in RHB with respect to RH. The two genotypes also showed differences in the expression patterns of several carotenoid and isoprenoid transcripts, compatible with a feed-back regulation of these transcripts. Abamine SG - an inhibitor of CCD enzymes - decreased the levels of both isoprenoid and non-isoprenoid volatiles in RHB fruits, indicating a complex regulation of volatile production.
Differential expression of ccd4 is likely to be the major determinant in the accumulation of carotenoids and carotenoid-derived volatiles in peach fruit flesh. More in general, dioxygenases appear to be key factors controlling volatile composition in peach fruit, since abamine SG-treated 'Redhaven Bianca' fruits had strongly reduced levels of norisoprenoids and other volatile classes. Comparative functional studies of peach carotenoid cleavage enzymes are required to fully elucidate their role in peach fruit pigmentation and aroma.
Among Rosaceae, peach (Prunus persica L. Batsch) is an appealing model crop, because of its economical value, small genome, rapid generation time and several Mendelian traits (i.e. flesh/leaf/flower color, smooth/fuzzy skin, clingstone/freestone, normal/dwarf growth habit) still to be functionally characterized [1, 2]. Peaches are appreciated for their visual, nutritional and organoleptic features, partially contributed by carotenoids, sugars, acids and volatile organic compounds (VOCs), which vary as a function of genetic, developmental and post-harvest factors [[3–5] and references therein]. In particular, carotenoid accumulation in the mesocarp determines the difference between yellow- and white-fleshed genotypes, the latter being generally characterized by a peculiar and more intense aroma. Flesh color is a Mendelian trait (white genotype dominant over yellow ), associated with the Y locus that has been mapped on the linkage group 1 of the Prunus map  but which has not been yet functionally characterized from the molecular or enzymatic point of view. Natural mutations, originating flesh color chimera with irregular yellow and white distribution, have long been observed in peach .
Carotenoids are a widespread class of compounds having important functions across living organisms, whose accumulation shows striking phylum- and genotype-specific regulation . Following the formation of the first carotenoid phytoene from the general isoprenoid pathway, the pathway bifurcates after lycopene with respect to the ring type, giving rise to carotenes and xanthophylls with either β-β or β-ε rings (Figure 1, Additional File 1). In addition to their roles in plants as photosynthetic accessory pigments and colorants, carotenoids are also precursors to norisoprenoids (also called apocarotenoids). Norisoprenoids are commonly found in flowers, fruits, and leaves of many plants  and possess aromatic properties together with low odor thresholds (e.g., β-ionone), thus having a strong impact on fruit and flower aroma even at low levels . An increasing number of dioxygenase enzymes that specifically cleave carotenoid compounds to form volatile norisoprenoids, abscisic acid (ABA) and regulators of plant growth and development has been characterized. These enzymes have been referred to as carotenoid cleavage dioxygenases (CCDs) and 9-cis-epoxycarotenoid dioxygenases (NCEDs)  and represent a plant multienzyme family: Arabidopsis has nine CCD/NCED members, of which four have been classified as CCDs (AtCCD1, AtCCD4, AtCCD7 and AtCCD8) and the remaining as ABA-related NCEDs . Functional analysis of CCD enzymes determined that CCD7 and CCD8 are mostly related to the synthesis of norisoprenoid (apocarotenoid) plant hormones, while CCD1 and CCD4 are preferentially involved in volatile production, by using different carotenoid substrates with variable specificity and cleavage site, which probably contributes to the diversity of norisoprenoids found in nature [[14–17] and references therein]. The synthesis of β-ionone, geranylacetone and 6-methyl-5-hepten-2-one in tomato fruits increases 10-20 fold during fruit ripening and these compounds were produced by the activity of ccd products . Silencing of ccd genes resulted in a significant decrease of the β-ionone content of tomato ripe fruits and petunia flowers [18, 19], and increased pigmentation of potato tubers and Chrysanthemum flowers [20, 21]. CCDs were also implied to be involved in the formation of important apocarotenoid aroma compounds in melon fruit . Furthermore, comparative genetics studies have indicated that carotenoid pigmentation patterns have profound effects on the norisoprenoid and monoterpene aroma volatile compositions of tomato and watermelon fruits . Many norisoprenoids strongly contribute to peach fruit aroma, and their levels increase during fruit ripening . The partial purification and biochemical characterization of β-carotene degrading enzyme(s) from nectarine skin extracts was associated with the formation of C13 norisoprenoids . The study of two NCED-encoding genes from peach showed their differential expression, suggesting a functional relation to ABA formation during fruit ripening . The recent synthesis of specific carotenoid dioxygenase inhibitors [25, 26] enables to assess the role of such enzymes in vivo, not only in ABA biosynthesis but also in fruit VOC metabolism.
In order to improve our knowledge of carotenoid and VOC biosynthesis in peach fruit and determine the factor(s) controlling carotenoid accumulation in peach flesh, the two cultivars 'Redhaven' (RH; yellow-fleshed) and its white-fleshed bud sport mutant 'Redhaven Bianca' (RHB)  were investigated. Carotenoid accumulation, VOC content and transcript levels of relevant carotenoid biosynthetic genes were studied at five stages of fruit development (Figure 2). The effect of the carotenoid dioxygenase inhibitor, abamine SG, on fruit VOC composition at full ripening was also studied.
Fruit phenotype during ripening
At S1-S3 stages, whole fruits of RHB (Figure 2A) and RH (Figure 2B) look similar, while the final ripening stages Br and S4, carotenogenesis is well established in RH fruits, and differences in flesh color between yellow-fleshed RH and white-fleshed RHB fruits become dramatic. The yellow pigmentation visible along the suture at full ripening stage (Figure 2A, S4) indicates that RHB is a bud sport chimera mutant, in which the mutation does not involve the L-I apical cell layer, that originates the epidermis and several cells layers at the suture of the ovary wall . The flesh color phenotype has remained stable throughout several cycles of clonal propagation, pointing out the stability of the chimera in RHB (Liverani A., unpublished data). The mutation is transmitted to the progeny in a Mendelian fashion, and is associated to the Y locus controlling fruit flesh color. RHB is heterozygous for this locus (Yy), originating either 1:3 of yellow- to white-fleshed seedlings when selfed or 1:1 of yellow- to white-fleshed progenies when crossed with yellow (yy) peach accessions (Liverani A., unpublished data). The yellow strip is not observed in white-fleshed RHB progenies (Liverani A., unpublished data), because the gametes does not originate from the L1 cell layer but from the fully-mutated L2 layer. The major fruit quality traits and skin color parameters of the two cultivars were also measured, showing statistically significant differences only for soluble solids content (Additional File 2).
Carotenoid composition of RHB and RH fruits
In carotenoid-containing fruits, the massive biosynthesis of such compounds is generally associated with late ripening stages and plastid transition from chloroplasts into chromoplasts. At early ripening stages S1 and S2, fruits of RHB and RH had similar total carotenoid levels and accumulated only a few carotenoid compounds, dominated by the presence of lutein and β-carotene (Figure 3; Table 1). From the S3 stage, RH mesocarp accumulated increasing amounts of carotenoids that peaked at the S4 stage to provide the solid yellow flesh color, while carotenoid content in RHB flesh remained low (Figure 3). Detailed HPLC analysis revealed the presence of specific carotenoid compounds in the two genotypes (Table 1), some of which rather uncommon and present in RH only, whose main chemical structures are illustrated (Additional File 1). Until stage S3 (RH) and Br (RHB), lutein and its (Z)-isomers were the major carotenoids in fruits of both genotypes, accounting for over 50% of the total carotenoid pool. Other major carotenoids at early stages were β-carotene, relatively abundant in fruits of both cultivars, and neochrome epimers, which accumulated only in RH fruits. From stage S3, not only did RH fruits have higher carotenoid levels, but also a range of carotenoid compounds much wider than RHB (Table 1). At full ripening S4 stage, xanthophylls made up the majority of carotenoids - zeaxanthin was the main carotenoid in RHB, while antheraxanthin, luteoxanthin and zeaxanthin were the most abundant compounds in RH fruits (Table 1).
Gene expression analyses
Relative transcript levels of three genes involved in isoprenoid metabolism [1-deoxy-d-xylulose 5-phosphate synthase (dxs), 4-(cytidine 5'-diphospho)-2-C-methyl-d-erythritol kinase (cmk) and 4-hydroxy-3-methylbut-2-enyl diphosphate reductase (hdr)] and twelve genes involved in carotenoid biosynthesis and cleavage [phytoene synthase (psy), phytoene desaturase (pds), ζ-carotene desaturase (zds), lycopene β-cyclase (lcy-b), lycopene ε-cyclase (lcy-e), carotene β-hydroxylase (chy-b), carotene ε-hydroxylase (chy-e), zeaxanthin epoxidase (zep), two carotenoid cleavage dioxygenases (ccd1 and ccd4), and two 9-cis-epoxycarotenoid dioxygenases (nced1 and nced2)] (Figure 1) were measured in RH and RHB mesocarp at S1, S2, S3, Br and, S4 stages by reverse transcription quantitative real-time PCR (RT-qPCR).
The isoprenoid pathway genes, dxs and cmk, had very low transcript levels throughout fruit development in both cultivars. hdr showed a sharp peak of expression at the S2 stage which declined at later stages in the yellow-fleshed RH, while its expression remained high in RHB (Figure 4A). Similarly, early carotenoid pathway genes psy and zds showed a peak at the S3 stage in RH, while in RHB these transcripts showed a constant increase until the S4 stage (Figure 4B). Among later pathway genes (lcy-b , lcy-e, chy-b, chy-e and zep), only chy-b was strongly up-regulated in RHB (Figure 4C), while, ccd and nced expression was generally low in both genotypes, with the exception of ccd4, which was significantly up-regulated in RHB at late ripening stages, its transcript level being 13-fold higher than that in RH at the S4 stage (Figure 4D).
Hierarchical clustering analysis (HCA) of gene expression data clustered the ripening stages consistently with their chronological order in joint analysis of both genotypes (Additional File 3A) and in RH alone (Additional File 3B). In RH, a clear co-regulation of genes encoding enzymes closely positioned in the pathway (dxs and cmk; psy, pds and zds; lcy-b and lcy-e; nced1, ccd1, ccd4 and, surprisingly, chy-b) was observed (Additional File 3C). Instead, in the RHB mutant the majority of these co-regulation clusters was broken, with the exception of ccd4 and chy-b genes which remained co-regulated (Additional File 3C).
Levels of different VOCs were studied in RHB and RH during fruit ripening by GC-MS. In total, 41 VOCs were detected, assigned to aromatic and branched chain amino acid-related, fatty acid-related, furan-related, lactone, monoterpene and norisoprenoid classes, quantified and underwent further analyses (Table 2).
The two genotypes had a similar, ripening-associated accumulation of total VOCs starting from the S3 stage, while early S1 and S2 stages were characterized by very low volatile content (Additional File 4). Detailed analysis pointed out differences in the accumulation of the distinct VOC pools in the two genotypes (Figure 5). Furan-related compounds accumulated at the highest levels in both genotypes, followed by norisoprenoids and fatty acid-related compounds, whose maximum levels were about 5-fold lower than those of the furans (compare Figures 5D, C and 5H). The other classes accumulated at lower absolute levels, with maxima in the range of hundreds of μg/g fresh weight. The two genotypes displayed similar ripening-associated patterns for aromatic and branched chain amino acid-, fatty acid-, and furan-related classes, with a peak at the S3/Br stages and a more or less pronounced decline at later ripening stage(s) (Figures 5A-D). Total lactone and monoterpene contents displayed a different pattern, with a strong up-regulation only at final S4 ripening stage in the two cultivars (Figure 5E-F). At the S4 stage, all the six above-mentioned VOC classes had higher levels in RH fruits (Figure 5A-F).
A remarkable exception was the norisoprenoid pool, which accumulated in RHB fruits at levels higher than those of RH from S3 stage on. Norisoprenoid pattern in RHB fruits peaked at Br and was constant through S4 stage, while in RH fruits it displayed a linear increase from S3 stage (Figure 5H). In particular, the three identified norisoprenoids 3-hydroxy-5,6-epoxide-β-ionone, 3-hydroxy-β-damascone and 4-hydroxy-3,5,6-trimethyl-4-(3-oxo-1-butenyl)-2-cyclo-hexen-1-one were responsible for the higher total norisoprenoid levels in fruits of RHB, with an almost 3-fold difference at the S4 stage (Figure 5 H1), when the typical floral scent of white-fleshed peach fruits reaches its maximum. At any ripening stage, the level of each identified norisoprenoid compound was always higher in the white-fleshed RHB than in RH (Additional File 5). Instead, the less prominent unidentified norisoprenoids had a similar accumulation pattern in the two genotypes, with higher levels in RH fruits (Figure 5 H2).
PCA was performed on the whole GC-MS dataset (41 major VOCs) to provide a more intuitive visualization of data and to discriminate the different ripening stages in the two varieties with respect to VOC composition. A preliminary PCA analysis was carried out including all five stages, and resulted in a poor separation of most samples (Figure 6A), with the exception of RHB-S1, RH-S2, RH-S3 and RH-Br. Principal components 1 and 2 explained 67% and 21% of the total variability, respectively (Figure 6A). A narrower analysis was then carried out by excluding the S1 and S2 samples, which allowed the complete discrimination of the six late ripening samples of both genotypes (Figure 6B). In this closer analysis, the new calculated principal components 1 and 2 accounted for 76% and 13% of the total variability, respectively (Figure 6B).
Effect of a carotenoid dioxygenase inhibitor on VOC production
The effect of abamine SG treatment was assessed in RHB ripe fruits, injected once a week from the S3 stage. As expected, abamine SG injection resulted in a drastic reduction of the levels of both identified and unknown norisoprenoids (Figure 7). Unexpectedly, this treatment also down-regulated the content of the other VOCs. The strongest reduction was observed for furan-related, monoterpene and lactone pools, while the total content of aromatic amino acid- and fatty acid-related compounds was the least affected by abamine SG treatment (Figure 7).
Peach fruits contain carotenoid compounds with significant antioxidant capacity and claimed beneficial health effects. The enzymatic cleavage of these compounds results in the production of volatile norisoprenoids (apocarotenoids). Our study investigated the expression of carotenoid genes, the carotenoid content and the volatile composition in the wild type yellow-fleshed RH and its white-fleshed RHB during fruit ripening.
Key regulatory steps and regulation mechanisms controlling isoprenoid and carotenoid flux in many species, also through biotechnology-based approaches to carotenoid manipulation, have been extensively reviewed [9, 29]. Differential accumulation of carotenoids in RH and RHB became evident at stage S3, and reached its maximum at stage S4, when RH fruits accumulated approximately 10-fold more carotenoids than RHB (Figure 3), composed mainly of β-ring carotenoids (Table 1). Accordingly, the two cultivars showed strikingly different developmental patterns of expression of several isoprenoid and carotenoid genes (Figure 4). However, several differential regulation events appeared to be the effect, rather than the cause, of the differential carotenoid accumulation. For instance, in RHB the up-regulation of hdr, psy, and zds genes at late (Br and S4) stages appears to be the result of a feed-back repression by either the lower carotenoid levels in this genotype or their cleavage products, acting on the transcript levels of these genes. Examples of feed-back regulation of carotenoid gene transcripts have been described in plants where carotenoid biosynthesis has been altered through the use of inhibitors or metabolic engineering [30–32]. A similar higher expression of early carotenoid genes was found in a white-fleshed apricot cultivar with respect to an orange-fleshed genotype , possibly suggesting common regulation mechanisms of early carotenoid gene expression based on feedback regulation mediated by carotenoid end products. The generally low expression at late Br/S4 stages of the studied carotenoid genes in the yellow-fleshed RH might point out a control of metabolite flux based on steady-state gene/enzyme expression rather than up- or down regulation of transcription.
On the other hand, in RHB fruits the strong up-regulation at late stages of chy-b and ccd4 genes is negatively correlated with the accumulation of β-ring carotenoids, and positively correlated with that of identified β-ring norisoprenoids. The biochemical function of the CHY-b and CCD4 gene products is compatible with the observed phenotype of RHB: CHY-b funnels carotenoids into the β-xanthophyll branch, producing substrates cleaved by carotenoid dioxygenases, including CCD4. Furthermore, chy-b is a negative regulatory step in several plant systems: inhibition of its expression through transgenosis or natural genetic variation results in higher β-carotenoid levels in potato tubers and maize endosperm [34, 35]. Similarly, the levels of ccd4 transcript negatively correlate with carotenoid levels (see below).
Aubert et al. [4, 5] described the presence of several C13 norisoprenoids, mainly derived from the cleavage of β-ring xanthophylls like in the present study, in yellow- white-fleshed nectarines. Unlike RH and RHB, the two cultivars were not isogenic, but similarly to the present case, the white-fleshed nectarine 'Vermeil' had higher norisoprenoid contents than the yellow-fleshed 'Springbright'. However, in two other studies, the β-ionone levels detected in a number of white-fleshed genotypes were lower than those of several unrelated yellow-fleshed accessions [36, 37]. Such examples show that flesh color per se is not sufficient to determine the levels of norisoprenoid VOCs, confirming that volatile composition strongly depends on genetic factors other than carotenoid levels. In our study, RHB ripe fruits had higher norisoprenoid content than that of RH, derived from carotenoid cleavage. As a likely consequence of the redirection of metabolite flux towards the synthesis of nor/isoprenoid compounds in RHB fruits, the levels of all other VOC pools were lower than those in RH fruits.
Among aroma-related carotenoid dioxygenases, both CCD1 and CCD4 enzymes cleave carotenoids at the 9,10 and 9',10' positions and have a major role in the formation of β-ionone and other fruit and flower norisoprenoids [18, 19, 38–41]: the β-carotene degradation displayed by yellow nectarine skin extracts  most likely corresponds to CCD1 and/or CCD4 activity. Compared with CCD1s, the CCD4 enzymes were characterized more recently from several crops [40–43], and shown to have a major impact on organ pigmentation: ccd4 expression was higher in white-fleshed potato and in white-flowered Chrysanthemum genotypes than in their respective yellow-pigmented counterparts, and RNAi-mediated knockout boosted carotenoid accumulation [20, 21]. The small CCD4 family contains at least two forms of genes with different structure, expression patterns and genome position, and their encoded enzymes show different substrate specificity [41, 42]. The existence of plastid target peptides and the demonstrated plastid localization of CCD4 enzymes  allow these enzymes direct access to the carotenoid substrates, suggesting that they start the carotenoid degradation and norisoprenoid synthesis pathway. On the other hand, the lack of correlation between the patterns of ccd1 expression (Figure 4D) and norisoprenoid content (Figure 5H, H1 and H2) in the flesh of both cultivars reflects a situation common to other fruits like grape, melon, and tomato [16, 18, 39]. This is a likely consequence of the different localization of CCD1s and their substrates, since carotenoids are accumulated in plastids but known CCD1s lack plastid transit peptides [14, 38], and/or substrate preference, since each carotenoid dioxygenase can accept different carotenoids. Our data suggest that because of their different subcellular localization, CCD1s only contribute to volatile production, while CCD4s are likely to control also carotenoid degradation. As to other dioxygenase-encoding genes studied, nced1 and nced2 expression patterns were up-regulated during ripening (Additional File 3, D), in agreement with the reported physiological increase of ABA content at late fruit ripening stages , with a profile similar to that of ethylene-related accS and accO genes (data not shown).
The transcriptional control of peach norisoprenoid content by specific dioxygenase genes is similar to that of fatty acid pathway-related genes and enzymes, whose expression patterns generally show a strong positive correlation with those of the corresponding volatiles [44, 45]. The strong down-regulation of norisoprenoids in ripe RHB fruits treated with abamine SG, a carotenoid dioxygenase inhibitor applied at late steps of fruit development, confirmed the hypothesis that their production is mediated by carotenoid cleavage dioxygenases. Although abamine SG was reported to be a specific carotenoid cleavage dioxygenase and ABA biosynthesis inhibitor, more potent than the previously synthesized abamine [25, 26], the authors did not investigate pathways other than that leading to ABA formation under stress conditions. Therefore, the observed unexpected reduction of other VOC classes could be due to the non-specific inhibition of other dioxygenase enzymes active in different aroma pathways, suggesting that abamine SG is a useful tool to investigate VOC metabolism.
If CCD4 is confirmed to be the major factor responsible for carotenoid degradation, then the mutation that generated the RHB cultivar is likely to be a gain-of-function mutation, restoring ccd4 function. This is compatible with the dominant nature of the white flesh phenotype over the yellow one. Since the mutation originating RHB is associated with the Y locus (cf unpublished data of progeny tests, Results section), a crucial step of future research is to determine whether the Y locus, mapped on LG1, is linked to the ccd4 gene. The recently released draft of the peach genome sequence (http://www.rosaceae.org/peach/genome) will enable future comprehensive research towards the functional characterization of the peach CCD/NCED family, required to confirm the presented data on ccd4 and to elucidate the function of each member in different peach genotypes.
This study presented a comprehensive molecular and biochemical research of the carotenoid and VOC metabolism in peach fruit, driven by developmental and genetic cues, and pointed out the central role of carotenoid cleavage dioxygenases, namely the product of ccd4, in flesh color and peach aroma formation. By taking advantage of a wild type-mutant system contrasted in fruit flesh color, we provided new information for understanding the mechanisms controlling carotenoid biosynthesis, flesh color and volatile content in peach, which could be useful for the optimization of these important fruit quality traits. The completion of peach genome will improve existing databases [, http://www.rosaceae.org] and boost genetic and molecular research to confirm the role and interactions of carotenoid cleavage dioxygenases with other factors controlling fruit pigmentation and aroma, stimulating new comprehensive research on peach fruit quality traits [47, 48].
Plant material and sampling
Fruits of cv RH and RHB were harvested from trees grown in an experimental field located near Forlì (Po Valley, Italy; 44.161° N, 12.088° E) at five different ripening stages according with the growth curve (Figure 2): S1 (about 35 days after pollination, dap), S2 (about 50 dap), S3 (about 90 dap), Breaker (Br; about 115 dap), and S4 (about 122 dap) (Figure 2). For molecular and biochemical analyses, two replicates of four representative fruits of each stage were sampled, peeled, cut into 0.5-cm slices and the mesocarp was immediately frozen in liquid nitrogen and stored at -80 °C. In RHB, the tissues near the suture were discarded since were not affected by the mutation originating the white-fleshed phenotype. For visual and organoleptic quality scoring (Additional File 2), S4 fruits underwent the following measurements immediately after harvest: flesh firmness (FTA penetrometer, Turroni, Italy, equipped with an 8-mm plunger tip), soluble solid content (Smart1 automatic refractometer; Atago Co., Tokyo, Japan), titratable acidity (titration of 10 ml peach juice with NaOH 0.2 M until pH 8.2), skin color (Chroma Meter CR-200 reflectance colorimeter; Minolta Camera Co., Osaka, Japan).
For abamine SG treatments, 100 μl of an abamine SG solution (0,1 mM in DMSO) were injected with a hypodermical needle in different points of RHB fruits once a week from S3 stage until S4. DMSO-injected fruits were taken as controls. Mesocarp samples were collected at S4 stage as previously described, and used for GC-MS analyses.
Experiments complied with the MIQE guidelines , and relevant information is contained in the manuscript. Total RNA (two independent replicates for each sample) was isolated from frozen fruit tissue as reported . First strand cDNA was synthesized from 1 μg total RNA in 30 μL with oligo-d(T)17 and Superscript III (Invitrogen, Milan, Italy), according to manufacturer's instructions. cDNA concentration in the RT mix was quantified using a ND-1000 UV spectrophotometer (Nanodrop, Wilmington, USA). First strand cDNA (10 ng) was used as template for RT-qPCR assays, carried out with primers for rps28 and carotenoid genes. The rps28 gene was chosen as reference gene, based on preliminary tests (data not shown). All reactions were performed using the Applied Biosystems 7900HT Real Time PCR system and Platinum® SYBR® Green RT-qPCR SuperMix-UDG plus ROX (Applera Italia, Monza, Italy), following the manufacturer's procedures. No template and no amplification controls were included in each experiment.
Three RT-qPCR runs were carried out for each cDNA and gene to serve as technical replicates. PCR conditions were: 5 min at 95°C, followed by 45 cycle at 95°C for 15 s and at 58°C for 60 s. At the end of the PCR, melting curve analysis was carried out to check for the presence of the correct amplicons only. Relative transcript abundance was quantified using the relative standard curve method described in the ABI PRISM 7900 HT manual, and the data was normalized against the quantity of the reference rps28 transcript. Serial 10-fold dilution of each gene fragment were used to calculate the standard curve and measure the amplification efficiency for each target and reference gene. Sequences of peach dxs, cmk, hdr, psy, pds, zds, lcy-b, lcy-e, chy-b, chy-e, zep, ccd1, nced1 and nced2 genes, and the reference gene rps28, were obtained from NCBI database and . A 1.84-kb putative peach ccd4 gene was retrieved from the recently published peach draft genome  (Phytozome v6.0 search tool at http://www.phytozome.org/peach) by searching with a Malus ccd4 sequence (EU327777). The identified peach ccd4 coding sequence share 80% identity with its Malus ortholog. All RT-qPCR primers for studied genes were designed with Primer Express (Applera Italia, Monza, Italy) software (Additional File 6).
Analysis of VOCs
Frozen fruit tissue (15 g) was ground in liquid N2 into a fine powder. The volatile components were extracted with 100 ml methyl tert-butyl ether (MTBE) by vigorous shaking on a shaker apparatus overnight. Iso-butyl benzene (5 μg) was added as internal standard. The upper MTBE layer was separated, dried with sodium sulfate, and concentrated to a volume of 0.5 ml under a N2 stream. Samples were kept at 4°C until analysis. A 1-μl aliquot of the concentrated MTBE extract was injected into a GC-MSD (splitless mode). Results were an average of 6-8 replicate measurements.
GC-MS analysis was carried out using an Agilent GC-MSD system (CA, USA) EI Scan, equipped with a RTx-5sil MS column (injector temperature: 250°C; detector temperature: 280°C, 70 eV). Oven temperature profile was: 50 °C (1 min), 5° C/min to 280° C, 280° (5 min). Gas flow: 0.8 ml/min. Mass range: 41 to 500 m/z. Compound identification was performed by comparing their relative retention indices and mass spectra with those of authentic standards or with those found in the literature and supplemented with NIST 98 and QuadLib 1607 GC-MS libraries. A mixture of straight-chain alkanes (C7-C23) was injected into the column under the above-mentioned conditions for retention index calculation .
Isolation and HPLC analysis of carotenoids
Plant material was homogenized in the presence of methanol and extracted three-times with MeOH. The methanolic extracts were combined and the carotenoid content of this solutions was transferred in a separatory funnel into toluene:hexane (1:1) mixture; evaporated and dissolved in diethyl-ether. After methanolic extractions, plant materials were then extracted once with diethyl-ether. The resulting extracts (ethereal solution of methanolic extracts + ethereal extract) were combined and saponified in a heterogeneous phase with 30% KOH/MeOH overnight (the 30% KOH/MeOH solution was layered on the lower part of the ethereal total extracts). After this process the reaction mixture was washed to alkali-free in a separatory funnel, dried over anhydrous Na2SO4, evaporated and the residues were dissolved in EtOH preparing the corresponding stock solutions for quantitative analysis. All of these procedures were carried out in semi-darkness and in N2-atmosphere [54–56].
The quantitative UV/VIS spectra of the corresponding stock solutions were recorded by a Jasco V-530 spectrophotometer. The total carotenoid content of the samples was calculated according to the Lambert-Beer's law (average molar extinction coefficient: 100,000; average molar mass of carotenoids: 600). The HPLC separation was performed with a Dionex Softron instrument (Germering, Germany) equipped with a Dionex P680 gradient pump, a Dionex PDA-100 detector and an end-capped C18 column (250 × 4.6 mm internal diameter; Merck LiChrospher 100 RP-18; 5 μm) thermostated at 22 ˚C. Elution was performed using 12% H2O/MeOH (A), MeOH (B), 50% Acetone/MeOH (C) at a flow rate of 1.25 ml min-1. The gradient program was: 100% A (0-2 min); 80% A and 20% B (2-10 min); 50% A and 50% B (10-18 min); 100% B (18-27 min); 100% C (27-34 min); 100% B (34-43 min); 100% A (43-56 min). Data acquisition was performed at 450 nm detection wavelength by Chromeleon 6.70 software.
GC-MS data were submitted to principal component analysis (PCA) using Systat 11 software (http://www.systat.com). The data set was made up of data from eight repetitions of each ripening stage of RH and RHB. The variable set was made of the major 41 volatile aroma compounds. PCA involves a mathematical procedure that transforms a number of possibly correlated variables into a smaller number of uncorrelated variables called principal components. The first principal component accounts for as much of the variability in the data as possible, and each succeeding component accounts for as much of the remaining variability as possible. Gene expression data were submitted to hierarchical clustering analysis (HCA) using Systat 11. The HCA data set was made up from the means from independent experiments, for every stage of fruit development, of both genotypes. Means from independent RT-qPCR experiments were subjected to one-way ANOVA and Tukey's pairwise comparisons, and fruit quality data to t test, carried out using PAST (http://folk.uio.no/ohammer/past/).
Folta KM, Gardiner SE: Genetics and genomics of Rosaceae. Plant Genetics and Genomics volume 6: Crops and Models., New York: Springer; 2009.
Bassi D, Monet R: Botany and taxonomy. The peach: botany, production and uses. Edited by: Layne DR, Bassi D. Wallingford: CAB International; 2008, 1-36. full_text.
Gil MI, Tomas-Barberan FA, Hess-Pierce B, Kader AA: Antioxidant capacities, phenolic compounds, carotenoids, and vitamin C contents of nectarine, peach, and plum cultivars from California. J Agric Food Chem. 2002, 50: 4976-4982. 10.1021/jf020136b.
Aubert C, Ambid C, Baumes R, Günata Z: Investigation of bound aroma constituents of yellow-fleshed nectarines (Prunus persica L. cv. Springbright). Changes in bound aroma profiles during maturation. J Agric Food Chem. 2003, 51: 6280-6286. 10.1021/jf034613h.
Aubert C, Günata Z, Ambid C, Baumes R: Changes in physicochemical characteristics and volatile constituents of yellow- and white-fleshed nectarines during maturation and artificial ripening. J Agric Food Chem. 2003, 51: 3083-3091. 10.1021/jf026153i.
Connors CH: Some notes on the inheritance of unit characters in the peach. Proc Am Soc Hortic Sci. 1920, 16: 24-36.
Bliss FA, Arulsekar S, Foolad MR, Becerra V, Gillen AM, Warburton ML, Dandekar AM, Kocsisne GM, Mydin KK: An expanded genetic linkage map of Prunus based on an interspecific cross between almond and peach. Genome. 2002, 45: 520-529. 10.1139/g02-011.
Yeager AF, Meader EM: A flesh-color chimera in the peach. J Hered. 1956, 47: 77-78.
Cazzonelli CI, Pogson BJ: Source to sink: regulation of carotenoid biosynthesis in plants. Trends Plant Sci. 2010, 15: 266-274. 10.1016/j.tplants.2010.02.003.
Winterhalter P, Rouseff RL: Carotenoid-derived aroma compounds. ACS Symposium series 802. Washington DC: American Chemical Society; 2002.
Winterhalter P, Schreier P: C13-Norisoprenoid glycosides in plant tissues: an overview on their occurrence, composition and role as flavour precursors. Flav Frag J. 1994, 9: 281-287. 10.1002/ffj.2730090602.
Giuliano G, Al-Babili S, von Lintig J: Carotenoid oxygenases: cleave it or leave it. Trends Plant Sci. 2003, 8: 145-149. 10.1016/S1360-1385(03)00053-0.
Tan BC, Joseph LM, Deng WT, Liu L, Li QB, Cline K, McCarty DR: Molecular characterization of the Arabidopsis 9-cis epoxycarotenoid dioxygenase gene family. Plant J. 2003, 35: 44-56. 10.1046/j.1365-313X.2003.01786.x.
Auldridge ME, Block A, Vogel JT, Dabney-Smith C, Mila I, Bouzayen M, Magallanes-Lundback M, DellaPenna D, McCarty DR, Klee HJ: Characterization of three members of the Arabidopsis carotenoid cleavage dioxygenase family demonstrates the divergent roles of this multifunctional enzyme family. Plant J. 2006, 45: 982-993. 10.1111/j.1365-313X.2006.02666.x.
Auldridge ME, McCarty DR, Klee HJ: Plant carotenoid cleavage oxygenases and their apocarotenoid products. Curr Op Plant Biol. 2006, 9: 1-7. 10.1016/j.pbi.2006.03.005.
Ibdah M, Azulay Y, Portnoy V, Wasserman B, Bar E, Meir A, Burger Y, Hirschberg J, Schaffer AA, Katzir N, Tadmor Y, Lewinsohn E: Functional characterization of CmCCD1, a carotenoid cleavage dioxygenase from melon. Phytochem. 2006, 67: 1579-1589. 10.1016/j.phytochem.2006.02.009.
Vogel JT, Tan BC, McCarty DR, Klee HJ: The carotenoid cleavage dioxygenase 1 enzyme has broad substrate specificity, cleaving multiple carotenoids at two different bond positions. J Biol Chem. 2008, 283: 11364-11373. 10.1074/jbc.M710106200.
Simkin AJ, Schwartz SH, Auldridge M, Taylor MG, Klee HJ: The tomato carotenoid cleavage dioxygenase 1 genes contribute to the formation of the flavor volatiles β-ionone, pseudoionone, and geranylacetone. Plant J. 2004, 40: 882-892. 10.1111/j.1365-313X.2004.02263.x.
Simkin AJ, Underwood BA, Auldridge M, Loucas HM, Shibuya K, Schmelz E, Clark DG, Klee HJ: Circadian regulation of the PhCCD1 carotenoid cleavage dioxygenase controls emission of beta-ionone, a fragrance volatile of petunia flowers. Plant Physiol. 2004, 136: 3504-3514. 10.1104/pp.104.049718.
Campbell R, Ducreux LJM, Morris WL, Morris JA, Suttle JC, Ramsay G, Bryan GJ, Hedley PE, Taylor MA: The metabolic and developmental roles of carotenoid cleavage dioxygenase4 from potato. Plant Physiol. 2010, 154: 656-664. 10.1104/pp.110.158733.
Ohmiya A, Kishimoto S, Aida R, Yoshioka S, Sumitomo K: Carotenoid cleavage dioxygenase (CmCCD4a) contributes to white color formation in chrysanthemum petals. Plant Physiol. 2006, 142: 1193-1201. 10.1104/pp.106.087130.
Lewinsohn E, Sitrit Y, Bar E, Azulay Y, Meir A, Zamir D, Tadmor Y: Carotenoid pigmentation affects the volatile composition of tomato and watermelon fruits, as revealed by comparative genetic analyses. J Agric Food Chem. 2005, 20: 3142-3148. 10.1021/jf047927t.
Baldermann S, Naim M, Fleischmann P: Enzymatic carotenoid degradation and aroma formation in nectarines (Prunus persica). Food Res Int. 2005, 38: 833-836. 10.1016/j.foodres.2005.02.009.
Zhang M, Leng P, Zhang G, Li X: Cloning and functional analysis of 9-cis-epoxycarotenoid dioxygenase (NCED) genes encoding a key enzyme during abscisic acid biosynthesis from peach and grape fruits. J Plant Physiol. 2009, 166: 1241-1252. 10.1016/j.jplph.2009.01.013.
Han SY, Kitahata N, Sekimata K, Saito T, Kobayashi M, Nakashima K, Yamaguchi-Shinozaki K, Shinozaki K, Yoshida S, Asami T: A novel inhibitor of 9-cis-epoxycarotenoid dioxygenase in abscisic acid biosynthesis in higher plants. Plant Physiol. 2004, 135: 1574-1582. 10.1104/pp.104.039511.
Kitahata N, Han SY, Noji N, Saito T, Kobayashi M, Nakano T, Kuchitsu K, Shinozaki K, Yoshida S, Matsumoto S, Tsujimoto M, Asami T: A 9-cis-epoxycarotenoid dioxygenase inhibitor for use in the elucidation of abscisic acid action mechanisms. Bioorg Med Chem. 2006, 14: 5555-5561. 10.1016/j.bmc.2006.04.025.
Conte L, Della Strada G, Fideghelli C, Insero O, Liverani A, Moser L, Nicotra A: Redhaven Bianca. Monografia di cultivar di pesche, nettarine e percoche. Edited by: Conte L, Della Strada G, Fideghelli C, Insero O, Liverani A, Moser L, Nicotra A. 1994, Roma: Mi.R.A.A.F.-I.S.F, 94-
Dermen H, Stewart RN: Ontogenetic study of floral organs of peach (Prunus persica) utilizing cytochimeral plants. Am J Bot. 1972, 60: 283-291. 10.2307/2441220.
Farré G, Sanahuja G, Naqvi S, Bai C, Capell T, Zhu C, Christou P: Travel advice on the road to carotenoids in plants. Plant Sci. 2010, 179: 28-48.
Giuliano G, Bartley GE, Scolnik PA: Regulation of carotenoid biosynthesis during tomato development. Plant Cell. 1993, 5: 379-387. 10.1105/tpc.5.4.379.
Diretto G, Al-Babili S, Tavazza R, Scossa F, Papacchioli V, Migliore M, Beyer P, Giuliano G: Transcriptional-metabolic networks in β-carotene-enriched potato tubers: the long and winding road to the Golden phenotype. Plant Physiol. 2010, 154: 899-912. 10.1104/pp.110.159368.
Rodríguez-Villalón A, Gas E, Rodríguez-Concepción M: Phytoene synthase activity controls the biosynthesis of carotenoids and the supply of their metabolic precursors in dark-grown Arabidopsis seedlings. Plant J. 2009, 60: 424-435.
Marty I, Bureau S, Sarkissian G, Gouble B, Audergon JM, Albagnac G: Ethylene regulation of carotenoid accumulation and carotenogenic gene expression in colour-contrasted apricot varieties (Prunus armeniaca). J Exp Bot. 2005, 56: 1877-1886. 10.1093/jxb/eri177.
Diretto G, Welsch R, Tavazza R, Mourgues F, Pizzichini D, Beyer P, Giuliano G: Silencing of beta-carotene hydroxylase increases total carotenoid and beta-carotene levels in potato tubers. BMC Plant Biol. 2007, 7: 11-10.1186/1471-2229-7-11.
Yan J, Kandianis CB, Harjes CE, Bai L, Kim EH, Yang X, Skinner DJ, Fu Z, Mitchell S, Li Q, Fernandez MG, Zaharieva M, Babu R, Fu Y, Palacios N, Li J, DellaPenna D, Brutnell T, Buckler ES, Warburton ML, Rocheford T: Rare genetic variation at Zea mays crtRB1 increases β-carotene in maize grain. Nat Genet. 2010, 42: 322-327. 10.1038/ng.551.
Eduardo I, Chietera G, Bassi D, Rossini L, Vecchietti A: Identification of key odor volatile compounds in the essential oil of nine peach accessions. J Sci Food Agric. 2009, 90: 1146-1154. 10.1002/jsfa.3932.
Wang Y, Yang C, Li S, Yang L, Wang S, Zhao J, Jiang Q: Volatile characteristics of 50 peaches and nectarines evacuate by HP-SPME with GC-MS. Food Chem. 2009, 116: 356-364. 10.1016/j.foodchem.2009.02.004.
Bouvier F, Suire C, Mutterer J, Camara B: Oxidative remodeling of chromoplast carotenoids: identification of the carotenoid dioxygenase CsCCD and CsZCD genes involved in Crocus secondary metabolite biogenesis. Plant Cell. 2003, 15: 47-62. 10.1105/tpc.006536.
Mathieu S, Terrier N, Procureur J, Bigey F, Gunata Z: A carotenoid cleavage dioxygenase from Vitis vinifera L.: functional characterization and expression during grape berry development in relation to C13-norisoprenoid accumulation. J Exp Bot. 2005, 56: 2721-2731. 10.1093/jxb/eri265.
Huang FC, Molnár P, Schwab W: Cloning and functional characterization of carotenoid cleavage dioxygenase 4 genes. J Exp Bot. 2009, 60: 3011-3022. 10.1093/jxb/erp137.
Rubio A, Rambla JL, Santaella M, Gomez MD, Orzaez D, Granell A, Gómez-Gómez L: Cytosolic and plastoglobule-targeted carotenoid dioxygenases from Crocus sativus are both involved in beta-ionone release. J Biol Chem. 2008, 283: 24816-24825. 10.1074/jbc.M804000200.
Ahrazem O, Trapero A, Gómez MD, Rubio-Moraga A, Gómez-Gómez L: Genomic analysis and gene structure of the plant carotenoid dioxygenase 4 family: A deeper study in Crocus sativus and its allies. Genomics. 2010, 96: 239-250. 10.1016/j.ygeno.2010.07.003.
Yamamizo C, Kishimoto S, Ohmiya A: Carotenoid composition and carotenogenic gene expression during Ipomoea petal development. J Exp Bot. 2010, 61: 709-719. 10.1093/jxb/erp335.
Carbone F, Mourgues F, Biasioli F, Gasperi F, Märk TD, Rosati C, Perrotta G: Development of molecular and biochemical tools to investigate fruit quality traits in strawberry elite genotypes. Mol Breed. 2006, 18: 127-142. 10.1007/s11032-006-9017-2.
Zhang B, Shen JY, Wei WW, Xi WP, Xu CJ, Ferguson I, Chen K: Expression of genes associated with aroma formation derived from the fatty acid pathway during peach fruit ripening. J Agric Food Chem. 2010, 58: 6157-6165. 10.1021/jf100172e.
Lazzari B, Caprera A, Vecchietti A, Merelli I, Barale F, Milanesi L, Stella A, Pozzi C: Version VI of the ESTree db: an improved tool for peach transcriptome analysis. BMC Bioinf. 2008, 9 (Suppl 2): S9-10.1186/1471-2105-9-S2-S9.
Boudehri K, Bendahmane A, Cardinet G, Troadec C, Moing A, Dirlewanger E: Phenotypic and fine genetic characterization of the D locus controlling fruit acidity in peach. BMC Plant Biol. 2009, 9: 59-10.1186/1471-2229-9-59.
Vizoso P, Meisel LA, Tittarelli A, Latorre M, Saba J, Caroca R, Maldonado J, Cambiazo V, Campos-Vargas R, Gonzalez M, Orellana A, Silva H: Comparative EST transcript profiling of peach fruits under different post-harvest conditions reveals candidate genes associated with peach fruit quality. BMC Genomics. 2009, 10: 423-10.1186/1471-2164-10-423.
Bustin SA, Beaulieu JF, Huggett J, Jaggi R, Kibenge FSB, Olsvik PA, Penning LC, Toegel S: MIQE précis: Practical implementation of minimum standard guidelines for fluorescence-based quantitative real-time PCR experiments. BMC Mol Biol. 2010, 11: 74-10.1186/1471-2199-11-74.
López-Gómez R, Gómez-Lim MA: A method for extracting intact RNA from fruit rich in polysaccharides using ripe mango mesocarp. HortSci. 1992, 27: 440-442.
Cecchi F: Studio dell'espressione genica degli enzimi preposti alla biosintesi dei carotenoidi nel corso della maturazione di frutti di pesco (Prunus persica). BS Thesis, Viterbo University (Italy), Food Science and Technology Department;2004.
International Peach Genome Initiative. http://www.peachgenome.org.
Davidovich-Rikanati R, Lewinsohn E, Bar E, Iijima Y, Pichersky E, Sitrit Y: Overexpression of the lemon basil α-zingiberene synthase gene increases both mono- and sesquiterpene contents in tomato fruit. Plant J. 2008, 56: 228-238. 10.1111/j.1365-313X.2008.03599.x.
Molnár P, Szabolcs J: Alkaline permanganate oxidation of carotenoid epoxides and furanoids. Acta Chim Acad Sci Hung. 1979, 99: 155-173.
Schiedt K, Liaaen-Jensen S: Isolation and analysis. Carotenoids Vol. 1A: Isolation and analysis. Edited by: Britton G, Liaaen-Jensen S, Pfander H. Basel-Boston-Berlin: Birkhäuser Verlag; 1995, 109-144.
Molnár P, Kawase M, Motohashi N: Isolation, crystallization and handling of carotenoids. Functional polyphenols and carotenoids with antioxidative action. Edited by: Motohashi N. Kerala: Chem. Pharm. Sci. (CPS); 2005, 111-131.
Pichersky E, Noel JP, Dudareva N: Biosynthesis of plant volatiles: Nature's diversity and ingenuity. Science. 2006, 311: 808-811. 10.1126/science.1118510.
Financial support by CARFLAVO and FIRB-Parallelomics projects to FB, FM, GG, AL and CR is acknowledged. The Hungarian authors were supported by the grants OTKA K 76 176 (Hungarian National Research Foundation) and PTE AOK KA-34039-35/2009 (Research Found of the Faculty of Medicine, University of Pécs). Authors wish to thank dr. Tadao Asami (University of Tokyo) for the kind gift of abamine SG, and profs. Péter Molnár and József Deli (University of Pécs) for their collaborative support in isolation of carotenoids and in HPLC investigations. The International Peach Genome Initiative (IPGI) is acknowledged for early online access to the draft genome sequence, which enabled the analysis of ccd4 gene expression.
FB did the samplings, organized the work, carried out molecular experiments and statistical analyses, and contributed to manuscript writing; EB carried out VOC extractions and GC-MS analyses; FM designed and participated to molecular experiments and data analysis, and contributed to manuscript writing; GH and ET performed carotenoid extractions, HPLC and data analyses; GG contributed early research design and molecular data, and to manuscript writing; AL and ST tutored FB, helped in field experiments and data analysis, and contributed to manuscript writing; EL designed GC-MS analyses and interpreted their data, and contributed to manuscript writing; CR designed and coordinated the research, contributed to statistical analyses and wrote the manuscript. All authors read and approved the final manuscript.
Einat Bar, Fabienne Mourgues, Györgyi Horváth, Erika Turcsi contributed equally to this work.
Electronic supplementary material
Additional File 1: Structures of the main carotenoids identified in RHB and RH fruits during ripening. Carotenoid composition is reported in Table 1. (PPT 138 KB)
Additional File 2: Major quality traits of RHB and RH fruits measured at S4 stage. SSC: soluble solids content (expressed in °Brix). TA: titratable acidity (expressed in meq NaOH/l). Firmness was measured with an 8-mm diameter tip (expressed in kg/cm2). Skin color parameters a* (red chromatic coordinate), b* (yellow chromatic coordinate) and L* (brightness) were recorded at two fruit cheeks (opposite equatorial points). Values are average measurements of ten representative fruits. Different letters indicate significant differences among mean values (t test; p ≤ 0.05). (DOC 29 KB)
Additional File 3: Hierarchical clustering analysis of carotenoid gene expression in RH and RHB genotypes. A: joint analysis of RHB and RH data. B: RHB data only. C: RH data only. Each cell corresponds to the relative expression value (Log-transformed) according to the color scale on the right. For enzyme abbreviations and fruit development stages, see text and Methods, respectively. (PPT 213 KB)
Additional File 4: Total VOC content in RHB and RH mesocarp during fruit ripening. RH: solid black squares. RHB: open squares. Values ± SD are in ng/g fresh weight. (PPT 120 KB)
Additional File 5: Accumulation patterns of identified norisoprenoids in RHB and RH mesocarp during fruit ripening. RH: solid black symbols. RHB: open symbols. Values are in ng/g fresh weight. (PPT 125 KB)
Additional File 6: Sequences of RT-qPCR primers used in this work. for experimental conditions, see Methods. (DOC 44 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Brandi, F., Bar, E., Mourgues, F. et al. Study of 'Redhaven' peach and its white-fleshed mutant suggests a key role of CCD4 carotenoid dioxygenase in carotenoid and norisoprenoid volatile metabolism. BMC Plant Biol 11, 24 (2011). https://doi.org/10.1186/1471-2229-11-24
- Flesh Color
- Carotenoid Accumulation
- Peach Fruit