Next Article in Journal
Mutual Modulation of the Activities of Human CYP2D6 and Four UGTs during the Metabolism of Propranolol
Next Article in Special Issue
Research Progress on Anthocyanin-Mediated Regulation of ‘Black’ Phenotypes of Plant Organs
Previous Article in Journal
Axonal Protection by Oral Nicotinamide Riboside Treatment with Upregulated AMPK Phosphorylation in a Rat Glaucomatous Degeneration Model
Previous Article in Special Issue
Exploring Genetic and Epigenetic Changes in Lingonberry Using Molecular Markers: Implications for Clonal Propagation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

QTL Mapping and Genome-Wide Association Study Reveal Genetic Loci and Candidate Genes Related to Soluble Solids Content in Melon

Jiangsu Yanjiang Institute of Agricultural Sciences, Jiangsu Academy of Agricultural Sciences, Nantong 226012, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Curr. Issues Mol. Biol. 2023, 45(9), 7110-7129; https://doi.org/10.3390/cimb45090450
Submission received: 1 August 2023 / Revised: 21 August 2023 / Accepted: 25 August 2023 / Published: 26 August 2023
(This article belongs to the Special Issue Molecular Breeding and Genetics Research in Plants)

Abstract

:
Melon (Cucumis melo L.) is an economically important Cucurbitaceae crop grown around the globe. The sweetness of melon is a significant factor in fruit quality and consumer appeal, and the soluble solids content (SSC) is a key index of melon sweetness. In this study, 146 recombinant inbred lines (RILs) derived from two oriental melon materials with different levels of sweetness containing 1427 bin markers, and 213 melon accessions containing 1,681,775 single nucleotide polymorphism (SNP) markers were used to identify genomic regions influencing SSC. Linkage mapping detected 10 quantitative trait loci (QTLs) distributed on six chromosomes, seven of which were overlapped with the reported QTLs. A total of 211 significant SNPs were identified by genome-wide association study (GWAS), 138 of which overlapped with the reported QTLs. Two new stable, co-localized regions on chromosome 3 were identified by QTL mapping and GWAS across multiple environments, which explained large phenotypic variance. Five candidate genes related to SSC were identified by QTL mapping, GWAS, and qRT-PCR, two of which were involved in hydrolysis of raffinose and sucrose located in the new stable loci. The other three candidate genes were involved in raffinose synthesis, sugar transport, and production of substrate for sugar synthesis. The genomic regions and candidate genes will be helpful for molecular breeding programs and elucidating the mechanisms of sugar accumulation.

1. Introduction

Melon (Cucumis melo L.; 2n = 2x = 24) is an economically important Cucurbitaceae crop grown around the globe. Based on differences of ovary pubescence, melon cultivars are divided into two subspecies: C. melo ssp. agrestis and C. melo ssp. melo [1,2]. The sweetness is an important quality trait of melons that affects consumer preferences, and the soluble solids content (SSC) is a key index of melon sweetness. Identifying genomic regions and genes controlling SSC will be helpful for molecular breeding.
Linkage analysis is an effective approach for detecting quantitative trait loci (QTLs) related to complex agronomic traits using biparental segregating populations, and several QTLs were found to regulate the accumulation of soluble solids content in melon [3,4,5,6,7,8,9,10,11,12,13,14]. However, a considerable number of QTLs for soluble solids content detected through diverse backgrounds have wide intervals, and these QTLs are difficult to apply to marker-assisted selection (MAS). Due to the rapid growth of sequencing technology, GWAS has proven a powerful approach in detecting genetic variations of important agronomic traits in melon [15,16,17,18,19]. Previously, some key loci related to soluble solids content were detected via GWAS in tomato [20,21], apricot [22], peach [23], and watermelon [24]. Leida et al. [25] identified a stable region on chromosome 3 associated with sugar content via association analysis using 251 single nucleotide polymorphism (SNP) markers. Wang et al. [26] detected 11 SNPs associated with soluble solids content in melon via GWAS using genotyping-by-sequencing (GBS).
Association mapping does not require a biparental segregating population and is useful for analyzing multiple alleles [27], but the population substructure and allele frequency could bring about false-positive results [28]. Combining association analysis and linkage analysis has been successfully used for dissecting complex quantitative traits in rice [29], maize [30], and soybean [31]. Wei et al. [29] identified eight candidate genes for rice seedlings in response to high temperature stress using linkage analysis and GWAS. Xu et al. [30] dissected the genetic basis of zinc deficiency tolerance in maize at the seedling stage via linkage mapping and GWAS, and three candidate genes were identified as being responsible for Zn transport. Hu et al. [31] identified a gibberellin 3b-hydroxylase gene related to soybean 100-seed weight through joint linkage mapping and GWAS which enhanced photosynthesis and increased seed yield in soybean. At present, association analysis and linkage analysis have been used to dissect the complex agronomic traits of melon. Nimmakayala et al. [32] detected several QTLs for fruit firmness through GWAS and QTL mapping. Perpiñá et al. [33] detected several stable QTLs related to SSC using linkage analysis with introgression lines and association analysis in backcross populations. Oren et al. [34] detected a major netting QTL on melon chromosome 2 by combining linkage analysis and association analysis. Du et al. [35] identified a candidate genomic region associate with the fruit surface groove of melon via linkage analysis and association analysis. These genetic loci identified by GWAS and linkage analysis will make MAS breeding more accurate and effective. Based on the rapid growth of sequencing technology, joint GWAS and linkage analysis will play an increasingly important role in gene mining for important agronomic traits in melon. Raffinose and stachyose are the core sugars of cucurbit crops transported in the phloem [36]. Based on the rapid advances of molecular biology and genomics, many genes related to sugar metabolism and accumulation during fruit growth have been reported such as alpha-galactosidase (AGA) [37,38], sugars will eventually be exported transporter (SWEET) [37,38,39], vacuolar sugar transporter (VST) [40], sucrose transporter (SUT) [39,41], tonoplast sugar transporter (TST) [42,43,44], hexose transporter [45], tonoplast H+/sugar antiporters [46], trehalose-6-phosphate (T6P) [47,48], sucrose synthase (SUS), sucrose phosphate synthase (SPS), and invertase [49].
Previous studies have shown that the integration of linkage analysis and GWAS provided a higher power and resolution to investigate complex quantitative traits. In the present study, we constructed a recombinant inbred line (RIL) population containing 146 lines derived from two oriental melon materials with different levels of sweetness. Then, a high-resolution genetic map of the RIL population was constructed to detect QTLs related to SSC in melon. In addition, a total of 213 melon accessions were utilized to detect the significance of SNPs associated with SSC via GWAS. Five candidate genes were identified by QTL mapping, GWAS, and qRT-PCR, which played important roles in synthesis and hydrolysis of raffinose, sucrose hydrolysis, and sugar transport. The genomic regions and candidate genes will be helpful for molecular breeding programs and elucidating the mechanisms of sugar accumulation.

2. Materials and Methods

2.1. Plant Materials and Phenotypic Data Analysis

The plant material contained an RIL population and a GWAS panel. The RIL population, consisting of 146 RILs, was derived from two oriental melon cultivars (1214 with low sweetness and 1228 with high sweetness). The RIL population was grown in a greenhouse at Xueyao Experiment Station (32°14′ N, 120°64′ E) of Jiangsu Yanjiang Institute of Agricultural Sciences (Nantong, China) in the spring of 2020, 2021, and 2022. The GWAS panel, consisting of 213 melon accessions that included 145 accessions of ssp. agrestis, 67 accessions of ssp. melo, and 1 Mapao melon accession, was planted in a greenhouse at Xueyao Experiment Station in the spring of 2020 and 2021 and the spring/autumn of 2021. A randomized block design with three replicates was utilized to evaluate the plant materials. Five plants were grown for each melon germplasm. One fruit per plant was harvested at 30–35 days after pollination, and the SSC of fruit flesh was measured with a hand-held digital PAL-1 refractometer (Atago). ANOVA for SSC was calculated in R using the aov function. Based on the phenotypic data of two populations for multiple environments, the best linear unbiased predictions (BLUP) were calculated by the R package lme4.

2.2. DNA Extraction and Whole-Genome Resequencing

