Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

A Multidisciplinary Approach Providing New Insight into Fruit Flesh Browning Physiology in Apple (Malus x domestica Borkh.)

  • Mario Di Guardo ,

    Contributed equally to this work with: Mario Di Guardo, Alice Tadiello

    Affiliation Genomics and Crop Biology Department Research and Innovation Centre, Fondazione Edmund, Mach, San Michele all’Adige (Trento), Italy

  • Alice Tadiello ,

    Contributed equally to this work with: Mario Di Guardo, Alice Tadiello

    Affiliation Genomics and Crop Biology Department Research and Innovation Centre, Fondazione Edmund, Mach, San Michele all’Adige (Trento), Italy

  • Brian Farneti,

    Affiliation Department of Fruit Trees & Woody Plant Science, University of Bologna, Bologna, Italy

  • Giorgia Lorenz,

    Affiliation Genomics and Crop Biology Department Research and Innovation Centre, Fondazione Edmund, Mach, San Michele all’Adige (Trento), Italy

  • Domenico Masuero,

    Affiliation Genomics and Crop Biology Department Research and Innovation Centre, Fondazione Edmund, Mach, San Michele all’Adige (Trento), Italy

  • Urska Vrhovsek,

    Affiliation Genomics and Crop Biology Department Research and Innovation Centre, Fondazione Edmund, Mach, San Michele all’Adige (Trento), Italy

  • Guglielmo Costa,

    Affiliation Department of Fruit Trees & Woody Plant Science, University of Bologna, Bologna, Italy

  • Riccardo Velasco,

    Affiliation Genomics and Crop Biology Department Research and Innovation Centre, Fondazione Edmund, Mach, San Michele all’Adige (Trento), Italy

  • Fabrizio Costa

    fabrizio.costa@fmach.it

    Affiliation Genomics and Crop Biology Department Research and Innovation Centre, Fondazione Edmund, Mach, San Michele all’Adige (Trento), Italy

Correction

6 Jan 2014: Di Guardo M, Tadiello A, Farneti B, Lorenz G, Masuero D, et al. (2014) Correction: A Multidisciplinary Approach Providing New Insight into Fruit Flesh Browning Physiology in Apple (Malus x domestica Borkh.). PLOS ONE 9(1): 10.1371/annotation/72dfe3df-4cfd-49a8-93db-979bdd8794c3. https://doi.org/10.1371/annotation/72dfe3df-4cfd-49a8-93db-979bdd8794c3 View correction

Abstract

In terms of the quality of minimally processed fruit, flesh browning is fundamentally important in the development of an aesthetically unpleasant appearance, with consequent off-flavours. The development of browning depends on the enzymatic action of the polyphenol oxidase (PPO). In the ‘Golden Delicious’ apple genome ten PPO genes were initially identified and located on three main chromosomes (2, 5 and 10). Of these genes, one element in particular, here called Md-PPO, located on chromosome 10, was further investigated and genetically mapped in two apple progenies (‘Fuji x Pink Lady’ and ‘Golden Delicious x Braeburn’). Both linkage maps, made up of 481 and 608 markers respectively, were then employed to find QTL regions associated with fruit flesh browning, allowing the detection of 25 QTLs related to several browning parameters. These were distributed over six linkage groups with LOD values spanning from 3.08 to 4.99 and showed a rate of phenotypic variance from 26.1 to 38.6%. Anchoring of these intervals to the apple genome led to the identification of several genes involved in polyphenol synthesis and cell wall metabolism. Finally, the expression profile of two specific candidate genes, up and downstream of the polyphenolic pathway, namely phenylalanine ammonia lyase (PAL) and polyphenol oxidase (PPO), provided insight into flesh browning physiology. Md-PPO was further analyzed and two haplotypes were characterised and associated with fruit flesh browning in apple.

Introduction

Fruit quality features, represented by the biochemical and physical properties making a fruit edible and appreciated by consumers, are nowadays considered a major priority in several breeding programmes worldwide. For this reason, in the last decade the scientific community has initiated a series of extensive studies addressed at examining the genetic and molecular mechanisms responsible for controlling the fundamental physiological processes ultimately leading to fruit quality [1-5]. In the modern system of fruit distribution and marketing, fruit quality also needs to be guaranteed during storage, allowing high quality standards to be maintained. In addition to this, pre-prepared fresh cut fruit is rapidly becoming interesting, given increased consumer demand for fresh and convenient food [6,7]. However, the fruit processing procedure can encounter serious problems, which need to be prevented in order to ensure high quality and at the same time avoid substantial loss of fruit. One of the most important problems occurring during fruit processing is flesh browning, which is undesirable due to the aesthetically unpleasant appearance and consequential off-flavour [8,9].

The apple’s susceptibility to flesh browning is thought to be the result of complex interplay between the polyphenoloxidase (PPO) enzyme and polyphenol content [10]. PPO is a bi-copper metalloenzyme showing two conserved copper-binding domains, CuA and CuB, which interact with molecular oxygen and phenolic substrates [11]. The PPO enzyme catalyses the formation of quinones from phenols through the hydroxylation of monophenols to o-diphenols (such as chlorogenic acid and (-)-epicatechin) and their subsequent dehydrogenation to o-quinones [12]. O-quinones are highly reactive, as they can undergo self-polymerisation or react with amines and thiol groups to form brown pigments (Figure 1). In the cell, after a PPO is encoded in the nucleus it is transported to the cytoplasm, where proPPO formed is stored in the plastid. The enzyme is thus physically separated from its substrate, which is located in the vacuole [12-14]. The loss of subcellular structural compartmentation, following damage or cutting, leads the PPO enzyme to come into contact with phenolic compounds. The quinones produced by oxidation rapidly condense to generate relatively insoluble brown polymers (melanins), ultimately detected as the browning reaction. For the conversion of o-dihydroxyphenols into o-benzoquinones, PPO can also use O2 as a second substrate (catecholase activity; [13]).

thumbnail
Figure 1. Diagram of the biosynthetic pathway for polyphenol compounds, highlighting the cascade of compounds produced from phenylalanine and the action of the PPO enzyme.

https://doi.org/10.1371/journal.pone.0078004.g001

Initial browning can also occur during storage in a controlled atmosphere, causing serious fruit loss. The importance of browning prevention depends on the fact that the symptoms in whole fruit are not usually externally visible, and thus the problem can only be detected after purchase, influencing consumer confidence. Of several strategies adopted to prevent fruit flesh browning, the application of low O2, together with a high concentration of CO2 are the most common, although an excessive level of CO2 may also affect fruit quality negatively [15]. For quality control in minimally processed fruit, another system for reducing flesh browning is control of PPO activity. This enzyme can indeed be inactivated by ascorbic acid or its derivates, which reduce phenoxyl radicals and quinone forms of phenolics back to precursor phenols, through a coupled oxidation/reduction reaction [16,17]. However, the effect of ascorbic acid is temporary, as it is further irreversibly oxidized into dehydroascorbic acid, enabling browning to occur again [18].

Due to this mechanism, development of functional markers for the selection of new apple accessions characterised by a minimum browning phenotype would be a valuable approach, also taking into account the fact that several chemicals (such as sulphate) and physical prevention strategies are no longer accepted by the FDA agency [15,19]. To date, genetic control of flesh browning has only been partially assessed using the QTL mapping approach, revealing the presence of two intervals located on linkage groups 3 and 17 [20,21]. However, these two QTLs have not yet been shown to be associated with any candidate genes involved in the flesh browning metabolism. In the apple, the only PPO gene positioned to date on a genetic map is represented by [22], but its genomic location on chromosome 5 has not yet been associated with any browning phenotype.

In this survey, two independent apple populations were employed to examine and characterise fruit flesh browning in apple. A new gene, here called Md-PPO, was identified and associated with the different browning rate. To better understand this physiological mechanism, both transcript and metabolite assessment was carried out on a set of apple samples developing browning. Finally, a new hypothesis regarding the physiological regulation of flesh browning is discussed here.

Materials and Methods

Plant materials

To investigate fruit flesh browning in apple and detect QTL regions involved in control of this phenotype, two independent full-sib populations were employed. The first progeny came from the controlled cross between ‘Fuji’ and ‘Pink Lady’ (POP_1), while the second (POP_2) was created by crossing the apple cv. ‘Golden Delicious’ with ‘Braeburn’. Both populations were planted in the same experimental orchard in 2003 and both belong to the breeding programme currently underway at FEM (Foundation Edmund Mach), located in northern Italy. All the plant material was grafted onto M9 rootstock and maintained by making use of standard technical management for pruning and disease control.

