Genotyping by Sequencing for SNP-Based Linkage Map Construction and QTL Analysis of Chilling Requirement and Bloom Date in Peach [Prunus persica (L.) Batsch]

Douglas Gary Bielenberg; Bradley Rauh; Shenghua Fan; Ksenija Gasic; Albert Glenn Abbott; Gregory Lynn Reighard; William R. Okie; Christina Elizabeth Wells

doi:10.1371/journal.pone.0139406

Abstract

Low-cost, high throughput genotyping methods are crucial to marker discovery and marker-assisted breeding efforts, but have not been available for many ‘specialty crops’ such as fruit and nut trees. Here we apply the Genotyping-By-Sequencing (GBS) method developed for cereals to the discovery of single nucleotide polymorphisms (SNPs) in a peach F₂ mapping population. Peach is a genetic and genomic model within the Rosaceae and will provide a template for the use of this method with other members of this family. Our F₂ mapping population of 57 genotypes segregates for bloom time (BD) and chilling requirement (CR) and we have extensively phenotyped this population. The population derives from a selfed F₁ progeny of a cross between ‘Hakuho’ (high CR) and ‘UFGold’ (low CR). We were able to successfully employ GBS and the TASSEL GBS pipeline without modification of the original methodology using the ApeKI restriction enzyme and multiplexing at an equivalent of 96 samples per Illumina HiSeq 2000 lane. We obtained hundreds of SNP markers which were then used to construct a genetic linkage map and identify quantitative trait loci (QTL) for BD and CR.

Citation: Bielenberg DG, Rauh B, Fan S, Gasic K, Abbott AG, Reighard GL, et al. (2015) Genotyping by Sequencing for SNP-Based Linkage Map Construction and QTL Analysis of Chilling Requirement and Bloom Date in Peach [Prunus persica (L.) Batsch]. PLoS ONE 10(10): e0139406. https://doi.org/10.1371/journal.pone.0139406

Editor: Yuepeng Han, Wuhan Botanical Garden of Chinese Academy of Sciences, CHINA

Received: July 23, 2015; Accepted: September 11, 2015; Published: October 2, 2015

This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This research was supported by the Israel United States Binational Agricultural Research and Development program project US-3746-05R to AGA, GLR and DGB and USDA-NRI CSREES grant 2007-35304-17896 to DGB. The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

High-throughput sequencing platforms such as the Illumina HiSeq® have inspired the development of methods for cost-efficient, genome-wide genotyping of numerous individuals in genetic mapping and population scale studies [1]. These methods achieve cost efficiency through deep multiplexing large numbers of individuals on a sequencing lane and targeting a limited percentage of the genome associated with restriction enzyme loci distributed across all chromosomes. One such method, Genotyping-by-Sequencing (GBS), was developed for low-cost, high-throughput genotyping in maize and has subsequently been used in a number of crop species [2]. GBS is particularly attractive because of the relatively simple bench protocol involved in generating bar-coded samples, which are then combined into a single library for sequencing on the Illumina platform [2]. Our aim was to evaluate the effectiveness with which GBS could be applied to a perennial tree fruit species such as peach [Prunus persica (L) Batsch] [3, 4]. Successful application of GBS to peach would suggest that the method is readily applicable to other economically important forest, fruit, and nut species within the Rosaceae.

Chilling requirement (CR) refers to the minimum duration of cold exposure required before dormant buds will bloom in response to bud break-inducing conditions [5]. CR contributes strongly to bloom date (BD), although BD is also affected by an endogenous heat requirement [6]. Growers and breeders must select cultivars whose CR and BD closely match local climatic conditions in order to avoid crop losses due to late frosts or poor bud break due to insufficient chilling [7]. CR and BD therefore impose a constraint on the introduction and spread of new cultivars with superior agronomic performance and marketability.

CR and BD are known to be quantitative genetic traits whose phenotypes vary widely among peach cultivars [8, 9]. Since CR and BD can only be evaluated after individuals reach reproductive maturity (3–4 years for peach), identification of genetic markers linked with CR and BD phenotypes would save time and resources by allowing selection of genotypes at the seedling stage. As a test of the GBS method in Prunus, we genotyped a F₂ peach mapping population segregating for CR and BD for which we have multiple years of phenotypic data. Here we present the successful use of GBS in peach and the identification of multiple QTLs for CR and BD in our F₂ population.

Methods

Plant material

A peach F₂ population of 57 genotypes was developed at the USDA-ARS Southern Fruit and Tree Nut Research Laboratory (Byron, GA, U.S.A.) by crossing the high chill requirement cultivar Hakuho with the low chill requirement cultivar UFGold and selfing its F₁ hybrid. The female grandparent ‘Hakuho’ is a white-fleshed commercial peach cultivar originating from Japan that requires approximately 900 chilling hours for spring bloom [9]. The male grandparent ‘UFGold’ is a yellow-fleshed commercial cultivar released by the University of Florida (Gainesville, FL, U.S.A.) breeding program [10] that requires approximately 400 chilling hours for spring bloom.

F₂ seeds were stratified, germinated and planted in a greenhouse in fall 2003, then transplanted in spring 2004 to Clemson University’s Musser Fruit Research Center (Seneca, SC, U.S.A.). Replicate clones of each tree were produced by rooting current-year stem cuttings during the 2005 growing season and replicate plantings were established on their own roots at the Musser Fruit Research Center in spring 2006. Replicates of each grandparent were also included in these plantings, although the F₁ tree was no longer available. ‘Hakuho’ was propagated from the original tree used in the cross. Twenty ‘UF Gold’ trees were obtained from a commercial nursery grafted onto ‘Nemaguard’ rootstock.

Phenotyping

Bloom date was visually scored as the date at which 50% of the floral buds on an individual tree reached full bloom stage. Bloom date of all F₂ individuals and grandparents was scored in the original 2004 planting each spring from 2006 through 2012. Observations of bloom progression were made every 2–3 days from the onset of floral bud break. Bloom date was recorded as the number of days from January 1^st of each year.

Chilling requirement of the F₂ population and grandparents was determined in two successive winters (2008/2009 and 2009/2010) using the procedures detailed in our previous publication [8]. In brief, temperature-recording data loggers were placed in the canopy of replicate trees and the average temperature was recorded at 10 minute intervals from mid-October through full bloom in late March. Sampling was performed approximately every 100 hours of accumulated chill time below 7.2°C beginning at 200 hours [11].

