Flesh texture and stone adherence are important factors of the overall quality of the peach [Prunus persica (L.) Batsch] fruit, contributing to the consumer’s satisfaction, fruit nutritional features, and postharvest behavior   .
Peaches are classified on the basis of their different flesh softening behavior during ripening. In Melting (M) flesh, a marked and steep loss of firmness (“melting” phase) occurs in the final stages of ripening. The melting process is delayed in Slow Softening (SS) fruit and undetectable in Non Melting (NM) and Stony Hard (SH) ones      . The SS phenotype is very interesting, retaining flesh firmness on the tree for longer time than standard M, with full development of sensory qualities, and remarkable keeping quality appreciated by both growers and consumers. The spontaneous occurrence of melting, although a few days later than standard M, allows to group SS fruit into the M phenotype, since the physiological basis of SS is due to a delay in ethylene production  . Unfortunately, this phenotype, like the SH one, is of very difficult assessment on the tree when scoring segregating progenies  .
Fruit flesh texture is affected by biochemical and physiological factors with a major role played by hydrolytic enzymes and other proteins that cooperate to modify the composition and architecture of the cell wall polysaccharides   . In particular, peach flesh melting has been related to a strong increase in the expression of Pp-endopolygalacturonase (Pp-endoPG) and accumulation of an active Pp-endoPG protein. In NM fruit, Pp-endoPG is absent or detected at levels lower than in M ones     . Therefore, Pp-endoPG is considered as a candidate gene for peach flesh softening behavior  .
In peach, availability of different molecular markers allowed to gain information on some 760 Quantitative Trait Loci and Mendelian Trait Loci linked to horticultural and physiological traits such as tree development, pest and disease resistance, flowering, ripening, seed and fruit quality   .
The availability of the peach genome sequence  and the development of SNP genotyping resources   offer the opportunity to study at the molecular level the development and inheritance of different phenotypic traits, as recently reported for skin blush  .
A high degree of allelic diversity has been identified in the peach Pp-endoPG sequence    . We previously identified specific SNPs in the open-reading frame of Pp-endoPG, one of which allowed to determine the configuration at the Melting flesh locus in the NM (Oro A) and M (Bolero) model cultivars, and in a few M and NM F1 seedlings from their cross   . Moreover, in Bolero (M), Oro A (NM), and in the cultivars Big Top (SS), Yumyeong (SH), and Ghiaccio (SH), representative of different flesh textures, variants of the Pp-endoPG gene characterized by peculiar InDels and SNPs have been individuated. The Pp-endoPG locus on chromosome 4 contains two sequences, i.e., Fa-Ppa006839m (GenBank ID: 18781156) and Fb-Ppa006857m (GenBank ID: 18779267), that coincide with Pp-endoPGm (GenBank ID: DQ659240.1) and Pp-endoPG_M (GenBank ID: DQ659241.1)    . Fa-Ppa006839m (Pp-endoPGm), the only sequence retrieved in Oro A, shows 97 % identity with its variants Pp-endoPGSH (GenBank ID: HQ891822.1 of SH Yumyeong and GenBank ID: HQ891821.1 of SH Ghiaccio) and Pp-endoPGBT (GenBank ID: HQ891820.1 of SS Big Top), for the presence in the latter ones of one 17-bp intronic deletion. Moreover, the Big Top-specific Pp-endoPGBT sequence shows a peculiar SNP (bp 348) in Exon 1, that does not induce any change in the deduced amino acidic sequence  . Fb-Ppa006857m (Pp-endoPG_M), present in Bolero together with Fa-Ppa006839m, is by 34 bp shorter than Pp-endoPGm, for the presence of one additional 17-bp intronic deletion. On this basis, amplification of a selected (1455-1892 bp) Pp-endoPG genomic sequence including these InDels generates fragments specifically referable to Fa-Ppa006839m (Pp-endoPGm; 437 bp), to its variants Pp-endoPGSH and Pp-endoPGBT (both of 420 bp), or to Fb-Ppa006857m (Pp-endoPG_M; 403 bp)  . Fa-Ppa006839m has recently been proposed as the main responsible for the Melting trait  .
In the present work the InDel and the CAPS (Cleaved Amplified Polymorphic Sequence) markers developed from the described polymorphisms have been exploited to screen a total of 85 accessions (commercial cultivars, F1 offsprings from selected crosses, and three haploids) to describe them for the configuration at the endoPG locus and evaluate a potential relationship with fruit flesh texture phenotype.
2. Materials and Methods
2.1. Plant Material
The study was conducted on commercial cultivars or on F1 offsprings of different crosses. Three haploid accessions (P VIN 1 1N from Vineland, P RRL 1 1N from Rutgers Red Leaf, P LOV 3 1N from Lovell) were also considered (Table 1).
Plant material was obtained from a peach germplasm collection grown at the Experimental Orchard “Zabina” (Castel San Pietro Terme, Italy; lat. 44˚23'52''N; long. 11˚35'22''E), under the weather and climate conditions detailed in Table 2. Trees were grafted on GFF677 rootstock planted at 3.5 × 4.5 m and trained as open vase. Seedlings were planted on their own roots with a spacing of 1 m within and 4 m between rows and trained as slender spindle (one stem with short lateral scaffolds). Pruning was performed yearly and standard cultural practices were applied.
2.2. Fruit Flesh Texture Phenotype Scoring
Flesh texture was scored in the orchard by trained personnel on fruit harvested at commercial ripening (onset of the veraison stage), as identified in each
Table 1. List of the Prunus persica accessions used, with fruit flesh texture phenotype, pedigree, estimated length of endoPG-derived amplicons, and asset at the endoPG locus as hypothesizable by InDel analysis or, where available, reported in the literature. M: Melting; NM: Non Melting; SH: Stony Hard; SS: Slow Softening. NC, not classifiable. OP, open pollination. -, lack of fruit production.
Table 2. Climate and weather conditions recorded in the orchard during the time period of fruit growing season. n. d.: not determined.
cultivar by visual assessment of the skin ground color   , and sensory determination of approximate firmness.
Assessment of the latter parameter by human inspectors’ grading is coherent with instrumented measurements  , and represents a non-destructive method capable to yield information on the time course of fruit texture evolution.
Discrimination between M and SS phenotypes was based upon sensory scoring of flesh texture changes during 4 d of postharvest at room temperature. Due to stunted growth and low fruit set, in a few F1 offsprings from selected crosses, fruit were not phenotypically classifiable; haploids also had stunted growth and did not produce fruit at all, so the fruit flesh phenotype remained unascertained.
2.3. Extraction of Genomic DNA
Genomic DNA from young leaves was used as template for the analysis of the amplification profiles based on major InDel polymorphisms in Intron 3 of the Pp-endoPG sequence   , and for CAPS analysis  based on the SNP348 of the Pp-endoPGBT variant of Fa-Ppa006839 m  .
Extraction was conducted according to Geuna et al.  . Frozen leaf samples (100 mg fresh weight) were powdered in liquid N2 by a mortar and a pestle and mixed with 600 µL of preheated (65˚C) extraction buffer [0.2 M Tris (tris hydroxymethyl aminomethane)-HCl, pH 8.5, 10 mM EDTA (ethylenediaminetetraacetic acid), 0.3 M LiCl, 1.5 % (v/v) SDS (sodium dodecyl sulphate), 1 % (w/v) Nonidet™ P-40, 1 mM DTT (1,4-dithiothreitol)]. Samples were incubated at 65˚C for 20 min and centrifuged (15,000 g, 4˚C, 20 min). The aqueous solution was extracted in one volume of (in sequence) phenol, phenol/chloroform/isoamyl alcohol (25/24/1 by volume), and chloroform, with centrifugation (15,000 g, 4˚C, 20 min) after each passage. DNA was precipitated by adding 0.1 volume of 3 M Na-acetate, pH 5.2 plus 0.6 volume of 2-propanol, followed by incubation at 4˚C for 1 h and centrifugation (15,000 g, 4˚C, 30 min). Pellets were washed with 1 mL of 70% (v/v) ethanol, air-dried and resuspended in 100 µL of double-distilled H2O. Three microliters of these extracts were analyzed for quality and yield on 1% (w/v) agarose gel in 1 ´ TBE buffer (89 mM Tris-Borate, 2 mM Na2-EDTA, pH 8.2 plus 1 µg mL−1 ethidium bromide; low ionic strength), also used as electrophoresis buffer. Lambda DNA/Hind III genomic DNA and 1 kb DNA Plus Ladder (Invitrogen Life Technologies, Monza, Italy) were run as standards of concentration and size. Nucleic acids were visualized under UV light. The quantified DNA was treated with 10 mg mL−1 RNase A (Invitrogen Life Technologies) at 37˚C for 1 h to remove any RNA contamination and precipitated as described above.
2.4. InDel Analysis
A ~440 bp sequence of the Pp-endoPG gene region encompassing a few bases of Exon 3 and a wide sequence of Intron 3 comprising the major InDels  was amplified with proper primers (Pp-endoPG_InDelFor: 5’-GTGCCCTGGTCAGGTAAG- 3’ ; Pp-endoPG_InDelRev: 5’-GGCTAAGCTACGATGAAGTC-3’) in a MyCycler Thermal Cycler (Bio-Rad Laboratories Srl, Segrate, Italy). The PCR mix contained 20 ng genomic DNA, 0.3 mM dNTPs (deoxynucleotides), 0.3 µM of each primer, 1 × GoTaq Reaction Buffer®, 1 U Go Taq® DNA Polymerase (Promega, Segrate, Italy) and double-distilled H2O to a final volume of 25 mL. The conditions of PCR reaction were: one denaturation cycle (94˚C, 4 min), 35 cycles comprehensive of denaturation (94˚C, 30 s), annealing (62˚C, 30 s) and extension phases (72˚C, 45 s) and one final extension cycle (72˚C, 5 min).
The amplification products were separated on 3% (w/v) agarose gel and visualized (DirectLoad™ 50 bp DNA Step Ladder; Sigma-Aldrich) with ethidium bromide under UV light.
2.5. CAPS Analysis
Primer pairs (Pp-endoPGFor: 5’-ATGGCGAACCGTAGAAGCCTCT-3’; Pp-endoPGBTRev: 5’-CCACAAGCAACGCCTTCTATCC-3’) were designed to amplify the 1-972 bp region of Pp-endoPG including the SNP348 of interest  , which determined the polymorphic restriction site. Amplification was conducted as described above for InDel analysis and was followed by digestion of the polymorphic fragments [1 μg of amplification products, 5 U of BstXI restriction enzyme (Promega), 37˚C, 90 min]. The reaction products, separated on 3% (w/v) agarose gels, were visualized by ethidium bromide (TrackIt™ 1 Kb Plus DNA Ladder; Invitrogen).
3.1. InDel Analysis of Different Peach Accessions
Amplification of the cited genomic sequence of Pp-endoPG including the InDels of interest confirmed the presence of two amplicons of apparent length of ~440 bp and ~410 bp in Bolero, consistent with the expected lengths of 437 bp and 403 bp for Fa-Ppa006839m (Pp-endoPGm) and Fb-Ppa006857m (Pp-endoPG_M)   .
The M cultivars Contender, Maycrest, Springbelle, Springcrest, Springred, Glohaven, plus two (BO 96016136 and BO 96016094) out of the seven offsprings from the cross Contender (M) × Ambra (M) showed a Bolero-like amplification pattern with two fragments. In Ambra and in five offsprings from Contender × Ambra (BO 96016015, BO 96016018, BO 96016023, BO 96016165, and BO 96016208) an additional amplicon of ~420 bp was observed (Figure 1).
Figure 2 shows the results obtained in NM, SH, and SS accessions. In the NM Oro A, Andross  , Ionia, and BO 82010554 (offspring of NM Jungerman × NM Loadel), only one amplicon of ~440 bp was retrieved. In the SH accessions Yumyeong, Ghiaccio, Helena Cling, and D 41-62, as well as in SS Big Top, amplification generated only one fragment of ~420 bp, confirming the presence of the Pp-endoPGSH-or Pp-endoPGBT-like variants of Fa-Ppa006839m, respectively  .
The amplification pattern of SS Alitop was similar to that of Bolero. SH 7-28 and the SS Honey Gold, Ruby Rich, Vista Rich, BO 96028059 (F1 from M Springred × SS Big Top) and BO 96013046 (F1 from M Bolero × SS Rich Lady) yielded three amplicons (~440 bp, ~420 bp, and ~410 bp). The two haploids P VIN 1 1N and P LOV 3 1N yielded amplification patterns similar to that of Oro A (a single fragment of ~440 bp), whereas the pattern of P RRL 1 1N was similar to that of Bolero (two fragments of ~440 bp and ~410 bp; Figure 3).
Figure 1. InDel analysis of Pp-endoPG in accessions with M fruit flesh texture. Positions and lengths of DNA markers (DirectLoad™ 50 bp DNA Step Ladder; Sigma-Aldrich) are shown on the left. Twenty micrograms DNA were loaded per lane; the amplicons were separated on 3% (w/v) agarose gels. One representative gel from three independent experiments.
Figure 2. InDel analysis of Pp-endoPG in accessions with NM, SH or SS fruit flesh texture. Positions and lengths of DNA markers (DirectLoad™ 50 bp DNA Step Ladder) are shown on the left. Twenty micrograms DNA were loaded per lane; the amplicons were separated on 3% (w/v) agarose gels. One representative gel from three independent experiments.
3.2. Ambra × Big Top and Max 7 × Big Top Crosses: InDel Analysis of F1 Offsprings
Figure 4 shows the amplification patterns obtained in the M Ambra and SS Big Top parents and in 23 F1 seedlings from their cross (from BO 02002001 to BO 02002006, BO 02002008, BO 02002009, from BO 02002020 to BO 02002026, from BO 02002011 to BO 02002015, from BO 02002017 to BO 02002019). The amplification profiles were similar to that of the Big Top parent (one amplicon of ~420 bp) for eight of them, and to that of the Ambra parent (three amplicons of ~440 bp, ~420 bp, and ~410 bp) for 12. Three accessions showed the same pattern as Bolero (~440 bp and ~410 bp; compare to Figure 1).
Figure 5 shows the amplification profiles obtained in the M Max 7 and SS Big Top parents and in 27 F1 seedlings from their cross (from BO 02004001 to BO 02004027). In Max 7, as well as in eight offsprings, amplification yielded three fragments (~440 bp, ~420 bp, and ~410 bp). Amplification profiles similar to that of Big Top (one fragment of ~420 bp) were found in 13 seedlings, whereas six seedlings yielded Bolero-like (~440 bp and ~410 bp fragments) patterns (compare to Figure 1).
Figure 3. InDel analysis of Pp-endoPG in three haploid accessions. Positions and lengths of DNA markers (50 bp DNA Step Ladder; Sigma-Aldrich) are shown on the left. Twenty micrograms DNA were loaded per lane; the amplicons were separated on 3% (w/v) agarose gel. One representative gel from three independent experiments.
Figure 4. InDel analysis of Pp-endoPG in M Ambra, SS Big Top and in 23 F1 seedlings. Positions and lengths of DNA markers (DirectLoad™ 50 bp DNA Step Ladder) are shown. NC, not classifiable. Twenty micrograms DNA were loaded per lane; amplicons were separated on 3% (w/v) agarose gel. One representative gel from three independent experiments.
Figure 5. InDel analysis of Pp-endoPG in M Max 7, SS Big Top, and in 27 F1 seedlings. Positions and lengths of DNA markers (DirectLoad™ 50 bp DNA Step Ladder) are shown on the left. NC, not classifiable. Twenty micrograms DNA were loaded per lane; amplicons were separated on 3% (w/v) agarose gel. One representative gel from three independent experiments.
3.3. Cleaved Amplified Polymorphic Sequence (CAPS) Analysis
In the Pp-endoPGBT variant of Fa-Ppa006839m present in Big Top, a peculiar, silent SNP348 in Exon 1 originates a polymorphic sequence recognized by the BstXI restriction enzyme (Figure 6(b)), suitable for distinguishing between the Pp-endoPGBT and Pp-endoPGSH variants of Fa-Ppa006839m (Pp-endoPGm)  .
The SNP348 was exploited to broaden CAPS analysis to a few additional cultivars/accessions phenotypically scored as SS (Alitop, Helena Cling, Honey Gold, Ruby Rich, Vista Rich, BO 96028059, BO 96013046) or SH (D 41-62, 7-28) (Table 1). BstXI confirmed to be ineffective in Oro A and in Bolero, as it was in all the accessions scored as SH. Concerning the accessions scored as SS, Helena Cling showed, like Big Top, two digestion fragments; Honey Gold, Ruby Rich, Vista Rich, BO 96028059, and BO 96013046 showed two Big Top-like restriction fragments plus an undigested one, whereas BstXI was ineffective in Alitop (Figure 6(c)).
The F allele at the endoPG locus on peach chromosome 4 contains two gene sequences, i.e. Fa (Ppa006839m) and Fb (Ppa006857m), at short distance (32-34 kbp) from each other. Mutations and deletions of these genes determine several allelic variants: F, resulting from Fa plus Fb; f, resulting from mutation of Fa plus complete deletion of Fb; f1, resulting from Fa plus complete deletion of Fb; fnull, resulting from deep mutation or complete deletion of both Fa and Fb. The different alleles F, f, f1, fnull at the endoPG locus contribute to the fruit phenotype for stone adhesion and flesh texture   . In particular, it has been recently suggested that Fa-Ppa006839m plays a driving role in the determination of the Melting/Non Melting trait, whereas Fb-Ppa006857m seems involved in the determination of the Clingstone/Freestone trait, although it has been speculated that Fb-Ppa006857m could exert a pleiotropic effect on flesh melting through a negative feedback control on the transcription of Fa-Ppa006839m  .
From reference peach genotypes, peculiar endoPG sequences have been isolated. In particular, in NM Oro A has been retrieved a single sequence (Pp-endoPGm) coincident with Ppa006839m, that does not present any deletion
Figure 6. CAPS analysis exploiting the SNP348 of the Pp-endoPGBT variant of Fa-Ppa006839m. (a) Schematic representation of the structure of the Pp-endoPG gene. Exons and introns are indicated by black solid bars and by solid lines, respectively; positions of the considered InDels of Intron 3 are also indicated. Arrows and vertical dotted lines define the region (1-972 bp) amplified by the primers used. (b) Lengths (~972 bp) of the amplicons obtained from the Fa-Ppa006839m variants of: NM Oro A and M Bolero (Pp-endoPGm), SH Yumyeong and SH Ghiaccio (Pp-endoPGSH), and SS Big Top (Pp-endoPGBT); position of the peculiar SNP348 of the Pp-endoPGBT variant and predicted lengths (348 bp and 626 bp) of the fragments obtainable by BstXI digestion (scissors symbol) are also indicated (modified from  ). The Pp-endoPG_M amplicon from Fb-Ppa006857m of Bolero is also reported. (c) Restriction patterns obtained in different accessions before (-) and after (D) digestion with BstXI. Positions and lengths of DNA markers (1 Kb Plus DNA Ladder, Invitrogen) are shown on the right. Expected lengths (bp) of the undigested amplicons and of the products of BstXI digestion are shown on the left. Twenty micrograms DNA were loaded per lane; 3% (w/v) agarose. One representative gel from three independent experiments.
in Intron 3, as reported also for the Fa and f1 endoPG selected sequences. The allelic variants (Pp-endoPGBT/SH) of Fa-Ppa006839m, isolated from SS Big Top and SH Yumyeong and SH Ghiaccio, present a specific 17-bp deletion in Intron 3, apparently coincident with the deletion of the f sequence. In M Bolero has been isolated a sequence (Pp-endoPG_M) coincident with Fb-Ppa006857m and presenting two 17-bp deletions in Intron 3 as reported for Fb    . Figure 7 summarizes the structure of the different alleles at the endoPG locus.
Genotype-specific Pp-endoPG polymorphisms (InDels and SNPs;   ) were exploited in the present work as tools to describe 85 peach accessions (Table 1). InDel analysis confirmed the co-segregation of Fa-Ppa006839m (Pp-endoPGm) and Fb-Ppa006857m (Pp-endoPG_M). In fact, the longest (~440 bp) amplicon from Fa-Ppa006839m was always accompanied by the shortest (~410 bp) one from Fb-Ppa006857m, as indicated by the simultaneous presence of three amplicons in several accessions and of both amplicons in the haploid P RRL 1 1N
Figure 7. Schematic picture of the structure of the Pp-endoPG locus with the different allelic variants and the discussed InDel mutations, as deduced from     . C: Clingstone; F: Freestone. Exons are indicated by black solid bars and introns by solid lines.
(Figures 1-5). In the doubled haploid P LOV2 2N (from open pollination of Lovell) used for genome sequencing both the Fa-Ppa006839m and the Fb-Ppa006857m sequences are present   , while only the longest amplicon (~440 bp) was retrieved in the haploid P LOV3 1N (Figure 3) suggesting that Lovell is heterozygous (F/f1) at the endoPG locus.
The individuated InDel polymorphisms in the endoPG sequences seem to distinguish the NM trait from the M one. In all the NM accessions (Oro A, Andross, Ionia, and BO 82010054) only the longest amplicon (~440 bp; Figure 2) referable to a Fa-Ppa006839m-like variant is present. This result confirms what reported for NM phenotypes as determined by the presence of the f1 allele, in homozygosis (f1/f1) or in heterozygosis with fnull (f1/fnull)    (Figure 7; Table 1).
Most (31) of the 64 accessions scored as M or SS (all grouped in a general Melting class   ) were characterized by presence of all the three amplicons (~440 bp, ~420 bp and ~410 bp), referable respectively to Fa-Ppa006839m (no deletion), its Pp-endoPGBT/SH-like variants (one 17-bp deletion), and Fb-Ppa006857m (two 17-bp deletions). Eighteen accessions showed a pattern with the two amplicons of ~440 bp and ~410 bp, and 15 showed only the amplicon of ~420 bp (Figures 1-4; Table 1). From these amplification profiles, it seems reasonable to deduce that the M and SS accessions can present all the three combinations of amplicons.
The InDel polymorphism, when scored in a cross progeny, appears useful to deduce the allelic combination at the endoPG locus of the parental genotypes, as in the case of Big Top. In fact, a few F1 individuals of the crosses Ambra × Big Top or Max 7 × Big Top, as well as Alitop [(Flavortop × Snow Queen) × Big Top], lacked the amplicon of ~420 bp (Figure 2, Figure 4 and Figure 5). These results can only be explained by hypothesizing that Big Top is heterozygous at the endoPG locus for the presence of a “null” allele, in addition to the f allele composed only by the Pp-endoPGBT variant of Fa-Ppa006839m. The literature proposes for SS Big Top the allelic combinations f/f or f/fnull   ; our results seem to support the latter one (Figure 4 and Figure 5; Table 1).
The coherence of the allelic combinations at the endoPG locus inferred by the InDel analysis with, where available, those reported in the literature (Table 1) suggests that the InDel polymorphism may represent a tool for a simple genotyping of peach accessions. Nevertheless, this polymorphism does not allow to discriminate between the peculiar gene sequences individuated in SS or SH genotypes.
CAPS analysis exploiting the SNP individuated in Exon 1 of the endoPG sequence of SS Big Top integrates the InDel results allowing to distinguish between amplicons of ~420 bp derived from the Pp-endoPGBT-like variant of Fa-Ppa006839m and the Pp-endoPGSH-like one  . All the SS genotypes tested in CAPS analysis proved sensitive to BstXI producing upon cleavage two digestion fragments. In a few accessions, the contemporary presence of an undigested fragment indicated that their allelic asset was heterozygous (F/f; Figure 6(c) and Figure 2). For Ruby Rich, this asset confirms the functional genotype (F/f) proposed in the literature   (Table 1). Lack of effect of BstXI in SS Alitop further confirms the inheritance of a possible ‘null’ Pp-endoPG variant from the Big Top parent. Conversely, the SH genotypes tested were not sensitive to BstXI digestion (Figure 6(c)).
In the climacteric peach fruit the evolution of texture characteristics is regulated by ethylene. The SH trait (hdhd) is related to absence of ethylene production for a mutation in the Pp-ACS1 gene and altered indole-3-acetic acid levels due to regulation of PpYUC11     . In SH 7-28, homozygous-recessive for the stony hard gene (hdhd)  , InDel analysis shows presence of the ~420 bp amplicon (Figure 2) attributable, by CAPS analysis, to a Pp-endoPGSH-like sequence (Figure 6(c)). The contemporary presence of the other two amplicons allows to assign this accession, as far as it concerns Pp-endoPG, to the Stony Hard melting category (hdhd M-) consistent to Haji et al.  .
Therefore, it can be concluded that although CAPS analysis of the specific Pp-endoPGBT/SH variants seems to allow a preliminary distinction of SS from SH genotypes, conclusive phenotypic attribution of a specific accession will necessarily require the assessment of other activities like fruit ethylene production.
The two molecular markers (InDel and CAPS) based upon Pp-endoPG polymorphisms tested in the present work may represent useful tools for simple genotyping of peach accessions as far as it concerns the characteristics and different allelic combinations at the endoPG locus.
In particular, the use of the Pp-endoPG InDel marker yields results consistent with the existing knowledge on the polymorphic structure of the endoPG locus in peach, and allows to individuate NM accessions. The CAPS marker developed on the basis of the peculiar SNP348 of SS Big Top integrates the results of InDel analysis by discriminating the Pp-endoPG variants in SS and SH accessions.
Future development of the work may foresee the testing of a broader number of NM, as well as SS and SH accessions, in order to confirm the potential use of the described InDel and CAPS markers in Marker-Assisted Selection peach breeding programs.
This work was funded by the Italian Ministry of University and Research (MIUR) PRIN-COFIN 2006 grant number 2006078852, and by the MAS.PES Project (public-private breeding program supported by Italian growers’ organizations and the University of Milan, Italy). The authors thank S. Foschi and M. Lama for their technical assistance in field operations.
 Bruhn, C.M. (1995) Consumer and Retailer Satisfaction with the Quality and Size of California Peaches and Nectarines. Journal of Food Quality, 18, 241-256.
 Bassi, D. and Monet, R. (2008) 1. Botany and Taxonomy. In: Layne, D.R. and Bassi, D., Eds., The Peach: Botany, Production and Uses, CAB International, Wallingford, 1-36.
 Yoshida, M. (1976) Genetic Studies on the Fruit Quality of Peach Varieties. III. Texture and Keeping Quality. Bulletin of the Fruit Tree Research Station, 3, 1-16. (In Japanese with English Abstract)
 Peace, C.P., Crisosto, C.H. and Gradziel, T.M. (2005) Endopolygalacturonase: A Candidate Gene for Freestone and Melting Flesh in Peach. Molecular Breeding, 16, 21-31.
 Morgutti, S., Negrini, N., Nocito, F.F., Ghiani, A., Bassi, D. and Cocucci, M. (2006) Changes in Endopolygalacturonase Levels and Characterization of a Putative Endo-PG Gene During Fruit Softening in Peach Genotypes with Nonmelting and Melting Flesh Fruit Phenotypes. New Phytologist, 171, 315-328.
 Ghiani, A., Negrini, N., Morgutti, S., Baldin, F., Nocito, F.F., Spinardi, A., et al. (2011) Melting of “Big Top” Nectarine Fruit: Some Physiological, Biochemical, and Molecular Aspects. Journal of the American Society for Horticultural Science, 136, 61-68.
 Byrne, D.H., Raseira, M.B., Bassi, D., Piagnani, M.C., Gasic, K., Reighard, G.L., et al. (2012) Peach. In: Badenes, M.L. and Byrne, D.H., Eds., Fruit Breeding, Series Handbook of Plant Breeding 8, Springer, New York, 505-569.
 Redgwell, R.J. and Fischer, M. (2002) 3. Fruit Texture, Cell Wall Metabolism and Consumer Perceptions. In: Knee, M., Ed., Fruit Quality and Its Biological Basis, Sheffield Academic Press, Sheffield, 46-88.
 Negi, P.S. and Kanda, A.K (2008) 8. Structural Deterioration of the Produce: The Breakdown of Cell Wall Components. In: Paliyath, G., Murr, D.P., Handa, A.K. and Lurie, S., Eds., Postharvest Biology and Technology of Fruits, Vegetables, and Flowers, Wiley-Blackwell, Ames, 162-194.
 Lester, D.R., Speirs, G., Orr, G. and Brady, C.J. (1994) Peach (Prunus persica) Endo-PG cDNA Isolation and mRNA Analysis in Melting and Non-Melting Peach Cultivars. Plant Physiology, 105, 225-231.
 Callahan, A.M., Scorza, R., Bassett, C., Nickerson, M. and Abeles, F.B. (2004) Deletions in an Endopolygalacturonase Gene Cluster Correlate with Non-Melting Flesh Texture in Peach. Functional Plant Biology, 31, 159-168.
 Peace, C.P., Callahan, A., Ogundiwin, E.A., Potter, D., Gradziel, T.M., Bliss, F.A., et al. (2007) Endopolygalacturonase Genotypic Variation in Prunus. In: Litz, R.E. and Scorza, R., Eds., International Symposium on Biotechnology of Temperate Fruit Crops and Tropical Species, Acta Horticulturae, 738, 639-646.
 Zou, X., Shi, C., Austin, R.S., Merico, D., Munholland, S., Marsolais, F., et al. (2014) Genome-Wide Single Nucleotide Polymorphism and Insertion-Deletion Discovery through Next-Generation Sequencing of Reduced Representation Libraries in Common Bean. Molecular Breeding, 33, 769-778.
 Agarwal, M., Shrivastava, N. and Padh, H. (2008) Advances in Molecular Marker Techniques and their Applications in Plant Science. Plant Cell Reports, 27, 617-631.
 Chavez, D.J., Beckman, T.G., Werner, D.J. and Chaparro, J.X. (2014) Genetic Diversity in Peach (Prunus persica (L.) Batsch) at the University of Florida: Past, Present and Future. Tree Genetics & Genomes, 10, 1399-1417.
 Salazar, J.A., Ruiz, D., Campoy, J.A., Sánchez-Pérez, R., Crisosto, C.H., Martínez-García, P.J., et al. (2014) Quantitative Trait Loci (QTL) and Mendelian Trait Loci (MTL) Analysis in Prunus: A Breeding Perspective and Beyond. Plant Molecular Biology Reporter, 32, 1-18.
 Verde, I., Abbott, A.G., Scalabrin, S., Jung, S., Su, S., Marroni, F., et al. (2013) The High-Quality Draft Genome of Peach (Prunus persica) Identifies Unique Patterns of Genetic Diversity, Domestication and Genome Evolution. Nature Genetics, 45, 487-494.
 Ahmad, R., Parfitt, D.E., Fass, J., Ogundiwin, E., Dhingra, A. and Crisosto, C.H. (2011) Whole Genome Sequencing of Peach (Prunus persica L.) for SNP Identification and Selection. BMC Genomics, 12, 569.
 Verde, I., Bassil, N., Scalabrin, S., Gilmore, B., Lawley, C.T., Gasic, K., et al. (2012) Development and Evaluation of a 9K SNP Array for Peach by Internationally Coordinated SNP Detection and Validation in Breeding Germplasm. PLoS ONE, 7, e35668.
 Frett, T.J., Reighard, G.L., Okie, W.R. and Gasic, K. (2014) Mapping Quantitative Trait Loci Associated with Blush in Peach (Prunus persica (L.) Batsch). Tree Genetics & Genomes, 10, 367-381.
 Gu, C., Wang, L., Wang, W., Zhou, H., Ma, B., Zheng, H., et al. (2016) Copy Number Variation of a Gene Cluster Encoding Endopolygalacturonase Mediates Flesh Texture and Stone Adhesion in Peach. Journal of Experimental Botany, 67, 1993-2005.
 Crisosto, C.H. and Valero, D. (2008) 22. Harvesting and Postharvest Handling of Peaches for the Fresh Market. In: Layne, D.R. and Bassi, D., Eds., The Peach: Botany, Production and Uses, CAB International, Wallingford, 575-596.
 Slaughter, D.C., Crisosto, C.H., Hasey, J.K. and Thompson, J.F. (2006) Comparison of Instrumental and Manual Inspection of Clingstone Peaches. Applied Engineering in Agriculture, 22, 883-889.
 Lukowitz, W., Gillmor, C.S. and Scheible, W. (2000) Positional Cloning in Arabidopsis. Why It Feels Good to Have a Genome Initiative Working for You. Plant Physiology, 123, 795-806.
 Geuna, F., Maitti, C., Digiuni, S. and Banfi, R. (2004) A Method for Extracting Genomic DNA Suitable for Medium-Throughput Applications from Plant Tissues Rich in Contaminants. Plant Molecular Biology Reporter, 22, 87.
 Warburton, M.L. and Bliss, F.A. (1996) Genetic Diversity in Peach (Prunus persica L. Batsch) Revealed by Randomly Amplified Polymorphic DNA (RAPD) Markers and Compared to Inbreeding Coefficients. Journal of the American Society for Horticultural Science, 121, 1012-1019.
 Iezzoni, A., Weebadde, C., Luby, J., Chengyan, Y., van de Weg, E., Fazio, G., et al. (2010) RosBREED: Enabling Marker-Assisted Breeding in Rosaceae. Acta Horticulturae, 859, 389-394.
 Tatsuki, M., Haji, T. and Yamaguchi, M. (2006) The Involvement of 1-Aminocyclopropane-1-Carboxylic Acid Synthase Isogene, Pp-ACS1, in Peach Fruit Softening. Journal of Experimental Botany, 57, 1281-1289.
 Hayama, H., Tatsuki, M., Ito, A. and Kashimura, Y. (2006) Ethylene and Fruit Softening in the Stony Hard Mutation in Peach. Postharvest Biology and Technology, 41, 16-21.
 Pan, L., Zeng, W., Niu, L., Lu, Z., Liu, H., Cui, G., et al. (2015) PpYUC11, a Strong Candidate Gene for the Stony Hard Phenotype in Peach (Prunus persica L. Batsch), Participates in IAA Biosynthesis during Fruit Ripening. Journal of Experimental Botany, 66, 7031-7044.
 Haji, T., Yaegaki, H. and Yamaguchi, M. (2005) Inheritance and Expression of Fruit Texture Melting, Non-Melting and Stony Hard in Peach. Scientia Horticulturae, 105, 241-248.