Construction and Comparative Analyses of Highly Dense Linkage Maps of Two Sweet Cherry Intra-Specific Progenies of Commercial Cultivars

Despite the agronomical importance and high synteny with other Prunus species, breeding improvements for cherry have been slow compared to other temperate fruits, such as apple or peach. However, the recent release of the peach genome v1.0 by the International Peach Genome Initiative and the sequencing of cherry accessions to identify Single Nucleotide Polymorphisms (SNPs) provide an excellent basis for the advancement of cherry genetic and genomic studies. The availability of dense genetic linkage maps in phenotyped segregating progenies would be a valuable tool for breeders and geneticists. Using two sweet cherry (Prunus avium L.) intra-specific progenies derived from crosses between ‘Black Tartarian’ × ‘Kordia’ (BT×K) and ‘Regina’ × ‘Lapins’(R×L), high-density genetic maps of the four parental lines and the two segregating populations were constructed. For BT×K and R×L, 89 and 121 F1 plants were used for linkage mapping, respectively. A total of 5,696 SNP markers were tested in each progeny. As a result of these analyses, 723 and 687 markers were mapped into eight linkage groups (LGs) in BT×K and R×L, respectively. The resulting maps spanned 752.9 and 639.9 cM with an average distance of 1.1 and 0.9 cM between adjacent markers in BT×K and R×L, respectively. The maps displayed high synteny and co-linearity between each other, with the Prunus bin map, and with the peach genome v1.0 for all eight LGs (LG1–LG8). These maps provide a useful tool for investigating traits of interest in sweet cherry and represent a qualitative advance in the understanding of the cherry genome and its synteny with other members of the Rosaceae family.