On each sampling date three branches (>40cm in length with floral buds) were harvested from each of three replicate trees from each genotype. Branch cuttings were brought to a 25°C greenhouse at Clemson University, recut under water, placed in a 1% Floralife solution (Floralife, Inc., Walterboro, SC, U.S.A.) and maintained under a 16 hour photoperiod. Branches were recut and fresh solution supplied after seven days. At 14 d, % floral bud break was scored on all cuttings. A genotype’s CR was considered to be fulfilled when 50% of the flower buds on all cuttings had opened sufficiently for the petals to be visible.

Genotyping

Genomic DNA for SSR marker screening was isolated as described previously [8]. Genomic DNA for SNP identification was isolated from powdered, freeze-dried leaf tissue of all F₂ individuals and grandparents using the Dellaporta et al. miniprep method [12]. DNA quality was assessed by 260 nm/280 nm absorbance ratios with a Biophotometer 6131 (Eppendorf, Hauppauge, NY). DNA was quantified using the QuantiFluor dsDNA labeling system (Promega, Madison, WI) with a TBS Mini-Fluorometer (Turner Biosystems, Sunnyvale, CA).

A set of 370 simple sequence repeat (SSR) markers from Prunus species were screened for polymorphism in the F₂ mapping population. Markers from all eight major peach linkage groups were selected based on known locations on the T x E Prunus reference map [13] and the peach ‘bin map’ [14]. SSR screening was performed as described in Fan et al. [8]. SSR origins and references are as described in Fan et al. [8].

We followed the ‘Genotyping By Sequencing’ (GBS) method of Elshire et al. [2] to generate ApeKI-associated DNA fragments for sequencing on the Illumina HiSeq® 2000 platform. Ninety-six double-stranded forward adaptors each with a unique barcode and a single common double-stranded reverse adaptor were created from a set of 194 single-stranded oligonucleotides (IDT, Coralville, IA, U.S.A.). Each adaptor contained a three base overhang for ligation with ApeKI digested DNA. The ApeKI compatible barcode set was that published in Elshire et al. [2].

Adaptors were ligated to restriction-digested DNA from F2 individuals and grandparents following the methods described in [2]. Pooled, amplified libraries were sent to the Cornell University Biotechnology Resource Center for single-end sequencing on the Illumina HiSeq® 2000 platform. Grandparental samples were run in triplicate (three separately barcodes samples), whereas F₂ individuals were run once. Ninety-six barcoded samples were run on a single Illumina lane, including some samples from a separate experiment. Samples were re-sequenced if less thant 100,000 reads with identifiable barcodes were obtained from the initial run.

Preliminary inspection of the raw fastq files showed low recovery of five expected barcodes. Using Homertools [15] and basic Linux commands, we extracted barcode sequences from all reads and determined the frequency with which each barcode appeared in the sequence file (S1 File). This analysis revealed a systematic sequencing error that produced a single ‘N’ in the seventh position of several of the barcodes. Due to the robust design of the barcode sequences, we were able to unambiguously assign these approximately 1.5 million mis-barcoded reads to their proper samples.

Processing of sequenced reads

Sequenced reads were processed using default parameters of the TASSEL 3.0 and 4.0 GBS pipeline obtained from the Maize Genetics and Diversity Lab (www.maizegenetics.net) at Cornell University [16]. Assembled scaffolds of the peach genome v1.0 were downloaded from the Genome Database for Rosaceae and used as the reference sequence for alignment of sequenced reads with Bowtie v2.1 [17, 18]. SNP calls for all genotypes were exported as vcf files for filtering with vcftools [19]. SNP calls within a genotype were filtered, retaining only those with a minimum depth of five reads and a genotype quality score of 98 [20]. Filtered vcf files were imported to the TASSEL GUI for visualization of results prior to exporting the SNP calls as a spreadsheet [21].

SNPs were named according to scaffold and base pair position within the peach genome v1.0 build. SNP names contain the scaffold number, an underscore, and eight characters denoting the base position. Leading zeroes in the base position were retained, and a SNP at base 1,234,567 of scaffold 2 would therefore have been named ‘2_01234567’.

Genetic map construction and QTL discovery

SNP data were converted to the ‘a, h, b’ codes with the female ‘Hakuho’ grandparent conferring the ‘a’ genotype. Because the F₁ individual of this cross was deceased, SNPs in the F₂ population were only used if the grandparental source of the alleles could unambiguously be assigned (i.e., one grandparent was homozygous at the SNP locus).

Genetic mapping of the F₂ population was performed using JoinMap® 4.1 [22]. Highly similar markers (>0.95) were excluded from the data set to reduce calculation time. Where applicable, SSRs were preferentially retained from groups of SSR and SNP markers sharing similar segregation. Remaining markers were grouped using the ‘independence LOD’ function in JoinMap® with default settings. Marker groupings were manually verified by inspecting the grouped markers for agreement with known physical locations (SNPs) and known linkage group locations on the ‘T×E’ almond × peach reference map (SSRs) [13, 14, 17]. Marker order and distances were calculated using the Maximum Likelihood Mapping function with default settings. Segregation distortion of individual markers was calculated using the χ² test in JoinMap. Markers with significant segregation distortion (p<0.01) were excluded if surrounding markers did not also show significant segregation distortion of a similar direction and magnitude.

QTLs were detected with MapQTL® 6 [23]. Genome-wise LOD significance thresholds were determined independently for each trait using the ‘Permutation Test’ function. The ‘Automatic Cofactor Selection’ tool was used iteratively to identify the strongest marker cofactors on each linkage group for each trait. Resulting cofactors were included in the search for QTLs that exceeded the LOD significance threshold. LOD curves were plotted with the MQM QTL detection algorithm.

Results

Sequencing and identification of SNPs

Illumina sequencing of 63 pooled, barcoded samples (1 sample per F2 individual, 3 samples per grandparent) generated approximately 1.5 million single-end 100 bp reads per sample. The final mean read number per genotype was 1,843,261 (+/- 1,092,631 SD) after merging genotypes with multiple barcodes and resequencing samples with low initial read numbers.

