Dynamic characteristics of sugar accumulation and related enzyme activities in sweet and non-sweet watermelon fruits

Sugar content largely determines watermelon fruit quality. We compared changes in sugar accumulation and activities of carbohydrate enzymes in the flesh (central portion) and mesocarp of elite sweet watermelon line 97103 (Citrullus lanatus subsp. vulgaris) and exotic non-sweet line PI296341-FR (C. lanatus subsp. lanatus) to elucidate the physiological and biochemical mechanisms of sugar accumulation in watermelon fruit. The major translocated sugars, raffinose and stachyose, were more unloaded into sweet watermelon fruit than non-sweet fruit. During the fruit development, acid α-galactosidase activity was much higher in flesh of 97103 than in mesocarp of 97103, in flesh and mesocarp of PI296341-FR fruit. Insoluble acid invertase activity was higher in 97103 flesh than in 97103 mesocarp, PI296341-FR flesh or mesocarp from 18 days after pollination (DAP) to 34 DAP. Changes in soluble acid invertase activity in 97103 flesh were similar to those in PI296341-FR flesh and mesocarp from 18 DAP to full ripening. Sucrose synthase and sucrose phosphate synthase activities in 97103 flesh were significantly higher than those in 97103 mesocarp and PI296341-FR fruits from 18 to 34 DAP. Only insoluble acid invertase and sucrose phosphate synthase activities were significantly positively correlated with sucrose content in 97103 flesh. Therefore, phloem loading, distribution and metabolism of major translocated sugars, which are controlled by key sugar metabolism enzymes, determine fruit sugar accumulation in sweet and non-sweet watermelon and reflect the distribution diversity of translocated sugars between subspecies.     

Non-synonymous single nucleotide polymorphisms in the watermelon eIF4E gene are closely associated with resistance to <Emphasis Type="Italic">Zucchini yellow mosaic virus</Emphasis>

Zucchini yellow mosaic virus (ZYMV) is one of the most economically important potyviruses infecting cucurbit crops worldwide. Using a candidate gene approach, we cloned and sequenced eIF4E and eIF(iso)4E gene segments in watermelon. Analysis of the nucleotide sequences between the ZYMV-resistant watermelon plant introduction PI 595203 (Citrullus lanatus var. lanatus) and the ZYMV-susceptible watermelon cultivar ‘New Hampshire Midget’ (‘NHM’) showed the presence of single nucleotide polymorphisms (SNPs). Initial analysis of the identified SNPs in association studies indicated that SNPs in the eIF4E, but not eIF(iso)4E, were closely associated to the phenotype of ZYMV-resistance in 70 F₂ and 114 BC_1R progenies. Subsequently, we focused our efforts in obtaining the entire genomic sequence of watermelon eIF4E. Three SNPs were identified between PI 595203 and NHM. One of the SNPs (A241C) was in exon 1 and the other two SNPs (C309A and T554G) were in the first intron of the gene. SNP241 which resulted in an amino acid substitution (proline to threonine) was shown to be located in the critical cap recognition and binding area, similar to that of several plant species resistance to potyviruses. Analysis of a cleaved amplified polymorphism sequence (CAPS) marker derived from this SNP in F₂ and BC_1R populations demonstrated a cosegregation between the CAPS-2 marker and their ZYMV resistance or susceptibility phenotype. When we investigated whether such SNP mutation in the eIF4E was also conserved in several other PIs of C. lanatus var. citroides, we identified a different SNP (A171G) resulting in another amino acid substitution (D71G) from four ZYMV-resistant C. lanatus var. citroides (PI 244018, PI 482261, PI 482299, and PI 482322). Additional CAPS markers were also identified. Availability of all these CAPS markers will enable marker-aided breeding of watermelon for ZYMV resistance.     

Georeferenced phylogenetic analysis of a global collection of wild and cultivated Citrullus species

