Comparative transcriptomic analysis of two Cucumis melo var. saccharinus germplasms differing in fruit physical and chemical characteristics

Background Hami melon (Cucumis melo var. saccharinus) is a popular fruit in China because of its excellent taste, which is largely determined by its physicochemical characteristics, including flesh texture, sugar content, aroma, and nutrient composition. However, the mechanisms by which these characteristics are regulated have not yet been determined. In this study, we monitored changes in the fruits of two germplasms that differed in physicochemical characteristics throughout the fruit development period. Results Ripe fruit of the bred variety ‘Guimi’ had significantly higher soluble sugar contents than the fruit of the common variety ‘Yaolong.’ Additionally, differences in fruit shape and color between these two germplasms were observed during development. Comparative transcriptome analysis, conducted to identify regulators and pathways underlying the observed differences at corresponding stages of development, revealed a higher number of differentially expressed genes (DEGs) in Guimi than in Yaolong. Moreover, most DEGs detected during early fruit development in Guimi were associated with cell wall biogenesis. Temporal analysis of the identified DEGs revealed similar trends in the enrichment of downregulated genes in both germplasms, although there were differences in the enrichment trends of upregulated genes. Further analyses revealed trends in differential changes in multiple genes involved in cell wall biogenesis and sugar metabolism during fruit ripening. Conclusions We identified several genes associated with the ripening of Hami melons, which will provide novel insights into the molecular mechanisms underlying the development of fruit characteristics in these melons. Supplementary Information The online version contains supplementary material available at 10.1186/s12870-022-03550-8.


Background
Hami melon (Cucumis melo var. saccharinus) is among the most important diploid crops within the family Cucurbitaceae, the different cultivars of which are characterized by highly variable fruit traits, such as esh color, sugar contents, and shape (Wang et al. 2018) [1]. Production of Hami melons in Guangxi Province, China, started in the middle to late 1990s (Nunez-Palenius et al. 2008). However, given the small differences in day and night temperatures, the Guangxi environment is not particularly conducive to sugar accumulation in Hami melon, as it is di cult to achieve the desirable 14% soluble solid content (Bi et al. 2003). In addition, high temperatures and humidity often contribute to the development of a range of diseases and pest infestations, thereby reducing yields and resulting in heavy economic losses, notably as a consequence of gummy stem blight caused by Mycosphaerella melonis. To address these problems, 'Guimi-12', a new high-quality Guangxi germplasm, characterized by high yields and disease resistance, has been selected and popularized. Guimi-12 differs from the common Hami melon germplasm ('Yaolong') with respect to shape, size, esh texture, sugar content, aroma, and nutrient composition, which are acknowledged to be important considerations in breeding of fruit-bearing plants and are typically the focus of studies seeking to identify and characterize the associated regulatory transcription factors and molecular mechanisms (Guo et al. 2015). At present, however, little is known about the differences between these two Hami melon varieties with respect to gene expression patterns and regulatory pathways.
To date, breeding programs for Hami melon have tended to focus primarily on the selection of traits associated with fruit thickness and sweetness (Yu et al. 2012). Pericarp hardening is generally associated with the development of secondary cell wall structures, although the underlying mechanisms have yet to be ascertained. In this regard, however, it has been established that a transcription network comprising NAC and MYB transcription factors is highly conserved in the regulation of pericarp formation and secondary cell wall development in Arabidopsis (Dardick et al. 2010). With respect to sweetness, it has been established that sucrose is the predominant sugar in ripe melon fruit (Burger et al., 2003;Villanueva et al., 2004). α-Galactosidase (α-Gal) has been identi ed as the initial enzyme in the decomposition of stachyose and its activity has been shown to be closely associated with the accumulation of fructose and glucose during fruit expansion (Alan et al. 2007). Consistently, acidic α-Gal has been found to have higher activity in the sweet tissue of watermelon fruit than that in unsweet tissue, thereby highlighting its signi cance in the accumulation of sugar in watermelon (Colla et al. 2006). Numerous metabolic pathways have been shown to contribute to sugar accumulation in watermelon fruit (Braun et al. 2012), with 62 sugar synthesis enzymes and 76 sugar transporter genes associated with these pathways having been identi ed (Guo et al. 2013), 13 and 14 genes of which, respectively, have been shown to be differentially expressed during fruit development ). In addition, it has also been found that the patterns of sugar accumulation during fruit maturation differ between oriental and other types of melon (Saladie et al. 2015). These observations thus indicate that differences in the genes and pathways associated sugar metabolism probably determine the differing sucrose contents of the fruits, and we accordingly speculate that sugar contents in the fruit of Guimi-12 Hami melon might also differ from those in Yaolong Hami melon.
Transcriptomes are collections of all transcripts produced by a species or speci c types of organs, tissues, and cells, and include the number of transcripts, expressional dynamics at speci c developmental stages, post transcriptional modi cation, and regulation of non-coding RNAs (Ashburner et al. 2000). By conducting transcriptomic studies, we can perform more comprehensive analyses of gene expression, structure, and function in speci c species, and thereby gain a more complete understanding of the molecular mechanism underlying speci c physiological characteristics (Pruitt et al. 2005).
Compared with the whole genome, properties at the transcriptome level tend to be more temporally and spatially discrete, as has been reported in numerous studies (Trapnell et al. 2010;Grassi et al. 2013). Thus, transcriptome sequencing technology can provide a new approach to the study of gene function and physiological characteristics in Hami melon. In this study, we sought to comparatively analyze the transcriptomes of two contrasting Hami melon germplasms, namely, Guimi-12 and Yaolong, throughout the fruit ripening process, with the aim of identifying those genes playing roles in the development and ripening of Hami melon fruit, and thereby provide further insights for identifying key pathways and regulators associated with the physical and chemical characteristics of Hami melon fruit.