Introduction
All cherry species belong to the Cerasus subgenus of the Prunus genus, within the Rosaceae family. Sweet cherry (Prunus avium L.) is an economically important crop that includes cherry trees cultivated for human consumption and wild cherry trees used for their wood, also called mazzards [1]; [2]. The majority of cultivated cherry trees belongs to the diploid (2n = 2x = 16) sweet cherry and allotetraploid (2n = 4x = 32) sour cherry (P. cerasus L.) species. Sweet cherry and the tetraploid ground cherry (P. fruticosa Pall., 2n = 4x = 32) are believed to be the parental species that gave rise to sour cherry [1]; [2]; [3]; [4].
Traditionally, the main sweet cherry breeding objectives have been to select for large, attractive and good-flavoured fruits, a short juvenile phase, abundant and consistent yields, reduced susceptibility to fruit cracking, self-compatibility and improved resistance or tolerance to diseases such as bacterial canker induced by Pseudomonas mors pv. prunorum and P. syringae [3]; [5]; [6]; [7]; [8].
Due to rapid climate change and the noteworthy reduction of chilling accumulation observed and/or predicted for several temperate zones [9], concern about climatic adaptation in temperate fruit trees has arisen [10]. Sweet cherries are cultivated commercially around the world in temperate, Mediterranean and even subtropical regions. In order to break bud dormancy, sweet cherries need to meet minimum chilling requirements in the autumn and winter [11]. Most fruit tree species are long-lived perennials, lasting more than 20 years, which implies that successful cultivars must be able to perform well despite changing climatic conditions. For this reason, climatic adaptation has become one of the major breeding objectives for many fruit crops and climatic adaptation of the parental plant material precedes breeding for commercial fruit qualities [12]. However, the long juvenile period (four to five years) in sweet cherry represents a drawback for quick and efficient breeding. Reducing the time needed to breed temperate fruit trees will be important to develop new cultivars that are able to face the challenges associated with temperate fruit production in a changing global environment. Marker-assisted breeding approaches based on genetic maps have the potential to save time and resources needed to select new sweet cherry cultivars.
The recent release of the peach genome v1.0 by the International Peach Genome Initiative [13]; [14], (GDR database: http://www.rosaceae.org/peach/genome) and the high synteny among Prunus species [15] will facilitate the characterization of genes responsible for agronomically important traits within this family. In cherry, maps from an inter-specific cross, P. avium 'Napoleon' 6 P. nipponica, have been developed [16]. From intraspecific crosses, Dirlewanger et al. [17] have developed a map from elite cultivars 'Regina' 6 'Lapins', whereas Olmstead et al. [18] have constructed a map from the great-grandparent of 'Lapins' (cultivar 'Emperor Francis'), and the wild forest cherry 'New York 549. More recently, Cabrera et al. [19] presented a consensus sweet cherry map constructed from four progenies using Rosaceae Conserved Orthologous Set (RosCOS) SNP and SSR markers. Despite the potential usefulness of genetic linkage maps for cherry, highly dense linkage maps for commercial sweet cherry cultivars have not yet been constructed. With the exception of the SNP-based map that we have recently reported [19], the majority of the cherry linkage maps have been based on low throughput molecular markers such as SSRs, AFLPs and isoenzymes. High throughput SNP genotyping technology has become available for some plant species. BeadChip or Infinium II technologies have been developed by consortiums or commercially for apple [20], peach [21], cherry, maize, tomato and Vitis [22], minimizing dramatically genotyping costs and facilitating the production of high density SNP maps. The RosBREED cherry 6K SNP array v1 for use with the Illumina InfiniumH system developed by the RosBREED consortium [23] for both diploid sweet cherry and allotetraploid sour cherry provides a new avenue for genome scanning capability in sweet and sour cherry. This publicly available genomics resource represents a significant advance for marker-locus-trait discovery and further research about the genome structure and diversity in this diploid-tetraploid crop set [23].
In this study, we develop highly dense linkage maps of two intraspecific progenies generated from the crosses of commercial cultivars: 'Black Tartarian'6'Kordia' (designated as BT6K) and 'Regina 6 Lapins' (designated as R6L), carried out in Chile and France, respectively, using the new RosBREED cherry 6K SNP array v1. To our knowledge, these are the first genetic maps using this new SNP genotyping resource developed in sweet cherry. The maps were compared with the Prunus bin map [24] and the peach genome v1.0 to assess synteny and co linearity between the studied progenies and the bin map and the peach genome v1.0. These maps provide a new framework for sweet cherry genomic information for the international community and new opportunities to discover QTLs associated with agronomical traits in this and related species, based upon the haplotype segregation in the progeny.
The field studies performed did not involve endangered or protected species, but cultivated sweet cherry trees. In addition, all the trees used in this study belong to public institutes and do not need any specific permit to be used for the analyses described in this paper.

DNA Extraction
Leaf material was frozen in liquid nitrogen and stored at 280uC for later use. Genomic DNA was extracted from the frozen tissue using the DNeasyH plant kit (Qiagen) according to the manufacturer's instructions. Genomic DNA was quantified using Accublue TM dsDNA quantification kit (BIOTIUM) according to the manufacturer's instructions for the BT6K progeny, whereas for the R6L progeny, genomic DNA was quantified using spectrophotometry Nanoview (GE Healthcare) and fluorimetry Quant-iT TM PicogreenH (Invitrogen) according to the manufacturer's instructions. Fifteen ml of DNA with a concentration between 50 ng/ml -75 ng/ml was used for subsequent analyses.

Genetic Marker Analysis
SNPs. Both sweet cherry mapping populations and their respective parents were genotyped using the RosBREED cherry 6K SNP array v1 that contains 5,696 SNPs of which 76% and 24% of the SNPs were chosen to target the sweet and sour cherry genomes, respectively [23]. SNPs were obtained from re-sequencing of a detection panel containing 16 sweet and eight sour cherry accessions. The accessions were founders, intermediate ancestors or important parents used in U.S. breeding programs [23]. SNPs were detected using SoapSNP [28] (SoapSNP: http://soap. genomics.org.cn/soapsnp.html) and validated by GoldenGateH assay [29]. Genotype differences were recorded in the iSCAN platform and SNP genotypes were determined using Genome Studio Genotyping Module (Version 1.8.4, Illumina, [30]; [31]; [32]) as described in Peace et al. [23]. All DNA samples were above the GenCall threshold of 0.15 and were therefore all used in further analyses [33]. Initial clustering was done using Gentrain2, a GenomeStudio build-in clustering algorithm ( [30]; [31]; [32]). Following the clustering by Gentrain2, all SNPs were visually examined for appropriateness of clustering, cluster separation, number of clusters, presence of null alleles and paralogs along with a check on whether the progenies follow expected genetic ratios and presence of genotyping errors. A SNP was considered ''failed'' if that have overlapping clusters or ambiguous clusters which cannot be improved by even manual clustering, more than 3 clusters suggesting presence of paralogs and very low call frequency, and positive outlier in terms of Mendelian inheritance error rates (e.g. parent-child or parent-parent-child) [33]. The failed SNPs were not used in the mapping process.
The RosBREED cherry 6K SNP array v1 markers used in this work were deposited in NCBI's dbSNP repository available at www.ncbi.nlm.nih.gov/projects/SNP [34] and each SNP was given a unique accession number that starts with the prefix ss (SNPs NCBI ss# database names). More information associated with these SNPs is available at the Genome Database for Rosaceae (GDR; www.rosaceae.org [35]). The location of the markers on the physical map of the peach genome v1.0 as well as their positions on the genetic maps of BT, K, BT6K, L, R and R6L are available in Table S1.
The MAF (minor allele frequency) distribution of SNPs heterozygous in each of the four mapping parents ('Black Tartarian', 'Kordia', 'Regina' and 'Lapins') was determined and then compared with a germplasm collection of 269 sweet cherry selections used by Peace et al. [23]. Similarly, the distribution and physical spacing (Mbp) of SNP heterozygocity along the eight cherry chromosomes for the four linkage mapping parents was compared with polymorphic SNP markers identified in the germplasm collection of the 269 sweet cherry selections [23].
Firstly, parental maps were constructed using only heterozygous markers for one parent (class ,nn 6 np. and ,lm 6 ll., less ambiguous in tracing origin of allele than ,hk 6 hk. class). 'Suspect Linkage' (recombination frequency .0.6) and 'Genotype Probabilities' (-Log10 (P) .2.0) tools were used to identify misgrouped markers and double recombination events, respectively. Spurious double crossovers were considered as missing value while revising the map. Chromosome painting was performed for each linkage group of the four parents to verify the double recombination events. Double recombination events in narrow genetic distances were omitted from subsequent mapping analyses. The x 2 (chi-square) goodness-of-fit test was used to evaluate the discrepancy with the expected segregation ratios. In a first step, all segregating markers were used for initial mapping. In a second step, isolated markers showing segregation distortions different from expected Mendelian ratios, with a probability lower than p = 0.01 were excluded from further analyses. The grouping and ordering of markers was done using JoinMapH 4.0's build algorithm for cross pollinated species, a likelihood based weighted least square procedure (with the squares of the LODs as weights) as described by Stam [38]. The inspection of proper assignment of a marker to a group was done by calculating the Strongest Cross Link (SCL) parameter. The mapping procedure implemented in JoinMap consisted of iterative process: (a) adding loci one by one starting from most informative and each added marker locus's best position is searched by comparing the goodness-of-fit of the calculated map for each tested position (b) when the goodness of fit measure decreases sharply (i.e. difference known as jump) or gives negative estimates in the map the locus is removed and process is continued till first round is completed [38]. Second and third round accommodate the marker loci that are previously removed without the constraints of maximum allowed jump level or negative values. Maps were calculated iterating different LOD values in order to minimize the Mean Chi-square contributions. LOD thresholds of six and three were used for grouping and mapping, respectively. Kosambi's mapping function [39] was used to convert recombination frequency into map distance. Although peach physical map based on peach physical genome v1.0 provided a reference for comparisons of the linkage map output, no fixed order was forced based on peach physical positions. MapChart [40] was used to draw the linkage maps and make group-wise comparisons between maps. Secondly, an integrated map of each progeny was undertaken including heterozygous markers in both parents (,hk 6 hk.) to decrease large gaps in parental LGs. New parental maps were developed including (,hk 6 hk.) markers as previously described and integrated using the ''Combine Groups for Map Integration'' function of JoinMapH 4.0 to produce the combined maps.
The MAF distribution of SNPs heterozygous in each of the four mapping parents was similar to the distribution observed within a germplasm collection of 269 sweet cherry selections used by Peace et al [23] (Figure 1). A similar proportion of SNPs for each of the MAF categories were observed. The majority of the heterozygous Table 1. Origin and characteristics of the four parental cultivars used to perform the two mapping progenies.

Cultivar Parental cultivars Country of origin Main characters Reference
Black Tartarian U 6 U United Kingdom Medium chill, early flowering, self-incompatible (S1S2), low firmness, low sugar content, small fruit, acid fruit.
Intra-specific segregating population linkage maps. To reduce gaps in parental linkage maps, two linkage maps were constructed for each intra-specific segregating population including heterozygous markers in both parents (,hk6hk. markers). Each map consisted of eight Prunus LGs (Figure 2a, 2b; Figure S4). The maps covered 752.9 and 639.9 cM, containing 723 and 687  markers, respectively (Table 3). Each LG contained multiple markers, ranging from 55 to 133 markers. The average cM distance between markers was 1.1 for BT6K and 0.9 cM for R6L.
LG1 was the longest group in both maps covering 171.9 cM with 133 markers in BT6K, and 136.2 cM with 98 markers in R6L.
LG5 was the shortest group in both maps, covering 67.4 cM with 95 markers in BT6K, and 56.8 cM with 91 markers in R6L. Maximum gap size in each chromosome ranged from 4.2 to 12.3 (Table 3). In BT6K the largest gap (12 cM) was found in LG5 between the ss490554476 and ss490554529 markers; whereas in R6L it was in LG1 (12.3 cM), between the markers ss490547111 and ss490547271. A total of 3.8% (28/723) and 3.9% (27/687) skewed markers (chi-square probability p,0.01) were mapped in BT6K and R6L, respectively. Similarly to the parental linkage maps, clusters of skewed markers were found on LG1 (16/28 in BT6K), LG2 (5/27 in R6L), LG3 (12/28 in BT6K) and on LG6 (20/27 in R6L) (Figure 2a, 2b; Figure S4). This last region on LG6 corresponds to the S locus location in Prunus [14]. In the R6L progeny, 'Regina' genotype for the S alleles being S 1 S 3 and 'Lapins' being S 1 S 4 (Table 1), all the individuals inherited the S 4' from 'Lapins' as the S 1 is incompatible in an S 1 style (Figure 3). Using the markers framing the S locus described in Peace et al. ( [23]), the segregation of a 6 cM segment containing the S locus was followed in the two mapping progenies. For three individuals, a recombination was detected between the S 1 , S 3 alleles from 'Regina'. Within the BT6K progeny, 92.1% (82/89) individuals were characterized by non-recombinant haplotypes including the four combinations: S 1 S 6 (25/89), S 1 S 3 (24/89), S 2 S 6 (19/89) and S 2 S 3 (14/89). The other seven individuals were detected with one recombination event only in the 'Black Tartarian' parent ( Figure 3).
Comparative analyses among the linkage maps. As parental maps were constructed with markers heterozygous for only one parent of the cross (i.e. either ,nn 6np. or ,lm 6ll.), markers were in common between just one parent of each cross and the two parents of the other cross. The comparison of the parental maps was based upon the position of 277 markers: 88 common makers between 'Black Tartarian' and 'Regina', 74 between 'Black Tartarian' and 'Lapins', 87 between 'Kordia' and 'Regina', and 28 between 'Kordia' and 'Lapins' (Figure S3). A perfect co-linearity of the common markers between the parental linkage maps compared was found.
Comparison between BT6K and R6L linkage maps was performed using a total 379 common markers, representing 52.4% of the markers mapped in BT6K and 55.2% of the markers mapped in R6L ( Figure S4 Table 3. Number of markers, and map size and density of each parent map and of the two consensus maps (BT6K) and (R6L). LG1 LG2 LG3 LG4 LG5 LG6 LG7 LG8 Total conserved between the two sweet cherry linkage maps constructed and only small rearranged orders were found ( Figure S4).
Comparison between the sweet cherry genetic linkage maps and the physical peach map. Anchored markers between the intra-specific sweet cherry maps (BT6K and R6L) and the peach genome revealed that the relative position of 91.8% (944/1028) SNP markers are conserved between sweet cherry and peach (Table S1), demonstrating a high level of synteny between these two species. Only a small number of markers differed in their relative position when compared to the peach physical map ( Table 4). The SNPs in this table have been ordered cronologically for the position (cM) of the SNP markers on the sweet cherry parental maps. The majority of SNP markers with the LGs are oriented such that the number increases, similar to the orientation seen on the peach physical map. The regions represented in this table have differences between the orientation of the SNP markers on the LG, when compared to the peach physical map. Black boxes represent regions of the LG outside of the region in question. Green boxes represent regions in which the orientation of the SNP markers are inverted in relationship to the rest of the LG as well as the peach physical map. Marker position (cM) of the parental maps compared to the physical peach position v1.0, revealed a small percentage of local discrepancies (8.2%). There was only partial co-linearity of markers within the distal regions of LG1 and LG4; as well as the proximal regions of LG5, LG6 and LG7 of sweet cherry and peach (Table 4; Figure S3; Table S1). In region ,42.3 Mb to ,46.6 Mb of LG1, a total of 45 markers differed in their relative position on the parental maps when compared to the peach genome ('Black Tartarian' = 25, 'Kordia' = 3, 'Regina' = 12, and 'Lapins' = 10), five of these markers demonstrate a conservation of their relative position in at least two sweet cherry parental maps. In the region between ,26.7 Mb and ,26.8 Mb on LG4, there were local discrepancies in two markers on the 'Regina' map and one local discrepancy on the 'Kordia' map. The relative distance of the local discrepancy detected in 'Kordia' is only 0.7 cM relative to the nearest marker. This difference may be due to imprecisions associated with the small population size used in creating this map (89 individuals). However, the marker was conserved between 'Kordia' and 'Regina' maps. It was not possible to confirm this discrepancy in the other two parental maps, due to a limited number of polymorphic markers in this region in the other two cultivars. Similarly, on LG5, a discrepancie were observed between ,0.7 Mb and ,2.7 Mb for eight and seven markers in 'Black Tartarian' and 'Lapins', respectively. Three markers also showed discrepancies in 'Kordia' between ,3.5 Mb to ,4.8 Mb. On LG6, a discrepancy between ,0.0 Mb and ,1.3 Mb was observed only in the 'Regina' map (one marker). Other discrepancy in the region between ,1.3 Mb and ,4.3 Mb in the same linkage group was found for 13 and five markers in 'Kordia' and 'Lapins', respectively. On LG7, in a region from ,0.0 Mb to ,4.5 Mb, three markers in both 'Regina' and 'Lapins' maps were placed in an ''inverted order'' when compared with the peach physical map. However, in 'Black Tartarian' only two of the three markers were found having discrepancy with the peach genome. This ''inverted order'' was not found in 'Kordia' map. Other region on LG7 has discrepancy in two of the maps ('Black Tartarian' and 'Kordia' maps), from ,5.5 Mb to ,6.6 Mb for two markers in each parent (Table 4; Figure S3).
Only partial co-linearity of markers within the proximal region on LG5, LG6 and LG7; and distal region on LG1 and LG4 of intra specific sweet cherry maps (BT6K and R6L) and peach were detected, similar to what we observed in the parental maps. On the other hand, 1% (10/1.028) of the SNP markers were located on different LGs in the BT6K and R6L maps than predicted based upon the peach genome (v1.0) [14] (Table 5). A comparison between the position of these ten markers in the peach genome and the sweet cherry high density linkage maps revealed a clustering of four of the SNP markers in two narrow zones on LG8 in both progenies, situated at 45.0 cM and 62.6 cM in BT6K and at 44.5 cM and ,58.4 cM in R6L (Table 5).
Syntenic regions between the sweet cherry high density linkage maps and the Prunus bin map [24] were determined by comparing the location of RosCOS markers on the maps ( Figure S5). The relative positions of the markers in the BT6K and R6L intra specific sweet cherry maps (55 and 51 RosCOS markers, respectively) and the Prunus bin map are highly conserved, with the exception of one marker which is differential positions in both the BT6K and R6L maps when compared to the Prunus bin map ( Figure S5, marker ss490559286 located on LG6). Interestingly, this marker is located in the region of LG6 in which there are discrepancies between the 'Kordia' map and the peach genome.
There is also an almost perfect conservation of the relative map positions of these RosCOS markers when comparing the BT6K and R6L intra specific sweet cherry maps with the consensus map reported by Cabrera [19]. Of a total of 39 and 45 RosCOS markers conserved between the map reported by Cabrera [19] and the BT6K and R6L intra specific sweet cherry maps, respectively, there is only one RosCOS marker (RosCOS1551) located on LG5 whose relative position is not conserved in the R6L map.

Discussion
The number of heterozygous markers observed using the RosBREED cherry 6K SNP array v1 in each of the four parents (515-634, Table 2) represent 9-11% of all SNPs on the array (28-35% of 1,825 SNPs reported to be polymorphic among the 269 sweet cherry accessions evaluated by Peace et al. [23]). These results are consistent with the heterozygosity estimation of 400-700 markers for any given sweet cherry cultivar suggested by Peace et al. [23]. The majority of these heterozygous SNPs were SNPs that would more likely be heterozygous given that the MAFs for the majority of these SNPs were .0.2. Additionally, these results sustain the choice of the accessions used for the detection panel in the RosBREED cherry 6K SNP array v1 to efficiently represent cherry breeding germplasm from different origins ( Table 1) Scalabrin, personal communication) as it was already reported for peach [21]. These regions correspond to low recombination frequency as shown in Figure S1; S2.
A linkage map based on an intra-specific cross between 'Regina' and 'Lapins' has been previously reported [17]. Using the progeny of this intra-specific cross, as well as a newly developed intraspecific segregating population produced from a cross between 'Black Tartarian' 6 'Kordia', we have constructed two highdensity sweet cherry linkage maps. Over 650 SNPs have been placed on these high-density linkage maps, using the RosBREED cherry 6K SNP array v1. The vast majority of these SNP markers have not previously been positioned on sweet cherry genetic linkage maps. Additionally, analyses of the SNP markers in the intra-specific segregating population enabled us to develop linkage maps specific for the parental lines used in these crosses. These results demonstrate that the RosBREED cherry 6K SNP array v1 may be used to efficiently conduct genome-wide genotyping of sweet cherry cultivars and their progeny.
The length of four parental sweet cherry cultivars (719.4 cM 'Black Tartarian', 788 cM 'Kordia', 619.4 cM 'Regina' and 610.1 cM 'Lapins') are similar to those of other sweet cherry cultivars that have been published previously: 799.4 cM [19]; 711 cM in 'EF' and 565.8 cM in 'NY' [18]. Similarly, the length of the BT6K and R6L maps (752.9 cM and 639.9 cM respectively) were similar to that seen in the interspecific 'Napoleon (P. avium) x P. nipponica consensus map (680 cM) [16] as well as the sweet cherry consensus linkage map developed by Cabrera et al. [19] (779.4 cM). The genetic distance determined for LG1, LG4, LG5, LG7 and LG8 are equivalent to those reported in the Prunus bin map ( Figure S5). However, LG2, LG3 and LG6 are longer in BT6K and R6L when compared to the Prunus bin map.
A high co-linearity between BT6K and R6L maps was observed and only minor rearrangements of markers in small regions of the map were detected. Due to the relatively small sizes of the progeny population, minor variation in marker order between the BT6K and R6L high density linkage maps may be due to the lack of precision in the positioning of markers located less than 2 cM from each other. Alternatively, these minor variations between the maps may be due to imprecisions in the mapping of SNP markers that are heterozygous in both parents (,hk 6 hk.). The distance between the flanking markers and the heterozygous marker may alter the exact positioning of the heterozygous marker.
Among the large number of SNPs located in the high density sweet cherry linkage maps, the segregation of a small number of markers were skewed (28/723 BT6K and 27/687 R6L).  Interestingly, most of these markers showed skewed inheritance in both segregating populations and mapped to similar regions of LG1 and LG6. The clustering of the skewed markers suggests that there may be zones within the LGs that contain deleterious genes. The cluster on the top of LG6 (,4 Mb in peach v1.0 genome) in both mapping populations coincides with the position of the peach locus for male sterility [15]; [41]. In R6L, there was also a clustering of skewed markers near the bottom of LG6 (20 Mb in peach v1.0 genome) a zone that coincides with the position of the Prunus self-incompatibility (S) locus [15]. The S locus genotypes for 'Regina' and 'Lapins' are S 1 S 3 and S 1 S 4' , respectively, therefore, progeny of this cross contain only S 1 S 4' and S 3 S 4' allele combinations. Pollen with the S 1 genotype cannot germinate and fertilization cannot occur [42], resulting in the absence of S 1 S 1 and S 1 S 3 allele combinations in the progeny. This explains the skewed segregation of the markers surrounding the S locus located near the distal end of LG6. By analyzing haplotypes in the region of the S locus, it was possible to predict the S-allele combination for all individuals of the two progenies ( Figure 3). This is an example of application offered by the use of this cherry 6K SNP array that will be useful in breeding programs. The high co-linearity of the maps between the two non-related intra-specific progenies increases the likelihood that the marker order conserved in these maps is indicative of the sweet cherry genome organization. Additionally, the resulting maps (four parents and two progenies) showed high co-linearity with the Prunus bin map ( Figure S5). Previous genomic comparison studies have shown a high synteny level within the Prunus genus [15]. As expected, we found a high synteny between the four sweet cherry parents studied and the peach genome (v1.0), supporting the peach can be use as the model species for other Prunus members [13]; [43]; [44]. In addition; the reported discrepancies found on LG1, LG6 and LG7-from ,0.0 Mb to ,4.5 Mb-coincide with minor assembly errors which have been detected by the International  LG4, LG5 and part of LG7. The divergences could be attributable to mapping inaccuracy due to the low number of individuals analysed. Moreover, each of these divergences have only been found on two of four parental maps. Thus, our results confirm not only the high synteny of peach and sweet cherry, but also the high reliability of the peach genome sequence for comparative studies within the Prunus species [13]; [43]; [44].
In addition to the discrepancies within the linkage-group mentioned above, a few markers were not were not located on the same linkage group in both sweet cherry and peach. One marker in BT6K and six markers in the R6L maps mapped to different LGs then predicted based upon the peach genome. In addition, three markers, ss490555352, ss490556354 and ss490550342 were mapped to LG8 in both sweet cherry high density maps, whereas these markers are located on chromosome 6 of the peach genome. These four markers are located at ,45 cM and ,60 cM on LG8, suggesting that this region may have undergone a duplication event during the divergence between peach and sweet cherry. Sweet and sour cherry are considered the cultivated Prunus fruit species most distant from peach, as predicted by their divergent evolutionary origin [43]. Further studies would be required in order to verify the local discrepancies between the sweet cherry and peach genomes. This is the first report of a linkage map using the RosBREED cherry 6K SNP array v1. As the four cultivars used in this work derive from different origins, the polymorphism observed for their progenies from genome-scanning with the RosBREED cherry 6K SNP array v1 indicates its utility for sweet cherry breeding programs with diverse germplasms. The high proportion of monomorphic SNP markers (Table 2) agrees with the monomorphic rate (58%) obtained in the 6K SNP array. The monomorphic rate of this array was higher than those found for the recently published SNP arrays on peach [21], and apple [20] probably due to the lower sequencing depth used in the 6K SNP array [23]. However, in this study we mapped ,700 SNPs in each of the unrelated sweet cherry progenies examined.
The large number of common markers mapped in both progenies (379 markers) provides an excellent framework for comparative QTL discovery in sweet cherry. In conclusion, we have constructed two high-density sweet cherry linkage maps from two progenies using the RosBREED cherry 6K SNP array v1 which are highly co-linear with the previously published consensus sweet cherry linkage map, the Prunus bin map and the peach genome. The high density of heterozygous markers on these maps has the potential to not only enable high-resolution identification of QTLs in sweet cherry, but also to predict the genotypes at a specific locus linked to an agronomical interesting characters, as illustrated for the self-incompatibility S locus.  Figure S4 Comparison of the two consensus sweet cherry highly dense linkage maps of two intraspecific progenies (BT6K and R6L). Anchored markers are indicated by connecting lines and are represented in green. Markers in black are unique to each map. Distance between markers is represented in cM. Markers grouped in a different LG in comparison with peach genome v1.0 are underlined. Skewed markers mapped are represented by asterisks to indicate distortion level (* for p,0.1; ** p,0.05; *** p,0.01; **** p,0.005; ***** p,0.001; ****** p,0.0005; ******* p,0.0001). (TIF) Figure S5 Comparison between the sweet cherry highly dense linkage maps (BT6K and R6L) and the Prunus bin map. Positions of anchor loci between maps are indicated by connecting lines and are represented in green. Distance between markers is represented in cM. Markers grouped in a different LG in comparison with peach genome v1.0 are underlined. Skewed markers mapped are represented by asterisks to indicate distortion level (* for p,0.1; ** p,0.05; *** p,0.01; **** p,0.005; ***** p,0.001; ****** p,0.0005; ******* p,0.0001). Bin map markers are situated at the bottom of their corresponding bin. (TIF)