Search

The Cucurbitaceae (cucurbit) family consists of about 1,000 species in 95 genera, including many economically important and popular fruit and vegetable crops. During the past several years, reference genomes have been generated for >20 cucurbit species, and variome and transcriptome profiling data have been rapidly accumulated for cucurbits. To efficiently mine, analyze and disseminate these large-scale datasets, we have developed an updated version of Cucurbit Genomics Database. The updated database, CuGenDBv2 (http://cucurbitgenomics.org/v2), currently hosts 34 reference genomes from 27 cucurbit species/subspecies belonging to 10 different genera. Protein-coding genes from these genomes have been comprehensively annotated by comparing their protein sequences to various public protein and domain databases. A novel ‘Genotype’ module has been implemented to facilitate mining and analysis of the functionally annotated variome data including SNPs and small indels from large-scale genome sequencing projects. An updated ‘Expression’ module has been developed to provide a comprehensive gene expression atlas for cucurbits. Furthermore, synteny blocks between any two and within each of the 34 genomes, representing a total of 595 pair-wise genome comparisons, have been identified and can be explored and visualized in the database.

The mutation of ClZISO identified in EMS-induced watermelon leads to photosensitive flesh in watermelon. Watermelon (Citrullus lanatus) has a colorful flesh that attracts consumers and benefits human health. We developed an ethyl-methanesulfonate mutation library in red-fleshed line '302' to create new flesh color lines and found a yellow-fleshed mutant which accumulated ζ-carotene. The initial yellow color of this mutant can be photobleached within 10 min under intense sunlight. A long-term light-emitting diode (LED) light treatment turned flesh color from yellow to pink. We identified this unique variation as photosensitive flesh mutant ('psf'). Using bulked segregant analysis, we fine-mapped an EMS-induced G-A transversion in 'psf' which leads to a premature stop codon in 15-cis-ζ-carotene isomerase (ClZISO) gene. We detected that wild-type ClZISO is expressed in chromoplasts to catalyze the conversion of 9,15,9'-tri-cis-ζ-carotene to 9,9'-di-cis-ζ-carotene. The truncated ClZISO

Cucurbitaceae is one of the most genetically diverse plant families in the world. Many of them are important vegetables or medicinal plants and are widely distributed worldwide. The rapid development of sequencing technologies and bioinformatic algorithms has enabled the generation of genome sequences of numerous important Cucurbitaceae species. This has greatly facilitated research on gene identification, genome evolution, genetic variation, and molecular breeding of cucurbit crops. So far, genome sequences of 18 different cucurbit species belonging to tribes Benincaseae, Cucurbiteae, Sicyoeae, Momordiceae, and Siraitieae have been deciphered. This review summarizes the genome sequence information, evolutionary relationships, and functional genes associated with important agronomic traits (e.g. fruit quality). The progress of molecular breeding in cucurbit crops and prospects for future applications of Cucurbitaceae genome information are also discussed.

NAC (NAM, ATAF1/2, and CUC2) transcription factors play important roles in fruit ripening and quality. The watermelon genome encodes 80 NAC genes, and 21 of these NAC genes are highly expressed in both the flesh and vascular tissues. Among these genes, ClNAC68 expression was significantly higher in flesh than in rind. However, the intrinsic regulatory mechanism of ClNAC68 in fruit ripening and quality is still unknown. In this study, we found that ClNAC68 is a transcriptional repressor and that the repression domain is located in the C-terminus. Knockout of ClNAC68 by the CRISPR-Cas9 system decreased the soluble solid content and sucrose accumulation in mutant flesh. Development was delayed, germination was inhibited, and the IAA content was significantly decreased in mutant seeds. Transcriptome analysis showed that the invertase gene ClINV was the only gene involved in sucrose metabolism that was upregulated in mutant flesh, and expression of the indole-3-acetic acid-amido synthetase gene ClGH3.6 in the IAA signaling pathway was also induced in mutant seeds. EMSA and dual-luciferase assays showed that ClNAC68 directly bound to the promoters of ClINV and ClGH3.6 to repress their expression. These results indicated that ClNAC68 positively regulated sugar and IAA accumulation by repressing ClINV and ClGH3.6. Our findings provide new insights into the regulatory mechanisms by which NAC transcription factors affect fruit quality and seed development.

Fruit characteristics of sweet watermelon are largely the result of human selection. Here we report an improved watermelon reference genome and whole-genome resequencing of 414 accessions representing all extant species in theCitrullusgenus. Population genomic analyses reveal the evolutionary history ofCitrullus, suggesting independent evolutions inCitrullus amarusand the lineage containingCitrullus lanatusandCitrullus mucosospermus. Our findings indicate that different loci affecting watermelon fruit size have been under selection during speciation, domestication and improvement. A non-bitter allele, arising in the progenitor of sweet watermelon, is largely fixed inC. lanatus. Selection for flesh sweetness started in the progenitor ofC. lanatusand continues through modern breeding on loci controlling raffinose catabolism and sugar transport. Fruit flesh coloration and sugar accumulation might have co-evolved through shared genetic components including a sugar transporter gene. This study provides valuable genomic resources and sheds light on watermelon speciation and breeding history.

To understand sex determination in watermelon (Citrullus lanatus), a spontaneous gynoecious watermelon mutant, XHBGM, was selected from the monoecious wild type XHB. Using map-based cloning, resequencing and fluorescence in situ hybridization analysis, a unique chromosome translocation between chromosome 2 and chromosome 3 was found in XHBGM. Based on the breakpoint location in chromosome 2, a putative C2H2 zinc finger transcription factor gene, ClWIP1 (gene ID Cla008537), an orthologue of the melon gynoecy gene CmWIP1, was disrupted. Using clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated system 9 to edit ClWIP1, we obtained gynoecious watermelon lines. Functional studies showed that ClWIP1 is expressed specifically in carpel primordia and is related to the abortion of carpel primordia in early floral development. To identify the cellular and metabolic processes associated with ClWIP1, we compared the shoot apex transcriptomes of two gynoecious mutants and their corresponding wild types. Transcriptome analysis showed that differentially expressed genes related to the ethylene and cytokinin pathways were upregulated in the gynoecious mutants. This study explores the molecular mechanism of sex determination in watermelon and provides a theoretical and technical basis for breeding elite gynoecious watermelon lines.

[EN] The Cucurbitaceae family (cucurbit) includes several economically important crops, such as melon, cucumber, watermelon, pumpkin, squash and gourds. During the past several years, genomic and genetic data have been rapidly accumulated for cucurbits. To store, mine, analyze, integrate and disseminate these large-scale datasets and to provide a central portal for the cucurbit research and breeding community, we have developed the Cucurbit Genomics Database (CuGenDB; http://cucurbitgenomics.org) using the Tripal toolkit. The database currently contains all available genome and expressed sequence tag (EST) sequences, genetic maps, and transcriptome profiles for cucurbit species, as well as sequence annotations, biochemical pathways and comparative genomic analysis results such as synteny blocks and homologous gene pairs between different cucurbit species. A set of analysis and visualization tools and user-friendly query interfaces have been implemented in the database to facilitate the usage of these large-scale data by the community. In particular, two new tools have been developed in the database, a `SyntenyViewer¿ to view genome synteny between different cucurbit species and an `RNA-Seq¿ module to analyze and visualize gene expression profiles. Both tools have been packed as Tripal extension modules that can be adopted in other genomics databases developed using the Tripal system., USDA National Institute of Food and Agriculture Specialty Crop Research Initiative [2015-51181-24285]; US-Israel Binational Agricultural Research and Development Fund [IS-3333-02, IS-3877-06CR and IS-4223-09C]; USDA Agricultural Research Service, and by SNC Laboratoire ASL, de Ruiter Seeds B.V., Enza Zaden B.V., Gautier Semences S.A., Nunhems B.V., Rijk Zwaan B.V., Sakata Seed Inc, Semillas Fito S.A., Seminis Vegetable Seeds Inc, Syngenta Seeds B.V., Takii and Company Ltd, Vilmorin and Cie S.A. and Zeraim Gedera Ltd, all of them as part of the support to the International Cucurbit Genomics Initiative (ICuGI). Funding for open access charge: USDA National Institute of Food and Agriculture.

Abscisic acid (ABA) is a critical regulator of seed development and germination. β-glucosidases (BGs) have been suggested to be contributors to increased ABA content because they catalyze the hydrolysis of ABA-glucose ester to release free ABA. However, whether BGs are involved in seed development is unclear. In this study, a candidate gene, ClBG1, in watermelon was selected for targeted mutagenesis via the CRISPR/Cas9 system. Seed size and weight were significantly reduced in the Clbg1-mutant watermelon lines, which was mainly attributed to decreased cell number resulting from decreased ABA levels. A transcriptome analysis showed that the expression of 1015 and 1429 unique genes was changed 10 and 18 days after pollination (DAP), respectively. Cytoskeleton- and cell cycle-related genes were enriched in the differentially expressed genes of wild type and Clbg1-mutant lines during seed development. Moreover, the expression of genes in the major signaling pathways of seed size control was also changed. In addition, seed germination was promoted in the Clbg1-mutant lines due to decreased ABA content. These results indicate that ClBG1 may be critical for watermelon seed size regulation and germination mainly through the modulation of ABA content and thereby the transcriptional regulation of cytoskeleton-, cell cycle- and signaling-related genes. Our results lay a foundation for dissecting the molecular mechanisms of controlling watermelon seed size, a key agricultural trait of significant economic importance.

Red-fleshed watermelons (Citrullus lanatus) that accumulate lycopene in their flesh cells have been selected and domesticated from their pale-fleshed ancestors. However, the molecular basis of this trait remains poorly understood. Using map-based cloning and transgenic analysis, we identified a lycopene β-cyclase (ClLCYB) gene that controls the flesh color of watermelon. Down-regulation of ClLCYB caused the flesh color to change from pale yellow to red, and ClLCYB overexpression in the red-fleshed line caused the flesh color to change to orange. Analysis of ClLCYB single-nucleotide polymorphisms using 211 watermelon accessions with different flesh colors revealed that two missense mutations between three haplotypes (ClLCYB(red), ClLCYB(white), and ClLCYB(yellow)) were selected and largely fixed in domesticated watermelon. Proteins derived from these three ClLCYB haplotypes were localized in plastids to catalyze the conversion of lycopene to β-carotene and showed similar catalytic abilities. We revealed that ClLCYB protein abundance, instead of ClLCYB transcript level, was negatively correlated with lycopene accumulation. Different amounts of ClLCYB protein degradation among the ClLCYB haplotypes were found in ClLCYB transgenic Arabidopsis (Arabidopsis thaliana) lines. After treatment with the proteasome inhibitor MG132, the concentration of ClLCYB(red) increased noticeably compared with other ClLCYB proteins. These results indicate that natural missense mutations within ClLCYB influence ClLCYB protein abundance and have contributed to the development of red flesh color in domesticated watermelon.

The Cucurbita genus contains several economically important species in the Cucurbitaceae family. Here, we report high-quality genome sequences of C. maxima and C. moschata and provide evidence supporting an allotetraploidization event in Cucurbita. We are able to partition the genome into two homoeologous subgenomes based on different genetic distances to melon, cucumber, and watermelon in the Benincaseae tribe. We estimate that the two diploid progenitors successively diverged from Benincaseae around 31 and 26 million years ago (Mya), respectively, and the allotetraploidization happened at some point between 26 Mya and 3 Mya, the estimated date when C. maxima and C. moschata diverged. The subgenomes have largely maintained the chromosome structures of their diploid progenitors. Such long-term karyotype stability after polyploidization has not been commonly observed in plant polyploids. The two subgenomes have retained similar numbers of genes, and neither subgenome is globally dominant in gene expression. Allele-specific expression analysis in the C. maxima × C. moschata interspecific F1 hybrid and their two parents indicates the predominance of trans-regulatory effects underlying expression divergence of the parents, and detects transgressive gene expression changes in the hybrid correlated with heterosis in important agronomic traits. Our study provides insights into polyploid genome evolution and valuable resources for genetic improvement of cucurbit crops.

Although it has been reported previously that ethylene plays a critical role in sex determination in cucurbit species, how the andromonoecy that carries both the male and hermaphroditic flowers is determined in watermelon is still unknown. Here we showed that the watermelon gene 1-aminocyclopropane-1-carboxylate synthase 4 (CitACS4), expressed specifically in carpel primordia, determines the andromonoecy in watermelon. Among four single nucleotide polymorphism (SNPs) and one InDel identified in the coding region of CitACS4, the C364W mutation located in the conserved box 6 was co-segregated with andromonoecy. Enzymatic analyses showed that the C364W mutation caused a reduced activity in CitACS4. We believe that the reduced CitACS4 activity may hamper the programmed cell death in stamen primordia, leading to the formation of hermaphroditic flowers.

Grafting is an important agricultural technique widely used to improve plant growth, yield, and adaptation to either biotic or abiotic stresses. However, the molecular mechanisms underlying grafting-induced physiological processes remain unclear. Watermelon (Citrullus lanatus L.) is an important horticultural crop worldwide. Grafting technique is commonly used in watermelon production for improving its tolerance to stresses, especially to the soil-borne fusarium wilt disease. In the present study, we used high-throughput sequencing to perform a genome-wide transcript analysis of scions from watermelon grafted onto bottle gourd and squash rootstocks. Our transcriptome and digital gene expression (DGE) profiling data provided insights into the molecular aspects of gene regulation in grafted watermelon. Compared with self-grafted watermelon, there were 787 and 3485 genes differentially expressed in watermelon grafted onto bottle gourd and squash rootstocks, respectively. These genes were associated with primary and secondary metabolism, hormone signaling, transcription factors, transporters, and response to stimuli. Grafting led to changes in expression of these genes, suggesting that they may play important roles in mediating the physiological processes of grafted seedlings. The potential roles of the grafting-responsive mRNAs in diverse biological and metabolic processes were discussed. Obviously, the data obtained in this study provide an excellent resource for unraveling the mechanisms of candidate genes function in diverse biological processes and in environmental adaptation in a graft system.

How sugar transporters regulate sugar accumulation in fruits is poorly understood and particularly so for species storing high-concentration Suc. Accumulation of soluble sugars in watermelon (Citrullus lanatus) fruit, a major quality trait, had been selected during domestication. Still, the molecular mechanisms controlling this quantitative trait are unknown. We resequenced 96 recombinant inbred lines, derived from crossing sweet and unsweet accessions, to narrow down the size of a previously described sugar content quantitative trait locus, which contains a putative Tonoplast Sugar Transporter gene (ClTST2). Molecular and biochemical analyses indicated that ClTST2 encodes a vacuolar membrane protein, whose expression is associated with tonoplast uptake and accumulation of sugars in watermelon fruit flesh cells. We measured fruit sugar content and resequenced the genomic region surrounding ClTST2 in 400 watermelon accessions and associated the most sugar-related significant single-nucleotide polymorphisms (SNPs) to the ClTST2 promoter. Large-scale population analyses strongly suggest increased expression of ClTST2 as a major molecular event in watermelon domestication associated with a selection sweep around the ClTST2 promoter. Further molecular analyses explored the binding of a sugar-induced transcription factor (SUSIWM1) to a sugar-responsive cis-element within the ClTST2 promoter, which contains the quantitative trait locus (QTL) causal SNP. The functional characterization of ClTST2 and its expression regulation by SUSIWM1 provide novel tools to increase sugar sink potency in watermelon and possibly in other vegetable and fruit crops.

Interactions between watermelon and a green fluorescent protein (GFP)-tagged isolate of Fusarium oxysporum f.sp. niveum race 1 (Fon-1) were studied to determine the differences in infection and colonization of watermelon roots in cultivars resistant to and susceptible to Fusarium wilt. The roots of watermelon seedlings were inoculated with a conidial suspension of the GFP-tagged isolate, and confocal laser scanning microscopy was used to visualize colonization, infection and disease development. The initial infection stages were similar in both the resistant and susceptible cultivars, but the resistant cultivar responded differentially after the pathogen had penetrated the root. The pathogen penetrated and colonized resistant watermelon roots, but further fungal advance appeared to be halted, and the fungus did not enter the taproot, suggesting that resistance is initiated postpenetration. However, the tertiary and secondary lateral roots of resistant watermelon also were colonized, although not as extensively as susceptible roots, and the hyphae had penetrated into the central cylinder of lateral roots forming a dense hyphal mat, which was followed by a subsequent collapse of the lateral roots. The initial infection zone for both the wilt-susceptible and wilt-resistant watermelon roots appeared to be the epidermal cells within the root hair zone, which the fungus penetrated directly after forming appressoria. Areas where secondary roots emerged and wounded root tissue also were penetrated preferentially.

Watermelon (Citrullus lanatus (Thunb.) Matsum. & Nakai) is a non-climacteric fruit. The modern sweet-dessert watermelon is the result of years of cultivation and selection for fruits with desirable qualities. To date, the mechanisms of watermelon fruit ripening, and the role of abscisic acid (ABA) in this process, has not been well understood. We quantified levels of free and conjugated ABA contents in the fruits of cultivated watermelon (97103; C. lanatus subsp. vulgaris), semi-wild germplasm (PI179878; C. lanatus subsp. mucosospermus), and wild germplasm (PI296341-FR; C. lanatus subsp. lanatus). Results showed that ABA content in the fruits of 97103 and PI179878 increased during fruit development and ripening, but maintained a low steady state in the center flesh of PI296341-FR fruits. ABA levels in fruits were highest in 97103 and lowest in PI296341-FR, but no obvious differences in ABA levels were observed in seeds of these lines. Examination of 31 representative watermelon accessions, including different C. lanatus subspecies and ancestral species, showed a correlation between soluble solids content (SSC) and ABA levels in ripening fruits. Furthermore, injection of exogenous ABA or nordihydroguaiaretic acid (NDGA) into 97103 fruits promoted or inhibited ripening, respectively. Transcriptomic analyses showed that the expression levels of several genes involved in ABA metabolism and signaling, including Cla009779 (NCED), Cla005404 (NCED), Cla020673 (CYP707A), Cla006655 (UGT) and Cla020180 (SnRK2), varied significantly in cultivated and wild watermelon center flesh. Three SNPs (-738, C/A; -1681, C/T; -1832, G/T) in the promoter region of Cla020673 (CYP707A) and one single SNP (-701, G/A) in the promoter of Cla020180 (SnRK2) exhibited a high level of correlation with SSC variation in the 100 tested accessions. Our results not only demonstrate for the first time that ABA is involved in the regulation of watermelon fruit ripening, but also provide insights into the evolutionary mechanisms of this phenomenon.

CRISPR/Cas9 system can precisely edit genomic sequence and effectively create knockout mutations in T0 generation watermelon plants. Genome editing offers great advantage to reveal gene function and generate agronomically important mutations to crops. Recently, RNA-guided genome editing system using the type II clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein 9 (Cas9) has been applied to several plant species, achieving successful targeted mutagenesis. Here, we report the genome of watermelon, an important fruit crop, can also be precisely edited by CRISPR/Cas9 system. ClPDS, phytoene desaturase in watermelon, was selected as the target gene because its mutant bears evident albino phenotype. CRISPR/Cas9 system performed genome editing, such as insertions or deletions at the expected position, in transfected watermelon protoplast cells. More importantly, all transgenic watermelon plants harbored ClPDS mutations and showed clear or mosaic albino phenotype, indicating that CRISPR/Cas9 system has technically 100% of genome editing efficiency in transgenic watermelon lines. Furthermore, there were very likely no off-target mutations, indicated by examining regions that were highly homologous to sgRNA sequences. Our results show that CRISPR/Cas9 system is a powerful tool to effectively create knockout mutations in watermelon.

Powdery mildew caused by Podosphaera xanthii is an important disease of melon, and race 2F is the predominant race in most areas of China. Resistance to P. xanthii race 2F in melon K7-1 was controlled by a dominant gene, designated Pm-2F, in a 106-member population of recombinant inbred lines derived from K7-1× susceptible K7-2. Using bulked segregant analysis with molecular markers, we have identified two polymorphic simple sequence repeats (SSR) to determine that Pm-2F is located on linkage group II. Comparative genomic analyses using mapped SSR markers and the cucumber genome sequence showed that the melon chromosomal region carrying Pm-2F is homologous to a 288,223 bp genomic region on cucumber chromosome (chr) 1. The SSR markers on chr 1 of cucumber, SSR02734, SSR02733 and CS27 were found linked with Pm-2F. Comparative mapping showed that two SSR markers (SSR02734 and CMBR8) flanked the Pm-2F locus and two nucleotide binding site-leucine-rich repeat resistance genes were identified in the collinear region of cucumber. A cleaved amplified polymorphic sequence (CAPS) marker was developed from the sequence of resistance genes and it delimits the genomic region carrying Pm-2F to 0.8 cM. The evaluation of 165 melon accessions and 13 race differential lines showed that the newly developed CAPS (CAPS-Dde I) marker can be used as a universal marker for effective marker assisted selection in melon powdery mildew resistance breeding. The putative resistance gene cluster provides a potential target site for further fine mapping and cloning of Pm-2F.

Watermelon, Citrullus lanatus Thunb. Matsum. & Nakai is an important vegetable crop worldwide. Due to its narrow genetic base, detection and utilization of the genetic variations, cultivar identification and increasing genetic diversity are some important tasks for watermelon breeders. Molecular markers, especially microsatellites or simple sequence repeats (SSRs) are playing increasingly important roles for these purposes. In the present study, a core set of 23 highly informative SSR markers was developed for watermelon genetic diversity analysis. Based on whole genome sequencing of 17 watermelon inbred lines, we identified 3.9 million single nucleotide polymorphisms (SNPs) which were used to construct a SNP-based dendrogram for the 17 lines. Meanwhile, from the sequenced genome, 13,744 SSRs were developed, of which 704 were placed on a high-resolution watermelon linkage map. To develop the core set SSR markers, 78 of the 704 mapped SSRs were selected as the candidate markers. Using the SNP-based dendrogram as calibration, 23 SSR markers evenly distributed across the genome were identified as the core marker set for watermelon genetic diversity analysis. Each marker was able to detect 2–7 alleles with polymorphism information content values ranging from 0.45 to 0.82. The dendrograms of 17 watermelon lines based on SNPs, the base set of 78 SSRs and the core set of 23 SSRs were highly consistent. The utility of this core set SSRs was demonstrated in 100 commercial watermelon cultivars and elite lines, which could be placed into six clusters that were largely consistent with previous classification based on morphology and parentage data. This core set of SSR markers should be very useful for genotyping and genetic variation analysis in watermelon.

Transcriptome profiling of watermelon during its incompatible interactions with Fusarium oxysporum f.sp. niveum (FON) was performed using an Agilent custom microarray, which contains 15,000 probes representing approximately 8,200 watermelon genes. A total of 24, 275, 596, 598, and 592 genes showed significant differential expression in FON-infected plant roots, as compared with mock-inoculated roots, at 0.5, 1, 3, 5 and 8 days post inoculation (dpi), respectively. Bioinformatics analysis of these differentially expressed genes revealed that during the incompatible interaction between watermelon and FON, the expression of a number of pathogenesis-related (PR) genes, transcription factors, signalling/regulatory genes, and cell wall modification genes, was significantly induced. A number of genes for transporter proteins such as aquaporins were down-regulated, indicating that transporter proteins might contribute to the development of wilt symptoms after FON infection. In the incompatible interaction, most genes involved in biosynthesis of jasmonic acid (JA) were expressed stronger and more sustained than those in a compatible interaction in FON-infected tissues. Similarly, genes associated with shikimate-phenylpropanoid-lignin biosynthesis were also induced during the incompatible interaction, but expression of these genes were not changed or repressed in the compatible interaction. Those results demonstrate that JA biosynthesis and shikimate-phenylpropanoid-lignin pathways might play important roles in watermelon against FON infection and thus provides new insights in understanding the molecular basis and signalling network in watermelon plants in response to FON infection. We also performed confocal imaging of watermelon roots infected with the green fluorescent protein (GFP)-tagged FON1 to revealed histological characteristics of the infection.

Watermelon is an important vegetable crop throughout the world and its worldwide harvested area is 21.7% of that of all vegetables. The genetic diversity of cultivated watermelon is very narrow, which is a major limiting factor of its breeding. Therefore it is necessary to develop new resources that can be used to widen watermelon genetic basis. We have initiated the Watermelon Whole Genome Sequencing project and currently have generated Solexa sequences which represent more than 89.9X coverage of the genome of watermelon cv. '97103' and covered 95% of the whole genome and 98% of the genic region. Further genome sequencing of '97103' using the traditional Sanger technology, comparative genome sequencing of watermelon cv. '90-4304' (10X coverage Solexa sequences), genome annotation is currently underway. In addition, using Roche 454 GS-FLX massively parallel pyrosequencing technology, we have generated 715,059 ESTs with an average length of 250 bp from fruits of '97103' during different developmental stages. These sequences are being assembled and functionally annotated. SSR markers and pathway prediction have been analyzed preliminarily. A comprehensive digital expression map during watermelon fruit development is being created. Furthermore, we have constructed a SSH library from root tissues of watermelon cv. 'PI296341-FR' infected by Fusarium wilt racel and randomly sequenced approximately 3,895 SSH clones. Sequence analysis indicated that 28.3% of the unique ESTs were related to defense response, signal transduction and interaction with the environment. GFP gene is used as the marker to detect the infection process of Fusarium wilt in root tissue. Three hundred twenty-four SSR and EST-SSR markers consist of the watermelon RIL SSR genetic map. The relationships between four Chromosomes and the genome are identified by FISH.

Watermelon grows throughout the world, and is one of the most important cucurbit crops. The draft genome sequence of the East-Asia watermelon cultivar ‘97103’ was published in 2013. The genome sequence allowed the prediction of 23,440 protein coding genes. Comparative genomics analysis showed that the 11 watermelon chromosomes are derived from a seven-chromosome paleohexaploid eudicot ancestor. Twenty watermelon accessions representing three different C. lanatus subspecies were re-sequenced, which produced numerous haplotypes and revealed the extent of genetic diversity and population structure of watermelon germplasm. Preferentially selected genomic regions were identified and several disease resistance genes were found to be lost during domestication. Integrative genomic and transcriptomic analyses identified genes critical to valuable fruit quality traits. The draft watermelon genome sequence represents an important resource for plant research and crop genetic improvement, and also support further evolutionary and comparative genomics studies of the Cucurbitaceae.

Fusarium wilt (FW) caused by Fusarium oxysporum f. sp. niveum (FON) is the major soilborne disease of watermelon (Citrullus lanatus L.). The development and deployment of resistant cultivars is generally considered to be an effective approach to control FW. In this study, an F8 population consisting of 103 recombinant inbred lines derived from a cross between the cultivar 97103 and a wild accession PI 296341-FR was used for FON race 1 and race 2 fungal inoculations. One major QTL on chromosome 1 for FON race 1 resistance was detected with a logarithm of odds of 13.2 and explained phenotypic variation R2 = 48.1 %; two QTLs of FON race 2 resistance on chromosomes 9 and 10 were discovered based on the high-density integrated genetic map we constructed. The nearest molecular marker should be useful for marker-assisted selection of FON race 1 and race 2 resistance. One receptor kinase, one glucan endo-1,3-β-glucosidase precursors and three acidic chitinase located in the FON-1 QTL genomic region. In Qfon2.1 QTL region, one lipoxygenase gene, five receptor-like kinases and four glutathione S-transferase genes are discovered. One arginine biosynthesis bifunctional protein, two receptor kinase proteins and one lipid-transfer protein located in Qfon2.2 QTL region. Based on SNP analysis by using 20 re-sequenced accessions of watermelon and 231-plant F2 population generated from Black Diamond × Calhoun Grey, we developed a SNP marker Chr1SNP_502124 for FON-1 detection. Electronic supplementary material The online version of this article (doi:10.1007/s11032-015-0375-5) contains supplementary material, which is available to authorized users.

Watermelon, Citrullus lanatus, is an important cucurbit crop grown throughout the world. Here we report a high-quality draft genome sequence of the east Asia watermelon cultivar 97103 (2n = 2x = 22) containing 23,440 predicted protein-coding genes. Comparative genomics analysis provided an evolutionary scenario for the origin of the 11 watermelon chromosomes derived from a 7-chromosome paleohexaploid eudicot ancestor. Resequencing of 20 watermelon accessions representing three different C. lanatus subspecies produced numerous haplotypes and identified the extent of genetic diversity and population structure of watermelon germplasm. Genomic regions that were preferentially selected during domestication were identified. Many disease-resistance genes were also found to be lost during domestication. In addition, integrative genomic and transcriptomic analyses yielded important insights into aspects of phloem-based vascular signaling in common between watermelon and cucumber and identified genes crucial to valuable fruit-quality traits, including sugar accumulation and citrulline metabolism.

Background 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. Results 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. Conclusion 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.

As part of our ongoing efforts to sequence and map the watermelon (Citrullus spp.) genome, we have constructed a high density genetic linkage map. The map positioned 234 watermelon genome sequence scaffolds (an average size of 1.41 Mb) that cover about 330 Mb and account for 93.5% of the 353 Mb of the assembled genomic sequences of the elite Chinese watermelon line 97103 (Citrullus lanatus var. lanatus). The genetic map was constructed using an F(8) population of 103 recombinant inbred lines (RILs). The RILs are derived from a cross between the line 97103 and the United States Plant Introduction (PI) 296341-FR (C. lanatus var. citroides) that contains resistance to fusarium wilt (races 0, 1, and 2). The genetic map consists of eleven linkage groups that include 698 simple sequence repeat (SSR), 219 insertion-deletion (InDel) and 36 structure variation (SV) markers and spans ∼800 cM with a mean marker interval of 0.8 cM. Using fluorescent in situ hybridization (FISH) with 11 BACs that produced chromosome-specifc signals, we have depicted watermelon chromosomes that correspond to the eleven linkage groups constructed in this study. The high resolution genetic map developed here should be a useful platform for the assembly of the watermelon genome, for the development of sequence-based markers used in breeding programs, and for the identification of genes associated with important agricultural traits.

BackgroundCucumber,Cucumis sativusL., is an economically and nutritionally important crop of theCucurbitaceaefamily and has long served as a primary model system for sex determination studies. Recently, the sequencing of its whole genome has been completed. However, transcriptome information of this species is still scarce, with a total of around 8,000 Expressed Sequence Tag (EST) and mRNA sequences currently available in GenBank. In order to gain more insights into molecular mechanisms of plant sex determination and provide the community a functional genomics resource that will facilitate cucurbit research and breeding, we performed transcriptome sequencing of cucumber flower buds of two near-isogenic lines, WI1983G, a gynoecious plant which bears only pistillate flowers, and WI1983H, a hermaphroditic plant which bears only bisexual flowers.ResultUsing Roche-454 massive parallel pyrosequencing technology, we generated a total of 353,941 high quality EST sequences with an average length of 175bp, among which 188,255 were from gynoecious flowers and 165,686 from hermaphroditic flowers. These EST sequences, together with ~5,600 high quality cucumber EST and mRNA sequences available in GenBank, were clustered and assembled into 81,401 unigenes, of which 28,452 were contigs and 52,949 were singletons. The unigenes and ESTs were further mapped to the cucumber genome and more than 500 alternative splicing events were identified in 443 cucumber genes. The unigenes were further functionally annotated by comparing their sequences to different protein and functional domain databases and assigned with Gene Ontology (GO) terms. A biochemical pathway database containing 343 predicted pathways was also created based on the annotations of the unigenes. Digital expression analysis identified ~200 differentially expressed genes between flowers of WI1983G and WI1983H and provided novel insights into molecular mechanisms of plant sex determination process. Furthermore, a set of SSR motifs and high confidence SNPs between WI1983G and WI1983H were identified from the ESTs, which provided the material basis for future genetic linkage and QTL analysis.ConclusionA large set of EST sequences were generated from cucumber flower buds of two different sex types. Differentially expressed genes between these two different sex-type flowers, as well as putative SSR and SNP markers, were identified. These EST sequences provide valuable information to further understand molecular mechanisms of plant sex determination process and forms a rich resource for future functional genomics analysis, marker development and cucumber breeding.

Background Pumpkin (Cucurbita maxima Duch.) is an economically important crop belonging to the Cucurbitaceae family. However, very few genomic and genetic resources are available for this species. As part of our ongoing efforts to sequence the pumpkin genome, high-density genetic map is essential for anchoring and orienting the assembled scaffolds. In addition, a saturated genetic map can facilitate quantitative trait locus (QTL) mapping. Results A set of 186 F2 plants derived from the cross of pumpkin inbred lines Rimu and SQ026 were genotyped using the genotyping-by-sequencing approach. Using the SNPs we identified, a high-density genetic map containing 458 bin-markers was constructed, spanning a total genetic distance of 2,566.8 cM across the 20 linkage groups of C. maxima with a mean marker density of 5.60 cM. Using this map we were able to anchor 58 assembled scaffolds that covered about 194.5 Mb (71.7 %) of the 271.4 Mb assembled pumpkin genome, of which 44 (183.0 Mb; 67.4 %) were oriented. Furthermore, the high-density genetic map was used to identify genomic regions highly associated with an important agronomic trait, dwarf vine. Three QTLs on linkage groups (LGs) 1, 3 and 4, respectively, were recovered. One QTL, qCmB2, which was located in an interval of 0.42 Mb on LG 3, explained 21.4 % phenotypic variations. Within qCmB2, one gene, Cma_004516, encoding the gibberellin (GA) 20-oxidase in the GA biosynthesis pathway, had a 1249-bp deletion in its promoter in bush type lines, and its expression level was significantly increased during the vine growth and higher in vine type lines than bush type lines, supporting Cma_004516 as a possible candidate gene controlling vine growth in pumpkin. Conclusions A high-density pumpkin genetic map was constructed, which was used to successfully anchor and orient the assembled genome scaffolds, and to identify QTLs highly associated with pumpkin vine length. The map provided a valuable resource for gene cloning and marker assisted breeding in pumpkin and other related species. The identified vine length QTLs would help to dissect the underlying molecular basis regulating pumpkin vine growth. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2312-8) contains supplementary material, which is available to authorized users.

Carotenoid pigments confer photoprotection and visual attraction and serve as precursors for many important signaling molecules. Herein, the orange-fruited phenotype of a tomato elite inbred line resulting from sharply reduced carotenoid levels and an increased β-carotene-to-lycopene ratio in fruit was shown to be controlled by a single recessive gene, oft3. BSA-Seq combined with fine mapping delimited the oft3 gene to a 71.23 kb interval on chromosome 4, including eight genes. Finally, the oft3 candidate gene SlIDI1, harboring a 116 bp deletion mutation, was identified by genome sequence analysis. Further functional complementation and CRISPR–Cas9 knockout experiments confirmed that SlIDI1 was the gene underlying the oft3 locus. qRT–PCR analysis revealed that the expression of SlIDI1 was highest in flowers and fruit and increased with fruit ripening or flower maturation. SlIDI1 simultaneously produced long and short transcripts by alternative transcription initiation and alternative splicing. Green fluorescent protein fusion expression revealed that the long isoform was mainly localized in plastids and that an N-terminal 59-amino acid extension sequence was responsible for plastid targeting. Short transcripts were identified in leaves and fruit by 5’ RACE and in fruit by 3’ RACE, which produced corresponding proteins lacking transit peptides and/or putative peroxisome targeting sequences, respectively. In SlIDI1 mutant fruit, SlBCH1 transcription involved in β-carotenoid catabolism was obviously suppressed, which may be responsible for the higher β-carotene-to-lycopene ratio and suggested potential feedback regulatory mechanisms involved in carotenoid pathway flux.

Background Modern watermelon (Citrullus lanatus L.) cultivars share a narrow genetic base due to many years of selection for desirable horticultural qualities. Wild subspecies within C. lanatus are important potential sources of novel alleles for watermelon breeding, but successful trait introgression into elite cultivars has had limited success. The application of marker assisted selection (MAS) in watermelon is yet to be realized, mainly due to the past lack of high quality genetic maps. Recently, a number of useful maps have become available, however these maps have few common markers, and were constructed using different marker sets, thus, making integration and comparative analysis among maps difficult. The objective of this research was to use single-nucleotide polymorphism (SNP) anchor markers to construct an integrated genetic map for C. lanatus. Results Under the framework of the high density genetic map, an integrated genetic map was constructed by merging data from four independent mapping experiments using a genetically diverse array of parental lines, which included three subspecies of watermelon. The 698 simple sequence repeat (SSR), 219 insertion-deletion (InDel), 36 structure variation (SV) and 386 SNP markers from the four maps were used to construct an integrated map. This integrated map contained 1339 markers, spanning 798 cM with an average marker interval of 0.6 cM. Fifty-eight previously reported quantitative trait loci (QTL) for 12 traits in these populations were also integrated into the map. In addition, new QTL identified for brix, fructose, glucose and sucrose were added. Some QTL associated with economically important traits detected in different genetic backgrounds mapped to similar genomic regions of the integrated map, suggesting that such QTL are responsible for the phenotypic variability observed in a broad array of watermelon germplasm. Conclusions The integrated map described herein enhances the utility of genomic tools over previous watermelon genetic maps. A large proportion of the markers in the integrated map are SSRs, InDels and SNPs, which are easily transferable across laboratories. Moreover, the populations used to construct the integrated map include all three watermelon subspecies, making this integrated map useful for the selection of breeding traits, identification of QTL, MAS, analysis of germplasm and commercial hybrid seed detection.