Characterization of transcriptome dynamics during watermelon fruit development: sequencing, assembly, annotation and gene expression profiles
© Guo et al; licensee BioMed Central Ltd. 2011
Received: 5 July 2011
Accepted: 21 September 2011
Published: 21 September 2011
Cultivated watermelon [Citrullus lanatus (Thunb.) Matsum. & Nakai var. lanatus] is an important agriculture crop world-wide. The fruit of watermelon undergoes distinct stages of development with dramatic changes in its size, color, sweetness, texture and aroma. In order to better understand the genetic and molecular basis of these changes and significantly expand the watermelon transcript catalog, we have selected four critical stages of watermelon fruit development and used Roche/454 next-generation sequencing technology to generate a large expressed sequence tag (EST) dataset and a comprehensive transcriptome profile for watermelon fruit flesh tissues.
We performed half Roche/454 GS-FLX run for each of the four watermelon fruit developmental stages (immature white, white-pink flesh, red flesh and over-ripe) and obtained 577,023 high quality ESTs with an average length of 302.8 bp. De novo assembly of these ESTs together with 11,786 watermelon ESTs collected from GenBank produced 75,068 unigenes with a total length of approximately 31.8 Mb. Overall 54.9% of the unigenes showed significant similarities to known sequences in GenBank non-redundant (nr) protein database and around two-thirds of them matched proteins of cucumber, the most closely-related species with a sequenced genome. The unigenes were further assigned with gene ontology (GO) terms and mapped to biochemical pathways. More than 5,000 SSRs were identified from the EST collection. Furthermore we carried out digital gene expression analysis of these ESTs and identified 3,023 genes that were differentially expressed during watermelon fruit development and ripening, which provided novel insights into watermelon fruit biology and a comprehensive resource of candidate genes for future functional analysis. We then generated profiles of several interesting metabolites that are important to fruit quality including pigmentation and sweetness. Integrative analysis of metabolite and digital gene expression profiles helped elucidating molecular mechanisms governing these important quality-related traits during watermelon fruit development.
We have generated a large collection of watermelon ESTs, which represents a significant expansion of the current transcript catalog of watermelon and a valuable resource for future studies on the genomics of watermelon and other closely-related species. Digital expression analysis of this EST collection allowed us to identify a large set of genes that were differentially expressed during watermelon fruit development and ripening, which provide a rich source of candidates for future functional analysis and represent a valuable increase in our knowledge base of watermelon fruit biology.
Watermelon [Citrullus lanatus (Thunb.) Matsum. & Nakai var. lanatus] belongs to the Cucurbitaceae family which includes several other important vegetable crops such as melon, cucumber, squash and pumpkin. It produces large edible fruits that serve as an important component in human diets throughout the world  and its farming accounts for ~7% of the world's total area devoted to vegetable production according to FAO statistics . Its production in the U.S. alone reached 4 billion pounds in 2010 with a net market value of half billion U.S. dollars. The quality of watermelon fruits consists of many factors including fruit shape and size, rind thickness and color, flesh texture and color, aroma, flavor, sugar content, carotenoid and flavonoid composition, and nutrient composition . During the development and ripening process, watermelon fruits undergo many biochemical and physiological changes including size expansion, fruit softening, and accumulation of sugars, pigments, and flavor and aromatic volatiles [4, 5]. Most of these traits are controlled by multiple QTLs and pose a significant challenge to traditional breeding [6, 7].
Currently genomics and functional genomics resources of watermelon that are publicly available are very limited. This lack of extensive genomics and functional genomics resources, combined with the narrow genetic diversity among watermelon cultivars, is one of the major limiting factors in watermelon research and breeding. However, this situation will soon be changed due to the recent advent of next-generation sequencing (NGS) technologies such as Roche/454 and Illumina/Solexa sequencing platforms. The extremely high throughput and relatively low cost of these sequencing technologies have offered unique opportunities to study genomics and functional genomics in non-model organisms. Using the NGS technologies, currently the genome of watermelon, which has an estimated size of 425 Mb , is being sequenced by the International Watermelon Genomics Initiative. The genome sequencing of cucumber, a closely-related cucurbit species, was completed , and the genome of melon, another closely-related cucurbit species, is being sequenced under the Spanish Genomics Initiative (MELONOMICS). Complementary to whole genome sequencing, which still requires huge effort and investment, large-scale transcriptome sequencing has proved to be efficient and cost-effective for gene discovery and gene function and expression analysis. Recently, several large expressed sequence tag (EST) datasets have been generated in cucurbit species. These include approximately 1.2 million, 350,000 and 500,000 ESTs generated from melon , cucumber  and Cucurbita pepo, respectively, using the Roche/454 sequencing technologies, and an additional ~127,000 ESTs generated from melon using the traditional Sanger sequencing approach . However, currently only around 12,000 watermelon ESTs are available in GenBank. Thus it's very important to expand the transcript catalog of watermelon in order to facilitate gene discovery, functional analysis, molecular breeding, and comparative genomics analysis of watermelon and its closely-related species.
Fruit is a major component of the human diet contributing a large portion of vitamins, minerals, antioxidants and fiber. Fruit development and ripening is a complex process influenced by numerous factors including light, hormones, temperature and genotype. Ripening associated events are brought about by developmentally and physiologically regulated changes in gene expression which ultimately lead to alterations in color, texture, flavor, and aroma of fruit. Fruit can be physiologically classified as climacteric or non-climacteric depending on the presence or absence of a burst in respiration at the onset of ripening . Tomato, a climacteric fruit, has served as the model system for fruit development and ripening and attracted extensive studies to understand molecular mechanisms of the development and ripening of its fleshy fruits . However, very little information is currently available at the molecular level regarding fruit development and ripening of watermelon, a non-climacteric and economically important fruit crop. To date, there is only one report describing expression profiles of a very small set of genes (832) during watermelon fruit development through microarray and qRT-PCR analysis .
In order to significantly expand the transcript catalog of watermelon and gain more insights into the molecular basis of watermelon fruit development, we performed large-scale transcriptome sequencing of watermelon fruits using the Roch/454 massively parallel pyrosequencing technology. A total of 577,023 ESTs were obtained and assembled into 75,068 unigenes which were further extensively annotated. Digital expression analysis of these ESTs identified more than 3,000 unigenes that were differentially expressed during watermelon fruit development. This information, coupled with profiles of several interesting metabolites that are important to fruit quality, provided novel insights into the biology of watermelon fruit development and ripening process.
Results and discussion
Sequencing and assembly of watermelon fruit transcriptome
Statistics of watermelon Roche/454 ESTs
No. of reads
Average length (bp)
Total bases (bp)
No. of reads in contigs
No. of reads in singletons
Distribution of number of ESTs in watermelon unigenes
No. EST members
No. ESTs in unigene
Our assembly included approximately 8,000 ESTs generated from a fruit normalized and subtracted library reported in Levi et al. . In the assembly, these ESTs were distributed in 4,616 unigenes, among which 962 (20.8%) were not captured by our 454 deep transcriptome sequencing. This indicated that although the 454 deep sequencing generated a large number of novel unigenes (more than 70,000), sequencing the transcriptome to higher depth is required to discover more rare genes.
Functional annotation, comparative genomics and pathway analysis
Statistics of watermelon unigenes with sequence matches against public protein databases
We further compared watermelon unigenes against pfam domain database . A total of 11,454 watermelon unigenes contained at least one pfam domain and 1,475 distinct pfam domains were represented by the unigenes. The most abundant pfam domains found in the collection of watermelon unigenes were PF01439 (metallothionein; 343 unigenes), PF00179 (ubiquitin-conjugating enzyme; 278), PF01200 (ribosomal protein S28e; 223), PF06522 (B12D protein; 186), and PF00025 (ADP-ribosylation factor family; 183). We also found that 465 unigenes contained transcription factor domains. The most represented pfam transcription factor domains were PF00319 (SRF-type transcription factor; 50), PF00249 (Myb-like DNA-binding domain; 42), PF00847 (AP2 domain; 37), PF00046 (homeobox domain; 27), and PF01486 (K-box region; 27).
To further identify active biochemical pathways during watermelon fruit development and ripening, we predicted biochemical pathways from our EST collection using the Pathway Tools . A watermelon metabolic pathway database has been constructed and is freely available at the Cucurbit Genomics Database . The database contained 318 metabolite pathways which were predicted from 2,968 enzyme-coding unigenes. Most major plant biochemical pathways such as calvin cycle, cellulose biosynthesis, ethylene biosynthesis, glycolysis II and IV, gluconeogenesis, and sucrose degradation, and several important secondary metabolite biosynthesis pathways including the citrulline-nitric oxide cycle, trans-lycopene biosynthesis, beta-carotene biosynthesis and flavonoid biosynthesis, were well covered by our EST collection.
Watermelon simple sequence repeats (SSRs)
During the past several years, SSR markers have been extensively used in cucurbit species such as cucumber and melon for construction of high-density genetic maps and identification of QTLs associated with economically important traits [20–22]. However, due to the limited resources of watermelon sequences, SSR markers in watermelon have been scarce. We screened watermelon unigenes for the presence of di-, tri-, tetra-, penta- and hexa-nucleotide SSR motifs and were able to predict 5,195 SSRs from 4,668 watermelon unigenes, among which 2,265, 2,709, 115, 57, and 49 were di-, tri-, tetra-, penta- and hexa-nucleotide SSR motifs, respectively. The most frequent SSR motif was AG/CT (1,616; 31.1%), followed by AAG/CTT (1,300; 25%), AT/AT (519, 10%), and AAT/ATT (465; 9%). Primer pairs were designed for SSR motifs that had sufficient flanking sequences. The complete list of SSR motifs and their corresponding primer pair information are provided in Additional file 2.
The SSRs identified in the present study provide a rich resource of valuable molecular markers in watermelon. However, polymorphism of these SSRs will of course need to be tested in specific resource populations.
Dynamic changes of gene expression and metabolite profiles during fruit development
Fruit development and ripening is a genetically programmed event that is defined by a series of biochemical and physiological changes that ultimately alter the fruit color, aroma, texture and its nutritional values. Researches toward elucidating the molecular basis of fruit development and ripening have been mainly performed on tomato, a typical climacteric fruit and model organism for fruit ripening, and numerous ripening-related genes have been isolated in tomato . However, current available information regarding the physiology, biochemistry and molecular biology of fruit development and ripening of watermelon, a non-climacteric fruit, is very limited. To gain more insights into molecular mechanisms of watermelon fruit development and ripening, we performed comprehensive digital expression analysis using ESTs generated from the four stages described above.
Validation of microarray results by qRT-PCR
Carotenoid and sugar metabolism pathway genes that are differentially expressed during watermelon fruit development
Expression (reads per million)
phytoene synthase 1
lycopene beta cyclase
Our digital expression analysis identified a total of 3,023 differentially expressed genes with at least two-fold difference in expression levels during watermelon fruit development and a false discovery rate (FDR) < 0.01 (Additional file 3). From this list of differentially expressed genes, we identified significantly enriched GO terms (FDR < 0.01), which serve as indicators of significantly altered biological processes during watermelon fruit development (Additional file 4). Among them are a number of interesting biological processes that are known to be associated with fruit development, including GO terms related to responses to different abiotic/biotic stresses, secondary metabolic process, organic acid metabolic process, and flavonoid biosynthetic/metabolic processes. We then further identified significantly altered biochemical pathways during watermelon fruit development; among which are several important pathways that are known to affect fruit quality, such as pathways of cellulose biosynthesis, suberin biosynthesis, sucrose biosynthesis and degradation, and starch biosynthesis and degradation (Additional file 5).
We further classified differentially expressed genes into different categories according to their expression patterns: 1) genes that are highly expressed in each of the four stages; 2) genes that are highly expressed in early stages; 3) genes that are highly expressed in late stages (Additional file 3). We also generated profiles of three major sugars (sucrose, fructose and glucose) that determine the sweetness of fruit and three carotenoids (lycopene, β-carotene and lutein) that play a critical role in fruit coloration and contribute significantly to fruit phyto-nutrient values. An integrative analysis of gene expression and metabolite profiles was performed to gain a better understanding of important fruit quality-related traits in watermelon.
Cell wall-related genes are highly expressed in immature white fruits
Immature watermelon fruits undergo a burst of cell division and later continue cell expansion to form large vacuolated cells that make up the majority of the flesh tissue . Cell expansion involves changes in cell wall structure and continuous accumulation (in the vacuoles) of carbohydrates, organic acids, and different compounds needed to retain the osmotic pressure and flow of water into the expanding cells . Through digital expression analysis, we identified a number of cell wall related genes that were expressed significantly higher in immature white fruits, including PRPs (proline-rich proteins), fasciclin-like arabinogalactan proteins (FLAs), xyloglucan endotransglycosylases (XETs) (Additional file 3). PRPs are a group of cell wall proteins characterized by their proline and hydroproline-rich repetitive peptides. They have been reported to be involved in cell wall formation and cell expansion in Arabidopsis, carrot and cotton [29–31]. FLAs are a subclass of arabinogalactan proteins (AGPs) that contain putative cell adhesion domains known as fasciclin domains. In eukaryotes, fasciclin domain-containing proteins are involved in cell adhesion and cell expansion [32–35]. During cell expansion and elongation, the cell wall continually undergoes temporary loosening followed by rapid reinforcement of wall structure. XETs are unique enzymes in plants that are capable of modulating the chemistry of the matrix and therefore performing both of these functions . XETs have proved to catalyze the formation of covalent linkages between xyloglucans and cellulosic substrates and between xyloglucans and (1,3; 1,4)-β-D-glucans, and influence cell wall strength, flexibility and porosity, and cell expansion [36–38].
Several other genes involved in cell division and expansion were also identified as highly expressed genes in immature white fruits, including early nodulin-like proteins (ENODs), S-adenosyl methionine decarboxylase (SAMDC) and auxin-repressed protein ARP1. ENODs have been reported to function in cell-to-cell signaling, cell differentiation, tissue development, and signal transduction pathways [39, 40]. SAMDC is the key gene involved in the biosynthesis pathway of polyamines (PAs), which are related to the cell growth and division in the early fruit development of many plant species [41–44]. Auxin is a plant growth hormone with many roles in cell division and enlargement, differentiation, and vascular bundle formation [45, 46]. The expression of auxin responsive genes was also up-regulated in the rapid expansion phase of tomato fruit and in the fast fiber cell elongation stage of cotton [47–49]. These genes might play important roles in contributing to the fast growth of the early stage watermelon fruits.
Initiation of pigment biosynthesis in white-pink flesh fruits
It is worth noting that a cinnamoyl-CoA reductase (CCR; WMU10737) was highly expressed in the fast-growing fruit (immature white and white-pink flesh). CCR catalyses the first step in the biosynthesis pathway of monolignols, the monomeric units of lignins. Lignified tissues play important roles in structural support and water/nutrient conduction . Consistent with its role in lignin biosynthesis during plant development, expression of the CCR gene was shown to be high in tissues undergoing active lignification, i.e. vascular cambium and differentiating xylem. Development of vascular system in fleshy fruits is a critical component of fruit growth and expansion, because fruits serve as the nutrient sink and the vascular system provides the framework for water, nutrients, and sugars to be transported from the vegetative parts . Thus, continued growth and expansion of watermelon fruits at early stages could rely on the establishment of the vascular system and similar connection between fruit development and vascular system formation has been reported in other fleshy fruit species like strawberry .
Rapid accumulation of sugars and secondary metabolites in mature fruits
Watermelon fruit at red flesh and over-ripe stages becomes much crispier and sweeter, and its flesh becomes red with accumulation of volatile compounds that gives watermelon its distinct aroma and flavor (Figure 1 and Additional file 1). Besides color, sweetness is another important trait of fruit quality and there has been enormous interest from growers and breeders to improve the sugar content of watermelon. Its sugar composition, as in many cucurbit fruits, is mainly determined by sucrose, fructose and glucose levels . In the present study, we measured the content of these three sugars during watermelon fruit development. At immature white and white-pink flesh stages, fructose and glucose are the predominant sugars and their levels start to decline as the fruit matured. On the contrary, the content of sucrose was low in immature white and white-pink flesh fruits and then had a sharp increase in mature fruits, in which sucrose replaced fructose and glucose as the determining factor of sugar content (Figure 5B). Previous reports indicated that the sharp increase in sucrose level in a high-sucrose-accumulating watermelon cultivar was a result of fructose and glucose reduction, and sucrose made up about 70% of the total reducing sugars in mature fruit . Our data confirmed the similar trend of sugar metabolism and, in the cultivar we examined, sucrose accounted for approximately 50% of total soluble sugar in over-ripe fruits, compared to 3-7% in fruits at immature and white-pink flesh stages (Figure 5B). We further analyzed sugar metabolism-related biochemical pathways and identified a sucrose synthase and a sucrose-phosphate synthase that were differentially expressed during watermelon fruit development (Table 5). These two genes were involved in the biosynthesis of sucrose and were highly expressed in fruits at red flesh and over-ripe stages. Correlation between sugar content and expression profiles of the two sugar metabolism related genes indicated that these two genes might play important roles in regulating sugar content in watermelon fruits.
A number of genes involved in the accumulation of second metabolites contributing to fruit flavor and aroma were found to be highly expressed in fruits at red flesh and/or over-ripe stages. They included ascorbate peroxidase (WMU23118), 9-cis-epoxycarotenoid dioxygenase (WMU05613), benzoquinone reductase (WMU23598), quinone reductase (WMU39499), flavonoid 3',5'-hdyroxylase (WMU39890), and squalene synthase (WMU27722 and WMU15567). Identification of these genes provided a rich source for further dissection of molecular mechanisms that governing fruit flavor and aroma, the important traits of fruit quality that currently are not well understood.
Watermelon is an important fruit crop and becomes a useful model for the research of non-climacteric fruits. However, genetic and genomic resources for watermelon are very scarce, which are among the major limiting factors in watermelon research and breeding. In the present study, using the Roche/454 massively parallel pyrosequencing technology, we have generated more than half million watermelon ESTs from four stages of watermelon fruit development. These ESTs were de novo assembled and extensively annotated. They represent the significant expansion of the watermelon transcript catalog and provide a comprehensive material basis for future functional and expression analysis of genes of interest. The availability of these ESTs will also facilitate the annotation of the watermelon genome that is currently being sequenced. SSRs were also identified from these ESTs, which provides a valuable resource for the development of molecular markers that can be further used to facilitate the watermelon breeding program and cloning genes of interest. Integrative analysis of digital gene expression and metabolite profiles provided novel insights into molecular mechanisms of watermelon fruit development and regulatory mechanisms of several import traits of watermelon fruit quality.
Seeds of watermelon inbred line (Citrullus lanatus (Thunb.) Matsum. & Nakai var. lanatus cv 97103) were germinated and grown in greenhouse with nutrition pots containing a soil mixture (peat:sand:pumice, 1:1:1, v/v/v). Flowers were hand-pollinated and tagged. Center fruit flesh samples were collected at stages of 10, 18, 28, and 34 DAP, respectively. Tissues were immediately frozen in liquid nitrogen and stored at -80°C till use.
cDNA preparation and sequencing
Total RNA was extracted from watermelon fruit flesh samples using the TRIzol Reagent (Invitrogen, USA). mRNA was purified from the total RNA using the Oligotex mRNA Midi Kit (QIAGEN, Germany). Double-strand cDNA was then synthesized using the SMART cDNA Library Construction kit (Clontech, USA) following the manufacturer's instructions. PCR products of cDNA were further purified using the QIAquick PCR Purification Kit (QIAGEN, Germany) and checked for quality using the Agilent 2100 Bioanalyzer. Approximately 10 ug cDNA from each of the four fruit flesh samples was used for sequencing on a Roche/454 GS-FLX Titanium instrument. A half-plate sequencing run was performed for each sample at the Cornell University Life Sciences Core Laboratories Center following manufacturer's instructions.
cDNA sequence processing, assembly, annotation and comparative genomics analysis
The raw 454 sequence files in SFF format were base called using the Pyrobayes base caller . In addition, around 12,000 watermelon ESTs were collected from GenBank. All these sequences were then processed to remove low quality regions and adaptor sequences using programs LUCY  and SeqClean . The resulting sequences were then screened against the NCBI UniVec database, E. coli genome sequences, and watermelon ribosomal RNA sequences, to remove possible contaminations of these sequences. Sequences shorter than 100 bp were discarded. The processed 454 and GenBank sequences were assembled into unigenes using the iAssembler program .
The resulting watermelon unigene sequences were compared against GenBank non-redundant protein (nr) and UniProt (TrEMBL and SwissProt) databases using the BLAST program with a cutoff E-value of 1e-5. The unigene sequences were also translated into proteins using ESTScan  and the translated protein sequences were then compared to pfam domain databases. Gene ontology (GO; ) terms were assigned to each unigene based on the GO terms annotated to its corresponding homologues in the UniProt databases and domains in pfam database . GO annotations of watermelon unigenes were then mapped to a list of plant-specific GO slim ontology . Annotations of unigenes were then formatted into the PathoLogic format and used to predict watermelon biochemical pathways using the Pathway Tools .
For comparative genomics analysis, watermelon unigenes were compared to protein databases of cucumber (version 2), a closely-related cucurbit species, Arabidopsis (TAIR version 10), a model system of dicot plants, and rice (RGAP 6.1), a model system of monocot plants, using the BLAST program with an E-value cutoff of 1e-5. Cucumber, Arabidopsis and rice protein databases were obtained from the following links, respectively:
Watermelon SSR identification
SSRs in watermelon unigene sequences were identified using the MISA program . The minimum repeat number was six for dinucleotide and five for tri-, tetra-, penta- and hexa-nucleotide. Primer pairs flanking each SSR loci were designed using the Primer3 program .
Identification of differentially expressed genes
Following cDNA sequence assembly, EST copy numbers of each unigene in each of the four stages of fruits were derived and used as an approximate estimation of gene expression levels in the corresponding tissues following normalization to RPM (reads per million sequenced reads). Significance of differential gene expression during watermelon fruit development was determined using the R statistic  and the resulting raw p values were corrected for multiple tests using the False Discovery Rate (FDR; ). Genes with expression changes no less than two during fruit development and FDRs less than 0.01 were identified as differentially expressed genes. GO terms enriched in the set of differentially expressed genes were identified using GO::TermFinder , requiring p values adjusted for multiple testing to be less than 0.01. Significantly altered pathways during fruit development were identified using the Plant MetGenMAP system .
Determination of the content of sugars, soluble solids and carotenoids, and flesh firmness
Sugars including sucrose, glucose and fructose were determined according to protocols described in Bethke et al.  and Yativ et al. , with minor modifications. Two hundred milligrams of frozen watermelon flesh samples were ground to a fine powder prior to being extracted for 1 h in 10 ml of 50% ethanol at 80°C, then centrifuged at 3000 g for 10 min. This step was repeated one more time and the supernatants were then dissolved to a volumetric flask (25 ml) as extracts. Two milliliter extracts were centrifuged at 3000 g for 10 min. The supernatants (1 ml) were filtered through a 0.45-mm HPLC nylon filter (Membrana, Germany). Sugars were separated in an analytical HPLC system (Pump System LC-10ATVP, Shimadzu, Japan) fitted with a Shodex Asahipak NH2P-50 4E column (4.6 × 250 mm, Shodex, Japan) using a refractive-index detector (RID-10A, Shimadzu, Japan). Soluble solid content (SSC) was measured using juices extracted from the center flesh of watermelon fruits with a digital refractometer (Digital Refractometer PR-1, Atago, Japan). The content of three carotenoids (lycopene, β-carotene and lutein) was determined according to the protocol described in Fraser et al. . Individual carotenoids were separated by HPLC on a Waters Nova-pak C18 column. Flesh firmness was measured on the center of watermelon flesh, using a pressure tester (FT327, Breuzzi, Italy). All the above experiments were performed using two biological replicates.
Two sugar metabolism genes (WMU23179, sucrose-phosphate synthase and WMU23817, sucrose synthase) were selected for qRT-PCR analysis. qRT-PCR reactions were run in a Roche 480 Real-Time PCR system. Each 20 μL reaction consisted of 5 μL cDNA (~20 ng/uL), 5 μL of primer mix (2 μM of each forward and reverse primer), 10 μL of lightcycler 480 SYBR Green I Master Mix. Cycling conditions were 95°C for 3 min, followed by 45 cycles of: 95°C for 10 sec, 58°C for 20 sec, and 72°C for 30 sec. Melting curves were performed at the end of each reaction run to detect primer dimers. The 18S rRNA gene was used as the internal control. Quantification was achieved by normalizing the number of target gene copies to a reference 18S rRNA gene by using the comparative Ct method . The ΔCt was calculated by subtracting the average Ct value of each tissue type from the average Ct values of 18S rRNA. The ΔΔCt was calculated by subtracting the ΔCt of each of the four fruit stages from the ΔCt of the fruit tissue. The formula 2 ^-(ΔΔCt) was used to calculate a relative fold change between the four fruit developing stages. Primer pairs used for WMU23179, WMU23817, and 18S rRNA were F:5'-TAACCTAGTGGTTGTGGCTGGAGA-3' and R:5'- TGGCCATTCAGGTTGTAGGTCT-3', F:5'-TCGCATCAAGAGCTCAAGCACTCA-3' and R:5'- GCACTTGCATCACCTGGTTTCCAT-3', and F:5'-AGCCTGAGAAACGGCTACCACATC-3' and R:5'-ACCAGACTCGAAGAGCCCGGTAT, respectively.
This research was supported by grants from the Ministry of Science and Technology of the People's Republic of China (2009BADB8B02, 2010DFB33740, 2010AA10A107, 30972015 and 31171980), Ministry of Agriculture of the People's Republic of China (CARS-26, 2008-Z42) and Beijing Municipal Science & Technology Commission, China (D08070500690803, D111100001311002, KJCX201101010, Z09090501040902 and 5100001) to Y.X. and H.Z., United States-Israel Binational Agricultural Research and Development Fund (IS-4223-09C) to Z.F.
- Erickson DL, Smith BD, Clarke AC, Sandweiss DH, Tuross N: An Asian origin for a 10,000-year-old domesticated plant in the Americas. Proc Natl Acad Sci USA. 2005, 102: 18315-18320. 10.1073/pnas.0509279102.PubMed CentralPubMedView Article
- FAO Statistics. [http://faostat.fao.org]
- Wechter WP, Levi A, Harris KR, Davis AR, Fei Z, Katzir N, Giovannoni JJ, Salman-Minkov A, Hernandez A, Thimmapuram J, et al: Gene expression in developing watermelon fruit. BMC Genomics. 2008, 9: 275-10.1186/1471-2164-9-275.PubMed CentralPubMedView Article
- Tadmor Y, King S, Levi A, Davis A, Meir A, Wasserman B, Hirschberg J, Lewinsohn E: Comparative fruit colouration in watermelon and tomato. Food Res Int. 2005, 38: 837-841. 10.1016/j.foodres.2004.07.011.View Article
- Yativ M, Harary I, Wolf S: Sucrose accumulation in watermelon fruits: genetic variation and biochemical analysis. J Plant Physiol. 2010, 167: 589-596. 10.1016/j.jplph.2009.11.009.PubMedView Article
- Hashizume T, Shimamoto I, Hirai M: Construction of a linkage map and QTL analysis of horticultural traits for watermelon [Citrullus lanatus (THUNB.) MATSUM & NAKAI] using RAPD, RFLP and ISSR markers. Theor Appl Genet. 2003, 106: 779-785.PubMed
- Gusmini G, Wehner TC: Qualitative inheritance of rind pattern and flesh color in watermelon. J Hered. 2006, 97: 177-185. 10.1093/jhered/esj023.PubMedView Article
- Arumuganathan K, Earle E: Nuclear DNA content of some important plant species. Plant Mol Biol Rep. 1991, 9: 208-218. 10.1007/BF02672069.View Article
- Huang S, Li R, Zhang Z, Li L, Gu X, Fan W, Lucas WJ, Wang X, Xie B, Ni P, et al: The genome of the cucumber, Cucumis sativus L. Nat Genet. 2009, 41: 1275-1281. 10.1038/ng.475.PubMedView Article
- Portnoy V, Diber A, Pollock S, Karchi H, Lev S, Tzuri G, Harel-Beja R, Forer R, Portnoy VH, Lewinsohn E, et al: Use of non-normalized, non-amplified cDNA for 454-based RNA sequencing of fleshy melon fruit. Plant Gen. 2011, 4: 36-46. 10.3835/plantgenome2010.11.0026.View Article
- Guo S, Zheng Y, Joung JG, Liu S, Zhang Z, Crasta OR, Sobral BW, Xu Y, Huang S, Fei Z: Transcriptome sequencing and comparative analysis of cucumber flowers with different sex types. BMC Genomics. 2010, 11: 384-10.1186/1471-2164-11-384.PubMed CentralPubMedView Article
- Blanca J, Canizares J, Roig C, Ziarsolo P, Nuez F, Pico B: Transcriptome characterization and high throughput SSRs and SNPs discovery in Cucurbita pepo (Cucurbitaceae). BMC Genomics. 2011, 12: 104-10.1186/1471-2164-12-104.PubMed CentralPubMedView Article
- Clepet C, Joobeur T, Zheng Y, Jublot D, Huang M, Truniger V, Boualem A, Hernandez-Gonzalez ME, Dolcet-Sanjuan R, Portnoy V, et al: Analysis of expressed sequence tags generated from full-length enriched cDNA libraries of melon. BMC Genomics. 2011, 12: 252-10.1186/1471-2164-12-252.PubMed CentralPubMedView Article
- Lelievre JM, Latche A, Jones B, Bouzayen M, Pech JC: Ethylene and fruit ripening. Physiol Plant. 1997, 101: 727-739. 10.1111/j.1399-3054.1997.tb01057.x.View Article
- Giovannoni JJ: Genetic regulation of fruit development and ripening. Plant Cell. 2004, 16: S170-S180. 10.1105/tpc.019158.PubMed CentralPubMedView Article
- Cucurbit Genomics Database. [http://www.icugi.org]
- Levi A, Davis A, Hernandez A, Wechter P, Thimmapuram J, Trebitsh T, Tadmor Y, Katzir N, Portnoy V, King S: Genes expressed during the development and ripening of watermelon fruit. Plant Cell Rep. 2006, 25: 1233-1245. 10.1007/s00299-006-0163-0.PubMedView Article
- Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K, et al: The Pfam protein families database. Nucleic Acids Res. 2010, 38: D211-D222. 10.1093/nar/gkp985.PubMed CentralPubMedView Article
- Karp PD, Paley S, Romero P: The Pathway Tools software. Bioinformatics. 2002, 18: S225-S232. 10.1093/bioinformatics/18.suppl_1.S225.PubMedView Article
- Ren Y, Zhang Z, Liu J, Staub JE, Han Y, Cheng Z, Li X, Lu J, Miao H, Kang H, et al: An integrated genetic and cytogenetic map of the cucumber genome. PLoS One. 2009, 4: e5795-10.1371/journal.pone.0005795.PubMed CentralPubMedView Article
- Harel-Beja R, Tzuri G, Portnoy V, Lotan-Pompan M, Lev S, Cohen S, Dai N, Yeselson L, Meir A, Libhaber SE, 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. 10.1007/s00122-010-1327-4.PubMedView Article
- Díaz A, Fergani M, Formisano G, Ziarsolo P, Blanca J, Fei Z, Staub JE, Zalapa JE, Cuevas HE, 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-10.1186/1471-2229-11-111.PubMed CentralPubMedView Article
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009, 10: 57-63. 10.1038/nrg2484.PubMed CentralPubMedView Article
- Fei Z, Tang X, Alba RM, White JA, Ronning CM, Martin GB, Tanksley SD, Giovannoni JJ: Comprehensive EST analysis of tomato and comparative genomics of fruit ripening. Plant J. 2004, 40: 47-59. 10.1111/j.1365-313X.2004.02188.x.PubMedView Article
- Wang SM: Understanding SAGE data. Trends Genet. 2007, 23: 42-50. 10.1016/j.tig.2006.11.001.PubMedView Article
- Kang B, Zhao W, Hou Y, Tian P: Expression of carotenogenic genes during development and ripening of watermelon fruit. Scientia Horticulturae. 2010, 124: 368-375. 10.1016/j.scienta.2010.01.027.View Article
- Sedgley M, Newbury HJ, Possingham JV: Early fruit development in the watermelon: anatomical comparison of pollinated, auxin-induced parthenocarpic and unpollinated fruits. Ann Bot-London. 1977, 41: 1345-1355.
- Carrari F, Fernie AR: Metabolic regulation underlying tomato fruit development. J Exp Bot. 2006, 57: 1883-1897. 10.1093/jxb/erj020.PubMedView Article
- John ME, Keller G: Characterization of mRNA for a proline-rich protein of cotton fiber. Plant Physiol. 1995, 108: 669-676. 10.1104/pp.108.2.669.PubMed CentralPubMedView Article
- Fowler TJ, Bernhardt C, Tierney ML: Characterization and expression of four proline-rich cell wall protein genes in Arabidopsis encoding two distinct subsets of multiple domain proteins. Plant Physiol. 1999, 121: 1081-1092. 10.1104/pp.121.4.1081.PubMed CentralPubMedView Article
- Holk A, Klumpp L, Scherer GF: A cell wall protein down-regulated by auxin suppresses cell expansion in Daucus carota (L.). Plant Mol Biol. 2002, 50: 295-305. 10.1023/A:1016052613196.PubMedView Article
- Johnson KL, Jones BJ, Bacic A, Schultz CJ: The fasciclin-like arabinogalactan proteins of Arabidopsis. A multigene family of putative cell adhesion molecules. Plant Physiol. 2003, 133: 1911-1925. 10.1104/pp.103.031237.PubMed CentralPubMedView Article
- Shi H, Kim Y, Guo Y, Stevenson B, Zhu JK: The Arabidopsis SOS5 locus encodes a putative cell surface adhesion protein and is required for normal cell expansion. Plant Cell. 2003, 15: 19-32. 10.1105/tpc.007872.PubMed CentralPubMedView Article
- Lee KJ, Sakata Y, Mau SL, Pettolino F, Bacic A, Quatrano RS, Knight CD, Knox JP: Arabinogalactan proteins are required for apical cell extension in the moss Physcomitrella patens. Plant Cell. 2005, 17: 3051-3065. 10.1105/tpc.105.034413.PubMed CentralPubMedView Article
- Ma H, Zhao J: Genome-wide identification, classification, and expression analysis of the arabinogalactan protein gene family in rice (Oryza sativa L.). J Exp Bot. 2010, 61: 2647-2668. 10.1093/jxb/erq104.PubMed CentralPubMedView Article
- Eckardt NA: Inside the matrix: crystal structure of a xyloglucan endotransglycosylase. Plant Cell. 2004, 16: 792-793. 10.1105/tpc.160411.PubMed CentralView Article
- Hrmova M, Farkas V, Lahnstein J, Fincher GB: A Barley xyloglucan xyloglucosyl transferase covalently links xyloglucan, cellulosic substrates, and (1,3;1,4)-beta-D-glucans. J Biol Chem. 2007, 282: 12951-12962. 10.1074/jbc.M611487200.PubMedView Article
- Nishikubo N, Takahashi J, Roos AA, Derba-Maceluch M, Piens K, Brumer H, Teeri TT, Stalbrand H, Mellerowicz EJ: Xyloglucan endo-transglycosylase-mediated xyloglucan rearrangements in developing wood of hybrid aspen. Plant Physiol. 2011, 155: 399-413. 10.1104/pp.110.166934.PubMed CentralPubMedView Article
- Fruhling M, Schroder G, Hohnjec N, Puhler A, Perlick AM, Kuster H: The promoter of the Vicia faba L. gene VfEnod12 encoding an early nodulin is active in cortical cells and nodule primordia of transgenic hairy roots of Vicia hirsuta as well as in the prefixing zone II of mature transgenic V. hirsuta root nodules. Plant Sci. 2000, 160: 67-75. 10.1016/S0168-9452(00)00362-9.PubMedView Article
- Khan JA, Wang Q, Sjölund RD, Schulz A, Thompson GA: An early nodulin-like protein accumulates in the sieve element plasma membrane of Arabidopsis. Plant Physiol. 2007, 143: 1576-1589. 10.1104/pp.106.092296.PubMed CentralPubMedView Article
- Theiss C, Bohley P, Voigt J: Regulation by polyamines of ornithine decarboxylase activity and cell division in the unicellular green alga Chlamydomonas reinhardtii. Plant Physiol. 2002, 128: 1470-1479. 10.1104/pp.010896.PubMed CentralPubMedView Article
- Liu J, Nada K, Pang X, Honda C, Kitashiba H, Moriguchi T: Role of polyamines in peach fruit development and storage. Tree Physiol. 2006, 26: 791-798.PubMedView Article
- Liu JH, Honda C, Moriguchi T: Involvement of polyamine in floral and fruit development. Japan Agricultural Research Quarterly. 2006, 40: 51-58.View Article
- Trenor M, Perez-Amador MA, Carbonell J, Blazquez MA: Expression of polyamine biosynthesis genes during parthenocarpic fruit development in Citrus clementina. Planta. 2010, 231: 1401-1411. 10.1007/s00425-010-1141-x.PubMedView Article
- Fukuda H: Tracheary element differentiation. Plant Cell. 1997, 9: 1147-1156. 10.1105/tpc.9.7.1147.PubMed CentralPubMedView Article
- Buchanan BB, Gruissem W, Jones RL: Biochemistry and molecular biology of plants. 2000, Rockville, MD: American Society of Plant Physiologists
- Lemaire-Chamley M, Petit J, Garcia V, Just D, Baldet P, Germain V, Fagard M, Mouassite M, Cheniclet C, Rothan C: Changes in transcriptional profiles are associated with early fruit tissue specialization in tomato. Plant Physiol. 2005, 139: 750-769. 10.1104/pp.105.063719.PubMed CentralPubMedView Article
- Gou JY, Wang LJ, Chen SP, Hu WL, Chen XY: Gene expression and metabolite profiles of cotton fiber during cell elongation and secondary cell wall synthesis. Cell Res. 2007, 17: 422-434.PubMed
- Mounet F, Moing A, Garcia V, Petit J, Maucourt M, Deborde C, Bernillon S, Le Gall G, Colquhoun I, Defernez M, et al: Gene and metabolite regulatory network analysis of early developing fruit tissues highlights new candidate genes for the control of tomato fruit composition and development. Plant Physiol. 2009, 149: 1505-1528. 10.1104/pp.108.133967.PubMed CentralPubMedView Article
- Di Mascio P, Kaiser S, Sies H: Lycopene as the most efficient biological carotenoid singlet oxygen quencher. Arch Biochem Biophys. 1989, 274: 532-538. 10.1016/0003-9861(89)90467-0.PubMedView Article
- Fraser PD, Truesdale MR, Bird CR, Schuch W, Bramley PM: Carotenoid biosynthesis during tomato fruit development (evidence for tissue-specific gene expression). Plant Physiol. 1994, 105: 405-413.PubMed CentralPubMed
- Pecker I, Gabbay R, Cunningham FJ, Hirschberg J: Cloning and characterization of the cDNA for lycopene beta-cyclase from tomato reveals decrease in its expression during fruit ripening. Plant Mol Biol. 1996, 30: 807-819. 10.1007/BF00019013.PubMedView Article
- Lacombe E, Van Doorsselaere J, Boerjan W, Boudet AM, Grima-Pettenati J: Characterization of cis-elements required for vascular expression of the cinnamoyl CoA reductase gene and for protein-DNA complex formation. Plant J. 2000, 23: 663-676. 10.1046/j.1365-313x.2000.00838.x.PubMedView Article
- Aharoni A, Keizer LC, Van Den Broeck HC, Blanco-Portales R, Munoz-Blanco J, Bois G, Smit P, De Vos RC, O'Connell AP: Novel insight into vascular, stress, and auxin-dependent and -independent gene expression programs in strawberry, a non-climacteric fruit. Plant Physiol. 2002, 129: 1019-1031. 10.1104/pp.003558.PubMed CentralPubMedView Article
- Quinlan AR, Stewart DA, Stromberg MP, Marth GT: Pyrobayes: an improved base caller for SNP discovery in pyrosequences. Nat Methods. 2008, 5: 179-181. 10.1038/nmeth.1172.PubMedView Article
- Chou HH, Holmes MH: DNA sequence quality trimming and vector removal. Bioinformatics. 2001, 17: 1093-1104. 10.1093/bioinformatics/17.12.1093.PubMedView Article
- SeqClean program. [http://compbio.dfci.harvard.edu/tgi/software]
- iAssembler program. [http://bioinfo.bti.cornell.edu/tool/iAssembler]
- Iseli C, Jongeneel CV, Bucher P: ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc Int Conf Intell Syst Mol Biol. 1999, 138-148.
- Gene Ontology. [http://www.geneontology.org]
- Camon E, Magrane M, Barrell D, Lee V, Dimmer E, Maslen J, Binns D, Harte N, Lopez R, Apweiler R: The Gene Ontology Annotation (GOA) Database: sharing knowledge in Uniprot with Gene Ontology. Nucleic Acids Res. 2004, 32: D262-D266. 10.1093/nar/gkh021.PubMed CentralPubMedView Article
- Plant specific GO slims. [http://www.geneontology.org/GO.slims.shtml]
- MISA program. [http://pgrc.ipk-gatersleben.de/misa]
- Primer3 program. [http://frodo.wi.mit.edu]
- Stekel DJ, Git Y, Falciani F: The comparison of gene expression from multiple cDNA libraries. Genome Res. 2000, 10: 2055-2061. 10.1101/gr.GR-1325RR.PubMed CentralPubMedView Article
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995, 57: 289-300.
- Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, Sherlock G: GO:TermFinder-open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics. 2004, 20: 3710-3715. 10.1093/bioinformatics/bth456.PubMed CentralPubMedView Article
- Joung JG, Corbett AM, Fellman SM, Tieman DM, Klee HJ, Giovannoni JJ, Fei Z: Plant MetGenMAP: an integrative analysis system for plant systems biology. Plant Physiol. 2009, 151: 1758-1768. 10.1104/pp.109.145169.PubMed CentralPubMedView Article
- Bethke P, Sabba R, Bussan A: Tuber water and pressure potentials decrease and sucrose contents increase in response to moderate drought and heat stress. Am J Potato Res. 2009, 86: 519-10.1007/s12230-009-9109-8.View Article
- Bièche I, Laurendeau I, Tozlu S, Olivi M, Vidaud D, Lidereau R, Vidaud M: Quantitation of MYC gene expression in sporadic breast tumors with a real-time reverse transcription-PCR assay. Cancer Res. 1999, 59: 2759-2765.PubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.