The young leaves of RIL population and GWAS panel were utilized to extract DNA using the CTAB method. The genomic DNA samples were fragmented to a size of 350 bp for sequencing on the Illumina HiSeq PE150 platform of Novogene, China. The melon reference genome HS_ZJU [50] and its annotation were downloaded online (http://www.ncbi.nlm.nih.gov, accessed on 19 July 2021). BWA [51] was utilized to map the clean reads to the melon reference genome using the command ‘mem -t 4 -k 32 -M’. SNP calling was performed using SAMtools [52] with a Bayesian approach. The high-quality SNPs were performed using Genome Analysis Toolkit (GATK) software 4.0.4.0 [53]. ANNOVAR software (Version 20130520) [54] was utilized to annotate the SNPs identified in this study.

2.3. Linkage Analysis

The sequencing depths of parents (1214 and 1228) and RILs were 20-fold and 5-fold, respectively. The RILs were genotyped by polymorphic SNPs with an aa × bb segregation pattern between the parental materials. The recombinant breakpoints were evaluated by the sliding window approach [55]. The genotype was performed with a sliding window of 15 SNPs. The high-resolution genetic map of the RIL population was constructed with bin markers using JoinMap 4.1 (https://www.kyazma.nl/index.php/JoinMap/, accessed on 13 May 2022). The ALLMAPS program [56] was utilized to visualize the relationships between the physical and genetic locations. The composite interval mapping (CIM) method of Windows QTL Cartographer 2.5 software (https://brcwebportal.cos.ncsu.edu/qtlcart/WQTLCart.htm, accessed on 7 June 2022) was utilized to detect QTLs. The permutation test of MapQTL6.0 (https://www.kyazma.nl/index.php/MapQTL/, accessed on 7 June 2022) was utilized to calculate the LOD threshold values for the SSC.

2.4. Association Analysis

The sequencing depths of 211 melon accessions were 10-fold, and two accessions (1214 and 1228) were sequenced with a depth of 20-fold. The sequencing depth ≥ 6, missing rate ≤ 0.5, and minor allele frequency ≥ 0.05 were considered standard in identifying SNPs. Finally, a total of 1,681,775 SNPs were used for association analysis. The software TreeBeST 1.9.2 (https://treesoft.sourceforge.net/treebest.shtml, accessed on 2 April 2022) was utilized to create an individual-based neighbor-joining (NJ) tree using p-distance with 1000 bootstrap replications. The phylogenetic tree was visualized using MEGA6.0 (http://www.megasoftware.net/, accessed on 2 April 2022). The program ADMIXTURE (version 1.23) [57] was utilized to calculate the population genetics structure using an expectation maximization algorithm. To evaluate the LD decay in the association population, the software PopLDdecay 3.40 [58] was utilized to calculate the degree of linkage disequilibrium coefficient (r2) between pairwise SNPs with the command ‘-n -dprime -minMAF 0.05’. The GEMMA 0.98.1 software package [59] was utilized to perform association analysis, and the significant P-value threshold was approximately 1 × 10−6.

2.5. Candidate Genes Identification and qRT-PCR Analysis

The co-detected regions by linkage analysis and GWAS, and the ±100 kb regions around the significant SNPs located within the reported QTL regions were considered significant regions to search for candidate genes. Annotations of these genes were analyzed following the melon reference genome ‘HS_ZJU’, and genes related to sugar metabolism and accumulation were considered candidate genes. The fruits of the two parents (1228 and 1214) harvested at 20DAP, 25DAP, and 30DAP were used to extract RNA for qRT-PCR analysis. The cDNA was synthesized using All-in-One First-Strand cDNA Synthesis SuperMix for qPCR (TransGen, Beijing, China). Gene expression was performed on the ABI7500 Real-Time PCR System (Applied Biosystems, Singapore) using Green qPCR SuperMix (TransGen, Beijing, China). The data of relative gene expression was analyzed according to the 2−ΔΔCt method with the 18s rRNA used as the internal control. The primers used in this study are shown in Table S8.

3. Results

3.1. Phenotypic Variation of Soluble Solids Content in Two Populations

In this study, 146 RILs and 213 melon accessions were used. The SSC of RILs was evaluated in the Nantong greenhouse in the spring of 2020, 2021, and 2022 (Table S9). The SSC of the 213 melon accessions was evaluated in the spring of 2020 and 2022 and the spring/autumn of 2021 (Table S10). The descriptive statistics including the mean, range, standard deviation, kurtosis, skewness, and coefficient of variation (CV) for SSC of two panels are presented in Table 1.
The SSC of RILs ranged from 5.67% to 14.50% (with a mean of 10.43%) in spring 2020, from 5.83% to 14.07% (with a mean of 11.01%) in spring 2021, and from 7.15% to 14.80% (with a mean of 11.29%) in spring 2022. The SSC of 213 melon accessions ranged from 3.00% to 17.30% (with a mean of 10.01%) in spring 2020, from 3.00% to 17.43% (with a mean of 9.86%) in spring 2021, from 3.55% to 15.50% (with a mean of 8.97%) in autumn 2021, and from 3.80% to 18.98% (with a mean of 10.68%) in spring 2022. The CV of the RILs ranged from 12.93% to 14.96%, while the CV of the accessions ranged from 26.17% to 30.21%. The skewness and kurtosis in Table 1 and the histograms in Figure 1a,b show that the data of SSC have a normal distribution. The melon accessions (G), environment (E), and the genotype × environment interaction (G × E) had a significant effect (p < 0.01) on the SSC (Table 1). The broad-sense heritability (H2) for SSC were 0.65 and 0.72 in two populations (Table 1).
The fruits of the two parents (1228 and 1214), were harvested at five growth stages to investigate the changes in the SSC (Figure 1c). With fruit development, the SSC of 1228 and 1214 accumulated continuously and peaked at 30 days after pollination (DAP), but the SSC of 1214 was significantly lower than that of 1228 at each stage (Figure 1c). The different levels of SSC led to the differences in sweetness between these two cultivars.

3.2. Population Sequencing and Linkage Map Construction

In order to construct a high-resolution genetic map, 146 RILs together with two parents were re-sequenced on the Illumina HiSeq PE150 platform. In total, 9.05 Gb (19.96-fold genome coverage) and 10.34 Gb (23.04-fold genome coverage) clean bases were obtained for two parents (1214 and 1228), respectively (Table S1). A total of 364 Gb clean bases were obtained for the 146 RILs with high quality (Q20 ≥ 95.42%, Q30 ≥ 88.72%), and the average sequencing depth for each individual was 6.7-fold (Table S1). A total of 380,864 SNPs were detected between the parents, and 111,500 SNPs were identified with the aa×bb segregation pattern. Then, a total of 1427 bin markers were obtained to construct a high-resolution genetic map (Figure 2a). The total length of the bin map was 1254.34 cM, and the average interval between the adjacent markers was 0.88 cM (Table 2). The genetic and physical positions of the bin markers on the 12 chromosomes were compared with each other, and a high level of collinearity was detected between the linkage map and the reference genome of melon cultivar ‘HS_ZJU’ (Figure 2b). However, there were several chromosomal locations of the bins that displayed inconsistencies with the genetic map.

3.3. Identification of QTLs for Soluble Solids Content

Under three growth environments and BLUP, 14 QTLs related to SSC were identified using Windows QTL Cartographer 2.5 software with the bin map. These 14 QTLs were scattered on chromosomes 1, 2, 3, 4, 6 and 8, as follows: 1 QTL on chromosome 1, 4 QTLs on chromosome 2, 4 QTLs on chromosome 3, 2 QTLs on chromosome 4, 2 QTLs on chromosome 6, and 1 QTL on chromosome 8 (Figure 3, Table S3). The physical distances of these QTLs ranged from 0.39 to 3.15 Mb, and the average distances was 0.88 Mb. The LOD value of these QTLs ranged from 3.2 to 5.7, and the phenotypic variation explained (PVE) values ranged from 6.7% to 13.8%. All the QTL showed negative additive effects, indicating that alleles from the high-SSC parent 1228 contributed to the greater SSC. Based on the physical position of flanking markers, the 14 QTLs were clustered into 10 common QTLs. Of the 10 common QTLs, one QTL was repeatedly detected across three environments, and two QTLs were detected across two environments. In addition, seven common QTLs were mapped to the same loci of previously reported QTLs (Table S3) and confirmed the reliability of the genetic map.

3.4. GWAS for Soluble Solids Content

We re-sequenced 213 melon accessions at an average sequencing depth of 9.7× on the Illumina HiSeq PE150 platform, generating 941.6 Gb resequencing data of high quality (Q20 ≥ 94.43%, Q30 ≥ 86.92%, Table S2). The raw sequence data reported in this paper have been deposited in the Genome Sequence Archive at the National Genomics Data Center, China National Center for Bioinformation/Beijing Institute of Genomics, Chinese Academy of Sciences, under accession number CRA011939 and are publicly accessible at https://ngdc.cncb.ac.cn/gsa (released on 24 July 2023). The resequencing data were mapped to the melon genome ‘HS_ZJU’, and 1,681,775 high-confidence SNPs (missing rate ≤ 0.5; minor allele frequency ≥ 0.05) were detected. Of these, 99,547, 51,623, 304,931, 1,029,614, and 96,114 SNPs were located in upstream regions, exons, introns, intergenic regions, and downstream regions, respectively (Table S4). The markers were scattered on the melon genome, with the lowest (107,238) and highest (174,458) number of markers present on chromosome 9 and chromosome 1, respectively (Figure 4a). Additionally, the highest SNP frequency was 5.04 SNPs/kb on chromosome 7. The linkage disequilibrium (LD) level was calculated as the physical distance at which r2 decreased to half of the maximum value, and LD decay was estimated at 103 kb (Figure 4b). From the phylogenetic tree constructed with the high-confidence SNPs (Figure 4c), the 213 accessions were divided into two subpopulations: subpopulation I primarily consisted of ssp. agrestis and subpopulation II consisted of ssp. melo. The result was further supported by a population structure plot (Figure 4d). In addition, the ‘Mapao’ melon accession was located within ssp. agrestis, and this result was consistent with a previous study [17].
Based on the data of four growth environments and BLUP, GWAS of 213 melon accessions was performed with 1,681,775 SNPs using GEMMA software. At a significance level of p < 10−6, a total of 211 SNPs significant associated with SSC was found (Figure 5; Table S5). Among these SNPs, s6-26722606 had the lowest p-value (−log10 p = 10.44), which was detected in 2022SPR. The two chromosomes with the highest number of markers were chromosome 9 with 61 SNPs and chromosome 12 with 34 SNPs. The intersections based on the SNPs detected across four environments and BLUP are presented in Figure 6 as the UpSet plot [60]. The UpSet plot showed the number of significant markers detected in each environment and BLUP, and also showed the number of the same markers detected across different environments (Figure 6). In Figure 6, 29 significant SNPs were detected coherently across three environments, 10 significant SNPs were detected coherently across four environments, and 1 significant SNP was detected coherently across all environments. Compared to the SNP loci detected by GWAS and the QTLs detected in previous studies, 138 SNPs overlapped with the reported QTLs (Table S6).

3.5. Co-Detected Regions by QTL Mapping and GWAS

Linkage analysis and association analysis could be complementary methods for QTLs detection. The results of linkage analysis and GWAS were analyzed together to identify stable regions significantly associated with SSC. Five QTLs identified by linkage analysis and 13 SNPs identified by GWAS were co-localized on chromosome 3, 4, and 6 (Table 3). QTL qSSC3-1 and five SNPs, and qSSC3-2 and four SNPs were co-localized on chromosome 3. QTL qSSC4-1 and one SNPs, and qSSC4-2 and two SNPs were co-localized on chromosome 4. QTL qSSC6-2 and one SNP were co-localized on chromosome 6. The two co-localized regions on chromosome 3 were identified across multiple environments with a large PVE (>10%) and likely included genes affecting the accumulation of SSC.

3.6. Identification of Candidate Genes Related to SSC

The genes located within QTL intervals and the LD regions of significant SNPs, which were related to sugar metabolism and accumulation, were considered as candidate genes. As a result of qSSC3-1 and qSSC3-2 being co-detected by linkage analysis and GWAS across multiple environments, candidate genes related to SSC were searched for in these two stable regions on chromosome 3. A total of 113 genes were detected within the two co-localized intervals (28.5–28.8 Mb and 28.9–29.3 Mb) on chromosome 3, which consisted of 46 and 67 genes, respectively (Table S7). Of the 46 genes in qSSC3-1, the invertase gene CWINV3 (MELO 06759) may be a candidate gene for SSC as a key enzyme for converting sucrose into glucose and fructose. The expression of CWINV3 was significantly different between the high sweetness and low sweetness oriental melon accessions at different fruit growth stages [61]. Meanwhile, we found an alpha-galactosidase (MELO 06810) in qSSC3-2, which was related to the production of galactose and sucrose during phloem unloading [62]. Yang et al. [50] reported that MELO 06810 had a parallel shift of amino acid in Cucurbitaceae via positive selection gene analysis, which may play an important role in phloem transport.
Based on the LD decay of 103 kb (Figure 4b), the ±100 kb regions around 138 significant SNPs that overlapped with the reported QTLs (Table S6) were considered significant regions to search for candidate genes. In these regions, nine candidate genes, which included sucrose phosphate synthase (SPS), trehalose 6-phosphate phosphatase (TPP), raffinose synthase (RS), sucrose synthase (SS), UDP-glucose 4-epimerase (UGE), sugar transport, and sugars will eventually be exported transporter (SWEET), were identified as being related to sugar metabolism and accumulation (Table 4). Quantitative real-time PCR was used to estimate the relative expression of 11 candidate genes at the stage of SSC rapid accumulation (Figure 1c). The results showed that five candidate genes had different expression levels between the two parents at the stage of 20DAP, 25DAP, and 30DAP (Figure 7). The expression levels of MELO21359, MELO06759, MELO22262, and MELO06810 were significantly higher in 1228 than 1214 at three stages. MELO18502 had higher expression levels in 1228 than 1214 at 20DAP and 25DAP, while no significant difference was observed at 30DAP. MELO22262 and MELO06810 shared a pattern of higher expression level at 20DAP, and decreased expression at 25DAP and 30DAP. MELO21359 had the higher expression level at 30DAP, while MELO06759 had the higher expression levels at 25DAP and 30DAP. These five genes with different expression levels between two parents are the presumed candidates controlling SSC in melon.

4. Discussion

4.1. Phenotypic Data over Multiple Years and High-Density Markers Improved the Accuracy of QTL Mapping and Gwas

Linkage mapping based on biparental segregating populations and GWAS based on a natural population are complementary approaches to detect genetic variations for important agronomic traits. Linkage mapping can only evaluate limited variations between two parents in a single population. However, Association analysis can investigate all the allelic variations in natural populations. For QTL mapping and GWAS, phenotypic data over multiple years could reduce environmental factors to identify genetic loci that consistently affect the trait. The sweetness of melon has a significant effect on its quality and consumer appeal, and SSC is an important index of melon sweetness. In this study, phenotypic data of the two melon populations were investigated for three years. In addition, two cultivars with different SSC levels were utilized to understand the accumulation process. We found that SSC accumulated quickly from 20 DAP to 30 DAP, and a similar accumulation pattern for sugar has been observed in previous studies [50,61]. Higher marker density could also enhance the accuracy of linkage analysis and GWAS. In this study, 1427 bin markers generated via genome resequencing were used to construct a high-resolution genetic linkage map with an average interval between the adjacent markers of 0.88 cM. For GWAS, 1.68 million high-confidence SNPs were used in this study, leading to a high precision of GWAS to investigate the complex traits.

4.2. QTL Intervals and SNP Loci for SSC Identified via Linkage Mapping and GWAS

In this study, 10 QTLs (with an average distance of about 0.88 Mb) related to SSC were identified by linkage mapping using the high-resolution genetic linkage map. Seven QTLs were mapped to the same loci of previously reported QTLs (Table S3) and confirmed the reliability of our study. Based on genome resequencing, 1.68 million marks were utilized to detect the loci related to SSC via GWAS. A total of 211 SNPs significantly associated with SSC were found by GWAS with a significance level of p < 10−6. Compared to the SNP loci detected by GWAS and the QTLs detected in previous studies, 138 SNPs overlapped with the reported QTLs (Table S6). These results indicate that the loci identified in the present study were reliable for candidate genes mining.
In this study, we found that 5 QTLs detected by linkage analysis were co-localized with 13 significant SNPs detected by GWAS. The genomic region overlapped by qSSC4-2, S4_10141036, and S4_10275303 on chromosome 4, and qSSC6-2 and S6_32656543 on chromosome 6 (Table 3, Table S3 and Table S6) were also identified within the QTLs identified in previous studies [4,10]. Due to limited parental materials used in QTL mapping, not all of the genetic variations in melon have been identified, and populations derived from different materials, as in previous studies, could determine new QTLs or confirm previous QTLs [8]. In previous studies, the high sweetness parental materials utilized to construct segregated populations were almost always muskmelon. Saladie et al. [63] showed that the patterns of sugar accumulation were different between the oriental melon and muskmelon during fruit development. In this study, two oriental melon inbreeds with different soluble solids content were utilized to derive a segregated population for QTL mapping, and two stable QTLs, qSSC3-1 and qSSC3-2, identified across multiple environments might be new QTLs associated with SSC.

4.3. Further Analysis of Candidate Genes Differentially Expressed between Parents at Different Fruit Growth Stages

Comparing the QTL intervals and SNP loci identified by linkage analysis and GWAS, the co-detected loci and ±100 kb regions around the significant SNPs within the reported QTLs were considered reliable regions to search for candidate genes. According to the annotation information, the genes related to sugar metabolism and accumulation were considered candidate gene controlling soluble solids content. The relative expression levels of eleven genes predicted to participate in the accumulation of SSC were investigated by qRT-PCR at different fruit stages. Finally, five genes with different expression levels between two parents during different stages were considered putative candidates controlling SSC.
In sweet melon, sucrose, fructose, and glucose are the major soluble sugars, and sucrose is the main factor affecting total sugar content [19,49,64]. The cell wall invertases (CWINV) hydrolyzing sucrose into fructose and glucose play an important role in sugar accumulation [65]. Zanor et al. [66] reported that silencing the cell wall invertase LIN5 could lead to a linear decrease in the total soluble solids content in tomato. Wang et al. [45] showed that a reduction in LIN5 activity caused a decrease of fructose, glucose, and SSC in tomato. A previous study revealed that the expression of CWINV3 was significantly higher in the high sweetness rather than low sweetness oriental melon accessions at different fruit growth stages [61]. In the present study, the same expression pattern of CWINV3 was observed (Figure 7a) and the result indicated that CWINV3 might be involved in SSC accumulation.
Alpha-galactosidases (AGA) that hydrolyze stachyose and raffinose play an important role in mediating source–sink communication during fruit development [62,67,68]. A previous study showed that ClAGA2 controls fruit raffinose hydrolysis in watermelon, and that the sucrose content of fruits increased significantly in the overexpressed plants [67]. Based on the positive selection analysis in Cucurbitaceae genomes including C. lanatus, C. maxima, C. pepo, C. sativus, L. siceraria, and C. melo (HS, DHL92, Payzawat) with A thaliana, S. lycopersicum, M. domestica, and O. sativa, a parallel shift of an amino acid in alpha-galactosidase (MELO 06810) was identified in Cucurbitaceae, which may play an important role in raffinose hydrolysis [50]. In this study, the expression of MELO 06810 was significantly higher in 1228 than 1214 at different fruit growth stages (Figure 7b). Therefore, the alpha-galactosidase (MELO 06810) located in qSSC3-2 is considered a candidate gene related to SSC accumulation.
Raffinose could be hydrolyzed into sucrose and then transported and stored in tissues or cells. Raffinose synthase is a key enzyme for raffinose production, and overexpression of raffinose synthase significantly increased raffinose content [69]. A previous study revealed that the expression of raffinose synthase (MELO18502) was significantly higher in the high sweetness rather than low sweetness oriental melon accessions at different fruit growth stages [61]. In this study, the expression of MELO18502 was significantly higher in 1228 than 1214 at 20DAP and 25DAP(Figure 7c).
Additionally, UDP-glucose 4′-epimerase (UGE) could catalyze the interconversion of UDP-glucose and UDP-galactose, and UDP-glucose may further serve as a substrate for the synthesis of sucrose. A previous study reported that overexpression of UGE could improve the accumulation soluble sugars [70]. In this study, the expression of MELO21359 was significantly higher in 1228 than 1214 at three stages (Figure 7d).
Sugar transport proteins (STP) played an important role on sugar accumulation via the transport of glucose from the apoplast into plant cells [39,71]. A previous study showed that the soluble solids content increased in the overexpressed plants [45]. In this study, the expression of MELO22262 was significantly higher in 1228 than 1214 at three stages (Figure 7e).
Sugar metabolism and accumulation are complex biological processes that require the synergistic action of key genes. In this study, five candidate genes related to SSC were identified by QTL mapping, GWAS, and qRT-PCR. These five genes play important roles in the synthesis and hydrolysis of raffinose, sucrose hydrolysis, sugar transport, and production of substrate for sugar synthesis. High gene expression levels of these five genes may collectively promote sugar accumulation.

5. Conclusions

The sweetness of melon is a significant factor in fruit quality and consumer appeal, and SSC is a key index of melon sweetness. In this study, a recombinant inbred population including 146 lines derived from two oriental melon materials with different sweetness levels and an association population consisting of 213 melon accessions were utilized to detect genomic regions influencing SSC. By comparing the results of linkage analysis and GWAS, we identified seven QTLs and 138 significant SNPs overlapped with the reported QTLs. In addition, we identified two new stable loci on chromosome 3 via QTL mapping and GWAS across multiple environments. Five candidate genes related to SSC were identified by QTL mapping, GWAS, and qRT-PCR, two of which were involved in raffinose and sucrose hydrolysis located in the new stable loci. The other three candidate genes were involved in raffinose synthesis sugar transport, and production of substrate for sugar synthesis. Synergistic action of these five candidate genes may cause the higher SSC. The identification of genomic regions and candidate genes that influence SSC will provide useful information and insights for molecular breeding and offer an opportunity for sugar accumulation research.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cimb45090450/s1, Table S1: The sequence data summary for RILs; Table S2: Sequence information on the genomes of 213 melon accessions; Table S3: QTL analyses of soluble solids content in melon; Table S4: Results of filtered SNP annotation; Table S5: The SNPs with significant association signals (−log10(p) ≥ 6) for SSC; Table S6: GWAS signals overlapped with the reported QTLs; Table S7: Annotations of the genes underlying qSSC3-1, and qSSC3-2 in melon; Table S8: Primer sequences for qRT-PCR analysis; Table S9: Phenotypic data of SSC for RILs; Table S10: Phenotypic data of SSC for the accessions.

Author Contributions

Conceptualization, H.Y. and K.W.; Formal analysis, H.Y. and M.W.; Validation, H.Y., M.W. and H.Z.; Investigation, H.Y., L.F. and X.W.; writing—original draft preparation, H.Y.; Writing—Review and editing, H.Y., M.W. and K.W.; visualization, H.Y. and L.F.; Supervision, K.W.; Project administration, M.W. and K.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of Jiangsu Province, China (BK20201205), the National Natural Science Foundation of China (31902037), Jiangsu Province Seed Industry Revitalization Project (JBGS [2021]075), and the Earmarked Fund for Jiangsu Agricultural Industry Technology System (JATS [2022]147).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets analyzed in this study are available in the manuscript text and Supplementary Materials. The raw sequence data reported in this paper have been deposited in the Genome Sequence Archive at the National Genomics Data Center, China National Center for Bioinformation/Beijing Institute of Genomics, Chinese Academy of Sciences, under accession number CRA011939 and are publicly accessible at https://ngdc.cncb.ac.cn/gsa (released on 24 July 2023).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jeffrey, C. A review of the Cucurbitaceae. Bot. J. Linn. Soc. 1980, 81, 233–247. [Google Scholar] [CrossRef]
  2. Pitrat, M.; Hanelt, P.; Hammer, K. Some Comments on Infraspecific Classification of Cultivars of Melon. Acta Hortic. 2000, 510, 29–36. [Google Scholar] [CrossRef]
  3. Monforte, A.J.; Oliver, M.; Gonzalo, M.J.; Alvarez, J.M.; Dolcet-Sanjuan, R.; Arús, P. Identification of quantitative trait loci involved in fruit quality traits in melon (Cucumis melo L.). Theor. Appl. Genet. 2003, 108, 750–758. [Google Scholar] [CrossRef] [PubMed]
  4. Eduardo, I.; Arús, P.; Monforte, A.J.; Obando, J.; Fernández-Trujillo, J.P.; Martínez, J.A.; Alarcón, A.L.; Álvarez, J.M.; 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. [Google Scholar] [CrossRef]
  5. Paris, M.K.; Zalapa, J.E.; McCreight, J.D.; Staub, J.E. Genetic dissection of fruit quality components in melon (Cucumis melo L.) using a RIL population derived from exotic × elite US Western Shipping germplasm. Mol. Breed. 2008, 22, 405–419. [Google Scholar] [CrossRef]
  6. Harel-Beja, R.; Tzuri, G.; Portnoy, V.; Lotan-Pompan, M.; Lev, S.; Cohen, S.; Dai, N.; Yeselson, L.; Meir, A.; Libhaber, S.E.; et al. A genetic map of melon highly enriched with fruit quality QTLs and EST markers, including sugar and carotenoid metabolism genes. Theor. Appl. Genet. 2010, 121, 511–533. [Google Scholar] [CrossRef]
  7. Diaz, A.; Fergany, M.; Formisano, G.; Ziarsolo, P.; Blanca, J.; Fei, Z.; Staub, J.E.; Zalapa, J.E.; Cuevas, H.E.; Dace, G.; et al. A consensus linkage map for molecular markers and Quantitative Trait Loci associated with economically important traits in melon (Cucumis melo L.). BMC Plant Biol. 2011, 11, 111. [Google Scholar] [CrossRef]
  8. Ramamurthy, R.K.; Waters, B.M. Identification of fruit quality and morphology QTLs in melon (Cucumis melo) using a population derived from flexuosus and cantalupensis botanical groups. Euphytica 2015, 204, 163–177. [Google Scholar] [CrossRef]
  9. Baloch, A.; Baloch, A.W.; Liu, S.; Gao, P.; Baloch, M.J.; Wang, X.; Davis, A.; Ali, M.; Luan, F. Linkage map construction and qtl analysis of fruit traits in melon (Cucumis melo L.) based on caps markers. Pak. J. Bot. 2016, 48, 1579–1584. [Google Scholar]
  10. Argyris, J.M.; Díaz, A.; Ruggieri, V.; Fernández, M.; Jahrmann, T.; Gibon, Y.; Picó, B.; Martín-Hernández, A.M.; Monforte, A.J.; Garcia-Mas, J. QTL Analyses in Multiple Populations Employed for the Fine Mapping and Identification of Candidate Genes at a Locus Affecting Sugar Accumulation in Melon (Cucumis melo L.). Front. Plant Sci. 2017, 8, 1679. [Google Scholar] [CrossRef]
  11. Amanullah, S.; Liu, S.; Gao, P.; Zhu, Z.; Zhu, Q.; Fan, C.; Luan, F. QTL mapping for melon (Cucumis melo L.) fruit traits by assembling and utilization of novel SNPs based CAPS markers. Sci. Hortic. 2018, 236, 18–29. [Google Scholar] [CrossRef]
  12. Galpaz, N.; Gonda, I.; Shem-Tov, D.; Barad, O.; Tzuri, G.; Lev, S.; Fei, Z.; Xu, Y.; Mao, L.; Jiao, C.; et al. Deciphering genetic factors that determine melon fruit-quality traits using RNA -Seq-based high-resolution QTL and eQTL mapping. Plant J. 2018, 94, 169–191. [Google Scholar] [CrossRef] [PubMed]
  13. Pereira, L.; Ruggieri, V.; Pérez, S.; Alexiou, K.G.; Fernández, M.; Jahrmann, T.; Pujol, M.; Garcia-Mas, J. QTL mapping of melon fruit quality traits using a high-density GBS-based genetic map. BMC Plant Biol. 2018, 18, 324. [Google Scholar] [CrossRef] [PubMed]
  14. Pereira, L.; Domingo, M.S.; Argyris, J.; Mayobre, C.; Valverde, L.; Martín-Hernández, A.M.; Pujol, M.; Garcia-Mas, J. A novel introgression line collection to unravel the genetics of climacteric ripening and fruit quality in melon. Sci. Rep. 2021, 11, 11364. [Google Scholar] [CrossRef] [PubMed]
  15. Gur, A.; Tzuri, G.; Meir, A.; Sa’ar, U.; Portnoy, V.; Katzir, N.; Schaffer, A.A.; Li, L.; Burger, J.; Tadmor, Y. Genome-Wide Linkage-Disequilibrium Mapping to the Candidate Gene Level in Melon (Cucumis melo). Sci. Rep. 2017, 7, 9770. [Google Scholar] [CrossRef] [PubMed]
  16. Zhao, G.; Lian, Q.; Zhang, Z.; Fu, Q.; He, Y.; Ma, S.; Ruggieri, V.; Monforte, A.J.; Wang, P.; Julca, I.; et al. A comprehensive genome variation map of melon identifies multiple domestication events and loci influencing agronomic traits. Nat. Genet. 2019, 51, 1607–1615. [Google Scholar] [CrossRef]
  17. Liu, S.; Gao, P.; Zhu, Q.; Zhu, Z.; Liu, H.; Wang, X.; Weng, Y.; Gao, M.; Luan, F. Resequencing of 297 melon accessions reveals the genomic history of improvement and loci related to fruit traits in melon. Plant Biotechnol. J. 2020, 18, 2545–2558. [Google Scholar] [CrossRef]
  18. Pavan, S.; Marcotrigiano, A.; Ciani, E.; Mazzeo, R.; Zonno, V.; Ruggieri, V.; Lotti, C.; Ricciardi, L. Genotyping-by-sequencing of a melon (Cucumis melo L.) germplasm collection from a secondary center of diversity highlights patterns of genetic variation and genomic features of different gene pools. BMC Genom. 2017, 18, 59. [Google Scholar] [CrossRef]
  19. Oren, E.; Dafna, A.; Tzuri, G.; Halperin, I.; Isaacson, T.; Elkabetz, M.; Meir, A.; Saar, U.; Ohali, S.; La, T.; et al. Pan-genome and multi-parental framework for high-resolution trait dissection in melon (Cucumis melo). Plant J. 2022, 112, 1525–1542. [Google Scholar] [CrossRef]
  20. Tripodi, P.; Soler, S.; Campanelli, G.; Díez, M.J.; Esposito, S.; Sestili, S.; Figàs, M.R.; Leteo, F.; Casanova, C.; Platani, C.; et al. Genome wide association mapping for agronomic, fruit quality, and root architectural traits in tomato under organic farming conditions. BMC Plant Biol. 2021, 21, 481. [Google Scholar] [CrossRef]
  21. Wang, Y.; Shi, C.; Ge, P.; Li, F.; Zhu, L.; Wang, Y.; Tao, J.; Zhang, X.; Dong, H.; Gai, W.; et al. A 21-bp InDel in the promoter of STP1 selected during tomato improvement accounts for soluble solid content in fruits. Hortic. Res. 2023, 10, uhad009. [Google Scholar] [CrossRef] [PubMed]
  22. Ferik, F.; Ates, D.; Ercisli, S.; Erdogan, A.; Orhan, E.; Tanyolac, M.B. Genome-wide association links candidate genes to fruit firmness, fruit flesh color, flowering time, and soluble solid content in apricot (Prunus armeniaca L.). Mol. Biol. Rep. 2021, 49, 5283–5291. [Google Scholar] [CrossRef] [PubMed]
  23. Cao, K.; Zhou, Z.; Wang, Q.; Guo, J.; Zhao, P.; Zhu, G.; Fang, W.; Chen, C.; Wang, X.; Wang, X.; et al. Genome-wide association study of 12 agronomic traits in peach. Nat. Commun. 2016, 7, 13246. [Google Scholar] [CrossRef]
  24. Guo, S.; Zhao, S.; Sun, H.; Wang, X.; Wu, S.; Lin, T.; Ren, Y.; Gao, L.; Deng, Y.; Zhang, J.; et al. Resequencing of 414 cultivated and wild watermelon accessions identifies selection for fruit quality traits. Nat. Genet. 2019, 51, 1616–1623. [Google Scholar] [CrossRef]
  25. Leida, C.; Moser, C.; Esteras, C.; Sulpice, R.; Lunn, J.E.; De Langen, F.; Monforte, A.J.; Picó, B. Variability of candidate genes, genetic structure and association with sugar accumulation and climacteric behavior in a broad germplasm collection of melon (Cucumis melo L.). BMC Genet. 2015, 16, 28. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, X.; Ando, K.; Wu, S.; Reddy, U.K.; Tamang, P.; Bao, K.; Hammar, S.A.; Grumet, R.; McCreight, J.D.; Fei, Z. Genetic characterization of melon accessions in the U.S. National Plant Germplasm System and construction of a melon core collection. Mol. Hortic. 2021, 1, 11. [Google Scholar] [CrossRef]
  27. Yang, X.; Yan, J.; Zheng, Y.; Yu, J.; Li, J. Reviews of association analysis for quantitative traits in plants. Acta Agron. Sin. 2007, 33, 523. [Google Scholar]
  28. Yu, J.; Buckler, E.S. Genetic association mapping and genome organization of maize. Curr. Opin. Biotechnol. 2006, 17, 155–160. [Google Scholar] [CrossRef]
  29. Wei, Z.; Yuan, Q.; Lin, H.; Li, X.; Zhang, C.; Gao, H.; Zhang, B.; He, H.; Liu, T.; Jie, Z.; et al. Linkage analysis, GWAS, transcriptome analysis to identify candidate genes for rice seedlings in response to high temperature stress. BMC Plant Biol. 2021, 21, 85. [Google Scholar] [CrossRef]
  30. Xu, J.; Ni, Z.; Chen, F.; Fu, X.; Yu, F. Integrated linkage mapping and genome-wide association study to dissect the genetic basis of zinc deficiency tolerance in maize at seedling stage. Crop. J. 2022, 10, 1807–1818. [Google Scholar] [CrossRef]
  31. Hu, D.; Li, X.; Yang, Z.; Liu, S.; Hao, D.; Chao, M.; Zhang, J.; Yang, H.; Su, X.; Jiang, M.; et al. Downregulation of a gibberellin 3β-hydroxylase enhances photosynthesis and increases seed yield in soybean. New Phytol. 2022, 235, 502–517. [Google Scholar] [CrossRef] [PubMed]
  32. Nimmakayala, P.; Tomason, Y.R.; Abburi, V.L.; Alvarado, A.; Saminathan, T.; Vajja, V.G.; Salazar, G.; Panicker, G.K.; Levi, A.; Wechter, W.P.; et al. Genome-Wide Differentiation of Various Melon Horticultural Groups for Use in GWAS for Fruit Firmness and Construction of a High Resolution Genetic Map. Front. Plant Sci. 2016, 7, 1437. [Google Scholar] [CrossRef] [PubMed]
  33. Perpiñá, G.; Esteras, C.; Gibon, Y.; Monforte, A.J.; Picó, B. A new genomic library of melon introgression lines in a cantaloupe genetic background for dissecting desirable agronomical traits. BMC Plant Biol. 2016, 16, 154. [Google Scholar] [CrossRef]
  34. Oren, E.; Tzuri, G.; Dafna, A.; Meir, A.; Kumar, R.; Katzir, N.; Elkind, Y.; Freilich, S.; Schaffer, A.A.; Tadmor, Y.; et al. High-density NGS-based map construction and genetic dissection of fruit shape and rind netting in Cucumis melo. Theor. Appl. Genet. 2020, 133, 1927–1945. [Google Scholar] [CrossRef] [PubMed]
  35. Du, X.; Liu, H.; Zhu, Z.; Liu, S.; Song, Z.; Xia, L.; Zhao, J.; Luan, F.; Liu, S. Identification of Candidate Chromosome Region Related to Melon (Cucumis melo L.) Fruit Surface Groove Trait Through Biparental Genetic Mapping and Genome-Wide Association Study. Front. Plant Sci. 2022, 13, 828287. [Google Scholar] [CrossRef]
  36. Zhang, C.; Yu, X.; Ayre, B.G.; Turgeon, R. The Origin and Composition of Cucurbit “Phloem” Exudate. Plant Physiol. 2012, 158, 1873–1882. [Google Scholar] [CrossRef] [PubMed]
  37. Ko, H.-Y.; Ho, L.-H.; Neuhaus, H.E.; Guo, W.-J. Transporter SlSWEET15 unloads sucrose from phloem and seed coat for fruit and seed development in tomato. Plant Physiol. 2021, 187, 2230–2245. [Google Scholar] [CrossRef]
  38. Li, Y.; Liu, H.; Yao, X.; Wang, J.; Feng, S.; Sun, L.; Ma, S.; Xu, K.; Chen, L.-Q.; Sui, X. Hexose transporter CsSWEET7a in cucumber mediates phloem unloading in companion cells for fruit development. Plant Physiol. 2021, 186, 640–654. [Google Scholar] [CrossRef]
  39. Wen, S.; Neuhaus, H.E.; Cheng, J.; Bie, Z. Contributions of sugar transporters to crop yield and fruit quality. J. Exp. Bot. 2022, 73, 2275–2289. [Google Scholar] [CrossRef]
  40. Ren, Y.; Sun, H.; Zong, M.; Guo, S.; Ren, Z.; Zhao, J.; Li, M.; Zhang, J.; Tian, S.; Wang, J.; et al. Localization shift of a sugar transporter contributes to phloem unloading in sweet watermelons. New Phytol. 2020, 227, 1858–1871. [Google Scholar] [CrossRef]
  41. Peng, Q.; Cai, Y.; Lai, E.; Nakamura, M.; Liao, L.; Zheng, B.; Ogutu, C.; Cherono, S.; Han, Y. The sucrose transporter MdSUT4.1 participates in the regulation of fruit sugar accumulation in apple. BMC Plant Biol. 2020, 20, 191. [Google Scholar] [CrossRef] [PubMed]
  42. Ren, Y.; Guo, S.; Zhang, J.; He, H.; Sun, H.; Tian, S.; Gong, G.; Zhang, H.; Levi, A.; Tadmor, Y.; et al. A Tonoplast Sugar Transporter Underlies a Sugar Accumulation QTL in Watermelon. Plant Physiol. 2017, 176, 836–850. [Google Scholar] [CrossRef] [PubMed]
  43. Cheng, J.; Wen, S.; Xiao, S.; Lu, B.; Ma, M.; Bie, Z. Overexpression of the tonoplast sugar transporter CmTST2 in melon fruit increases sugar accumulation. J. Exp. Bot. 2017, 69, 511–523. [Google Scholar] [CrossRef] [PubMed]
  44. Huang, W.; Hu, B.; Liu, J.; Zhou, Y.; Liu, S. Identification and Characterization of Tonoplast Sugar Transporter (TST) Gene Family in Cucumber. Hortic. Plant J. 2020, 6, 145–157. [Google Scholar] [CrossRef]
  45. Wang, Z.; Wei, X.; Yang, J.; Li, H.; Ma, B.; Zhang, K.; Zhang, Y.; Cheng, L.; Ma, F.; Li, M. Heterologous expression of the apple hexose transporter Md HT 2.2 altered sugar concentration with increasing cell wall invertase activity in tomato fruit. Plant Biotechnol. J. 2019, 18, 540–552. [Google Scholar] [CrossRef] [PubMed]
  46. Zhu, L.; Li, B.; Wu, L.; Li, H.; Wang, Z.; Wei, X.; Ma, B.; Zhang, Y.; Ma, F.; Ruan, Y.-L.; et al. MdERDL6-mediated glucose efflux to the cytosol promotes sugar accumulation in the vacuole through up-regulating TSTs in apple and tomato. Proc. Natl. Acad. Sci. USA 2020, 118, e2022788118. [Google Scholar] [CrossRef]
  47. Lunn, J.E.; Delorge, I.; Figueroa, C.M.; Van Dijck, P.; Stitt, M. Trehalose metabolism in plants. Plant J. 2014, 79, 544–567. [Google Scholar] [CrossRef]
  48. Zhang, Z.; Deng, Y.; Song, X.; Miao, M. Trehalose-6-phosphate and SNF1-related protein kinase 1 are involved in the first-fruit inhibition of cucumber. J. Plant Physiol. 2014, 177, 110–120. [Google Scholar] [CrossRef]
  49. Dai, N.; Cohen, S.; Portnoy, V.; Tzuri, G.; Harel-Beja, R.; Pompan-Lotan, M.; Carmi, N.; Zhang, G.; Diber, A.; Pollock, S.; et al. Metabolism of soluble sugars in developing melon fruit: A global transcriptional view of the metabolic transition to sucrose accumulation. Plant Mol. Biol. 2011, 76, 1–18. [Google Scholar] [CrossRef]
  50. Yang, J.; Deng, G.; Lian, J.; Garraway, J.; Niu, Y.; Hu, Z.; Yu, J.; Zhang, M. The Chromosome-Scale Genome of Melon Dissects Genetic Architecture of Important Agronomic Traits. iScience 2020, 23, 101422. [Google Scholar] [CrossRef]
  51. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows—Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed]
  52. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. Genome Project Data Processing S: The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed]
  53. McKenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M.; et al. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef]
  54. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef] [PubMed]
  55. Huang, X.; Feng, Q.; Qian, Q.; Zhao, Q.; Wang, L.; Wang, A.; Guan, J.; Fan, D.; Weng, Q.; Huang, T.; et al. High-throughput genotyping by whole-genome resequencing. Genome Res. 2009, 19, 1068–1076. [Google Scholar] [CrossRef] [PubMed]
  56. Tang, H.; Zhang, X.; Miao, C.; Zhang, J.; Ming, R.; Schnable, J.C.; Schnable, P.S.; Lyons, E.; Lu, J. ALLMAPS: Robust scaffold ordering based on multiple maps. Genome Biol. 2015, 16, 3. [Google Scholar] [CrossRef] [PubMed]
  57. Alexander, D.H.; Novembre, J.; Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009, 19, 1655–1664. [Google Scholar] [CrossRef]
  58. Zhang, C.; Dong, S.-S.; Xu, J.-Y.; He, W.-M.; Yang, T.-L. PopLDdecay: A fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics 2019, 35, 1786–1788. [Google Scholar] [CrossRef]
  59. Zhou, X.; Stephens, M. Efficient multivariate linear mixed model algorithms for genome-wide association studies. Nat. Methods 2014, 11, 407–409. [Google Scholar] [CrossRef]
  60. Lex, A.; Gehlenborg, N.; Strobelt, H.; Vuillemot, R.; Pfister, H. UpSet: Visualization of Intersecting Sets. IEEE Trans. Vis. Comput. Graph. 2014, 20, 1983–1992. [Google Scholar] [CrossRef]
  61. Cheng, H.; Kong, W.; Tang, T.; Ren, K.; Zhang, K.; Wei, H.; Lin, T. Identification of Key Gene Networks Controlling Soluble Sugar and Organic Acid Metabolism During Oriental Melon Fruit Development by Integrated Analysis of Metabolic and Transcriptomic Analyses. Front. Plant Sci. 2022, 13, 830517. [Google Scholar] [CrossRef]
  62. Ohkawa, W.; Kanayama, Y.; Daibo, N.; Sato, T.; Nishiyama, M.; Kanahama, K. Metabolic process of the 14C-sugars on the translocation pathways of cucumber plants. Sci. Hortic. 2010, 124, 46–50. [Google Scholar] [CrossRef]
  63. Saladié, M.; Cañizares, J.; Phillips, M.A.; Rodriguez-Concepcion, M.; Larrigaudière, C.; Gibon, Y.; Stitt, M.; Lunn, J.E.; Garcia-Mas, J. Comparative transcriptional profiling analysis of developing melon (Cucumis melo L.) fruit from climacteric and non-climacteric varieties. BMC Genom. 2015, 16, 440. [Google Scholar] [CrossRef]
  64. Shin, A.-Y.; Kim, Y.-M.; Koo, N.; Lee, S.M.; Nahm, S.; Kwon, S.-Y. Transcriptome analysis of the oriental melon (Cucumis melo L. var. makuwa) during fruit development. PeerJ 2017, 5, e2834. [Google Scholar] [CrossRef] [PubMed]
  65. Jin, Y.; Ni, D.-A.; Ruan, Y.-L. Posttranslational Elevation of Cell Wall Invertase Activity by Silencing Its Inhibitor in Tomato Delays Leaf Senescence and Increases Seed Weight and Fruit Hexose Level. Plant Cell 2009, 21, 2072–2089. [Google Scholar] [CrossRef]
  66. Zanor, M.I.; Osorio, S.; Nunes-Nesi, A.; Carrari, F.; Lohse, M.; Usadel, B.; Kühn, C.; Bleiss, W.; Giavalisco, P.; Willmitzer, L.; et al. RNA Interference of LIN5 in Tomato Confirms Its Role in Controlling Brix Content, Uncovers the Influence of Sugars on the Levels of Fruit Hormones, and Demonstrates the Importance of Sucrose Cleavage for Normal Fruit Development and Fertility. Plant Physiol. 2009, 150, 1204–1218. [Google Scholar] [CrossRef] [PubMed]
  67. Ren, Y.; Li, M.; Guo, S.; Sun, H.; Zhao, J.; Zhang, J.; Liu, G.; He, H.; Tian, S.; Yu, Y.; et al. Evolutionary gain of oligosaccharide hydrolysis and sugar transport enhanced carbohydrate partitioning in sweet watermelon fruits. Plant Cell 2021, 33, 1554–1573. [Google Scholar] [CrossRef]
  68. Liu, H.; Liu, X.; Zhao, Y.; Nie, J.; Yao, X.; Lv, L.; Yang, J.; Ma, N.; Guo, Y.; Li, Y.; et al. Alkaline α-galactosidase 2 (CsAGA2) plays a pivotal role in mediating source–sink communication in cucumber. Plant Physiol. 2022, 189, 1501–1518. [Google Scholar] [CrossRef] [PubMed]
  69. Liu, Y.; Li, T.; Zhang, C.; Zhang, W.; Deng, N.; Dirk, L.M.A.; Downie, A.B.; Zhao, T. Raffinose positively regulates maize drought tolerance by reducing leaf transpiration. Plant J. 2023, 114, 55–67. [Google Scholar] [CrossRef]
  70. Tang, Y.; Wang, M.; Cao, L.; Dang, Z.; Ruan, N.; Wang, Y.; Huang, Y.; Wu, J.; Zhang, M.; Xu, Z.; et al. OsUGE3-mediated cell wall polysaccharides accumulation improves biomass production, mechanical strength, and salt tolerance. Plant Cell Environ. 2022, 45, 2492–2507. [Google Scholar] [CrossRef]
  71. Bavnhøj, L.; Paulsen, P.A.; Flores-Canales, J.C.; Schiøtt, B.; Pedersen, B.P. Molecular mechanism of sugar transport in plants unveiled by structures of glucose/H+ symporter STP10. Nat. Plants 2021, 7, 1409–1419. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Frequency distribution and changes in soluble solids content. (a) Frequency distribution of SSC for 213 melon accessions in four environments. (b) Frequency distribution of SSC for RILs in three environments. SPR represents spring and AUT represents autumn. Arrows indicate the soluble solids content for the parents, 1214(P1) and 1228(P2). (c) SSC in the fruit during the five ripening stages. ** Significant at p ≤ 0.01.
