Mapping the Flavor Contributing Traits on "Fengwei Melon" (Cucumis melo L.) Chromosomes Using Parent Resequencing and Super Bulked-Segregant Analysis

We used a next-generation high-throughput sequencing platform to resequence the Xinguowei and Shouxing melon cultivars, the parents of Fengwei melon. We found 84% of the reads (under a coverage rate of “13×”) placed on the reference genome DHL92. There were 2,550,000 single-nucleotide polymorphisms and 140,000 structural variations in the two genomes. We also identified 1,290 polymorphic genes between Xinguowei and Shouxing. We combined specific length amplified fragment sequencing (SLAF-seq) and bulked-segregant analysis (super-BSA) to analyze the two parents and the F2 extreme phenotypes. This combined method yielded 12,438,270 reads, 46,087 SLAF tags, and 4,480 polymorphic markers (average depth of 161.81×). There were six sweet trait-related regions containing 13 differential SLAF markers, and 23 sour trait-related regions containing 48 differential SLAF markers. We further fine-mapped the sweet trait to the genomic regions on chromosomes 6, 10, 11, and 12. Correspondingly, we mapped the sour trait-related genomic regions to chromosomes 2, 3, 4, 5, 9, and 12. Finally, we positioned nine of the 61 differential markers in the sweet and sour trait candidate regions on the parental genome. These markers corresponded to one sweet and eight sour trait-related genes. Our study provides a basis for marker-assisted breeding of desirable sweet and sour traits in Fengwei melons.


Introduction
The melon (Cucumis melo L.) is an important horticultural crop, and has diverse phenotypes and genotypes that result in variations in fruit quality including sugar and acid content, aroma, and fruit shape [1]. The principal component of fruit quality in sweet melons is the content of sugar, primarily sucrose [2][3][4]. Details of the metabolic pathways involved in sucrose accumulation in melon fruit have been previously described [5]. The organic acid content in most melon fruit is very low (pH >6.0) [6]. However, some varieties accumulate high levels of acid (pH <5) [7][8][9][10][11], and some breeders of melons have combined high acidity and high sugar traits to develop dessert cultivars [12][13][14]. The hereditary characteristics of citric acid content and titratable acid have been studied in melon fruit [7]. A single major QTL for pH has been found to co-localize with major QTLs for citric and malic acids [8]. Furthermore, the linkage map position for pH genes has been reported and an SSR marker was identified as being tightly associated with pH [10]. A previous study showed in bred melons that the high sugar and low pH traits were inherited independently of each other, and that the acidity trait was particularly associated with citric acid levels [12]. However, the relationship between genomic variations and the sour trait phenotype has not been reported in melon fruit.
Fengwei melon is a hybrid obtained from the Shouxing and Xinguowei varieties, and has a unique sweet and sour taste that results from the combined traits of its parents. Shouxing [Cucumis melo L. ssp. melo cover. ameri (Pang) Greb] is a crossbred from multiple parents, while Xinguowei [Cucumis melo L. ssp. melo convar. ameri (Pang) Greb, a sour tasting Hami melon line bred using high-dose Co 60 γ-ray irradiation] has a unique sweet and sour taste [13].
Over the past few years, various genetic and genomic resources for melon have become available. New genetic maps have been reported for numerous traits such as fruit size and shape, climacteric ripening, netting, color, various metabolites, total soluble solids, and pH [8,9,15,16]. Resequencing technology and bioinformatics tools have helped to unravel the genetic variations and diversity on a genome-wide scale [17]. Garcia-Mas et al. [18] sequenced the complete genome of the double haploid line, DHL92, and provided valuable information about this important cucurbit. Blanca et al. [19] resequenced the transcriptomes of 67 melon genotypes using a high-throughput SOLiD™ system. Bulked-segregant analysis (BSA) is used for fine-mapping of phenotypic mutants [20] and is substantially enhanced by whole genome sequencing [21]. The use of sequence-based markers such as restriction site-associated DNA (RAD) markers, helps to achieve high-density discovery and genotyping in large populations [22]. However, these technologies usually fail in the presence of repetitive genomic sequences that usually are not useful in mapping experiments. Sun et al. [23] developed an enhanced reduced representation library (RRL) sequencing method, known as specific length amplified fragment sequencing (SLAF-seq) technology, that combines locus-specific amplification and high-throughput sequencing. This technique avoids the selection of repetitive genomic sequences and is highly accurate, low-cost, and fast. Therefore, the present study aimed to localize important genes related to the sweet and sour traits in Fengwei melon (the F1 hybrid of the Shouxing and Xinguowei varieties). We used parental resequencing and comparative genomics to compare whole genome variations and functional genes between the two parental lines, in order to provide references for the gene mapping of sweet and sour traits. We used a combination of SLAF-seq and BSA technology (super-BSA) to help identifying the genomic regions affecting these traits in pools of extreme phenotypes from the F 2 population (Fengwei melon ). Our specific objectives were: 1) to obtain sweetness and sourness candidate genes through bioinformatics mapping of candidate regions closely correlated with sweet and sour traits; and 2) to perform a combined analysis of genes in these regions and SLAF differential markers.