The geographical origin of watermelon (Citrullus lanatus) remains debated. While a first hypothesis suggests the center of origin to be West Africa, where the endemic sister species C. mucosospermus thrives, a second hypothesis suggests northeastern Africa where the white‐fleshed Sudanese Kordophan melon is cultivated. In this study, we infer biogeographical and haplotype genealogy for C. lanatus, C. mucosospermus, C. amarus, and C. colocynthis using noncoding cpDNA sequences (trnT‐trnL and ndhF‐rpl32 regions) from a global collection of 135 accessions. In total, we identified 38 haplotypes in C. lanatus, C. mucosospermus, C. amarus, and C. colocynthis; of these, 21 were found in Africa and 17 appear endemic to the continent. The least diverse species was C. mucosospermus (5 haplotypes) and the most diverse was C. colocynthis (16 haplotypes). Some haplotypes of C. mucosospermus were nearly exclusive to West Africa, and C. lanatus and C. mucosospermus shared haplotypes that were distinct from those of both C. amarus and C. colocynthis. The results support previous findings that revealed C. mucosospermus to be the closest relative to C. lanatus (including subsp. cordophanus). West Africa, as a center of endemism of C. mucosospermus, is an area of interest in the search of the origin of C. lanatus. This calls for further historical and phylogeographical investigations and wider collection of samples in West and northeastern Africa.     

Development of an integrated 200K SNP genotyping array and application for genetic mapping,genome assembly improvement and genome wide association studies in pear (Pyrus)

Pear (Pyrus; 2n = 34), the third most important temperate fruit crop, has great nutritional and economic value. Despite the availability of many genomic resources in pear, it is challenging to genotype novel germplasm resources and breeding progeny in a timely and cost‐effective manner. Genotyping arrays can provide fast, efficient and high‐throughput genetic characterization of diverse germplasm, genetic mapping and breeding populations. We present here 200K AXIOM^® PyrSNP, a large‐scale single nucleotide polymorphism (SNP) genotyping array to facilitate genotyping of Pyrus species. A diverse panel of 113 re‐sequenced pear genotypes was used to discover SNPs to promote increased adoption of the array. A set of 188 diverse accessions and an F₁ population of 98 individuals from ‘Cuiguan’ × ‘Starkrimson’ was genotyped with the array to assess its effectiveness. A large majority of SNPs (166 335 or 83%) are of high quality. The high density and uniform distribution of the array SNPs facilitated prediction of centromeric regions on 17 pear chromosomes, and significantly improved the genome assembly from 75.5% to 81.4% based on genetic mapping. Identification of a gene associated with flowering time and candidate genes linked to size of fruit core via genome wide association studies showed the usefulness of the array in pear genetic research. The newly developed high‐density SNP array presents an important tool for rapid and high‐throughput genotyping in pear for genetic map construction, QTL identification and genomic selection.     

A high resolution genetic map anchoring scaffolds of the sequenced watermelon genome

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.     

The high‐quality genome of Brassica napus cultivar ‘ZS11’ reveals the introgression history in semi‐winter morphotype

Allotetraploid oilseed rape (Brassica napus L.) is an agriculturally important crop. Cultivation and breeding of B. napus by humans has resulted in numerous genetically diverse morphotypes with optimized agronomic traits and ecophysiological adaptation. To further understand the genetic basis of diversification and adaptation, we report a draft genome of an Asian semi‐winter oilseed rape cultivar ‘ZS11’ and its comprehensive genomic comparison with the genomes of the winter‐type cultivar ‘Darmor‐bzh’ as well as two progenitors. The integrated BAC‐to‐BAC and whole‐genome shotgun sequencing strategies were effective in the assembly of repetitive regions (especially young long terminal repeats) and resulted in a high‐quality genome assembly of B. napus ‘ZS11’. Within a short evolutionary period (~6700 years ago), semi‐winter‐type ‘ZS11’ and the winter‐type ‘Darmor‐bzh’ maintained highly genomic collinearity. Even so, certain genetic differences were also detected in two morphotypes. Relative to ‘Darmor‐bzh’, both two subgenomes of ‘ZS11’ are closely related to its progenitors, and the ‘ZS11’ genome harbored several specific segmental homoeologous exchanges (HEs). Furthermore, the semi‐winter‐type ‘ZS11’ underwent potential genomic introgressions with B. rapa (A_r). Some of these genetic differences were associated with key agronomic traits. A key gene of A03.FLC3 regulating vernalization‐responsive flowering time in ‘ZS11’ was first experienced HE, and then underwent genomic introgression event with A_r, which potentially has led to genetic differences in controlling vernalization in the semi‐winter types. Our observations improved our understanding of the genetic diversity of different B. napus morphotypes and the cultivation history of semi‐winter oilseed rape in Asia.     