The TASSEL GBS pipeline initially identified 9,998 SNP loci distributed across all major scaffolds of the peach genome. Prior to the creation of a linkage map, data were filtered to remove SNPs with low read support. Within a genotype, data for individual SNPs were retained only if they possessed a minimum read depth of five and a minimum quality score of 98. After depth and quality filtering, SNPs for which data were missing in more than five genotypes were removed from further analysis. SNPs were also removed in grandparental genotype data were missing, if both grandparents shared a heterozygous genotype or if the minor allele frequency was less than 0.20. These filtering procedures resulted in 410 final SNP loci.

To estimate the effect of sequencing depth on SNP identification, we included the grandparental genotypes in triplicate in the sequencing pool (i.e. each grandparental genotype was tagged with three separate barcodes). This increased the combined sequencing depth of the grandparents threefold relative to the F₂ individuals. Deeper sequencing of the grandparents (‘Hakuho’: 4.2 million reads, ‘UFGold’ 4.6 million reads) resulted in the identification of 4,833 final SNPs, approximately ten times more SNPs than were identified in the analysis of the F₂ population. It should be noted that the decreased number of samples in this analysis (two grandparents vs. 57 F2s) results in fewer SNPs being discarded due to frequency of missing data across all individuals. The distribution of these grandparental SNPs on the physical map is shown in S1 Fig.

Linkage mapping

Of the initial 370 SSR markers screened, thirty-seven were identified as polymorphic in the F₂ population. SSRs were combined with SNP markers identified from the TASSEL GBS pipeline and used to create a genetic linkage map. Following removal of SNPs with highly similar segregation or locally discontinuous segregation distortion, a linkage map was calculated using 201 SNPs and 33 SSR markers with an average intermarker distance of 2.85 cM. The resulting linkage map comprised eight linkage groups and a total map distance of 666.1 centiMorgans (Fig 1). One SNP from scaffold_10 and two SNPs from scaffold_4 of the peach genome v1.0 mapped to linkage group 3 (Fig 1). Two SNPs from scaffold_16 of the peach genome v1.0 mapped to linkage group 2 (Fig 1). Placement of these five SNPs is in agreement with the recent corrected assembly of the peach genome [17]. No SNPs were identified in the first several million bp of linkage groups 5, 6, and 8 (Fig 1).

Download:

Fig 1. Linkage map of F₂ mapping population segregating for chilling requirement (CR) and bloom date (BD).

Linkage group numbering corresponds to the peach reference genetic map. Left side of bars are marker distances in cM. Labels on the right side of bars are the marker name. Markers with names followed by an asterix have significantly distorted genotypic ratios (*, p<0.05; **, p<0.01; ***, p<0.005; ****, p<0.001; *****, p<0.0005; ******, p<0.0001). Italicized labels are SSR markers, all other markers are SNPs and are named for the scaffold and bp position on the scaffold as found in the Prunus persica genome (v1.0). QTL intervals are shown to the left of each LG. Bars and whiskers on QTL represent the 1 LOD and 2 LOD confidence intervals from the QTL peak, respectively. Chilling requirement and bloom date QTLs are represented by solid and hatched bars, respectively. A QTL is named as qXXYa—ZZZZ, with ‘XX’ being the trait abbreviation, ‘Y’ the number of the linkage group, ‘a’ the letter to specify different QTLs for the same trait in one linkage group (G), and ‘ZZZZ’ the year in which the trait was phenotyped.

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

Significant segregation distortion was observed for markers in the lower half of linkage group 1 (Fig 1). Segregation distortion in this region was caused by underrepresentation of genotypes homozygous for the allele derived from the high chill female grandparent ‘Hakuho’ and overrepresentation of genotypes homozygous for the allele derived from the low chill male grandparent ‘UFGold’.

CR and BD phenotyping

F₂ individuals segregated for both chilling requirement (CR) and bloom date (BD). CR ranged from as little as 300–400 chilling hours to >1100 chilling hours (Fig 2). BD varied strongly by year (Fig 3). The interval between the earliest and latest blooming F₂ individuals varied with year from as little as 9 days (2010) to as long as 37 days (2007). Duration of the interval between earliest and latest bloom was reflective of year to year variation in chilling accumulation. Interruptions of chilling accumulation by warm periods promoted earlier bloom of low chill genotypes while consistent cold weather compressed bloom of all genotypes into a shorter interval (S2 and S3 Figs). Despite year to year variation in BD and earliest and latest BDs in the population, relative order of genotype BDs was highly correlated across years (Table 1). BD and CR phenotypes were also highly correlated (Table 1).

Download:

Fig 2. Histogram of chilling requirement (CR) phenotypes of F₂ progeny in 2008/2009 and 2009/2010.

Bins are 100 chill hour intervals. Letters indicate grandparental phenotypes (H, ‘Hakuho’; U, ‘UFGold’).

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

Download:

Fig 3. Histogram of bloom date phenotypes of F₂ progeny in seven seasons of observations.

Bloom date is measured from January 1^st. Letters indicate grandparental phenotypes (H, ‘Hakuho’; U, ‘UFGold’).

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

Download:

Table 1. Spearman’s rank order correlation coefficients for observed chilling requirement (CR) and bloom date (BD).

https://doi.org/10.1371/journal.pone.0139406.t001

QTL detection

Eight QTLs for CR were detected in 2008/2009 and two QTLs for CR were detected in the 2009/2010 (Fig 1, Table 2). Individual CR QTLs accounted for 4.0 to 27.8 percent of the phenotypic variation across years (Table 2). Nineteen QTLs were detected for BD across seven years of observations (Fig 1, Table 3). Individual BD QTLs accounted for 7.6 to 44.6 percent of phenotypic variation across years (Table 3). No CR or BD QTLs were identified on linkage groups 3 or 6 in any year.

Download:

Table 2. Chilling requirement QTL detected in 2008/2009 and 2009/2010.

https://doi.org/10.1371/journal.pone.0139406.t002

Download:

Table 3. Bloom date QTL detected by year from 2006–2012.

QTL are named as in Table 2.

https://doi.org/10.1371/journal.pone.0139406.t003

Discussion

GBS is an effective genotyping method for peach