Fruit Flesh Browning Phenotyping

Fruit flesh browning in apple was assessed in the two populations at two specific stages in the fruit ripening process. Apples from POP_1 were assessed at harvest, while for POP_2 flesh browning was also measured after two months’ cold storage, to estimate the impact of low temperature in the control of this disorder. For convenience in terms of the genotyping procedure, both progenies were represented by the two parental cultivars, plus 94 seedlings each (to set a 96 well plate). The browning phenotype was assessed by measuring the flesh colour from two halves of a cut apple with a Konica Minolta digital colorimeter CM2600D. Five whole apples/genotype were considered, collecting ten measurements/samples at the end. The use of this instrument allowed digital characterisation of the flesh colour, based on the tristimulus coordinates L*, a* and b*. This colour space was chosen with regards to the most popular RGB (red, green and blue) system because it provides a more uniform colour difference in relation to human perception. For this reason the L* a* b* is the colour space generally employed in photoelectric measurements. The L* value spans from 0 to 100, corresponding to black (0) and white (100), while a* and b* indexes do not have specific numerical limits. Negative numbers for a* indicate a green shade while positive values point to magenta. For b*, negative and positive values are for blue and yellow respectively. Analysis of flesh browning was performed by measuring the flesh colour at three specific moments: after cutting (T0), after 30 minutes (T30) and after 60 minutes (T60) on the same sample. Furthermore, the absolute value for L*, a* and b* taken at these three stages (T0, T30 and T60) and the delta values (∆) between T30 and T0, T60 and T0 and T60 and T30 were considered. The percentage variation, represented by ((Ti - T0)/T0)x100, where Ti is the L*, a* or b* value after 30 or 60 minutes and T0 is the corresponding index measured after cutting, was also calculated. Finally, apple flesh browning was examined using twenty-four new colour parameters.

PPO identification

