Genomic rearrangements and signatures of breeding in the allo-octoploid strawberry as revealed through an allele dose based SSR linkage map
© van Dijk et al.; licensee BioMed Central Ltd. 2014
Received: 19 December 2013
Accepted: 18 February 2014
Published: 1 March 2014
Breeders in the allo-octoploid strawberry currently make little use of molecular marker tools. As a first step of a QTL discovery project on fruit quality traits and resistance to soil-borne pathogens such as Phytophthora cactorum and Verticillium we built a genome-wide SSR linkage map for the cross Holiday x Korona. We used the previously published MADCE method to obtain full haplotype information for both of the parental cultivars, facilitating in-depth studies on their genomic organisation.
The linkage map incorporates 508 segregating loci and represents each of the 28 chromosome pairs of octoploid strawberry, spanning an estimated length of 2050 cM. The sub-genomes are denoted according to their sequence divergence from F. vesca as revealed by marker performance. The map revealed high overall synteny between the sub-genomes, but also revealed two large inversions on LG2C and LG2D, of which the latter was confirmed using a separate mapping population. We discovered interesting breeding features within the parental cultivars by in-depth analysis of our haplotype data. The linkage map-derived homozygosity level of Holiday was similar to the pedigree-derived inbreeding level (33% and 29%, respectively). For Korona we found that the observed homozygosity level was over three times higher than expected from the pedigree (13% versus 3.6%). This could indicate selection pressure on genes that have favourable effects in homozygous states. The level of kinship between Holiday and Korona derived from our linkage map was 2.5 times higher than the pedigree-derived value. This large difference could be evidence of selection pressure enacted by strawberry breeders towards specific haplotypes.
The obtained SSR linkage map provides a good base for QTL discovery. It also provides the first biologically relevant basis for the discernment and notation of sub-genomes. For the first time, we revealed genomic rearrangements that were verified in a separate mapping population. We believe that haplotype information will become increasingly important in identifying marker-trait relationships and regions that are under selection pressure within breeding material. Our attempt at providing a biological basis for the discernment of sub-genomes warrants follow-up studies to streamline the naming of the sub-genomes among different octoploid strawberry maps.
KeywordsGenomic rearrangement Inversion Fragaria Polyploid Haplotype Homozygosity MAS Selection Breeding signature
Cultivated strawberry (Fragaria x ananassa) is an important soft fruit species that is grown worldwide. Strawberry is a vegetatively propagated outbred species derived from the hybridisation of two new world species (Fragaria chiloensis and Fragaria virginiana) in the 18th century . As a member of the Rosacaea family, it shares ancestry with a variety of important food and ornamental crops such as apple, pear, peach and rose. Despite its economic importance and membership in a well-studied family, strawberry breeding to date rarely incorporates the use of molecular marker resources due to its complex, allo-octoploid genetic composition . Because of this complexity, there are only a limited number of studies where clear marker-trait relationships for major genes/QTLs were identified ([3–8]).
The first comprehensive molecular genetic maps in strawberry were developed for the diploid wild species Fragaria vesca[9–12]. This effort culminated in the completion of the draft sequence of the diploid Fragaria vesca clone ‘Hawai 4’ in late 2010 , which provided the rosaceous community a highly valuable tool for further genomic research.
Soon after the first genetic map of diploid strawberry was published, similar studies were initiated for the octoploid strawberry, resulting in the completion of several (partial) genetic maps [6, 14–22]. These mapping studies also conclusively revealed that the octoploid strawberry showed genome-wide disomic inheritance [1, 2] and could therefore be classified as a full allopolyploid.
The origins of the homoeologues (or sub-genomes) in allo-octoploid strawberry have not been studied as extensively as those of other allopolyploid crops such as bread wheat and Cotton [23, 24]. Molecular genetic studies revealed that the chloroplast DNA of octoploid strawberry is most closely related to that of the diploid species Fragaria vesca (subsp. bracteata) [25, 26]. In another study on nuclear genes, it was confirmed that part of the genome was clearly related to Fragaria vesca, and another part was related to the wild diploid Fragaria innumae, leading to the hypothesis that the octoploid genome originated from the fusing of unreduced gametes of two tetraploid species from quite distinct genetic backgrounds. To date, no convention exists for naming homoeologues, and none of the octoploid mapping studies have incorporated information on the origins of the different homoelogues in the naming of the linkage groups. For this reason, the assignment of homoelogue letters is not consistent between the different octoploid maps.
The need to obtain complete haplotype information from microsatellites utilised in polyploids motivated van Dijk et al.  to develop the MADCE (Microsatellite Allele Dose Configuration & Establishment) methodology for determining the allelic configuration of allopolyploid plant species . This method essentially converts any allopolyploid genome into a diploid genome regarding the software and methodologies that can be employed for genetic analysis.
In this study, we created a highly comprehensive genetic SSR linkage map of the octoploid strawberry using MADCE. We used this map to differentiate homoeologues based on their efficiency in amplifying F. vesca-derived markers and to discover genomic rearrangements among the diploid sub-genomes (homoeologues). This map provides the genetic makeup of the two parental varieties and their levels of homozygosity and haplotype sharing. Finally, we made comparisons of the cultivated strawberry genetic map to the physical reference map of the wild diploid F. vesca.
Genomic DNA was extracted according to a modified version of the Fulton et al.  mini-prep protocol. Briefly, 1 g of young, folded leaves were harvested. The leaves were freeze-dried and ground to powder in a 2 ml tube. To this tube, 700 μL of warm (65°C) containing 2% CTAB buffer was added, and the contents were mixed by vortexing and incubated for 10 min. Next, 700 μL of chloroform:isoamyl alcohol (24:1) was added. The mixture was centrifuged at room temperature at 10,000 g for 2 min. Next, 600 μL of the top phase was transferred to a fresh tube. Isopropanol (480 μL) was added, and the sample was mixed and then centrifuged at 10,000 g for 2 min at room temperature. The supernatant was discarded, and the pellet was washed with 500 μL of 70% ethanol, left for 2 min and then centrifuged at 10,000 g for 2 min. The supernatant was discarded by pipetting, and the pellet was resuspended in 400 μL of Tris EDTA. LiCl (135 μL, 8 M) was added to remove RNA and polysaccharides, and the mixture was incubated for 30 min at −20°C. After incubation, the mixture was centrifuged at room temperature at 10,000 g for 2 min, and the supernatant was transferred to a fresh tube. Isopropanol (320 μL) was added, and the mixture was incubated at −20°C for 30 min. The mixture was then centrifuged at 10,000 g for 5 min, and the supernatant was discarded. The pellet was then washed and centrifuged twice with 500 μL of ethanol (70%), and then the dried pellet was dissolved in 50 μL of TE.
A total of 186 primer combinations from a variety of sources [6, 10–12, 15, 18, 22, 30–43] were used for the construction of the linkage map. These primers were selected to obtain genome-wide coverage of 10-20 cM intervals with the least possible number of markers. The parameters considered were the length of the SSR repeat, the polymorphism level between our mapping parents and, when available, mapping information from previous publications. Complete information on these primer combinations can be found in Additional file 1: Table S1.
The first approximately 50 primer combinations used in this study were directly labelled with fluorescent dyes (6-FAM, NED or HEX). Subsequent PCR reactions were performed with indirect fluorescent labelling  using a universal 17 bp 5’ end tail sequence (AACAGGTATGACCATGA) on the forward primer, which matched a universal fluorescently labelled primer (6-FAM, HEX or ROX) . All reverse primers had a GTTT tail  on the 5′ end to minimise stutter formation. The PCR mixture was composed of 1 X Goldstar PCR buffer, 0.05 μM unlabelled forward primer with a tail, 0.2 μM unlabelled reverse primer, 0.2 μM labelled universal primer, 0.3 U of Goldstar Taq polymerase (Eurogentec Nederland B.V., Maastricht, The Netherlands) and 10 ng of DNA in a total reaction volume of 20 μL. In the case of directly labelled primers, the same mixture was used except that the forward (directly) labelled primer and reverse primer were both present at 0.2 μM concentrations. The PCR conditions were one cycle at 94°C for 3 min followed by 35 cycles at 94°C for 30 s, 50°C for 30 s and 72°C for 2 min and a final extension cycle at 72°C for 10 min for both labelling methods.
Fluorescently labelled amplicons were separated and detected using an ABI capillary automated sequencing platform (Initially ABI 3700 and later ABI 3730, Perkin Elmer Biosystems, Foster City, Calif.). The output from the ABI platform was analysed with either Genotyper 3.6 (ABI3700) or Genemapper 4.0 (ABI3730) software. Peaks corresponding to alleles were identified, and their bin ranges were defined. Next, for each sample, the software automatically identified the presence of alleles (peaks) and the area under the peak. Allele detection was checked manually and adjusted where necessary. The allelic data (size and area) for each individual (parents and progeny) were transferred to an Excel spread sheet. The analysis of the data followed the MADCE procedure for establishing allelic configurations in allopolyploid populations , which allowed us to estimate allele dose and to identify pairs of homologous alleles for each of the sub-genomes.
Construction of linkage maps
The construction of linkage maps followed the same procedure as described by van Dijk et al. . Briefly, during data analysis, the alleles were first assigned to homologous groups on the assumption that alleles shared between parents are most likely to originate from the same sub-genome, unless the data indicated otherwise. This approach allowed the definition of so-called bridge markers that link the two pairs of parental homologs. These markers are of type <hkxhk> and <efxeg> (Annotation of JoinMap® 3.0 and following versions for Cross Pollinating systems). Early data consistency checks were performed using allelic pairs as described previously by Sargent et al. . Next, linkage maps were created for each parent separately using JoinMap® 4.0 (Kyazma B.V.)  applying the regression approach and Kosambi mapping function. These separate parental maps were compared to each other to match the parental maps belonging to the same homoeologue based on the already-identified <hkxhk>, <efxeg> markers similar to the method of Barrett et al. . This information was used for increasing the number of integrated loci by converting <lmxll> and <nnxnp> markers from the same primer pair into <abxcd> markers, as well as for the validation of the previously identified <hkxhk> and <efxeg> loci. After this data check, integrated maps were created when possible. JoinMap® output was imported into Excel to check for possible genotyping errors (double recombinants) through a graphical genotyping approach . Putative double recombination events were checked up to the level of the original ABI output. The map was regarded as final when the latest corrections did not result in new putatively erroneous double recombination events, which typically required one or two rounds of corrections. The map positions of loci where both parents were homozygous were added later by imputing them from the relative positions of their homoeologous loci. Primer pairs that amplified heterologous chromosomes were never imputed and were only shown on Linkage Groups (LGs) to which their amplicons mapped. The phase information generated by JoinMap® was used to establish the parental haplotypes. Drawings of the linkage maps were first created with the software packages MapChart  and later finalised in Adobe Illustrator CS5 (Adobe Systems, San Jose, CA).
Denotation of sub-genomes
The assignment of a homoeologue letter (A, B, C or D) to a linkage group was based on the amplification efficiency of the F. vesca-derived primer pairs. The efficiency was expressed as the proportion of amplified alleles observed for all F. vesca primer pairs on a linkage group over the total numbers of alleles that were possible (amplified and null alleles).
Loci for which it was uncertain whether null alleles occurred (e.g., due to homozygosity on multiple LGs) were not included in the calculation for these LGs. Loci that amplified from heterologous chromosomes were only used for the efficiency calculations for the LGs on which they mapped.
Physical map locations of the microsatellites used in this study were obtained by blasting the SSR primer sequences to the F. vesca pseudo-chromosome assembly v 1.1 [13, 50]. When no clear hits were found, we used full-length sequences of the marker, when available. In the visualisations, the physical positions of the microsatellites in mega-base pairs were multiplied by three to better fit the scale of the genetic maps. The octoploid genetic map was represented by the homoeologue that had a good density of segregating markers and showed few inconsistencies in marker order with the other homoeologues. The diploid genetic map of Sargent et al.  was chosen as it had most primer pairs in common with our map. The genetic positions of the CO–and CX-series of markers , were imputed using data from a recent diploid genetic map . The maps were completed in Mapchart  and finalised in Adobe Illustrator CS5.
Estimation of homozygosity levels and haplotype sharing between parents
Holiday and Korona both have ancestors that occur multiple times in the known parts of their pedigree (Figure 1), due to which, part of their genomes are likely to be homozygous by descent. In addition, Holiday and Korona are likely to have shared haplotypes, as they have some ancestors in common. The theoretically expected level of homozygosity was derived from a numeric relationship matrix obtained with FlexQTL . This matrix consists of doubled kinship coefficients. The observed levels of homozygosity and haplotype-sharing were estimated using our linkage map. For this map, we identified genetic regions that had multiple (3 or more) adjacent loci where the alleles were identical, within parents (for homozygosity estimation) or between parents (for haplotype sharing/kinship estimates). Because multiple adjacent loci were used, these identical by state (IBS) regions were assumed to be identical by descent (IBD). The genetic length covered by such regions was assessed and totalled. To calculate the homozygosity levels, we divided the genetic size of the homozygous stretches by the genetic size of the genome. Briefly, the linkage map-derived kinship coefficients were calculated using Gillois identity states  for each genomic region (in cM) and their associated Jacquard condensed coefficients of identity . On our linkage map, several identity states can be distinguished. First, there are areas where no haplotype is found in common; these areas have a kinship coefficient of 0 (Gillois identity state Δ9). Next, areas with one haplotype in common between the parents have a kinship coefficient of 0.25 (state Δ8). Areas with two different haplotypes in common between the parents have a kinship coefficient of 0.5 (state Δ7). Areas where a haplotype is homozygous in one parent and the same haplotype is heterozygous in the other parent have a kinship coefficient of 0.5 (state Δ3 and Δ5). Areas where a haplotype is homozygous in both parents have a kinship coefficient of 1 (state Δ1). As an example, when a chromosome of 60 cM has identity states Δ8, Δ3, Δ1 and Δ9 on areas of 15, 10, 5, and 30 cM, respectively, its total kinship coefficient amounts to 0.23 ((15 cM*0.25 + 10 cM*0.5 + 5 cM*1 + 30 cM*0)/60 cM).
Duplicated microsatellite analysis
To investigate the underlying causes of multi-locus targeting of microsatellites we performed a BLAST search of these sequences against the Fragaria vesca reference genome (cutoff value 1*E−10). We checked whether the reference genome annotation showed a transposable element identified by LTRHarvest  overlapping the location to which the marker was BLASTed. We then used 4 kb of flanking sequence from the most significant hit and performed a BLAST search against the nucleotide collection from NCBI to establish whether the microsatellite was present within a gene.
Global mapping results
Denotation of sub-genomes
Amplification efficiency of Fragaria vesca -derived SSR primer pairs
nr of vesca derived primer pairs
Amplified alleles/total alleles
Total vesca efficiency in %
Genomic organisation of homoeologues: collinearity and re-arrangements
The overall collinearity between the homoeologues of a chromosome was very high. There were many small-scale divergences of only 1-2 cM (Additional file 2: Figure S1), but these divergences are likely due to mapping or scoring errors that were overlooked in the error checking, missing values or the presence of less informative <hkxhk> markers. In some cases, the differences in marker order were caused by the integration of the two parental maps. An example of this discrepancy can be observed in the order of markers EMFn185 and UDF067 (LG6, Figure 2). For LG6A, UDF067 occurred before EMFn185, whereas for the other homoeologues, it did not. The nearest locus for which both parents segregated was EMFn123. The distance from EMFn185 to EMFn123 was based solely on recombination events within Holiday, and the distance from UDF067 to EMFn123 was based solely on recombination events within Korona. A difference in the recombination frequency between the parents for that small region is the likely cause of the altered marker order.
Large rearrangement on chromosome 2
A second putative rearrangement was found on LG2C and occurs within the homoeologous region of the former inversion (Figure 2). Here, a large gap appeared between marker loci UFFxa02C07 (at 4 cM) and CFVCT031 (at 31 cM), and close linkage was observed between CFVCT031 and ARSFL031, whereas for linkage groups 2A and 2B, these three markers showed the opposite pattern (Figure 2). Unfortunately, it was not possible to verify whether this result occurred due to an inversion, a translocation or simply a large difference in the recombination rate because the markers that are normally located between CFVCT031 and ARSFL031 were not informative, being either homozygous or impossible to discern for LG2C. In any case, the size of the rearrangement is smaller than that for LG2D, due to the difference in the position of the homoeologous loci for UFFxa03B05.
Homozygosity & heterozygosity
Haplotype sharing (kinship)
Pedigree based kinship coefficients for Holiday and Korona and their common ancestors
List of duplicated microsatellites in the octoploid Holiday x Korona map
HxK nr of distinct alleles
Transposable element (LTRharvest)
Inside gene sequence
Large gene family
Similarity of best_scaffold 4 k region
Fragaria x ananassa beta-1,3-glucanase (BG2-2) gene
Within intron of populus trichocarpa cytochrome P450 (CYP721)
Not within region with apparant function
Fragaria x ananassa beta-1,3-glucanase (BG2-2) gene, complete cds
Not within region with apparant function
Not within region with apparant function
Not within region with apparant function
Not within region with apparant function
Not within region with apparant function
Not within region with apparant function
Within intron of ferredoxin-dependent glutamate synthase
Within intron of adiponectin receptor protein 1-like (LOC100777701), mRNA
Within intron of chromatin remodeling complex subunit (CHR923)
Within intron of molybdopterin cofactor sulfurase (ABA3)
Pyrus pyrifolia genes for F-box proteins and S2-RNase, complete cds, haplotype: S2
Not within region with apparant function
Finally, for nine primer pairs (CFVCT018, CFVCT023, CFVCT032, EMFn181, EMFn198, EMFn225, EMFn230, EMFv019 and UDF004), we could not find any putative explanation for their targeting of heterologous loci. Two of these (EMFn198 and EMFn230) did not yield sufficiently specific hits in the reference genome to do further analysis. Another two (EMFn181 and EMFn225) corresponded to loci of varying positions among the four homoeologues of a chromosome and finally, two pairs (CFVCT018 and EMFn198) corresponded to loci of varying position within a homoeologue. This result strongly indicates that markers EMFn181, EMFn225 CFVCT018 and EMFn198 represent mobile elements, which is consistent with the lack of adjacent markers showing similar behaviour.
Comparison to the diploid genome
In this study, we used the MADCE method  to develop an integrated genetic map of two octoploid strawberry cultivars that, for the first time for a polyploid plant species, included comprehensive haplotype information, even for homozygous regions and areas with null alleles. The benefits of having such extensive haplotype information are discussed in the following sections.
Map length and marker density
The map length of 2050 cM (corrected for homozygosity) is largely in line with the results of previous studies [6, 15–17, 22]. The marker density of one marker per every 3.6 cM does not provide improvements over some previously reported linkage maps in octoploid strawberry [16, 22]. However, the use of MADCE allowed us to maximise the number of segregating loci per primer pair, allowing for better comparisons of marker order retention between homoeologues. In addition, the very precise pinpointing of homozygous regions showed that marker saturation in these areas for obtaining a better resolution would be futile, as was also indicated previously by Sargent et al. .
Denotation of sub-genomes
No convention exists regarding the differentiation and notation of sub-genomes. Here, we distinguished the sub-genomes based on their sequence divergence from F. vesca using molecular markers. This approach is in contrast to those of previous studies where technical features such as number of loci and map length were used to distinguish the homoeologues [6, 14–22]. These parameters are affected by the level of homozygosity and may thus largely vary with the parents that were used for crossing. We believe that our approach, though rudimentary, is biologically more relevant and that in the near future the link between the octoploid sub-genomes and their diploid and tetraploid ancestors will become better specified, as has occurred for bread wheat .
Because strawberry is fully diploidized and for the ease of distinction, we denoted the sub-genomes with the letters A–D. However, if in the near future octoploid strawberries are shown to have originated from two tetraploid species of different origins, as suggested by Rousseau-Gueutin et al. , and as supported by our data on LGs 1, 5 and 7, then we may adopt an identical homoeologue naming convention.
Currently, the ‘Holiday’ x ‘Korona’ map is used for the mapping of thousands of SNP markers from the recently released Axiom® Strawberry Genotyping Array (also called International Strawberry 90 K SNP array or IStraw90), which was generated in a joint effort of the USDA-SCRI RosBREED project (http://www.rosbreed.org) and Affymetrix Ltd (Santa Clara, CA, USA). Now that large scale SNP arrays have come into use for strawberry, the use of a similar naming convention for octoploid maps should become straightforward due to the expected high number of common markers.
Genomic organization of homoeologues: collinearity & re-arrangements
Homoeologue collinearity and re-arrangements
The marker order between homoeologous linkage groups was highly collinear, as observed in previous studies [6, 15–17, 22]. Inconsistencies spanning very small distances are more likely to be attributable to scoring errors, missing values, the use of markers that are not equally informative, and differences in recombination rates between the single parent maps than actual genomic rearrangements.
We described for the first time a large inversion in a linkage map of octoploid strawberry. This inversion on LG2D spanned almost 30 cM and was verified in an independent mapping population. Because we could not trace this inversion in the linkage maps of Fragaria vesca and Fragaria bucharica[10, 11, 42, 51, 57], it should derive from one of the other ancestors of Fragaria x ananassa, or, less likely, may have occurred after polyploidisation. Apart from the LG2D inversion, we also found evidence that LG2C may contain a rearrangement, although the evidence was less clear. It would be interesting to further investigate these rearrangements using different progenies, higher marker densities and fluorescent in-situ hybridisation (FISH, similar to Tang et al. ), as the rearrangements could reveal interesting insights into the relationship between the octoploid strawberry and its diploid relatives.
Breeding signatures in mapping parents
The information on homozygosity and haplotype sharing generated by our genetic map revealed interesting features for both of the parental varieties, which allows us to hypothesise the possibility of breeding signatures. On a genome-wide level, the homozygosity in Holiday was found to be similar to what was theoretically expected (approximately 30%), whereas for Korona, the inbreeding level was much higher than expected (13% vs. 3.6%). Holiday has a history of heavy inbreeding, and Korona does not. It is likely that normalisation occurs for homozygosity, where extreme levels of homozygosity and extreme levels of heterozygosity are not favoured during selection. This trend may be especially true for high levels of homozygosity, as such levels are known to lead to inbreeding depression in strawberry [59, 60]. In crop species, certain traits of high agronomic value, such as plant size, adaptability and vigour, frequently favour heterozygous states . Conversely, traits such as fruit firmness, shape and size are often inherited recessively or occur due to the additive nature of alleles and thereby favour homozygous states . The differences in phenotype between the parents corroborate this hypothesis. The heavily inbred cultivar Holiday was purposefully bred for high fruit firmness and skin toughness, whereas Korona is a popular garden variety because of its adaptability and taste. However, it suffers from soft fruits and irregular fruit shape. Both varieties exhibit large fruit size and yield. There is a high probability that the differences and similarities in the distribution of homozygosity along the genome reflect the phenotype differences and similarities of the parental lines. Genes controlling fruit firmness, shape and skin vulnerability may therefore be located in areas where Holiday is homozygous and Korona is heterozygous, whereas genes controlling traits that favour agronomic performance are more likely to be located in regions where both parents are heterozygous.
IBD of haplotypes (kinship)
For Holiday and Korona, linkage map-derived kinship is more than twice as high as expected from their pedigrees. This result could indicate positive selection of these shared genetic regions. The two varieties with the theoretically largest contribution to the kinship between Holiday and Korona are Howard 17 and Aberdeen. These two varieties have been used extensively as parents in early 20th century strawberry breeding  and are therefore present in the pedigrees of many modern varieties. It is very likely that certain genomic regions of these founders are under positive selection, which would result in a higher than expected level of kinship in their descendants. Another explanation for the relatively high level of kinship could be the presence of close common ancestry among the founders of Holiday and Korona. The distribution of the shared haplotypes appears to be non-random. Certain chromosomes, such as 7A and 7D, were found to have almost three times the average haplotype-sharing, whereas all homoeologues of chromosome 4 had virtually no shared haplotypes. This result could be a coincidence but may also be due to positive or purifying selection for specific regions. Tracing the shared haplotypes between Holiday and Korona over a pedigree to their founders as well as their descendants could further clarify which haplotypes are under strong positive or purifying selection. These could be interesting for marker assisted breeding, even without knowing the associated trait(s).
Nineteen (10%) of the 186 microsatellite markers tested mapped to two or more heterologous chromosomes. Six of these had previously been reported as duplicated [6, 15, 16], and some were suggested to be remnants of a putative ancient chromosomal duplication event . Our findings did not support the presence of ancient duplication events, as we could not find a clear pattern where the same heterologous chromosome segments were consistently being amplified by multiple duplicated microsatellites. The occurrence of several of these markers in known transposable elements and large gene families provides further evidence against the hypothesis of duplication events.
Comparison to the diploid Fragariagenome
The comparative mapping revealed a generally high level of collinearity between the octoploid genetic map and the diploid physical and genetic map. This result is in line with that of previous studies [6, 15–17, 22]. However, we did find some discrepancies, which in most cases showed that the octoploid genetic map had better collinearity with the diploid Fv × Fb genetic map  than with the diploid Fragaria vesca pseudo chromosomes (v1.1). This result indicates that most of the divergence between our genetic maps and the physical map are due to limitations of the physical map. A possible explanation for this result could be that some of the scaffolds have been mapped and oriented with a BIN set comprised of a relatively low number of individuals. Knowledge of the identity of erroneously placed or erroneously oriented scaffolds as provided by this study may be of great help when fine-mapping genes of interest. However, it is not always possible to pinpoint which of the discrepancies are due to errors in the genetic maps or pseudo-chromosomes or due to real rearrangements. We could therefore not positively confirm nor reject the rearrangements observed by Sargent et al.  on LGs 1, 3 and 4.
The MADCE approach enabled the full assessment of marker haplotypes for sets of SSR-loci across the four homoeologues. It also enabled the identification of genomic rearrangements, and the discernment of homoeologues based on their similarity to the F. vesca genome. Moreover, we were able to assess the level and distribution of homozygosity and haplotype-sharing, which could indicate breeding signatures. The availability of haplotype information is crucial to go from mapping population-derived QTLs to marker-assisted selection in breeding germplasm. Haplotype information will also prove to be a valuable tool in several other aspects of strawberry breeding, such as parent selection, the verification of pedigree information and IBD analysis. New technologies such as SNP arrays and Genotyping by Sequencing (GBS) could speed up the availability of such information in allopolyploids, if coupled with appropriate methodologies to discern the different sub-genomes. We hope that with this study we have provided a significant step towards the availability of such comprehensive genetic information in strawberry.
The Socrates-Erasmus Student Exchange Programme is acknowledged for the fellowship awarded to Hulya Yilmaz Temel. Tom Davis, University of New Hampshire, is acknowledged for sharing sequence information for the design of the primers TDFM1, −2, and −6, and Paul Arens, Wageningen UR Plant Breeding is acknowledged for the design of these specific primers. Funding for this research was provided through the Technological Top Institute Green Genetics, project number 1C004RP.
- Hancock JF, Sjulin TM, Lobos GA: Strawberries. Temperate Fruit Crop Breeding. Netherlands: Springer, 2008:393-437.View ArticleGoogle Scholar
- Whitaker VM: Applications of molecular markers in strawberry. J Berry Res. 2011, 1 (3): 115-127.Google Scholar
- Lerceteau-Kohler E, Guerin G, Denoyes-Rothan B: Identification of SCAR markers linked to Rca2 anthracnose resistance gene and their assessment in strawberry germplasm. Theor Appl Genet. 2005, 111 (5): 862-870. 10.1007/s00122-005-0008-1.View ArticlePubMedGoogle Scholar
- Haymes KM, Henken B, Davis TM, Van de Weg WE: Identification of RAPD markers linked to a Phytophthora fragariae resistance gene (Rpf1) in the cultivated strawberry. Theor Appl Genet. 1997, 94 (8): 1097-1101. 10.1007/s001220050521.View ArticleGoogle Scholar
- Gaston A, Perrotte J, Lerceteau-Kohler E, Rousseau-Gueutin M, Petit A, Hernould M, Rothan C, Denoyes B: PFRU, a single dominant locus regulates the balance between sexual and asexual plant reproduction in cultivated strawberry. J Exp Bot. 2013, 64 (7): 1837-1848. 10.1093/jxb/ert047.View ArticlePubMedGoogle Scholar
- Spigler RB, Lewers KS, Johnson AL, Ashman TL: Comparative mapping reveals autosomal origin of sex chromosome in octoploid Fragaria virginiana. J Hered. 2010, 101: S107-S117. 10.1093/jhered/esq001.View ArticlePubMedGoogle Scholar
- Zorrilla-Fontanesi Y, Rambla JL, Cabeza A, Medina JJ, Sanchez-Sevilla JF, Valpuesta V, Botella MA, Granell A, Amaya I: Genetic analysis of strawberry fruit aroma and identification of O-Methyltransferase FaOMT as the locus controlling natural variation in mesifurane content. Plant Physiol. 2012, 159 (2): 851-10.1104/pp.111.188318.PubMed CentralView ArticlePubMedGoogle Scholar
- Haymes KM, Van de Weg WE, Arens P, Maas JL, Vosman B, Den Nijs APM: Development of SCAR markers linked to a Phytophthora fragariae resistance gene and their assessment in European and North American strawberry genotypes. J Am Soc Hortic Sci. 2000, 125 (3): 330-339.Google Scholar
- Davis TM, Yu H: A linkage map of the diploid strawberry, Fragaria vesca. J Hered. 1997, 88 (3): 215-221. 10.1093/oxfordjournals.jhered.a023091.View ArticleGoogle Scholar
- Sargent DJ, Davis TM, Tobutt KR, Wilkinson MJ, Battey NH, Simpson DW: A genetic linkage map of microsatellite, gene-specific and morphological markers in diploid Fragaria. Theor Appl Genet. 2004, 109 (7): 1385-1391. 10.1007/s00122-004-1767-9.View ArticlePubMedGoogle Scholar
- Sargent DJ, Clarke J, Simpson DW, Tobutt KR, Arus P, Monfort A, Vilanova S, Denoyes-Rothan B, Rousseau M, Folta KM, Bassil NV, Battey NH: An enhanced microsatellite map of diploid Fragaria. Theor Appl Genet. 2006, 112 (7): 1349-1359. 10.1007/s00122-006-0237-y.View ArticlePubMedGoogle Scholar
- Cipriani G, Pinosa F, Bonoli M, Faedi W: A new set of microsatellite markers for Fragaria species and their application in linkage analysis. J Horticult Sci Biotechnol. 2006, 81 (4): 668-675.Google Scholar
- Shulaev V, Sargent DJ, Crowhurst RN, Mockler TC, Folkerts O, Delcher AL, Jaiswal P, Mockaitis K, Liston A, Mane SP, Burns P, Davis TM, Slovin JP, Bassil NV, Hellens RP, Evans C, Harkins T, Kodira C, Desany B, Crasta OR, Jensen RV, Allan AC, Michael TP, Setubal JC, Celton JM, Rees DJG, Williams KP, Holt SH, Rojas JJR, Chatterjee M: The genome of woodland strawberry (Fragaria vesca). Nat Genet. 2011, 43 (2): 109-U151. 10.1038/ng.740.PubMed CentralView ArticlePubMedGoogle Scholar
- Lerceteau-Kohler E, Guerin G, Laigret F, Denoyes-Rothan B: Characterization of mixed disomic and polysomic inheritance in the octoploid strawberry (Fragaria x ananassa) using AFLP mapping. Theor Appl Genet. 2003, 107 (4): 619-628. 10.1007/s00122-003-1300-6.View ArticlePubMedGoogle Scholar
- Rousseau-Gueutin M, Lerceteau-Kohler E, Barrot L, Sargent DJ, Monfort A, Simpson D, Arus P, Guerin G, Denoyes-Rothan B: Comparative genetic mapping between octoploid and diploid Fragaria species reveals a high level of colinearity between their genomes and the essentially disomic behavior of the cultivated octoploid strawberry. Genetics. 2008, 179 (4): 2045-2060. 10.1534/genetics.107.083840.PubMed CentralView ArticlePubMedGoogle Scholar
- Sargent D, Passey T, Šurbanovski N, Lopez Girona E, Kuchta P, Davik J, Harrison R, Passey A, Whitehouse A, Simpson D: A microsatellite linkage map for the cultivated strawberry Fragaria x ananassa suggests extensive regions of homozygosity in the genome that may have resulted from breeding and selection. Theor Appl Genet. 2012, 124 (7): 1229-1240. 10.1007/s00122-011-1782-6.View ArticlePubMedGoogle Scholar
- Sargent DJ, Fernandez-Fernandez F, Ruiz-Roja JJ, Sutherland BG, Passey A, Whitehouse AB, Simpson DW: A genetic linkage map of the cultivated strawberry (Fragaria x ananassa) and its comparison to the diploid Fragaria reference map. Mol Breed. 2009, 24 (3): 293-303. 10.1007/s11032-009-9292-9.View ArticleGoogle Scholar
- Spigler RB, Lewers KS, Main DS, Ashman TL: Genetic mapping of sex determination in a wild strawberry, Fragaria virginiana, reveals earliest form of sex chromosome. Heredity. 2008, 101 (6): 507-517. 10.1038/hdy.2008.100.View ArticlePubMedGoogle Scholar
- Weebadde CK, Wang D, Finn CE, Lewers KS, Luby JJ, Bushakra J, Sjulin TM, Hancock JF: Using a linkage mapping approach to identify QTL for day-neutrality in the octoploid strawberry. Plant Breed. 2008, 127 (1): 94-101.Google Scholar
- Viruel MA: An SSR and RFLP linkage map for the octoploid strawberry (Fragaria x ananassa). Plant, Animal and Microbe genomes X conference. San Diego, California, USA, 2002:660Google Scholar
- Zorrilla-Fontanesi Y, Cabeza A, Dominguez P, Medina JJ, Valpuesta V, Denoyes-Rothan B, Sanchez-Sevilla JF, Amaya I: Quantitative trait loci and underlying candidate genes controlling agronomical and fruit quality traits in octoploid strawberry (Fragaria x ananassa). Theor Appl Genet. 2011, 123 (5): 755-778. 10.1007/s00122-011-1624-6.View ArticlePubMedGoogle Scholar
- Isobe SN, Hirakawa H, Sato S, Maeda F, Ishikawa M, Mori T, Yamamoto Y, Shirasawa K, Kimura M, Fukami M, Hashizume F, Tsuji T, Sasamoto S, Kato M, Nanri K, Tsuruoka H, Minami C, Takahashi C, Wada T, Ono A, Kawashima K, Nakazaki N, Kishida Y, Kohara M, Nakayama S, Yamada M, Fujishiro T, Watanabe A, Tabata S: Construction of an integrated high density simple sequence repeat linkage map in cultivated strawberry (Fragaria x ananassa) and its applicability. DNA Res. 2013, 20 (1): 79-92. 10.1093/dnares/dss035.PubMed CentralView ArticlePubMedGoogle Scholar
- Kerby K, Kuspira J: The phylogeny of the polyploid wheats Triticum-Aestivum (Bread Wheat) and Triticum-Turgidum (Macaroni Wheat). Genome. 1987, 29 (5): 722-737. 10.1139/g87-124.View ArticleGoogle Scholar
- Reinisch AJ, Dong J, Brubaker CL, Stelly DM, Wendel JF, Paterson AH: A detailed Rflp map of cotton, Gossypium-Hirsutum X Gossypium Barbadense–chromosome organization and evolution in a disomic polyploid genome. Genetics. 1994, 138 (3): 829-847.PubMed CentralPubMedGoogle Scholar
- Njuguna W, Liston A, Cronn R, Ashman TL, Bassil N: Insights into phylogeny, sex function and age of Fragaria based on whole chloroplast genome sequencing. Mol Phylogenet Evol. 2013, 66 (1): 17-29. 10.1016/j.ympev.2012.08.026.View ArticlePubMedGoogle Scholar
- Potter D, Luby JJ, Harrison RE: Phylogenetic relationships among species of Fragaria (Rosaceae) inferred from non-coding nuclear and chloroplast DNA sequences. Syst Bot. 2000, 25 (2): 337-348. 10.2307/2666646.View ArticleGoogle Scholar
- van Dijk T, Noordijk Y, Dubos T, Bink MCAM, Meulenbroek EJ, Visser RGF, Van de Weg WE: Microsatellite Allele Dose and Configuration Establishment (MADCE): an integrated approach for genetic studies in allopolyploids. BMC Plant Biol. 2012, 12: 25-10.1186/1471-2229-12-25.PubMed CentralView ArticlePubMedGoogle Scholar
- Voorrips RE, Bink MCAM, van de Weg WE: Pedimap: software for the visualization of genetic and phenotypic data in pedigrees. J Hered. 2012, 103 (6): 903-907. 10.1093/jhered/ess060.PubMed CentralView ArticlePubMedGoogle Scholar
- Fulton TM, Chunwongse J, Tanksley SD: Microprep protocol for extraction of DNA from tomato and other herbaceous plants. Plant Mol Biol Rep. 1995, 13 (3): 207-209. 10.1007/BF02670897.View ArticleGoogle Scholar
- Ashley MV, Wilk JA, Styan SMN, Craft KJ, Jones KL, Feldheim KA, Lewers KS, Ashman TL: High variability and disomic segregation of microsatellites in the octoploid Fragaria virginiana Mill. (Rosaceae). Theor Appl Genet. 2003, 107 (7): 1201-1207. 10.1007/s00122-003-1370-5.View ArticlePubMedGoogle Scholar
- Bassil NV, Gunn M, Folta K, Lewers K: Microsatellite markers for Fragaria from ‘Strawberry Festival’ expressed sequence tags. Mol Ecol Notes. 2006, 6 (2): 473-476. 10.1111/j.1471-8286.2006.01278.x.View ArticleGoogle Scholar
- Bassil NV, Njuguna W, Slovin JP: EST-SSR markers from Fragaria vesca L. cv. yellow wonder. Mol Ecol Notes. 2006, 6 (3): 806-809. 10.1111/j.1471-8286.2006.01351.x.View ArticleGoogle Scholar
- Cipriani G, Testolin R: Isolation and characterization of microsatellite loci in Fragaria. Mol Ecol Notes. 2004, 4 (3): 366-368. 10.1111/j.1471-8286.2004.00655.x.View ArticleGoogle Scholar
- Gil-Ariza DJ, Amaya I, Botella MA, Blanco JM, Caballero JL, Lopez-Aranda JM, Valpuesta V, Sanchez-Sevilla JF: EST-derived polymorphic microsatellites from cultivated strawberry (Fragaria x ananassa) are useful for diversity studies and varietal identification among Fragaria species. Mol Ecol Notes. 2006, 6 (4): 1195-1197. 10.1111/j.1471-8286.2006.01489.x.View ArticleGoogle Scholar
- Hadonou AM, Sargent DJ, Wilson F, James CM, Simpson DW: Development of microsatellite markers in Fragaria, their use in genetic diversity analysis, and their potential for genetic linkage mapping. Genome. 2004, 47 (3): 429-438. 10.1139/g03-142.View ArticlePubMedGoogle Scholar
- Hadonou M, Sargent D, Walden R, Simpson D: Characterisation of Fragaria vesca single sequence repeats (SSR) markers. Acta Hort. 99-102. 649Google Scholar
- James CM, Wilson F, Hadonou AM, Tobutt KR: Isolation and characterization of polymorphic microsatellites in diploid strawberry (Fragaria vesca L.) for mapping, diversity studies and clone identification. Mol Ecol Notes. 2003, 3 (2): 171-173. 10.1046/j.1471-8286.2003.00365.x.View ArticleGoogle Scholar
- Keniry A, Hopkins CJ, Jewell E, Morrison B, Spangenberg GC, Edwards D, Batley J: Identification and characterization of simple sequence repeat (SSR) markers from Fragaria x ananassa expressed sequences. Mol Ecol Notes. 2006, 6 (2): 319-322. 10.1111/j.1471-8286.2005.01215.x.View ArticleGoogle Scholar
- Lewers KS, Styan SMN, Hokanson SC, Bassil NV: Strawberry GenBank-derived and genomic simple sequence repeat (SSR) markers and their utility with strawberry, blackberry, and red and black raspberry. J Am Soc Hortic Sci. 2005, 130 (1): 102-115.Google Scholar
- Monfort A, Vilanova S, Davis TM, Arus P: A new set of polymorphic simple sequence repeat (SSR) markers from a wild strawberry (Fragaria vesca) are transferable to other diploid Fragaria species and to Fragaria x ananassa. Mol Ecol Notes. 2006, 6 (1): 197-200. 10.1111/j.1471-8286.2005.01191.x.View ArticleGoogle Scholar
- Oyant LHS, Crespel L, Rajapakse S, Zhang L, Foucher F: Genetic linkage maps of rose constructed with new microsatellite markers and locating QTL controlling flowering traits. Tree Genetics Genomes. 2008, 4 (1): 11-23.View ArticleGoogle Scholar
- Sargent DJ, Rys A, Nier S, Simpson DW, Tobutt KR: The development and mapping of functional markers in Fragaria and their transferability and potential for mapping in other genera. Theor Appl Genet. 2007, 114 (2): 373-384.View ArticlePubMedGoogle Scholar
- Yan Z, Denneboom C, Hattendorf A, Dolstra O, Debener T, Stam P, Visser PB: Construction of an integrated map of rose with AFLP, SSR, PK, RGA, RFLP, SCAR and morphological markers. Theor Appl Genet. 2005, 110 (4): 766-777. 10.1007/s00122-004-1903-6.View ArticlePubMedGoogle Scholar
- Schuelke M: An economic method for the fluorescent labeling of PCR fragments. Nat Biotechnol. 2000, 18 (2): 233-234. 10.1038/72708.View ArticlePubMedGoogle Scholar
- Brownstein MJ, Carpten JD, Smith JR: Modulation of non-templated nucleotide addition by taq DNA polymerase: primer modifications that facilitate genotyping. Biotechniques. 1996, 20 (6): 1004-PubMedGoogle Scholar
- Van Ooijen JW: In Joinmap 4 ®. 2006, Wageningen, Netherlands: Software for the calculation of genetic linkage maps in experimental populations. Kyazma BVGoogle Scholar
- Barrett B, Griffiths A, Schreiber M, Ellison N, Mercer C, Bouton J, Ong B, Forster J, Sawbridge T, Spangenberg G, Bryan G, Woodfield D: A microsatellite map of white clover. Theor Appl Genet. 2004, 109 (3): 596-608.View ArticlePubMedGoogle Scholar
- Young ND, Tanksley SD: Restriction fragment length polymorphism maps and the concept of graphical genotypes. Theor Appl Genet. 1989, 77 (1): 95-101. 10.1007/BF00292322.View ArticlePubMedGoogle Scholar
- Voorrips RE: MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002, 93 (1): 77-78. 10.1093/jhered/93.1.77.View ArticlePubMedGoogle Scholar
- Sargent DJ, Kuchta P, Girona EL, Zhang H, Davis TM, Celton J-M, Marchese A, Korbin M, Folta KM, Shulaev V, Simpson DW: Simple sequence repeat marker development and mapping targeted to previously unmapped regions of the strawberry genome sequence. Plant Genome. 2011, 4 (3): 165-177. 10.3835/plantgenome2011.05.0014.View ArticleGoogle Scholar
- Sargent DJ, Davis TM, Simpson DW, Folta KM, Gardiner SE: Strawberry (Fragaria spp.) Structural Genomics Genetics and Genomics of Rosaceae. 6th edition. New York: Springer, 2009:437-456.View ArticleGoogle Scholar
- Bink MCAM, Jansen J, Madduri M, Voorrips RE, Durel CE, Kouassi AB, Laurens F, Mathis F, Gessler C, Gobbin D, Rezzonico F, Patocchi A, Kellerhals M, Boudichevskaia A, Dunemann F, Peil A, Nowicka A, Lata B, Stankiewicz-Kosyl M, Jeziorek K, Pitera E, Soska A, Tomala K, Evans KM, Fernandez-Fernandez F, Guerra W, Korbin M, Keller S, Lewandowski M, Plocharski W: Bayesian QTL analyses using pedigreed families of an outcrossing species, with application to fruit firmness in apple. Theor Appl Genet. Advance online publication. doi:10.1007/s00122-014-2281-3Google Scholar
- Gillois M: La relation d’identité genetique. PhD thesis. University of Paris: Faculté des Sciences 1964.Google Scholar
- Jacquard A: The genetic structure of populations. New York: Springer-Verlag 1974.View ArticleGoogle Scholar
- Ellinghaus D, Kurtz S, Willhoeft U: LTRharvest, an efficient and flexible software for de novo detection of LTR retrotransposons. BMC Bioinformatics. 2008, 9: 18-10.1186/1471-2105-9-18.PubMed CentralView ArticlePubMedGoogle Scholar
- Rousseau-Gueutin M, Gaston A, Ainouche A, Ainouche ML, Olbricht K, Staudt G, Richard L, Denoyes-Rothan B: Tracking the evolutionary history of polyploidy in Fragaria L. (strawberry): new insights from phylogenetic analyses of low-copy nuclear genes. Mol Phylogenet Evol. 2009, 51 (3): 515-530. 10.1016/j.ympev.2008.12.024.View ArticlePubMedGoogle Scholar
- Sargent DJ, Cipriani G, Vilanova S, Gil-Ariza D, Arus P, Simpson DW, Tobutt KR, Monfort A: The development of a bin mapping population and the selective mapping of 103 markers in the diploid Fragaria reference map. Genome. 2008, 51 (2): 120-127. 10.1139/G07-107.View ArticlePubMedGoogle Scholar
- Tang XM, Szinay D, Lang C, Ramanna MS, van der Vossen EAG, Datema E, Lankhorst RK, de Boer J, Peters SA, Bachem C, Stiekema W, Visser RGF, de Jong H, Bai YL: Cross-species bacterial artificial chromosome-fluorescence in situ hybridization painting of the tomato and potato chromosome 6 reveals undescribed chromosomal rearrangements. Genetics. 2008, 180 (3): 1319-1328. 10.1534/genetics.108.093211.PubMed CentralView ArticlePubMedGoogle Scholar
- Rho IR, Woo JG, Jeong HJ, Jeon HY, Lee CH: Characteristics of F-1 hybrids and inbred lines in octoploid strawberry (Fragaria x ananassa Duchesne). Plant Breed. 2012, 131 (4): 550-554. 10.1111/j.1439-0523.2012.01958.x.View ArticleGoogle Scholar
- Shaw DV: Comparison of ancestral and current generation inbreeding in an experimental strawberry breeding population. Theor Appl Genet. 1995, 90 (2): 237-241.View ArticlePubMedGoogle Scholar
- Hochholdinger F, Hoecker N: Towards the molecular basis of heterosis. Trends Plant Sci. 2007, 12 (9): 427-432. 10.1016/j.tplants.2007.08.005.View ArticlePubMedGoogle Scholar
- Bai YL, Lindhout P: Domestication and breeding of tomatoes: what have we gained and what can we gain in the future?. Ann Bot. 2007, 100 (5): 1085-1094. 10.1093/aob/mcm150.PubMed CentralView ArticlePubMedGoogle Scholar
- Darrow GM: The strawberry; history, breeding, and physiology. 1st edition. New York: Holt 1966.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.