We successfully employed the GBS method to detect SNPs in a peach F₂ mapping population. We employed the library preparation method as published with no modifications and were able to obtain a data set containing >1.5 million of ApeKI-associated reads per sample. TASSEL GBS pipeline 3.0 SNP calling appeared to be sensitive to low sequence depth in a species with greater heterozygosity than is typical of the cultivated cereals for which it was developed [20]. Nonetheless, after filtering all loci for read depth and genotype quality, mapping success with the remaining SNPs was greatly improved. Relatively stringent filtering of loci for minimum read depth, missing data and identifiable grandparental alleles reduced the number of SNP loci below that which has been typically reported in other species [24–26]. Relaxation of these filtering standards and the use of genotype imputation methods may increase the number of usable SNP loci obtained. The low background polymorphism of peach and the use of an F₂ population may have also decreased the number of SNP loci relative to studies on unrelated individuals or F₁ hybrids, which would likely possess greater allelic diversity and heterogeneity [24–26]. Finally, if increased SNP density were desired, reducing the number of individuals per sequencing lane would also increase the capture rate of SNPs in the population. We observed a nearly tenfold increase in detected SNPs when read depth was increased three-fold. GBS is, to some extent, a ‘tunable’ technique in which read depth can be adjusted to achieve a desired marker density.

Our ‘Hakuho’ × ‘UFGold’ population lacked detectable polymorphism in the distal ends of LGs 5 and 6 (Fig 1). These regions also lack observable polymorphisms between the grandparental genotypes (S1 Fig) and could represent either true monomorphic regions or regions with a very low density of ApeKI restriction sites. An in silico ApeKI digestion of the peach genome predicted that the greatest distance between two ApeKI sites in LG 5 and 6 would be approximately 75 and 25 kb, respectively. Since these two regions are several Mb in length, a complete absence of ApeKI sites is insufficient to explain the absence of SNPs. However, since the GBS procedure selects genomic fragments between 100 and 400 bp in length which are flanked by ApeKI sites, regions with more widely-spaced ApeKI sites would also appear to lack polymorphism.

Creation of merged or consensus linkage maps across independent populations is a valuable method for improving QTL localization. Creation of consensus maps relies on the existence of ‘anchor’ markers which have been genotyped across all populations being mapped. To assess the potential for an independent GBS experiment to find markers in common with SNPs from other research groups, we cross-referenced the physical locations of grandparental SNP loci against two published and publically-available peach SNP datasets. The first set of 6,557 SNPs (‘UCD’) was generated by GoldenGate^® technology to facilitate mapping of two F₂ families segregating for fruit quality characters [3]. A second set of 9,000 SNPs (‘IGA_SNP’) was generated by re-sequencing of community-selected genotypes and curated for genome-wide representation [4]. Our pool of 4,833 GBS grandparental SNPs shared 33 SNP loci in common with the UCD dataset and 102 SNP loci in common with the IGA_SNP dataset (S2 File). The low number and uneven genomic distribution of shared SNP markers suggests that GBS-derived SNPs from a single population may have limited utility as anchoring markers for the integration of independent linkage maps (S2 File). This suggests a continued need for established anchoring markers currently used by the peach community to facilitate comparative mapping [13]. Conversion of novel SNPs discovered in one mapping project to cleaved amplified polymorphic sequence (CAPS) markers would also be a rapid and cost effective method to generate anchoring markers between maps of interest.

Genetic control of CR and BD

We identified 29 QTLs for CR and BD in this F₂ population across the two years of CR phenotyping and seven years of BD observation. Fewer QTLs were observed for CR than for BD. This discrepancy is consistent with the role of CR as a major, but not exclusive, determinant of BD. Additionally, more frequent BD observations allowed finer discrimination between genotypes than did the <10 bins of chill hour bins into which F₂ genotypes were partitioned. 2010 and 2011 had notably fewer BD QTLs of all of the years (two and one, respectively). All other years had three or four QTL. The reduced number of QTL in these two years likely resulted from the relatively compact population variation in bloom in 2010 and 2011 (Fig 3). The temperatures in these two years was consistently cold through March and as a result most of the genotypes were saturated in their chilling accumulation. Most trees were therefore able to rapidly bloom with the occurrence of warm temperatures, reducing the differentiation between genotypes.

Across all years the three CR and BD QTLs with the greatest effect were those found on LG1, LG4, and LG7 (Tables 1 and 2). QTLs in the same regions have been identified in a number of peach mapping populations [8, 27, 28] and also appear to be conserved across other Prunus species [28–32].

Several QTLs regions overlap with locations of hypothesized gene candidates for the genetic control of CR and BD. One group of strong CR and BD QTLs at the end of LG1 corresponds to the genomic location of the peach DAM gene cluster responsible for the evg mutation [33, 34]. In three years, QTLs for BD were detected within 10 cM of the DAM gene cluster, suggesting that this genomic region may also contain other candidates (Fig 1). The CR QTL located on LG2 (qCR2-2008) spans a region of the genome that contains PpeMADS22, a peach homolog of ParSOC1, that has been implicated in control of CR for vegetative bud break in P. armeniaca [32, 35, 36].

Our ‘Hakuho’ × ‘UFGold’ population displayed a strong segregation distortion at the bottom of LG1 as previously observed in a study of the unrelated ‘Contender’ × ‘FLA92-2C’ population [8]. Identification of similar segregation distortion in an unrelated cross is consistent with a linkage between the low chill grandparent derived alleles and a locus that distorts segregation in this region. However, the seeds of both of these populations (‘Hakuho’ × ‘UFGold’ and ‘Contender’ × ‘FLA92-2C’) were harvested, stratified and germinated in the same season and location. This allows the possibility that environmental selection occurred at an early developmental stage, biasing the resulting populations in a similar direction.

Conclusions

GBS is a rapid and cost-competitive method for genotyping in peach. In our example, 96 individuals were multiplexed on a single lane of an Illumina HiSeq® 2000 machine to produce a dataset of several hundred SNPs for use in linkage mapping. Assuming a conservatively high sequencing cost of $2000 per lane, this is approximately $21 per sample excluding the initial expense of the adaptor oligonucleotides, which can be amortized over >200 sequencing runs. Although this cost is still prohibitive for large scale breeding programs that screen thousands of individuals annually, it is competitive for genotyping populations in trait mapping studies. As sequencing costs continue to decline, the cost competitiveness of this method will only improve.

Supporting Information

S1 Fig. Physical positions of SNP markers polymorphic between the ‘Hakuho’ and ‘UFGold’ grandparental genotypes.