To target the PPO gene set located over the apple genome assembly, BLAST searching was performed, using the only available and functionally well-characterised PPO gene for apple as the query, coded in the GeneBank (www.ncbi.nlm.nih.gov) as L29450 [12]. BLAST was performed using the internal resources of FEM (http://genomics.research.iasma.it) and the predicted gene set produced during the apple genome sequencing project [23]. The sequences retrieved were confirmed through gene annotation, performed using the Uniprot database (http://uniprot.org), and further compared with the one released by GDR (http://rosaceae.org). Gene distribution over the chromosomes was illustrated using MapChart [24].

Marker development

Two types of functional molecular markers related to the PPO elements identified here were newly designed: SSRs and SNPs. Initially a set of SSR motifs were discovered in silico in the genomic contigs of each PPO, using Sputnik software (http://espressosoftware.com/sputnik/index.html). Primer pairs (Table S1) flanking the microsatellite motif, were designed by Primer3 software (http://frodo.wi.mit.edu/primer3/) and further used for specific amplification. SSR markers were initially exploited on the four parental cultivars to search for polymorphism in a PCR mix with a final volume of 20 μL. The PCR conditions were as follow: 5 ng of DNA, 10 X Buffer, 0.25mM dNTPs, 0.075 µM of forward labelled and reverse primers and 1000 U of 5Prime® Taq polymerase. Temperature conditions were: initial denaturation at 94 °C for 150 s followed by 32 cycles at 94°C for 30 s, 58°C for 45 s, 72°C for 60 s, and a final extension at 72 °C for 5 min. Fragment analysis was performed by an ABI PRISM® 3730 capillary sequencer (Applied Biosystems by Life Technologies, Carlsbad, CA, USA) in a final volume of 0.3 µL of PCR product, 9.67 µL formamide and 0.03 µL of 500-LIZ denaturated for 3 min at 95 °C. Fragment sizing was carried out with PeakScanner® software (Applied Biosystems by Life Technologies, Carlsbad, CA, USA).

In addition to the SSR markers, a set of SNPs related to the PPO gene corresponding to L29450 was also created. Three pairs of primers were designed, in order to cover the full length of the entire gene: PPO_1for: CTTCTTGGTCTTGGAGGTCT and PPO_1rev: ATCGGAGCTTGTCGTAGAGA; PPO_2for: CCACAACTCATGGCTCTTCT and PPO_2rev: CTAACTCTGCTGTCTCGTTG; PPO_3for: GTTCTTTGGGAACCCGTACA and PPO_3rev: CATCAAACTTCACAGCCACG. Amplification of the three fragments was carried out in a final volume of 10 µL containing 5 ng of DNA, 10X buffer, 0.25 mM dNTPs, 1 µM forward and reverse primers and 1000U 5Prime® Taq polymerase. PCR was carried out using an ABI 2720 Thermal Cycler (Applied Biosystems by Life Technologies, Carlsbad, CA, USA) with the following thermal conditions: denaturation at 94 °C for 2 minutes, followed by 35 cycles of denaturation at 94 °C for 45 s, annealing at 65 °C for 45 s, extension at 72 °C for 2 min, finishing with a conclusive extension at 72 °C for 10 min. Amplicons were than purified by adding 4 µl of water and 1µl of ExoSAP to 1 µL of PCR product. The mix was then incubated for 45 minutes at 37 °C, followed by 15 minutes at 75 °C. In a second step, 2 µL of Buffer 5X, 1 µL of the forward primer (at 3.2 uM), 1 µL of Big dye terminator (Applied Biosystem by Life Technologies, Carlsbad, CA, USA) and 0.5 µL of water were added to the purification mix to reach a final volume of 10 µL. The sequencing reaction was performed for 2 minutes at 96 °C followed by 39 cycles of denaturation (96 °C for 10 s), annealing (50 °C for 5 s) and extension (60 °C for 4 minutes). The sequencing runs were carried out using an ABI PRISM® 3730 capillary sequencer (Applied Biosystems by Life Technologies, Carlsbad, CA, USA) and finally the SNPs were genotyped by re-sequencing the amplicons across the two populations. SNP scoring was carried out by analysing the sequences obtained using Pregap4 software (Staden package; http://staden.sourceforge.net).

QTL analysis

The new markers designed in this study were used to improve the current version of the two molecular maps. The POP_1 map was originally created for comprehensive examination of fruit texture in apple [25], while POP_2 was employed to anchor the contigs produced during the apple genome sequencing projects [23]. Marker integration in the context of these two maps was done using JoinMap 4.0 software [26], employing the Kosambi mapping function, while visual display of the linkage groups was carried out using MapChart software [24]. Finally, genetic and phenotypic data of 47 seedlings of the POP_1 were used for QTL analysis, in order to target the genomic regions associated with apple flesh browning. QTL intervals were detected using MapQTL®6 [27], and the Interval Mapping (IM) algorithm. To reduce residual variance and the effect of possible false positives, the markers with the highest LOD value were selected as cofactors in the subsequent MQM computation. The threshold for calling positive QTLs was established at a LOD value of 3, after running 1000 permutations.

To examine the genetic regulation of this phenotype, each QTL interval was anchored on the assembled apple genome and the sequences of the predicted gene were retrieved by using the Computational Web Resources of FEM (www.genomics.research.iasma.it). From the geneset underlying the QTL regions, the predicted aminoacid sequences were employed to perform gene annotation by interrogating the UniProt (Viridiplantae) database.

Finally, the effect of the Md-PPO haplotype was validated on two set of seedlings belonging to the POP_2, 64 measured at harvest and 58 after two months of cold storage.

RNA isolation and qPCR analysis

For transcript profiling, flesh samples from four replicates of the ‘Golden Delicious’ collected at T0 (immediately after cutting), T30 and T60 were used. Cortex tissues were initially frozen in liquid nitrogen and further ground to a fine powder. RNA extraction was performed using the Spectrum Plant total RNA kit (Sigma). RNA was quantified using a NanoDrop ND-8000 spectrophotometer (Thermo Scientific, USA), while its purity and integrity was assessed by an Agilent 2100 Bioanalyzer. The RNA isolated from the apple samples (T0, T30 and T60) was converted into cDNA by the ‘SuperScript VILO cDNA Synthesis Kit’ (Invitrogen). Prior to this, 1 µg of total RNA from each sample was pretreated with 2 Units of rDNAse I (DNA free kit, Ambion) and used as a starting template. To clarify the physiological regulation of flesh browning, the expression of two genes, Md-PAL and Md-PPO, was assessed by using the following specific primer pairs: Md-PALfor: AGACCCTCAATGCCTCAGAA; Md-PALrev: CAAGCCAGAACCAACAGCAG; Md-PPOfor: CCTACTCACAAAGCCCAAGC and Md-PPOrev: CCTCCAAGACCAAGAAGCAC.

Real Time PCR was performed using ACTIN as housekeeping gene (Md-ACTfor: TGACCGAATGAGCAAGGAAATTACT; Md-ACTrev:TACTCAGCTTTGGCAATCCACATC).

Transcript quantification, carried out using the ViiA7 instrument (Applied Biosystem), was performed by using the ‘KAPA SYBR FAST Universal qPCR kit’ (Kapa Biosystems). PCR thermal conditions were: incubation at 95°C for 20 sec, 40 cycles of 95°C 1 sec and 60°C 20 sec. Finally, a cycle at 95°C for 15 sec, 60°C for 1 min and 95°C for 15 sec was applied to determine the melting curve. The Ct results were obtained by averaging three independent normalised expression values for each sample, carried out using Q-gene software [28]. Relative gene expression was plotted as the mean of the normalised expression values of the triplicates.

Phenolic profiling

Phenols were extracted from the ground cortex tissues of ‘Golden Delicious’ collected at T0, T30 and T60, following the procedure reported in [29] and analysed as described in [30]. Briefly, 2 g of tissue powder collected from four replicates each experimental time, previously prepared for RNA isolation, were treated in sealed glass vials using 4 ml of water/methanol/chloroform solution (20:40:40). After vortexing for 1 min, the samples were mixed using an orbital shaker for 15 min at room temperature, and further centrifuged at 1000g (4 °C) for 10 min, after which the upper phases, made up of aqueous methanol extract, were collected. Extraction was repeated by adding another 2.4 ml of water/methanol (1:2) to the pellet and chloroform fractions. After the final centrifugation, the upper phases from the two extractions were combined and brought to the volume of 10 ml and filtered with a 0.2 μm PTFE filter prior to liquid chromatography-mass spectrometry analysis. Ultraperformance liquid chromatography was performed employing a Waters Acquity UPLC system (Milford, MA, USA) coupled to a Waters Xevo TQMS (Milford, MA, USA) working in ESI ionisation mode. Separation of the phenolic compounds was achieved on a Waters Acquity HSS T3 column 1.8 μm, 100 mm × 2.1 mm (Milford, MA, USA), kept at 40 °C, and with two solvents: A (water containing 0.1% formic acid) and B (acetonitrile containing 0.1% formic acid). The samples were eluted according to a linear gradient method described in detail by Vrhovsek et al. [30]. 2 μL of the final extraction were injected by an autosampler set at the temperature of 6 °C. Data were processed by using Waters MassLynx 4.1 and TargetLynx software. The compounds analysed belong to four chemical classes: hydroxycinnamic acids, dihydrochalcones, flavan-3-ols and flavonols.

Results and Discussion

Distribution of apple flesh browning traits

Flesh browning in apple is understood as the colour change in the fruit cortex after cutting. Analysis was performed immediately after cutting and after 30 and 60 minutes. Within 1 hour a significant variation in colour intensity was observed, as shown by the four parental cultivars (Figure S1 i-iv). Flesh browning was assessed in two populations (POP_1 and POP_2), and the distribution of the several parameters considered to examine the phenotype is illustrated in Figure 2. This shows the data distribution of the L*, a* and b* index absolute values, measured after 60 minutes of exposure to the air, in order to evaluate the maximum trait variability existing within each progeny. In both cases, all the browning sub-phenotypes showed quantitative distribution, which is the main prerequisite for a QTL mapping survey. However, some exceptions with skewed segregation were detected, such as L* at T60. It is also worth noting that analysis of data distribution revealed a transgressive type of segregation (particularly in the POP_1). A set of seedlings exceeding the phenotypic value of the two parents was indeed observed, making the efforts to identify molecular markers suitable for the anticipated selection of this trait worthwhile. This type of segregation is also reflected in the difference of the L*, a* and b* values measured over the three stages for the four parental cultivars (Figure S1_v).

thumbnail
Figure 2. Trait distribution for L*, a* and b* parameters measured after 60 minutes of exposure to air.

The panels on the left are for POP_1 (‘Fuji x Pink Lady’), while the panels on the right are for POP_2 (‘Golden Delicious x Braeburn’). The grey bars in each panel represent the data distribution observed, while the black line is the normal fit of data. The x-axis gives the value for the three indexes, while the y-axis plots the number of observations.

The position of the four parental varieties is indicated by four letters: F_’Fuji’, PL_’Pink Lady’, GD_’Golden Delicious’ and B_’Braeburn’.

https://doi.org/10.1371/journal.pone.0078004.g002

PPO Organisation in the Apple Genome

To identify the several PPO genes positioned within the apple genome assembly, the L29450 sequence was used as query, as to date it is the only PPO gene whose regulation has been experimentally validated for apple [12]. The predicted apple gene set was investigated using iterative nucleotide BLAST searching, and only genes showing an e-value of ≤ 2e-04 were further considered as possible polyphenol oxidase candidates. The resulting gene list was further annotated using BLASTp (adopting the amino acid sequences resulting from the predicted gene set) through interrogation of the Uniprot database. Of the initial set, ten genes were finally annotated as PPOs and located in three distinct chromosomes (Figure 3). In particular, one element was located at the bottom of chromosome 2 (MDP0000500159), five were clustered at the top of chromosome 5 (MDP0000609966, MDP0000222503, MDP0000207799, MDP0000173059 and MDP0000221498) while the last group, made up of four elements, was detected at the bottom of chromosome 10 (MDP0000744636, MDP0000234782, MDP0000709073 and MDP0000699845; available at GDR: http://www.rosaceae.org/gb/gbrowse/malus_x_domestica/ and FEM: http://genomics.research.iasma.it/gb2/gbrowse/apple). The number of PPOs identified for apple is also consistent with what has been reported so far for other species. In tomato, for instance, seven genes belonging to this family were positioned on chromosome 8 [31]. In this species, as observed for apple, the PPO elements are organised in clusters. With regards to tomato, apple shows three chromosomes characterized by the presence of PPO elements, but chromosome 2 contains only one element, and chromosomes 5 and 10 are homoeologous, due to the recent genome duplication occurring in this species [23]. Moreover, other authors have already presented the duplication of these two chromosomes in apple in previous publications [32-35]. Recently, a comprehensive survey performed on 25 sequenced genomes revealed that other species shared a similar number of PPO genes [36], such as Sorghum bicolour (8), Glycine max (11) and Populus trichocarpa (11). Beside this, it was also discovered that the PPO family can vary greatly in size across species. Indeed, Solanum tuberosum shows only 5 PPO elements [37], Zea mays 6 and Vitis vinifera 4. This difference was thought to be a consequence of duplication events, responsible, for instance, for the complete loss of this gene family in Arabidopsis thaliana [36,38]. The organisation of the cluster of PPO genes in apple is also supported by phylogenetic analysis carried out by Tran and colleagues, which suggests the tandem arrangement of PPO genes on chromosomes as a consequence of these recent genome duplication events. Of the total number of PPO genes found in the apple genome, the one encoded as MDP0000699845 was shown to be the most similar to L29450, initially used as the query, and for this reason it was further studied in more detail, to exploit its association with the fruit flesh browning rate.

thumbnail
Figure 3. Physical localization of the ten PPO genes discovered in the ‘Golden Delicious’ genome.

The position in Kb is plotted for each chromosome on the left , while the gene ID is shown on the right.

https://doi.org/10.1371/journal.pone.0078004.g003

Functional marker design and QTL mapping

The first set of functional markers designed was related to microsatellite motifs discovered on the contigs where the ten PPO genes were targeted. Initially, these SSR markers were tested for polymorphism, searching in the four parental cultivars (‘Fuji’, ‘Pink Lady’, ‘Golden Delicious’ and ‘Braeburn’). Of these markers, only two primer pairs provided interesting results: Md-PPO_SSR_ch5e and Md-PPO_SSR_ch10d (Figures S2 and S3). The position of these two markers, assigned in silico to chromosomes 5 and 10 (as reported by the Genome Browser), was further confirmed by genetic mapping of the two populations considered in this investigation (Figure S4 and S5). Md-PPO_SSR_ch05e was mapped at the top of LG (linkage group) 5 in both progenies (at 17.4 cM in POP_1 and 4.8 cM in POP_2). In the same way, Md-PPO_SSR_ch10d was also mapped at 64.4 cM and 80 cM from the top of linkage group 10 for POP_1 and POP_2 respectively. The different genetic positions observed could be attributed to different genome coverage of these two linkage maps, being higher in POP_2.

Because SNPs found within gene sequences are considered to be one of the most frequent and important causal events controlling phenotype variation [39-41], a set of SNP markers was also exploited. It is indeed known that substitution of a single base can modify the aminoacid sequence, leading to functional variation [42]. Furthermore, given that they are the most abundant type of molecular marker within the genome, SNP markers are currently recognised as those most valuable in finding associations with phenotypic traits [43]. As proof of this, several references have already reported the use of SNPs in trait association analysis [41,44-51]. SNP discovery was performed by selecting the gene MDP0000699845 as the main candidate, on which three primer pairs were designed to characterise the full-length sequence of the gene. On aligning the sequences read in the four parental cultivars, two SNPs were exploitable in POP_2 only, being homozygous in the parents of POP_1. The two SNPs, located respectively at 170 bp and 500 bp, were genotyped by sequencing all the seedlings in the ‘Golden Delicious x Braeburn’ population, showing 1:1 segregation (χ2 = 0.4). Both SNPs were heterozygous in the apple cv. ‘Braeburn’ (SNP170: CG and SNP500: AG), while in ‘Golden Delicious’ the allelotype configuration was always homozygous (GG) for both SNPs. The close proximity between these two SNPs generated identical segregation (high LD level), which enabled construction of a haplotype, further positioned on the POP_2 genetic map (Figure S4). As expected the Md-PPO gene was mapped by means of this haplotype on LG 10 at 1.7 cM from Md-PPO_SSR_ch10d. This distance may however be overestimated, due to the fact that the SNP set used to saturate the maps came from the sequencing of the ‘Golden Delicious’ genome [52]. Because of this, the SNPs were fully informative only for the ‘Golden Delicious’ parent (“ab x aa” and “ab x ab”), while in ‘Braeburn’ the segregation type “aa x ab” was completely absent, leading to reduced representation of recombination events for this paternal cultivar in the POP_2 map. Between the two SNPs identified by comparing ‘Golden Delicious’ and ‘Braeburn’, SNP170 caused an amino acid substitution in the protein sequence (from Glycine in ‘Golden Delicious’ to Glutamic Acid in ‘Braeburn’; Figure S6) that may theoretically contribute towards modifying the occurrence of flesh browning observed in this progeny.

Integration of these markers led to an improved version of the two maps, which were subsequently used for a QTL mapping survey carried out on POP_1 (composed by 481 markers for a total length of 1430.8 cM with an average distance between markers of 2.9 cM) and a further haplotype validation step performed on POP_2 (composed by 608 markers for a total length of 1204 cM with an average distance between markers of 1.9 cM). Within the POP_1 genome, six linkage groups out of seventeen (apple chromosome numbers), showed significant intervals for the presence of twenty-five QTLs associated with fruit flesh browning. Five of these were identified when a* and b* absolute values were measured at T0, T30 and T60. In this case it is interesting to note that at the three measuring times not a single QTL was assigned the L* absolute value parameter, which is related to colour brightness (black/white gradient). The distribution of the a* value led instead to identification of a QTL located on LG 16 (LOD: 3.85) when measured after cutting, and on LG 11 (LOD: 4.12) and 14 (LOD: 3.92) at the two stages (T30 and T60) of exposure to air (Figure 4 and Table 1). The change in LGs during the time-course was also detected for the b* colorimetric parameter. Indeed, after cutting the QTL was identified on LG 9 (LOD: 4.0), while after 30 and 60 minutes another set of QTLs was instead targeted on LG 11 (LOD: 3.21), 13 (LOD: 3.93) and 14 (LOD: 3.39). LOD profile comparison of the three times for each colorimetric value revealed a shift in the set of chromosomes associated with the QTLs from the first assessment (immediately after cutting), as compared to fruit exposed to air. It is worth noting that the QTLs identified after cutting (T0), are mainly located on chromosomes 9 and 16, not being detected at any other location over the genome. These two QTLs should thus be more associated with flesh colour properties than flesh browning, hypothesis supported by the recent discovery of the gene set controlling red flesh in the apple and located respectively on these two chromosomes, such as MdMYB10 [53,54] and MdLAR [55,56]. As browning is thought to be the intrinsic capacity of a fruit to change flesh colour after wounding, a delta (∆) value, calculated between two moments, would possibly be more effective in finding genomic regions involved in the control of browning development. When ∆ values (as well as the variation expressed as the percentage variation between T30-T0 and T60-T0) for the three colorimetric indexes (L*, a* and b*) were considered in the computation, two linkage groups, LG 10 and LG 11, were mainly shown to be involved, spanning from a minimum LOD value of 3.08 (26.1 % of expressed variance) to a maximum of 4.68 (36.8%; Table 1). The last group of targeted QTLs was instead located on chromosome 14, and was associated with the ∆ value for a* and b* between T60 and T30, the absolute value for b* at T30 and T60 and for a* at T60, with LOD values from 3.39 (28.3%) to 4.99 (38.6%). Of the chromosomes identified in this survey considered to be relevant in the control of the flesh browning, chromosomes 10, 11 and 14 are the most significant, due to the presence of several QTLs in clusters, associated with different colour parameters, whereas the other three (9, 13 and 16) are represented by only one QTL each. The QTL-LOD profile defined for LG 11 also suggested the presence of two intervals associated with fruit flesh browning, one located at 25 cM and one below, at approximately 50 cM from the top of the linkage group. It is moreover worth noting that interesting QTLs associated with the concentration of Ascorbic Acid (AsA) and its oxidized form Dehydroascorbate (DHA), playing an important role in control of fruit flesh browning, were recently mapped [57] on chromosomes 10 (containing the Md-PPO gene) and 11. Indeed, in the apple the concentration of DHA was correlated to susceptibility to flesh browning [20], while in pears the tendency to develop internal browning was linked to a decreased concentration of AsA [58]. This finding, together with the results discussed here, suggests that chromosome 10 is the best candidate for regulation of the flesh browning phenotype, due to the simultaneous presence of genes responsible for its occurrence as well as its prevention.

thumbnail
Figure 4. QTL-LOD profile targeted on 6 linkage groups of the POP_1 and related to several flesh browning components.

In particular:.

LG9: b*_T0 (solid red).

LG10: L*T30-0 (solid green), a*_T30-0 (solid violet), L*_T60-0 (solid dark green), a*_T60-0 (dashed black), L*T30-0% (dashed blue), L*_T60-0% (dashed pink).

LG11: a*_T30 (solid dark green), b*_T30 (solid red), a*_T60 (solid blue), L*_T30-0 (solid green), a*_T30-0 (solid violet), b*_T30-0 (solid light blue), a*_T60-0 (dashed black), L*_T30-0% (dashed blue), a*_T30-0% (dashed yellow), L*_T60-0% (dashed pink), a*_T60-0% (dashed light green).

LG13: b*_T30 (solid blue).

LG14: b*_T30 (solid black), a*_T60 (solid yellow), b*_T60 (solid pink), L*_T60-30 (dashed red), a*_T60-30 (dashed green).

LG16: a*_T0 (solid black).

The QTL analysis was performed by using the phenotypic data of 47 seedlings of the POP_1. For each seedling a total of five apples were considered (biological replicates), and for each fruit the browning was measured on the two cut halves (technical replicates).

https://doi.org/10.1371/journal.pone.0078004.g004

TraitLGLODMarker%Expmu_acmu_admu_bcmu_bd
a*_T0163.85GDsnp0004731.41.2246.971.422.43
b*_T094GDsnp0080332.410.7110.5112.3210.39
a*_T30114.12GDsnp0191233.22.594.123.213.36
b*_T30113.21GDsnp002252713.7415.9113.7714.47
133.93CH01B1231.916.5413.6814.1414.29
a*_T60113.52GDsnp0041829.22.944.793.563.8
143.39GDsnp0006228.33.673.095.033.69
b*_T60143.96GDsnp0188832.215.3214.0417.1315.29
L_T30-T0103.96GDsnp0200932.1-3.78-2.21-2.42-3.4
114.68GDsnp0024036.8-2.18-4.01-2.17-2.88
a_T30-T0103.08GDsnp0218326.12.221.281.522.32
114.17GDsnp0024033.61.222.411.541.73
b_T30-T0113.42GDsnp0003728.52.474.233.343.42
L_T60-T0104.04GDsnp0218332.7-4.78-2.67-3.29-4.19
a_T60-T0103.69GDsnp0218330.32.831.591.952.69
114.39GDsnp00240351.572.941.872.24
L_T60-T30144.04GDsnp0016532.7-0.06-0.75-1.72-0.47
a_T60-T30144.99GDsnp0016538.60.350.330.910.32
L_T30-T0%104.2GDsnp0200933.7-5.23-2.91-3.22-4.57
a_T30-T0%113.39GDsnp0003728.395.06201.69124.7119.44
L_T60-T0%103.93GDsnp0218332-6.47-3.62-4.4-5.58
a_T60-T0%113.39GDsnp0024028.3119.88242.25147.89151.07

Table 1. Summary of QTLs associated with flesh browning identified in POP_1.

The trait, linkage group (LG), maximum LOD value identified within the interval, closest marker to the LOD peak, percentage of expressed phenotypic variance (% Exp) and estimated mean of distribution for the four possible genotypes (mu_ac, mu_ad, mu_bc and mu_bd) are given for each QTL .
CSV
Download CSV

Finally, for better understanding the genetic control of flesh browning, the genomic intervals identified in this work were anchored on the genome assembly of ‘Golden Delicious’, in order to allow in silico annotation. Of the several genes identified over the six targeted chromosomes (Table S2), it is worth highlighting 52 elements (Table 2) mainly involved in cell wall metabolism and secondary metabolites (polyphenols [60-62]). In the first category, three gene families known to participate in cell wall disassembly were retrieved, namely pectinesterase, pectin lyase and polygalacturonase. Furthermore, the fifteen genes were mainly located on chromosomes 10 and 16, which are known to be hot spot regions for control of fruit texture physiology [25]. In particular, the presence of a polygalacturonase gene responsible for the texture variability observed in apple [48,51] has recently been validated on chromosome 10. The collocation of genes encoding cell wall degrading enzymes and QTLs discovered for fruit flesh browning has already been described in peach [59], where two browning QTLs co-mapped with candidates annotated as pectate lyase and expansin genes. The simultaneous presence of these two categories in the QTL intervals discovered for flesh browning can be explained by the fact that disassembly of internal cellular structure compartmentation is needed for flesh browning to occur. In this scenario, we can surmise the involvement of this class of genes to promote this process, facilitating the interplay between the polyphenol oxidase (stored in the plastid) and its phenolic substrate (stored in the vacuole).

Gene familyNumberLG_9LG_10LG_11LG_13LG_14LG_16
Glycosyltransferase10314110
pectinesterase9200403
Cellulose synthase5040001
Flavonoid o-methyltransferase5050000
Polygalacturonase5010103
Caffeic acid 3-O-methyltransferase5000050
Chalcone synthase2200000
Glucosyltransferase2002000
Leucoanthocyanidin dioxygenase2000200
flavonoid 3-o-glucosyltransferase1010000
Isoflavone reductase1000100
Trans-cinnamate 4-monooxygenase1000100
Polyphenol oxidase1000100
Cellulase 1000001
Isoflavonoid regulator1000001
Pectin lyase1000001

Table 2. Characterisation of 52 elements belonging to 16 categories mainly involved in cell wall and polyphenol metabolism.

The relative number and chromosome locations are given for each class. This list is a specific selection of the complete QTL annotation given in Table S2.
CSV
Download CSV

Finally, another PPO element was newly discovered on chromosome 13, not previously targeted during in silico genome analysis (because of its poor sequence similarity with L29450), increasing the final number of PPO genes present in the apple genome to eleven.

Md-PPO haplotype validation

The effect of the Md-PPO haplotype (defined by the two SNPs targeted within the MDP0000699845 full-length) was further validated in POP_2 (‘Golden Delicious’ x ‘Braeburn’) assessed at two specific ripening stages, at harvest and after two months’ cold storage, in order to evaluate the effect of low storage temperature on flesh browning development. Marker validation was performed by grouping the seedlings into two classes according to their haplotypes (given as “np” and “nn” in Figure 5), which showed significant differences for ∆ a* at T30-T0 and T60-T0 based on the LSD-ANOVA test (Pvalue ≤ 0.05). For the two ripening stages (harvest and after two months’ cold storage) both parameters were statistically significant, with the heterozygous haplotype associated with a reduced variation in flesh colour, thus less prone to developing fruit flesh browning after cutting. These results also suggest that the colour change observed during flesh browning can be mostly attributed to a variation in the a* parameter. This data is also in agreement with observations of significant changes in this colour index between T0 and T30 for the four apple cultivars, suggesting that fruit browning develops rapidly in the first 30 minutes (Figure S1). It is also interesting to note the differences observed during postharvest storage (2 months’ cold storage). The association observed between the rate of flesh browning and the Md-PPO haplotype was indeed confirmed after storage, but with a reduction of about 27% and 10% for the “np” and “nn” classes of seedling respectively. A reduced polyphenol content during cold storage [63] may reduce substrate availability for PPOs, resulting in a lower flesh browning rate. The difference observed between harvest and two months after storage was also confirmed by the P-values. At harvest, in fact, the value computed between the two seedling classes was more significant (P-value: 0.04) than the one calculated after storage (P-value: 0.05).

thumbnail
Figure 5. Md-PPO haplotype validation in the POP_2 progeny.

White bars represent individuals characterised by a heterozygous haplotype (“np”), while the grey bars show the homozygous category (“nn”). The standard error is reported in each bar.

On the x-axis, 30-0 and 60-0 are the percentage variation of ∆ a* calculated between T30-T0 and T60-T0 respectively (2). indicates the ∆ value calculated after two months’ cold storage.

As for the analysis carried out for the POP_1, five apples (biological replicates) were assessed for each genotypes, and for each fruit the flesh browning was measured of the two halves (technical replicates).

Asterisks show statistically significant comparison based on the LSD-ANOVA test (P-value ≤ 0.05).

https://doi.org/10.1371/journal.pone.0078004.g005

Expression profiling and polyphenolic characterisation during the development of fruit flesh browning in apple

To explain the physiological regulation leading to the occurrence of fruit flesh browning in more detail, two main genes were considered, Md-PAL (phenylalanine ammonia lyase) and Md-PPO. The first gene, Md-PAL, is responsible for the biotransformation of L-phenyalanine to ammonia and trans-cinnamic acid, which is the first step in the biosynthetic pathway of polyphenol compounds [64,65]. The second, Md-PPO, is the same element investigated here and genetically associated with the flesh browning observed in the two populations. The expression profile of these two genes over the time-course showed distinct temporal activation, also in agreement with their respective positions along the polyphenol cascade (Figure 1). Md-PAL, which is located upstream in this pathway, was progressively expressed throughout the time-course, already showing increased mRNA accumulation after 30 minutes’ (T30) exposure to air (Figure 6a), consistent with the different accumulation of the main polyphenol classes. At T30 a general accumulation was indeed observed for hydroxycinnamic acid, dihydrochalcone and flavan-3-ol classes (Figure 6 c, d, e), while an almost unchanged situation was observed for flavonols (Figure 6f). It is worth noting that individual compounds showed the same trend as their respective phenolic groups, as indicated in Table S3. Increased production of these compounds by wounding has already been observed in lettuce [66], but never clearly examined in apple. This enhanced production, coordinated by the activation of Md-PAL expression, can be considered a response of the defence mechanism (antioxidant protection) as well as signalling [67,68]. In this scenario we can hypothesize that PPO enzyme is eventually synthesized to oxidize the amount of polyphenols produced after wounding, theory supported by the late expression observed for Md-PPO. This gene, located downstream in this pathway, indeed showed basal and consistent expression between T0 and T30, which was crucial in the development of flesh browning, while at T60 its transcript accumulation increased by about 5 fold (Figure 6b).

thumbnail
Figure 6. Expression profile of the two genes involved in the polyphenolic pathway, PAL (Md-PAL, panel a) and PPO (Md-PPO, panel b).

In both graphs the y axis shows the mean normalised expression, graphing the three samples selected to monitor flesh browning (T0, T30 and T60).

Over the same time-course, the accumulation of the four main polyphenolic compounds, namely hydroxycinnamic acid (c), dihydrochalcones (d), flavan-3-ols (e) and flavonols (f), is shown below in the figure. The amount of these compounds is plotted on the y axis and expressed as μg/g of fresh weight (FW). Polyphenolic profiling was performed analyzing four replicates/experimental time.

For each bar the standard error is also reported.

Samples significantly different are shown using different letters following a LSD-ANOVA test (P-value ≤ 0.05).

https://doi.org/10.1371/journal.pone.0078004.g006

These functional dynamics were also validated by metabolite screening of polyphenol characterisation. At T60 the four main polyphenol classes, described in Figure 6 (c, d, e, f), showed an average concentration decrease of about 1.5-2 fold. The browning colouration occurring in the first 30 minutes may be possibly caused by PPO enzymes already available and stored in the plastid. After this activation, the wounded organs stimulate higher production of polyphenols as a defence signal, which is then maintained by a feed-back controlling mechanism, through the production of additional PPO enzymes devoted to their oxidation. In this system the fruit can regulate the signalling triggered by the wounding events, while the flesh browning appearing immediately after cutting seems to be more related to the initial genetically determined amount of PPO enzyme, together with its substrates.

Conclusion

This work offers new insight, shedding light on the genetic regulation of apple flesh browning. This flesh browning, which seriously decreases the quality of the final fruit in minimally processed fruit, seems to be related to a genetically determined accumulation of the PPO enzyme, which can differ greatly in the various apple cultivars. This theory is supported by the fact that Md-PPO gene is highly activated only in the last phase of the time-course designed here, while browning occurs much earlier, corresponding to the increased transcription of another gene, Md-PAL, known to stimulate the entire polyphenolic cascade. Md-PPO was then thought to be activated more to regulate the polyphenolic signalling system. However, this gene showed two distinct haplotypes, which may be responsible for the initial enzymatic amount stored in the cell (before wounding), responsible for the extent of final browning.

The availability of a marker associated with this phenomenon may represent a valid alternative to destructive methods for the selection of low browning accession, suitable for improving the quality of minimally processed apples.

Supporting Information

Figure S1.

Fruit flesh browning evolution in the four parental apple cultivars, ‘Fuji’(i), ‘Pink Lady’ (ii), ‘Golden Delicious’ (iii) and ‘Braeburn’ (iv). For each variety the three panels are a_T0 (after cutting), b_T30 (after 30 minutes) and c_T60 (after 60 minutes). The histograms below each panel give the digital colour measurements obtained by the colorimeter and expressed as L*, a* and b*. Also for the four parental cultivars, five apples were assessed for each experimental time (T0, T30 and T60), representing the biological replicates, on which two colour measurements were performed on the two sides of a cut apple (technical replicates). The letters given show the statistical significance following to the LSD-ANOVA test (P-value ≤ 0.05). In slide “v” the difference between the two parental cultivars for L*, a* and b* values measured at T0, T30 and T60 is shown for each progeny. Statistically significant differences (P value ≤ 0.05) are highlighted with asterisks. The four parents are indicated as follows: ‘Fuji’ (green), ‘Pink Lady’ (pink), ‘Golden Delicious’ (blue) and ‘Braeburn’ (red). For each bar the standard error is also visualized. Below each histogram the actual value is reported.

https://doi.org/10.1371/journal.pone.0078004.s001

(PPT)

Figure S2.

Allelic polymorphism and segregation of the microsatellite marker MdPPO_SSR_ch5e. In the figure the first two rows refer to ‘Golden Delicious’ and ‘Braeburn’.

https://doi.org/10.1371/journal.pone.0078004.s002

(PPT)

Figure S3.

Allelic polymorphism and segregation of the microsatellite marker MdPPO_SSR_ch10d. In the figure the first two rows refer to ‘Braeburn’ and ‘Golden Delicious’.

https://doi.org/10.1371/journal.pone.0078004.s003

(PPT)

Figure S4.

POP_1 (‘Fuji x Pink Lady’) genetic map. The red text highlights the genetic position of the two SSR markers, respectively positioned on chromosome 5 (MdPPO_SSR_ch5e) and 10 (MdPPO_SSR_ch10d).

https://doi.org/10.1371/journal.pone.0078004.s004

(PPT)

Figure S5.

POP_2 (‘Golden Delicious x Braeburn’) genetic map. The red text highlights the genetic position of the two SSR markers, respectively positioned on chromosome 5 (MdPPO_SSR_ch5e) and 10 (MdPPO_SSR_ch10d).

https://doi.org/10.1371/journal.pone.0078004.s005

(PPT)

Figure S6.

Nucleotide and aminoacid sequence of Md-PPO. Section “a” of the figure shows the Md-PPO nucleotide sequence of ‘Golden Delicious’ (>PPO_GD; accession number: KF561111) and ‘Braeburn’ (>PPO_Br; accession number: KF561112). In the two sequences the two SNPs are highlighted in bold black text. Section “b” shows the alignment of the predicted aminoacids in ‘Golden Delicious’ (GD_prot) and Braeburn (Br_prot). The aminoacid substitution conferred by SNP170 is shown in bold black text.

https://doi.org/10.1371/journal.pone.0078004.s006

(DOC)

Table S1.

List of microsatellite marker primers designed on the contigs, related to the ten targeted PPO genes. The name, primer sequences, number of the contig and chromosome are given for each marker.

https://doi.org/10.1371/journal.pone.0078004.s007

(DOC)

Table S2.

Annotation of the gene underlying the QTL intervals identified on chromosomes 9, 10, 11, 13, 14 and 16.

https://doi.org/10.1371/journal.pone.0078004.s008

(XLS)

Table S3.

Trend for each individual phenolic compound. The mean concentration and respective standard deviation are given for each phenolic compound.

https://doi.org/10.1371/journal.pone.0078004.s009

(XLS)

Acknowledgments

We would like to thank Pierluigi Magnago and his team for field maintenance and the care of plant materials and Marco Fontanari and Giacomo Bazzanella for their help in sample collection.

Author Contributions

Conceived and designed the experiments: FC. Performed the experiments: MDG AT BF GL DM. Analyzed the data: MDG UV BF FC. Contributed reagents/materials/analysis tools: GC RV. Wrote the manuscript: FC MDG AT BF UV.

References

  1. 1. Giovannoni JJ (2001) Molecular biology of fruit maturation and ripening. Annu Rev Plant Physiol Plant Mol Biol 52: 725-749. doi:https://doi.org/10.1146/annurev.arplant.52.1.725. PubMed: 11337414.
  2. 2. Liebhard R, Kellerhals M, Pfammatter W, Jertmini M, Gessler C (2003) Mapping quantitative physiological traits in apple (Malus x domestica Borkh.). Plant Mol Biol 52: 511-526. doi:https://doi.org/10.1023/A:1024886500979. PubMed: 12956523.
  3. 3. Newcomb RD, Crowhurst RN, Gleave AP, Rikkerink EHA, Allan AC et al. (2006) Analyses of Expressed Sequence Tags from Apple. Plant Physiol 141: 147-166. doi:https://doi.org/10.1104/pp.105.076208. PubMed: 16531485.
  4. 4. Park S, Sugimoto N, Larson MD, Beaudry R, van Nocker S (2006) Identification of Genes with Potential Roles in Apple Fruit Development and Biochemistry through Large-Scale Statistical Analysis of Expressed Sequence Tags. Plant Physiol 141: 811-824. doi:https://doi.org/10.1104/pp.106.080994. PubMed: 16825339.
  5. 5. Costa F, Alba R, Schouten H, Soglio V, Gianfranceschi L et al. (2010) Use of homologous and heterologous gene expression profiling tools to characterize transcription dynamics during apple fruit maturation and ripening. BMC Plant Biol 10: 299. PubMed: 20973957.
  6. 6. Ohlsson T (1994) Minimal processing — preservation methods of the future — an overview. Trends Food Sci Technol 5: 341–344. doi:https://doi.org/10.1016/0924-2244(94)90210-0.
  7. 7. Gunes G, Watkins CB, Hotchkiss JH (2001) Physiological responses of fresh-cut apple slices under high CO2 and low O2 partial pressures. Postharvest Biol Technol 22: 197–204. doi:https://doi.org/10.1016/S0925-5214(01)00083-7.
  8. 8. Coseteng MY, Lee CY (1987) Changes in apple polyphenoloxidase and polyphenol concentrations in relation to degree of browning. J Food Sci 52(4): 985 - 989. doi:https://doi.org/10.1111/j.1365-2621.1987.tb14257.x.
  9. 9. Murata M, Tsurutani M, Tomita M, Homma S, Kaneko K (1995) Relationship between Apple Ripening and Browning: Changes in Polyphenol Content and Polyphenol Oxidase. J Agric Foods Chemestry 43(5).
  10. 10. Amiot MJ, Tacchini M, Aubert S, Nicolas J (1992) Phenolic compounds and browning susceptibility of various apple cultivars at maturity. J Food Sci 57(4): 958-962. doi:https://doi.org/10.1111/j.1365-2621.1992.tb14333.x.
  11. 11. Demeke T, Morris CF (2002) Molecular characterization of wheat polyphenol oxidase (PPO). Theor Appl Genet 104: 813-818. doi:https://doi.org/10.1007/s00122-001-0847-3. PubMed: 12582641.
  12. 12. Boss PK, Gardner RC, Janssen B-J, Ross GS (1995) An apple polypgenol oxidase cDNA is up-regulated in wounded tissues. Plant Mol Biol 27: 429-433. doi:https://doi.org/10.1007/BF00020197. PubMed: 7888632.
  13. 13. Martinez MV, Whitaker JR (1995) The biochemistry and control of enzymatic browning. Trends Food Sci Technol 6: 195-200. doi:https://doi.org/10.1016/S0924-2244(00)89054-8.
  14. 14. Murata M, Nishimura M, Murai N, Haruta M, Homma S et al. (2001) A transgenic apple callus showing reduced polyphenol oxidase activity and lower browning potential. Biosci Biotechnol Biochemestry 65(2): 383-388. doi:https://doi.org/10.1271/bbb.65.383.
  15. 15. Gorny JR, Hess-Pierce B, Cifuentes RA, Kader AA (2002) Quality changes in fresh-cut pear slices as affected by controlled atmospheres and chemical preservatives. Postharvest Biol Technol 24: 271-278. doi:https://doi.org/10.1016/S0925-5214(01)00139-9.
  16. 16. Nicolas JJ, Richard-Forget FG, Goupy PM, Amiot MJ, Aubert SY (1994) Enzymatic browning reactions in apple and apple products. Crit Rev Food Sci Nutr 34: 109-157. doi:https://doi.org/10.1080/10408399409527653. PubMed: 8011143.
  17. 17. Chow Y-N, Louarme L, Bonazzi C, Nicolas J, Billaud C (2011) Apple polyphenoloxidase inactivation during heating in the presence of ascorbic acid and chlorogenic acid. Foods Chemestry 129: 761-767.
  18. 18. Serpen A, Gökmen V (2007) Reversible degradation kinetics of ascorbic acid under reducing and oxidizing condition. Food Chem 104: 721-725. doi:https://doi.org/10.1016/j.foodchem.2006.11.073.
  19. 19. Buta JG, Moline HE, Spaulding DW, Wang CY (1999) Extending Storage Life of Fresh-Cut Apples Using Natural Products and Their Derivatives. J Agric Food Chem 1: 1–6. PubMed: 10563838.
  20. 20. Davey MW, Kenis K, Keulemans J (2006) Genetic control of fruit vitamin C contents. Plant Physiol 142: 343-351. doi:https://doi.org/10.1104/pp.106.083279. PubMed: 16844833.
  21. 21. Kenis K, Keulemans J, Davey MK (2008) Identification and stability of QTLs for fruit quality traits in apple. Tree Genet Genomes 4: 647-661. doi:https://doi.org/10.1007/s11295-008-0140-6.
  22. 22. Cova V, Perini D, Soglio V, Komjanc M, van de Weg WE et al. (2012) Exploiting expressed sequence tag databases for mapping markers associated with fruit development and fruit quality in apple. Mol Breeding 29: 699-715. doi:https://doi.org/10.1007/s11032-011-9584-8.
  23. 23. Velasco R, Zharkikh A, Affourtit J, Dhingra A, Cestaro A et al. (2010) The genome of the domesticated apple (Malus x domestica Borkh.). Nat Genet 42(10): 833-839. doi:https://doi.org/10.1038/ng.654. PubMed: 20802477.
  24. 24. Voorrips RE (2002) MapChart: Software for the graphical presentation of linkage maps and QTLs. J Hered 93: 77-78. doi:https://doi.org/10.1093/jhered/93.1.77. PubMed: 12011185.
  25. 25. Longhi S, Moretto M, Viola R, Velasco R, Costa F (2012) Comprehensive QTL mapping survey dissects the complex fruit texture physiology in apple (Malus x domestica Borkh.). J Exp Bot 63(3): 1107-1121. doi:https://doi.org/10.1093/jxb/err326. PubMed: 22121200.
  26. 26. Van Ooijen JW (2006) JoinMap® 4, Software for the Calculation of Genetic Linkage Maps in Experimental Populations. Wageningen, Netherlands: Kyazma, B.V.
  27. 27. Van Ooijen JW (2009) MAPQTL® 6, Software for the mapping of quantitative trait loci in experimental populations of diploid species. Wageningen, Netherlands: Kyazma, B.V.
  28. 28. Muller PY, Janovjak H, Miserez AR, Dobbie Z (2002) Processing of gene expression data generated by quantitative real-time RT-PCR. BioTechniques 32: 1372–1379. PubMed: 12074169.
  29. 29. Theodoridis G, Gika H, Franceschi P, Caputi L, Arapitsas P et al. (2012) LC-MS based global metabolite profiling of grapes: solvent extraction protocol optimization. Metabolomics 8: 175−185. doi:https://doi.org/10.1007/s11306-011-0389-x.
  30. 30. Vrhovsek U, Masuero D, Gasperotti M, Franceschi P, Caputi L et al. (2012) A Versatile Targeted Metabolomics Method for the Rapid Quantification of Multiple Classes of Phenolics in Fruits and Beverages. J Agric Food Chem 60: 8831−8840. doi:https://doi.org/10.1021/jf2051569. PubMed: 22468648.
  31. 31. Newman SM, Eannetta NT, Yu H, Prince JP, de Vicente MC et al. (1993) Organization of the tomato polyphenol oxidase gene family. Plant Mol Biol 21: 1035-1051. doi:https://doi.org/10.1007/BF00023601. PubMed: 8098228.
  32. 32. Liebhard R, Gianfranceschi L, Koller B, Ryder CD, Tarchini R et al. (2002) Development and characterisation of 140 new microsatellites in apple (Malus x domestica Borkh.). Mol Breeding 10: 217-241. doi:https://doi.org/10.1023/A:1020525906332.
  33. 33. Liebhard R, Koller B, Gianfranceschi L, Gessler C (2003) Creating a saturated reference map for the apple (Malus x domestica Borkh.) genome. Theor Appl Genet 106: 1497-1508. PubMed: 12677403.
  34. 34. Yamamoto T, Kimura T, Saito T, Kotobuki K, Matsuta N et al. (2004) Genetic linkage maps of Japanese and European pears aligned to the apple consensus map. Acta Hort 663: 51-56.
  35. 35. Folta KM, Gardiner SE (2009) Genetics and genomics of rosaceae. Springer Verlag, ISBN 978-0-387-77490-9. doi:https://doi.org/10.1007/978-0-387-77491-6.
  36. 36. Tran TL, Taylor JS, Constabel CP (2012) The polyphenol oxidase gene family in land plants: Lineage-specific duplication and expansion. BMC Genomics 13: 395. doi:https://doi.org/10.1186/1471-2164-13-395. PubMed: 22897796.
  37. 37. Thygesen PW, Dry IB, Robinson SP (1995) Polyphenol oxidase in potato. A multigene family that exhibits differential expression patterns. Plant Physiol 109: 525-531. doi:https://doi.org/10.1104/pp.109.2.525. PubMed: 7480344.
  38. 38. Van der Hoeven R, Ronning C, Giovannoni J, Martin G, Tanksley S (2002) Deduction about the number, organization and evolution of genes in the tomato genome based on analysis of a large expressed sequence tag collection and selective genomic sequencing. Plant Cell 14: 1441-1456. doi:https://doi.org/10.1105/tpc.010478. PubMed: 12119366.
  39. 39. Jorde LB (1995) Linkage disequilibrium as a gene-mapping tool. Am J Hum Genet 56: 11-14. PubMed: 7825565.
  40. 40. Jorde LB (2000) Linkage disequilibrium and the search for complex disease genes. Genome Res 10: 1435-1444. doi:https://doi.org/10.1101/gr.144500. PubMed: 11042143.
  41. 41. Rafalski A (2002) Application of single nucleotide polymorphisms in crop genetics. Curr Opin Plant Biol 5: 94–100. doi:https://doi.org/10.1016/S1369-5266(02)00240-6. PubMed: 11856602.
  42. 42. Krawczak M, Ball EV, Fenton I, Stenson PD, Abeysinghe S et al. (2000) Human gene mutation database-A biomedical information and research resource. Hum Mutat 15: 45-51. doi:https://doi.org/10.1002/(SICI)1098-1004(200001)15:1. PubMed: 10612821.
  43. 43. Schlötterer C (2004) The evolution of molecular markers – just a matter of fashion? Nat Rev Genet, 5: 63-69. PubMed: 14666112.
  44. 44. Fridman E, Pleban T, Zamir D (2000) A recombination hotspot delimits a wild-species quantitative trait locus for tomato sugar content to 484 bp within an invertase gene. Proc Natl Acad Sci USA 97(9): 4718-4723. doi:https://doi.org/10.1073/pnas.97.9.4718. PubMed: 10781077.
  45. 45. Gupta PK, Roy JK, Prasad M (2001) Single nucleotide polymorphisms: A new paradigm for molecular marker technology and DNA polymorphism detection with emphasis on their use in plants. Curr Sci 80(4): 524-535.
  46. 46. Wilson LM, Whitt SR, Ibanez AM, Rocheford TR, Goodman MM et al. (2004) Dissection of Maize Kernel Composition and Starch Production by Candidate Gene Association. The Plant Cell, doi/10.1105/tpc. p. 104.025700.
  47. 47. Lu Y, Yan J, Guimarães CT, Taba S, Hao Z et al. (2009) Molecular characterization of global maize breeding germplasm based on genome-wide single nucleotide polymorphisms. Theor Apple Genet 120: 93-115. doi:https://doi.org/10.1007/s00122-009-1162-7. PubMed: 19823800.
  48. 48. Costa F, Peace CP, Stella S, Serra S, Musacchi S et al. (2010) QTL dynamics for fruit firmness and softening around an ethylene-dependent polygalacturonase gene in apple (Malusxdomestica Borkh.). J Exp Bot 61(11): 3029-3039. doi:https://doi.org/10.1093/jxb/erq130. PubMed: 20462945.
  49. 49. Emanuelli F, Battilana J, Costantini L, Le Cunff et al Boursiquot JM. (2010) A candidate gene association study on muscat flavor in grapevine (Vitis vinifera. p. L.). BMC Plant Biology 10, 241.
  50. 50. Dunemann F, Ulrich D, Malysheva-Otto L, Weber WE, Longhi S et al. (2012) Functional allelic diversity of the apple alcohol acyl-transferase gene MdAAT1 associated with fruit ester volatile contents in apple cultivars. Mol Breeding 29: 609-625. doi:https://doi.org/10.1007/s11032-011-9577-7.
  51. 51. Longhi S, Hamblin MT, Trainotti L, Peace CP, Velasco R et al. (2013) A candidate gene based approach validates Md-PG1 as the main responsible for a QTL impacting fruit texture in apple (Malus x domestica Borkh.). BMC Plant Biol 13: 37. doi:https://doi.org/10.1186/1471-2229-13-37. PubMed: 23496960.
  52. 52. Micheletti D, Troggio M, Zharkikh A, Costa F, Malnoy M et al. (2011) Genetic diversity of the genus Malus and implications for linkage mapping with SNPs. Tree Genet Genomes 7: 857-868. doi:https://doi.org/10.1007/s11295-011-0380-8.
  53. 53. Chagné D, Carlisle CM, Blond C, Volz RK, Whitworth CJ et al. (2007) Mapping a candidate gene (MdMYB10)for red flesh and foliage colour in apple. BMC Genomics 8: 212. doi:https://doi.org/10.1186/1471-2164-8-212. PubMed: 17608951.
  54. 54. Espley RV, Hellens RP, Putterill J, Stevenson DE, Kutty-Amma S et al. (2007) Red colouration in apple fruit is due to the activity of the MYB transcription factor, MdMYB10. Plant J 49: 414-427. doi:https://doi.org/10.1111/j.1365-313X.2006.02964.x. PubMed: 17181777.
  55. 55. Chagné D, Krieger C, Rassam M, Sullivan M, Fraser J et al. (2012) QTL and candidate gene mapping for polyphenolic composition in apple fruit. BMC Plant Biol 12: 12. doi:https://doi.org/10.1186/1471-2229-12-12. PubMed: 22269060.
  56. 56. Kumar S, Chagné D, Bink MCAM, Volz RK, Whitworth C et al. (2012) Genomic selection for fruit quality traits in apple (Malus x domestica Borkh.). PLOS_ONE 7(5): 1-10.
  57. 57. Mellidou I, Chagné D, Laing WA, Keulemans J, Davey MW (2012) Allelic variation in paralogs of GDP-L-Galactose phosphoryòase is a major determinant of vitamin C concentration in apple fruit. Plant Physiol 160: 1613-1629. doi:https://doi.org/10.1104/pp.112.203786. PubMed: 23001142.
  58. 58. Franck C, Baetens M, Lammertyn J, Verboven P, Davey MW et al. (2003) Ascorbic acid concentration in cv. Conference pears during fruit development and postharvest storage. Journal of Agricultural Food Chemestry 51, 4757-4763..
  59. 59. Ogundiwin EA, Peace CP, Gradziel TM, Parfitt DE, Bliss FA et al. (2009) A fruit quality gene map of Prunus. BMC Genomics 10: 587. doi:https://doi.org/10.1186/1471-2164-10-587. PubMed: 19995417.
  60. 60. Treutter D (2001) Biosynthesis of phenolic compounds and its regulation in apple. Plant Growth Regul 34: 71-89. doi:https://doi.org/10.1023/A:1013378702940.
  61. 61. Pannala AS, Chan TS, O’Brien PJ, Rice-Evans CA (2001) Flavonoid B-ring chemistry and antioxidant activity: fast reaction kinetics. Biochem Biophys Res Commun 282: 1161–1168. doi:https://doi.org/10.1006/bbrc.2001.4705. PubMed: 11302737.
  62. 62. Regev-Shoshani G, Shoseyov O, Bilkis I, Kerem Z (2003) Glycosylation of resveratrol protects it from enzymic oxidation. Biochem J 374: 157–163. doi:https://doi.org/10.1042/BJ20030141. PubMed: 12697026.
  63. 63. Spanos GA, Wrolstad RE, Heatherbell DA (1990) Influence of processing and storage on the phenolic composition of apple juice. J, Agric Food Chem 38: 1572-1579. doi:https://doi.org/10.1021/jf00097a031.
  64. 64. MacDonald MJ, D’Cunha GB (2007) A modern view of phenylalanine ammonia lyase. Biochemestry Cell Biol 85(3): 273-282. doi:https://doi.org/10.1139/O07-018.
  65. 65. Bagal UR, Leebens-Mack JH, Lorenz WW, Dean JFD (2012) The phenylalanine ammonia lyase (PAL) gene family shows a gymnosperm-specific lineage. BMC Genomics 13: 51. doi:https://doi.org/10.1186/1471-2164-13-51. PubMed: 22296850.
  66. 66. Choi Y-j, Tomás-Barberán , Saltveit ME (2002) Wound-induced phenolic accumulation and browning in lettuce (Lactuca sativa L.) leaf tissue is reduced by exposure to n-alcohols. Postharvest Biol Technol 37: 47-55.
  67. 67. Constabel CP, Bergey DR, Ryan CA (1995) Systemin activates synthesis of wound-inducible tomato leaf polyphenol oxidase via the octadecanoid defense signaling pathway. Proc Natl Acad Sci USA 92: 407-411. doi:https://doi.org/10.1073/pnas.92.2.407. PubMed: 7831300.
  68. 68. Thipyapong P, Steffens JC (1997) Tomato Polyphenol Oxidase (Differential Response of the Polyphenol Oxidase F Promoter to Injuries and Wound Signals). Plant Physiol 1158(2): 409-418.