Exploring the differential mechanisms of carotenoid biosynthesis in the yellow peel and red flesh of papaya
BMC Genomics volume 20, Article number: 49 (2019)
Red-fleshed papaya is a good material to study the different carotenoids accumulation mechanism in the peel and flesh. Although the peel and flesh of papaya closely integrated into one body, the flesh coloration changing from white to red, while the exocarp coloration changing from green to yellow. In this study, the major carotenoids accumulation and the expression patterns of key carotenoid biosynthesis pathway genes in the process of papaya fruit ripening were studied, and the carotenoid biosynthetic pathways in the yellow peel and red flesh of papaya were investigated.
The carotenoid composition in papaya flesh and peel were different. The major carotenoids were lutein and β-carotene in the peel, while lycopene in the flesh. The accumulation of carotenoids, including lycopene, β-carotene, and β-cryptoxanthin were considered to cause the orange-red color of papaya cv. ‘Daqing No.10’ flesh. The color of peel changed from green to yellow because of the fast degradation of chlorophyll and the appearance of carotenoids such as lutein and β-carotene. Thirteen genes that encode enzymes in the carotenoid biosynthetic pathway were detected in papaya fruit transcriptome: two phytoene synthase (PSY1, PSY2), two phytoene desaturase (PDS1, PDS2), one ζ-carotene desaturase (ZDS), four lycopene cyclase (CYCB, LCYB1, LCYB2, LCYE), one β-carotene hydroxylase (CHYB), one carotene ε-monooxygenase (LUT1), one violaxanthin de-epoxidase (VDE), and one zeaxanthin epoxidase (ZEP). The results of RNA-Seq and RT-qPCR showed the expression of carotenoid biosynthetic pathway genes was consistent with the change of carotenoid content. Carotenoid biosynthetic pathways in the yellow peel and red flesh of papaya were analysed based on the major carotenoids accumulation and the expression patterns of key carotenoid biosynthesis pathway genes. There was only a β-branch of carotenoid biosynthesis in the flesh of papaya, while there were both α- and β-branch of carotenoid biosynthesis in papaya peel. In the process of papaya fruit ripening, the α-branch was inhibited and the β-branch was enhanced in the peel.
The differential carotenoid accumulation and biosynthesis pathway genes expression in peel and flesh, lay a foundation for further study and provide further insights to control fruit color and improve fruit quality and appearance.
Carotenoids are one of the most essential components for human nutrition and health, mainly due to their pro-vitamin A and antioxidant proportion. Humans are not able to synthesize carotenoids and depend entirely on natural sources or dietary supplements. Dietary carotenoids are thought to provide health benefits in decreasing the risk of disease, particularly certain cancers and eye disease [1, 2]. The most common dietary carotenoids are α-carotene, β-carotene, β-cryptoxanthin, lycopene, lutein, and zeaxanthin; α-carotene, β-carotene, and β-cryptoxanthin are precursors of vitamin A, while lycopene, lutein, zeaxanthin have no vitamin A activity but can also act as antioxidants. Besides the nutritional value of carotenoids, they are the main coloring pigments of the yellow, orange, and red colors of fruits: mango (Mangifera indica L.) , tomato (Lycopersicon esculentum Mill.) , orange (Citrus) , banana (Musa spp.) , papaya (Carica papaya L.) , and loquat (Eriobotrya japonica L.) .
The biosynthetic pathway of carotenoids in high plants has been gradually clarified by biochemical, classical genetic and molecular genetic analysis. Plant carotenoids are isoprenoid-derived molecules generally synthesized and located in plastids. Two Geranylgeranyl pyrophosphate (GGPP) molecules condense via phytoene synthase (PSY), forming phytoene; phytoene forms lycopene via phytoene desaturase (PDS) and ζ-carotene desaturase (ZDS); then the biosynthetic pathway of carotenoids splits after biosynthesis of lycopene into the α-branch and the β-branch: ① in α-branch, lycopene cyclase (LCY) catalyzes the cyclization of lycopene to form α-carotene; α-carotene is converted into lutein by β-carotene hydroxylase (CHYB) and carotene ε-monooxygenase (LUT1); ② in β-branch, LCY catalyzes the cyclization of lycopene to form β-carotene, and CHYB catalyzes the conversion of β-carotene to β-cryptoxanthin and β-cryptoxanthin to zeaxanthin [9, 10]. Many carotenoid synthesis genes including PSY, PDS, ZDS, and LCY have been cloned successfully and genetically modified plants with increased carotenoids levels have been obtained such as rice, tomato, potato, and canola [11, 12].
Papaya is an important tropical and sub-tropical fruit crop which is known for its high nutritional values like vitamins A and vitamins C [13, 14]. There are two types of papaya, red-fleshed and yellow-fleshed. The major carotenoid in the pulp of red-fleshed papaya is lycopene, while the major carotenoids in the yellow-fleshed papaya are β-carotene and β-cryptoxanthin [15, 16]. It has been documented that a 2-bp insertion is present in the chromoplast specific lycopene β-cyclase (CYCB) gene resulting in a recessive loss-of-function mutation, so as to accumulate a large amount of lycopene which is responsible for the red flesh . Pigments in the flesh of papaya are known, but pigments in the peel of papaya remain unknown. The color of exocarp turns from green to yellow and the color of the flesh changes from white to red during the fruit ripening process in red-fleshed papaya. The differential accumulation of carotenoids and differential expression of carotenoid synthesis genes were studied in order to understand the different carotenoids accumulation mechanism in the peel and flesh of papaya.
Papaya coloration during papaya fruit ripening
C. papaya cv. ‘Da Qing No.10’ is a red-fleshed papaya. During fruit ripening process, the most visible change of papaya is the color of exocarp turning from green to yellow and the color of flesh changing from white to orange-red (Fig. 1). As seen in Fig. 1a, at green stage, the peel is dark green and the flesh is white; at color break stage, most of the peel is still green, but there are two yellow stripes appear on the top of fruit, as well as the flesh is red; at half yellow stage, the color of peel is half green and half yellow; at full yellow stage, the color of peel is full yellow; the flesh and peel closely integrated into one body under the anatomical lens. The peel is not as circumscribed from the flesh. Under a microscope, the flesh cells are bigger than exocarp cells in FY red-fleshed papaya. Yellow round-shaped elements are abundant and well-distributed in peel cells with larger yellow round-shaped and several red elongated contours crystals present in flesh cells (Fig. 1b, c). Thus, the color and shape of chromoplasts have a clear difference between flesh and peel. These findings illustrated the differences in carotenoid composition between peel and flesh.
Chlorophyll change and carotenoid accumulation
The content of chlorophyll and major carotenoids was determined during the ripening of papaya fruit. The content of chlorophyll in the peel decreased during papaya ripening, and there was almost no chlorophyll detected in the flesh. The color of green stage in papaya was green because of the high contents of chlorophyll, after that the content of chlorophyll decreased sharply during the process of papaya ripening (Fig. 2). A large amount of lutein, β-carotene, as well as small amount of β-cryptoxanthin, zeaxanthin, lycopene were detected in the peel. The content of β-cryptoxanthin, zeaxanthin, and lycopene increased, while and the content of lutein and β-carotene decreased in the peel during papaya ripening. The presence of lutein, β-carotene, β-cryptoxanthin, zeaxanthin and the simultaneous rapid loss of chlorophyll caused a color shift from green to yellow during fruit ripening.
Almost no chlorophyll, lutein, and zeaxanthin were detected in the flesh, but a large amount of lycopene, some β-carotene, and a small amount of β-cryptoxanthin were detected in the flesh. In the process of fruit ripening, there was a sustained increase in carotenoids (lycopene, β-carotene, β-cryptoxanthin) content. Therefore, the flesh color changed from white to red. Except for lycopene, the other carotenoids (lutein, β-carotene, β-cryptoxanthin, zeaxanthin) content in peels were higher than that in fleshes. This illustrates that exocarp is an important source of carotenoids.
The results of RNA-Seq are shown in Table 1 and Fig. 3. Table 1 showed the clean reads number and the gene number of the four samples. After removing low-quality reads, adaptor sequences, and rRNA reads, we obtained 44,431,496 (PE1), and 41,487,992 (PE2), 44,812,308 (FL1), 52,110,632 (FL2) clean reads. Clean reads were then mapped to the papaya reference genome. There were 18,941 (PE1), 18,242 (PE2), 18,315 (FL1), and 17,907 (FL2) expressed genes detected, with the genes in the FL2 being the least. More expressed genes at the green stage were than those at color break stage. Figure 3 showed the Differentially Expressed Genes (DEGs) in different samples. Compared with FL1, 2663 genes were up-regulated and 4501 genes were down-regulated in FL2; in papaya peel, there were 2317 genes up-regulated and 3749 genes down-regulated. This demonstrates that a large number of genes may participate in the process of fruit ripening. Comparing with the peel, the number of down-regulated genes is more than the up-regulated genes’ number in the flesh.
The pathway enrichment analysis was carried out using the KEGG pathway database. We focused on the analysis of carotenoid biosynthetic pathway. Comparing with FL1, the expression of PDS, ZDS, CYCB (chromoplast specific lycopene β-cyclase), CHYB, and VDE (Violaxanthin de-epoxidase) increased, while the expression of LCYB2 decreased in FL2 (Additional file 1: Figure S1). Also, the expression of most carotenoid biosynthetic pathway genes (PDS, ZDS, CYCB, LUT1, CHYB, VDE, and ZEP (zeaxanthin epoxidase)) in the peel increased with papaya ripening, while the expression of LCYE decreased (Additional file 2: Figure S2). Comparing of carotenoid biosynthetic genes in the peel and flesh of papaya, we found that LCYE was the only gene that had been enriched in FL2 by KEGG pathway enrichment. The expression level of LCYE in FL2 was significantly lower than that in PE2 (Additional file 3: Figure S3). This indicated that the differential expression of LCYE may be the key reason for the difference of carotenoid components between flesh and peel in color break stage. As previously described, the major carotenoid in the peel was lutein, but it was not detected in the flesh. There was no lutein in the flesh because of non-detectable expression of LCYE.
Screening carotenoid biosynthesis related genes in the transcriptome
We analyzed the genes involved in carotenoid biosynthesis. Thirteen papaya genes that encode enzymes in the carotenoid biosynthetic pathway were detected in papaya transcriptome: two phytoene synthase (PSY1, PSY2), two phytoene desaturase (PDS1, PDS2), one ζ-carotene desaturase (ZDS), four lycopene cyclase (CYCB, LCYB1, LCYB2, LCYE), one β-carotene hydroxylase (CHYB), one carotene ε-monooxygenase (LUT1), one violaxanthin de-epoxidase (VDE), and one zeaxanthin epoxidase (ZEP) (Additional file 4: Table S1). The heat map of the carotenoid biosynthesis related genes was drawn according to FPKM (Fragments Per Kilobase of transcript per Million mapped reads) values (Fig. 4). LCYB2 and LCYE were closely related in a cluster, while the other ten genes are in a big cluster (Fig. 4). Except for LCYB2 and LCYE, the expression of most carotenoid biosynthetic pathway genes (PDS, ZDS, CYCB, LUT1, CHYB, VDE, and ZEP) increased during papaya ripening.
Phytoene synthase (PSY) is a key regulator in the carotenoid biosynthetic pathway ; Two PSY genes were detected in the transcriptome. The expression of PSY1 increased in the flesh and peel, and no expression of PSY2 was detected, thereby demonstrating that PSY1 was active in the carotenoid biosynthetic pathway in papaya. Similar to papaya, PSY1 plays an important role in tomato fruit tissue pigmentation too . Two PDS genes were detected, and the expression patterns were the same, but the expression abundance between the two genes was different. The expression of PDS1 in FL2 is 394.46, while that of PDS2 is 7.9 (Additional file 4: Table S1). Because the abundant gene transcripts account for more enzymatic activity participating in carotenoid synthesis, PSY1 and PDS1 are the major genes involved in the synthesis of carotenoid in fruit. The expression of two genes (ZDS and CHYB) increased with papaya ripening. Four lycopene cyclase genes were detected: LCYB1, LCYB2, CYCB, and LCYE. LCYB1 and CYCB had been cloned and characterized previously [20, 21]. We didn’t analyze CYCB because there is a 2-bp insertion in CYCB gene in red-fleshed papaya, CYCB is a functional gene in yellow-fleshed papaya, but a pseudo gene in red-fleshed papaya . The expression of the other three lycopene cyclase genes was different. This showed that they may play different roles in carotenoid synthesis.
RT-qPCR analysis of transcription levels of five carotenoid biosynthetic pathway genes
The expression of five key carotenoid biosynthetic pathway genes PSY1, LCYB1, LCYB2, LCYE, CHYB, in different ripening papaya fruits and different organs was studied using RT-qPCR analysis. The transcript levels of PSY1 increased first and then decreased with papaya ripening (Fig. 5b). The expression pattern of PSY1 was similar in the flesh and peel. LCYB1 also expressed in both flesh and peel, and the transcript levels increased first then decreased too. The expression of LCYB2 and LCYE in the peel were higher than that in the flesh, and the transcript levels were decreased with papaya ripening. The transcript level of CHYB was higher in the peel than that in the flesh, and the CHYB’s expression presented a trend of rising. Figure 5b also shows the results of the relative transcript levels of the five carotenoid biosynthesis-related genes in different organs (root, leaf, stem, and flower). PSY1 expressed in leaf, stem, and ripening fruit; LCYB1 expressed in leaf, stem, fruit, flower, while almost no expression was detected in the root; LCYB2 has a higher expression in leaf, stem, and flower, a lower expression in fruit. LCYE transcripts predominate in green color tissues, such as green peel, leaf, and stem.
Carotenoid biosynthetic pathway in the yellow peel and red flesh of papaya
Based on the major carotenoids accumulation and the expression patterns of key carotenoid biosynthesis pathway genes, the carotenoid biosynthetic pathway in the yellow peel and red flesh of papaya were speculated (Fig. 6). There were both α- and β-branch of carotenoid biosynthesis in the peel of papaya, while there was only β-branch in the flesh. A large amount of lutein, β-carotene, as well as small amount of β-cryptoxanthin, zeaxanthin, lycopene were in the peel. The most abundant carotenoid in the flesh was lycopene, followd by β-carotene and β-cryptoxanthin. In the peel, the key carotenoid biosynthesis pathway genes: PSY1, PDS1, ZDS, LCYB1, CHYB, LUT1, ZEP, and VDE were up-regulated during papaya fruit ripening, while LCYE and LCYB2 were down-regulated; in the flesh, PSY1, PDS1, ZDS, LCYB1, CHYB, ZEP, and VDE were also up-regulated during papaya fruit ripening.
Although the flesh and peel of red-fleshed papaya fruit are closely integrated into one body, the coloration of the flesh and peel are quite different: the color of exocarp turning from green to yellow and the color of flesh changing from white to red. So exploring the differential mechanisms of carotenoid biosynthesis in the yellow peel and red flesh of papaya is an interesting research point. The color transformation is directly related to the change of the carotenoid component. Some studies showed that the concentration of carotenoids was significantly higher in the peel than that in the flesh such as Citrus iyo fruit  and loquat fruit . The coloration of loquat peel due to the transient increase of carotenoids and degradation of chlorophylls . In papaya fruit, carotenoids (lycopene, β-carotene, β-cryptoxanthin) accumulates extensively in the flesh, therefore the flesh color changes from white to red; although the content of lutein and β-carotene decreased in the peel, the rapid loss of chlorophyll leads to the appearance of carotenoids causes a color shift from green to yellow. Microscopic observation also proved that the carotenoid composition between the peel and the flesh are different. Yellow round-shaped elements are abundant in peel cells, while larger yellow round-shaped and several red elongated contours crystals present in flesh cells. This result is consistent with previous research that β-carotene is deposited in globular and tubular elements in a lipid-dissolved form whereas lycopene is deposited in a solid crystalline form [23, 24].
The synthesis of carotenoid in peel and flesh of papaya had been analysed base on the distribution of carotenoid components and the expression of carotenoid biosynthesis pathway genes. The increase in the content of β-cryptoxanthin and zeaxanthin was accompanied by the decrease in the content of lutein and β-carotene, indicating that there were both α- and β-branch in the peel and the β-branch was enhanced in the process of papaya fruit ripening. However, there is only a β-branch in the flesh because no lutein and expression of LCYE were detected, LCYE was the only gene that had been enriched in FL2 by KEGG pathway enrichment. It expressed abundantly in the green peel, but hardly expressed in the flesh. LCYE is the key reason for the differential carotenoid accumulation in the flesh and peel of papaya. Except for LCYE, three lycopene cyclases had been detected in papaya fruit transcriptome, LCYB1, CYCB, and LCYB2. LCYB1 has been considered as a chloroplast-specific lycopene β-cyclase as well as CYCB is a chromoplast-specific lycopene β-cyclase [20, 21]. Although the gene of CYCB was named as LCY-β2 in previous studies, CYCB was used in our paper [20, 25]. LCYB1 was supposed to be a chloroplast-specific lycopene β-cyclase, but it had higher expression level in red flesh and the downstream products of lycopene β-cyclase like β-carotene and β-cryptoxanthin were also detected, which means LCYB1 expresses in non-chlorophyllous tissue as well. This is not consistent with previous studies . LCYB2 is a newly discovered gene of lycopene cyclase in papaya, the expression of LCYB2 in peel was higher than that in flesh, and the transcript levels had been decreasing with papaya fruit ripening. However, the function of LCYB2 is still unclear in papaya.
The reason for the conversion from α-branch to β-branch in the peel of papaya during fruit ripening is not clear. One of the key roles of carotenoids in plants is protecting the photosynthetic apparatus against photoinhibition and photodamage . The xanthophyll cycle transforms the excitation energy into heat and thereby prevents the formation of damaging ROS. It is an important photoprotective process in plants . Violaxanthin de-epoxidase (VDE) and zeaxanthin epoxidase (ZEP) are the two enzymes in the xanthophyll cycle. In excess light conditions, VDE catalyzes the conversion of violaxanthin to zeaxanthin via antheraxanthin, whereas ZEP catalyzes the reverse reaction . The expression of ZEP and VDE were up-regulated in color break stage of papaya, this may due to chlorophyll degradation enhanced the activity of the xanthophyll cycle to avoid severe photodamage under strong illumination.
By above knowable, the differential expression of the key carotenoid biosynthesis pathway genes in the peel and flesh of papaya result in the differential accumulation of carotenoids. But little is known about the regulatory mechanisms of the genes involved in ɑ- and β-branch. Transcriptional factors (TFs) play important roles in a variety of cellular and developmental processes, such as fruit ripening, senescence, hormone signalling, and various biotic and abiotic stress responses. For example, several TFs have been reported to be involved in the ripening of banana , tomato , and papaya . Recently, some TFs were found to play essential roles in regulating the carotenoid biosynthesis in plants. An R2R3-MYB represses the transformation of α- and β-branch carotenoids by negatively regulating expression of CrBCH2 and CrNCED5 in flavedo of Citrus reticulate . CpNAC1 is a positive regulator of carotenoid biosynthesis during papaya fruit ripening . The protein of CpNAC2 physically interacts with CpEIN3a, which act as a transcriptional activator of CpPDS4 and CpCHY-b by directly binding to their promoters has been proposed . Several TFs including MYB, ERF, bZIP, bHLH, WRKY, and NAC had been selected according to the expression patterns with papaya fruit ripening in our research (data not given). Some TFs expression decreased in both flesh and peel, while some TFs were gradually increased, paralleling the increase of the carotenoids content; some TFs have higher expression levels in the flesh, and some TFs have much higher expression levels in the peel. The differential expression of these TFs suggests that they may be involved in the regulation of carotenoid biosynthesis genes transcription in flesh or peel. The functions of these TFs should be studied in the future.
Carotenoid biosynthetic pathway in the yellow peel and red flesh of papaya was revealed on the basis of the major carotenoids accumulation and the expression patterns of key carotenoid biosynthesis pathway genes. There is only a β-branch of carotenoid biosynthesis in the flesh of papaya, while there are both α- and β-branch in papaya peel. However, additional genomics wide analysis are required to determine and validate. This work elucidated the differential anabolic pathway of carotenoids in peel and flesh of papaya, and provide further insights to control fruit color and improve fruit quality and appearance.
Papaya fruits (C. papaya cv. ‘Daqing No.10’) in different development stages (GS, Green stage; CB, Color break stage; HY, Half yellow stage; FY, Full yellow stage) were harvested from Longzhijia fruit and vegetable farmers’ specialized cooperatives in November in Zhangzhou, China. The healthy fruits of similar size, shape, and maturity were divided into four groups. Five papaya in the same stage were peeled, seeds were removed, and the peel and flesh were cut into some pieces respectively. The pieces were mixed and weighed, frozen in liquid nitrogen, and stored at − 80 °C. Finally, eight types of samples were obtained, namely FL1 (the flesh of GS), FL2 (the flesh of CB), FL3 (the flesh of HY), FL4 (the flesh of FY), PE1 (the peel of GS), PE2 (the peel of CB), PE3 (the peel of HY), PE4 (the peel of FY). Tissue samples (root, stem, leaf, and flower) were collected for reverse transcription quantitative PCR (RT-qPCR) analysis. For chromoplast observation, small pieces of flesh and peel of FY stage papaya were obtained respectively. All pieces were cut with a sharp razor blade to thin sections. Then the ultrathin sections were put on a clean slide and examined with an Olympus IX51 (Tokyo, Japan).
Spectrophotometric determination of chlorophyll content
Chlorophyll was extracted from 200 mg of ground papaya fruit materials using 25 mL of 80% acetone. The absorbance was determined at 663 nm and 646 nm with a UV-Vis spectrophotometer equipped with 1.0 cm quartz cells. Based on the equation of total chlorophyll concentration CT = 7.18 × A663 + 17.32 × A646, total chlorophyll content was determined .
UPLC analysis of major carotenoids
An modified method for the extraction of major carotenoids from papaya fruit was used in the study . Briefly, the sample (peel or flesh) was lyophilized using a freeze dryer at − 40 °C for two days, followed by mechanical grinding with liquid nitrogen, and stored at − 20 °C. papaya fruit (250 mg) was added into 2.5 mL of acetone (1:10). Ultrasonic-assisted extraction was used to extract carotenoids for 1 h at 240 W at 20 °C, and then centrifuged for 10 min at 10,000 rpm at 4 °C. The liquid supernatant was filtered by 0.22 μm organic filter, and stored at − 20 °C for later use.
Carotenoid pigments analysis was referenced to Yan  and Lee . Chromatography was carried out with a Waters Ultra-high performance liquid chromatography. Samples were measured by a C18 carotenoid column (100 mm × 2.1 i.d., 1.7 μm) from Waters. The eluent was methyl tertiary butyl ether and methanol (V/V = 30:70) with the type of isocratic elution. And each eluent contained 0.01% BHT and 0.05% TEA (triethylamine) as modifiers in order to prevent the degradation of carotenoids on the column. The following rate was 0.5 mL/min, column temperature was 25 °C, and the injection volume was 5 μL, analyzing with UV at 445 nm. The lutein, β-carotene, lycopene standards were obtained from Solarbio (Beijing Solarbio Science & Technology Co., Ltd., Beijing, China). The β-cryptoxanthin standard was from Sigma (USA).
RNA extraction and RNA-Seq
Total RNA was extracted from the papaya fleshes or peels using an RNA extraction kit (Dongsheng Biotech, Guangzhou, China). Total RNA of four samples (PE1, PE2, FL1, FL2) was used to prepare cDNA libraries using the Illumina Dynabeads® mRNA DIRECT™ Kit. Then, the cDNA libraries were used for paired-end 125 sequencing using an Illumina Hiseq2500 at Genedenovo Biotechnology Co., Ltd. (Guangzhou, China). In total, four sets of raw reads were obtained, and all sequencing data were deposited in the NCBI Sequence Read Archive (SRA). The raw reads were filtered to remove “dirty” data, including low-quality reads, adaptor sequences, and rRNA reads, to generate “clean” reads. FPKM (Fragments Per Kilobase of transcript per Million mapped reads) was used to measure the transcript abundance of each gene. Those with a fold-change of ≥2 and a false discovery rate (FDR) < 0.05 were considered significantly differentially expressed genes (DEGs) . All of the genes were annotated using the reference papaya genome database (Carica papaya Version1.181: CpGDB181(JGI)) , NCBI non-redundant (Nr) database, the Kyoto Encyclopedia of genes and genomes (KEGG) pathway database. The KEGG enrichment analysis was performed with a Qvalue cut off of 0.05.
Reverse transcription quantitative PCR analysis
For RT-qPCR, oligonucleotide primers were designed according to each gene’s 3′-untranslated region with DNAMAN (Table 2). CpActin was used as the reference gene. RT-qPCR was carried out using SYBR Green-based PCR assay in a Bio-Rad CFX96 Real-Time PCR System (Bio-Rad, USA). Each reaction mix contained 1.0 μL of cDNAs, SYBR Premix ExTaq™ 10 μL, PCR forward primer (10 μmol·L− 1) 0.5 μL, PCR reverse primer (10 μmol·L− 1) 0.5 μL, ddH2O 8.0 μL, to a final volume of 20 μL. The PCR conditions were 95 °C for 3 min, followed by 40 cycles of 95 °C for 15 s, 56 °C for 30 s, and 72 °C for 20 s. Each RT-qPCR analysis was performed in triplicate, and the mean was used for RT-qPCR analysis. The relative expression of the genes was calculated according to the method of 2−△△Ct , and SPSS was used to analyze the data.
Color break stage
Chromoplast-specific lycopene β-cyclase
Differentially expressed Genes
Differential gene expression
The flesh of GS
The flesh of CB
The flesh of HY
The flesh of FY
Fragments Per Kilobase of transcript per Million mapped reads
Full yellow stage
Half yellow stage
Kyoto Encyclopedia of Genes and Genomes
The peel of GS
The peel of CB
The peel of HY
The peel of FY
Reverse transcription quantitative PCR
Sequence Read Archive
Ultra-high Performance Liquid Chromatography
Meyers KJ, Mares JA, Jr RPI, Truitt B, Liu Z, Millen AE, et al. Genetic evidence for role of carotenoids in age-related macular degeneration in the carotenoids in age-related eye disease study (CAREDS). Invest Ophthalmol Vis Sci. 2014;55(1):587.
Sharoni Y, Linnewiel-Hermoni K, Khanin M, Salman H, Veprik A, Danilenko M, et al. Carotenoids and apocarotenoids in cellular signaling related to cancer: a review. Mol Nutr Food Res. 2012;56(2):259–69.
Zhou Y, Tang Y, Tan X, Guo J. Effects of exogenous ABA, GA3 and cell-wall-degrading enzyme activity, carotenoid content in ripening mango fruit. Acta Phytophisiologica Sinica. 1996;22:421–6.
Abushita AA, Daood HG, Biacs PA. Change in carotenoids and antioxidant vitamins in tomato as a function of varietal and technological factors. J Agric Food Chem. 2000;48(6):2075–81.
Kato M, Ikoma Y, Matsumoto H, Sugiura M, Hyodo H, Yano M. Accumulation of carotenoids and expression of carotenoid biosynthetic genes during maturation in citrus fruit. Plant Physiol. 2004;134(2):824–37.
Davey MW, Bergh IVD, Markham R, Swennen R, Keulemans J. Genetic variability in Musa fruit provitamin a carotenoids, lutein and mineral micronutrient contents. Food Chem. 2009;115(3):806–13.
Schweiggert RM, Steingass CB, Mora E, Esquivel P, Carle R. Carotenogenesis and physico-chemical characteristics during maturation of red fleshed papaya fruit (Carica papaya L.). Food Res Int. 2011;44(5):1373–80.
Hadjipieri M, Georgiadou EC, Marin A, Diazmula HM, Goulas V, Fotopoulos V, et al. Metabolic and transcriptional elucidation of the carotenoid biosynthesis pathway in peel and flesh tissue of loquat fruit during on-tree development. BMC Plant Biol. 2017;17(1):102.
Giuliano G, Aquilani R, Dharmapuri S. Metabolic engineering of plant carotenoids. Trends Plant Sci. 2000;5(10):406–9.
Hirschberg J. Carotenoid biosynthesis in flowering plants. Curr Opin Plant Biol. 2001;4(3):210–8.
Paine JA, Shipton CA, Chaggar S, Howells RM, Kennedy MJ, Vernon G, et al. Improving the nutritional value of Golden Rice through increased pro-vitamin a content. Nat Biotechnol. 2005;23(4):482.
Rosati C, Diretto G, Giuliano G. Biosynthesis and engineering of carotenoids and apocarotenoids in plants: state of the art and future prospects. Biotechnol Genet Eng Rev. 2009;26(1):139–62.
Krishna K, Paridhavi M, Patel JA. Review on nutritional, medicinal and pharmacological properties of papaya (Carica papaya Linn.). Natural product radiance. 2008;7(4):364–73.
Yamamoto HY. Biochemistry of the violaxanthin cycle in higher plants. Pure Appl Chem. 1979;51(3):639–48.
Wall MM. Ascorbic acid, vitamin a, and mineral composition of banana (Musa sp.) and papaya (Carica papaya) cultivars grown in Hawaii. J. Food Compos. Anal. 2006;19:434–45.
Saengmanee P, Burns P, Chaisan T, Thaipong K, Siriphanich J. Genetic diversity of genes involved in the carotenoid pathway of Carica papaya L. and their expression during fruit ripening. J. Plant Biochem. Biotechnol. 2017;27(1):1–10.
Blas AL, Ming R, Liu Z, Veatch OJ, Paull RE, Moore PH, et al. Cloning of the papaya chromoplast-specific lycopene beta-cyclase, CpCYC-b, controlling fruit flesh color reveals conserved microsynteny and a recombination hot spot. Plant Physiol. 2010;152(4):2013–22.
Welsch R, Beyer P, Hugueney P, Kleinig H, Lintig JV. Regulation and activation of phytoene synthase, a key enzyme in carotenoid biosynthesis, during photomorphogenesis. Planta. 2000;211(6):846–54.
Cunningham FX, Gantt E. Genes and enzymes of carotenoid biosynthesis in plants. Annu Rev Plant Physiol Plant Mol Biol. 1998;49(1):557–83.
Devitt LC, Fanning K, Dietzgen RG, Holton TA. Isolation and functional characterization of a lycopene β-cyclase gene that controls fruit colour of papaya (Carica papaya L.). J Exp Bot. 2010;61(1):33–9.
Skelton RL, Yu Q, Srinivasan R, Manshardt R, Moore PH, Ming R. Tissue differential expression of lycopene beta-cyclase gene in papaya. Cell Res. 2006;16(8):731–9.
Wang W, Xu J, Xu C. Developmental changes in external color, pigment content and composition in Citrus iyo fruit. Acta Horticulturae Sinica. 2006;33(3):461–5.
Schweiggert RM, Mezger D, Schimpf F, Steingass CB, Carle R. Influence of chromoplast morphology on carotenoid bioaccessibility of carrot, mango, papaya, and tomato. Food Chem. 2012;135(4):2736–42.
Lu P, Wang C, Yin T, Zhong S, Grierson D, Chen K, et al. Cytological and molecular characterization of carotenoid accumulation in normal and high-lycopene mutant oranges. Sci Rep. 2017;7(1):761–76.
Chan-León AC, Estrella-Maldonado H, Dubé P, Fuentes GO, Espadas-Gil F, Talavera CM, et al. The high content of β-carotene present in orange-pulp fruits of Carica papaya L. is not correlated with a high expression of the CpLCY-β2 gene. Food Res Int. 2017;100(Pt 2):45–56.
Armstrong GA, Hearst JE. Carotenoids 2: genetics and molecular biology of carotenoid pigment biosynthesis. FASEB J. 1996;10(2):228–37.
Huner NPA, Öquist G, Sarhan F. Energy balance and acclimation to light and cold. Trends Plant Sci. 1998;3(6):224–30.
Shan W, Kuang J, Chen L, Xie H, Peng H, Xiao Y, et al. Molecular characterization of banana NAC transcription factors and their interactions with ethylene signalling component EIL during fruit ripening. J Exp Bot. 2012;63(14):5171–87.
Zhu M, Chen G, Zhou S, Tu Y, Wang Y, Dong T, et al. A new tomato NAC (NAM/ATAF1/2/CUC2) transcription factor, SlNAC4, functions as a positive regulator of fruit ripening and carotenoid accumulation. Plant Cell Physiol. 2014;55(1):119–35.
Fu CC, Han YC, Qi XY, Shan W, Chen JY, Lu WJ, et al. Papaya CpERF9 acts as a transcriptional repressor of cell-wall-modifying genes CpPME1/2 and CpPG5 involved in fruit ripening. Plant Cell Rep. 2016;35(11):1–12.
Zhu F, Luo T, Liu C, Wang Y, Yang H, Yang W, et al. An R2R3-MYB transcription factor represses the transformation of α- and β-branch carotenoids by negatively regulating expression of CrBCH2 and CrNCED5 in flavedo of Citrus reticulate. New Phytol. 2017;216(1):178–92.
Fu CC, Han YC, Fan ZQ, Chen JY, Chen WX, Lu WJ, et al. The papaya transcription factor CpNAC1 modulates carotenoid biosynthesis through activating phytoene desaturase genes CpPDS2/4 during fruit ripening. J Agric Food Chem. 2016;64(27):5454–63.
Fu CC, Han YC, Kuang JF, Chen JY, Lu WJ. Papaya CpEIN3a and CpNAC2 cooperatively regulate carotenoid biosynthesis related genes CpPDS2/4, CpLCY-e, and CpCHY-b during fruit ripening. Plant Cell Physiol. 2017;58(12):2155–65.
Lichtenthaler HK, Wellburn AR. Determination of total carotenoids and chlorophylls a and b of leaf in different solvents. Analysis. 1983;11(5):591–2.
Yan J, Cai Z, Shen Z, Zhang B, Ruijuan MA, Mingliang YU. Extraction and analytical methods of carotenoids in fruit of yellow flesh peach. Journal of Fruit Science. 2015;32:1267–74.
Lee HS. Characterization of carotenoids in juice of red navel orange (Cara Cara). J. Agric. Food Chem. 2001;49(5):2563–8.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.
Ming R, Hou S, Feng Y, Yu Q, Dionne-Laporte A, Saw JH, et al. The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus). Nature. 2008;452(7190):991–6.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25(4):402–8.
Naik PS, Chanemougasoundharam A, Khurana SMP, Kalloo G. Genetic manipulation of carotenoid pathway in higher plants. Curr Sci. 2003;85(25):1423–30.
We sincerely thank Qi Xin Chen (Genedenovo Biotechnology Co., Ltd., Guangzhou, Beijing) for the technical assistance.
Innovation and Technology Foundation in Fujian Agriculture and Forestry University, China (KFA17355A).
Science Foundation for Doctorate Research in Hebei Normal University of Science and Technology, China (2018YB023).
Availability of data and materials
The raw data had been submitted to the NCBI SRA database, the accessions number is SRA688610.
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.
Figure S1. KEGG graph of carotenoid biosynthetic pathway (FL1-vs-FL2). 18.104.22.168 indicates PDS (evm.TU.supercontig_157.3, fold 4.4); 22.214.171.124 indicates ZDS (evm.TU.supercontig_117.67, fold 4.3); CrtL-b indicates CYCB (evm.TU.supercontig_195.16, fold 1.0); CruA/P indicates LCYB2 (evm.TU.supercontig_132.5, fold − 1.6); CrtR indicates CHYB (evm.TU.supercontig_107.106, fold 3.2); 126.96.36.199 indicates VDE (evm.TU.supercontig_51.78, fold 1.4). (DOCX 51 kb)
Figure S2. KEGG graph of carotenoid biosynthetic pathway (PE1-vs-PE2). 188.8.131.52 indicates PDS (evm.TU.supercontig_157.3, fold 3.9); 184.108.40.206 indicates ZDS (evm.TU.supercontig_117.67, fold 5.3); CrtL-e indicates LCYE (evm.TU.supercontig_28.134, fold − 2.5); CrtL-b indicates CYCB (evm.TU.supercontig_195.16, fold 5.4); LUT1(evm.TU.supercontig_5.131, fold 1.1); CrtR indicates CHYB (evm.TU.supercontig_107.106, fold 3.9); 220.127.116.11 indicates VDE (evm.TU.supercontig_51.78, fold 1.6); 1.14.1390 indicates ZEP (evm.TU.supercontig_55.146, fold 1.4). (DOCX 50 kb)
Figure S3. KEGG graph of carotenoid biosynthetic pathway (PE2-vs-FL2). CrtL-e indicates LCYE (evm.TU.supercontig_28.134, fold − 3.9). (DOCX 50 kb)
Table S1. Carotenoid biosynthesis-related genes. (DOCX 16 kb)
About this article
Cite this article
Shen, Y.H., Yang, F.Y., Lu, B.G. et al. Exploring the differential mechanisms of carotenoid biosynthesis in the yellow peel and red flesh of papaya. BMC Genomics 20, 49 (2019). https://doi.org/10.1186/s12864-018-5388-0