Figure 1. Frequency distribution and changes in soluble solids content. (a) Frequency distribution of SSC for 213 melon accessions in four environments. (b) Frequency distribution of SSC for RILs in three environments. SPR represents spring and AUT represents autumn. Arrows indicate the soluble solids content for the parents, 1214(P1) and 1228(P2). (c) SSC in the fruit during the five ripening stages. ** Significant at p ≤ 0.01.
Cimb 45 00450 g001
Figure 2. Construction of the linkage map (a) and collinearity analysis between the linkage map and melon genome ‘HS_ZJU’ (b).
Figure 2. Construction of the linkage map (a) and collinearity analysis between the linkage map and melon genome ‘HS_ZJU’ (b).
Cimb 45 00450 g002
Figure 3. Distribution of QTLs related to SSC on the melon chromosomes. Gray bars represent the chromosomes and colored bars represent the QTLs for SSC. Bar length indicates the physical distance of QTLs.
Figure 3. Distribution of QTLs related to SSC on the melon chromosomes. Gray bars represent the chromosomes and colored bars represent the QTLs for SSC. Bar length indicates the physical distance of QTLs.
Cimb 45 00450 g003
Figure 4. Distribution of the SNPs and the population structure of 213 melon accessions. (a) Distribution of SNPs and nucleotide diversity across the melon genome ‘HS_ZJU’. Different colors indicate the number of SNPs within the 0.1 Mb window size. (b) Genome-wide average LD decay for 213 melon accessions. (c) Phylogenetic tree of 213 melon accessions. (d) Population structure of melon accessions estimated by ADMIXTURE with k = 2.
Figure 4. Distribution of the SNPs and the population structure of 213 melon accessions. (a) Distribution of SNPs and nucleotide diversity across the melon genome ‘HS_ZJU’. Different colors indicate the number of SNPs within the 0.1 Mb window size. (b) Genome-wide average LD decay for 213 melon accessions. (c) Phylogenetic tree of 213 melon accessions. (d) Population structure of melon accessions estimated by ADMIXTURE with k = 2.
Cimb 45 00450 g004
Figure 5. GWAS signals for SSC estimated across four environments and BLUP. Manhattan plot and quantile–quantile (Q–Q) plot estimated in (a) spring 2020 (E1), (b) spring 2021 (E2), (c) autumn 2021 (E3), (d) spring 2022 (E4), and (e) BLUP. The gray lines on the Y-axis denote the significance threshold (−log10 p > 6).
Figure 5. GWAS signals for SSC estimated across four environments and BLUP. Manhattan plot and quantile–quantile (Q–Q) plot estimated in (a) spring 2020 (E1), (b) spring 2021 (E2), (c) autumn 2021 (E3), (d) spring 2022 (E4), and (e) BLUP. The gray lines on the Y-axis denote the significance threshold (−log10 p > 6).
Cimb 45 00450 g005aCimb 45 00450 g005b
Figure 6. Intersection of significant SNPs across four different environments (E1, E2, E3, and E4) and BLUP.
Figure 6. Intersection of significant SNPs across four different environments (E1, E2, E3, and E4) and BLUP.
Cimb 45 00450 g006
Figure 7. Relative expression levels of five candidate genes. (a) MELO06759; (b) MELO06810; (c) MELO18502; (d) MELO21359; and (e) MELO22262; ** Significant at p ≤ 0.01.
Figure 7. Relative expression levels of five candidate genes. (a) MELO06759; (b) MELO06810; (c) MELO18502; (d) MELO21359; and (e) MELO22262; ** Significant at p ≤ 0.01.
Cimb 45 00450 g007
Table 1. Descriptive statistics and variance analysis for SSC of the two panels in multiple environments.
Table 1. Descriptive statistics and variance analysis for SSC of the two panels in multiple environments.
PopulationEnvironment aParentsMean ± SE (%)SD bRange (%)CV c (%)SkewnessKurtosisG dE dG × E dH2 e
12141228
RILs2020 SPR7.4312.4010.43 ± 0.131.565.67–14.5014.96−0.150.13******0.65
2021 SPR9.2613.4811.01 ± 0.131.555.83–14.0714.07−0.350.46
2022 SPR9.6214.3511.29 ± 0.131.467.15–14.8012.93−0.34−0.01
Accessions2020 SPR--10.01 ± 0.202.623.00–17.3026.17−0.39−0.17******0.72
2021 SPR--9.86 ± 0.202.873.00–17.4329.01−0.17−0.57
2021 AUT--8.97 ± 0.212.713.55–15.5030.210.16−0.57
2022 SPR--10.68 ± 0.233.123.80–18.9829.21−0.03−0.75
a SPR represents spring and AUT represents autumn; b SD represents the standard deviation; c CV represents the coefficient of variation; d G, E, and G × E represent the effect for the genotype, environment, and genotype × environment interaction, respectively; e broad-sense heritability. ** Significant at p ≤ 0.01.
Table 2. Information of the novel constructed genetic linkage map for melon.
Table 2. Information of the novel constructed genetic linkage map for melon.
LGsNumber of
Bins
Map Length (cM)Average
Distance (cM)
Max Gap
(cM)
<5 cM>5 cM
LG01148133.750.904.271480
LG0210984.730.783.001090
LG03119105.390.893.841190
LG049690.160.943.84960
LG05152126.170.833.421520
LG06158128.210.814.271580
LG07114105.840.933.841140
LG08109101.160.934.271090
LG09118105.970.903.001180
LG104358.461.3614.26421
LG11149133.920.903.841490
LG1211280.590.723.841120
total14271254.340.8814.2614261
Table 3. Co-detected loci by QTL mapping and GWAS.
Table 3. Co-detected loci by QTL mapping and GWAS.
Chr.EnvironmentLODPosition Interval (cM)Physical Position (bp)SNP Loci−log10(p)
2020SPR2021SPR2021AUT2022SPRBLUP
32020SPR4.595.2–104.328,509,689–28,965,70528,672,844 6.01
32022SPR5.593.7–99.528,446,470–28,845,99828,694,3926.276.05 6.27
3BLUP5.795.6–102.928,509,689–28,947,99828,711,4966.47 6.51
28,714,701 6.19
28,724,7766.08 6.67
3BLUP5.2102.9–105.128,947,998–29,341,84929,267,7826.68 6.048.497.65
29,283,426 6.14
29,289,1576.22
29,330,7887.77 8.337.887.88
42021SPR3.520.9–28.62,334,788–3,020,2262,491,4346.74 6.80
4BLUP3.451.7–54.28,961,531–10,360,27910,141,0367.13 6.63
10,275,303 7.06
62021SPR3.595.5–10132,217,543–32,935,48832,656,543 6.56
Table 4. Candidate genes related to SSC within the QTL/SNP regions.
Table 4. Candidate genes related to SSC within the QTL/SNP regions.
GeneChromosomeStartEndFunctional Annotation
MELO03698chr214,397,61514,406,002Sucrose-phosphate synthase
MELO06759chr328,826,32928,829,358Beta-fructofuranosidase, insoluble isoenzyme CWINV3
MELO06810chr329,119,43629,126,262Alpha-galactosidase
MELO12562chr66,109,5126,114,363Trehalose 6-phosphate phosphatase
MELO18502chr816,852,49716,859,087Raffinose synthase
MELO20547chr914,794,45514,800,212Sucrose synthase
MELO21359chr923,772,89323,777,340UDP-glucose 4′-epimerase
MELO22262chr105,145,5185,147,793Sugar transport protein 10-like
MELO27354chr1220,598,70620,606,584Trehalose-6-phosphate synthase
MELO27634chr1223,004,02423,006,015Sugar transporter SWEET12-like
MELO27635chr1223,013,80023,015,623Sugar transporter SWEET12-like
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yan, H.; Wang, K.; Wang, M.; Feng, L.; Zhang, H.; Wei, X. QTL Mapping and Genome-Wide Association Study Reveal Genetic Loci and Candidate Genes Related to Soluble Solids Content in Melon. Curr. Issues Mol. Biol. 2023, 45, 7110-7129. https://doi.org/10.3390/cimb45090450

AMA Style

Yan H, Wang K, Wang M, Feng L, Zhang H, Wei X. QTL Mapping and Genome-Wide Association Study Reveal Genetic Loci and Candidate Genes Related to Soluble Solids Content in Melon. Current Issues in Molecular Biology. 2023; 45(9):7110-7129. https://doi.org/10.3390/cimb45090450

Chicago/Turabian Style

Yan, Honglang, Kang Wang, Manman Wang, Lulu Feng, Huimin Zhang, and Xiaoyun Wei. 2023. "QTL Mapping and Genome-Wide Association Study Reveal Genetic Loci and Candidate Genes Related to Soluble Solids Content in Melon" Current Issues in Molecular Biology 45, no. 9: 7110-7129. https://doi.org/10.3390/cimb45090450

APA Style

Yan, H., Wang, K., Wang, M., Feng, L., Zhang, H., & Wei, X. (2023). QTL Mapping and Genome-Wide Association Study Reveal Genetic Loci and Candidate Genes Related to Soluble Solids Content in Melon. Current Issues in Molecular Biology, 45(9), 7110-7129. https://doi.org/10.3390/cimb45090450

Article Metrics

Back to TopTop