Physical map of peach chromosomes (vertical bars) with horizontal marks indicating position of SNPs between grandparental genotypes.

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

(TIF)

S2 Fig. Cumulative chill hours below 7 degrees C in the seven winters of bloom date observation.

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

(TIF)

S3 Fig. Cumulative chill hours between 0–7 degrees C in the seven winters of bloom date observation.

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

(TIF)

S1 File. Scripts for extracting barcode sequences from raw fastq files.

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

(PDF)

S2 File. List of peach SNPs in common between this study, IGA_SNP, and UCD datasets.

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

(XLSX)

Acknowledgments

The authors thank Kathy Brock for the bloom date phenotyping, the many people who assisted with sample collection for chilling requirement phenotyping, and staff of the Clemson University Musser Fruit Research Center for their diligent care of trees.

Author Contributions

Conceived and designed the experiments: DGB AGA BR. Performed the experiments: DGB SF BR WRO. Analyzed the data: DGB BR CEW. Contributed reagents/materials/analysis tools: DGB AGA CEW KG GLR. Wrote the paper: DGB.

References

1. Rowe HC, Renaut S, Guggisberg A. RAD in the realm of next-generation sequencing technologies. Molecular Ecology. 2011;20(17):3499–502. pmid:WOS:000294224400004.
- View Article
- PubMed/NCBI
- Google Scholar
2. Elshire RJ, Glaubitz JC, Sun Q, Poland JA, Kawamoto K, Buckler ES, et al. A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity Species. PLoS ONE. 2011;6(5). pmid:WOS:000290224800026.
- View Article
- PubMed/NCBI
- Google Scholar
3. Martinez-Garcia PJ, Parfitt DE, Ogundiwin EA, Fass J, Chan HM, Ahmad R, et al. High density SNP mapping and QTL analysis for fruit quality characteristics in peach (Prunus persica L.). Tree Genet Genomes. 2013;9(1):19–36. pmid:WOS:000313816600002.
- View Article
- PubMed/NCBI
- Google Scholar
4. Verde I, Bassil N, Scalabrin S, Gilmore B, Lawley CT, Gasic K, et al. Development and Evaluation of a 9K SNP Array for Peach by Internationally Coordinated SNP Detection and Validation in Breeding Germplasm. PLoS ONE. 2012;7(4). pmid:WOS:000305339200093.
- View Article
- PubMed/NCBI
- Google Scholar
5. Dennis FG. Problems in standardizing methods for evaluating the chilling requirements for the breaking of dormancy in buds of woody plants. Hortscience. 2003;38(3):347–50. pmid:ISI:000183523700005.
- View Article
- PubMed/NCBI
- Google Scholar
6. Okie WR, Blackburn B. Increasing Chilling Reduces Heat Requirement for Floral Budbreak in Peach. Hortscience. 2011;46(2):245–52. pmid:WOS:000287336000019.
- View Article
- PubMed/NCBI
- Google Scholar
7. Luedeling E, Brown PH. A global analysis of the comparability of winter chill models for fruit and nut trees. International Journal of Biometeorology. 2011;55(3):411–21. pmid:WOS:000289561600013.
- View Article
- PubMed/NCBI
- Google Scholar
8. Fan S, Bielenberg DG, Zhebentyayeva TN, Reighard GL, Okie WR, Holland D, et al. Mapping quantitative trait loci associated with chilling requirement, heat requirement and bloom date in peach (Prunus persica). New Phytologist. 2010;185(4):917–30. pmid:20028471
- View Article
- PubMed/NCBI
- Google Scholar
9. Okie WR. Handbook of peach and nectarine varieties: performance in the southeastern United States and index of names (SuDoc A 1.76:714): United States Department of Agriculture, Agricultural Research Service; 1998.
10. Sherman WB, Lyrene PM. 'UFGold' peach. Fruit Var J. 1997;51(2):76–7. pmid:WOS:A1997XB95700003.
- View Article
- PubMed/NCBI
- Google Scholar
11. Weinberger JH. Chilling Requirements of Peach Varieties. Proceedings of the American Society for Horticultural Science. 1950;56:122–8.
- View Article
- Google Scholar
12. Dellaporta SL, Wood J, Hicks JB. A plant DNA minipreparation: version II. Plant Mol Biol Rep. 1983;1(4):19–22.
- View Article
- Google Scholar
13. Dirlewanger E, Graziano E, Joobeur T, Garriga-Caldere F, Cosson P, Howad W, et al. Comparative mapping and marker-assisted selection in Rosaceae fruit crops. Proc Natl Acad Sci U S A. 2004;101(26):9891–6. pmid:ISI:000222405600067.
- View Article
- PubMed/NCBI
- Google Scholar
14. Howad W, Yamamoto T, Dirlewanger E, Testolin R, Cosson P, Cipriani G, et al. Mapping with a few plants: Using selective mapping for microsatellite saturation of the Prunus reference map. Genetics. 2005;171(3):1305–9. pmid:WOS:000233967200036.
- View Article
- PubMed/NCBI
- Google Scholar
15. Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple Combinations of Lineage-Determining Transcription Factors Prime cis-Regulatory Elements Required for Macrophage and B Cell Identities. Mol Cell. 2010;38(4):576–89. pmid:WOS:000278448100012.
- View Article
- PubMed/NCBI
- Google Scholar
16. Glaubitz JC, Casstevens TM, Lu F, Harriman J, Elshire RJ, Sun Q, et al. TASSEL-GBS: A High Capacity Genotyping by Sequencing Analysis Pipeline. PLoS ONE. 2014;9(2). pmid:WOS:000332396200176.
- View Article
- PubMed/NCBI
- Google Scholar
17. Verde I, Abbott AG, Scalabrin S, Jung S, Shu SQ, Marroni F, et al. The high-quality draft genome of peach (Prunus persica) identifies unique patterns of genetic diversity, domestication and genome evolution. Nature Genetics. 2013;45(5):487–U47. pmid:WOS:000318158200009.
- View Article
- PubMed/NCBI
- Google Scholar
18. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–U54. pmid:WOS:000302218500017.
- View Article
- PubMed/NCBI
- Google Scholar
19. Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and VCFtools. Bioinformatics. 2011;27(15):2156–8. pmid:WOS:000292778700023.
- View Article
- PubMed/NCBI
- Google Scholar
20. Hyma KE, Barba P, Wang M, Londo JP, Acharya CB, Mitchell SE, et al. HetMappS: Heterozygous Mapping Strategy for High Resolution Genotyping-by-Sequencing Markers. PLoS ONE. 2015;10(8). Epub August 5, 2015. pmid:26244767
- View Article
- PubMed/NCBI
- Google Scholar
21. Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007;23(19):2633–5. pmid:WOS:000250673800021.
- View Article
- PubMed/NCBI
- Google Scholar
22. Van Ooijen JW. JoinMap® 4, Software for the caluculation of genetic linkage maps in experimental populations. Wageningen, Netherlands: Kyazma B.V.; 2006.
23. Van Ooijen JW. MapQTL ® 6, Software for the mapping of quantitative trait loci in experimental populations of diploid species.: Kyazma B.V., Wageningen, Netherlands.; 2009.
24. Barba P, Cadle-Davidson L, Harriman J, Glaubitz JC, Brooks S, Hyma K, et al. Grapevine powdery mildew resistance and susceptibility loci identified on a high-resolution SNP map. Theor Appl Genet. 2014;127(1):73–84. pmid:WOS:000330436500007.
- View Article
- PubMed/NCBI
- Google Scholar
25. Gardner KM, Brown P, Cooke TF, Cann S, Costa F, Bustamante C, et al. Fast and Cost-Effective Genetic Mapping in Apple Using Next-Generation Sequencing. G3-Genes Genomes Genetics. 2014;4(9):1681–7. pmid:WOS:000342570600013.
- View Article
- PubMed/NCBI
- Google Scholar
26. Ward JA, Bhangoo J, Fernandez-Fernandez F, Moore P, Swanson JD, Viola R, et al. Saturated linkage map construction in Rubus idaeus using genotyping by sequencing and genome-independent imputation. BMC Genomics. 2013;14. pmid:WOS:000315033800001.
- View Article
- PubMed/NCBI
- Google Scholar
27. Romeu JF, Monforte AJ, Sanchez G, Granell A, Garcia-Brunton J, Badenes ML, et al. Quantitative trait loci affecting reproductive phenology in peach. BMC Plant Biology. 2014;14. pmid:WOS:000335065200001.
- View Article
- PubMed/NCBI
- Google Scholar
28. Dirlewanger E, Quero-Garcia J, Le Dantec L, Lambert P, Ruiz D, Dondini L, et al. Comparison of the genetic determinism of two key phenological traits, flowering and maturity dates, in three Prunus species: peach, apricot and sweet cherry. Heredity. 2012;109(5):280–92. pmid:WOS:000310098100003.
- View Article
- PubMed/NCBI
- Google Scholar
29. Sanchez-Perez R, Dicenta F, Martinez-Gomez P. Inheritance of chilling and heat requirements for flowering in almond and QTL analysis. Tree Genet Genomes. 2012;8(2):379–89. pmid:WOS:000301846800014.
- View Article
- PubMed/NCBI
- Google Scholar
30. Castede S, Campoy JA, Quero-Garcia J, Le Dantec L, Lafargue M, Barreneche T, et al. Genetic determinism of phenological traits highly affected by climate change in Prunus avium: flowering date dissected into chilling and heat requirements. New Phytologist. 2014;202(2):703–15. pmid:WOS:000333060500038.
- View Article
- PubMed/NCBI
- Google Scholar
31. Campoy JA, Ruiz D, Egea J, Rees DJG, Celton JM, Martinez-Gomez P. Inheritance of Flowering Time in Apricot (Prunus armeniaca L.) and Analysis of Linked Quantitative Trait Loci (QTLs) using Simple Sequence Repeat (SSR) Markers. Plant Mol Biol Rep. 2011;29(2):404–10. pmid:WOS:000288712400014.
- View Article
- PubMed/NCBI
- Google Scholar
32. Olukolu BA, Trainin T, Fan SH, Kole C, Bielenberg DG, Reighard GL, et al. Genetic linkage mapping for molecular dissection of chilling requirement and budbreak in apricot (Prunus armeniaca L.). Genome. 2009;52(10):819–28. pmid:WOS:000271832300001.
- View Article
- PubMed/NCBI
- Google Scholar
33. Jiménez S, Reighard GL, Bielenberg DG. Gene expression of DAM5 and DAM6 is suppressed by chilling temperatures and inversely correlated with bud break rate. Plant Molecular Biology. 2010;73(1–2):157–67. pmid:20143130
- View Article
- PubMed/NCBI
- Google Scholar
34. Bielenberg DG, Wang Y, Li Z, Zhebentyayeva T, Fan S, Reighard GL, et al. Sequencing and annotation of the evergrowing locus in peach [Prunus persica (L.) Batsch] reveals a cluster of six MADS-box transcription factors as candidate genes for regulation of terminal bud formation. Tree Genetics and Genomes. 2008;4(3):495–507.
- View Article
- Google Scholar
35. Trainin T, Bar-Ya'akov I, Holland D. ParSOC1, a MADS-box gene closely related to Arabidopsis AGL20/SOC1, is expressed in apricot leaves in a diurnal manner and is linked with chilling requirements for dormancy break. Tree Genet Genomes. 2013;9(3):753–66. pmid:WOS:000318879000011.
- View Article
- PubMed/NCBI
- Google Scholar
36. Wells CE, Vendramin E, Tarodo SJ, Verde I, Bielenberg DG. A genome-wide analysis of MADS-box genes in peach Prunus persica (L.) Batsch. BMC Plant Biology. 2015;15. pmid:WOS:000349410600001.
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Rowe HC, Renaut S, Guggisberg A. RAD in the realm of next-generation sequencing technologies. Molecular Ecology. 2011;20(17):3499–502. pmid:WOS:000294224400004.
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Elshire RJ, Glaubitz JC, Sun Q, Poland JA, Kawamoto K, Buckler ES, et al. A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity Species. PLoS ONE. 2011;6(5). pmid:WOS:000290224800026.
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Martinez-Garcia PJ, Parfitt DE, Ogundiwin EA, Fass J, Chan HM, Ahmad R, et al. High density SNP mapping and QTL analysis for fruit quality characteristics in peach (Prunus persica L.). Tree Genet Genomes. 2013;9(1):19–36. pmid:WOS:000313816600002.
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Verde I, Bassil N, Scalabrin S, Gilmore B, Lawley CT, Gasic K, et al. Development and Evaluation of a 9K SNP Array for Peach by Internationally Coordinated SNP Detection and Validation in Breeding Germplasm. PLoS ONE. 2012;7(4). pmid:WOS:000305339200093.
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Dennis FG. Problems in standardizing methods for evaluating the chilling requirements for the breaking of dormancy in buds of woody plants. Hortscience. 2003;38(3):347–50. pmid:ISI:000183523700005.
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. Okie WR, Blackburn B. Increasing Chilling Reduces Heat Requirement for Floral Budbreak in Peach. Hortscience. 2011;46(2):245–52. pmid:WOS:000287336000019.
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref7] 7. Luedeling E, Brown PH. A global analysis of the comparability of winter chill models for fruit and nut trees. International Journal of Biometeorology. 2011;55(3):411–21. pmid:WOS:000289561600013.
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref8] 8. Fan S, Bielenberg DG, Zhebentyayeva TN, Reighard GL, Okie WR, Holland D, et al. Mapping quantitative trait loci associated with chilling requirement, heat requirement and bloom date in peach (Prunus persica). New Phytologist. 2010;185(4):917–30. pmid:20028471
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref9] 9. Okie WR. Handbook of peach and nectarine varieties: performance in the southeastern United States and index of names (SuDoc A 1.76:714): United States Department of Agriculture, Agricultural Research Service; 1998.