Plants
The Shouxing parental line [Cucumis melo L. ssp. melo convar. ameri (Pang.) Greb] is a crossbred from multiple parents and its fruit has a sweet taste. The other parental line used in this study, Xinguowei [designated guowei (sour taste Hami melon line bred from the 76-2 line using high-dose Co 60 γ-ray irradiation [13]) × PI140471], is a sour-tasting variety. Fengwei melon is the F 1 hybrid obtained using Shouxing as the female parent and Xinguowei as the male parent. Its fruit has a unique sour and sweet taste. We obtained the F 2 segregating population by self-pollinating the F 1 hybrid (Fengwei melon).
Ten plants of each of the female parent, male parent, and F 1 and F 2 populations were grown at the experimental field of Xinjiang Academy of Agricultural Science (Xinjiang Turpan, China) in an open field during spring 2012 (F 2 individuals, n = 479) and autumn 2012 (F 2 individuals, n = 500), and in a greenhouse in Xinjiang Turpan during autumn 2012 (F 2 individuals, n = 514). No endangered or protected species were involved in this study.
Total soluble solids and pH of the fruit A digital refractometer (Atago Co, Ltd, Tokyo, Japan) was used to measure the total soluble solids (TSS). A digital pH meter (Ohaus, Shanghai, China) was used to measure the pH of the juice [5]. The juice used in these experiments was extracted from a single fruit harvested from each plant. About 5 g of ripe flesh was taken from the equatorial portion of the fruit after removing the rind and seed cavity. The flesh was squeezed by hand, and the juice was collected and centrifuged to remove solid matter.

DNA extraction
DNA was extracted from young leaf tissue obtained from 10 female parent plants, 10 male parent plants, and 50 plants from each of the three extreme phenotypes from the F 2 population, as described by Murray and Thompson [24].

