- Research article
- Open Access
Analysis of ripening-related gene expression in papaya using an Arabidopsis-based microarray
BMC Plant Biology volume 12, Article number: 242 (2012)
Papaya (Carica papaya L.) is a commercially important crop that produces climacteric fruits with a soft and sweet pulp that contain a wide range of health promoting phytochemicals. Despite its importance, little is known about transcriptional modifications during papaya fruit ripening and their control. In this study we report the analysis of ripe papaya transcriptome by using a cross-species (XSpecies) microarray technique based on the phylogenetic proximity between papaya and Arabidopsis thaliana.
Papaya transcriptome analyses resulted in the identification of 414 ripening-related genes with some having their expression validated by qPCR. The transcription profile was compared with that from ripening tomato and grape. There were many similarities between papaya and tomato especially with respect to the expression of genes encoding proteins involved in primary metabolism, regulation of transcription, biotic and abiotic stress and cell wall metabolism. XSpecies microarray data indicated that transcription factors (TFs) of the MADS-box, NAC and AP2/ERF gene families were involved in the control of papaya ripening and revealed that cell wall-related gene expression in papaya had similarities to the expression profiles seen in Arabidopsis during hypocotyl development.
The cross-species array experiment identified a ripening-related set of genes in papaya allowing the comparison of transcription control between papaya and other fruit bearing taxa during the ripening process.
Papaya (Carica papaya L.) is an important crop cultivated in tropical and subtropical areas and the ripe fruit has a soft and sweet pulp with high amounts of pro-vitamin A and antioxidants . Papaya is a typical climacteric fruit, with striking colour changes, a rapid rise in ethylene production, and substantial pulp softening; it also responds to exogenous ethylene and 1-MCP applications [2, 3]. The physico-chemical changes during papaya ripening are dependent on the expression of specific genes, and the identification of ripening-related genes involved in the activation of biochemical steps relevant for fruit quality is of both scientific and commercial interest.
In order to understand the network of ripening genes in fleshy fruits, transcriptome studies are valuable tools. In the case of fruit such as tomato, microarrays have been used extensively [4, 5]. However, for less well studied fruits, transcriptome analyses are based on ‘home-made’ microarrays, such as the μPEACH1.0 array , or classical transcript profiling by Differential Display-PCR or cDNA-AFLP [7–9]. With the development of high-throughput sequencing, several species have had their genome sequenced including the Hawaiian variety of papaya fruit . Commercial oligo-chips are not currently available for these organisms and comprehensive RNA sequencing can still be costly often prohibiting routine experiments. However, a cross-species (XSpecies) microarray is an alternative approach that has been successfully used to study the transcriptomes of non-model organisms [11, 12].
Papaya is a member of family Caricaceae within the Brassicales, the same order as the ‘model plant’ Arabidopsis thaliana, which has been the object of many microarray experiments based on commercially available oligo-chips. Because the two species are relatively closely related, the use of Arabidopsis arrays to hybridize RNA from papaya should provide information on the transcriptome changes during ripening in papaya fruits.
In the present study we report the use of RNA from unripe and ripe papaya to probe the Affymetrix Arabidopsis GeneChip ATH1-121501 to profile ripening–related gene expression in papaya. The expression pattern of a number of genes likely to be related to fruit quality was validated by quantitative real-time PCR, and the data from papaya cross-species microarray was compared to microarray data from tomato (a climacteric fruit) and grape (a non-climacteric fruit). A comparative biology approach was then used to compare the putative proteins from papaya and protein sequences from Arabidopsis and other fleshy fruits in order to obtain information on the differences between these organisms in respect to evolutionary role in fruit ripening. The expression of transcription factors was divergent amongst three species, and XSpecies data indicated transcription factors (TFs) that may be involved in the control of papaya ripening. Comparison of the expression patterns of ripening-related TFs and down-stream effectors such as cell wall genes between papaya, fleshy fruits and Arabidopsis indicated both common and unique features in these higher regulatory networks governing ripening.
Transcriptome characteristics of papaya fruit ripening: Use of XSpecies microarray
A probe-masking strategy utilising hybridization of papaya genomic DNA was used to identify probes with low or non-specific hybridization. The number of probe-pairs retained for analysis decreased rapidly (Additional file 1) as the DNA hybridization threshold was increased. In comparison the number of probe-sets reduced at a slower rate, which was consistent with the results obtained with other species [11, 12]. The number of differentially expressed genes (fold change >1.25, p<0.05) was calculated at each threshold and the mask value of 75 returned the highest number of differentially expressed putative genes (414 probe-sets) (Additional file 2). The hierarchical clustering of the log2 values of these probe-sets intensities resulted in the discrimination of eight main clusters (Additional file 3) with different expression patterns. Clusters II, III, IV, V and VI, with 208 probe-sets, were up-regulated, while 205 probe-sets from clusters I, VII and VIII were down-regulated during ripening. Clusters I, II and III was composed by genes with high levels of expression while clusters VI and VIII enclosed genes with the lowest levels of expression during papaya ripening. Clusters IV, V and VII enclosed genes with the highest differences in gene expression when log2 intensities were compared between unripe and ripe fruit.
Papaya probe-sets were separated into gene categories using the PageMan software and up-regulated probe-sets were organized according to the correspondent Arabidopsis Gene Ontology (Additional file 4) function (Figure 1). As was expected in a climacteric fruit, there were a wide range of genes up-regulated during fruit ripening, including those involved in primary metabolism (especially carbohydrate degradation) and energy transport and these findings are consistent with the high demand associated with ripening and the peak in CO2 production (as it was previously reported [3, 8]) (Additional file 5). Genes encoding proteins involved in lipids, protein, and hormone metabolism, as well as cell signalling, signal transduction, and the response to biotic and abiotic stress were also up-regulated (Table 1).
Comparative genomics of fruit ripening in papaya vs. Tomato and grape
In order to visualize the probable cell functions altered in papaya and to compare them to those of another climacteric fruit (tomato) and a non-climacteric fruit (grape), the bins of differentially expressed probe-sets for papaya (fold change >1.25, p<0.10; Additional file 6), tomato and grape (fold change >2.0, authors’ statistical cut-off) were run using the MapMan software . The aim was identify sets of expression profiles that might be conserved between or unique to climacteric and non-climacteric fruit. Data chosen for analysis focused on the main changes associated with the ripening process, such as the comparison between red and mature green tomato and post-véraison and pre-vérasion grape.
Figure 2 shows a schematic Venn diagram of up-regulated probe-sets from papaya XSpecies, tomato and grape microarrays obtained from MapMan analysis (Additional file 7). Despite differences between climacteric and non-climacteric ripening, the proportion of up-regulated genes varied amongst species, with a similar number of over-expressed genes between ripe papaya and grapes.
A schematic overview of ripening regulation and some of its cellular responses (Figure 3) showed a high proportion of up-regulated transcription factors (TFs) in papaya (61%) (Additional file 8). When probe-sets from papaya experiments were analysed in terms of GO molecular function (Additional file 4), 103 probe-sets (25%) accounted for genes responsible for nucleotides and nucleic acids binding, translation initiation factors and ribosomes constituents. Table 2 shows the main up-regulated TFs in the three species regarding the three different ripening experiments.
The cellular response to ripening in these three fruits was also different, especially to biotic and abiotic stress and redox regulation (Figure 3). While tomato responses were the increase of genes related to abiotic stress and decrease of those related to biotic stress, grape responded in an opposite way (Table 3). Nonetheless, papaya responded to ripening in a similar way to tomato, with abiotic genes being expressed in the same time course than tomato ones.
Validation of gene expression and comparative biology between papaya and fleshy fruit organs
In order to validate the XSpecies experiment, the expression of 21 probe-sets that were satisfactory aligned to Arabidopsis genes were analysed by qPCR throughout papaya ripening. Papaya genes were named according to the Arabidopsis functions and grouped together, as shown in Figure 4. Quantitative analysis revealed A-GAL, EXP, FRUCT, LIP3, ACX, ARP, ENOL, HSP70 and ERF2 were induced during ripening, while ARF18, CELL, MEV, LIP, PME1, PME2, PL, RAP2, XET, ERD3, ERD4 and AXR were decreased during ripening. Contrasting results between qPCR and XSpecies experiments were observed for three genes (ARF18, EXP and PME1).
Comparative analysis of the translated amino acids sequences for all genes analysed from qPCR experiments indicated some papaya putative proteins were more similar to those from species with fleshy fruits than to the Arabidopsis ones (Figure 5). In general, proteins from primary and secondary metabolism (FRUCT, LIP, MEV) and transcription factors (ERF2, ARF18, AXR1, ARP) appeared to be phylogenetically closer to those from species with fleshy fruits, whereas abiotic stress (ERD3, ERD4) and cell wall metabolism (A GAL, PL, PME1, PME2, CELL, XTH) proteins are more closely related to those from Arabidopsis. Contrasting results were observed for proteins with high degree of conservation, such as ENOL, RAP2, HSP, ACX and PL, where homology was determined primarily by non-conserved regions and less by conserved domains. Because cell wall related genes shared greater homology with Arabidopsis, a comparison of gene expression in ripe papayas was conducted using a published Arabidopsis microarray study of gene expression in 5 or 11-day-old hypocotyls (ratio 11-day/5-day) . These data (Table 4) also indicate that this set of genes had comparable gene expression levels and might act on cell wall polymers in a similar way.
In this study we investigated the transcriptome of ripening papaya [7, 8] using fruits at two contrasting physiological stages and the Affymetrix ATH1-121501 GeneChip for Arabidopsis in an XSpecies microarray . The number of differentially expressed probe-sets was consistent with the results obtained with other species [11, 12] but was fewer than those from regular microarray experiments, especially with regarding comparative analyses between unripe and ripe fruits [15, 16]. Although XSpecies results have less power than a regular microarray experiment and technologies such as RNAseq performed by Next Generation Sequencing (NGS), the technique was able to identify important genes related to metabolic processes involved in papaya ripening. Surprisingly, despite papaya and Arabidopsis being members of the order Brassicales, only cell wall-related genes shared a greater homology among the genes studied. This may reflect shared metabolic pathways which differ from other fleshy fruit. Fleshy fruits showed greater homology between transcription factors which might indicate shared transcriptional regulatory networks.
Deduced roles of ripening-related genes in papaya
The high energy demand at climacteric ripening was evidenced by the up-regulation of genes related to Krebs and TCA cycles (AT5G08300, AT2G29990, AT3G24200 and AT5G51060), cytochromes (AT2G05180, AT5G06900, AT1G01280), and hexose metabolism (AT2G36530 and AT1G50390). Similarly, genes involved in lipid metabolism, were also affected. These included cell membrane lipases (AT3G62590, AT1G53920 and AT1G28590), and up-regulation of genes responsible for the synthesis of unsaturated fatty acids (AT4G30950 and AT3G15850). Down-regulation of the MEV gene (AT2G38700) suggests the main precursors for volatiles production and carotenoid biosynthesis could be produced by the MEP pathway inside the plastids during papaya ripening . Fabi et al.  identified another transcript related to MEV gene, suggesting mevalonate accumulation prior ripening.
When the probe-sets were analysed in terms of GO molecular function (Additional file 4), a high number of genes responsible for nucleotide and nucleic acid binding, translation initiation and ribosomes constituents were observed, indicating control of ripening at both the level of transcription and translation. Five members of the MYB family, which is important for ripening since MYB genes are known to regulate secondary metabolism and colour accumulation , were differentially expressed during ripening. Since colour of papaya pulp is correlated with accumulation of carotenoids , papaya MYB TFs might show distinct role in other metabolism during ripening.
Three closely-related ERF/AP2 genes had diverse expression throughout papaya ripening, with an AP2 gene (AT4G39780) and an ERF gene (AT1G15360) being up-regulated (Table 2). While some genes from AP2/ERF family are related to improving Arabidopsis resistance to abiotic stress , others have been shown to be responsible for controlling tomato ripening [20, 21]. Specific APETALA2 genes (so called SlAP2a to SlAP2e genes) are expressed throughout normal tomato ripening and SlAP2a gene may balance the activities of positive ripening regulators as a negative feedback loop. Moreover, a previous study showed an APETALA2 gene subclade IIIc (COLD BIND FACTOR II gene) is involved in regulation of apple pulp softening during cold storage and/or ethylene treatment . In fact, papaya ERF2 and RAP2 putative proteins share higher similarity to diverse APETALA2 proteins from apples and tomato (Figure 6) rather than the SlAP2a protein from tomato.
Three genes members of auxin signalling pathway (ARF18, AXR1 and ARP), a hormone that impairs ripening , were identified in our experiments. Whereas ARP was the only up regulated gene during ripening, ARF18 was phylogenetically close to SlARF family, the auxin response factors from S. lycopersicum (Figure 5). These data suggest that ripening in papaya involves auxin signalling, in common with other fruits [23, 24].
Many different TFs are involved in control of ripening in fleshy fruits. In tomato, master regulators include MADS-box (SEP4 like, RIN, TDR4, TAG1, TAGL1), SBP-box (CNR), HB box (LeHB 1) and NAC genes . Strawberry ripening involves a SEP1/2 like gene (FaMADS9) , and for banana a SEP3 like gene (MaMADS2) . In bilberry, accumulation of anthocyanins is controlled by an SQUAMOSA class MADS-box TF VmTDR4, orthologous to the TDR4 gene in tomato . Papaya shows range of TFs which are homologous to these master regulators and show ripening-related changes in gene expression (Table 2). Despite each species seeming to have specific sets of ripening-related transcription factors, their exact functions in papaya fruits remain to be elucidated.
Cell wall-related genes
Papaya pulp softening shows a remarkable change during ripening, and this is thought to be due to the activities of cell wall hydrolases . Seven putative cell wall-related genes were differentially expressed: four related to pectin hydrolysis (α-galactosidase - AT3G26380; pectate lyase - AT4G24780; pectin methylesterases - AT5G48450 and AT4G12420), two related to cellulose hydrolysis and rearrangement (cellulase - AT4G28320; xyloglucan endotransglycosylase - AT2G06850) and one related to cell expansion (expansin - AT2G39700). Mostly genes related to PGs and expansin were up-regulated while those related to PMEs, PLs and glucosidases were down-regulated.
The up-regulation of an α-galactosidase gene (A GAL) observed in our experiments is consistent with previous data reported by Soh et al.  and Nogueira et al.  on ripening of papaya fruit. On the other hand, changes in pectate lyases (PL) were detected, but the reduced levels contrast with the expected role of the enzyme during ripening. In other fruits, such as bananas  and mangoes PL is thought to contribute to the pectin disassembly that leads to the pulp softening of ripe fruits. Pectin methylesterases were also apparently down-regulated during ripening. Solubilisation of pectin in ripe papayas may be governed principally by the action of polygalacturonases on previously de-esterified pectin chains [29, 34] in addition to the action of α-galactosidases. Cellulose and hemicellulose processing in the papaya pulp could not be inferred from our results, since the genes related to hydrolysis and rearrangements of cellulose (cellulase and xyloglucan endotransglycosylase, respectively) were down-regulated during papaya ripening, in contrast to the up-regulation in ripe tomato and strawberry [35, 36]; respectively). Gaete-Eastman et al.  observed expansin gene expression was up-regulated during ripening in mountain papaya and it was inversely correlated to pulp texture. This pattern of expression is in agreement with the expansin gene from ‘Golden’ papaya (Figure 7). Papaya EXP is related to the gene products encoded by FaEXP4, PcEXP5 and PpEXP2 from strawberry, pear and peach, respectively [38–40]; Figure 6). All of these genes are up-regulated during fruit ripening, suggesting conservation of gene action in the softening of a range of fleshy fruits species. However in some instances cell wall changes in papaya resembled events in non-fruit tissues. An Arabidopsis microarray study  investigating gene expression during 5 or 11-day-old hypocotyls (ratio 11-day/5-day) showed a transcription profile, at least for some of the cell wall-related genes, that was highly similar to that found in ripening papaya (Table 4). This reveals interesting characteristics of cell wall remodelling which may have evolved from the same ancestral Brassicales.
Response to stress and plant defence
No representative differences in expression of probe-sets related to ‘response to stress’ (six probe-sets up-regulated and six probe-sets down-regulated) were observed, and only two out of six ‘plant defence response’ probe-sets were up-regulated. However, genes from other GO classes may be considered based on their putative functions. Acyl-CoA oxidases (ACXs) are members of ‘lipid metabolism’ GO, and have a key role in the jasmonic acid (JA) biosynthesis, an important compound in protection against pathogens and insects . The up-regulated papaya ACX gene could contribute to the defence system during ripening, since it showed a 5-fold increase in transcript abundance. Papaya ACX protein shared high similarity with that from tomato, where a peroxisomal ACX protein is critical for the β-oxidation through JA biosynthesis and systemic wound signalling, indicating this gene might also be responsible for JA biosynthesis in papaya . Two putative papaya ERD3 and ERD4 genes (AT2G39750 and AT4G04340, respectively) were down-regulated. These proteins are methyltransferases that respond positively to abiotic stresses such as cold treatment and prolonged mild osmotic stress .
Analysis of the papaya transcriptome reveals that like tomato there is overexpression of heat-shock protein (HSPs) genes. Heat-shock proteins can work as chaperones in protein folding; under stress conditions they can re-establish normal protein conformation  and degrade damaged or misfolded peptides . Genes normally associated with response to abiotic stress are often seen expressed in ripening fruits [8, 46].
Chitinases and thaumatins were over-expressed in grapes, conferring plant protection against saprophytic organisms , but down-regulated in tomato instead , none of these genes were differentially expressed in papaya. Only some genes related to protein inhibition and those of so-called leucine repeat proteins (LRRs), which are very important in defence against pathogenic fungi  were induced. Regarding redox regulation, ascorbic acid (Vitamin C) accumulates in these three fruits [49–51], but only papaya cross-species microarray revealed genes related to ascorbic acid biosynthesis: the up regulated monodehydroascorbate reductase and an L-galactose dehydrogenase genes.
The heterologous hybridization microarray was successfully applied in the study of transcripts changes associated to papaya fruit, a commercially relevant crop from a non-model organism. The ripening of papaya represents a time-course of cellular metabolism changes, characterized by differential expression of numerous genes involved in primary metabolism, hormonal signalling, transcriptional regulation, abiotic stress and cell wall metabolism. Among the genes identified are transcription factors (TFs) of the MADS-box, NAC and AP2/ERF gene families, which are master regulators in other fruits indicating conservation of function for ripening control genes across different taxa. Moreover, data revealed cell wall-related gene expression was more similar to Arabidopsis hypocotyl development profiles than those from other fleshy fruit perhaps revealing characteristics of cell wall remodelling mechanisms specific to the Brassicales.
The phylogenetic relationship between transcription factors of fleshy fruits might indicate shared transcriptional regulatory networks. Although data presented in the manuscript is not enough to indicate how climacteric and non-climacteric fruits evolved, and some genes specific to papaya or those specifically expressed in fleshy organs might have been overlooked with the heterologous microarray, other approaches, e.g. RNAseq, are likely to contribute to new perspectives on how different fleshy fruits respond to ripening.
Papaya fruit (Carica papaya L. cv. ‘Golden’) were harvested at colour break to ¼ yellow (around 150 days after anthesis). One replicate of each unripe and ripe fruit samples for microarray assay was obtained from a previous study in 2007 . The other two replicates of unripe and ripe samples for microarray and also for quantitative PCR analyses were obtained in 2010. Soon after harvest, the initial respiration and ethylene levels were determined . Each experimental sample (unripe papaya and successive time points during the ripening process) comprised 12 individual fruits. After removal of the peel and seeds, the sliced pulp of each fruit was frozen in liquid N2 and stored at −80°C.
Genomic DNA and total RNA extractions
Genomic DNA and total RNA extraction protocols were the same as previously described . Total RNA was purified using the “RNeasy® Kit” (Qiagen). Nucleic acids were quantified in the NanoDrop® ND-1000 (Nanodrop Technologies©) and gel analyzed to verify their integrity . Total RNA were also analyzed with Agilent 2100 Bioanalyser” (Agilent Technologies©) to confirm its integrity before each hybridization.
Genomic DNA hybridization, probe-selection and cRNA hybridization
Genomic DNA was labelled with the Bioprime DNA labelling System kit (Invitrogen™), hybridized with Affymetrix GeneChip ATH1-121501 (Affymetrix) for 16 hours at 45°C using standard hybridization protocols (Affymetrix©) and analyzed by computed scanning. Using a perl script  and a range of user-defined threshold values (from 0 to 500), chip definition files (CDF) were created with papaya genomic DNA hybridization data (mask files). Probe-pairs were retained for analysis if their signal values were greater than the defined threshold. Probe-sets were retained if they contained more than one probe-pair. Total RNA (5 μg) from papaya pulp was reverse-transcribed to generate first strand cDNA containing a 5’-T7 RNA polymerase promoter sequence. Double stranded cDNA was synthesized using standard protocol, and the resulting samples were in vitro transcribed by T7 DNA polymerase using biotinylated nucleotides to generate complementary RNAs (cRNAs). Purified cRNA (15 μg) were heat-fragmented and hybridized to ATH1-121501 for 16 hours at 45°C. The complete protocol has already been published . For each fruit stage (unripe and ripe), a triplicate of hybridization was undertaken and all the hybridizations have been submitted to GEO (http://www.ncbi.nlm.nih.gov/geo; accession number GSE38105).
Microarray Analysis Suite (MAS Version 5.0; Affymetrix) was used to generate .CEL files for each RNA hybridization. These files were loaded into GeneSpring version 7.2 (Agilent Technologies) software using the Robust Multichip Average (RMA) pre-normalization algorithm . The computer files that were being loaded in GeneSpring were filtered using the CDF files generated using the genomic DNA hybridization. For each replicate array, each probe-set signal value from ripe samples was compared to the probe-set signal value of unripe samples to give gene expression ratios. Differentially expressed genes were identified using one-way ANOVA with a Benjamini and Hochberg false discovery rate multiple testing correction. The differentially expressed genes had their log2 signal intensities computed from the replicate chips and hierarchical clustering was carried out using the EPCLUST software  with a complete linkage algorithm and the Euclidian distance on normalized vectors of length 1 (chord distance) as parameters.
Comparative analyses using PageMan and MapMan
Normalized gene expression data were subjected to analysis of functional categories using PageMan and MapMan functional categories [56, 57]. Using the Wilcoxon test it was possible to assume whether significantly more genes in ripe vs. unripe point were up-regulated when normalized to their average expression. Expression data (log2 fold) from papaya microarray (ripe X unripe) was loaded and analysed in MapMan software for visualization of the cellular pathway  against Arabidopsis mapping (Ath_AFFY_ATH1_TAIR8_Jan2010). In the same way, expression data (log2 fold) from Solanum lycopersicum cultivar MicroTom (ripe X mature green; ) and Vitis vinifera cv. Pinot Noir (post-veráison X pre-véraison; ) were loaded and analysed in MapMan software, using data from Ozaki et al.  against S. lycopersicum mapping (Slyc_AFFY_SGN_BUILD2_070709) and Pilati et al.  against V. vinifera mapping (Vvin_AFFY_09), respectively. Log2 fold values of differentially expressed probe-sets from three species were compared and a Venn diagram enclosing only the log2 values, not the probe-set IDs (since three different platforms were used) was created. For papaya up-regulated probe-sets, the Arabidopsis gene annotations and functional classifications were analysed using the gene ontology (GO) function of the GeneSpring software.
Quantitative analysis of gene expression by real time-PCR (qPCR)
For the validation of putative papaya genes identified by the XSpecies hybridization, the Arabidopsis genes represented in the hybridized chips were aligned individually against the WGS database using the BLASTN tool (score ≥ 100 and e-value ≤ 1e-30 as cut-off values). Papaya putative coding sequences that were satisfactory aligned to Arabidopsis genes were evaluated following the ‘Minimum Information for Publication of Quantitative Real-Time PCR Experiments – MIQE’  and also according to the parameters: (1) pair of primers with melting temperature of 60°C and absence of primers-dimers and hairpins; (2) amplicon size between 75 and 200 bp; (3) amplicon evaluation with the mfold program ; and (4) amplicon alignment with an unique papaya WGS sequence (in order to not amplify two closely-related genes or a duplicated one). Twenty-one putative genes satisfied the above mentioned criteria and had their expression levels quantified. Primers were designed using Primer 3 (v.0.4.0) tool , and the sequences are shown in Additional file 9. As an internal controls, the putative actin gene located on chromosome LG9 contig 1059 (GenBank accession no. ABIM01001059) was used with sense (5′-CGT GAC CTT ACT GAT CAC TTG-3’) and reverse (5’-GTC AAG GGC AAT GTA AGA CAG-3’) primers in combination with the 18S rRNA (GenBank accession no. U42514) with sense (5′-AAA CGG CTA CCA CAT CCA AG-3’) reverse (5’-CGA AGA GCC CGG TAT TGT TAG GG-3’) primers. After on-column digestion of DNA with DNase RNase-free (NucleoSpin® - Macherey Nagel©), first-strand cDNA was synthesized with random primers from 1 μg total RNA using the ImProm-II Reverse Transcription System (Promega). For primer testing and identity confirmation, the fragments from preliminary PCR were cloned and sequenced. Real-time PCR was performed using the Platinum SYBR Green qPCR Supermix UDG (Invitrogen) in a “Rotor-Gene 3000 four channel Multiplexing System” (Corbett Research). The amplification was carried out under the following conditions; 50°C for 2 min followed by an initial denaturation step at 95°C for 2 min, 40 cycles at 95°C for 15 s, 60°C for 30 s, and 72°C for 30 s. Non Template Controls (NTCs) and melting curve analyses of amplicons were analysed for all experiments. The threshold cycle (Ct) values of the four technical replicate reactions were averaged using the Rotor-Gene 3000 software and quantification was performed using the relative standard curve method . Samples used in qPCR experiments were a mixture of, at least, 32 fruits from 2010 second and third biological replicates  (Additional file 5). Developmental parameters included ripening-related events such as CO2 and ethylene production as well as pulp softening . The results of the standard curves calculation are in Additional file 10. Data were analysed against the first day after harvest by one-way ANOVA, and means were compared using the Tukey test at α<0.05. Statistical analysis was carried out using OriginPro version 8 (OriginLab®).
Comparative biology analysis between fleshy fruit organisms
In order to identify genes similarities shared between some fleshy fruit organs, a comparative biology analysis was done using qPCR tested papaya genes. The validated probe-sets were aligned to WGS database of papaya genome and the putative corresponding proteins were identified. Sequential TBlastN analyses were done using the putative papaya proteins and the WGS database for fleshy fruit organisms such as tomato (S. lycopersicum), strawberry (F. vesca), peach (P. persica) and grape (V. vinifera), and the corresponding proteins were also obtained for each organism. If otherwise stated, all sequences were identified from this manner. Together with Arabidopsis proteins and some other proteins, they were aligned using ClustalW  and phylogenetic trees were drawn using the Phylodendron© software . The Genbank ID and sequences of proteins used in comparative biology analysis are in Additional file 11.
Basic Local Alignment Search Tool
Polymerase chain reaction
Whole genome shotgun
- XSpecies Microarray:
Cross Species Microarray.
Fabi JP, Peroni FHG, Gomez MLPA: Papaya, mango and guava fruit metabolism during ripening: postharvest changes affecting tropical fruit nutritional content and quality. Fresh Produce. 2010, 4: 56-66.
Manenoi A, Bayogan ERV, Thumdee S, Paull RE: Utility of 1-methylcyclopropene as a papaya postharvest treatment. Post Biol Technol. 2007, 44: 55-62. 10.1016/j.postharvbio.2006.11.005.
Fabi JP, Cordenunsi BR, Barreto GPM, Mercadante AZ, Lajolo FM, Nascimento JRO: Papaya fruit ripening: response to ethylene and 1-methylcyclopropene (1-MCP). J Agric Food Chem. 2007, 55: 6118-6123. 10.1021/jf070903c.
Moore S, Payton P, Wright M, Tanksley S, Giovannoni JJ: Utilization of tomato microarrays for comparative gene expression analysis in the Solanaceae. J Exp Botany. 2005, 56: 2885-2895. 10.1093/jxb/eri283.
Jiang F, Zheng X, Chen J: Microarray analysis of gene expression profile induced by the biocontrol yeast Cryptococcus laurentii in cherry tomato fruit. Gene. 2009, 430: 12-16. 10.1016/j.gene.2008.09.036.
Trainotti L, Tadiello A, Casadoro G: The involvement of auxin in the ripening of climacteric fruits comes of age: the hormone plays a role of its own and has an intense interplay with ethylene in ripening peaches. J Exp Bot. 2007, 58: 3299-3308. 10.1093/jxb/erm178.
Fabi JP, Lajolo FM, Nascimento JRO: Cloning and characterization of transcripts differentially expressed in the pulp of ripening papaya. Sci Hort. 2009, 121: 159-165. 10.1016/j.scienta.2009.01.036.
Fabi JP, Mendes LRBC, Lajolo FM, Nascimento JRO: Transcript profiling of papaya fruit reveals differentially expressed genes associated with fruit ripening. Plant Sci. 2010, 179: 225-233. 10.1016/j.plantsci.2010.05.007.
Cin VD, Barbaro E, Danesin M, Velasco HMR, Ramina A: Fruitlet abscission: a cDNA-AFLP approach to study genes differentially expressed during shedding of immature fruits reveals the involvement of a putative auxin hydrogen symporter in apple (Malus domestica L. Borkh). Gene. 2009, 442: 26-36. 10.1016/j.gene.2009.04.009.
Ming R, Hou S, Feng Y, et al: The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus). Nature. 2008, 452: 991-997. 10.1038/nature06856.
Hammond JP, Broadley MR, Craigon DJ, Higgins J, Emmerson Z, Townsend H, White PJ, May ST: Using genomic DNA-basedprobe-selection to improve the sensitivity of high-densityoligonucleotide arrays when applied to heterologous species. Plant Methods. 2005, 9 (1): 10.
Bagnaresi P, Moschella A, Beretta O, Vitulli F, Ranalli P, Perata P: Heterologous microarray experiments allow the identification of the early events associated with potato tuber cold sweetening. BMC Genomics. 2008, 9: 176-10.1186/1471-2164-9-176.
Thimm O, Blasing O, Gibon Y, Nagel A, Meyer S, Kruger P, Selbig J, Muller LA, Rhee SY, Stitt M: MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004, 37: 914-939. 10.1111/j.1365-313X.2004.02016.x.
Jamet E, Roujol D, San-Clement H, Irshad M, Soubigou-Taconnat L, Renou JP, Pont-Lezica R: Cell wall biogenesis of Arabidopsis thaliana elongating cells: transcriptomics complements proteomics. BMC Genomics. 2009, 10: 505-10.1186/1471-2164-10-505.
Pilati S, Perazzolli M, Malossini A, Cestaro A, Demattè L, Fontana P, Dal Ri A, Viola R, Velasco R, Moser C: Genome-wide transcriptional analysis of grapevine berry ripening reveals a set of genes similarly modulated during three seasons and the occurrence of an oxidative burst at vèraison. BMC Genomics. 2007, 8: 428-10.1186/1471-2164-8-428.
Ozaki S, Ogata Y, Suda K, Kurabayashi A, Suzuki T, Yamamoto N, Iijima Y, Tsugane T, Fujii T, Konishi C, Inai S, Bunsupa S, Yamazaki M, Shibata D, Aoki K: Coexpression analysis of tomato genes and experimental verification of coordinated expression of genes found in a functionally enriched coexpression module. DNA Res. 2010, 17: 105-116. 10.1093/dnares/dsq002.
Rodríguez-Concepción M: Supply of precursors for carotenoid biosynthesis in plants. Arch Bioch Biophysics. 2010, 504: 118-122. 10.1016/j.abb.2010.06.016.
Allan AC, Roger PH, Laing WA: MYB transcription factors that colour our fruit. T Plant Sci. 2008, 13: 99-102. 10.1016/j.tplants.2007.11.012.
Sharabi-Schwager M, Amnon Lers , Alon Samach , Charles L, Porat G: Overexpression of the CBF2 transcriptional activator in Arabidopsis delays leaf senescence and extends plant longevity. J Exp Botany. 2010, 61: 261-273. 10.1093/jxb/erp300.
Chung MY, Vrebalov V, Alba R, Lee J, McQuinn R, Chung JD, Klein P, Giovannoni J: A tomato (Solanum lycopersicum) APETALA2/ERF gene, SlAP2a, is a negative regulator of fruit ripening. Plant J. 2010, 64: 936-947. 10.1111/j.1365-313X.2010.04384.x.
Karlova R, Rosin FM, Busscher-Lange J, Parapunova V, Do PT, Fernie AR, Fraser PD, Baxter C, Angenent GC, de Maag RA: Transcriptome and metabolite profiling show that APETALA2a is a major regulator of tomato fruit ripening. Plant Cell. 2011, 23: 923-941. 10.1105/tpc.110.081273.
Taken E, Ireland H, Gunaseelan K, Karunairetnam S, Wang D, Schultz K, Bowen J, Atkinson RG, Johnston JW, Putterill J, Hellens RP, Schaffer RJ: The role of ethylene and cold temperature in the regulation of the apple POLYGALACTURONASE1 gene and fruit softening. Plant Physiol. 2010, 153: 294-305. 10.1104/pp.109.151092.
Li J, Dai X, Zhao Y: A role for auxin response factor 19 in auxin and ethylene signaling in Arabidopsis. Plant Physiol. 2006, 140: 899-908. 10.1104/pp.105.070987.
De Jong M, Wolters-Arts M, Feron R, Mariani C, Vriezen WH: The Solanum lycopersicum auxin response factor 7 (SlARF7) regulates auxin signaling during tomato fruit set and development. Plant J. 2009, 57: 160-170. 10.1111/j.1365-313X.2008.03671.x.
Klee HJ, Giovannoni JJ: Genetics and control of tomato fruit ripening and quality attributes. Annu Rev Genetics. 2011, 45: 41-59. 10.1146/annurev-genet-110410-132507.
Seymour GB, Ryder CD, Cevik V, Hammond JP, Popovich A, King GJ, Vrebalov J, Giovannoni JJ, Manning K: A SEPALLATA gene is involved in the development and ripening of strawberry (Fragaria x ananassa duch.) fruit, a non-climacteric tissue. J Exp Botany. 2011, 62: 1179-1188. 10.1093/jxb/erq360.
Elitzur T, Vrebalov J, Giovannoni JJ, Goldschmidt EE, Friedman H: The regulation of MADS-box gene expression during ripening of banana and their regulatory interaction with ethylene. J Exp Botany. 2010, 61: 1523-1535. 10.1093/jxb/erq017.
Jaakola L, Poole M, Jones MO, Kämäräinen-Karppinen T, Koskimäki JJ, Hohtola A, Häggman H, Fraser PD, Manning K, King GJ, Thomson H, Seymour GB: A SQUAMOSA MADS box gene involved in the regulation of anthocyanin accumulation in bilberry fruits. Plant Physiol. 2010, 153: 1619-1629. 10.1104/pp.110.158279.
Fabi JP, Cordenunsi BR, Seymour GB, Lajolo FM, Nascimento JRO: Molecular cloning and characterization of a ripening-induced polygalacturonase related to papaya fruit softening. Plant Physiol Biochem. 2009, 47: 1075-1081. 10.1016/j.plaphy.2009.08.002.
Soh CP, Ali ZM, Lazan H: Characterisation of a α-galactosidase with potential relevance to ripening related texture changes. Phytochem. 2006, 67: 242-254. 10.1016/j.phytochem.2005.09.032.
Nogueira SB, Labate CA, Gozzo FC, Pilau EJ, Lajolo FM, Nascimento JRO: Proteomic analysis of papaya fruit ripening using 2DE-DIGE. J Proteomics. 2012, 75: 1428-1439. 10.1016/j.jprot.2011.11.015.
Marín-Rodríguez MC, Smith DL, Manning K, Orchard J, Seymour GB: Pectate lyase gene expression and enzyme activity in ripening banana fruit. Plant Mol Biol. 2003, 51: 851-857. 10.1023/A:1023057202847.
Chourasia A, Sane VA, Nath P: Differential expression of pectate lyase during ethylene-induced postharvest softening of mango (Mangifera indica var. Dashehari). Physiol Plant. 2006, 128: 546-555. 10.1111/j.1399-3054.2006.00752.x.
Shiga TM, Fabi JP, Nascimento JRO, Petkowicz CLO, Vriesmann LC, Lajolo FM, Cordenunsi BR: Changes in cell wall composition associated to the softening of ripening papaya: evidence of extensive solubilisation of large molecular mass galactouronides. J Agric Food Chem. 2009, 57: 7064-7071. 10.1021/jf900073b.
Schroder R, Wegrzyn TF, Sharma NN, Atkinson RG: LeMAN4 Endo-beta-mannanase from ripe tomato fruit can act as a mannan transglycosylase or hydrolase. Planta. 2006, 224: 1091-1102. 10.1007/s00425-006-0286-0.
Opazo MC, Figueroa CR, Henríquez J, Herrera R, Bruno C, Valenzuela PD, Moya-León MA: Characterization of two divergent cDNAs encoding xyloglucan endotransglycosylase/hydrolase (XTH) expressed in fragaria chiloensis fruit. Plant Sci. 2010, 179: 479-488. 10.1016/j.plantsci.2010.07.018.
Gaete-Eastman C, Figueroa CR, Balbontína C, Moya M, Atkinson RG, Herrera R, Moya-León MA: Expression of an ethylene-related expansin gene during softening of mountain papaya fruit (Vasconcellea pubescens). Post Biol Technol. 2009, 53: 58-65. 10.1016/j.postharvbio.2009.03.007.
Dotto MC, Matrínez GA, Civello PM: Expression of expansin genes in strawberry varieties with contrasting fruit firmness. Plant Physiol Biochem. 2006, 44: 301-307. 10.1016/j.plaphy.2006.06.008.
Hisawa K, Rose JKC, Nakano R, Inaba A, Kubo Y: Differential expression of seven α-expansin genes during growth and ripening of pear fruit. Physiol Plant. 2003, 117: 564-572. 10.1034/j.1399-3054.2003.00064.x.
Hayama H, Moriguchi AI, Kashimura Y: Identification of a new expansin gene closely associated with peach fruit softening. Post Biol Technol. 2003, 29: 1/10
Kang JH, Wang L, Giri A, Baldwin IT: Silencing threonine deaminase and JAR4 in Nicotiana attenuata impairs jasmonic acid-isoleucine-mediated defenses against Manduca sexta. Plant Cell. 2006, 18: 3303-3320. 10.1105/tpc.106.041103.
Li C, Schilmiller AL, Liu G, Lee GI, Jayanty S, Sageman C, Vrebalov J, Giovannoni JJ, Yagi K, Kobayashi Y, Howe GA: Role of beta-oxidation in jasmonate biosynthesis and systemic wound signaling in tomato. Plant Cell. 2005, 17: 971-986. 10.1105/tpc.104.029108.
Skirycz A, De Bodt S, Obata T, De Clercq I, Claeys H, De Rycke R, Andriankaja M, Van Aken O, Van Breusegem F, Fernie AR, Inze D: Developmental stage specificity and the role of mitochondrial metabolism in the response of Arabidopsis leaves to prolonged mild osmotic stress. Plant Physiol. 2010, 152: 226-244. 10.1104/pp.109.148965.
Wang W, Vinocur B, Shoseyov O, Altman A: Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci. 2004, 9: 244-252. 10.1016/j.tplants.2004.03.006.
Faurobert M, Mihr C, Bertin N, Pawlowski T, Negroni L, Sommerer N, Causse M: Major proteome variations associated with cherry tomato pericarp development and ripening. Plant Physiol. 2007, 143: 1327-1346. 10.1104/pp.106.092817.
Manganaris GA, Rasori A, Bassi D, Geuna F, Ramina A, Tonutti P, Bonghi C: Comparative transcript profiling of apricot (Prunus armeniaca L.) fruit development and on-tree ripening. Tree Genetics Genomes. 2011, 7: 609-616. 10.1007/s11295-010-0360-4.
Rotter A, Camps C, Lohse M, Kappel C, Pilati S, Hren M, Stitt M, Coutos-Thévenot P, Moser C, Usadel B, Delrot S, Gruden K: Gene expression profiling in susceptible interaction of grapevine with its fungal pathogen eutypa lata: extending MapMan ontology for grapevine. BMC Plant Biol. 2009, 9: 104-10.1186/1471-2229-9-104.
De Lorenzo G, D’Ovidio R, Cervone F: The role of polygalacturonase-inhibiting proteins (pgips) in defense against pathogenic fungi. Annu Rev Phytopathol. 2001, 39: 313-335. 10.1146/annurev.phyto.39.1.313.
Ioannidi E, Kalamaki MS, Engineer C, Pateraki I, Alexandrou D, Mellidou I, Giovannonni JJ, Kanellis AK: Expression profiling of ascorbic acid-related genes during tomato fruit development and ripening and in response to stress conditions. J Exp Botany. 2009, 60: 663-678. 10.1093/jxb/ern322.
Cruz-Rus E, Botella MA, Valpuesta V, Gomez-Jimenez MC: Analysis of genes involved in L-ascorbic acid biosynthesis during growth and ripening of grapeberries. J Plant Physiol. 2010, 167: 739-748. 10.1016/j.jplph.2009.12.017.
Ali A, Muhammad MTM, Sijam K, Siddiqui Y: Effect of chitosan coatings on the physicochemical characteristics of eksotika II papaya (Carica papaya L.) fruit during cold storage. Food Chem. 2011, 124: 620-626. 10.1016/j.foodchem.2010.06.085.
Sambrook J, Fritsch EF, Maniatis T: Molecular Cloning: A Laboratory Manual. 2000, Cold Spring Harbor, New York USA: Cold Spring Harbor Laboratory, 3
The perl script was downloaded from The Nottingham Arabidopsis Stock Centre (NASC) website (http://arabidopsis.info/), specifically at the XSpecies section. (http://affymetrix.arabidopsis.info/xspecies/).
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4: 249-264. 10.1093/biostatistics/4.2.249.
The software used for hierarchical clustering was accessed on line from Expression Profiler website. Analysis and clustering of gene expression and sequence data - available at http://www.bioinf.ebc.ee/EP/EP.
Usadel B, Nagel A, Steinhauser D, Gibon Y, Blaesing OE, Redestig H, Sreenivasulu N, Krall L, Hannah MA, Poree F, Fernie AR, Stitt M: PageMan an interactive ontology tool to generate, display, and annotate overview graphs for profiling experiments. BMC Bioinforma. 2006, 7: 535-10.1186/1471-2105-7-535.
Usadel B, Nagel A, Thimm O, Redestig H, Blaesing OE, Palacios-Rojas N, Selbig J, Hannemann J, Piques MC, Steinhauser D, Scheible WR, Gibon Y, Morcuende R, Weicht D, Meyer S, Stitt M: Extension of the visualization tool MapMan to allow statistical analysis of arrays, display of corresponding genes, and comparison with known responses. Plant Physiol. 2005, 138: 1195-1204. 10.1104/pp.105.060459.
Expression data from Solanum lycopersicum cultivar MicroTom was downloaded from Tomato Functional Genomics Database website available at http://ted.bti.cornell.edu.
Expression data from Vitis vinifera cultivar Pinot Noir was downloaded from BMC website available at http://www.biomedcentral.com/1471-2164/8/428, referred to additional file 3 from reference 15 cited above (http://www.biomedcentral.com/content/supplementary/1471-2164-8-428-s3.xls).
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT: The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clinical Chem. 2009, 4: 611-622.
Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595. Program available at http://mfold.bioinfo.rpi.edu/cgi-bin/dna-form1.cgi
Primers for Real Time-PCR experiments were designed using Primer 3 (v.0.4.0) tool available at http://fokker.wi.mit.edu/primer3/input.htm
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: 2002-2007.
Thompson JD, Higgis DG, Gibson TJ: CLUSTALW: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Phylogenetic trees for comparative biology analysis were drawn using the Phylodendron© software available at http://iubio.bio.indiana.edu/treeapp
The authors thank the NASC team for helping with XSpecies hybridizations, Dr. Mervin Poole for assistance with the qPCR analyses and Caliman Agrícola S/A for providing fruits. The Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) supported this work through a grant (Project 02/12452-9) and scholarship to J.P.F. (03/00932-9 and 07/56515-8).
The authors declare that they have no competing interests.
JPF wrote the article and performed the research. GBS and JRON corrected the manuscript and helped in its preparation. NSG, MRB and STM helped in data analysis and manuscript preparation. FML and BRC helped in manuscript preparation. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Carica papaya genomic DNA (gDNA) hybridization on probe- pairs and probe- sets from the ATH1- 121501 chip. This additional figure describes how different values of hybridization intensity threshold (masks) affect the probe-sets and probe-pairs retained after papaya genomic hybridization in A.thaliana commercial chip. (PDF 134 KB)
Additional file 2: -sets identified in papaya fruit. This table is a full version (edited in Excel program) of Table 1 with crude data (19,286 probe-sets) and treated data (414 probe-sets). It describes the principal characteristics of probe-sets retained after threshold cut-off and comparison between ripe X unripe fruit, that includes systematic and common names of probe sets, normalized and raw values, fold time changes and p values, as well as GO identities for biological process, molecular function and cellular component. (XLS 3 MB)
Additional file 3: 121501 probe- sets hybridized with papaya RNA. This figure describes the hierarchical clustering of the 414 probe-sets identified by the XSpecies microarray technique, showing different clusters of gene expression. (DOCX 664 KB)
Additional file 4: sets identified in papaya fruit. This table is a shorter version of Additional file 2, also edited in Excel program. However this table shows the probe-sets collected for generation of GO plot figure by names (Figure 1A and B) and also describes the parameters of Genespring software from statistical analyses. (XLSX 167 KB)
Additional file 5: This table describes the main parameters that were used for classifying as unripe (green) and ripe (yellow) papayas. This includes respiration (CO2 production), ethylene production and pulp firmness. (PDF 90 KB)
Additional file 6: sets identified in papaya fruit. This table shows, in a reduced form, the probe-sets that were differentially expressed according to fold changes between ripe X unripe, but with a statistical p-value cut-off of 0.10, which returned 1,091 differentially expressed probe-sets. (XLS 108 KB)
Additional file 7: grapes and papayas in response to ripening. This figure describes the cell function transcripts from tomato, grape and papaya fruits when analysed with the MapMan software.These values (and identities) were the same used to structure the Venn diagram in Figure 2. (PDF 605 KB)
Additional file 8: tomato and grape. This table is a simplified version of Table 2 (edited in Excel program) and Figure 3. The table only returns the results from MapMan analysis of overview ripening regulation of papaya fruit, but there is no comparison of transcription factors between three species. (XLS 37 KB)
Additional file 10: PCR primers. This table describes the efficiency values, as well the standard curves for all primers used in Real Time-PCR experiments. (DOCX 18 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Fabi, J.P., Seymour, G.B., Graham, N.S. et al. Analysis of ripening-related gene expression in papaya using an Arabidopsis-based microarray. BMC Plant Biol 12, 242 (2012). https://doi.org/10.1186/1471-2229-12-242
- Heterologous microarray
- Papaya ripening
- Quantitative gene expression
- Whole genome shotgun
- Transcript profiling