[ref10] 10. Sherman WB, Lyrene PM. 'UFGold' peach. Fruit Var J. 1997;51(2):76–7. pmid:WOS:A1997XB95700003.
View Article
PubMed/NCBI
Google Scholar

[35] View Article

[36] PubMed/NCBI

[37] Google Scholar

[ref11] 11. Weinberger JH. Chilling Requirements of Peach Varieties. Proceedings of the American Society for Horticultural Science. 1950;56:122–8.
View Article
Google Scholar

[39] View Article

[40] Google Scholar

[ref12] 12. Dellaporta SL, Wood J, Hicks JB. A plant DNA minipreparation: version II. Plant Mol Biol Rep. 1983;1(4):19–22.
View Article
Google Scholar

[42] View Article

[43] Google Scholar

[ref13] 13. Dirlewanger E, Graziano E, Joobeur T, Garriga-Caldere F, Cosson P, Howad W, et al. Comparative mapping and marker-assisted selection in Rosaceae fruit crops. Proc Natl Acad Sci U S A. 2004;101(26):9891–6. pmid:ISI:000222405600067.
View Article
PubMed/NCBI
Google Scholar

[45] View Article

[46] PubMed/NCBI

[47] Google Scholar

[ref14] 14. Howad W, Yamamoto T, Dirlewanger E, Testolin R, Cosson P, Cipriani G, et al. Mapping with a few plants: Using selective mapping for microsatellite saturation of the Prunus reference map. Genetics. 2005;171(3):1305–9. pmid:WOS:000233967200036.
View Article
PubMed/NCBI
Google Scholar