Plant materials
The germplasm 'Guimi' was provided by Vegetable Research Institute, Guangxi Academy of Agricultural Sciences, Nanning, China. The germplasm 'Yaolong' was purchased from Fuyou Seedlings Co. LTD, Hainan, China. Both germplasms were cultivated in a greenhouse, in which the temperature was maintained between 25℃ to 30℃ during the day and above 15℃ at night. Flowers were manually pollinated, and three to ve biological replicate fruits were harvested from selected plants at 10,20,30, and 40 DAP. These fruits were frozen in liquid nitrogen and stored at -80℃ until used for analyses.

Determination of physical and chemical characteristics
Among the physical characteristics measured, the size and weight of fruit were measured at 5,10,15,20,25,30,35, and 40 DAP. The Brix value of fruit was determined using an ATAGO PAL-2 refractometer, as previously described (Abbas et al. 2020). Each measurement was conducted for least three replicates.
The soluble sugar contents (total soluble solids, glucose, fructose, and sucrose) were analyzed spectrophotometrically, as described previously (Zhu et al. 2017).
RNA isolation, cDNA library preparation, and RNA-Seq Total RNA was extracted from the pulp of frozen fruit collected at each time point using a TRIzol reagent kit (Invitrogen, Carlsbad, CA, USA) based on the manufacturer's protocol, and the quality of the isolated total RNA was evaluated using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) and RNase-free agarose gel electrophoresis. mRNA was enriched, fragmented, and reverse transcribed to rst-strand cDNA using Oligo(dT), fragmentation buffer, and random primers. Second-strand cDNA was subsequently synthesized using DNA polymerase I, RNase H, dNTPs, and buffer. Finally, 24 cDNA libraries were constructed and sequenced using the Illumina Novaseq6000 platform with pair-end 150bp mode.