The pomegranate (Punica granatum L.) draft genome dissects genetic divergence between soft‐ and hard‐seeded cultivars

Complete and highly accurate reference genomes and gene annotations are indispensable for basic biological research and trait improvement of woody tree species. In this study, we integrated single‐molecule sequencing and high‐throughput chromosome conformation capture techniques to produce a high‐quality and long‐range contiguity chromosome‐scale genome assembly of the soft‐seeded pomegranate cultivar ‘Tunisia’. The genome covers 320.31 Mb (scaffold N50 = 39.96 Mb; contig N50 = 4.49 Mb) and includes 33 594 protein‐coding genes. We also resequenced 26 pomegranate varieties that varied regarding seed hardness. Comparative genomic analyses revealed many genetic differences between soft‐ and hard‐seeded pomegranate varieties. A set of selective loci containing SUC8‐like, SUC6, FoxO and MAPK were identified by the selective sweep analysis between hard‐ and soft‐seeded populations. An exceptionally large selective region (26.2 Mb) was identified on chromosome 1. Our assembled pomegranate genome is more complete than other currently available genome assemblies. Our results indicate that genomic variations and selective genes may have contributed to the genetic divergence between soft‐ and hard‐seeded pomegranate varieties.     

Use of VeraCode 384-plex assays for watermelon diversity analysis and integrated genetic map of watermelon with single nucleotide polymorphisms and simple sequence repeats

Watermelon (Citrullus lanatus var. lanatus) is one of the most important vegetable crops in the world. Molecular markers have become the tools of choice for resolving watermelon taxonomic relationships and evolution. Increased numbers of single nucleotide polymorphism (SNP) markers together with simple sequence repeat (SSR) markers would be useful for phylogenetic analyses of germplasm accessions and for linkage mapping for marker-assisted breeding with quantitative trait loci and single genes. We aimed to construct a genetic map based on SNPs (generated by Illumina Veracode multiplex assays for genotyping) and SSR markers and evaluate relationships inferred from SNP genotypes between 130 watermelon accessions collected throughout the world. We incorporated 282 markers (232 SNPs and 50 SSRs) into the linkage map. The genetic map consisted of 11 linkage groups spanning 924.72 cM with an average distance of 3.28 cM between markers. Because all of the SNP-containing sequences were assembled with the whole-genome sequence draft for watermelon, chromosome numbers could be readily assigned for all the linkage groups. We found that 134 SNPs were polymorphic in 130 watermelon accessions chosen for diversity studies. The current 384-plex SNP set is a powerful tool for characterizing genetic relatedness and for developing medium-resolution genetic maps.     

Construction of a watermelon BAC library and identification of SSRs anchored to melon or Arabidopsis genomes