[49] View Article

[50] PubMed/NCBI

[51] Google Scholar

[ref15] 15. Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple Combinations of Lineage-Determining Transcription Factors Prime cis-Regulatory Elements Required for Macrophage and B Cell Identities. Mol Cell. 2010;38(4):576–89. pmid:WOS:000278448100012.
View Article
PubMed/NCBI
Google Scholar

[53] View Article

[54] PubMed/NCBI

[55] Google Scholar

[ref16] 16. Glaubitz JC, Casstevens TM, Lu F, Harriman J, Elshire RJ, Sun Q, et al. TASSEL-GBS: A High Capacity Genotyping by Sequencing Analysis Pipeline. PLoS ONE. 2014;9(2). pmid:WOS:000332396200176.
View Article
PubMed/NCBI
Google Scholar

[57] View Article

[58] PubMed/NCBI

[59] Google Scholar

[ref17] 17. Verde I, Abbott AG, Scalabrin S, Jung S, Shu SQ, Marroni F, et al. The high-quality draft genome of peach (Prunus persica) identifies unique patterns of genetic diversity, domestication and genome evolution. Nature Genetics. 2013;45(5):487–U47. pmid:WOS:000318158200009.
View Article
PubMed/NCBI
Google Scholar

[61] View Article

[62] PubMed/NCBI

[63] Google Scholar

[ref18] 18. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–U54. pmid:WOS:000302218500017.
View Article
PubMed/NCBI
Google Scholar

[65] View Article

[66] PubMed/NCBI

[67] Google Scholar

[ref19] 19. Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and VCFtools. Bioinformatics. 2011;27(15):2156–8. pmid:WOS:000292778700023.
View Article
PubMed/NCBI
Google Scholar

[69] View Article

[70] PubMed/NCBI

[71] Google Scholar

[ref20] 20. Hyma KE, Barba P, Wang M, Londo JP, Acharya CB, Mitchell SE, et al. HetMappS: Heterozygous Mapping Strategy for High Resolution Genotyping-by-Sequencing Markers. PLoS ONE. 2015;10(8). Epub August 5, 2015. pmid:26244767
View Article
PubMed/NCBI
Google Scholar

[73] View Article

[74] PubMed/NCBI

[75] Google Scholar

[ref21] 21. Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007;23(19):2633–5. pmid:WOS:000250673800021.
View Article
PubMed/NCBI
Google Scholar

[77] View Article

[78] PubMed/NCBI

[79] Google Scholar

[ref22] 22. Van Ooijen JW. JoinMap® 4, Software for the caluculation of genetic linkage maps in experimental populations. Wageningen, Netherlands: Kyazma B.V.; 2006.

[ref23] 23. Van Ooijen JW. MapQTL ® 6, Software for the mapping of quantitative trait loci in experimental populations of diploid species.: Kyazma B.V., Wageningen, Netherlands.; 2009.

[ref24] 24. Barba P, Cadle-Davidson L, Harriman J, Glaubitz JC, Brooks S, Hyma K, et al. Grapevine powdery mildew resistance and susceptibility loci identified on a high-resolution SNP map. Theor Appl Genet. 2014;127(1):73–84. pmid:WOS:000330436500007.
View Article
PubMed/NCBI
Google Scholar

[83] View Article

[84] PubMed/NCBI