Global and differential gene expression analysis
To obtain high-quality clean reads, reads containing adapter and poly-N sequences or low-quality bases were ltered out using fastp (version 0.18.0) (Chen et al., 2018) with default parameters. Thereafter, the rRNA mapped reads were removed following the alignment of clean reads to a ribosomal RNA (rRNA) database using Bowtie2(version2.28) (Langmead and Salzberg, 2012). As a Cucumis melo reference genome, we used the RefSeq assembly accession GCF_000313046.1. The remaining cleaned reads were aligned to the reference genome with HISAT2 (version 2.4) using "-rna-strandness RF" (Kim et al., 2015), and genes were assembled using StringTie (version 1. For analysis of differential gene expression between the Yaolong and Guimi groups, we used DESeq2 software. Differentially expressed genes (DEGs) were identi ed based on the criteria of a false discovery rate (FDR) < 0.05 and an absolute value of log2(fold change) ≥ 1.

Temporal analysis
To determine trends in the changes in gene expression during C. melo fruit development, we performed temporal analysis for the DEGs of each variety using Short Time-series Expression Miner (STEM) software with default parameters (Ernst et al. 2006). Genes were pro led by clustering based on the corresponding P-values and those with a P-value ≤ 0.05 were considered differentially expressed, These DEGs were subjected to GO term and KEGG pathway enrichment analyses using hypergeometric distribution tests. GO terms and KEGG pathways with Q-values equal to or less than 0.05 were considered functional annotations.

Results
Variations in the physical and chemical characteristics of melon fruit during ripening For comparisons of fruits of the two assessed Hami melon varieties, we examined the physiochemical characteristics of fruit collected at different stages of development (Fig. 1). The color of the fruit peel of the two varieties was found to be similar at 5 and 10 days after pollination (DAP). Thereafter, however, the peel of Yaolong fruit turned a darker green compared with that of Guimi from 15 to 20 DAP. With continued growth, Guimi melons developed as medium-sized oval fruits with a yellow skin. In both varieties, the juicy esh underwent a change in color from white to yellow during ripening.
Trends in the changes in physical and chemical characteristics of the fruits, including weight, size, and soluble sugar contents, are shown in Fig. 2. Similar trends in the two varieties were observed with respect to the weight and size of fruit during repining. Interestingly, fruit weight increased rapidly from 35 to 40 DAP, whereas there were only slight increases in size. At 40 DAP, the soluble sugar contents of Guimi fruit had reached a value 12, which was higher than that obtained for Yaolong, although differences in the brix values of the two varieties were found to be non-signi cant.

Global analysis of the RNA-Seq data
To determine gene expression patterns, we performed RNA-seq analysis using the C. melo reference genome. After ltering out rRNAs and low-quality reads, a total of 111 million reads were mapped to the reference genome (Table S1). For these clean reads, we obtained average mapped reads per sample greater that 90%. The one exception in this the regard was the low (74.44%) alignment ratio obtained for Yaolong at 40 DAP. In total, we detected 21,172 expressed genes among the melon fruit samples.
To determine differences in the expression of genes between two close sampling times, we identi ed the DEGs in Guimi and Yaolong based on the threshold criteria of a log2 fold change ≥ 1 and FDR ≤ 0.05 (Table S2). For this purpose, we de ned the earliest time point of paired groups as the control sample for subsequent measurements. The number of DEGs in Guimi was found to be markedly higher than that in Yaolong at 20 and 30 DAP, while DEGs being identi ed in Guimi less than Yaolong (223 vs 95) at 40 DAP ( Fig. 3A). These nding would thus appear to indicate that during fruit development the constituents of Guimi fruit undergo more pronounced changes than those of Yaolong fruit. Over time, however, there was a gradual reduction in the number of DEGs in both melon varieties during the fruit development, indicting a corresponding reduction in the speed of fruit growth and that the fruit was fully mature at 40 DAP. For both Guimi and Yaolong, we also detected perturbations in gene expression at the same time points. The number of genes differentially expressed between Guimi and Yaolong initially increased, reaching a peak at 20 DAP, after which the number declined to 458 at 40 DAP (Fig. 3B).
For the purpose of functional annotation, identi ed DEGs in the two melon varieties were assigned to GO terms and KEGG pathways. We accordingly found that for Yaolong melon, the majority of enriched GO terms between 10 and 20 DAP were classi ed in the category biological process, including 'histone lysine methylation', 'peptidyl-lysine methylation', and 'DNA alkylation', whereas in the case of Guimi melon, the most signi cantly enriched GO terms were 'polysaccharide metabolism', 'cell wall biogenesis', and 'external encapsulating structure organization' (Fig. S1 and S2). For the second comparison (20 DAP vs. 30 DAP), the DEGs in Guimi melon were found to be signi cant enriched in the processes of 'response to acid chemical', 'plant-type cell wall biogenesis', and 'response to chemical', whereas for Yaolong, the three most enriched biological processes were 'hemicellulose metabolism', xylan metabolism', and 'cell wall polysaccharide metabolism'. Comparing between 30 and 40 DAP samples revealed that the DEGs of Guimi melon were notably enriched with respect to 'amine metabolism', 'cellular amine metabolism', and 'jasmonic acid metabolism', whereas those of Yaolong melon were found to be enriched in the cellular component categories 'external encapsulating structure', 'cell periphery', and 'vacuole'. With regards to KEGG pathway annotation, we found that in the 10 to 20 DAP comparison, DEGs in Yaolong were enriched with respect to 'DNA replication', 'ABC transporters', and ' avone and avonol biosynthesis', whereas 'metabolic pathways', 'biosynthesis of secondary metabolites', and 'phenylpropanoid biosynthesis' were found to be enriched with Guimi DEGs (Fig. S3 and S4). Interestingly, for the 20 to 30 DAP comparison, we observed that for both melon varieties, the three most signi cantly DEG-enriched KEGG pathways were 'biosynthesis of secondary metabolites', 'metabolic pathways' and 'phenylpropanoid biosynthesis'. For the third comparison between 30 and 40 DAP, Guimi DEGs were found to be enriched in pathways such as 'valine, leucine, and isoleucine biosynthesis', 'biosynthesis of secondary metabolites' and 'cyanoamino acid metabolism', whereas pathways enriched with Yaolong DEGs included 'biosynthesis of secondary metabolites' and 'metabolic pathways'.
Comparison of trends in temporal gene expression during melon fruit ripening To gain further insights into the change trends in gene expression during fruit development, we clustered a total of 4,731 DEGs from Guimi melon and 3,198 DEGs from Yaolong into 38 pro les using the STEM algorithm. Among these, 2,897 Guimi DEGs were signi cantly clustered into the following six pro les: two upregulated pro les (Pro les 17 and 12), three downregulated pro les (Pro les 0, 7, and 2), and one biphasic expression pattern pro le (Pro le 18) (Fig. 4A). Similarly, 2,217 Yaolong DEGs were classi ed into the following six pro les based on P values ≤ 0.05: two upregulated patterns, one biphasic expression pattern, and three downregulated patterns (Fig. 4B).
To systematically investigate the biological functions of candidate genes, we extracted DEGs from upand downregulated cluster groups for further GO term and KEGG pathway analyses. GO analysis revealed that Yaolong DEGs assigned to pro les 17 and 19 were separately categorized into 18 biological processes and two molecular processes, whereas Guimi DEGs in pro le 12 were enriched with respect to 110 major functions in the Biological process, Cellular component, and Molecular function categories (Fig. S5). However, we detected no signi cant enrichment of the Guimi DEGs clustered in pro le 17. The On the basis of KEGG analysis, we identi ed nine KEGG pathways enriched with downregulated Guimi DEGs including 'metabolic pathways' (ko01100), 'biosynthesis of secondary metabolites' (ko01110), 'phenylpropanoid biosynthesis' (ko00940), and 'biosynthesis of various secondary metabolites -part 2' (ko00998) ( Figure S7). In contrast, only two pathways were enriched with upregulated Guimi DEGs, namely, 'carbon xation in photosynthetic organisms' (ko00710) and 'galactose metabolism' (ko00052).
These results thus indicate that most of the DEGs regulated during fruit development appear to be associated with the functioning of metabolic pathways.

Analysis of cell wall biogenesis during melon fruit ripening
Fifty-one DEGs in Guimi, clustered in Pro le 0 (n = 28) and Pro le 2 (n = 23), were found to be signi cantly associated with cell wall biogenesis, showing downregulated gene expression patterns, whereas in Yaolong, 48 DEGs showing downregulated trends, clustered into Pro les 0, 2, and 7, showed similar associations. Fig. 5 shows the differences in the expression trends of these DEGs between Guimi and Yaolong. We observed 33 common DEGs in cell wall biogenesis-related pro les of Guimi and Yaolong, a majority of which showed similar patterns of expression regulation in the two melon varieties at different stages of fruit development and ripening (Fig. 5A). For example, gradual reductions in the levels of GUX3 expression were observed during fruit development in both germplasms ( Fig. 5B and 5C). A similar reduction in expression was detected for ODO1, which peaked during the early stages of fruit development, after which there was a slight reduction, which became more pronounced prior to maturity. Although melon type-speci c DEGs involved in cell wall biogenesis showed an overall downregulated pattern, differences in the changes in the direction of gene expression were still different between consecutive stages of development. For example, in Guimi melon, UAM1 was downregulated from stage 1 to 2, then constantly expressed at a stable level from stage 2 to 3, prior to undergoing a decline in the mature stages. These ndings thus indicate differences in the expression patterns of melon cultivarspeci c DEGs associated with cell wall synthesis during fruit development.

Analysis of sugar metabolism during melon fruit ripening
Some of the genes involved in sugar metabolism were found to be differentially expressed during melon fruit development, including those associated with 'galactose metabolism', 'starch and sucrose metabolism', 'fructose and mannose metabolism', and 'amino sugar and nucleotide sugar metabolism' (Fig. 6A). However, in the case of Yaolong, only Pro le 0 DEGs were found to be signi cantly associated with the starch and sucrose metabolism pathway. Contrastingly, in Guimi, 'amino sugar and nucleotide sugar metabolism' was identi ed as being enriched with Pro le 2 DEGs. For a majority of the DEGs associated with the 'amino sugar and nucleotide sugar metabolism' pathway, including MUR4, UGD1, and UPTG2, the lowest levels of expression were detected in ripe fruit (Fig. 6B). Difference in gene expression were also detected during the different stages of fruit development. For example, there was a marked reduction in the expression of GAUT6 during in the middle phase of fruit development, followed by a slight increase in pre-mature fruit, before declining to the lowest levels in mature fruit. Although the change trend in UGD1 expression in young fruit was similar to that observed for GAUT6, it differed in that it was expressed at a constantly stable level in the pre-mature fruit. In Yaolong, a total of 11 Pro le 0 DEGs were found to be signi cantly associated with the starch and sucrose metabolism pathway. In contrast to trends in the expression of sugar metabolism-related genes in Guimi, nine of these DEGs were characterized by gradual downregulated expression trends during fruit development and ripening (Fig.  6C). For example, the lowest level of At4g02290 expression in Yaolong was detected at stage T3. Conversely, however, whereas INV1 was characterized by a downregulated expression trend in young fruit, high levels of expression were detected during the expansion stage, prior to a subsequent decline in expression, reaching the lowest values in ripe fruit.

Discussion
The trends we observed regarding changes in the physical and chemical characteristics of developing Guimi and Yaolong melon, including weight, size, and soluble sugar contents (Fig. 2), are presume to in uence the avor, taste, storage, and uses of Hami melon (Zhu et al. 2017). Fruit ripening is a complex genetically programmed process, which in Hami melon has not been extensively studied with respect to the underlying molecular mechanisms. To gain a better understanding of the associated gene regulation, we undertook RNA-seq of the Guimi and Yaolong varieties, on the basis of which, we identi ed 21,172 genes that were expressed among melon fruit samples during different stages development. Of these genes, the numbers of differentially expressed genes in Guimi were found to be markedly higher than those in Yaolong at 20, 30, and 40 DAP (Fig. 3), which would thus tend to indicate that Guimi fruit undergoes more pronounced changes during the course of development, which are assumed to be related to its improved characteristics compared with Yaolong. In maize, analyses at the whole-genome level revealed minimal difference between the drought-tolerant maize inbred line Han21 and drought-sensitive maize inbred line B73, with respect gene expression levels and the numbers of genes and transcription factors. However, when the authors assessed those genes that were differentially expressed between these two lines, they found that in response to drought stress, the number of DEGs in Han21 (2,832) was signi cantly higher than that in B73 (917) (Zenda et al. 2019;Thirunavukkarasu et al. 2017).
In the present study, we identi ed polysaccharide metabolism, cell wall biogenesis, and external encapsulating structure organization as the most signi cantly enriched GO terms among Guimi melon DEGs between 10 and 20 DAP (Fig. S1 and S2). Similarities in differential gene expression patterns between 10 and 20 DAF and 30 and 40 DAF during fruit development have been reported previously, indicating that gene expression is signi cantly altered during the period from 20 to 30 DAF, with implications for fruit development (Zhang et al. 2016). On the basis of our observations, we speculate that differences between Yaolong and Guimi with respect to cell wall biogenesis start to manifest at an early point in fruit development. Moreover, in Guimi melon, we detected differences with respect to a number of metabolic pathways, notably, 'metabolic pathways', 'biosynthesis of secondary metabolites' and 'phenylpropanoid biosynthesis', between 10 and 20 DAP (Figs. S3 and S4), which are probably indicative of the early modi cation of metabolic regulation during the ripening stage. In this regard, it has been established that plant polysaccharide metabolic pathways play a prominent role in mediating fruit ripening (Yamaguchi et al. 1977).
Melons are consumed for their sweet fruit (Burger et al. 2007), a trait that can be attributed to the accumulation of sucrose during the late stage of fruit development (Schaffer et al. 1987;Hubbard et al. 1989;Schaffer et al. 1996), and in this context, we found that DEGs associated with plant polysaccharide metabolic pathways showed different gene expression patterns during the fruit ripening process, thereby implying a role in the regulation of fruit ripening.
In total, we detected 4,731 and 3,198 genes that were differentially expressed in the Guimi and Yaolong varieties, respectively, during fruit development, which were clustered into 38 pro les based on STEM analysis. Three downregulated pro les (Pro les 0, 2, and 7) were obtained for both Guimi and Yaolong melon, whereas upregulated pro les were notably more prominent in Guimi melon (Fig. 4). In Guimi melon, we identi ed a number of pathways enriched with upregulated DEGs, namely, 'carbon xation in photosynthetic organisms' and 'galactose metabolism' (Fig. S8) Cell wall metabolism is one of the most important aspects of fruit ripening and is closely associated to esh texture (Sun et al. 2012), which is a signi cant quality attribute de ning commercial fruit quality (Guo et al. 2013). Plant cell walls play a signi cant role in determining the shape and elasticity of fruit, and changes in wall structure and composition are considered primary factors contributing to changes in fruit texture (Baldwin et al. 1989). Fruit texture is an important factor affecting fruit quality, storage and transportation characteristics, and disease resistance, and changes in pulp cell wall components are a primary factor contributing to changes in fruit essence (Hobson et al. 1981), whereas fruit softening is positively correlated with the degradation of cell wall components caused by cell wall hydrolases (Rasmussen et al. 1975). Studies have shown that mutations in several key genes involved in plant hormone synthesis and signal transduction can lead to a downregulated expression of cell wall-related genes. For example, the content of aldose and the degree of pectin methyl esteri cation have been found to be reduced in gibberellin gal-3 and Gai mutants (Derbyshire et al. 2007). The expression of XTH (xyloglucan endotransglucosylase hydrolase) in Arabidopsis hypocotyls has been shown to be induced by gibberellin (Liu et al. 2007), and in the present study, we established that gibberellin metabolic process (GO:0009685), cytoskeleton (GO:0005856), and cell wall organization or biogenesis (GO:0071554) were the most signi cantly enriched functions in Guimi melon DEG pro les (Fig. 4), thereby indicating that the genes involved in cell wall metabolism might play essential roles in determining fruit ripening. Moreover, we also found that the levels of GUX3 expression underwent a gradual reduction during fruit development. It is currently known that Gux, a glucuronosyltransferase-related gene, plays an important role in the formation of xylan side chains (Bromley et  Sugars are essential components of melon fruit quality and serve as prominent signals in the regulation of fruit ripening (Jia et al. 2013;Wu et al. 2014). The soluble sugars found in fruits are primarily sucrose, fructose, and glucose (Lv et al. 2000;), the accumulation of which is closely associated with fruit ripening. Sugars serve as regulatory signals for the expression of numerous genes, among which, sucrose is considered the most effective signal molecule. For example, it has been demonstrated that sucrose functions as a signal molecule that interacts with abscisic acid signaling to regulate strawberry fruit ripening (Jia et al. 2013). In the present study, we identi ed a number of genes involved in sugar metabolism showing differential expression patterns during fruit development, including those associated with 'galactose metabolism', 'starch and sucrose metabolism', 'fructose and mannose metabolism', and 'amino sugar and nucleotide sugar metabolism' (Fig. 6). In the case of several of these genes, we detected marled reductions in expression levels during the middle phase of fruit development, followed by a slight increase in the pre-mature fruit stage, and subsequent decline to the lowest recorded levels of expression in the mature fruit (Fig. 6). The contents of sugar differ according to organ and tissue, and the types of sugar also vary, as in, for example, the vascular bundles of tomato pedicels and fruit stalks, which mainly contain sucrose, whereas in other parts of the fruit, the accumulated sugars are primarily glucose and fructose. Moreover, the contents of glucose and fructose in the pericarp tissue and vascular bundles of tomato fruit differ signi cantly from those in the pectin and septum (Qi et al. 2001). Similarly, whereas the predominant sugar in grapes is glucose, peaches and apricots accumulate sucrose. These differences are generally determined by the type of sugar metabolism in fruit and the activity of related enzymes (Wang et al. 2009). It has also been established that there are considerable differences with respect sugar type and contents among the different varieties of the same plant, as has been demonstrated in a study of nine white pear varieties . In the early stages of fruit development, it has been established the sucrose transporter Pu SUT and β-glucomannases-pubglu1, Pu bglu2, and Pu bglu4-are highly expressed . Furthermore, there have been numerous studies that have examined expression of the sucrose acid invertase gene. For example, Qin et al. (2006) have characterized the cell wall acid invertase gene Mr ivr1 in Myrica rubra, the semi quantitative RT-PCR expression analysis of which revealed highest expression during the early stage of fruit development and relatively low expression in mature fruit. Thus, consistent with the ndings of the present study, it has been widely established that the expression of multiple genes associated glucose metabolism differs among different species and varieties of fruit-bearing plants, and also at different stages of fruit ripening.

Conclusion
In this study, we performed comparative transcriptomic analyses of two varieties of Hami melon (Guimi and Yaolong), with the aim of characterizing the genes and associated pathways controlling the physical and chemical properties of fruits. We successfully identi ed a number of candidate genes and pathways associated with polysaccharide metabolism, cell wall processes, biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, and carbon xation, which probably contribute to regulating the size and avor traits of these fruits. Comparing gene expression patterns at different stages of fruit development in Guimi melons will provide further insights into the regulatory mechanisms underlying the physical and chemical characteristics of melon fruit during the progression of development.