Resequencing and detection of genome variations in the two parents
Paired-end (PE) and mate-pair Solexa libraries were prepared, according to the manufacturer's instructions (Illumina, Hayward, CA, USA). Shouxing and Xinguowei DNAs were randomly sheared. After electrophoresis, DNA fragments of the desired length were gel-purified. Adaptor ligation and DNA cluster preparation were performed, and the resulting DNA libraries were subjected to Solexa sequencing using an Illumina GAIIx (Illumina) platform. Low-quality reads (<20), reads with adaptor sequences, and duplicate reads were filtered out, and the remaining high-quality data were used for mapping, which was carried out with the Burrows-Wheeler alignment (BWA) software [25].
The localization of SNPs, indels, and SVs was based on annotations of gene models provided by reference genome databases (https://melonomics.net) [18], and was carried out using the snpEff software (http://snpeff.sourceforge.net/SnpEff_manual.html). The three types of polymorphism were annotated as genic (in gene regions) or intergenic (in non-gene regions).
Genic SNPs, indels, and SVs were classified as exonic or intronic according to their localization. The GeneWise software [30] was used to separate the exonic SNPs into synonymous or nonsynonymous types.

Annotation of differential gene function between Shouxing and Xinguowei
We compared the non-synonymous SNPs of the two parental lines and the SVs (including indels) in exons against the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, using the BLAST program with a cutoff E-value of 1×10 −5 .
Super-BSA pools, library construction, and high-throughput sequencing Five different DNA pools were prepared by mixing equal amounts of DNA from the female Shouxing bulk 'M' (sweet) and male Xinguowei bulk 'P' (sour). Fifty plants were used from each of the three extreme phenotypes from the F 2 population bulk (sweet, sour, and non-sweet non-sour). The library was constructed as described by Sun et al. [23], with small modifications. The GC content, repeated sequences, and genetic characteristics of the DNA pools were analyzed using SLAFPredict (Biomarker, Beijing, China). The genomic DNA pools were digested using the XhoI and MseI restriction enzymes, followed by PCR amplification, fragment amplification, fragment selection, fragment extraction and amplification, and fragment sequencing using the Illumina GAIIx system.

SLAF marker development and polymorphism analysis
All SLAF pair-end reads with clear index information were clustered based on sequence similarity, as detected by BLAT v.34 [31] (-tileSize = 11, -stepSize = 11, -minScore = 30). Sequences with over 90% identity were grouped in one SLAF locus, as described by Sun et al. [23]. Melon is a diploid species, and one locus contains at most four SLAF tags; therefore, groups containing more than four tags were filtered out as repetitive SLAFs. In this study, SLAFs with a sequence depth <164 were defined as low-depth SLAFs and were filtered out. SLAFs with 2, 3, or 4 tags were identified as polymorphic SLAFs and considered to be potential markers.

Fine mapping of the sweet and sour traits
All markers were identified based on the parental origin of alleles M and P, according to the sequencing depth. M sweet represented the depth for the sweet phenotype from the female line; P non-sweet non-sour represented the depth for the non-sweet non-sour phenotype from the male line; M non-sweet non-sour represented the depth for the non-sweet non-sour phenotype from the female parent; and P sour represented the depth for the sour phenotype from the male line. The following ratios were calculated: Ratio_sweet = M sweet /P non-sweet non-sour ; and Ratio_sour = P sour /M non-sweet non-sour . In the case that P non-sweet non-sour = 0, Ratio_sweet was set to 1000; in the case that M non-sweet non-sour = 0, Ratio_sour was set to 1000. The thresholds for association were set at a ratio 3. Although the ratio measurement employed in this analysis has not been widely used in previous studies, it was recently described to be a good approach to BSA analysis [32]. It can be inferred that a ratio 3 means an SNP-index ratio 0.5. The SNP-index indicates the proportion of reads harboring a SNP that is different from the reference sequence [33].
We performed fine mapping of the genomic region according to the sequences of the sweet and sour trait-related SLAF markers in the scaffold. We used BLAT v.34 [31] (-tileSize = 11, -stepSize = 11, -minScore = 30) for every SLAF marker (https://melonomics.net) to determine the SLAF position in the genetic map.

Combined analysis of regions associated with sweet and sour traits and parental resequencing data
The genes in regions associated with sweet and sour traits and the differential SLAF markers were matched with the parental resequencing data using BLAT v.34 [31] (-tileSize = 11, -step-Size = 11, -minScore = 30). We used the annotation in 'gff' files of the genome to determine whether the markers covered the gene. We determined the variation in traits by analyzing the SNPs between genes and genome where correlated markers were located; or by analyzing the SNPs between samples (whether the SNP variation induced genetic variation and influenced the variation in the traits). Finally, we extracted the corresponding genes according to the physical location of the markers (in genes) and that correlated with the sweet and sour traits through gene function annotation.

Variability in fruit TSS and pH in the two parents and the F 1 and F 2 populations
To estimate the variation in the TSS and pH of the ripe fruits, TSS and pH were measured in the parents and the F 1 (Fengwei melon) and three F 2 populations grown in Turpan during two growing seasons (open field in spring and autumn 2012, and greenhouse in autumn 2012). Results showed that there was significant variability in the TSS and pH traits between the parental, F 1 , and F 2 populations (p0.01; Welch's ANOVA test; Table 1). Fruits grown in the spring showed higher levels of TSS and pH than those grown in autumn. In addition, fruits grown in autumn had higher levels of TSS and pH when grown in a greenhouse than when grown in an open field. We used the greenhouse results from autumn 2012 in subsequent parent resequencing and super-BSA. Melon fruits with TSS 8 and pH >5.3 were considered sweet, while those with TSS <8 and pH 5.3 were considered sour. We characterized Shouxing as sweet (TSS >12, pH >6), Xinguowei as sour (TSS <8, pH <5), and F 1 Fengwei as sweet and sour (TSS >12, pH <5). The F 2 population produced four types of trait: sweet (TSS 8 and pH >5.3), sour (TSS <8 and pH 5.3), sweet and sour (TSS >12, pH <5) and non-sweet non-sour (TSS <8, pH >5.3). We selected the sweet (TSS >12, pH >6), sour (TSS <8, pH <5) and non-sweet non-sour (TSS <8, pH >5.3) pools from the F 2 lines for subsequent super-BSA.

Discovery of polymorphism within genic regions related to fruit flavor traits
We used a Solexa Genome Analyzer II for genomic resequencing of the Shouxing and Xinguowei parental lines. Shouxing yielded 29,882,152 101-bp short sequences while Xinguowei produced 30,603,252 short sequences. We used the BWA software to match the 24,796,210 and 25,746,516 reads in the two parental lines to the reference genome, DHL92 ( Table 2).
The SAMtools software [26] was used to search for SNPs in the genomes of the two parental lines. There were 1,278,396 SNPs between Shouxing and DHL92, and 1,268,988 SNPs between Xinguowei and DHL92 (Fig 1). According to the differences in the nucleotide substitutions, the SNPs were either transitions (C/T or G/A, shown in red) or transversions (C/G, T/A, A/C or G/T, shown in black) (Fig 1). For both parental lines, the proportion of transitions (Ts) was higher than that of transversions (Tv). The Ts/Tv ratios for Shouxing and Xinguowei, with reference to the DHL92 genome, were 2.44 and 2.42, respectively. There were 171,072 and 176,345 genic SNPs, respectively. Of the genic SNPs in Shouxing and Xinguowei, 48,642 and 49,070 were located in exons, respectively, and included 26,670 and 26,773 non-synonymous SNPs, respectively (Fig 2A and 2B). Of the 1,152,860 differential SNPs between Shouxing and Xinguowei, 22,322 were non-synonymous. We also used the Pindel [27] and BreakDancer [28] software to detect indels and SVs in the genomes of the two parental lines (Fig 1). We identified 26,802 indels and 10,109 SVs in Shouxing melon, and 28,546 indels and 12,968 SVs in Xinguowei melon. There were 16,703 and 17,903 genic SVs (including indels), and 1,186 and 1,218 exonic SVs (including indels), in Shouxing and Xinguowei, respectively (Fig 3A and 3B). Of the 64,265 and 71,337 SVs (including indels) detected in these two parental lines, INS and DEL accounted for 94.98% and 95.54%, respectively. In Shouxing and Xinguowei, other types of SV and CTX accounted for 4.58% and 4.04%, respectively; In DELs and INS for 0.23% and 0.22%, respectively; CTX for 0.15% and 0.13%, respectively; and INV for 0.06% and 0.07%, respectively. To verify the accuracy of the detection of SNPs and SVs, we randomly selected variations in 400 SNPs and 100 SVs (INS and DEL) of 100 to 300 bp in length for validation in the two inbred lines by PCR and sequencing (S1, S2, S3 and S4 Tables). Among these, 394 SNPs and 92 SV loci could be amplified and sequenced. These data verified 98.75% of the SNPs and 92.0% of the SVs, indicating that our re-sequencing data were reliable.

Identification of pathways associated with fruit flavor traits
Previous studies have speculated that genetic variations in melon fruit may contribute to phenotypic differences in traits [8,9,15,16]. Therefore, we focused our analysis on non-synonymous SNPs, indels, and SVs in the exons of genes in the two parental lines. We performed a preliminary examination of the genes' functional annotations to assess the influence of genotypic variations on phenotypic variations. Pathway enrichment analysis for the 575 Shouxing polymorphic genes in 101 KEGG pathways and the 715 Xinguowei polymorphic genes in 110 KEGG pathways enabled the identification of several pathways contributing to fruit flavor traits. The pathways identified included starch and sucrose metabolism, fructose and mannose metabolism, galactose metabolism, citrate cycle metabolism, nicotinate and nicotinamide metabolism, folate biosynthesis, and vitamin B6 metabolism (S5 Table and S6 Table). We identified four genes in Xinguowei that were annotated in the pathway for nicotinate and nicotinamide metabolism: nudix hydrolase-19, nicotinamide mononucleotide adenylyltransferase-3, nicotinate phosphoribosyltransferase, and L-aspartate oxidase. We also annotated two genes as part of the folate pathway: folylpolyglutamate synthase and bifunctional dihydrofolate reductase-thymidylate synthase. Threonine synthase, whose encoding   gene was identified as polymorphic, participates in vitamin B6 metabolism, which is associated with fruit acid. The higher content of nicotinate, folate, and vitamin B6 in Xinguowei (data not shown) may be associated with variations in the genes of these pathways. In Xinguowei, three polymorphic genes were identified and annotated in the valine, leucine, and isoleucine pathways. These included: pyruvate dehydrogenase E1 component subunit beta-3, branched-chain-aminoacid aminotransferase-like protein-3, and 2-isopropylmalate synthase-2.
The three main soluble sugars in melons are sucrose, glucose and fructose. Sugar metabolism and accumulation determines the sweetness of the melon. We found polymorphisms between the two parents in the starch and sucrose pathway, fructose and mannose pathway, and galactose pathway. Organic acids in fruits are essential components that reflect fruit flavor characteristics. The citrate cycle pathway and its key steps are necessary for the accumulation of organic acids in fruits. Therefore, DNA polymorphisms in these genes could explain the differences in sweet and sour traits between the two inbred lines.
In both Xinguowei and Shouxing, we also identified genes encoding members of the lipoxygenase family, alcohol dehydrogenase family, and cadinene synthase, as well as genes involved in the synthesis of linalool derivatives or the nerolidol synthase of the linalool monomer.

Super-BSA, SLAF markers and polymorphism analysis
Super-BSA analysis was carried out to genotype the bulk pools of the two parental lines and three F 2 population: 50 individual plants with TSS 8 and pH >5.3, 50 with TSS <8 and pH 5.3, and 50 with TSS <8 and pH >5.3. The data from the above pools corroborated with each other to some extent, eliminating the need for biological replicates.
We obtained 12,438,270 reads (Table 3), and used the BLAT matching software to cluster these reads (http://www.blat.net) and map them to the reference DHL92 genome after correction. SLAFs were selected with depths larger than 10×, resulting in a total of 46,087 SLAF tags being obtained with an average depth of 161.81×. Fig 4 shows the even distribution of SLAF tags on the scaffolds (>5 Mb).

Association mapping of genes related to sweet and sour traits
We compared the 4,480 SLAF markers with the allelic frequencies of the parental lines and found 2,800 sweet trait-related SLAF markers and 2,711 sour trait-related SLAF markers from the parents. We identified 114 differential SLAF markers related to the sweet trait (Ratio_sweet 3; Fig  6) and 215 related to the sour trait (Ratio_sour 3; Fig 7). Examination of the distribution of the differential SLAF markers on the CM3.5_scaffold00001 suggested that they correlated most intensively with the sweet trait, showing a coverage area of 0.9 Mb containing 10 markers (marking density of 0.09 Mb/marker) (Fig 8). The markers in the CM3.5_scaffold00018 were most intensively correlated with the sour trait, presenting a coverage area of 2.75 Mb containing 14 markers (marking density of 0.2 Mb/marker) (Fig 9).
The presence of three or more differential SLAF markers (Ratio_sweet 3) in succession helped identifying the regions associated with the sweet and sour traits. There were six sweet trait-related regions containing 13 differential SLAF markers and 62 genes (Table 5), and 23 sour trait-related regions containing 48 differential SLAF markers and 185 genes ( Table 6).
To anchor the above 13 sweet trait-related and 48 sour trait-related SLAF markers in the scaffolds to the genetic map, we performed BLAT analysis of every SLAF marker in the scaffolds against the genetic map (https://melonomics.net/genetic_map/map_set_info). We mapped the 13 sweet trait-related SLAF markers to between 0.86 and 26.46 kb on chromosome 6; 3.54 and 3.87 kb on chromosome 10; 1.99 and 2.01 kb on chromosome 11; and 7.09 and   Table 8). We discarded 10 sour trait-related SLAF markers, each of which mapped to more than two different locations in the genome. We used the method described by Harel-Beja et al. [15] to measure the glucose and sucrose content (i.e. sweetness) of the melons, and the method of Cohen et al. [8] to measure the pH, citrate content and malic acid content (i.e. acidity). Based on the genetic mapping positions, we mapped sucrose QTLs, TSS QTLs and glucose QTLs to linkage groups LG2, LG3, LG4, LG5 and LG8 [14], while pH, citrate and malate QTLs were mapped to LG4, LG8, LG11 and LG12 [7,9,14,15].

Combined analysis of regions associated with sweet and sour traits and parental resequencing data
We studied the genes from the regions associated with sweet and sour traits, and found that some of the SLAF differential markers (13 sweet trait-related and 48 sour trait-related) matched the genome-wide genetic variation in the parental resequencing data. Some sweet (62) and sour (185) trait-related variations within the associated regions showed polymorphisms between the two parents, including SNPs (70, 247) and SVs (38, 106) ( Table 9). Comparative analysis of SLAF markers and parental resequencing data showed that only one marker (SLAF18745), located in one gene, N-acetylglucosaminyl transferase III (MELO3C011944T1), was correlated with the sweet trait. However, seven SLAF markers located in eight genes were associated with  (Table 10 and S7 Table). Therefore, there was rich variation in sweet and sour trait-related genes within the two parental lines compared with DHL92.

Discussion
The rapid development of sequencing technologies and bioinformatics tools has provided a starting point for unraveling the tremendous genetic diversity that occurs in the genome. Genetic variations on a genome-wide scale has been reported in several model organisms [17,29,[35][36][37]. Techniques such as parental genome resequencing and BSA have helped identifying markers linked to genes or QTLs. These markers can be further used in the construction of genetic maps and in high-throughput genotyping polymorphisms (SSRs, SVs, or indels). We resequenced the two parental lines (Shouxing and Xinguowei) of the Fengwei melon to uncover nearly two million SNPs, indels, and SVs, and identified 1290 polymorphic genes associated with sweet and sour traits. Some of these genes are involved in metabolic pathways for sucrose, organic acids, aromas, and vitamins. Our study provides the first report on genome-wide patterns of genetic variation in melon, and will be valuable to future genotype-phenotype studies and molecular breeding. Melon fruits accumulate various levels of soluble sugars, organic acids, vitamins, and aromatic volatiles [5], and these compounds affect fruit quality through complex networks of metabolic pathways that are active during fruit ripening. Transcriptome and gene cloning have identified the metabolic pathways for sugars that are related to sucrose accumulation in melon fruit [5]. Our results demonstrated that the two parents of Fengwei melon (Xinguowei and Shuoxing) showed polymorphisms of genes in several pathways involved in the biosynthesis of sucrose, fructose, mannose, and galactose. DNA polymorphisms in genes coding for components of these pathways may account for the differences in sweet taste between the two parents. Moreover, a transcriptome study that we have carried out has identified differential expressions of these pathways and genes in melon fruit (unpublished data).
Studies of organic acids in other fruit species have shown that the acid content of a fruit is determined by the balance between acid synthesis and degradation. The genes associated with the citrate cycle pathways have been documented to play a critical role in the accumulation of organic acids in fruits [38][39][40][41][42][43][44]. Our current findings reveal that the two parents showed polymorphisms of genes related to the citrate cycle pathways. DNA polymorphisms in genes involved in these pathways may contribute to the differences in sour taste between the two parents. Indeed, our transcriptome study showed differing expressions of the genes involved in the citrate cycle pathways in melon fruit (unpublished data). Recently, a pH-trait gene (CmPH) that has a major effect on fruit acidity has been cloned by map-based techniques and characterized in melon [45]. In future studies, the mechanisms of organic acid accumulation in melon fruit could be further explored by the continuous development of sour traits in melon fruit, and research into comparative genomics and transcription in fruits of other species. The lipoxygenase family of genes governs fruit aroma, which is an important trait contributing to fruit quality [46]. The release of aromatic substances signifies the maturity of the melon fruit. Aromatic substances are closely associated with the health and nutrition of humans and are indices for evaluating melon fruit quality. We identified genes annotated as encoding members of the lipoxygenase family in the two parental lines, including genes that play roles in the ripening and aging processes in melon fruit such as ethylene synthesis, after-ripening, softening, and aroma formation [47,48]. We also identified genes encoding the alcohol dehydrogenase family that participates in ethylene signal transduction and the biosynthesis of aromatic substances [49][50][51], cadinene synthase that catalyzes farnesyl pyrophosphate cyclization to  form (+)-δ-cadinene [52], and genes involved in the synthesis of linalool derivatives or the nerolidol synthase of the linalool monomer [53]. We also identified polymorphic genes in Xinguowei that were annotated in the metabolic pathways for valine, leucine, isoleucine, nicotinate, and nicotinamide. Valine, leucine, and isoleucine are precursors involved in the synthesis of most branched ester aromatic substances [54,55]. DNA polymorphisms in genes associated with the above pathways may also contribute to differences in taste between the two parents. Our transcriptome study also demonstrated differing expressions of these pathways in melon fruit (unpublished data). Sherman et al. [56] combined BSA and microarrays to map the pH trait to chromosome 8 in melon. Although Harel-Beja et al. [15] measured the sweetness in melons with respect to glucose and sucrose content, none of their experiments overlapped with ours. We mapped pH and acidity (i.e. sourness) traits to chromosome 12 (LG12), as also reported in a previous study by Cohen et al. [8]. We attribute the differences in the mapping data between our work and previous studies to the different accessions used. Compared with 'Dulce' and 'PI414723', the Hami melon inbred lines Shouxing and Guowei (sour-taste Hami melon lines bred from high-dose Co 60 γ-ray irradiation) may have different genes related to sweet and sour traits. Some of the crucial fruit quality traits such as flavor, aroma, and vitamin content, have not received enough attention. The polymorphic genes annotated in our study provide abundant resources for research concerning the molecular mechanisms associated with melon fruit traits.
Takagi et al. [57] used BSA in combination with high-throughput sequencing technologies, including QTL-seq and BSR-seq, to fine-map partial resistance to fungal rice blast disease and seedling vigor trait in rice. Trick et al. [58] used BSA and high-throughput sequencing to map Gpc-B1, a wheat quantitative trait locus associated with an increased grain content of protein, zinc, iron, and gl3, while Liu et al. [59] mapped the phenotypic mutants of maize using these techniques. We combined SLAF-seq and super-BSA to identify rapidly the genomic regions associated with sweet and sour traits, using members of the F 2 population with extreme phenotypes. The success of fine mapping the sweet and sour traits using the super-BSA technology provides strong technical support for its future application. In this method, high-density scanning of 10,000+ SNPs on large-scale bulk pools generated a massive amount of sequencing information that covered the whole genome. We compared the difference in the frequency of occurrence of different genes marked by a SNP in the two bulk pools and determined the molecular markers and fine mapping regions closely correlated with the traits. Thus, we directly used the sequencing results for molecular marker development. Moreover, our research provides a convenient way to develop trait-related functional molecular markers for molecular marker-assisted breeding.
Methods combining the efficiency of genome-wide variation mining and genomic mapping, including resequencing, super-BSA, QTL-seq [57], and BSR-seq [58,59], will  dramatically accelerate crop improvement in a cost-effective manner. The technologies that take full advantage of the rapidly declining cost of genome sequencing are expected to contribute to the on-going efforts aimed at addressing the world food security problem by reducing breeding time.
Supporting Information S1