[85] Google Scholar

[ref25] 25. Gardner KM, Brown P, Cooke TF, Cann S, Costa F, Bustamante C, et al. Fast and Cost-Effective Genetic Mapping in Apple Using Next-Generation Sequencing. G3-Genes Genomes Genetics. 2014;4(9):1681–7. pmid:WOS:000342570600013.
View Article
PubMed/NCBI
Google Scholar

[87] View Article

[88] PubMed/NCBI

[89] Google Scholar

[ref26] 26. Ward JA, Bhangoo J, Fernandez-Fernandez F, Moore P, Swanson JD, Viola R, et al. Saturated linkage map construction in Rubus idaeus using genotyping by sequencing and genome-independent imputation. BMC Genomics. 2013;14. pmid:WOS:000315033800001.
View Article
PubMed/NCBI
Google Scholar

[91] View Article

[92] PubMed/NCBI

[93] Google Scholar

[ref27] 27. Romeu JF, Monforte AJ, Sanchez G, Granell A, Garcia-Brunton J, Badenes ML, et al. Quantitative trait loci affecting reproductive phenology in peach. BMC Plant Biology. 2014;14. pmid:WOS:000335065200001.
View Article
PubMed/NCBI
Google Scholar

[95] View Article

[96] PubMed/NCBI

[97] Google Scholar

[ref28] 28. Dirlewanger E, Quero-Garcia J, Le Dantec L, Lambert P, Ruiz D, Dondini L, et al. Comparison of the genetic determinism of two key phenological traits, flowering and maturity dates, in three Prunus species: peach, apricot and sweet cherry. Heredity. 2012;109(5):280–92. pmid:WOS:000310098100003.
View Article
PubMed/NCBI
Google Scholar

[99] View Article

[100] PubMed/NCBI

[101] Google Scholar

[ref29] 29. Sanchez-Perez R, Dicenta F, Martinez-Gomez P. Inheritance of chilling and heat requirements for flowering in almond and QTL analysis. Tree Genet Genomes. 2012;8(2):379–89. pmid:WOS:000301846800014.
View Article
PubMed/NCBI
Google Scholar

[103] View Article

[104] PubMed/NCBI

[105] Google Scholar

[ref30] 30. Castede S, Campoy JA, Quero-Garcia J, Le Dantec L, Lafargue M, Barreneche T, et al. Genetic determinism of phenological traits highly affected by climate change in Prunus avium: flowering date dissected into chilling and heat requirements. New Phytologist. 2014;202(2):703–15. pmid:WOS:000333060500038.
View Article
PubMed/NCBI
Google Scholar

[107] View Article

[108] PubMed/NCBI

[109] Google Scholar

[ref31] 31. Campoy JA, Ruiz D, Egea J, Rees DJG, Celton JM, Martinez-Gomez P. Inheritance of Flowering Time in Apricot (Prunus armeniaca L.) and Analysis of Linked Quantitative Trait Loci (QTLs) using Simple Sequence Repeat (SSR) Markers. Plant Mol Biol Rep. 2011;29(2):404–10. pmid:WOS:000288712400014.
View Article
PubMed/NCBI
Google Scholar

[111] View Article

[112] PubMed/NCBI

[113] Google Scholar

[ref32] 32. Olukolu BA, Trainin T, Fan SH, Kole C, Bielenberg DG, Reighard GL, et al. Genetic linkage mapping for molecular dissection of chilling requirement and budbreak in apricot (Prunus armeniaca L.). Genome. 2009;52(10):819–28. pmid:WOS:000271832300001.
View Article
PubMed/NCBI
Google Scholar

[115] View Article

[116] PubMed/NCBI

[117] Google Scholar

[ref33] 33. Jiménez S, Reighard GL, Bielenberg DG. Gene expression of DAM5 and DAM6 is suppressed by chilling temperatures and inversely correlated with bud break rate. Plant Molecular Biology. 2010;73(1–2):157–67. pmid:20143130
View Article
PubMed/NCBI
Google Scholar

[119] View Article

[120] PubMed/NCBI

[121] Google Scholar

[ref34] 34. Bielenberg DG, Wang Y, Li Z, Zhebentyayeva T, Fan S, Reighard GL, et al. Sequencing and annotation of the evergrowing locus in peach [Prunus persica (L.) Batsch] reveals a cluster of six MADS-box transcription factors as candidate genes for regulation of terminal bud formation. Tree Genetics and Genomes. 2008;4(3):495–507.
View Article
Google Scholar

[123] View Article

[124] Google Scholar

[ref35] 35. Trainin T, Bar-Ya'akov I, Holland D. ParSOC1, a MADS-box gene closely related to Arabidopsis AGL20/SOC1, is expressed in apricot leaves in a diurnal manner and is linked with chilling requirements for dormancy break. Tree Genet Genomes. 2013;9(3):753–66. pmid:WOS:000318879000011.
View Article
PubMed/NCBI
Google Scholar

[126] View Article

[127] PubMed/NCBI

[128] Google Scholar

[ref36] 36. Wells CE, Vendramin E, Tarodo SJ, Verde I, Bielenberg DG. A genome-wide analysis of MADS-box genes in peach Prunus persica (L.) Batsch. BMC Plant Biology. 2015;15. pmid:WOS:000349410600001.
View Article
PubMed/NCBI
Google Scholar

[130] View Article

[131] PubMed/NCBI

[132] Google Scholar

Figures

Abstract

Introduction

Methods

Plant material

Phenotyping

Genotyping

Processing of sequenced reads

Genetic map construction and QTL discovery

Results

Sequencing and identification of SNPs

Linkage mapping

CR and BD phenotyping

QTL detection

Discussion

GBS is an effective genotyping method for peach

Genetic control of CR and BD

Conclusions

Supporting Information

S1 Fig. Physical positions of SNP markers polymorphic between the ‘Hakuho’ and ‘UFGold’ grandparental genotypes.

S2 Fig. Cumulative chill hours below 7 degrees C in the seven winters of bloom date observation.

S3 Fig. Cumulative chill hours between 0–7 degrees C in the seven winters of bloom date observation.

S1 File. Scripts for extracting barcode sequences from raw fastq files.

S2 File. List of peach SNPs in common between this study, IGA_SNP, and UCD datasets.

Acknowledgments

Author Contributions

References