- Research article
- Open Access
A consensus linkage map for molecular markers and Quantitative Trait Loci associated with economically important traits in melon (Cucumis meloL.)
BMC Plant Biology volume 11, Article number: 111 (2011)
A number of molecular marker linkage maps have been developed for melon (Cucumis melo L.) over the last two decades. However, these maps were constructed using different marker sets, thus, making comparative analysis among maps difficult. In order to solve this problem, a consensus genetic map in melon was constructed using primarily highly transferable anchor markers that have broad potential use for mapping, synteny, and comparative quantitative trait loci (QTL) analysis, increasing breeding effectiveness and efficiency via marker-assisted selection (MAS).
Under the framework of the International Cucurbit Genomics Initiative (ICuGI, http://www.icugi.org), an integrated genetic map has been constructed by merging data from eight independent mapping experiments using a genetically diverse array of parental lines. The consensus map spans 1150 cM across the 12 melon linkage groups and is composed of 1592 markers (640 SSRs, 330 SNPs, 252 AFLPs, 239 RFLPs, 89 RAPDs, 15 IMAs, 16 indels and 11 morphological traits) with a mean marker density of 0.72 cM/marker. One hundred and ninety-six of these markers (157 SSRs, 32 SNPs, 6 indels and 1 RAPD) were newly developed, mapped or provided by industry representatives as released markers, including 27 SNPs and 5 indels from genes involved in the organic acid metabolism and transport, and 58 EST-SSRs. Additionally, 85 of 822 SSR markers contributed by Syngenta Seeds were included in the integrated map. In addition, 370 QTL controlling 62 traits from 18 previously reported mapping experiments using genetically diverse parental genotypes were also integrated into the consensus map. Some QTL associated with economically important traits detected in separate studies mapped to similar genomic positions. For example, independently identified QTL controlling fruit shape were mapped on similar genomic positions, suggesting that such QTL are possibly responsible for the phenotypic variability observed for this trait in a broad array of melon germplasm.
Even though relatively unsaturated genetic maps in a diverse set of melon market types have been published, the integrated saturated map presented herein should be considered the initial reference map for melon. Most of the mapped markers contained in the reference map are polymorphic in diverse collection of germplasm, and thus are potentially transferrable to a broad array of genetic experimentation (e.g., integration of physical and genetic maps, colinearity analysis, map-based gene cloning, epistasis dissection, and marker-assisted selection).
Saturated genetic linkage maps (< 1 cM between markers) are required for the efficient and effective deployment of markers in plant breeding and genomic analysis. Linkage map applications include, but are not limited to: gene mapping, positional cloning, QTL analysis, MAS, epistasis dissection, linkage disequilibrium analysis, comparative genomics, physical and genetic map integration, and genome assembly. The construction of highly saturated maps is often a time-consuming process, especially if investigators are employing different parental stocks and markers are not easily transferable. Merged maps are attractive since their integration allows for an increase in marker density without the need of additional genotyping, increased marker portability (i.e., polymorphic markers can be used in more than one population), improved marker alignment precision (i.e., congruent anchor maker position), and broader inferential capabilities (i.e., cross-population prognostication). A number of integrated linkage maps have been developed in numerous economically important crop plants including grapevine (Vitis vinifera L.) , lettuce (Lactuca sativa L.) , maize (Zea mays L.) , red clover (Trifolium pratense L.) , ryegrass (Lolium ssp.) , wheat (Triticum aestivum L.) , among others.
The genome of melon (Cucumis melo L.; 2n = 2x = 24) is relatively small (450 Mb, ), consisting of 12 chromosomes. The first molecular marker-based melon map was constructed in 1996  using mainly restriction fragment length polymorphism (RFLP) markers and morphological traits, although the markers did not cover the predicted 12 melon chromosomes. This was comparatively late for a major crop species like melon that is among the most important horticultural crops in terms of worldwide production (25 millions of tons in 2009) and which production has been increased around 40% in the last ten years . Subsequently, the first linkage maps that positioned markers on 12 linkage groups (LG) were constructed few years later, using the F2 progeny of a cross between the Korean accession PI161375 and the melon type "Pinyonet Piel de Sapo"  and two Recombinant Inbred Line (RIL) populations derived from the crosses "Védrantais" × PI161375 and "Védrantais" × PI414723 . However, these maps had few markers in common and different LG nomenclature, making comparative mapping intractable. More recently, dense linkage maps have been constructed using Simple Sequence Repeat (SSR) [12–16] and Single Nucleotide Polymorphism (SNP) [17, 18] markers. Nevertheless, although these maps share common markers, they possess large numbers of map-specific markers that makes map-wide comparisons complicated.
Melon germplasm displays an impressive variability for fruit traits and response to diseases [19–22]. Recently, part of this variability has been genetically dissected by QTL analysis [18, 23–27]. Inter-population QTL comparisons among these maps are, however, difficult given the aforementioned technical barriers.
Databases integrating genomic, genetic, and phenotypic information have been well developed in some plant species such as the Genome Database for Rosaceae , SOL Genomics Network for Solanaceae  or Gramene , and provide powerful tools for genomic analysis. In 2005, the International Cucurbit Genomics Initiative (ICuGI)  was created to further genomic research in Cucurbitaceae species by integrating genomic information in a database (http://www.icugi.org). Thirteen private seed companies funded this project, which sought to construct an integrated genetic melon map through merging existing maps using common SSR markers as anchor points. We present herein an integrated melon map, including the position of QTL controlling economically important traits, to facilitate comparative mapping comparison and to create a dynamic genetic backbone for the placement of additional markers and QTL.
Results and discussion
Construction of the integrated map
Anchor molecular markers
Based on their previously observed even map distribution, polymorphism, and repeatability, 116 SSR markers and 1 SNP marker (Additional File 1) were chosen as anchor points to integrate the eight genetic maps (Table 1). Anchor marker segregation varied among maps, where the greatest number of polymorphic anchor markers were in IRTA (Institut de Recerca i Tecnologia Agroalimentáries, Barcelona, Spain)  and INRA (Institut National de la Recherche Agronomique, Montfavet Cedex, France)  maps containing 100 and 82 anchor polymorphic markers, respectively. The minimum number of anchor polymorphic markers was recorded in the NERCV (National Engineering Research Center for Vegetables, Beijing, China)  map (35 polymorphic markers). Most of the anchor markers were originally mapped in the IRTA population, that shared a common parent (the Korean line PI 161375) with the INRA population, while the other parent was an Occidental cultivar ("Piel de Sapo" and "Vedrantais" for IRTA and INRA populations, respectively), so it was actually expected that the proportion of markers that can be transferred successfully from IRTA to INRA populations is larger than to the any other studied population developed from different germplasm.
Molecular marker segregation analysis among individual maps
Considerable and significant skewed marker segregations (p < 0.005) were detected in seven genomic regions of the DHL-based IRTA map (Table 1). Although significant skewed segregations were also detected in a region on LG VIII of the F2-based IRTA map , on LGs I, IV, and VI in NIVTS (National Institute of Vegetable and Tea Science, Mie, Japan) map  and on LGs V, VII, VIII and X in the ARO (Agricultural Research Organization, Ramat Yishay 30095, Israel) map . No significant segregation distortion was detected in the other maps used herein (data not shown). The relatively high number of genomic regions with skewed segregation detected in the DHL-based map reinforces the hypothesis that such distortion likely originated from unintentional selection during the in vitro line development process . The low number of genomic regions showing skewed segregation in most melon maps contrasts with that reported in other crops such as lettuce , red clover , sorghum , and tomato . The degree of such distortion has been correlated to the extent of taxonomic divergence between mapping parents . The use of inter-specific hybrids in order to construct genetic maps is a common strategy to ensure the availability of a high number of polymorphic markers, and in such cases segregation distortion may be frequent . However, depending on the relative frequency and intensity, segregation distortion may not interfere on the map construction. Nevertheless, such distortion may hinder the transfer of economically important alleles during plant improvement. The comparatively low frequency of segregation distortion present in melon maps may be partially explained by the use of intra-specific crosses during population development. Given the infrequent occurrence of segregation distortion in melon, the introgression of novel, economically important alleles from exotic melon germplasm into elite modern cultivars should be relatively unimpeded.
Marker polymorphism and recombination rates among individual maps
The number of polymorphic markers for individual maps ranged from 168 (USDA-ARS, Vegetable Crops Research Unit, Department of Horticulture, Madison USA) to 713 (ARO) (Table 1). INRA and IRTA maps consisted of 12 LGs, coinciding with the basic chromosome number of melon, whereas the remaining maps consisted of more LGs (see http://www.icugi.org for further details). The number of common markers in pairwise individual map comparisons was quite variable, with a mean of 40 common markers among maps. Each individual map shared between 41 and 111 markers with at least one of the other maps (Table 1). Marker order and recombination rates among markers were very consistent among maps, where significant recombination rate heterogeneities (p < 0.001) were detected between only a few marker pairs (CMN22_85-CMTCN66 in LGIII, CMAGN75-CMGA15 in LG VII, and TJ2-TJ3 in LG VIII). Similar results have been found during genetic map integration in grapevine , but more frequent recombination rate differences have been reported among integrated maps in apple (Malus domestica Borkh) , Brassica ssp. , and lettuce . Differences in locus order and recombination rates may be attributed, in part, to bands that were scored as single alleles instead of duplicated loci or to evolutionary events (chromosomal rearrangements). Nevertheless, it must be concluded from the data presented that major chromosomal rearrangements have not occurred during the recent evolutionary history (i.e., domestication) of this species.
Consensus linkage map
The construction of the integrated map described herein involved two stages: 1) the building of a framework map by merging all the available maps (Table 1) using Joinmap 3.0 ; 2) the addition of subsequent markers using a "bin-mapping" approach .
Given the high co-linearity among melon maps, 1565 markers from all maps were initially employed for map integration. However, 258 (16%) of these markers could not be included in the final integrated map. This proportion was smaller than that obtained during map integration of lettuce (19.6% ), and larger than in the grapevine integrated map (8%, ). The markers segregating within each individual map were quite complementary, what made the inclusion of a large number of markers into the final merged map possible. For example, the IRTA_F2 map was constructed with an important proportion of RFLP markers that were not used in most of the other maps. However, this map had enough RFLP markers in common with the IRTA_LDH map, which has a good proportion of common markers with INRA (68) and NIVTS (70) maps, making possible to integrate the IRTA_F2 RFLP markers in the final map.
Given the congruency detected among melon maps, the inability to incorporate some previously mapped markers into the integrated map is likely due to the lack of sufficient linkage among markers in some genomic regions, especially in small LGs drawn from some individual maps where there was a paucity of common framework map markers.
The framework integrated map contained 1307 markers (110 SNPs, 588 SSRs, 252 AFLPs, 236 RFLPs, 89 RAPDs, 6 indels, 15 IMAs, and 11 morphological traits) spanning 1150 cM that were distributed across 12 LGs with a mean genetic distance between adjacent loci of 0.88 cM (Figures 1 and 2, Additional Files 2 and 3). Integrated map length was similar to previously published maps (Table 1). While the largest marker gap was 11 cM (on the distal ends of LG × and LG IV), the remaining gaps were less than 10 cM, and occurred mainly on the distal ends of LGs (Figures 1 and 2). These gaps are likely due to the lack of sufficient common anchor markers in some maps or slight inconsistencies (distance and/or order) among maps.
Bin-mapping subsequently resulted in the addition of 285 markers (225 SNPs, 52 SSRs, 3 RFLPs, and 5 indels) producing the final integrated map containing 1592 markers (640 SSRs, 335 SNPs, 252 AFLPs, 239 RFLPs, 89 RAPDs, 15 IMAs, 11 indels, and 11 morphological traits) with a mean marker density of 0.72 cM/marker (Table 2 Figures 1 and 2, Additional Files 2 and 3, http://www.icugi.org). One hundred and seventy-eight of these markers were developed, released, or mapped for the first time for the ICuGI Consortium. The marker saturation of this integrated map is far greater than previously published maps (Table 1), increasing dramatically the number of easily transferable markers from 200  to 3353 SNPs and from 386  to 640 SSRs. Noteworthy is the fact that 17 previously bin-mapped markers were positioned on the integrated map after being genotyped in several populations. In each case, these markers mapped to their predicted positions inferred by the bin mapping approach (Table 3), demonstrating the suitability of the bin mapping set  to quickly map new markers onto the melon reference map.
Marker distribution in the integrated map varied depending on the marker type. For instance, AFLP markers clustered mainly in certain regions of LGs I, II, III, V, VI, VIII, and × (Figures 1 and 2). AFLP clustering has been commonly reported (e.g., in saturated maps of lettuce , potato  or tomato ), and it is usually associated with heterochromatic regions near centromeres. Even though regions showing AFLP clustering are likely indicative of centromeric positions, comprehensive cytogenetic analyses would be necessary to demonstrate this association in melon. In contrast, SSR, SNP and RFLP markers were generally more evenly distributed throughout the genome. Similar conclusions can not be reached about the remaining markers (RAPDs, IMAs, indels and morphological traits) due to their low number. Nevertheless, SSR marker clustering was observed in LGs III, IV, VII, VIII, XI, and XII, involving mainly SSR markers originated from genomic libraries (e.g., CMBR-SSRs ), not from ESTs. This result might indicate that those SSRs are located in repetitive DNA regions as centromeres or telomeres. However, such SSR marker clusters did not overlap those of AFLPs, even though these clusters were in the same LG (i.e., LGs III and VIII), suggesting that SSR marker clustering may be due to reasons not associated with centromeric or telomeric regions.
Integration of QTL information
Eighteen previously reported melon-mapping experiments identified 370 QTL for 62 traits (Table 4 and Additional File 4), and these were aligned in the integrated map described herein. The distribution of these QTL varied from 18 on LG IV to 57 on LG VIII (Figures 3 and 4, Additional File 5). The number of QTLs defined per trait ranged from 1 (e.g., CMV, ETH, and FB) to 40 (FS), with QTL for FS, FW, and SSC being identified in 7, 5, and 5 of the previously reported 18 mapping experiments, respectively. The number of QTL experiments in melon must be considered modest when compared with other major species, with a significant number of the traits being genetically characterized in only one or two different mapping experiments, which thereby limits the meta-analysis of QTL in this species.
Even though additional studies would be necessary to draw definitive conclusions, the position of FS QTL tend to be more consistent among experiments than those for FW and SSC QTL, mapping on LG I in six out of seven works, and on LGs II, VI, VII, VIII, XI, and XII in at least three experiments. Clustering of FW and SSC QTL was, however, only observed in LGs VIII and XI, and in LGs II, III, and V, respectively. FS is a highly heritable trait in melon, whereas FW and SSC usually show a lower heritability . The differences in QTL detection among experiments might be partially explained by trait heritability differences. Another possible explanation is that the variability of FS among the germplasm used in the experimental crosses might be controlled by a low number of common QTL with large effects, whereas a higher number of QTL with lower effects and/or more allelic variability among them might be underling SSC and FW.
Utility of the integrated molecular and QTL map
The integrated map described herein dramatically enhances the development and utility of genomic tools (i.e., markers, map-based cloning and sequencing) over previous melon maps. A large proportion of the markers in the integrated map are SSRs and SNPs, which are easily transferable across laboratories. Moreover, the populations used to construct the integrated map include genotypes from the most important market class cultivars ("Charentais", "Cantaloup", "Hami melon", "Piel de Sapo" and "U. S. Western Shipper") in broad horticultural groups (cantalupensis, inodorus, and reticulatus), guaranteeing the future utility of the markers in a broad range of cultivars and experimental crosses. The high marker density of the map allows for the selection of specific markers to customize mapping and molecular breeding applications, such as fine mapping, the development of novel genetic stocks (e.g., nearly isogenic lines and inbred backcross lines), MAS, and hybrid seed production.
The positioning of economically important QTL in the integrated map and the standardization of trait nomenclature will facilitate comparative QTL analyses among populations of different origins to provide deeper insights into the genetic control of the diverse phenotypic variability observable in melon germplasm. For example, QTL for SSC on LG III co-localize with QTL associated with SUC, GLU, and SWEET, suggesting perhaps the existence of pleiotropic effects (Figures 3 and 4). The search of candidate genes is also facilitated, as presently little correlation has been detected between candidate gene and trait for ethylene production [45, 46], fruit flesh firmness , carotenoid content [13, 18], or sugar accumulation . These associations were studied in single population, which limits the possibility of identifying associations between candidate genes and QTL. Multi-population analysis is a more powerful approach for detecting QTL/candidate gene associations. For instance, two clusters of QTL involved in carotenoid accumulation and flesh color co-localized with carotenoid-related genes: CMCRTR and BOH_1 in LG VI and CMBCYC and LYCB in LG VIII (Figures 3 and 4), and as such become candidate genes for those QTL. Similar associations can been found between genes involved in polysaccharide metabolism and transport and clusters of QTL related to fruit sugar content on LGs II, III, V, VIII, and X. Likewise, associations have been detected between ethylene biosynthesis genes and groups of QTL with effects on fruit ripening on LG VIII.
Preliminary synteny analyses have been conducted between cucumber and melon based on the IRTA SNP and EST-SSR based melon map  and the cucumber genome sequence . A large number of EST-based markers (RFLPs, EST-SSRs, and SNPs) mapped in the integrated map will facilitate synteny studies with cucumber and other cucurbit species such as watermelon, squash, and pumpkins as genomic information on such species becomes available. Most cucurbit species display a myriad of variability for economically important vegetative (e. g., branch number, sex expression) and fruit (e.g. morphology, carotenes, sugars) traits. Comparative QTL mapping based on syntenic relationships will allow the evaluation of associations between the allelic constitution at the same genetic loci and the phenotypic variability among the different cucurbit species, as is the case with fruit size between pepper and tomato in Solanaceae family .
Eight molecular marker melon maps were integrated into a single map containing 1592 markers, with a mean marker density of 0.72 cM/marker, increasing dramatically the density over previously published maps in melon. The integrated map contains a large proportion of easily transferable markers (i.e. SSRs and SNPs) and putative candidate genes that control fruit ripening, flesh softening, and sugar and carotenoid accumulation. Moreover, QTL information for 62 traits from 18 different mapping experiments was integrated into the melon map that, together with the mapped candidate genes, may provide a suitable framework for QTL/candidate gene analysis. In summary, the integrated map will be a valuable resource that will prompt the Cucurbitaceae research community for next generation genomic and genetic studies. All the individual maps, the integrated map, marker and QTL information are available at ICuGI web site (http://www.icugi.org). Researchers interested in including their QTL data into the integrated map may contact the corresponding author.
Eight mapping populations derived from seven independent crosses were used to develop the integrated map (Table 1). Three crosses involved genotypes from the two C. melo subspecies (ssp. melo and ssp. agrestis), three of them between two C. melo ssp. melo cultivars and one cross between a C. melo ssp. melo cultivar and a breeding line derived from a cross between C. melo ssp. melo and C. melo ssp. agrestis cultivars. The C. melo ssp. melo genotypes represent the most important economically market classes (Charentais, Cantaloup, Hami melon, Piel de Sapo, and U. S. Western Shipper) belonging to horticultural groups inodorus, cantalupensis, and reticulatus (Table 1) according to the classification described by Pitrat et al. (2000) . Most of the mapping populations were RILs, where two were F2 and one was a double haploid line (DHL) population (Table 1).
Development of new genomic SSR markersNew genomic SSR marker (designated DE- and DM-) were developed by Syngenta seeds. DNA plasmid libraries were constructed using approximately 1 kb fragments of sheared total DNA. SSRs were targeted via 5'-biotinylated total LNA capture probes (12-16 bases long and containing 2, 3, or 4 base repeating units) (Proligo LLC--now IDT). These probes disrupted the double helix of the library DNA at the probe sequence and as a consequence the single strand subsequently formed a double helix with the LNA probe sequence. Streptavidin coated magnetic beads (Invitrogen M-280 Dynabeads) were then used to separate the targeted plasmids from the library. Beads were washed several times and the DNA was then eluted from the beads and transformed into electrocompetent Escherichia coli DH12S cells (Life Technologies, California, USA) which were grown up and plated on large Qubit plates. Resultant colonies were then picked using the Qubit, incubated in LB broth, purified and recovered DNA was Sanger sequenced. Proprietary programs selected sequences with SSRs and designed flanking primers.
A large proportion of molecular markers developed and/or mapped in previous works (Table 1) were positioned in the integrated map. Additionally, 196 unpublished markers described bellow were included in the merged map. Additional file 2 details the major properties of these markers. On one hand, Syngenta Seeds kindly released 822 SSR markers (see above) to the ICuGI mapping project that were polymorphic in either ARO and/or INRA mapping populations. Eighty-five of them (selected based on their position calculated in in-house built genetic maps by Syngenta, unpublished results) were mapped in the ARO population and subsequently included in the merged map.
On the other hand, new 9 SSRs, 5 indels, s and 27 SNPs were released by ARO group. These indels and SNPs were detected and genotyped according to Harel-Beja et al. (2010)  in genes associated with organic acid metabolism or transport (designated OAMG-organic acid melon genes) that were cloned by two methods: (1) from melon cDNA and gDNA by PCR using degenerate primers based on conserved protein sequences; (2) ICuGI database mining. All of them were incorporated to the ARO's map .
In contrast, MU- markers are EST-SSRs were developed from their respective EST contigs available at ICuGI web page and mapped by the NERCV group. Four SNPs (AF- and AB- markers) were released by NIVTS (National Institute of Vegetable and Tea Science, Mie, Japan) group and mapped in their respective map. Finally, the unpublished indel MC264 and the SSR marker TJ22 were included in the IRTA map [10, 15, 17].
Construction of the integrated map
Various combinations of RFLP, RAPD, IMA, AFLP, SSR, indel and SNP markers had previously been employed to genotyped individuals in each of the eight mapping populations (Table 1). In order to ensure a minimum number of common anchor points among markers, 116 SSR and 1 SNP markers evenly distributed through the melon genome according to two previous linkage maps [15, 16] (Additional File 1) were selected to be genotyped in the eight mapping populations. When possible, two markers per anchor-point position were chosen to maximize the probability of identifying polymorphisms in populations examined. Standard, published protocols were employed for SSR marker genotyping [13–16, 18].
Marker segregation distortion was investigated employing Joinmap 3.0 software  in each of the mapping populations used for map merging. Given the large number of maps and markers evaluated, marker distortion was considered significant at p < 0.005 and when adjacent linked markers also showed distortion at p < 0.01. The heterogeneity of recombination frequency (REC) between common markers among different maps was also evaluated with Joinmap 3.0 and declared significant at p < 0.001.
Initially, a map was constructed for each mapping population, where LGs were defined with the "group" command with a minimum LOD score of 4.0. Groups were then assigned to LGs by comparing their marker composition with the LGs defined in previous reference maps [11, 12, 15, 17]. Groups belonging to the same LG in different populations were then integrated with the "combine groups for map integration" module of Joinmap 3.0 using the following parameters: Kosambi's mapping function LOD > 2, REC < 0.4, goodness of fit jump threshold for removal of loci = 5, performing ripple after adding 1 locus and the third integration round = No. The resulting map was designated the "framework map" and was used in further marker integrations. To add markers mapped by bin mapping [15, 17], markers defining the bins in the IRTA map were identified on the framework map. The bins were redefined in the framework map and markers were located subsequently to their respective bins from the IRTA to the framework map..
Trait and QTL definition
Traits and QTL were selected from 17 published works and 1 unpublished work (Additional Files 4 and 5) by the collaborating project researchers. Crosschecking and evaluation of recording methods allowed for the unification of trait descriptions and common abbreviations were assigned accordingly (Additional File 4). QTL were defined following the directions of the Gramene database .
Nevertheless, QTL controlling the same trait expression were often defined in independent publications and/or in different mapping populations and, consequently, QTL characterized in those different populations may correspond to the same genetic locus. Therefore, each QTL was treated independently, making it possible to notice the number of times that a QTL is reported in a similar genomic location across independent experiments.
A specific identifier was assigned to each QTL, where the first letters designate the trait abbreviation, followed by a "Q" that stands for QTL, then a letter indicating a reference to a mapping experiment (publication) followed by a digit representing the LG to which the QTL maps, and then followed by a dot and a final digit that distinguishes different QTL from the same experiment on the same LG (Additional File 5). For example, the designation FDQJ2.2 stands for one of the QTL for FD (fruit diameter) reported in the experiment J and mapping in the LG II.
QTL were defined within a marker interval according to the information presented in the original publication from which it was taken or as a personal communication from a project collaborator. If a flanking marker defining a QTL was not included in the framework map during the merging process, then the next closely linked marker was chosen for representation in the integrated map. Where only a single marker was associated with a QTL, marker position was used as both the start and stop position of the QTL. For illustration purposes, graphic representation of a QTL's position was defined in the centre of a marker interval (Figures 3 and 4).
To provide visual images of their genomic positions, integrated markers and QTL were plotted using Mapchart 2.0 . Colour codes were used to identify marker types, traits, QTL, and candidate genes in order to facilitate visualization of the co-localization of possible QTL and candidate genes involved in similar processes across different mapping experiments.
Vezzulli S, Troggio M, Coppola G, Jermakow A, Cartwright D, Zharkikh A, Stefanini M, Grando MS, Viola R, Adam-Blondon AF, Thomas M, This P, Velasco R: A reference integrated map for cultivated grapevine (Vitis vinifera L.) from three crosses, based on 283 SSR and 501 SNP-based markers. Theor Appl Genet. 2008, 117: 499-511. 10.1007/s00122-008-0794-3.
Truco MJ, Antonise R, Lavelle D, Ochoa O, Kozik A, Witsenboer H, Fort SB, Jeuken MJW, Kesseli RV, Lindhout P, Michelmore RW, Peleman J: A high-density, integrated genetic linkage map of lettuce (Lactuca spp.). Theor Appl Genet. 2007, 115: 735-746. 10.1007/s00122-007-0599-9.
Falque M, Decousset L, Dervins D, Jacob AM, Joets J, Martinant JP, Raffoux X, Ribiere N, Ridel C, Samson D, Charcosset A, Murigneux A: Linkage mapping of 1454 new maize candidate gene loci. Genetics. 2005, 170: 1957-1966. 10.1534/genetics.104.040204.
Isobe S, Kölliker R, Hisano H, Sasamoto S, Wada T, Klimenko I, Okumura K, Tabata S: Construction of a consensus linkage map for red clover (Trifolium pratense L.). BMC Plant Biol. 2009, 9: 57-10.1186/1471-2229-9-57.
Studer B, Kölliker R, Muylle H, Asp T, Frei U, Roldán-Ruiz I, Barre P, Tomaszewski C, Meally H, Barth S, Skøt L, Armstead IP, Dolstra O, Lübberstedt T: EST-derived SSR markers used as anchor loci for the construction of a consensus linkage map in ryegrass (Lolium spp.). BMC Plant Biol. 2010, 10: 177-10.1186/1471-2229-10-177.
Somers DJ, Isaac P, Edwards K: A high-density microsatellite consensus map for bread wheat (Triticum aestivum L.). Theor Appl Genet. 2004, 109: 1105-1114. 10.1007/s00122-004-1740-7.
Arumanagathan K, Earle ED: Nuclear DNA content of some important plant species. Plant Mol Biol Rep. 1991, 9: 208-218. 10.1007/BF02672069.
Baudracco-Arnas S, Pitrat M: A genetic map of melon (Cucumis melo L.) with RFLP, RAPD, isozyme, disease resistance and morphological markers. Theor Appl Genet. 1996, 93: 57-64. 10.1007/BF00225727.
FAOSTAT: Food and Agricultural Organization of the United Nations. 2011, [http://faostat.fao.org/]
Oliver M, Garcia-Mas J, Cardus M, Pueyo N, Lopez-Sese A, Arroyo M, Gomez-Paniagua H, Arus P, de Vicente MC: Construction of a reference linkage map for melon. Genome. 2001, 44: 836-845. 10.1139/g01-073.
Perin C, Hagen LS, De Conto V, Katzir N, Danin-Poleg Y, Portnoy V, Baudracco-Arnas S, Chadoeuf J, Dogimont C, Pitrat M: A reference map of Cucumis melo based on two recombinant inbred line populations. Theor Appl Genet. 2002, 104: 1017-1034. 10.1007/s00122-002-0864-x.
Gonzalo MJ, Oliver M, Garcia-Mas J, Monfort A, Dolcet-Sanjuan R, Katzir N, Arus P, Monforte AJ: Simple-sequence repeat markers used in merging linkage maps of melon (Cucumis melo L.). Theor Appl Genet. 2005, 110: 802-811. 10.1007/s00122-004-1814-6.
Cuevas HE, Staub JE, Simon PW, Zalapa JE, McCreight JD: Mapping of genetic loci that regulate quantity of beta-carotene in fruit of US Western Shipping melon (Cucumis melo L.). Theor Appl Genet. 2008, 117: 1345-1359. 10.1007/s00122-008-0868-2.
Cuevas HE, Staub JE, Simon PW, Zalapa JE: A consensus linkage map identifies genomics regions controlling fruit maturity and beta-carotene-associated flesh color in melon (Cucumis melo L.). Theor Appl Genet. 2009, 119: 741-756. 10.1007/s00122-009-1085-3.
Fernandez-Silva I, Eduardo I, Blanca J, Esteras C, Pico B, Nuez F, Arus P, Garcia-Mas J, Monforte AJ: Bin mapping of genomic and EST-derived SSRs in melon (Cucumis melo L.). Theor Appl Genet. 2008, 118: 139-150. 10.1007/s00122-008-0883-3.
Fukino N, Yoshioka Y, Kubo N, Hirai M, Sugiyama M, Sakata Y, Matsumoto S: Identification of QTLs for resistance to powdery mildew and SSR markers diagnostic for powdery mildew resistance genes in melon (Cucumis melo L.). Theor Appl Genet. 2008, 118: 165-175. 10.1007/s00122-008-0885-1.
Deleu W, Esteras C, Roig C, Gonzalez-To M, Fernandez-Silva I, Gonzalez-Ibeas D, Blanca J, Aranda MA, Arus P, Nuez F, Monforte AJ, Picó MB, Garcia-Mas J: A set of EST-SNPs for map saturation and cultivar identification in melon. BMC Plant Biol. 2009, 9: 90-10.1186/1471-2229-9-90.
Harel-Beja R, Tzuri G, Portnoy V, Lotan-Pompan M, Lev S, Cohen S, Dai N, Yeselson L, Meir A, Libhaber SE, Avisar E, Melame T, van Koert P, Verbakel H, Hofstede R, Volpin H, Oliver M, Fougedoire A, Stalh C, Fauve J, Copes B, Fei Z, Giovannoni J, Ori N, Lewinsohn E, Sherman A, Burger J, Tadmor Y, SchaVer AA, Katzir N: A genetic map of melon highly enriched with fruit quality QTL and EST markers, including sugar and carotenoid metabolism genes. Theor Appl Genet. 2010, 121: 511-533. 10.1007/s00122-010-1327-4.
Stepansky A, Kovalski I, Perl-Treves R: Intraspecific classification of melons (Cucumis melo L.) in view of their phenotypic and molecular variation. Plant Syst Evol. 1999, 217: 313-332. 10.1007/BF00984373.
Monforte AJ, Eduardo I, Abad S, Arus P: Inheritance mode of fruit traits in melon: Heterosis for fruit shape and its correlation with genetic distance. Euphytica. 2005, 144: 31-38. 10.1007/s10681-005-0201-y.
Burger Y, Sa'ar U, Paris HS, Lewinsohn E, Katzir N, Tadmor Y, Schaffer AA: Genetic variability for valuable fruit quality traits in Cucumis melo. Isr J Plant Sci. 2006, 54: 233-242. 10.1560/IJPS_54_3_233.
Burger Y, Jonas-Levi A, Gurski E, Horev C, Saar U, Cohen R: Variation in antifungal activity in extracts from Momordica plants. Isr J Plant Sci. 2010, 58: 1-7. 10.1560/IJPS.58.1.1.
Perin C, Hagen LS, Giovinazzo N, Besombes D, Dogimont C, Pitrat M: Genetic control of fruit shape acts prior to anthesis in melon (Cucumis melo L.). Mol Genet Genomics. 2002, 266: 933-941. 10.1007/s00438-001-0612-y.
Monforte AJ, Oliver M, Gonzalo MJ, Alvarez JM, Dolcet-Sanjuan R, Arus P: Identification of quantitative trait loci involved in fruit quality traits in melon (Cucumis melo L.). Theor Appl Genet. 2004, 108: 750-758. 10.1007/s00122-003-1483-x.
Eduardo I, Arus P, Monforte AJ, Obando J, Fernandez-Trujillo JP, Martinez JA, Alarcon AL, Alvarez JM, van der Knaap E: Estimating the genetic architecture of fruit quality traits in melon using a genomic library of near isogenic lines. J Am Soc Hortic Sci. 2007, 132: 80-89.
Obando-Ulloa JM, Eduardo I, Monforte AJ, Fernandez-Trujillo JP: Identification of QTL related to sugar and organic acid composition in melon using near-isogenic lines. Sci Hortic. 2009, 121: 425-433. 10.1016/j.scienta.2009.02.023.
Boissot N, Thomas S, Sauvion N, Marchal C, Pavis C, Dogimont C: Mapping and validation of QTLs for resistance to aphids and whiteflies in melon. Theor Appl Genet. 2010, 121: 9-20. 10.1007/s00122-010-1287-8.
Genome Data Base for Rosaceae. [http://www.rosaceae.org/]
SOL Genomics Network. [http://solgenomics.net/]
Cucurbit Genomics Database. [http://www.icugi.org/]
Fang L, Yong X, Yue Z, Di C, Jian-ming F, Shao-gui G, Guo-yi G, Hong-ping Y, Ming-zhu W, Hai-ying Z: Construction of permanent genetic map and comparative analysis of Xinjiang Hami melon (Cucumis melo L. ssp. melo. convar. ameri (Pang.) Greb). Acta Hortic Sinica. 2009, 36: 1767-1774.
Gonzalo MJ, Claveria E, Monforte AJ, Dolcet-Sanjuan R: Parthenogenic haploids in melon: generation and molecular characterization of a doubled haploid line population. J Amer Soc Hortic Sci. 2011, 136: 145-154.
Mace ES, Rami JF, Bouchet S, Klein PE, Klein RR, Kilian A, Wenzl P, Xia L, Halloran K, Jordan DR: A consensus genetic map of sorghum that integrates multiple component maps and high-throughput Diversity Array Technology (DArT) markers. BMC Plant Biol. 2009, 9: 13-10.1186/1471-2229-9-13.
Villalta I, Reina-Sanchez A, Cuartero J, Carbonell EA, Asins MJ: Comparative microsatellite linkage analysis and genetic structure of two populations of F-6 lines derived from Lycopersicon pimpinellifolium and L. cheesmanii. Theor Appl Genet. 2005, 110: 881-894. 10.1007/s00122-004-1906-3.
Taylor DR, Ingvarsson PK: Common features of segregation distortion in plants and animals. Genetica. 2003, 117: 27-35. 10.1023/A:1022308414864.
Zamir D, Tadmor Y: Unequal segregation of nuclear genes in plants. Bot Gaz. 1986, 147: 355-358. 10.1086/337602.
N'Diaye A, de Weg WEV, Kodde LP, Koller B, Dunemann F, Thiermann M, Tartarini S, Gennari F, Durel CE: Construction of an integrated consensus map of the apple genome based on four mapping populations. Tree Genet Genomes. 2008, 4: 727-743. 10.1007/s11295-008-0146-0.
Nicolas SD, Le Mignon G, Eber F, Coriton O, Monod H, Clouet V, Huteau V, Lostanlen A, Delourme R, Chalhoub B, Ryder CD, Chevre AM, Jenczewski E: Homeologous recombination plays a major role in chromosome rearrangements that occur during meiosis of Brassica napus haploids. Genetics. 2007, 175: 487-503. 10.1534/genetics.106.062968.
Van Ooijen JW, Vorrips RE: JoinMap® Version 3.0, Software for thecalculation of genetic linkage map. Plant Research InternationalWageningen, The Netherlands; 2001.
Howad W, Yamamoto T, Dirlewanger E, Testolin R, Cosson P, Cipriani G, Monforte AJ, Georgi L, Abbott AG, Arus P: Mapping with a few plants: Using selective mapping for microsatellite saturation of the Prunus reference map. Genetics. 2005, 171: 1305-1309. 10.1534/genetics.105.043661.
van Os H, Andrzejewski S, Bakker E, Barrena I, Bryan GJ, Caromel B, Ghareeb B, Isidore E, de Jong W, van Koert P, Lefebvre V, Milbourne D, Ritter E, van der Voort JNAMR, Rousselle-Bourgeois F, van Vliet J, Waugh R, Visser RGF, Bakker J, van Eck HJ: Construction of a 10,000-Marker Ultradense Genetic Recombination Map of Potato: Providing a framework for accelerated gene isolation and a genomewide physical map. Genetics. 2006, 173: 1075-1087. 10.1534/genetics.106.055871.
Haanstra JPW, Wye C, Verbakel H, Meijer-Dekens F, van den Berg P, Odinot P, van Heusden AW, Tanksley S, Lindhout P, Peleman J: An integrated high density RFLP-AFLP map of tomato based on two Lycopersicon esculentum × L. pennellii F2 populations. Theor Appl Genet. 1999, 99: 254-271. 10.1007/s001220051231.
Ritschel PS, Lins TCL, Tristan RL, Buso GSC, Buso JA, Ferreira ME: Development of microsatellite markers from an enriched genomic library for genetic analysis of melon (Cucumis melo L.). BMC Plant Biol. 2004, 4: 1-14. 10.1186/1471-2229-4-1.
Perin C, Gomez-Jimenez M, Hagen L, Dogimont C, Pech JC, Latche A, Pitrat M, Lelievre JM: Molecular and genetic characterization of a non-climacteric phenotype in melon reveals two loci conferring altered ethylene response in fruit. Plant Physiol. 2002, 129: 300-309. 10.1104/pp.010613.
Moreno E, Obando JM, Dos-Santos N, Fernandez-Trujillo JP, Monforte AJ, Garcia-Mas J: Candidate genes and QTL for fruit ripening and softening in melon. Theor Appl Genet. 2008, 116: 589-602. 10.1007/s00122-007-0694-y.
Huang SW, Li RQ, Zhang ZH, Li L, Gu XF, Fan W, Lucas WJ, Wang XW, Xie BY, Ni PX, et al: The genome of the cucumber, Cucumis sativus L. Nature Genet. 2009, 41: 1275-U1229. 10.1038/ng.475.
Paran I, van der Knaap E: Genetic and molecular regulation of fruit and plant domestication traits in tomato and pepper. J Exp Bot. 2007, 58: 3841-3852. 10.1093/jxb/erm257.
Pitrat M, Hanel P, Hammer K: Some comments on infraspecific classification on cultivars of melon. Acta Hort. 2000, 510: 29-36.
Ni J, Pujar A, Youens-Clark K, Yap I, Jaiswal P, Tecle I, Tung CW, Ren L, Spooner W, Wei X, Avraham S, Ware D, Stein L, McCouch S: Gramene QTL database: development, content and applications. Database: J Biol Databases Curation. 2009, bap005.
Voorrips RE: MapChart: Software for the graphical presentation of linkage maps and QTL. J Hered. 2002, 93: 77-78. 10.1093/jhered/93.1.77.
This work was supported in part by SNC Laboratoire ASL, Ruiter Seeds B.V., Enza Zaden B.V., Gautier Semences S.A., Nunhems B.V., Rijk Zwaan B.V., Sakata Seed Inc, Semillas Fitó S.A., Seminis Vegetable Seeds Inc, Syngenta Seeds B.V., Takii and Company Ltd, Vilmorin & Cie S.A., and Zeraim Gedera Ltd (all of them as part of the support to the ICuGI); the grants AGL2009-12698-C02-02 from the Spanish "Ministerio de Ciencia e Innovación" to AJM. NK lab was supported in part by Research Grant Award No. IS-4223-09C from BARD, the United States - Israel Binational Agricultural Research and Development Fund, and in part by Israel Science Foundation Grant No. 386-06, De Ruiter Seeds, Enza Zaden, Keygene, Rijk Zwaan, Sakata Seed Corporation, Semillas Fitó, Syngenta Seeds and Vilmorin Clause & Cie. AD was supported by a JAE-Doc contract from "Consejo Superior de Investigaciones Científicas" (CSIC-Spain). MF was supported by a postdoctoral contract from CRAG. The research carried out at YX's laboratory was supported by Chinese funds (Grant No.2008-Z42(3), 5100001, 2010AA101907).
AJM coordinated the map integration study, provided the marker and QTL data of the IRTA mapping populations, performed the map merging, and drafted the manuscript. MF obtained additional genotype data for the IRTA mapping population. GF integrated QTL information into the merged map, AD assisted in the map merging, prepared tables, and graphic representations and helped to draft the manuscript. PZ and JB formatted the data for representation with C-maps for publication in the ICuGI web site. ZF is the responsible for the ICuGI web site. JES, JZ, and HC provided new marker and QTL data of the USDA-ARS mapping populations; JES assisted with manuscript editing. NF and SM provided new marker and QTL data of the NITVS mapping population and new SNP markers MO provided new marker mapping data for the ARO mapping population and GD developed the DE and DM SSR markers. CD, NB and MP provided new marker and QTL data of the INRA mapping population. RH and PK assisted with map merging construction. RHB, GL, VP, SC, AS, NK, provided new SSR and OGM markers, marker and QTL data of the ARO mapping population. YX and HYZ provided new SSR markers from melon ESTs, and also marker and QTL data of the NERCV mapping population. NF and SM provided the SSR markers used as anchor points for map integration and marker and QTL data of the NITVS mapping population. JGM was the coordinator of the ICuGI project and participated in the design of the study. All authors have read and approved the final manuscript