- Research article
- Open Access
Isolation of ripening-related genes from ethylene/1-MCP treated papaya through RNA-seq
BMC Genomics volume 18, Article number: 671 (2017)
Since papaya is a typical climacteric fruit, exogenous ethylene (ETH) applications can induce premature and quicker ripening, while 1-methylcyclopropene (1-MCP) slows down the ripening processes. Differential gene expression in ETH or 1-MCP-treated papaya fruits accounts for the ripening processes. To isolate the key ripening-related genes and better understand fruit ripening mechanisms, transcriptomes of ETH or 1-MCP-treated, and non-treated (Control Group, CG) papaya fruits were sequenced using Illumina Hiseq2500.
A total of 18,648 (1-MCP), 19,093 (CG), and 15,321 (ETH) genes were detected, with the genes detected in the ETH-treatment being the least. This suggests that ETH may inhibit the expression of some genes. Based on the differential gene expression (DGE) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, 53 fruit ripening-related genes were selected: 20 cell wall-related genes, 18 chlorophyll and carotenoid metabolism-related genes, four proteinases and their inhibitors, six plant hormone signal transduction pathway genes, four transcription factors, and one senescence-associated gene. Reverse transcription quantitative PCR (RT-qPCR) analyses confirmed the results of RNA-seq and verified that the expression pattern of six genes is consistent with the fruit senescence process. Based on the expression profiling of genes in carbohydrate metabolic process, chlorophyll metabolism pathway, and carotenoid metabolism pathway, the mechanism of pulp softening and coloration of papaya was deduced and discussed. We illustrate that papaya fruit softening is a complex process with significant cell wall hydrolases, such as pectinases, cellulases, and hemicellulases involved in the process. Exogenous ethylene accelerates the coloration of papaya changing from green to yellow. This is likely due to the inhibition of chlorophyll biosynthesis and the α-branch of carotenoid metabolism. Chy-b may play an important role in the yellow color of papaya fruit.
Comparing the differential gene expression in ETH/1-MCP-treated papaya using RNA-seq is a sound approach to isolate ripening-related genes. The results of this study can improve our understanding of papaya fruit ripening molecular mechanism and reveal candidate fruit ripening-related genes for further research.
Papaya (Carica papaya L.) is one of the most important fruit crops cultivated in tropical and sub-tropical areas and the ripe fruit has a soft and sweet pulp with high amounts of vitamin C, vitamin A, and carotenes. Consumption of fruits rich in vitamin C, carotene, and vitamin E have been associated with a reduced risk of colon cancer. In addition, papaya is a rich source of the digestive enzyme papain, which is used widely in the textile, food, animal feed, chemicals and pharmaceutical industries [1,2,3]. Global papaya production has grown significantly over the last few years, and papaya is now ranked as the third most popular tropical fruit, behind mango and pineapple . Because papaya is a typical climacteric fruit, it ripens very quickly with striking color changes, substantial pulp softening, and rotting after harvesting. One of the major problems faced by the global papaya industry is significant post-harvest losses throughout the marketing chain. Postharvest losses up to 75% have been reported for papaya fruits shipped from Hawaii to USA mainland . In the Southeast Asia region, post-harvest losses of papaya ranged from 30 to 60%.
Ethylene (ETH) plays a critical role in determining the timing of ripening. Today, supplemental ethylene is commonly used to speed up the ripening of bananas, avocados, mangos in the marketplace. A downside of this treatment is that ethylene shortens the shelf life of many fruits by hastening fruit ripening . Exogenous ethylene can prematurely induce greater endogenous ethylene production, and quicker ripening in climacteric fruit. Exogenous ethylene applications (100 μL·L−1) stimulated papaya skin degreening and yellowing, and flesh softening . Researchers have developed several ways to impair the ethylene signaling pathway, including inhibiting ethylene synthesis and perception. Inhibitors of ethylene perception include compounds that have a similar shape to ethylene, but are not able to elicit the ethylene response. A good example of ethylene perception inhibitors is 1-methylcyclopropene (1-MCP). The mechanism of action of 1-MCP involves its tight binding to the ethylene receptor in plants, thereby preventing the binding of ethylene and blocking the effects of ethylene . 1-MCP is used commercially to slow down the ripening of fruits, such as apples, kiwifruits, tomatoes, bananas, plums, persimmons, avocados. 1-MCP can also slow down the ripening of papaya fruit [9, 10]. This is the result of dynamic processes that involve in a complex series of molecular and biochemical changes under genetic regulation.
To better understand the mechanisms of papaya fruit ripening, numerous studies have focused on the analysis of transcript, protein, and metabolite levels in papaya fruits. Using an Arabidopsis thaliana microarray, 414 ripening-related genes were identified, and some transcription factors were found in papaya fruit . Twenty-seven protein spots showing differences in abundance during papaya ripening were successfully identified using the 2-DE analysis . Although some studies on papaya ripening have been conducted, little is known about the genetic control of ripening due to technical limitations. RNA-seq is a good method to examine the total RNA levels in different samples. Exogenous ethylene stimulates papaya ripening, while 1-MCP inhibits the ripening progression. Differential gene expression induced by the different treatments accounts for these observations. Therefore, screening for differentially expressed genes will be helpful for further elucidating the fruit ripening molecular mechanism. This study will analyze transcript levels in papaya fruits after 1-MCP or ethylene treatment to isolate the key ripening-related candidate genes, as well attempting to enhance our understanding of papaya fruit ripening molecular mechanism.
To identify candidate ripening-related genes in papaya, we conducted a large-scale transcriptome analysis of ETH-treated, 1-MCP-treated and untreated papaya fruit using an RNA-seq based approach. Exogenous ethylene accelerates climacteric fruit ripening, whereas 1-MCP slows it down. The effects of exogenous ethylene and 1-MCP treatments on skin color, total soluble solids (TSS), and pulp texture are shown in Fig. 1. In 1 day the pulp firmness of ETH-treated papaya changed from 14 to 3.93 105 Pa, whereas 1-MCP-treated fruits retained pulp firmness at high levels for the entire duration of the experiment (Fig. 1). 1-MCP also delayed papaya fruit skin coloration, and fruits treated with 1-MCP were still green 6 days after treatment (Fig. 1a). The TSS changing trend of 1-MCP-treated papaya was significantly different compared to that of the control and ETH-treatments. The TSS was highest at 24 h and then decreased in control and ETH-treatments, while TSS kept a low level and increased a little in 6 d in 1-MCP-treatments. Papaya fruits treated with 1-MCP showed a significant delay in the softening, color changing, and the formation of soluble solids, In contrast the ETH-treated papayas showed quick softening and color changing, with the differences in firmness and TSS at 24 h among three treatments reaching significant levels (Fig. 1b, c). Therefore, three samples collected at 24 h after treatments were used for RNA extraction, cDNA libraries construction, and sequencing.
The main results of RNA-seq are showed in Table 1 and Fig. 2. Table 1 shows the clean reads number and the gene number of the three samples. After removing low-quality reads, adaptor sequences, and rRNA reads, we obtained 43,873,036 (1-MCP), 65,149,940 (CG), and 33,805,002 (ETH) clean reads. Clean reads were then mapped to the papaya reference genome. There were 18,648 (1-MCP), 19,093 (CG), and 15,321 (ETH) genes (including 1127 (1-MCP), 1138 (CG), 1009 (ETH) new genes) were detected, which the genes were detected in the ETH-treatment was the least. The data showed that ETH inhibited the expression of some genes. Compared with the untreated papaya, 760 genes were up-regulated and 4753 genes were down-regulated in ETH-treatment; in the 1-MCP-treated papaya, there were only 608 genes up-regulated and 738 genes down-regulated (Fig. 2). This demonstrates that ETH inhibited more genes than 1-MCP. This may because ETH is an important hormone in fruit ripening, it can affect a lot of genes’ expression.
DEGs function annotation
Papaya genes were annotated using the Gene Ontology database and genes were classed into three categories (Additional file 1: Figure S1). Among the three categories, the terms of cellular progress and metabolic progress were observed to occur most frequently in the ontology of biological process; while the terms of cell and cell part were observed to occur most frequently in the ontology of cellular component; and the terms of binding, catalytic activity were observed to occur most frequently in the ontology of molecular function. In comparison with the control treatment, the terms of cell killing, reproductive process, reproduction, multi-organism process, extracellular region, and nutrient reservoir activity were up-regulated with ETH-treatment, while the terms of extracellular region, and nutrient reservoir activity were down-regulated with 1-MCP-treatment. This implied that ethylene may accelerate the reproductive process, the induction of death, and nutritious substrate storage. On the contrary, 1-MCP may inhibit the expression of genes in the extracellular region, or genes relating to nutritious substrate storage.
In addition, the pathway enrichment analysis was carried out using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. The results showed that DEGs participating in the pathways of pentose and glucuronate interconversions and the porphyrin and chlorophyll metabolism were significantly enriched in the ETH-treated fruits; while DEGs that participated in the pathways of plant-pathogen interaction, plant hormone signal transduction, and diterpenoid biosynthesis, pentose and glucuronate interconversions were significantly enriched in the 1-MCP-treated fruits (Fig. 3). Genes involved in the pentose and glucuronate interconversions pathway are shown in Additional file 2: Figure S2. The 22.214.171.124 in the red frame refers to polygalacturonase (PG), which is the key gene in the pentose and glucuronate interconversions pathway. This gene encodes for a very important enzyme that participates in papaya softening . PG was up-regulated significantly response to ETH treatment and down-regulated after 1-MCP treatment.
Ripening-related genes selection
According to the DEGs expression and the KEGG pathway enrichment, 53 genes were selected: 20 cell wall-related genes, 18 chlorophyll and carotenoid metabolism-related genes, four proteinases and their inhibitors, six plant hormone signal transduction pathway genes, four transcription factors, and one senescence-associated gene (Additional files 3: Table S1 and Fig. 4).
Among the cell wall-related genes, except pectinesterase/pectinesterase inhibitor (PMIS), the expression patterns of polygalacturonase (PG), beta-galactosidase (GAL-B), pectate lyase (PL), pectin methylesterase (PME), beta-glucosidase (GLU-B), xyloglucan endotransglucosylase (XTH), endoglucanase 8-like (EGase), endoxylanase (EXY1), beta-D-xylosidase 5 (XYL), sucrose synthase 4 (SUS4), sugar transporter (STP), UDP-galactose transporter 3 (UTR3), beta-fructofuranosidase (BFF), and expansin A (EXPA) were the same. Expression of these genes in the ETH-treatment was higher than that in the control treatment, and lowest in the 1-MCP treatment. PG, GAL-B, PL, and PME are involved in the degradation of pectin. GLU-B, XTH, EGase, EXY1, XYL may be involved in the degradation of cellulose and semi-cellulose. This suggested that ETH induced the expression of PG, GAL-B, PL, PME, GLU-B, XTH, EGase, EXY1, and XYL genes, then accelerated the degradation of pectin, cellulose, and hemi-cellulose, and therefore ethylene treatment results in the rapid softening of papaya fruit. Conversely, 1-MCP inhibited the expression of these genes and the fruits retained pulp firmness. These genes may play important roles in papaya cell wall degradation and fruit softening. PMIS is active against plant PMEs, PMIS showed a very low expression level in ETH-treated papaya. ETH induced the expression of SUS4, STP, UTR3, and BFF genes, promoted the synthesis of sucrose and fructose, and the transport of sugar. Expansins are non-enzymatic proteins found in the plant cell wall, with important roles in plant cell growth, fruit softening, and meristem function. In our study, an expansin A gene showed lower expression level in 1-MCP-treated papaya than that in ETH-treatment and control.
The color of the papaya fruit is determined largely by the presence of chlorophyll and carotenoid . Ethylene applications affect papaya fruit coloration causing it to change from green to yellow. Fig. 3 showed that DEGs involved in the pathways of the porphyrin and chlorophyll metabolism were enriched in the ETH-treatment, Additional file 5: Figure S3 and Additional file 3: Table S1 showed that ETH inhibited chlorophyll metabolism. The major chlorophyll biosynthesis enzymes, including ferretin 1 (FER1), hydroxymethylbilane synthase (HEMC), protoporphyrinogen oxidase (HEMG), magnesium-chelatase subunit chlH (ABAR), magnesium-chelatase subunit ChlI (CH-42), magnesium-protoporphyrin IX methyltransferase (CHLM), and protochlorophyllide oxidoreductase (PORA) were highly inhibited in ETH-treated samples, where chlorophyll biosynthesis was inhibited. Chlorophyll degradation enzymes such as NAD (P)-binding rossmann-fold superfamily protein (NYC1), chlorophyllase 2 (CLH2), pheophorbide a oxygenase, chloroplastic-like (ACD1) were also inhibited. This is indicative that the synthesis and degradation of chlorophyll both decreased in ETH-treated papaya. On the contrary, the expression of ferretin 1 in 1-MCP-treated papaya fruits was about two fold higher than that found in untreated papaya (Additional file 3: Table S1).
In the carotenoid biosynthesis pathway, lycopene beta/epsilon cyclase protein isoform 2 (LCY-B/E) and β-carotene hydroxylase (CHY-B) were enriched in KEGG pathway enrichment (Additional file 4: Figure S4). ETH treatment clearly up-regulated the expression of β-carotene hydroxylase gene, and decreased the expression of LCY-B/E (Additional file 3: Table S1). Other key carotenoid synthesis genes including phytoene synthase (PSY), phytoene desaturase 3 (PDS), zeta-carotene desaturase (ZDS), chromoplast-specific lycopene β-cyclase (CYC-B), lycopene beta-cyclase (LCY-B) were also detected, but were not significantly differentially expressed among treatments.
Protein degradation is one of the most significant features of senescence, and a lot of genes up-regulated during senescence are proteases . Several lines of research suggest that serine proteases are strongly associated with senescence in wheat [15, 16], and barley . Cysteine proteases (CPs) are the most abundant class of up-regulated proteases during natural or induced senescence [14, 18]. In our study, we identified one subtilisin-like serine endopeptidase family protein and one cysteine proteinases superfamily protein, which were up-regulated in response to ETH-treatment and down-regulated in 1-MCP-treated papaya. In addition, the expression of one serine protease inhibitor and one cysteine proteinase inhibitor were greatly inhibited by ETH (Additional file 3: Table S1).
KEGG pathway enrichment analysis of the DEGs indicated that genes encoding the plant hormone signaling components were significantly enriched in the 1-MCP-treated samples. Plant hormones regulate a wide range of processes, including the development and ripening of fruit. Ethylene is an important hormone controlling papaya fruit ripening. In ethylene signal transduction pathway, expression of the ethylene receptor ethylene-insensitive 4-like protein (EIN4) was inhibited by 1-MCP, and the expression of EIN3 binding F-box 2(EBF2) increased in response to 1-MCP-treatment. In Arabidopsis, overexpression of F-box genes display ethylene insensitivity . The high expression levels of EBF2 in the 1-MCP-treatment may reduce the sensitivity to ethylene. Genes in ABA, cytokinin, and auxin signal transduction pathway, with different expression levels in three treatments were also found.
We also selected several transcription factors and one senescence-associated gene through gene different expression. One NAC domain protein and two WRKY DNA-binding proteins showed an obvious decrease, the other NAC and senescence-associated gene increased in expression with application of ethylene.
Confirming genes expression using RT-qPCR
To confirm the accuracy and reproducibility of the transcriptome analysis results, 18 genes were selected for RT-qPCR validation. RNA samples from the 24 h ETH-treated papaya fruit were used as templates. Primers of the candidate genes are shown in Table 2. Except FER1 (evm.TU.supercontig_5.220) and EBF2 (evm.TU.supercontig_750.1), the expression profiles of the other 16 candidate genes determining using RT-qPCR data were consistent with the RNA-seq results (Fig. 5), thus confirming our transcriptome analysis.
Transcription levels of six selected genes during different processing
Six fruit ripening-related candidate genes were selected and their transcription levels in ETH-treated, 1-MCP-treated, and non-treated papaya fruits were analyzed by RT-qPCR (Fig. 6). RT-qPCR analysis revealed that the expression patterns of XTH, GLU-B, and PL were similar. The expression of these genes was highly induced by ethylene and inhibited by 1-MCP. In addition, the expression pattern is consistent with the fruit senescence process and fruit softening, indicating that the three genes are involved in fruit softening in papaya. PMIS inhibits the expression of PMEs, which are known to catalyze the demethoxylation of partial purification and characterization of pectinmethylesterase . Thus, during fruit ripening, PME may play an important role in determining the extent of pectin degradation. PMIS showed very low expression levels in ETH-treated papaya, while PMIS retained a high level in 1-MCP-treatment in 24 h and 3 d. This indicated that high levels of expression of PMIS inhibited the expression of PME in 1-MCP-treated papaya, thus preventing pectin degradation. The expression of one proteinase gene was also analyzed. CPSP showed extremely up-regulated and a high expression level maintained in ETH-treatment. The expression of CPSP increased during the process of fruit ripening and reached the highest level at 6d in 1-MCP-treated and control papayas. CHLM is one of the key chlorophyll biosynthesis enzymes, the expression of CHLM was maintained at a lower level in ETH-treated and control papayas than the 1-MCP-treatment.
To help in elucidating the molecular basis of ripening process in papaya, 24 h ETH-treated, 1-MCP-treated and untreated papaya fruits were analyzed by RNA-seq. 18,648 (1-MCP), 19,093 (CG), and 15,321 (ETH) genes were detected, with thousands of genes showing differential expression. According to the KEGG enrichment analysis and differential gene expression analysis, about 50 ripening-related genes were isolated. Because fruit color and firmness are the major fruit quality parameters, we focused on the analysis of genes related to fruit softening and coloration.
Pulp softening renders the fruit better for eating but harder for storage and transport. Firmness is one of the most important fruit quality features characterizing fruit quality. Papaya is a soft fruit with a short postharvest shelf life. The primary cell wall of fruits is a structure mainly composed of polysaccharides (pectin, cellulose, and hemi-cellulose). During papaya fruits ripening, pectin, cellulose, and hemi-cellulose depolymerisation increase [21, 22]; endo- and exo-polygalacturonase (PG), pectin methyl esterase (PME), glucanase, galactosidase, and xylanase have been detected in ripening papaya fruits [23,24,25]. In our study, several cell wall degrading candidate genes were selected, such as the genes of polygalacturonase, beta-galactosidase, pectin methyl esterase, pectate lyase, xyloglucan endotransglucosylase, endoxylanase, endoglucanase 8-like. Our results are consistent with the previous research.
Among the cell wall degrading candidate genes, PGs, GAL-B, PME, PL are important pectinolytic enzymes which are responsible for the solubilisation of pectins during papaya ripening. PGs and GAL-B play a central role in pectin solubilization during papaya fruit ripening [26, 27]. Three PG genes were detected, the expression patterns were same but the expression abundance among the three genes were different. The fpkm of Pg1 in ETH-treatment is 15,642.8, while that of Pg2 and Pg3 are only 12.09 and 32.85. There were more than 10 beta-galactosidase genes detected, only Gal-b (evm.TU.supercontig_93.34) was abundantly expressed in ethylene-treated papaya and was inhibited by 1-MCP, with the expression pattern being consistent with fruit softening . Because more abundant gene transcripts account for more enzymatic activity participating in softening, Pg1(evm.TU.supercontig_250.6) and Gal-b (evm.TU.supercontig_93.34) are the major cell wall degrading genes involved in the solubilization of pectins.
Several genes that could be involved in the degradation of cellulose and semi-cellulose were also selected, including Glu-b, Xth, Egase, Exy1, and Xyl. Endoglucanase and beta-glucosidase are key enzymatic component involved in cellulose hydrolysis . Exy1, Xyl , and Xth  are important genes participating in the degradation of hemicelluloses, thus playing an important role in plant cell wall remodeling. Some of the genes participating in fruit softening in other species have been identified, including the beta-glucosidase gene in sweet cherry , FaXyl1 in strawberry , Ma-Xth2 and Ma-Xth10 in apple, Ad-Xth4,5 and 7 in kiwifruit .
Cell wall disassembly during fruit ripening is a cooperative process involving the coordinated expression of multiple genes. The expression of sucrose synthase 4, sugar transporter, beta-fructofuranosidase genes were highly up-regulated in ETH-treated papaya, which had the highest TSS and lowest firmness. This suggested that ETH improves the rate of sugar synthesis, transport and degradation. We also found that the transcripts of the expansin A gene were abundant in the ETH-treatment sample. Although expansin is non-enzymatic protein, it is well known to play important roles in fruit softening [32, 33]. The papaya genome contains at least 15 expansin A genes , only the expression of expansin A (evm.TU.supercontig_1.419) was associated with firmness changes, indicating that this gene may be involved in the papaya fruit softening.
To confirm our hypotheses, the expression of Xth30, Glu-b, and Pl genes was measured in ETH-, 1-MCP-, and Non-treated papaya fruits by RT-qPCR. We found that the expression of three genes was ripening-related and correlated with the variation in softening patterns of different treatments. Therefore, we suggested that the high expression levels of Xth30, Glu-b, and Pl genes may be related to the degradation of the papaya fruit cell wall. In general, besides the cell wall-related genes that have been cloned (e.g. Pg, Gal-b, Pl, Xth32, and Exy1), several novel genes (e.g. Pme, Pmis, Glu-b, Egase, Xyl, Expa) were selected for further analysis. We illustrate that papaya fruit softening is a complex process with significant cell wall hydrolases, such as pectinases, cellulases, and hemicellulases participating in the process. Among all the cell wall-related genes, Pg1 was strongly induced by ETH and was highly inhibited by 1-MCP, and the expression abundance was very high. Therefore, Pg1 may play the most important role in the fruit softening of papaya.
Exogenous ethylene applications accelerate the climacteric fruit ripening process, with the most visible change being the color of exocarp turning from green to yellow during papaya fruit ripening. Papaya fruits treated with 1 μL·L−1 1-MCP retained the green color 6 days post treatment. The color of the papaya fruit is determined largely by the presence of chlorophyll and carotenoids . Typically, young fruits are green, due to the high contents of chlorophyll in young fruits. Then a steady increase in carotenoids and simultaneous rapid loss of chlorophyll causes a color shift from green to yellow during fruit ripening [35, 36]. We analyzed the genes involved in chlorophyll and carotenoid metabolism. ETH highly inhibited the expression of several chlorophyll biosynthesis enzyme genes such as Fer1, Hemc, Hemg, Abar, Ch-42, Chlm, and Pora, therefore inhibiting chlorophyll biosynthesis.
‘Da Qing No 7’ is a red-fleshed papaya. The pulp of full-ripe fruit is red and the peel is yellow. β-cryptoxanthin, lycopene, and β-carotene are the major carotenoids in these red-fleshed papaya fruits [37, 38]. Subsequent red-fleshed fruit maturation leads to a gradual accumulation of β-cryptoxanthin and total lycopene . Thus, we speculated that the accumulation of β-cryptoxanthin may be the main reason for the yellow peel, β-carotene hydroxylase (Chy-b) was enriched in the KEGG pathway, and the expression of Chy-b was highly inhibited by 1-MCP and boosted by ETH. β-carotene hydroxylase is a key enzyme in the pathway of carotenoid biosynthesis in plants, which catalyses the conversion of β-carotene to β-cryptoxanthin [40, 41]. 1-MCP inhibited the expression of Chy-b, while ETH increased the expression of this gene, illustrating that CHY-B is a key enzyme in β-cryptoxanthin synthesis and may be responsible for the yellow color of papaya fruit.
The other key carotenoid synthesis genes including Psy, Pds, Zds, and Cyc were also found to be expressed in all treatments. There were no significant differences in the expression of Psy, Pds, Zds, and Cyc among the three treatments. Lcy-b/e was enriched in the KEGG pathway. Cyclization of lycopene is a key branch point in the carotenoid biosynthetic pathway, the biosynthetic pathway of carotenoids splits at the level of lycopene into the α-branch and the β-branch . As we can see from Additional file 6 Figure S5, in α-branch, lycopene ε-cyclase catalyzes the cyclization of lycopene to form α-carotene; In β-branch, lycopene β-cyclase catalyzes the cyclization of lycopene to form β-carotene . The expression of Lcy-b/e in ETH-treated papaya has a significant lower transcriptional level than that in untreated papaya, the expression of the lycopene ε-cyclase genes were down-regulated by ethylene in orange too , this suggests that ethylene may inhibit the α-branch and β-branch may be the main synthetic pathway in the later stage of fruit ripening . Besides Lcy-b/e and Lcy-b, chromoplast-specific lycopene beta-cyclase (Cyc-b) was detected as well. A 2-bp insertion is present in the coding region of the recessive red flesh allele resulting in a frame-shift mutation and a premature stop codon, so as to accumulate a large amount of lycopene which is responsible for the red flesh . In ‘Da Qing No 7’, lycopene beta-cyclase (LCY-B, evm.TU.supercontig_5.198) is in β-branch and may be responsible for the biosynthesis of β-carotene, which are then converted to β-cryptoxanthin, thus contributing to the yellow peel of papaya. The expression pattern of Chy-b is consistent with the β-cryptoxanthin accumulation. Thus, we speculated that the accumulation of β-cryptoxanthin may be responsible for the yellow color of papaya fruit.
Fruit ripening is a complex, developmentally regulated process requiring the participation of numerous genes. Exogenous ethephon (0.5 g·L−1) treatment stimulated papaya ripening, while 1-MCP (1 μL·L−1) treatment inhibited ripening progression and made the fruits retain a “rubbery” texture . Differential gene expression induced by the different treatments accounts for these observations. Therefore, comparing the differential gene expression in ETH-treated and 1-MCP-treated papaya using RNA-seq is a sound approach to isolate papaya ripening-related genes. Numerous cell wall-related genes, color-related genes, proteinases, plant hormone signal transduction pathway genes, transcription factors, and a senescence-associated gene were isolated in our study. Among these genes, several genes had already been isolated and characterized as players in papaya ripening, corroborating our findings. These include the Pg and Gal-b, that play a central role in pectin solubilization during papaya fruit ripening [26, 27]. RT-qPCR were used to analyze the transcription levels of several selected genes under different treatments. The results also confirmed the RNA-seq data and verified that the expression pattern of the six analyzed genes is consistent with the fruit senescence process. These findings strongly suggest that these genes may be involved in the ripening of papaya fruit. Nevertheless, the specific role of these genes in fruit ripening remains unclear, and further work is required to accurately characterize the function of these genes in fruit ripening. The mechanisms of pulp softening and coloration of papaya were also deduced and discussed in the paper. We illustrate that papaya fruit softening is a complex process with significant cell wall hydrolases, such as pectinases, cellulases, and hemicellulases involved in the process. Pg1 may play the most important role in the fruit softening of papaya. Exogenous ethylene accelerates the coloration of papaya changing from green to yellow. This is likely due to the inhibition of chlorophyll biosynthesis and the α-branch of carotenoid metabolism. Chy-b may play an important role in the yellow color of papaya fruit. Our results deepen our understanding of the molecular mechanisms regulating papaya fruit ripening, and provide candidate fruit ripening-related genes for further research.
Plant materials and RNA preparation
Papaya fruits (C. papaya L. cv. ‘Daqing No.7’) at the mature-green stage were harvested from a local commercial plantation in Zhangzhou, China. The healthy fruits of similar size, shape, and maturity were separated into three groups and treated at 25 °C. Thirty-six fruits were incubated with 1 μL·L−1 of 1-methylcyclopropene (1-MCP) gas for 18 h in a sealed box; 36 fruits were dipped into 0.5 g·L−1 of ethephon solution for 3 min, then dried and put in a sealed box for 2 h; the 36 control fruits (Control Group, CG) were dipped into water for 3 min, then dried and put in a sealed box for 2 h, with three replicates. After treatments, all fruits were stored at 25 °C and allowed to ripen. Fruits were taken randomly at 6 h, 24 h, 3 d, and 6 d after treatments. The firmness of papaya was measured with a GY-3 firmness meter and the total soluble solids (TSS) were measured with a pocket refractometer. Three biological replicates were taken. Three fruits were peeled, seeds were removed, and the flesh was cut into some pieces. The pieces of papaya pulp were mixed, frozen in liquid nitrogen, and stored in - 80 °C. Total RNA was extracted from the papaya pulp using an RNA extraction kit (DongSheng Biological Technology Ltd., Guangzhou, P. R. China).
Sequencing, assembly and gene annotation
Total RNA of three samples (CG 24 h, ETH 24 h, 1-MCP 24 h) 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, three 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. The estimate of gene expression and identification of differentially expressed genes (DEGs) were conducted using a modified method described previously . 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 significant 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 were performed with a Qvalue cut off of 0.05.
For RT-qPCR, oligonucleotide primers were designed according to each gene’s 3′-untranslated region with DNAMAN (Table 2). CpEIF 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 30s, and 72 °C for 35 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.
Magnesium-chelatase subunit chlH
Pheophorbide a oxygenase, chloroplastic-like
Magnesium-chelatase subunit ChlI
Magnesium-protoporphyrin IX methyltransferase
Chromoplast-specific lycopene β-cyclase
Differentially expressed genes
Differential gene expression
Fragments Per Kilobase of transcript per Million mapped reads
Kyoto Encyclopedia of Genes and Genomes
Lycopene beta/epsilon cyclase protein isoform 2
NAD (P)-binding rossmann-fold superfamily protein
Protochlorophyllide oxidoreductase A
Reverse transcription quantitative PCR
Sucrose synthase 4
Total soluble solids
UDP-galactose transporter 3
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.
Fabi J, Peroni F, Gomez M. Papaya, mango and guava fruit metabolism during ripening: postharvest changes affecting tropical fruit nutritional content and quality. Fresh Produce. 2010;4:56–66.
Maisarah A, Asmah R, Fauziah O. Proximate analysis, antioxidant and antiproliferative activities of different parts of Carica papaya. Journal of Nutrition and Food Science. 2014;4:267.
Evans EA, Ballen FH. An overview of global papaya production, trade, and consumption. 2015. http://edis.ifas.ufl.edu/fe913. Accessed Aug 2015.
Paull RE, Nishijima W, Reyes M, Cavaletto C. Postharvest handling and losses during marketing of papaya (Carica papaya L.). Postharvest Biology & Technology. 1997;11(3):165–79.
Grierson D. Ethylene and the control of fruit ripening. In: Seymour GB, Poole M, Giovannoni JJ, Tucker GA, editors. The molecular biology and biochemistry of fruit ripening. New York City: Wiley; 2013. p. 43–6.
An J-F, Paull RE. Storage temperature and ethylene influence on ripening of papaya fruit. J Am Soc Horticultural Sci. 1990;115(6):949–53.
Serek M, Tamari G, Sisler E, Borochov A. Inhibition of ethylene-induced cellular senescence symptoms by 1-methylcyclopropene, a new inhibitor of ethylene action. Physiol Plantarum. 1995;94(2):229–32.
Huerta-Ocampo JA, Osuna-Castro JA, Lino-Lopez GJ, Barrera-Pacheco A, Mendoza-Hernández G, De León-Rodríguez A, et al. Proteomic analysis of differentially accumulated proteins during ripening and in response to 1-MCP in papaya fruit. J Proteome. 2012;75(7):2160–9.
Manenoi A, Bayogan ERV, Thumdee S, Paull RE. Utility of 1-methylcyclopropene as a papaya postharvest treatment. Postharvest Biology & Technology. 2007;44(1):55–62.
Fabi JP, Seymour GB, Graham NS, Broadley MR, May ST, Lajolo FM, et al. Analysis of ripening-related gene expression in papaya using an Arabidopsis-based microarray. BMC Plant Biol. 2012;12(1):242.
Fabi JP, Broetto SG, da Silva SLGL, Zhong S, Lajolo FM, JRO d N. Analysis of papaya cell wall-related genes during fruit ripening indicates a central role of polygalacturonases during pulp softening. PLoS One. 2014;9(8):e105685.
Yamamoto H. Comparison of the carotenoids in yellow-and red-fleshed Carica papaya. Nature. 1964;201(4923):1049–50.
Guo Y, Cai Z, Gan S. Transcriptome of Arabidopsis leaf senescence. Plant, Cell Envi. 2004;27(5):521–49.
Chauhan S, Srivalli S, Nautiyal A, Khanna-Chopra R. Wheat cultivars differing in heat tolerance show a differential response to monocarpic senescence under high-temperature stress and the involvement of serine proteases. Photosynthetica. 2009;47(4):536–47.
Roberts IN, Caputo C, Kade M, Criado MV, Barneix AJ. Subtilisin-like serine proteases involved in N remobilization during grain filling in wheat. Acta Physiol Plant. 2011;33(5):1997–2001.
Parrott DL, McInnerney K, Feller U, Fischer AM. Steam-girdling of barley (Hordeum vulgare) leaves leads to carbohydrate accumulation and accelerated leaf senescence, facilitating transcriptomic analysis of senescence-associated genes. New Phytol. 2007;176(1):56–69.
Bhalerao R, Keskitalo J, Sterky F, Erlandsson R, Björkbacka H, Birve SJ, et al. Gene expression in autumn leaves. Plant Physiol. 2003;131(2):430–42.
Guo H, Ecker JR. Plant responses to ethylene gas are mediated by SCF EBF1/EBF2-dependent proteolysis of EIN3 transcription factor. Cell. 2003;115(6):667–77.
Cameron RG, Grohmann K. Partial purification and thermal characterization of pectinmethylesterase from red grapefruit finisher pulp. J Food Sci. 1995;60(4):821–5.
Kondo S, Sae-Lee K, Kanlayanarat S. Xyloglucan and polyuronide in the cell wall of papaya fruit during development and storage. Acta Hortic. 2010;875:351–6.
Lazan H, Selamat M, Ali ZM. β-Galactosidase, polygalacturonase and pectinesterase in differential softening and cell wall modification during papaya fruit ripening. Physiol Plantarum. 1995;95(1):106–12.
Paull RE, Chen NJ. Postharvest variation in cell wall-degrading enzymes of papaya (Carica papaya L.) during fruit ripening. Plant Physiol. 1983;72(2):382–5.
Chen NJ, Paull RE. Endoxylanase expressed during papaya fruit ripening: purification, cloning and characterization. Functional Plant Biol. 2003;30(4):433–41.
Thumdee S, Manenoi A, Chen NJ, Paull RE. Papaya fruit softening: role of hydrolases. Tropical Plant Biol. 2010;3(2):98–109.
Fabi JP, Cordenunsi BR, Seymour GB, Lajolo FM, JRO d N. Molecular cloning and characterization of a ripening-induced polygalacturonase related to papaya fruit softening. Plant Physiol Biochemist. 2009;47(11):1075–81.
Othman R, Choo T, Ali ZM, Zainal Z, Lazan H. A full-length beta-galactosidase cDNA sequence from ripening papaya. Plant Physiol. 1998;118:1102.
Singhania RR, Patel AK, Sukumaran RK, Larroche C, Pandey A. Role and significance of beta-glucosidases in the hydrolysis of cellulose for bioethanol production. Bioresour Technol. 2013;127(1):500–7.
Martínez GA, Chaves AR, Civello PM. β-Xylosidase activity and expression of a β-xylosidase gene during strawberry fruit ripening. Plant Physiol Biochemist. 2004;42(2):89–96.
Atkinson RG, Johnston SL, Yauk YK, Sharma NN, Schröder R. Analysis of xyloglucan endotransglucosylase/hydrolase (XTH) gene families in kiwifruit and apple. Postharvest Biol Technol. 2009;51(2):149–57.
Gerardi C, Blando F, Santino A, Zacheo G. Purification and characterisation of a β-glucosidase abundantly expressed in ripe sweet cherry (Prunus avium L.) fruit. Plant Sci. 2001;160(5):795–805.
Brummell DA, Harpster MH, Civello PM, Palys JM, Bennett AB, Dunsmuir P. Modification of expansin protein abundance in tomato fruit alters softening and cell wall polymer metabolism during ripening. Plant Cell. 1999;11(11):2203–16.
Hayama H, Ito A, Moriguchi T, Kashimura Y. Identification of a new expansin gene closely associated with peach fruit softening. Postharvest Biology Technol. 2003;29(1):1–10.
Paull RE, Wu P, Chen NJ. Genomics of papaya fruit development and ripening. In: Ming R, Moore PH, editors. Genetics and genomics of papaya. New York: Springer; 2014. p. 244–5.
Sancho GG, Yahia EM, González-Aguilar GA. Identification and quantification of phenols, carotenoids, and vitamin C from papaya ( Carica papaya L., cv. Maradol) fruit determined by HPLC-DAD-MS/MS-ESI. Food Res Int. 2011;44(5):1284–91.
Vasquez-Caicedo AL, Heller A, Neidhart S, Carle R. Chromoplast morphology and beta-carotene accumulation during postharvest ripening of mango cv. ‘Tommy Atkins. J Agric Food Chemist. 2006;54(16):5769–76.
Chandrika UG, Jansz ER, Wickramasinghe S, Warnasuriya ND. Carotenoids in yellow-and red-fleshed papaya (Carica papaya L). J Sci Food Agric. 2003;83(12):1279–82.
Wall MM. Ascorbic acid, vitamin a, and mineral composition of banana (Musa sp.) and papaya (Carica papaya) cultivars grown in Hawaii. J Food Composition Analy. 2006;19(5):434–45.
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.
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. Current Opinion in Plant Biol. 2001;4(3):210–8.
Cunningham FX, Gantt E. Genes and enzymes of carotenoid biosynthesis in plants. Annu Rev Plant Physiol Plant Mol Biol. 1998;49(1):557–83.
Zhang L, Ma G, Shirai Y, Kato M, Yamawaki K, Ikoma Y, Matsumoto H. Expression and functional analysis of two lycopene β-cyclases from citrus fruits. Planta. 2012;236(4):1315–25.
Rodrigo MJ, Zacarias L. Effect of postharvest ethylene treatment on carotenoid accumulation and the expression of carotenoid biosynthetic genes in the flavedo of orange (Citrus sinensis L. Osbeck) fruit. Postharvest Biol Technol. 2007;43(1):14–22.
Blas AL, Ming R, Liu Z, Veatch OJ, Paull RE, Moore PH, et al. Cloning of the papaya chromoplast-specific lycopene β-cyclase, CpCYC-b, controlling fruit flesh color reveals conserved microsynteny and a recombination hot spot. Plant Physiol. 2010;152(4):2013–22.
Jeong J, Huber DJ. Suppression of avocado (Persea americana mill.) fruit softening and changes in cell wall matrix polysaccharides and enzyme activities: differential responses to 1-MCP and delayed ethylene application. J Am Soc Horticult Sci. 2004;129(5):752–9.
Audic S, Claverie J-M. The significance of digital gene expression profiles. Gen Res. 1997;7(10):986–95.
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.
Zhu X, Li X, Chen W, Chen J, Lu W, Chen L, et al. Evaluation of new reference genes in papaya for accurate transcript normalization under different experimental conditions. PLoS One. 2012;7(8):e44405.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔC T method. Methods. 2001;25(4):402–8.
We sincerely thank the Genedenovo Biotechnology Co., Ltd. for the technical assistance. We thank Zu-Fu Cheng for the technical assistance. We thank Jennifer Wai (UIUC) for helpful discussion on the data analysis. We thank Irene Lavagi (Plant Bio Press Editing) for editing the manuscript.
Teacher Training Fund in College of Horticulture, FAFU (FAFU2012YYPY05).
Chinese Scholarship Council, China (201407870004).
Availability of data and materials
The raw data had been submitted to the NCBI SRA database, the BioSample accessions number is SRA493837.
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.
Gene Ontology classification of genes. A, CG-vs-ETH; B, CG-vs-1-MCP. (DOCX 558 kb)
KEGG graph of pentose and glucuronate interconversions pathway (only part of the pictures were shown). A, CG-vs-ETH; B, CG-vs-1-MCP. Genes with a red frame are up-regulated DEGs, while down-regulated genes are inside a green frame; genes inside a half red half green frame belong to gene families containing both up- and down-regulated DEGs. The 126.96.36.199 in the frame refers to polygalacturonase (evm.TU.supercontig_250.6, evm.TU.supercontig_92.36). (DOCX 58 kb)
Selected ripening-related genes. (DOCX 23 kb)
KEGG graph of porphyrin and chlorophyll metabolism pathway (CG-vs-ETH). 188.8.131.52 indicates ferretin 1 (evm.TU.supercontig_5.220); 184.108.40.206 indicates hydroxymethylbilane synthase (evm.TU.supercontig_43.43); 220.127.116.11 indicates protoporphyrinogen oxidase (evm.TU.supercontig_130.4); 18.104.22.168 indicates magnesium-chelatase subunit ChlI (evm.TU.supercontig_7.16) and magnesium-chelatase subunit chlH (evm.TU.supercontig_92.51); 22.214.171.124 indicates magnesium-protoporphyrin IX methyltransferase (evm.TU.supercontig_13.164); 126.96.36.199 indicates protochlorophyllide oxidoreductase A (evm.TU.supercontig_80.29); 1.1.1294 indicates NAD (P)-binding rossmann-fold superfamily protein (evm.TU.supercontig_23.100); 188.8.131.52 indicates chlorophyllase 2 (evm.TU.supercontig_44.117); 1.14.1220 indicates pheophorbide a oxygenase, chloroplastic-like (evm.TU.supercontig_47.4). (DOCX 85 kb)
KEGG graph of carotenoid biosynthesis pathway. A, CG-vs-ETH; B, CG-vs-1-MCP. The CruA indicates lycopene beta/epsilon cyclase protein (evm.TU.supercontig_132.5); The CrtR and CrtZ indicate beta-carotene hydroxylase (evm.TU.supercontig_107.106). (DOCX 57 kb)
General pathway of carotenoid metabolism. (DOCX 26 kb)
About this article
Cite this article
Shen, Y.H., Lu, B.G., Feng, L. et al. Isolation of ripening-related genes from ethylene/1-MCP treated papaya through RNA-seq. BMC Genomics 18, 671 (2017). https://doi.org/10.1186/s12864-017-4072-0
- Carica papaya L.