A bacterial artificial chromosome (BAC) library was constructed for watermelon (Citrullus lanatus (Thunb.) Matsum. & Nakai var. lanatus) with an average insert-size of 106 kb, providing 21 haploid genome equivalents. The library was used to identify BAC clones that are anchored to probes evenly distributed on the genomes of melon or Arabidopsis. Twenty eight probes (representing 66% of the tested probes) from melon and 30 probes (65%) from Arabidopsis identified positive BAC clones. Two methods were implemented to identify SSRs from the positively hybridizing BAC clones. First, analysis of BAC end sequences revealed 37 SSRs. For the second method, pooled DNA of BACs identified by the melon probes was used to develop a shotgun library. The library was then screened with synthetic SSR oligonucleotides by hybridization. Sequence analysis of positively hybridizing shotgun clones revealed 142 different SSRs. Thirty eight SSRs were characterized using three watermelon cultivars, five plant introduction (PI) accessions of C. lanatus var lanatus and four PIs of C. lanatus var citroides. Of these, 36 (95%) were found to be polymorphic with up to six alleles per marker. Polymorphism information content values for polymorphic markers varied between 0.22 and 0.79 with an average of 0.53. The methods described herein will be valuable for the construction of a watermelon linkage map with SSRs evenly distributed on its genome that is anchored to the genomes of melon and Arabidopsis.     

Single nucleotide polymorphisms generated by genotyping by sequencing to characterize genome-wide diversity,linkage disequilibrium,and selective sweeps in cultivated watermelon

Background

A large single nucleotide polymorphism (SNP) dataset was used to analyze genome-wide diversity in a diverse collection of watermelon cultivars representing globally cultivated, watermelon genetic diversity. The marker density required for conducting successful association mapping depends on the extent of linkage disequilibrium (LD) within a population. Use of genotyping by sequencing reveals large numbers of SNPs that in turn generate opportunities in genome-wide association mapping and marker-assisted selection, even in crops such as watermelon for which few genomic resources are available. In this paper, we used genome-wide genetic diversity to study LD, selective sweeps, and pairwise F_ST distributions among worldwide cultivated watermelons to track signals of domestication.

Results

We examined 183 Citrullus lanatus var. lanatus accessions representing domesticated watermelon and generated a set of 11,485 SNP markers using genotyping by sequencing. With a diverse panel of worldwide cultivated watermelons, we identified a set of 5,254 SNPs with a minor allele frequency of ≥ 0.05, distributed across the genome. All ancestries were traced to Africa and an admixture of various ancestries constituted secondary gene pools across various continents. A sliding window analysis using pairwise F_ST values was used to resolve selective sweeps. We identified strong selection on chromosomes 3 and 9 that might have contributed to the domestication process. Pairwise analysis of adjacent SNPs within a chromosome as well as within a haplotype allowed us to estimate genome-wide LD decay. LD was also detected within individual genes on various chromosomes. Principal component and ancestry analyses were used to account for population structure in a genome-wide association study. We further mapped important genes for soluble solid content using a mixed linear model.

Conclusions

Information concerning the SNP resources, population structure, and LD developed in this study will help in identifying agronomically important candidate genes from the genomic regions underlying selection and for mapping quantitative trait loci using a genome-wide association study in sweet watermelon.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-767) contains supplementary material, which is available to authorized users.     

Resurgence of cucurbit downy mildew in the United States: Insights from comparative genomic analysis of Pseudoperonospora cubensis

Pseudoperonospora cubensis, the causal agent of cucurbit downy mildew (CDM), is known to exhibit host specialization. The virulence of different isolates of the pathogen can be classified into pathotypes based on their compatibility with a differential set composed of specific cucurbit host types. However, the genetic basis of host specialization within P. cubensis is not yet known. Total genomic DNA extracted from nine isolates of P. cubensis collected from 2008 to 2013 from diverse cucurbit host types (Cucumis sativus, C. melo var. reticulatus, Cucurbita maxima, C. moschata, C. pepo, and Citrullus lanatus) in the United States were subjected to whole‐genome sequencing. Comparative analysis of these nine genomes confirmed the presence of two distinct evolutionary lineages (lineages I and II) of P. cubensis. Many fixed polymorphisms separated lineage I comprising isolates from Cucurbita pepo, C. moschata, and Citrullus lanatus from lineage II comprising isolates from Cucumis spp. and Cucurbita maxima. Phenotypic characterization showed that lineage II isolates were of the A1 mating type and belonged to pathotypes 1 and 3 that were not known to be present in the United States prior to the resurgence of CDM in 2004. The association of lineage II isolates with the new pathotypes and a lack of genetic diversity among these isolates suggest that lineage II of P. cubensis is associated with the resurgence of CDM on cucumber in the United States.     

