An integrative “omics” approach identifies new candidate genes to impact aroma volatiles in peach fruit
© Sánchez et al.; licensee BioMed Central Ltd. 2013
Received: 21 December 2012
Accepted: 15 May 2013
Published: 23 May 2013
Ever since the recent completion of the peach genome, the focus of genetic research in this area has turned to the identification of genes related to important traits, such as fruit aroma volatiles. Of the over 100 volatile compounds described in peach, lactones most likely have the strongest effect on fruit aroma, while esters, terpenoids, and aldehydes have minor, yet significant effects. The identification of key genes underlying the production of aroma compounds is of interest for any fruit-quality improvement strategy.
Volatile (52 compounds) and gene expression (4348 genes) levels were profiled in peach fruit from a maturity time-course series belonging to two peach genotypes that showed considerable differences in maturation characteristics and postharvest ripening. This data set was analyzed by complementary correlation-based approaches to discover the genes related to the main aroma-contributing compounds: lactones, esters, and phenolic volatiles, among others. As a case study, one of the candidate genes was cloned and expressed in yeast to show specificity as an ω-6 Oleate desaturase, which may be involved in the production of a precursor of lactones/esters.
Our approach revealed a set of genes (an alcohol acyl transferase, fatty acid desaturases, transcription factors, protein kinases, cytochromes, etc.) that are highly associated with peach fruit volatiles, and which could prove useful in breeding or for biotechnological purposes.
Peach (Prunus persica L. Batsch) was definitively positioned as a new fruit model when its genome was sequenced and released in 2010 by an international initiative [1, 2]. This is further supported by the availability of several genetic and genomic tools, including molecular markers, genetics maps, transient fruit expression assays, microarrays, EST databases, and a 9K SNP array [1, 3–7], and by the fact that peach fruit is a drupe and therefore has a different physiology, anatomy, and metabolism from other “post-genomic” fruit crops, such as grape and tomato (berries), strawberry (an aggregate of achenes), orange (speridium), and apple (a pome-type fruit). Hence, the peach represents an excellent opportunity to isolate novel genes related to specific traits, like aroma volatiles. Nevertheless, there still exist many obstacles to discovering gene function in peach. For example, in-silico analyses are still quite restricted when compared to other species, like Arabidopsis, for which more exhaustive phenotypic and molecular database repositories are available. The use of a priori knowledge (e.g., based on co-expression data) to select genes for functional analyses is very limited in peach, which makes it necessary to develop in-house data in order to identify the candidate genes associated with important traits or physiological processes.
Underlying its emerging role as a fruit model is the fact that peach is an important food commodity with an estimated net worldwide production of 11 billion US$ . Aroma is one of the main attributes that affects fruit quality  and has been recognized as one of the main factors that affect peach prices in the market . For this reason, volatile organic compounds (VOCs), which define aroma and, in combination with sugars and organic acids, also contribute to fruit taste, have received a great deal of attention. More than 100 VOCs have been described in peach to date ( and references therein), of which about 25 of them appear to conform the typical peach aroma. In particular, γ- and δ-decalactone play a key role in association with C6 compounds, alcohols, esters, terpenoids, and phenolic volatiles . In addition to their contribution to fruit quality, peach volatiles are also important in the food and fragrance industry, where they are used as flavoring agents. Indeed, γ-decalactone is a sought-after industrial product that confers a “peach-like” odor  with an expanding annual world demand estimated at 10,000 Kg in 1997 . Despite the importance of lactones, their biosynthetic pathways in peach, and in plants in general, are still poorly understood . An early study suggested that epoxide hydrolases were involved in lactone production, since it was observed that nectarines (a glabrous mutation of peach) are able to produce an artificial lactone when infiltrated with a synthetic, radiolabeled epoxy acid . The analysis of EST libraries later showed that a homologous gene to epoxyde hydrolases was expressed in peach skin , although this gene has not been further characterized. Indeed, no gene involved in volatile production in peach has been reported to date. Most studies on genes related to peach aroma have focused on analyzing genes whose homologs are characterized in other plant species, i.e., literature-derived candidate genes. For example, Vecchietti et al.  analyzed an EST library to show that a set of candidate genes was expressed in peach fruit and could thus be related to the formation of different volatile compounds. Another study targeted certain members of the carotenoid cleavage dioxigenase gene family for an expression analysis of genotypes differing in carotenoid accumulation to support their involvement in the production of norisoprenoid volatiles in peach . The identification of QTLs for volatile compounds of peach was recently reported . The study in question also proposed putative candidate genes for minor contributors to peach aroma, such as linalool, p-menth-1-al, and nonanal, based on low-resolution co-localization of candidate genes within major QTL regions for these volatiles compounds . The use of “omics” technologies is needed to increase the resolution and accuracy of candidate-gene approaches in peach. Biosynthetic pathways for the main volatile compounds in peach are still poorly understood , and given the peach-specific nature of the volatiles involved, information from other model species may not be sufficient.
Integrating transcriptomic and metabolomic data by means of pair-wise correlation analyses is emerging as a promising approach for discovering novel gene functions, despite its currently being limited to the model plant Arabidopsis thaliana (reviewed in ) and tomato [21, 22]. On the other hand, integrative “omics” approaches have been widely employed for the purpose of studying diverse physiological processes of plants [23–27] rather than discovering candidate genes.
In order to identify the key genes underlying the production of aroma compounds in peach, we conducted an integrative analysis using in-house-developed transcriptomic and metabolomic data sets. Our reductionist approach, based on correlation analysis [Hierarchical Cluster Analysis (HCA), Pearson Correlation (PC), and Correlation Network Analysis (CNA)] between gene expression and metabolite accumulation levels, allowed us to propose candidate genes for aroma control in peach. Moreover, we cloned and expressed in yeast one of the candidate genes identified herein to analyze its enzymatic activity. Our results demonstrate that the encoded protein possesses ω-6 Oleate desaturase activity, which might participate in the biosynthesis of esters and lactones in peach fruit by generating a common precursor, linoleic acid. In addiction, our results provide some insights into the physiology and metabolism of peach, which are discussed. To the best of our knowledge, this is the first study to use non-targeted approaches to identify candidate genes in peach. The set of genes related to aroma compounds provided herein could prove useful for marker-assisted breeding and/or biotechnology improvements of peach fruit.
Fruit material and analysis
Peach fruit (genotypes ‘MxR_01’ and ‘Granada’) were harvested in July 2009 from a local commercial orchard situated in Murcia, Spain. ‘MxR_01’ is a freestone melting-flesh peach which was obtained through the IVIA (Instituto Valenciano de Investigaciones Agrarias) breeding program. ‘MxR_01’ is a seedling from the cross between ‘RedCandem’ and ‘Maruja’, a traditional Spanish germplasm. ‘Granada’ is a clingstone non-melting peach low chilling obtained from a Brazilian breeding program. Fruits corresponding to four maturity stages (S1, S2, S3, and S4) were harvested from one tree per genotype. For S1, S2, and S3, ten fruits from each maturity stage were collected and maturity parameters were immediately analyzed. For S4, “harvest ripe” (or ′ready to buy′) stage, 20 fruits were harvested and divided into two groups. One group (denominated S4) was immediately analyzed and the other (denominated S4+SL) was subjected to shelf-life conditions (2 days at 20°C, 85% RH) prior to the maturity analysis. S4+SL corresponded to “consumption ripe” (or “ready to eat”) stage. The maturity parameters (external color, flesh firmness, weight, total soluble solids (SSC), and ethylene and CO2 production) were analyzed as described in Sanchez et al. . A Principal Component Analysis was carried out with maturity parameters in order to identify the three most homogeneous fruits per maturity stage for each genotype (data not shown). These three fruits (biological replicates) were selected for volatile and microarray analyses.
Sample preparation and HS-SPME-GC-MS conditions
Frozen samples of fruit mesocarp were ground to powder in liquid nitrogen and used for volatile and microarray analyses as follows. Volatile compounds were analyzed from 500 mg of frozen tissue powder, as previously described . The volatile analysis was performed on an Agilent 6890N gas chromatograph coupled to a 5975B Inert XL MSD mass spectrometer (Agilent Technologies). For the chromatography and mass spectra conditions, see Sánchez et al. . A total of 36 commercial standards were used to confirm compound annotation; the VOCs confirmed are listed in Additional file 1: Figure S1. Volatiles were quantified relatively by means of the Multivariate Mass Spectra Reconstruction (MMSR) approach developed by Tikunov et al. . A detailed description of the quantification procedure is provided in Sánchez et al. .
RNA was extracted from 3 g of frozen tissue powder as described by Meisel et al. . RNA quantity and purity were determined spectrophotometrically with Nanodrop (Nanodrop Technologies Inc.; http://www.nanodrop.com/). RNA integrity was verified by agarose gel electrophoresis. After the quantity and quality checks, RNA were used for the microarray analyses and the qRT-PCR as described below.
Microarray hybridization and scanning
The microarray analysis was performed essentially as described in Ogundiwin et al. . For the microarray hybridization, the RNA from the samples and the reference pool (1μg of each sample) were amplified and aminoallyl-labeled using the MessageAmp™ II aRNA Kit (Ambion, http://www.ambion.com) and 5-(3-aminoallyl)-20-deoxyuridine-50-triphosphate (aa-dUTP, Ambion), following the manufacturer’s instructions. For each RNA sample, 7.5 μg of aminoallyl-labeled aRNA were re-suspended in 0.1 M Na2CO3 (pH 9.0) and labeled with either Cy5 or Cy3 Mono NHS Ester (CyTM Dye Postlabelling Reactive Dye Pack, Amersham), respectively. Samples were purified with MegaclearTM (Ambion) following the manufacturer’s instructions. Incorporation of Cy5 and Cy3 into probes was measured with a Nanodrop spectrofluorometer (Nanodrop Technologies Inc., http://www.nanodrop.com/). Microarray hybridization of samples and references to the ChillPeach microarray slides was performed manually using Telechem Hybridization Chambers (Corning), following the manufacturer’s instructions. After hybridization, slides were washed in 2x SSC, 0.1% SDS for 5 min at 42°C, 0.1x SSC, 0.1% SDS for 10 min at room temperature, 0.1x SSC for 5 min at room temperature 4 times, and 0.01x SSC for 5 min at room temperature 4 times. Arrays were drained by centrifugation at 528g for 2 min. Slides were scanned with a GenePix 4000B scanner (Axon Instruments) at 10 μm resolution, 100% laser power and with different PMT values to adjust the ratio to 1.0. Microarray images were analyzed and globally normalized using the GenePix 4.1 software (Axon Instruments). Only the spots with a background-subtracted intensity greater than 2-fold the mean background intensity in at least one channel were selected for the analysis. Data files were imported into Acuity 4.0 (Axon Instruments) and normalized by the Lowess normalization method. Finally, only the spots with valid values in at least 2 of the 3 analyzed hybridizations were considered for further analysis. The means and standard deviations of the values were calculated from each sample as log2 values and were later normalized to the median of the reference pool.
Real-time qRT-PCR analysis
One microgram of total RNA was used to synthesize first-strand cDNA using the SuperScript first-strand synthesis system for RT-PCR (Invitrogen). Two microliters of diluted cDNA (100 ng/μL) were used for qRT-PCR using the SYBR Green PCR master mix (Applied Biosystems), following the manufacturer’s recommendations, and an ABI Prism 7000 sequence detection system (Applied Biosystems). Each biological replicate was assayed in triplicate. Gene-specific oligonucleotide primers were designed using the Primer Express_ version 2.0 software (Applied Biosystems). Primer information is available in Additional file 2: Table S1. The expression levels for target genes were calculated in relation to a reference gene by the DDthreshold cycle method (DDCt, Applied Biosystems). From the microarray data, the gene (ID: PPN078E12, http://bioinfo.ibmcp.upv.es/genomics/ChillPeachDB), whose profile was corroborated as being constant throughout both time courses, was selected as a reference for the normalization of all the subsequent qRT-PCR analyses. The relative gene expression for each candidate gene was expressed by means of the DDCT method using the expression value of the reference gene to normalize expression and the sample from stage S1 in each time-course series as reference samples.
The Acuity 4.0 software (Axon Instruments) was used for: hierarchical cluster analysis, heatmap visualization, principal component analysis, Pearson correlation evaluation, and for the Student’s t-test for significant differences of volatile levels.
To detect differentially expressed genes between fruit at harvest (S4) and after shelf-life simulation (S4+SL) in both genotypes, data were analyzed with the SAM (Significance Analysis of Microarray) package . Statistical significance was assessed using a two-class (unpaired) SAM analysis, with a false discovery rate of 5% and a q-value of < 0.05.
Correlation network analyses were conducted with the Expression Correlation (http://www.baderlab.org/Software/ExpressionCorrelation) plug-in for the Cytoscape software . Network topological parameters were calculated with the NetworkAnalyzer plug-in . Networks were visualized with the Cytoscape software, v2.8.2 (http://www.cytoscape.org).
Venn diagrams were drawn with Microsoft PowerPoint.
Cloning and bioinformatics analysis of the peach candidate gene
The cDNA (from ‘Granada’ fruit at S4+SL) synthesized by qRT-PCR was used as a template for cloning the ORF of candidate gene PP1002E07. Coding sequences were amplified by PCR (FAD_F: ATGGGTGCCGGTGGAAGAAT and FAD_R: TTATAACTTATTATTGTACC) using Taq DNA polymerase (Biotools B and M Labs SA, Spain). PCR products were cloned in pCR®8 TOPO (Invitrogen, Spain), according to the manufacturer’s instructions, in order to create the pEntry-PpFAD_1B-6 vector. Cloned ORFs were verified by sequencing both DNA strands.
For the sequence analyses, different tools were used. Protein sequences were aligned by the clustalW method using the MegAlign software (DNAStar). Transmembrane domains were predicted with TMpred (http://www.ch.embnet.org) and protein localization with ProtComp v9.0 (http://linux1.softberry.com/berry.phtml).
For yeast expression vector construction, the coding sequence was PCR-amplified from the pEntry-PpFAD_1B-6 vector with forward primer FAD-Kozack_F (GGCATGGGTGCCGGTGGAAGAAT) and reverse primer FAD_R. FAD-Kozack_F includes three non-template nucleotides at the 5′ end (underlined) to improve protein synthesis in yeast. The PCR product was cloned in pYES2.1 TOPO® vectors (Invitrogen) by TA cloning following the manufacturer’s instructions. After checking the insert by sequencing, the plasmid was introduced into Saccharomyces cerevisiae strain W303-1A MATa [leu2-3112 trp1-1 can1-100 ura3-1 ade2-1 his3-11,15} for expression analysis as follows.
Expression of the peach candidate gene in yeast and fatty acid analysis
Yeast cells were grown on a rotary shaker at 200 rpm at 30°C in synthetic defined (SD) medium containing raffinose as a carbon source and mixtures of amino acids and nucleoside precursors for marker selection. The expression of the transgenes was induced by the addition of galactose to 2% (w/v) or 1% glucose as a negative control. After 3 days of induction at 20°C, cells were collected, washed with water, and dried. Fatty acids were extracted and methylated from cell pellets and the latter were silanized to obtain trimethylsilyl derivatives, as previously described . The methyl-ester derivatives were analyzed by gas chromatography (GC) with a Hewlett–Packard 6890 gas chromatograph (Palo Alto, CA, USA). The column, chromatographic, and detection conditions are described in Venegas-Calerón et al. .
As part of our final goal of identifying genes and sources of variability to improve peach quality, we undertook complementary genetics and genomics approaches. ‘MxR_01’ and ‘Granada’ peach genotypes vary greatly for important traits, such as melting/non-melting, freestone/clingstone, chilling requirement, aromas, and fruit flavour. To exploit this variability, an F1 population was developed to analyze quantitative trait loci, which will be presented elsewhere. Here we present together the analyses of gene expression and volatile accumulation during fruit maturity and ripening of fruit of the parental genotypes in order to find aroma-related genes.
Physiological and shelf-life ripening of peach genotypes
In order to characterize the ripening stages of the ‘Granada’ and ‘MxR_01’ peach genotypes, typical maturity parameters were evaluated and are presented in the supplementary data (Additional file 3: Figure S2). The results indicate that, as expected, peel ground color and weight increased with fruit ripening, whereas flesh firmness decreased, but with differences between varieties. Soluble solids content (SSC) was not affected throughout the study period, which indicates that this parameter is not a suitable indicator of ripening for either of the genotypes analyzed. Ethylene and CO2 production was also monitored, as the evolution of these compounds reflects the physiological ripening stage of fruits and could also underlie the differences observed in ripening. Both genotypes, ‘Granada’ and ‘MxR_01’, showed increased ethylene production at mature stages (S3 for ‘Granada’ and S4 for ‘MxR_01’), which is typical of climacteric fruit with differences in degrees (Additional file 3: Figure S2). Shelf-life ripening was also used to increase the complexity of our data set and to evaluate the effect on volatile and gene expression. Storage at 20°C for 2 days affected peel ground color, firmness, and ethylene production in both genotypes (Additional file 3: Figure S2). CO2 production decreased in both genotypes with shelf-life simulation, although differences were only significant for the ’MxR_01’ genotype. These results indicate that our postharvest treatment was effective in stimulating ripening off the tree in the fruit of both genotypes and that, all together, our samples represent different stages of fruit development and ripening that can be interrogated to obtain the network of interactions between transcripts and volatiles.
Non-melting (‘Granada’ genotype) and melting (‘MxR_01’ genotype) peaches showed different volatile evolution patterns during ripening
Clusters C1, C2, and C3 are formed by lactones and some lineal esters (Figure 1A). They showed similar overall trends in which compound levels increased in both genotypes, but with different profiles and fold changes (Figure 1B). In ‘Granada’, the volatiles in these clusters reached high levels at S3 and remained high for the remaining time series (mainly for C1 and C2), while in ‘MxR_01’, the increment observed throughout the time-course series was gradual (mainly for C2 and C3).
Clusters C4, C5, and C6 are formed mostly by terpenoid volatiles. They increased during ripening in ‘Granada’, while in ’MxR_01’ they showed no changes, or even decreased (C5) during ripening.
Ethanol and Ethyl acetate (C7) levels rose at a mature stage (S4) and even reached higher levels with shelf-life conditioning (64-fold as compared to S1) for ‘Granada’, while in ‘MxR_01’, its content increased only after shelf life reached a 2-fold difference as compared to S1.
Volatiles that decreased during ripening (II) also showed different evolution profiles according to genotype (Figure 1B). Cluster C10 is composed of aldehydes with eight and nine carbons (Octanal, Nonanal, and 2-Nonenal), which showed moderate changes during ripening (up to 1.5-fold changes in relation to S1, Figure 1B). Clusters C11 and C12 are formed by the volatiles derived from the catabolism of linoleic and linolenic acids, the so-called “green compounds” (Figure 1A). Both clusters displayed a similar decreasing trend in both time-course series (Figure 1B). The phenolic volatiles Benzyl chloride, Benzaldehyde, and Benzeneacetaldehyde (C13) decreased from S1 to S2, and remained at low levels until the end of the time series for ‘Granada’, while they showed no significant changes in the ‘MxR_01’ time series.
Several aroma-related volatiles accumulated at different levels in ‘Granada’ and ‘MxR_01’ when fruits reached the commercial mature stage (S4), which means that they could display very marked differences in aroma. As shown in the supplementary data (Additional file 4: Figure S3), 13 volatiles showed significant differences, with most (12) exhibiting higher levels in ‘Granada’ (the more aromatic variety). Compound γ-jasmolactone, which has a characteristic peach-like odor, showed a 39-fold higher level in ‘Granada’ when compared to ‘MxR_01’. Other compounds with pleasant aroma descriptions included 2-Hexen-1-ol, acetate (E) (fruity), γ-Hexalactone (coconut), and 3-Hexen-1-ol acetate (Z) (fruity), which also had higher levels (between 3-fold and 5-fold) in ‘Granada’ compared to ‘MxR_01’. The only compound with higher levels in ‘MxR_01’ at S4 was Benzeneacetaldehyde (4-fold higher than ‘Granada’ at S4), which has been described as confering a “green” aroma.
The ‘MxR_01’ genotype showed an enhanced shelf-life response
Due to commercial considerations (handling, long-distance transport, storage, etc.), peaches are generally harvested before their complete maturity, but the subsequent shelf life allows the fruit to ripen to the minimum quality threshold for consumers’ acceptance. To analyze the shelf-life response of the ‘MxR_01’ and ‘Granada’ genotypes, we compared the changes in volatile content and gene expression.
To analyze the effect of shelf-life treatment on the transcriptome, we conducted a PCA analysis with the 4348-gene data set and studied the differentially expressed genes between S4 and S4+SL in both genotypes. The first component (PC1, explaining 42% of variance) separated samples according to genotype (Figure 2A, right). The second component (PC2) explained 15% of variance and separated between the S4 and S4+SL replicates of the ‘MxR_01’ genotypes. While the ‘Granada’ samples remained close together in the first 2 components space, the ‘MxR_01’ samples separated according to treatment (Figure 2A, right), indicating that shelf-life has a greater impact on gene expression in the ‘MxR_01’ genotype than in the ‘Granada’ one. Accordingly, the direct analysis of the differentially expressed genes between S4 and S4+SL in ‘Granada’ and ‘MxR_01’ revealed a stronger gene expression response to the shelf-life in the ‘MxR_01’ genotype (Figure 2B, right). By taking a False Discovery Rate (FDR) and a q-value of < 0.05 as criteria, we found that for ‘Granada’, 13 genes were differentially expressed after treatment (nine up-regulated and four down-regulated), while the ‘MxR_01’ genotype fruits showed a drastic change in gene expression after shelf-life simulation, with 794 differentially expressed genes (451 up-regulated and 343 down-regulated). Only four genes (three up- and one down-regulated) showed the same trend in both genotypes (Figure 2B, right, Additional file 6: Table S3). The gene showing the largest changes after shelf-life in both genotypes (id: PPN044F04) had a calcium-binding EF hand domain. A gene related to the metabolism of amino acids (Aspartate aminotransferase) showed moderate changes in both genotypes, increasing after treatment. The third gene which increased after shelf-life (id: PPN064C01) is a homolog to a translation initiation factor, suggesting that protein synthesis could be stimulated by treatment in both genotypes. Only one gene decreased after shelf-life in both genotypes, an acid phosphatase belonging to class B. The ‘MxR_01’ genotype showed a higher decrease in the expression of this gene (3.2-fold) compared to ‘Granada’ (1.8-fold).
These results support that mostly ‘MxR_01’ fruits undergo significant molecular, and consequently physiological, changes during postharvest conditioning, and that very little happens to ‘Granada’. It is worth mentioning that despite the larger number of genes showing significant differences in expression after shelf-life in ‘MxR_01’ (794 genes vs. 13 genes in ‘Granada’), most involved only minor fold-change differences. Of the genes, 9% (74:794) showed a fold change of over 2 (between S4 and S4+SL) for the MxR genotype. For ‘Granada’, nine of the 13 differentially expressed genes showed fold differences over 2. In the supplementary data (Additional file 7: Table S4), a complete description of differentially expressed genes with a 2-fold change cutoff is provided.
Identification of the genes related to peach volatile compounds by hierarchical cluster analyses
In order to obtain a global view of the interrelationships between volatile accumulation and gene expression, a Pearson correlation analysis was conducted for all the possible volatile gene pairs (Additional file 8: Table S5). Volatiles showed wide variability in the range of the distribution of volatile-gene correlations (data not shown). For example, Hexanoic acid 2-ethyl- had correlations ranging from −0.42 to 0.52, while γ-Jasmolactone showed correlations ranging from −0.92 to 0.96. Thus, our data set contains some volatiles which have no clearly associated genes, while others show many strongly correlated genes.
Genes that highly correlate with volatile compounds proposed as candidate genes
Genes for volatile cluster C1
Acid phosphatase 1 precursor (EC 18.104.22.168) (Apase-1(1))
CAMP response element binding (CREB) protein
Tub family, putative
Genes for volatile cluster C2
Prunus_persica transcript; similar to CRHB10 [Ceratopteris richardii (Triangle waterfern)]
Cytochrome P450 monooxygenase CYP72A59
HCA, PC, CNA
Phosphatidylinositol transfer protein/retinal degeneration b protein
Genes for volatile cluster C3
Omega-6 fatty acid desaturase
Pyruvate decarboxylase 1
Vesicle-associated membrane protein 727
AJ533535 S3II Prunus persica cDNA clone PP_S3II_T5_SP6, mRNA sequence
Serine-threonine protein kinase, putative
HCA, PC, CNA
HCA, PC, CNA
HCA, PC, CNA
Malus_x_domestica transcript; similar to F17O7.2 [Arabidopsis thaliana]
Novel plant SNARE 13
Calcium binding protein
Genes for volatile cluster C4
Genes for volatile cluster C5
AJ827262 S3 Prunus persica cDNA clone S312E11, mRNA sequence
Putative pod-specific dehydrogenase SAC25
Putative pod-specific dehydrogenase SAC25
Genes for volatile cluster C7
PP_LEa0004N09f Peach developing fruit mesocarp
Genes for volatile cluster C11
Desaturase delta 9
Genes for volatile cluster C12
PU1_plate44_F06 PU1 Prunus persica cDNA weakly similar to putative protein (AL161573)
GTP-binding protein SAR1B
Putative acid cluster protein 33
Genes for volatile cluster C13
Cysteine protease CP1
Putatative tyrosine aminotransferase
Genes related to the main aroma compounds of peach
Although aroma is a complex trait, and each volatile contributes to defining overall peach aroma, some compounds appear to be more important than others according to sensory assessments (Additional file 1: Figure S1). For example, the aroma of two lactones, γ-Decalactone and γ-Jasmolactone, is described as peach-like . We have evaluated correlations between the levels of gene transcripts and the levels of these two aroma volatiles across the different samples in order to find genes related to them that were not identified in the previous analysis. Additional file 11: Table S7 provides the ten genes with the highest direct correlation values; in addition, the full list of genes with correlation values between the corresponding transcript and volatile levels is included in the supplementary data (Additional file 8: Table S5). As expected, some of the genes identified in the HCA were also found in this analysis. For example, genes related to ethylene signaling proved to be highly correlated with γ-decalactone (PPN054G06 and PPN004H06) and γ-Jasmolactone (PP1005G06). Cytochrome P450 (PPN070H11) was strongly correlated with both lactones. γ-Jasmolactone also showed high correlations with a protein involved in metabolite transport (PPN032F06) and with two protein kinases (PPN069F09 and PPN010B11). Nevertheless, two other genes that were not previously identified (PPN059A01 and PPN026D01) with no homology in Arabidopsis showed the highest correlations with γ-decalactone (r=0.90 and r=0.87, respectively). Moreover, γ-Jasmolactone was highly correlated with the genes showing no homology with Arabidopsis (PPN002B03 and PPN031G12), and which were not identified in previous analyses. Interestingly, some genes involved in endoplasmatic reticulum traffic were also associated with lactone production in peach fruit. A homolog to a Phosphatidylinositol transfer protein (PPN037A04) correlated well with γ-decalactone, while a novel plant SNARE 13 (PPN031F12) highly correlated with γ-Jasmolactone. Finally, γ-Jasmolactone showed a good correlation with the level of a transcript that encodes a protein related to the caleosin family (PP1004C02).
Correlation network analysis for discovering hub genes
As mentioned previously, in order to analyze interactions with genes, a correlation network was constructed by preselecting those genes with a cutoff of r > 0.85 in absolute values. To integrate the interactions between genes and VOCs, the two networks were merged to create a new one (Figure 4A, bottom). The resulting network consisted of 160 nodes (52 VOCs and 108 genes) and 1364 edges (correlations). The network clearly reflected the different distributions of the correlations between VOCs and genes observed. As the volatiles from C1, C2, and C3 showed the largest number of strong correlations with genes (both positive and negative), their nodes were highly connected to gene nodes (Figure 4A, bottom). Compound γ-jasmolactone showed strong correlations (r > 0.95) with a group of 14 genes (sub-cluster A, Figure 4B). Some members of this sub-cluster had already been selected as candidate genes for the metabolites of C3 (PPN069F09, PPN032F06, PPN031F12, PPN036E10, PPN031G12, PPN010B11, and PP1004C02) or for lactones of C2 (PPN070H11), as described in previous sections. Nevertheless, one gene, which our previous analysis did not reveal, was selected as a candidate gene (PPN051H01, Table 1). Sub-cluster A was related to a group of 20 genes (sub-cluster A.1) by strong positive correlations. It is interesting to note that some of the genes of sub-cluster A are involved in signal transduction processes (Additional file 12: Table S8). For example, there are two putative protein kinases (PPN069F09 and PPN010B11) that show high connectivity (degree=39 and 40, respectively), as well as two ethylene signaling genes (PP1005G06 and PPN004H06). A gene with homology to Asparagine synthetase (PPN079G06) is strongly and inversely correlated with γ-jasmolactone and could be considered a hub given its high connectivity (degree=37, Additional file 12: Table S8, Figure 4B). Volatile cluster C1 was highly related to sub-cluster B by strong correlations with 4-Methyl-5-penta-1,3-dienyltetrahydrofuran-2-one and 3-Hexen-1-ol acetate (Z)- (Figure 4B). Two transcription factors belong to this sub-cluster: one had already been selected (PP1009A02) and the other was a newly identified one (PPN066C05, Table 1). Furthermore, a gene that was not identified previously as showing homology to ripening-related proteins was selected (PPN066B05, Table 1). Sub-cluster B is highly interconnected with sub-cluster C, which consists of three genes with strong correlations with the lactones of C2 (Figure 4B). The members of this sub-cluster are: a gene related to Gibberellin metabolism (Gibberellin 2-oxidase, PP1004A06), a gene that is likely related to cell wall physiology (UDP-arabinose 4-epimerase 1, PPN062D06), and a gene with no homolog in Arabidopsis (PP1004G05), which had already been identified (Additional file 12: Table S8). The main sub-clusters that correlated with lactones and esters (A, B, and C) were highly interconnected to a group of 13 genes that formed sub-cluster D (Figure 4A bottom, Additional file 12: Table S8).
With the selected cutoff value (genes with a correlation higher than 0.85), the correlation network analysis failed to identify genes related to the other VOC clusters (C4 to C13). To gain insight into the genes associated with these volatile compounds, a new data set was composed by selecting genes after lowering the cutoff to > 0.8 for volatiles belonging to clusters C4 to C13, which allowed a new correlation network to be constructed (Additional file 13: Figure S5). The aromatic VOCs in C13 are related to a putative tyrosine (Tyr) aminotransferase (PPN018G03) through a direct correlation with Benzeneacetaldehyde (Additional file 13: Figure S5B) and, therefore, was selected as a candidate gene (Table 1). The VOCs from C4 are correlated with a group of 5 genes, two of which are related to hormone signaling. One is an Auxin-responsive protein (IAA16, PP1009D02) and the other (PP1000F05) belongs to a family of proteins regulated by gibberellins. These 5 genes correlated well with a group of 42 genes (sub-cluster E), and some were also associated with auxin (PPN044D01 and PP1005H05) and gibberellin (PPN037E06) signal transduction pathways (Additional file 13: Figure S5A). In more detail, we see that the lipid-derived compounds Furan, 2-pentyl-, and Hexanal inversely correlated with a lipid delta 9 desaturase homolog (PPN052H12), which, in turn, was highly correlated with a gene with no homology in Arabidopsis (PPN023E05). A BZIP-like transcription factor (PPN019F11) strongly correlated with PPN023E05 and with some genes of the sub-cluster E. Compound cis-Linaloloxide (belonging to C5) also strongly correlated with a gene that had a short-chain dehydrogenase/reductase (SDR) domain (PPN053G07). This gene formed a sub-cluster with two other genes, one which also had a short-chain dehydrogenase/reductase (SDR) domain (PPN078H04) and the other with no homologs in Arabidopsis thaliana (PPN029H12). The previous selection strategy also revealed these three genes as candidates (Table 1).
Validation of microarray data by qRT-PCR analysis
In order to validate the expression profile of the candidate genes identified by microarray analysis, gene-specific qRT-PCR analyses were conducted. We selected a subgroup of genes that we identified by their association with the main aroma-contributing volatiles, i.e., lactones and esters. Seven genes selected for clusters C1, C2, and C3 (PP1002E07, PPN070H11, PPN001H09, PPN032F06, PPN002B03, PPN059A01, and PPN066B05) were analyzed. Three genes that were inversely correlated with lactones/esters (PP1002D12, PPN079G06, and PPN008D02, Additional file 8: Table S5) were also included to obtain a better validation of our results. The Pearson correlation coefficient between the expression results obtained with both RNA profiling technologies was evaluated for all the genes studied. The ten genes analyzed by transcript-specific qRT-PCR analysis corroborated the microarray analysis, showing correlation coefficients ranging from 0.710 to 0.988 (Additional file 2: Table S1). By way of example, the comparison between the qRT-PCR and the microarray profiles of gene PPN001H09 is shown in Additional file 14: Figure S6.
The enzyme encoded by candidate gene PP1002E07 displays ω-6 Oleate desaturase activity
Non-melting (‘Granada’ genotype) and melting (‘MxR_01’ genotype) peaches show different maturity/ripening and shelf-life responses, which can be exploited to analyze co-regulation patterns
An “omics” approach based on integrative transcriptomics- and metabolomics-derived data was undertaken to identify those genes putatively involved in volatile production in peach fruit. The main aroma-contributing volatiles of peach increase during maturity and ripening following different patterns [35–38]. It is also known that the peach transcriptome undergoes significant re-organization during fruit ripening and maturity [39, 40]. Thus, our rationale was that the genes involved in aroma-related volatile production could be discovered by analyzing the co-regulation between gene expression and volatile accumulation during maturity and ripening. In addition, we hypothesized that the inclusion of different peach types could improve the robustness of metabolite-gene relationships since the gene-volatile correlation should hold across the genotype with different ripening characteristics. To this end, four maturity stages were analyzed for clingstone non-melting (‘Granada’) and freestone melting-flesh (‘MxR_01’) peaches (Additional file 3: Figure S2). To add further complexity to the data set and to increase the robustness of our model, one postharvest treatment was also included, since it is known that non-melting and melting peaches respond differently to shelf-life ripening . The volatile profiles analyzed in both time-course series confirmed that these genotypes show a different volatile evolution during maturity and ripening (Figure 1), suggesting that our hypothesis was correct. Each genotype also showed a different response to shelf-life simulation. Volatile profiles were more affected after treatment in ‘MxR_01’ as compared to ‘Granada’, although the lactones and esters in both genotypes increased after treatment (Figure 2), which is in agreement with previous studies [37, 38]. Correspondingly, ‘MxR_01’ also showed a more dramatic change in gene expression as compared to ‘Granada’ after the shelf-life simulation (Figure 2), indicating that the restructuration of volatile content after shelf-life could be, at least in part, determined by a modification in gene expression.
Therefore, our approach combined different genetic backgrounds with a development-regulated process and an artificial post-harvest treatment in order to analyze co-regulation patterns. This strategy differs substantially from other previously used integrative approaches. In tomato, the change in gene expression and metabolite accumulation throughout fruit development was also used to identify not only candidate genes, but also variability between different fruit tissues (mesocarp and locular), which was included to improve the robustness of the study . Carrera et al.  applied an alternative strategy consisting in profiling (transcripts and metabolites) an RIL population to find candidate genes that affect tomato fruit aroma by modeling the “omics” data obtained. In this case, system perturbation was given by the genomic region introgressed in each line. In other study, Arabidopsis plants were starved in order to affect the transcriptome and the metabolome so that they could analyze the co-regulation patterns and thereby identify the genes associated with different metabolites [42, 43]. The variation in gene expression and the volatile levels obtained in our study suggest that the strategy used herein could be a suitable alternative for the analysis of co-regulation patterns.
Lactone and ester production in peach fruit requires a tight regulation of lipid catabolism
A number of candidate genes for peach volatiles were identified using a combination of data analysis techniques based on correlations between gene expression and volatile accumulation data (HCA, PC, and CNA). The functional annotation indicated that, in some cases, this is due to a general effect on ripening. For example, the genes related to ethylene biosynthesis and perception (PPN004H06, PP1005G06, and PPN054G06) directly correlate with lactones (Figure 3, Figure 4, Additional file 11: Table S7), suggesting that they could be regulators of their biosynthesis. It is known that ethylene is the hormone that controls the ripening of climacteric fruits, such as peach , which makes it an unsuitable candidate for specific aroma control since its modification could alter the whole ripening syndrome. However, the possibility of there being downstream effectors specific for aroma production, which would be much more ideal candidates, cannot be ruled out, as different aspects of ripening appear to be controlled by specific transduction elements . Overall, these results are not only in agreement with Zhang et al. , who indicated that ethylene plays a regulatory role in aroma formation in peach, but also provide additional support for further research to search for the specific ethylene signal transduction elements involved in the specific activation for volatile production during ripening.
Unlike esters, the lactone biosynthesis pathway in plants is still incomplete. Nevertheless, it seems clear that lactone biosynthesis starts from fatty acids with the introduction of an O atom to form hydroxy fatty acids . Hydroxy fatty acids are then shortened by β-oxidation to form 4- or 5- hydroxy acids, which in turn result in the corresponding lactones after intramolecular esterification (Figure 6A, ). However, it is unclear how the introduction of the O atom is achieved. By infiltrating a synthetic radiolabeled epoxy acid into fruits, it has been demonstrated that nectarines (peach glabrous mutants) are able to produce an artificial lactone . Based on this observation, it was proposed that the introduction of the O atom is achieved by fatty acid epoxydation and that the subsequent breakdown of the epoxy group (by epoxyde hydrolase) forms hydroxy acids; unfortunately, no further evidence has been presented. Alternatively, the introduction of hydroxyl groups could be produced by the hydration of unsaturated fatty acids or the direct hydroxylation of fatty acids catalyzed by Cytochrome P450s (CYPs) or by other hydroxylases not related to CYPs (Figure 6A). To date, neither have the proposed endogenous fatty acid modification enzymes been cloned nor has their involvement in lactone biosynthesis in peach fruit been demonstrated.
For volatile cluster C1, an AAT homolog (PPN001H09) was selected as a putative candidate (Table 1, Additional file 10: Table S6) because it is highly associated with the lactone of C1 in the correlation network analysis (Figure 4). Moreover, the expression of this gene is moderately to highly correlated to the other lactones (from r=0.64 to 0.86, Additional file 8: Table S5). It is still unclear whether the last cyclation step leading to lactone biosynthesis requires an enzyme or not , but if this were the case, this gene could represent a novel type of AAT that produces cyclic esters (lactones). This AAT (PPN001H09) also correlated with lineal esters (r=0.69 to 0.75, Additional file 8: Table S5), making it therefore possible that this gene encodes a standard AAT. Nonetheless, a phylogenetic analysis has revealed that the AAT found here (PPN001H09) is not related to previously characterized fruit AATs (data not shown). As it is unlikely that this gene can be identified by candidate gene approaches based solely on sequence similarity to the genes reported in the bibliography, these results emphasize the advantage of using non-targeted analyses based on correlation rather than more targeted gene approaches based solely on homology with genes previously described as being associated with a given process. To date, no AAT from peach has been characterized, although an EST analysis has revealed that some AATs are expressed in fruit mesocarp . Moreover, a targeted approach has demonstrated that the expression of another AAT (ppAAT1) did not correlate with ester and lactone production during postharvest peach ripening . We can confirm this result as ppAAT1 is represented in our microarray (PPN018D08) and it poorly or moderately correlates with lactones and lineal esters (r between 0.27 and 0.60, Additional file 8: Table S5).
Correlation analyses revealed a CYP, PPN070H11, which is highly associated with the lactones in C2 (Table 1). CYPs represent the largest family of genes involved in plant metabolism, which include enzymes that are capable of epoxidizing or directly hydroxylating fatty acids. A CYP77A4 (belonging to the CYP77 subfamily) from Arabidopsis catalyzes the epoxidation of free oleic acid (18:1) to form 9,10 epoxyestearic acid , which has been proposed to be the precursor of γ-dodecalactone in peach . Other CYPs (belonging to the CYP71 family), which are capable of epoxydating unsaturated C18 fatty acids (18:2 and 18:3) as phosphatidylcholine esters, have been characterized in a Euphorbia species . Moreover, the position of the introduction of the hydroxyl group is crucial to producing a long hydroxy fatty acid, which, after β-oxidation, will produce the proper hydroxy acid molecule. The CYP703A2 (CYP703 subfamily) from Arabidopsis can catalyze the hydroxylation of saturated medium chain fatty acids (C10-C14) with a preference for hydroxylate at position C7 , and can therefore produce hydroxy fatty acids, which could be precursors for lactone generation after β-oxidation and cyclation. According to phylogenetic analyses (data not shown), the CYP identified here belongs to the CYP72A subfamily. CYP72A subfamily members are quite diverse, and most of them have unknown functions . The CYP72 clan is also associated with the metabolism of other fairly hydrophobic compounds besides fatty acids, such as isoprenoids. Other members have been associated with the catabolism of hormones (brassinosteroids and gibberellins) or with the biosynthesis of cytokinins . Thus, we cannot exclude the possibility that this candidate gene (CYP, PPN070H11) could be involved in lactone production by, for example, controlling a hormone metabolism instead of participating directly in lactone biosynthesis.
The HCA also identified a gene, PP1002E07, with a homology to fatty acid desaturase (FAD), as being highly associated with 2-Hexen-1-ol acetate (E)- (Additional file 9: Figure S4, Table 1). Nevertheless, this gene also correlated well with lactones (r= 0.65 to 0.89, Additional file 8: Table S5) and, in some cases, has correlation coefficients higher than those of 2-Hexen-1-ol, acetate (E)- (r=0.73, Additional file 8: Table S5). The Arabidopsis homolog of this gene is an ω-6 fatty acid desaturase, which catalyzes the reduction of oleic acid (18:1) esterified to the sn-2 position of the membrane lipid phosphatidylcholine (PC) to linoleic acid (18:2). Yet desaturase genes also share a good homology with other fatty acid modification enzymes, like hydroxylases and epoxidases . For example, a hydroxylase that can hydroxylate oleic acid (18:1) and is linked to PC to form ricinoleic acid (18:1-OH) has been described in castor plants . In addition, a bifunctional enzyme that can both catalyze hydroxylation and desaturate the same fatty acid has been reported in a Brassicacae species . Moreover, it has been demonstrated that a plant desaturase can be converted into a hydroxylase, and vice versa, by swapping a few specific amino acids by targeted mutagenesis . Since hydroxylase, desaturase, and epoxidase activity could lead to the production of hydroxy acids (Figure 6A), we expressed the functional protein encoded by PP1002E07 in yeast to clarify the reaction catalyzed by the putative encoded enzyme. Yeasts expressing the peach ORF gene PP1002E07 accumulated linoleic acid (18:2) when oleic acid (18:1) levels were reduced (Figure 5), indicating that the candidate gene has ω-6 Oleate desaturase activity. No ricinoleic acid was produced in the yeast expressing the candidate gene, suggesting that the encoded protein is not a bifunctional enzyme. Since no enzymes in plants have been shown to simultaneously desaturate and epoxidate fatty acid , our results suggest that the candidate gene identified herein is a monofunctional ω-6 Oleate desaturase.
Linoleic acid (18:2) may be further desaturated to linolenic acid (18:3), and both compounds could enter the Lipoxygenase (LOX) pathway to be catabolized into C6 and C9 alcohols, which are the substrates of AAT enzymes in the synthesis of a variety of esters (Figure 6A, ). Moreover, double bonds of unsaturated fatty acids could undergo the introduction of oxygen by hydratation to form hydroxy acids, which could be biosynthetic precursors of lactones (Figure 6A). In a previous report, we suggested that lipid-derived compounds and lactones are inversely regulated in peach fruit . An analysis of the current data set also revealed lactones, specifically those in C2 and C3, showing strong inverse correlations with lipid-derived compounds belonging to C11 and C12 (Figure 4A). This is in agreement with previous results indicating that lactone content increases during maturity and ripening, while lipid-derived compounds decrease [35–38]. The fact that ω-6 Oleate desaturase positively correlates with lactones, i.e., its levels increase during ripening, while lipid-derived compounds decrease, seems to suggest that the formation of lipid-derived compounds downstream of linoleic acid (18:2) is regulated during peach ripening, and that this pathway can feed into the precursors for lactone synthesis.
Beside CYPs and desaturase-like epoxidases, a third kind of enzyme, Peroxygenase (PXG), which can catalyze the epoxygenation of fatty acids, has been reported in plants. PXGs share no structural similarity with either peroxidase or cytochrome P450, but instead do so with caleosin, a small oil body-associated protein with both heme- and calcium-binding motifs . Here, a caleosin-related protein (PP1004C02), showing a strong correlation with g-jasmolactone, was indentified (Table 1, Additional file 10: Table S6). Recently, a PXG capable of epoxidating free oleic acid (18:1) was described in oat seeds . It is interesting to note that PXG has potential kinase phosphorylation sites , and our approach also identified kinases that correlate well with C3 (PPN010B11, PPN069F09).
PXG, like FAD and CYPs, is found in plant microsomes, which implies that it could be associated with the endoplasmatic reticulum in vivo [50, 59]. Interestingly, two SNARE genes associated with the endomembrane trafficking system (PP1004F06 and PPN031F12) also correlated with the C3 volatile cluster (Table 1, Additional file 10: Table S6), indicating the co-activation of both structural and metabolic enzymes.
It has been proposed that highly correlated metabolites in replicated experiments reflect the existence of a common regulatory mechanism in the system under study . In this sense, our approach identifies several transcription factors relating to volatile clusters C1 (PP1009A02, PPN066C05) and C2 (PPN010C01), which are tentatively involved in the regulation of volatile production during peach fruit ripening.
The highly stringent criteria used in our approach indicate that our candidate genes are unlikely to be the result of random correlations. In any case, further functional studies are required to unequivocally prove the implication of the candidate gene in controlling the proposed volatile accumulation.
Selection of candidate genes for other aroma-related compounds
Candidate genes for other minor aroma volatiles in peach have also been discovered. Ethyl acetate has been described to have a “fruity” aroma  and has also been proposed to be regulated independently of other esters in peach . Here we demonstrate that the HCA identifies an enolase (phosphopyruvate hydratase, PPN071H07), which is highly related to the accumulation of both Ethanol and its acetate ester, Ethyl acetate (Additional file 9: Figure S4, Table 1). Enolase catalyzes one of the final glycolysis steps (EC. 22.214.171.124), which generates Pyr as an end product (Figure 6B). As previously mentioned, Pyr is also converted into acetaldehyde by the action of Pyr decarboxylase . From acetaldehyde, the pathway splits into two branches, one producing ethanol and the other acetyl-CoA (Figure 6B, ). Ethyl acetate can be formed by the condensation of ethanol with Acetyl-CoA by a putative AAT, which is yet to be described. The correlation between ethanol and ethyl acetate with the enolase suggests that this gene is a good candidate for driving the production of glycolysis end products, i.e., Pyr, into aroma-contributing compounds.
The aroma of Benzeneacetaldehyde is described as “green” (Additional file 1: Figure S1) and could, therefore, confer immature fruit notes. Accordingly, the levels of this compound decreased during the ripening of the ‘Granada’ genotype (Figure 1) and after the shelf-life ripening in ‘MxR_01’ (Figure 2), which is in accordance with the shift of aroma from immature to ripe fruit. The correlation network analysis reveals that putative tyrosine (Tyr) aminotransferase (PPN018G03) correlates well (r=0.85) with this volatile, which is, in turn, related to other phenolic volatiles (Benzaldehyde and Benzyl chloride, Figure 1, Figure 4). Two pathways for the biosynthesis of aromatic volatiles have been described in different fruits (Figure 6C). In tomato, Phenylalanine (Phe) undergoes two sequential enzymatic steps: a decarboxylation step followed by a deamination step to form Benzeneacetaldehyde . In contrast, in melon, Phe first loses the amine groups by transamination reactions to then undergo a decarboxylation reaction . Aromatic acid transaminase (CmArAT1), which is involved in Benzeneacetaldehyde biosynthesis in melon, is capable of catalyzing the transamination of either Tyr or Phe . The gene identified herein (PPN018G03) is phylogenetically related to CmArAT1 (data not shown), suggesting that it may have both transaminase activities. Therefore, we hypothesize that Benzeneacetaldehyde in peach is produced by a similar transaminase pathway to that in melon (Figure 6C).
As previously mentioned, the catabolism of linoleic and linolenic acid by the LOX/HPL pathway produces the so-called lipid-derived compounds . As some of these volatiles are described to confer “green” aromas (Additional file 1: Figure S1), which is the typical aroma of freshly-cut grass, these volatiles have been traditionally associated with unripe fruits. Accordingly, we previously found that lipid-derived volatiles correlated highly among themselves and also with parameters that measure the peach ripening stage, indicating that high levels are linked to unripe fruit . In addition, some of these compounds have been described to have other aromas, which could even be unpleasant when present in fruit, these being “fermented”, “spicy”, “fatty”, or even “chemical” (Additional file 1: Figure S1). In the data set used herein, these compounds are grouped into two highly interrelated clusters (C11 and C12, Figure 1, Figure 4A). The correlation network analysis (Additional file 13: Figure S5) reveals that a transcription factor (PPN019F11) strongly correlates with the expression levels of a gene with no homolog in Arabidopsis (PPN023E05) which, in turn, is highly correlated to a lipid δ 9 desaturase (PPN052H12). This lipid desaturase strongly and inversely correlates with lipid-derived compounds Furan 2-pentyl- and Hexanal. It is possible that the lipid desaturase catalyzes the formation of unsaturated lipids, which are not substrates for the LOX/HPL pathway, thus driving the carbon allocation toward other products. If this were the case, this gene would be a candidate for reducing lipid-derived compound production. Four candidate genes for C12 volatiles were selected because they correlated with 2,4-Heptadienal (E,E)- in the HCA (Additional file 9: Figure S4). Fatty acids are stored as triacylglycerides, so acyl hydrolases should free them to enter the LOX/HPL pathway . One candidate gene is an esterase homolog (PPN033D11) that may be involved in the supply of fatty acids to the LOX/HPL pathway that produces lipid-derived volatiles, which makes it an interesting target for decreasing undesired aromas in peach fruit.
Cluster C4 is formed by two volatiles deriving from the terpenoid biosynthesis pathway and by four volatiles of unknown origin (Figure 1, Additional file 1: Figure S1). Our approach was able to identify the genes associated with C4 volatiles, although their functional annotations do not suggest a mechanism (Table 1, Additional file 10: Table S6). Nonetheless, it is quite likely that hormones are involved in the regulation of these compounds since gibberellin- and auxin-responsive proteins were identified (PP1000F05 and PP1009D02, respectively), as was aminopeptidase P (PPN017B08; Table 1, Additional file 10: Table S6). In Arabidopsis, it has been demonstrated that auxin signaling can be blocked by auxin-transport inhibitors, and that the specific Aminopeptidase P regulates hormone perception by degrading these inhibitors .
Our “omics” approach could be improved by including additional samples, i.e., fruits from other genotypes or those subjected to other treatments. For example, it has been demonstrated that low temperature storage and hot water treatments lower the volatile content of peach, while other quality parameters, e.g., flesh firmness, remain unaffected [46, 65]. Thus, if the volatile reduction is due to changes in gene expression, these treatments may be used to unlink both processes and therefore improve gene identification by co-regulation patterns. Although further experiments are required to unequivocally prove some of the links between the genes and volatile compounds found herein, we provide empirical data, which may prove useful to boost gene-function discovery studies in peach fruit.
In the present study, a hypothesis-free approach based on correlation analyses was applied to discover genes that are highly associated with the main aroma-contributing compounds of peach fruit. Consequently, a set of candidate genes for lactones, esters, terpenoids, phenolics and lipid-derived volatiles was identified. Some have functional annotations, which clearly suggest a possible role in controlling the production of related volatiles, which qualifies them as strong candidates to be selected for functional studies. In addition, for one of the candidate genes, we demonstrate that the encoded protein possesses enzymatic activity that produces a volatile peach precursor. Although our main interest was related to the practical implications of our results, the data provided herein not only support a previous hypothesis related to the physiology and metabolism of peach fruit, but also suggest novel ones. Although more studies are required to unequivocally prove the associations found herein, we provide experimental evidence and data analyses that support the role of a number of candidate genes in the control of volatile compounds, some of which may prove useful as targets for peach improvement and/or biotechnology industry interests.
We are grateful to Cristina Besada, PhD (Instituto Valenciano de Investigaciones Agrarias, IVIA, Spain) for her help with the fruit quality parameter analyses. We are also thankful to Cristina Martí and Clara Pons (Instituto de Biología Molecular y Celular de Plantas, IBMCP, Spain) for their advice on microarray analyses. Jesús Garcia Brunton, PhD for providing the fruits used in this study (Instituto Murciano de Investigación y Desarrollo Agrario, IMIDA, Spain). HS-SPME-GC-MS analyses were performed at the Metabolomic lab facilities at the IBMCP (CSIC) in Spain. GS has financial support from INTA (Instituto Nacional de Tecnología Agropecuaria, Argentine). This project has been funded by the Ministry of Economy and Competitivity grant AGL2010-20595.
- Arús P, Verde I, Sosinski B, Zhebentyayeva T, Abbott A: The peach genome. Tree Genetics & Genomes. 2012, 8 (3): 531-547.View ArticleGoogle Scholar
- Verde I, Abbott AG, Scalabrin S, Jung S, Shu S, Marroni F, Zhebentyayeva T, Dettori MT, Grimwood J, Cattonaro F, et al: The high-quality draft genome of peach (Prunus persica) identifies unique patterns of genetic diversity, domestication and genome evolution. Nat Genet. 2013, 45 (5): 487-494.View ArticlePubMedGoogle Scholar
- Ahmad R, Parfitt D, Fass J, Ogundiwin E, Dhingra A, Gradziel T, Lin D, Joshi N, Martinez-Garcia P, Crisosto C: Whole genome sequencing of peach (Prunus persica L.) for SNP identification and selection. BMC Genomics. 2011, 12 (1): 569-PubMed CentralView ArticlePubMedGoogle Scholar
- Verde I, Bassil N, Scalabrin S, Gilmore B, Lawley CT, Gasic K, Micheletti D, Rosyara UR, Cattonaro F, Vendramin E, et al: Development and Evaluation of a 9K SNP Array for Peach by Internationally Coordinated SNP Detection and Validation in Breeding Germplasm. PLoS One. 2012, 7 (4): e35668-PubMed CentralView ArticlePubMedGoogle Scholar
- Ogundiwin EA, Marti C, Forment J, Pons C, Granell A, Gradziel TM, Peace CP, Crisosto CH: Development of ChillPeach genomic tools and identification of cold-responsive genes in peach fruit. Plant Mol Biol. 2008, 68 (4–5): 379-397.View ArticlePubMedGoogle Scholar
- ESTree Consortium: Development of an oligo-based microarray (uPEACH 1.0) for genomics studies in peach fruit. Acta Hortic. 2005, 682: 263-268.View ArticleGoogle Scholar
- Spolaore S, Trainotti L, Casadoro G: A simple protocol for transient gene expression in ripe fleshy fruit mediated by Agrobacterium. J Exp Bot. 2001, 52 (357): 845-850.PubMedGoogle Scholar
- FAO: 2010, Accessed 16 Nov. 2012. http://faostat.fao.org,
- Bruhn CM, Feldman N, Garlitz C, Harwood J, Ivans E, Marshall M, Riley A, Thurber D, Williamson E: Consumer perceptions of quality: Apricots, cantaloupes, peaches, pears, strawberries, and tomatoes. Journal of Food Quality. 1991, 14 (3): 187-195.View ArticleGoogle Scholar
- Parker DD, Zilberman D, Moulton K: How quality relates to price in California fresh peaches. Calif Agr. 1991, 45 (2): 14-16.Google Scholar
- Sánchez G, Besada C, Badenes ML, Monforte AJ, Granell A: A non-targeted approach unravels the volatile network in peach fruit. PLoS One. 2012, 7 (6): e38992-PubMed CentralView ArticlePubMedGoogle Scholar
- Eduardo I, Chietera G, Bassi D, Rossini L, Vecchietti A: Identification of key odor volatile compounds in the essential oil of nine peach accessions. J Sci Food Agric. 2010, 90 (7): 1146-1154.View ArticlePubMedGoogle Scholar
- Dufosse L, Latrasse A, Spinnler H-E: Importance of lactones in food flavours: structure, distribution, sensory properties and biosynthesis. Sci Aliment. 1994, 14: 17-50. (in french)Google Scholar
- Gatfield IL: Biotechnological production of natural flavor materials. 1999, New York: Kluwer Academic / Plenum Publishing Corporation, 211-227.Google Scholar
- Schwab W, Davidovich-Rikanati R, Lewinsohn E: Biosynthesis of plant-derived flavor compounds. Plant J. 2008, 54 (4): 712-732.View ArticlePubMedGoogle Scholar
- Schöttler M, Boland W: Biosynthesis of Dodecano-4-lactone in ripening fruits: crucial role of an epoxide-hydrolase in enantioselective generation of aroma components of the nectarine (Prunus persica var. nucipersica) and the strawberry (Fragaria ananassa). Helvetica Chimica Acta. 1996, 79 (5): 1488-1496.View ArticleGoogle Scholar
- Vecchietti A, Lazzari B, Ortugno C, Bianchi F, Malinverni R, Caprera A, Mignani I, Pozzi C: Comparative analysis of expressed sequence tags from tissues in ripening stages of peach (Prunus persica L. Batsch). Tree Genetics & Genomes. 2009, 5: 377-391.View ArticleGoogle Scholar
- Brandi F, Bar E, Mourgues F, Horvath G, Turcsi E, Giuliano G, Liverani A, Tartarini S, Lewinsohn E, Rosati C: Study of ‘Redhaven’ peach and its white-fleshed mutant suggests a key role of CCD4 carotenoid dioxygenase in carotenoid and norisoprenoid volatile metabolism. BMC Plant Biol. 2011, 11 (1): 24-PubMed CentralView ArticlePubMedGoogle Scholar
- Eduardo I, Chietera G, Pirona R, Pacheco I, Troggio M, Banchi E, Bassi D, Rossini L, Vecchietti A, Pozzi C: Genetic dissection of aroma volatile compounds from the essential oil of peach fruit: QTL analysis and identification of candidate genes using dense SNP maps. Tree Genetics & Genomes. 2013, 9 (1): 189-204.View ArticleGoogle Scholar
- Saito K, Matsuda F: Metabolomics for functional genomics, systems biology, and biotechnology. Annu Rev Plant Biol. 2010, 61 (1): 463-489.View ArticlePubMedGoogle Scholar
- Carrera J, Rambla JL, Pons C, Jaramillo A, Elena SF, Granell A: Fine-tuning tomato agronomic properties by computational genome redesign. PLoS Comput Biol. 2012, 8 (6): e1002528-PubMed CentralView ArticlePubMedGoogle Scholar
- Mounet F, Moing A, Garcia V, Petit J, Maucourt M, Deborde C, Bernillon S, Le Gall G, Colquhoun I, Defernez M: Gene and metabolite regulatory network analysis of early developing fruit tissues highlights new candidate genes for the control of tomato fruit composition and development. Plant Physiol. 2009, 149 (3): 1505-1528.PubMed CentralView ArticlePubMedGoogle Scholar
- Lombardo VA, Osorio S, Borsani J, Lauxmann MA, Bustamante CA, Budde CO, Andreo CS, Lara MV, Fernie AR, Drincovich MF: Metabolic profiling during peach fruit development and ripening reveals the metabolic networks that underpin each developmental stage. Plant Physiol. 2011, 157 (4): 1696-1710.PubMed CentralView ArticlePubMedGoogle Scholar
- Amiour N, Imbaud S, Clément G, Agier N, Zivy M, Valot B, Balliau T, Armengaud P, Quilleré I, Cañas R, et al: The use of metabolomics integrated with transcriptomic and proteomic studies for identifying key steps involved in the control of nitrogen metabolism in crops such as maize. J Exp Bot. 2012, 63 (14): 5017-5033.View ArticlePubMedGoogle Scholar
- Carrari F, Baxter C, Usadel B, Urbanczyk-Wochniak E, Zanor MI, Nunes-Nesi A, Nikiforova V, Centero D, Ratzka A, Pauly M: Integrated analysis of metabolite and transcript levels reveals the metabolic shifts that underlie tomato fruit development and highlight regulatory aspects of metabolic network behavior. Plant Physiol. 2006, 142 (4): 1380-1396.PubMed CentralView ArticlePubMedGoogle Scholar
- Casati P, Campi M, Morrow D, Fernandes J, Walbot V: Transcriptomic, proteomic and metabolomic analysis of UV-B signaling in maize. BMC Genomics. 2011, 12 (1): 321-PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Y, Yang C, Li S, Yang L, Wang Y, Zhao J, Jiang Q: Volatile characteristics of 50 peaches and nectarines evaluated by HP-SPME with GC-MS. Food Chem. 2009, 116 (1): 356-364.View ArticleGoogle Scholar
- Tikunov Y, Lommen A, de Vos CHR, Verhoeven HA, Bino RJ, Hall RD, Bovy AG: A novel approach for nontargeted data analysis for metabolomics. Large-scale profiling of tomato fruit volatiles. Plant Physiol. 2005, 139 (3): 1125-1137.PubMed CentralView ArticlePubMedGoogle Scholar
- Meisel L, Fonseca B, Gonzalez S, Baeza-Yates R, Cambiazo V, Campos R, Gonzalez M, Orellana A, Retamales J, Silva H: A rapid and efficient method for purifying high quality total RNA from peaches (Prunus persica) for functional genomics analyses. Biol Res. 2005, 38 (1): 83-88.View ArticlePubMedGoogle Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci. 2001, 98 (9): 5116-5121.PubMed CentralView ArticlePubMedGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13 (11): 2498-2504.PubMed CentralView ArticlePubMedGoogle Scholar
- Assenov Y, Ramírez F, Schelhorn S-E, Lengauer T, Albrecht M: Computing topological parameters of biological networks. Bioinformatics. 2008, 24 (2): 282-284.View ArticlePubMedGoogle Scholar
- Venegas-Calerón M, Beaudoin F, Garcés R, Napier JA, Martínez-Force E: The sunflower plastidial ω3-fatty acid desaturase (HaFAD7) contains the signalling determinants required for targeting to, and retention in, the endoplasmic reticulum membrane in yeast but requires co-expressed ferredoxin for activity. Phytochemistry. 2010, 71 (10): 1050-1058.View ArticlePubMedGoogle Scholar
- Derail C, Hofmann T, Schieberle P: Differences in key odorants of handmade juice of yellow-flesh peaches (Prunus persica L.) induced by the workup procedure. J Agric Food Chem. 1999, 47 (11): 4742-4745.View ArticlePubMedGoogle Scholar
- Chapman GW, Horvat RJ, Forbus WR: Physical and chemical changes during the maturation of peaches (cv. Majestic). J Agric Food Chem. 1991, 39 (5): 867-870.View ArticleGoogle Scholar
- Visai C, Vanoli M: Volatile compound production during growth and ripening of peaches and nectarines. Sci Hortic. 1997, 70 (1): 15-24.View ArticleGoogle Scholar
- Zhang B, Shen J-y, Wei W-w, Xi W-p, Xu C-J, Ferguson I, Chen K: Expression of genes associated with aroma formation derived from the fatty acid pathway during peach fruit ripening. J Agric Food Chem. 2010, 58 (10): 6157-6165.View ArticlePubMedGoogle Scholar
- Aubert C, Gunata Z, Ambid C, Baumes R: Changes in physicochemical characteristics and volatile constituents of yellow- and white-fleshed nectarines during maturation and artificial ripening. J Agric Food Chem. 2003, 51 (10): 3083-3091.View ArticlePubMedGoogle Scholar
- Trainotti L, Bonghi C, Ziliotto F, Zanin D, Rasori A, Casadoro G, Ramina A, Tonutti P: The use of microarray μPEACH1.0 to investigate transcriptome changes during transition from pre-climacteric to climacteric phase in peach fruit. Plant Sci. 2006, 170 (3): 606-613.View ArticleGoogle Scholar
- 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 (12): 3299-3308.View ArticlePubMedGoogle Scholar
- Manganaris GA, Vasilakakis M, Diamantidis G, Mignani I: Diverse metabolism of cell wall components of melting and non-melting peach genotypes during ripening after harvest or cold storage. J Sci Food Agric. 2006, 86 (2): 243-250.View ArticleGoogle Scholar
- Hirai MY, Sugiyama K, Sawada Y, Tohge T, Obayashi T, Suzuki A, Araki R, Sakurai N, Suzuki H, Aoki K, et al: Omics-based identification of Arabidopsis Myb transcription factors regulating aliphatic glucosinolate biosynthesis. Proc Natl Acad Sci. 2007, 104 (15): 6478-6483.PubMed CentralView ArticlePubMedGoogle Scholar
- Hirai MY, Klein M, Fujikawa Y, Yano M, Goodenowe DB, Yamazaki Y, Kanaya S, Nakamura Y, Kitayama M, Suzuki H, et al: Elucidation of gene-to-gene and metabolite-to-gene networks in arabidopsis by integration of metabolomics and transcriptomics. J Biol Chem. 2005, 280 (27): 25590-25595.View ArticlePubMedGoogle Scholar
- Giovannoni J: Molecular biology of fruit maturation and ripening. Annu Rev Plant Biol. 2001, 52 (1): 725-749.View ArticleGoogle Scholar
- Klee HJ, Giovannoni JJ: Genetics and control of tomato fruit ripening and quality attributes. Annu Rev Genet. 2011, 45 (1): 41-59.View ArticlePubMedGoogle Scholar
- Zhang B, Xi W-p, Wei W-w, Shen J-y, Ferguson I, Chen K-s: Changes in aroma-related volatiles and gene expression during low temperature storage and subsequent shelf-life of peach fruit. Postharvest Biology and Technology. 2011, 60 (1): 7-16.View ArticleGoogle Scholar
- D’ Auria JC: Acyltransferases in plants: a good time to be BAHD. Curr Opin Plant Biol. 2006, 9 (3): 331-340.View ArticleGoogle Scholar
- Pesis E: The role of the anaerobic metabolites, acetaldehyde and ethanol, in fruit ripening, enhancement of fruit quality and fruit deterioration. Postharvest Biology and Technology. 2005, 37 (1): 1-19.View ArticleGoogle Scholar
- Sauveplane V, Kandel S, Kastner P-E, Ehlting J, Compagnon V, Werck-Reichhart D, Pinot F: Arabidopsis thaliana CYP77A4 is the first cytochrome P450 able to catalyze the epoxidation of free fatty acids in plants. FEBS J. 2009, 276 (3): 719-735.View ArticlePubMedGoogle Scholar
- Cahoon EB, Ripp KG, Hall SE, McGonigle B: Transgenic production of epoxy fatty acids by expression of a cytochrome p450 enzyme from euphorbia lagascaeseed. Plant Physiol. 2002, 128 (2): 615-624.PubMed CentralView ArticlePubMedGoogle Scholar
- Morant M, Jørgensen K, Schaller H, Pinot F, Møller BL, Werck-Reichhart D, Bak S: CYP703 is an ancient cytochrome p450 in land plants catalyzing in-chain hydroxylation of lauric acid to provide building blocks for sporopollenin synthesis in pollen. The Plant Cell Online. 2007, 19 (5): 1473-1487.View ArticleGoogle Scholar
- Nelson D, Werck-Reichhart D: A P450-centric view of plant evolution. Plant J. 2011, 66 (1): 194-211.View ArticlePubMedGoogle Scholar
- Shanklin J, Cahoon EB: Desaturation and related modifications of fatty acids. Annu Rev Plant Physiol Plant Mol Biol. 1998, 49 (1): 611-641.View ArticlePubMedGoogle Scholar
- van de Loo FJ, Broun P, Turner S, Somerville C: An oleate 12-hydroxylase from Ricinus communis L. is a fatty acyl desaturase homolog. Proc Natl Acad Sci. 1995, 92 (15): 6743-6747.PubMed CentralView ArticlePubMedGoogle Scholar
- Broun P, Boddupalli S, Somerville C: A bifunctional oleate 12-hydroxylase: desaturase from Lesquerella fendleri. Plant J. 1998, 13 (2): 201-210.View ArticlePubMedGoogle Scholar
- Broun P, Shanklin J, Whittle E, Somerville C: Catalytic plasticity of fatty acid modification enzymes underlying chemical diversity of plant lipids. Science. 1998, 282 (5392): 1315-1317.View ArticlePubMedGoogle Scholar
- Hanano A, Burcklen M, Flenet M, Ivancich A, Louwagie M, Garin J, Blée E: Plant seed peroxygenase is an original heme-oxygenase with an ef-hand calcium binding motif. J Biol Chem. 2006, 281 (44): 33140-33151.View ArticlePubMedGoogle Scholar
- Meesapyodsuk D, Qiu X: A peroxygenase pathway involved in the biosynthesis of epoxy fatty acids in oat. Plant Physiol. 2011, 157 (1): 454-463.PubMed CentralView ArticlePubMedGoogle Scholar
- Murphy DJ, Hernández-Pinzón I, Patel K: Role of lipid bodies and lipid-body proteins in seeds and other tissues. J Plant Physiol. 2001, 158 (4): 471-478.View ArticleGoogle Scholar
- Camacho D, de la Fuente A, Mendes P: The origin of correlations in metabolomics data. Metabol. 2005, 1: 53-View ArticleGoogle Scholar
- Guillot S, Peytavi L, Bureau S, Boulanger R, Lepoutre J-P, Crouzet J, Schorr-Galindo S: Aroma characterization of various apricot varieties using headspace-solid phase microextraction combined with gas chromatography–mass spectrometry and gas chromatography-olfactometry. Food Chem. 2006, 96 (1): 147-155.View ArticleGoogle Scholar
- Baldwin EA, Scott JW, Shewmaker CK, Schuch W: Flavor trivia and tomato aroma: biochemistry and possible mechanisms for control of important aroma components. HortScience. 2000, 35 (6): 1013-1022.Google Scholar
- Gonda I, Bar E, Portnoy V, Lev S, Burger J, Schaffer AA, Tadmor Y, Gepstein S, Giovannoni JJ, Katzir N, et al: Branched-chain and aromatic amino acid catabolism into aroma volatiles in Cucumis melo L. fruit. J Exp Bot. 2010, 61 (4): 1111-1123.PubMed CentralView ArticlePubMedGoogle Scholar
- Murphy AS, Hoogner KR, Peer WA, Taiz L: Identification, purification, and molecular cloning of n-1-naphthylphthalmic acid-binding plasma membrane-associated aminopeptidases from Arabidopsis. Plant Physiol. 2002, 128 (3): 935-950.PubMed CentralView ArticlePubMedGoogle Scholar
- Raffo A, Nardo N, Tabilio M, Paoletti F: Effects of cold storage on aroma compounds of white- and yellow-fleshed peaches. Eur Food Res Technol. 2008, 226 (6): 1503-1512.View ArticleGoogle Scholar