Rapid Detection of Colletotrichum lagenarium,Causal Agent of Anthracnose of Cucurbitaceous Crops,by PCR and Real‐time PCR

Rapid and accurate polymerase chain reaction (PCR) and real‐time PCR methods were developed for the detection of Colletotrichum lagenarium, the causal agent of anthracnose, in tissues of squash (Cucurbita moschata), watermelon (Citrullus lanatus), cucumber (Cucumis sativus) and muskmelon (Cucumis melo). PCR assays amplified different internal transcribed spacer sequences from C. lagenarium, so effectively detected this pathogen in infected tissues. PCR analysis with the primer co‐m‐337F1/R1 was able to differentiate C. lagenarium from other fungal pathogens, including Colletotrichum spp., Fusarium spp., Alternaria spp. and Didymella spp. An optimized real‐time PCR assay was developed to detect and monitor C. lagenarium in both infected plant tissues and soil samples. The sensitivity of real‐time PCR can detect down to 1 pg of DNA. Thus, PCR‐based analysis is a useful technique for rapid detection and diagnosis of C. lagenarium in infected plants or infested soils.     

An archaeobotanical contribution to the history of watermelon,<Emphasis Type="Italic"> Citrullus lanatus</Emphasis> (Thunb.) Matsum. &; Nakai (syn.<Emphasis Type="Italic"> C. vulgaris</Emphasis> Schrad.)

The discovery of several 5000-year old seeds of wild watermelon, Citrullus lanatus, at an archaeological site Uan Muhuggiag in southwest Libya, re-opens the debate on the origin, wild distribution and domestication history of this species. The seeds were found within a plant assemblage of wild seeds and fruits, associated with pottery and bones of domestic animals belonging to Neolithic pastoralists. The presumed wild progenitor of the modern cultivar C. lanatus is today found exclusively in a region centring on the Kalahari Desert. This new archaeobotanical record raises the possibility that this distribution was much more extensive in the past.     

Mitigation of saline conditions in watermelon with mycorrhiza and silicon application

Citrullus lanatus L. is critical vegetable for salinity stress. Arbuscular mycorrhizal fungi (AMF) and silicon treatments are known to help as bio-ameliorator of saline soils that can improve salinity tolerance in plants. But their combined effect has never been examined on watermelon therefore, present study investigated the effect of inoculation with the Arbuscular mycorrhizal fungi (AMF) along with silicon on the growth and yield parameters, antioxidant enzyme activities, pigment and mineral content of Citrullus lanatus L. plants grown during salt stress conditions. Outcomes from the study point out that salt stressed watermelon plants showed the best morphological and biochemical values when inoculated with Silicon (4 mM) + Glomus mosseae + Gigaspora gigantean. In addition, the plants inoculated by similar treatment demonstrated less osmotic activity, electrolyte leakage, as well as peroxide content. Treatments comprising Silicon (4 mM) with either Glomus mosseae and Gigaspora gigantean also performed significantly similar for most of the traits studied in the present investigation and better than the treatment only with either one of Glomus mosseae and Gigaspora gigantean. Antioxidant efficiency of melon was certainly appreciably enhanced after incubation with AMF and Si combination in salinity stress. Overall, the application of mycorrhiza and silicon can be considered to overcome the salinity stress in watermelon.     

Identification of a Novel 1033‐Nucleotide Deletion Polymorphism in the Prophage Region of ‘Candidatus Liberibacter asiaticus’: Potential Applications for Bacterial Epidemiology

The prophage/phage region in the genome of ‘Candidatus Liberibacter asiaticus’, an alpha‐proteobacterium associated with citrus Huanglongbing, included many valuable loci for genetic diversity studies. Previously, a mosaic genomic region (CLIBASIA_05640 to CLIBASIA_05650) was characterized, and this revealed inter‐ and intracontinental variations of ‘Ca. L. asiaticus’. In this study, 267 ‘Ca. L. asiaticus’ isolates collected from eight provinces in China were analysed with a primer set flanking the same mosaic region plus downstream sequence. While most amplicon sizes ranged from 1400 to 2000 bp, an amplicon of 550 bp (S550) was found in 14 samples collected from south‐western China. Sequence analyses showed that S550 was the result of a 1033 bp deletion which included the previously known mosaic region. The genetic nature of the deletion event remains unknown. The regional restriction of S550 suggests that the ‘Ca. L. asiaticus’ population from south‐western China is different from those in eastern China. The small and easy‐to‐detect S550 amplicon could serve as a molecular marker for ‘Ca. L. asiaticus’ epidemiology.     

A new genomic tool for walnut (Juglans regia L.): development and validation of the high‐density Axiom™ J. regia 700K SNP genotyping array

Over the last 20 years, global production of Persian walnut (Juglans regia L.) has grown enormously, likely reflecting increased consumption due to its numerous benefits to human health. However, advances in genome‐wide association (GWA) studies and genomic selection (GS) for agronomically important traits in walnut remain limited due to the lack of powerful genomic tools. Here, we present the development and validation of a high‐density 700K single nucleotide polymorphism (SNP) array in Persian walnut. Over 609K high‐quality SNPs have been thoroughly selected from a set of 9.6 m genome‐wide variants, previously identified from the high‐depth re‐sequencing of 27 founders of the Walnut Improvement Program (WIP) of University of California, Davis. To validate the effectiveness of the array, we genotyped a collection of 1284 walnut trees, including 1167 progeny of 48 WIP families and 26 walnut cultivars. More than half of the SNPs (55.7%) fell in the highest quality class of ‘Poly High Resolution’ (PHR) polymorphisms, which were used to assess the WIP pedigree integrity. We identified 151 new parent‐offspring relationships, all confirmed with the Mendelian inheritance test. In addition, we explored the genetic variability among cultivars of different origin, revealing how the varieties from Europe and California were differentiated from Asian accessions. Both the reconstruction of the WIP pedigree and population structure analysis confirmed the effectiveness of the Applied Biosystems? Axiom? J. regia 700K SNP array, which initiates a novel genomic and advanced phase in walnut genetics and breeding.     

Chromosomal‐level assembly of Takifugu obscurus (Abe, 1949) genome using third‐generation DNA sequencing and Hi‐C analysis

The Tetraodontidae family are known to have relatively small and compact genomes compared to other vertebrates. The obscure puffer fish Takifugu obscurus is an anadromous species that migrates to freshwater from the sea for spawning. Thus the euryhaline characteristics of T. obscurus have been investigated to gain understanding of their survival ability, osmoregulation, and other homeostatic mechanisms in both freshwater and seawater. In this study, a high quality chromosome‐level reference genome for T. obscurus was constructed using long‐read Pacific Biosciences (PacBio) Sequel sequencing and a Hi‐C‐based chromatin contact map platform. The final genome assembly of T. obscurus is 381 Mb, with a contig N50 length of 3,296 kb and longest length of 10.7 Mb, from a total of 62 Gb of raw reads generated using single‐molecule real‐time sequencing technology from a PacBio Sequel platform. The PacBio data were further clustered into chromosome‐scale scaffolds using a Hi‐C approach, resulting in a 373 Mb genome assembly with a contig N50 length of 15.2 Mb and and longest length of 28 Mb. When we directly compared the 22 longest scaffolds of T. obscurus to the 22 chromosomes of the tiger puffer Takifugu rubripes, a clear one‐to‐one orthologous relationship was observed between the two species, supporting the chromosome‐level assembly of T. obscurus. This genome assembly can serve as a valuable genetic resource for exploring fugu‐specific compact genome characteristics, and will provide essential genomic information for understanding molecular adaptations to salinity fluctuations and the evolution of osmoregulatory mechanisms.