Genome of 'Charleston Gray', the principal American watermelon cultivar, and genetic characterization of 1,365 accessions in the U.S. National Plant Germplasm System watermelon collection

Chromosome-level genome assembly and population genomics reveals crucial selection for subgynoecy development in chieh-qua

Chieh-qua is an important cucurbit crop and very popular in South China and Southeast Asia. Despite its significance, its genetic basis and domestication history are unclear. In this study, we have successfully generated a chromosome-level reference genome assembly for the chieh-qua 'A36' using a hybrid assembly strategy that combines PacBio long reads and Illumina short reads. The assembled genome of chieh-qua is approximately 953.3 Mb in size and is organized into 12 chromosomes, with contig N50 of 6.9 Mb and scaffold N50 of 68.2 Mb. Notably, the chieh-qua genome is comparable in size to the wax gourd genome. Through gene prediction analysis, we have identified a total of 24 593 protein-coding genes in the A36 genome. Additionally, approximately 56.6% (539.3 Mb) of the chieh-qua genome consists of repetitive sequences. Comparative genome analysis revealed that chieh-qua and wax gourd are closely related, indicating a close evolutionary relationship between the two species. Population genomic analysis, employing 129 chieh-qua accessions and 146 wax gourd accessions, demonstrated that chieh-qua exhibits greater genetic diversity compared to wax gourd. We also employed the GWAS method to identify related QTLs associated with subgynoecy, an interested and important trait in chieh-qua. The MYB59 (BhiCQ0880026447) exhibited relatively high expression levels in the shoot apex of four subgynoecious varieties compared with monoecious varieties. Overall, this research provides insights into the domestication history of chieh-qua and offers valuable genomic resources for further molecular research.

Mapping the genetic architecture of low-temperature stress tolerance in citron watermelon

Sweet-fleshed watermelon (Citrullus lanatus) is an important vegetable crop of the tropical origin. It is widely grown and consumed around the world for its hydration and nutritional quality values. Low-temperature stress can affect early planting, seedling establishment, and expansion of crop production to new areas. A collection of 122 citron watermelon (Citrullus amarus) accessions were obtained from the USDA's National Plant Germplasm Repository System gene bank in Griffin, GA. The accessions were genotyped using whole genome resequencing to generate single nucleotide polymorphisms (SNPs) molecular markers and screened under cold-stressed and non-stressed control conditions. Four low-temperature stress tolerance related traits including shoot biomass, vine length, maximum quantum efficiency of photosystem II, and chlorophyll content were measured under cold-stressed and non-stressed control treatment conditions. Correlation analysis revealed the presence of positive relationships among traits. Broad-sense heritability for all traits ranged from 0.35 to 0.73, implying the presence of genetic contributions to the observed phenotypic variation. Genomic regions underlying these traits across several citron watermelon chromosomes were identified. Four low-temperature stress tolerance related putative candidate genes co-located with the peak SNPs from genome-wide association study. These genomic regions and marker information could potentially be used in molecular breeding to accelerate genetic improvements for low-temperature stress tolerance in watermelon.

Global characteristics and trends of researches on watermelon: Based on bibliometric and visualized analysis

Watermelon is an important horticultural plant. A bibliometric analysis of the watermelon literature was carried out in order to analyze the research state, hotspots, and trends, as well as to highlight the overall watermelon research development from a holistic viewpoint. The summary of watermelon research is given via metrological analysis based on a set of indices using a newly built Bibliometrix R-package tool. This study gathered 6,632 documents indexed in the Core Collection of Web of Science (WoS) in the domain of watermelon from 1992 to 2022 using bibliometrix. The results indicated that the number of published articles showed an apparently upward trend. The United States was in the first place, with Plant Disease being the most productive journal. Levi A from the United States Department of Agriculture-Agricultural Research Service is the most prolific author, and Levi A is the most cited; The most frequently used keywords by authors are "growth", "resistance", "identification", "yield", "quality" "plants", "watermelon stomach" and "expression"; The most talked-about issues in this subject are resistance, yield, and quality, which highlight the crucial research areas. To effectively comprehend the turning moments for future research, it is useful to monitor the hotspots and frontiers of watermelon studies. The results highlight the future paths for study in the field of watermelon and provide useful information for researchers interested in the topic.

Development of a Gene-Based Marker Set for Orange-Colored Watermelon Flesh with a High β-Carotene Content

The fruit flesh of watermelons differs depending on the distinct carotenoid composition. Orange-colored flesh relates to the accumulation of β-carotene, which is beneficial to human health. Canary-yellow-fleshed OTO-DAH and orange-β-fleshed (orange-fleshed with high β-carotene) NB-DAH near-isogenic lines (NILs) were used to determine the genetic mechanism attributed to orange watermelon flesh. For genetic mapping, an F2 population was developed by crossing the two NILs. The segregation ratio of flesh color in the F2 population indicated that the orange-β flesh of the NB-DAH NIL was controlled by a single incompletely dominant gene. Through a comparative analysis of the whole-genome sequences of the parent lines and NILs, a major introgression region unique to the NB-DAH NIL was detected on Chr. 1; this was considered a candidate region for harboring genes that distinguish orange from canary-yellow and red flesh. Among the 13 genes involved in the carotenoid metabolic pathway in watermelons, only ClPSY1 (ClCG01G008470), which encodes phytoene synthase 1, was located within the introgression region. The genotyping of F2 plants using a cleaved amplified polymorphic sequence marker developed from a non-synonymous SNP in ClPSY1 revealed its relationship with orange-β flesh. The insights gained in this study can be applied to marker-assisted breeding for this desirable trait.

Genome-wide association study of soluble solids content, flesh color, and fruit shape in citron watermelon

Fruit quality traits are crucial determinants of consumers' willingness to purchase watermelon produce, making them major goals for breeding programs. There is limited information on the genetic underpinnings of fruit quality traits in watermelon. A total of 125 citron watermelon (Citrullus amarus) accessions were genotyped using single nucleotide polymorphisms (SNPs) molecular markers generated via whole-genome resequencing. A total of 2,126,759 genome-wide SNP markers were used to uncover marker-trait associations using single and multi-locus GWAS models. High broad-sense heritability for fruit quality traits was detected. Correlation analysis among traits revealed positive relationships, with the exception of fruit diameter and fruit shape index (ratio of fruit length to fruit diameter), which was negative. A total of 37 significant SNP markers associated with soluble solids content, flesh color, fruit length, fruit diameter, and fruit shape index traits were uncovered. These peak SNPs accounted for 2.1%-23.4% of the phenotypic variation explained showing the quantitative inheritance nature of the evaluated traits. Candidate genes relevant to fruit quality traits were uncovered on chromosomes Ca01, Ca03, Ca06, and Ca07. These significant molecular markers and candidate genes will be useful in marker-assisted breeding of fruit quality traits in watermelon.

Molecular Mapping of Putative Genomic Regions Controlling Fruit and Seed Morphology of Watermelon

The genetic regulatory basis of qualitative and quantitative phenotypes of watermelon is being investigated in different types of molecular and genetic breeding studies around the world. In this study, biparental F2 mapping populations were developed over two experimental years, and the collected datasets of fruit and seed traits exhibited highly significant correlations. Whole-genome resequencing of comparative parental lines was performed and detected single nucleotide polymorphism (SNP) loci were converted into cleaved amplified polymorphic sequence (CAPS) markers. The screened polymorphic markers were genotyped in segregating populations and two genetic linkage maps were constructed, which covered a total of 2834.28 and 2721.45 centimorgan (cM) genetic lengths, respectively. A total of 22 quantitative trait loci (QTLs) for seven phenotypic traits were mapped; among them, five stable and major-effect QTLs (PC-8-1, SL-9-1, SWi-9-1, SSi-9-1, and SW-6-1) and four minor-effect QTLs (PC-2-1 and PC-2-2; PT-2-1 and PT-2-2; SL-6-1 and SSi-6-2; and SWi-6-1 and SWi-6-2) were observed with 3.77-38.98% PVE. The adjacent QTL markers showed a good fit marker-trait association, and a significant allele-specific contribution was also noticed for genetic inheritance of traits. Further, a total of four candidate genes (Cla97C09G179150, Cla97C09G179350, Cla97C09G180040, and Cla97C09G180100) were spotted in the stable colocalized QTLs of seed size linked traits (SL-9-1 and SWi-9-1) that showed non-synonymous type mutations. The gene expression trends indicated that the seed morphology had been formed in the early developmental stage and showed the genetic regulation of seed shape formation. Hence, we think that our identified QTLs and genes would provide powerful genetic insights for marker-assisted breeding aimed at improving the quality traits of watermelon.

Genetic control of flowering time and fruit yield in citron watermelon

Flowering time and fruit yield are important traits in watermelon crop improvement. There is limited information on the inheritance and genomic loci underlying flowering time and yield performance, especially in citron watermelon. A total of 125 citron watermelon accessions were evaluated in field trials over two growing seasons for days to male and female flowers, fruit count, fruit weight, and fruit yield. The germplasm was genotyped with more than two million single-nucleotide polymorphism (SNP) markers generated via whole-genome resequencing. Trait mapping was conducted using a genome-wide association study (GWAS). Broad-sense heritability for all traits ranged from moderate to high, indicating that genetic improvement through breeding and selection is feasible. Significant marker-trait associations were uncovered for days to female flower (chromosomes Ca04, Ca05, Ca08, and Ca09), fruit count (on Ca02, Ca03, and Ca05), fruit weight (on Ca02, Ca06, Ca08, Ca10, and Ca11), and fruit yield on chromosomes Ca05, Ca07, and Ca09. The phenotypic variation explained by the significant SNPs ranged from 1.6 to 25.4, highlighting the complex genetic architecture of the evaluated traits. Candidate genes relevant to flowering time and fruit yield component traits were uncovered on chromosomes Ca02, Ca04, Ca05, Ca06, Ca09, and Ca11. These results lay a foundation for marker-assisted trait introgression of flowering time and fruit yield component traits in watermelons.

Identification of a novel locus C2 controlling canary yellow flesh color in watermelons

The flesh color of watermelon is an important trait that is determined by carotenoid composition and affects consumers' fruit desirability. Although a complete dominant control by C locus (Cllcyb) for canary yellow flesh (CY) over red flesh has been reported, red and CY colors frequently appear as a mixed pattern in the same flesh (incomplete canary yellow, ICY) in F1 and inbred lines carrying dominant C alleles. Therefore, we examined the genetic control of the mixed color pattern in ICY using whole-genome resequencing of three ICY (ICY group) and three CY inbred lines (CY group), as well as genetic linkage mapping of an F2 population. The segregation pattern in 135 F2 plants indicated that CY is controlled by a single locus (named C 2) dominant over ICY. The whole-genome resequencing of ICY and CY inbred lines revealed an ICY/CY-specific region of approximately 27.60-27.88 Mb on Chr. 2 that was polymorphic between the ICY and CY groups. Our genetic map, using nine cleaved amplified polymorphic sequence markers developed based on the single-nucleotide polymorphisms from the ICY/CY-specific region, confirmed that C 2 is located on Chr. 2 and cosegregated with the marker (M7) derived from a non-synonymous single-nucleotide polymorphism of the pentatricopeptide repeat (PPR) gene (ClPPR, Cla97C02G039880). Additionally, 27 watermelon inbred lines of ICY, CY, and red flesh were evaluated using previously reported Cllcyb (C locus)-based markers and our C 2 locus-linked ClPPR-based marker (M7). As a result, dominant alleles at the C 2 locus were required to produce CY, in addition to dominant alleles at the C locus, while a recessive homozygous genotype at the C locus gave the red flesh irrespective of the genotype at the C 2 locus. Using a ClPPR-based cleaved amplified polymorphic sequence developed in this study and Cllcyb-based markers, watermelon cultivars with CY, ICY, and red flesh could be successfully discerned, implying that the combined use of these markers will be efficient for marker-assisted selection of flesh color in watermelon breeding.

Genome-wide characterization of the PLATZ gene family in watermelon (Citrullus lanatus L.) with putative functions in biotic and abiotic stress response

Plant AT-rich sequence and zinc-binding (PLATZ) proteins are plant-specific transcription factors involved in growth, development, and stress responses. Here, we conducted a genome-wide characterization of the watermelon ClPLATZ family and examined its expression responsiveness to defense hormones and pathogen infection along with putative functions in biotic and abiotic stress responses. The watermelon genome contains 12 putative ClPLATZ genes, encoding proteins with a characteristic PLATZ domain, and their promoters contain various cis-elements related to plant growth, development, phytohormones and stress response. The ClPLATZ genes, except ClPLATZ6, are differentially expressed in response to defense hormones (e.g., salicylic acid and methyl jasmonate) and fungal infections caused by Fusarium oxysporum f. sp. niveum and Stagonosporopsis cucurbitacearum. Most ClPLATZ proteins interact with other proteins (viz., ClDP, ClRPT2a, and ClRPC53). Among ClPLATZ proteins, ClPLATZ8, 9, 10, and 11 are predominately localized in the nucleus. ClPLATZ3 and 8 positively, but ClPLATZ11 negatively regulate resistance against Pseudomonas syringe pv. tomato DC3000 in transgenic Arabidopsis lines. ClPLATZ8 and 11 positively regulate stress tolerance to NaCl and mannitol during seed germination in transgenic Arabidopsis. In conclusion, the characterization of the ClPLATZ family provides insights into the biological functions of ClPLATZ genes in growth, development, and stress response in watermelon. Further, the involvement of certain ClPLATZ genes in biotic and abiotic stress response in transgenic Arabidopsis suggests their potential application in engineering stress-tolerant crops.

Exploring the sorghum race level diversity utilizing 272 sorghum accessions genomic resources

Due to evolutionary divergence, sorghum race populations exhibit significant genetic and morphological variation. A k-mer-based sorghum race sequence comparison identified the conserved k-mers of all 272 accessions from sorghum and the race-specific genetic signatures identified the gene variability in 10,321 genes (PAVs). To understand sorghum race structure, diversity and domestication, a deep learning-based variant calling approach was employed in a set of genotypic data derived from a diverse panel of 272 sorghum accessions. The data resulted in 1.7 million high-quality genome-wide SNPs and identified selective signature (both positive and negative) regions through a genome-wide scan with different (iHS and XP-EHH) statistical methods. We discovered 2,370 genes associated with selection signatures including 179 selective sweep regions distributed over 10 chromosomes. Co-localization of these regions undergoing selective pressure with previously reported QTLs and genes revealed that the signatures of selection could be related to the domestication of important agronomic traits such as biomass and plant height. The developed k-mer signatures will be useful in the future to identify the sorghum race and for trait and SNP markers for assisting in plant breeding programs.

Characterization of the USDA Cucurbita pepo, C. moschata, and C. maxima germplasm collections

The Cucurbita genus is home to a number of economically and culturally important species. We present the analysis of genotype data generated through genotyping-by-sequencing of the USDA germplasm collections of Cucurbita pepo, C. moschata, and C. maxima. These collections include a mixture of wild, landrace, and cultivated specimens from all over the world. Roughly 1,500 - 32,000 high-quality single nucleotide polymorphisms (SNPs) were called in each of the collections, which ranged in size from 314 to 829 accessions. Genomic analyses were conducted to characterize the diversity in each of the species. Analysis revealed extensive structure corresponding to a combination of geographical origin and morphotype/market class. Genome-wide associate studies (GWAS) were conducted using both historical and contemporary data. Signals were observed for several traits, but the strongest was for the bush (Bu) gene in C. pepo. Analysis of genomic heritability, together with population structure and GWAS results, was used to demonstrate a close alignment of seed size in C. pepo, maturity in C. moschata, and plant habit in C. maxima with genetic subgroups. These data represent a large, valuable collection of sequenced Cucurbita that can be used to direct the maintenance of genetic diversity, for developing breeding resources, and to help prioritize whole-genome re-sequencing.

A comprehensive and conceptual overview of omics-based approaches for enhancing the resilience of vegetable crops against abiotic stresses

Abiotic stresses adversely affect the productivity and production of vegetable crops. The increasing number of crop genomes that have been sequenced or re-sequenced provides a set of computationally anticipated abiotic stress-related responsive genes on which further research may be focused. Knowledge of omics approaches and other advanced molecular tools have all been employed to understand the complex biology of these abiotic stresses. A vegetable can be defined as any component of a plant that is eaten for food. These plant parts may be celery stems, spinach leaves, radish roots, potato tubers, garlic bulbs, immature cauliflower flowers, cucumber fruits, and pea seeds. Abiotic stresses, such as deficient or excessive water, high temperature, cold, salinity, oxidative, heavy metals, and osmotic stress, are responsible for the adverse activity in plants and, ultimately major concern for decreasing yield in many vegetable crops. At the morphological level, altered leaf, shoot and root growth, altered life cycle duration and fewer or smaller organs can be observed. Likewise different physiological and biochemical/molecular processes are also affected in response to these abiotic stresses. In order to adapt and survive in a variety of stressful situations, plants have evolved physiological, biochemical, and molecular response mechanisms. A comprehensive understanding of the vegetable's response to different abiotic stresses and the identification of tolerant genotypes are essential to strengthening each vegetable's breeding program. The advances in genomics and next-generation sequencing have enabled the sequencing of many plant genomes over the last twenty years. A combination of modern genomics (MAS, GWAS, genomic selection, transgenic breeding, and gene editing), transcriptomics, and proteomics along with next-generation sequencing provides an array of new powerful approaches to the study of vegetable crops. This review examines the overall impact of major abiotic stresses on vegetables, adaptive mechanisms and functional genomic, transcriptomic, and proteomic processes used by researchers to minimize these challenges. The current status of genomics technologies for developing adaptable vegetable cultivars that will perform better in future climates is also examined.

High-quality genome assembly and genetic mapping reveal a gene regulating flesh color in watermelon (Citrullus lanatus)

The unique color and type characteristics of watermelon fruits are regulated by many molecular mechanisms. However, it still needs to be combined with more abundant genetic data to fine-tune the positioning. We assembled genomes of two Korean inbred watermelon lines (cv. 242-1 and 159-1) with unique color and fruit-type characteristics and identified 23,921 and 24,451 protein-coding genes in the two genomes, respectively. To obtain more precise results for further study, we resequenced one individual of each parental line and an F₂ population composed of 87 individuals. This identified 1,539 single-nucleotide polymorphisms (SNPs) and 80 InDel markers that provided a high-density genetic linkage map with a total length of 3,036.9 cM. Quantitative trait locus mapping identified 15 QTLs for watermelon fruit quality-related traits, including β-carotene and lycopene content in fruit flesh, fruit shape index, skin thickness, flesh color, and rind color. By investigating the mapping intervals, we identified 33 candidate genes containing variants in the coding sequence. Among them, Cla97C01G008760 was annotated as a phytoene synthase with a single-nucleotide variant (A → G) in the first exon at 9,539,129 bp of chromosome 1 that resulted in the conversion of a lysine to glutamic acid, indicating that this gene might regulate flesh color changes at the protein level. These findings not only prove the importance of a phytoene synthase gene in pigmentation but also explain an important reason for the color change of watermelon flesh.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

Primary mapping of quantitative trait loci regulating multivariate horticultural phenotypes of watermelon (Citrullus lanatus L.)

Watermelon fruits exhibit a remarkable diversity of important horticultural phenotypes. In this study, we initiated a primary quantitative trait loci (QTL) mapping to identify the candidate regions controlling the ovary, fruit, and seed phenotypes. Whole genome sequencing (WGS) was carried out for two differentiated watermelon lines, and 350 Mb (96%) and 354 Mb (97%) of re-sequenced reads covered the reference de novo genome assembly, individually. A total of 45.53% non-synonymous single nucleotide polymorphism (nsSNPs) and 54.47% synonymous SNPs (sSNPs) were spotted, which produced 210 sets of novel SNP-based cleaved amplified polymorphism sequence (CAPS) markers by depicting 46.25% co-dominant polymorphism among parent lines and offspring. A biparental F_2:3 mapping population comprised of 100 families was used for trait phenotyping and CAPS genotyping, respectively. The constructed genetic map spanned a total of 2,398.40 centimorgans (cM) in length and averaged 11.42 cM, with 95.99% genome collinearity. A total of 33 QTLs were identified at different genetic positions across the eight chromosomes of watermelon (Chr-01, Chr-02, Chr-04, Chr-05, Chr-06, Chr-07, Chr-10, and Chr-11); among them, eight QTLs of the ovary, sixteen QTLs of the fruit, and nine QTLs of the seed related phenotypes were classified with 5.32-25.99% phenotypic variance explained (PVE). However, twenty-four QTLs were identified as major-effect and nine QTLs were mapped as minor-effect QTLs across the flanking regions of CAPS markers. Some QTLs were exhibited as tightly localized across the nearby genetic regions and explained the pleiotropic effects of multigenic nature. The flanking QTL markers also depicted significant allele specific contributions and accountable genes were predicted for respective traits. Gene Ontology (GO) functional enrichment was categorized in molecular function (MF), cellular components (CC), and biological process (BP); however, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were classified into three main classes of metabolism, genetic information processing, and brite hierarchies. The principal component analysis (PCA) of multivariate phenotypes widely demonstrated the major variability, consistent with the identified QTL regions. In short, we assumed that our identified QTL regions provide valuable genetic insights regarding the watermelon phenotypes and fine genetic mapping could be used to confirm them.

Genetic mapping and molecular characterization of the delayed green gene dg in watermelon (Citrullus lanatus)

Leaf color mutants are common in higher plants that can be used as markers in crop breeding and are important tools in understanding regulatory mechanisms of chlorophyll biosynthesis and chloroplast development. Genetic analysis was performed by evaluating F₁, F₂ and BC₁ populations derived from two parental lines (Charleston gray with green leaf color and Houlv with delayed green leaf color), suggesting that a single recessive gene controls the delayed green leaf color. In this study, the delayed green mutant showed a conditional pale green leaf color at the early leaf development but turned to green as the leaf development progressed. Delayed green leaf plants showed reduced pigment content, photosynthetic, chlorophyll fluorescence parameters, and impaired chloroplast development compared with green leaf plants. The delayed green (dg) locus was mapped to 7.48 Mb on chromosome 3 through bulk segregant analysis approach, and the gene controlling delayed green leaf color was narrowed to 53.54 kb between SNP130 and SNP135 markers containing three candidate genes. Sequence alignment of the three genes indicated that there was a single SNP mutation (G/A) in the coding region of ClCG03G010030 in the Houlv parent, which causes an amino acid change from Arginine to Lysine. The ClCG03G010030 gene encoded FtsH extracellular protease protein family is involved in early delayed green leaf development. The expression level of ClCG03G010030 was significantly reduced in delayed green leaf plants than in green leaf plants. These results indicated that the ClCG03G010030 might control watermelon green leaf color and the single SNP variation in ClCG03G010030 may result in early delayed green leaf color development during evolutionary process.

The NF-Y Transcription Factor Family in Watermelon: Re-Characterization, Assembly of ClNF-Y Complexes, Hormone- and Pathogen-Inducible Expression and Putative Functions in Disease Resistance

Nuclear factor Y (NF-Y) is a heterotrimeric transcription factor that binds to the CCAAT cis-element in the promoters of target genes and plays critical roles in plant growth, development, and stress responses. In the present study, we aimed to re-characterize the ClNF-Y family in watermelon, examine the assembly of ClNF-Y complexes, and explore their possible involvement in disease resistance. A total of 25 ClNF-Y genes (7 ClNF-YAs, 10 ClNF-YBs, and 8 ClNF-YCs) were identified in the watermelon genome. The ClNF-Y family was comprehensively characterized in terms of gene and protein structures, phylogenetic relationships, and evolution events. Different types of cis-elements responsible for plant growth and development, phytohormones, and/or stress responses were identified in the promoters of the ClNF-Y genes. ClNF-YAs and ClNF-YCs were mainly localized in the nucleus, while most of the ClNF-YBs were localized in the cytoplasm of cells. ClNF-YB5, -YB6, -YB7, -YB8, -YB9, and -YB10 interacted with ClNF-YC2, -YC3, -YC4, -YC5, -YC6, -YC7, and -YC8, while ClNF-YB1 and -YB3 interacted with ClNF-YC1. A total of 37 putative ClNF-Y complexes were identified, e.g., ClNF-YA1, -YA2, -YA3, and -YA7 assembled into 13, 8, 8, and 8 ClNF-Y complexes with different ClNF-YB/-YC heterodimers. Most of the ClNF-Y genes responded with distinct expression patterns to defense hormones such as salicylic acid, methyl jasmonate, abscisic acid, and ethylene precursor 1-aminocyclopropane-1-carboxylate, and to infection by the vascular infecting fungus Fusarium oxysporum f. sp. niveum. Overexpression of ClNF-YB1, -YB8, -YB9, ClNF-YC2, and -YC7 in transgenic Arabidopsis resulted in an earlier flowering phenotype. Overexpression of ClNF-YB8 in Arabidopsis led to enhanced resistance while overexpression of ClNF-YA2 and -YC2 resulted in decreased resistance against Botrytis cinerea. Similarly, overexpression of ClNF-YA3, -YB1, and -YC4 strengthened resistance while overexpression of ClNF-YA2 and -YB8 attenuated resistance against Pseudomonas syringae pv. tomato DC3000. The re-characterization of the ClNF-Y family provides a basis from which to investigate the biological functions of ClNF-Y genes in respect of growth, development, and stress response in watermelon, and the identification of the functions of some ClNF-Y genes in disease resistance enables further exploration of the molecular mechanism of ClNF-Ys in the regulation of watermelon immunity against diverse pathogens.

Genotyping by Sequencing (GBS)-Based QTL Mapping for Bacterial Fruit Blotch (BFB) in Watermelon

Watermelon (Citrullus lanatus), an economically important and nutritionally rich Cucurbitaceous crop grown worldwide, is severely affected by bacterial fruit blotch (BFB). Development of resistant cultivar is the most eco-friendly, cost-effective, and sustainable way to tackle this disease. This requires wider understanding of the genetics of resistance to BFB. In this study, we identified quantitative trait loci (QTLs) associated with BFB resistance in an F₂ mapping population developed from BFB-resistant 'PI 189225' (Citrullus amarus) and -susceptible 'SW 26' (C. lanatus) genotypes based on the polymorphic markers identified by genotyping by sequencing (GSB). A linkage map covering a total genetic distance of 3377.1 cM was constructed. Two QTLs for BFB resistance, namely, ClBFB10.1 and ClBFB10.2, both located on chromosome 10 explaining 18.84 and 15.41% of the phenotypic variations, respectively, were identified. Two SNP-based high-resolution melting (HRM) markers WmBFB10.1 and WmBFB10.2 having high positive correlation with resistance vs. susceptible alleles were developed. The efficacy of the markers was validated in another F₂ population derived from SW34 × PI 189225. The highest phenotypic variation was found in the locus ClBFB10.2, which also contains three putative candidate genes for resistance to BFB. These findings will accelerate the development of BFB-resistant watermelon varieties via molecular breeding.

CuGenDBv2: an updated database for cucurbit genomics

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.

Diploid chromosome-level reference genome and population genomic analyses provide insights into Gypenoside biosynthesis and demographic evolution of Gynostemma pentaphyllum (Cucurbitaceae)

Gynostemma pentaphyllum (Thunb.) Makino is a perennial creeping herbaceous plant in the family Cucurbitaceae, which has great medicinal value and commercial potential, but urgent conservation efforts are needed due to the gradual decreases and fragmented distribution of its wild populations. Here, we report the high-quality diploid chromosome-level genome of G. pentaphyllum obtained using a combination of next-generation sequencing short reads, Nanopore long reads, and Hi-C sequencing technologies. The genome is anchored to 11 pseudo-chromosomes with a total size of 608.95 Mb and 26 588 predicted genes. Comparative genomic analyses indicate that G. pentaphyllum is estimated to have diverged from Momordica charantia 60.7 million years ago, with no recent whole-genome duplication event. Genomic population analyses based on genotyping-by-sequencing and ecological niche analyses indicated low genetic diversity but a strong population structure within the species, which could classify 32 G. pentaphyllum populations into three geographical groups shaped jointly by geographic and climate factors. Furthermore, comparative transcriptome analyses showed that the genes encoding enzyme involved in gypenoside biosynthesis had higher expression levels in the leaves and tendrils. Overall, the findings obtained in this study provide an effective molecular basis for further studies of demographic genetics, ecological adaption, and systematic evolution in Cucurbitaceae species, as well as contributing to molecular breeding, and the biosynthesis and biotransformation of gypenoside.

Mapping and functional verification of leaf yellowing genes in watermelon during whole growth period

Increasing light energy utilization efficiency is an effective way to increase yield and improve quality of watermelon. Leaf is the main place for photosynthesis, and the color of leaf is directly related to the change of photosynthesis. In addition, leaf yellowing can be used as a marker trait to play an important role in watermelon hybrid breeding and improve seed breeding. It can not only be used to eliminate hybrids at seedling stage, but also be used to determine seed purity. In this study, transcriptome analysis was first carried out using the whole growth period leaf yellowing watermelon mutant w-yl and inbred line ZK, and identified 2,471 differentially expressed genes (DEGs) in the comparison group w-yl-vs-ZK. Among the top 20 terms of the gene ontology (GO) enrichment pathway, 17 terms were related to photosynthesis. KEGG pathway enrichment analysis showed that the most abundant pathway was photosynthesis-antenna proteins. The F₂ population was constructed by conventional hybridization with the inbred line ZK. Genetic analysis showed that leaf yellowing of the mutant was controlled by a single recessive gene. The leaf yellowing gene of watermelon located between Ind14,179,011 and InD16,396,362 on chromosome 2 by using indel-specific PCR markers, with a region of 2.217 Mb. In the interval, it was found that five genes may have gene fragment deletion in w-yl, among which Cla97C02G036010, Cla97C02G036030, Cla97C02G036040, Cla97C02G036050 were the whole fragment loss, and Cla97C02G0360 was the C-terminal partial base loss. Gene function verification results showed that Cla97C02G036040, Cla97C02G036050 and Cla97C02G036060 may be the key factors leading to yellowing of w-yl leaves.

Identification of Candidate Genes for Rind Color and Bloom Formation in Watermelon Fruits Based on a Quantitative Trait Locus-Seq

Watermelon fruit rind color (RC) and bloom formation (BF) affect product value and consumer preference. However, information on the candidate gene(s) for additional loci involved in dark green (DG) RC and the genetic control of BF and its major chemical components is lacking. Therefore, this study aimed to identify loci controlling RC and BF using QTL-seq of the F₂ population derived by crossing 'FD061129' with light-green rind and bloom and 'SIT55616RN' with DG rind and bloomless. Phenotypic evaluation of the F₁ and 219 F₂ plants indicated the genetic control of two complementary dominant loci, G₁ and G₂, for DG and a dominant locus, Bf, for BF. QTL-seq identified a genomic region on Chr.6 for G₁, Chr.8 for G₂, and Chr.1 for Bf. G₁ and G₂ helped determine RC with possible environmental effects. Chlorophyll a-b binding protein gene-based CAPS (RC-m5) at G₁ matched the highest with the RC phenotype. In the 1.4 cM Bf map interval, two additional gene-based CAPS markers were designed, and the CAPS for a nonsynonymous SNP in Cla97C01G020050, encoding a CSC1-like protein, cosegregated with the BF trait in 219 F₂ plants. Bloom powder showed a high Ca²⁺ concentration (16,358 mg·kg^-1), indicating that the CSC1-like protein gene is possibly responsible for BF. Our findings provide valuable information for marker-assisted selection for RC and BF and insights into the functional characterization of genes governing these watermelon-fruit-related traits.

Genome-wide analysis of IQD proteins and ectopic expression of watermelon ClIQD24 in tomato suggests its important role in regulating fruit shape

The plant-specific IQ67 domain (IQD) is the largest class of calmodulin targets found in plants, and plays an important role in many biological processes, especially fruit development processes. However, the functional role of IQD proteins in the development of watermelon (Citrullus lanatus) shape remains unknown, as the IQD protein family in watermelon has not been systematically characterized. Herein, we elucidated the gene structures, chromosomal locations, evolutionary divergence, and functions of 35 IQD genes in the watermelon genome. The transcript profiles and quantitative real-time PCR analysis at different stages of fruit development showed that the ClIQD24 gene was highly expressed on 0 days after pollination. Furthermore, we found that the ectopic overexpression of ClIQD24 promoted tomato fruit elongation, thereby revealing the significance of ClIQD24 in the progression of watermelon shape. Our study will serve as a reference for further investigations on the molecular mechanisms underlying watermelon fruit shape formation.

A telomere-to-telomere gap-free reference genome of watermelon and its mutation library provide important resources for gene discovery and breeding

Watermelon, Citrullus lanatus, is the world's third largest fruit crop. Reference genomes with gaps and a narrow genetic base hinder functional genomics and genetic improvement of watermelon. Here, we report the assembly of a telomere-to-telomere gap-free genome of the elite watermelon inbred line G42 by incorporating high-coverage and accurate long-read sequencing data with multiple assembly strategies. All 11 chromosomes have been assembled into single-contig pseudomolecules without gaps, representing the highest completeness and assembly quality to date. The G42 reference genome is 369 321 829 bp in length and contains 24 205 predicted protein-coding genes, with all 22 telomeres and 11 centromeres characterized. Furthermore, we established a pollen-EMS mutagenesis protocol and obtained over 200 000 M1 seeds from G42 . In a sampling pool, 48 monogenic phenotypic mutations, selected from 223 M1 and 78 M2 mutants with morphological changes, were confirmed. The average mutation density was 1 SNP/1.69 Mb and 1 indel/4.55 Mb per M1 plant and 1 SNP/1.08 Mb and 1 indel/6.25 Mb per M2 plant. Taking advantage of the gap-free G42 genome, 8039 mutations from 32 plants sampled from M1 and M2 families were identified with 100% accuracy, whereas only 25% of the randomly selected mutations identified using the 97103v2 reference genome could be confirmed. Using this library and the gap-free genome, two genes responsible for elongated fruit shape and male sterility (ClMS1) were identified, both caused by a single base change from G to A. The validated gap-free genome and its EMS mutation library provide invaluable resources for functional genomics and genetic improvement of watermelon.

Genome sequencing of up to 6,000-yr-old Citrullus seeds reveals use of a bitter-fleshed species prior to watermelon domestication

Iconographic evidence from Egypt suggests that watermelon pulp was consumed there as a dessert by 4,360 BP. Earlier archaeobotanical evidence comes from seeds from Neolithic settlements in Libya, but whether these were watermelons with sweet pulp or other forms is unknown. We generated genome sequences from 6,000- and 3,300-year-old seeds from Libya and Sudan, and from worldwide herbarium collections made between 1824 and 2019, and analyzed these data together with resequenced genomes from important germplasm collections for a total of 131 accessions. Phylogenomic and population-genomic analyses reveal that (1) much of the nuclear genome of both ancient seeds is traceable to West African seed-use “egusi-type” watermelon (Citrullus mucosospermus) rather than domesticated pulp-use watermelon (Citrullus lanatus ssp. vulgaris); (2) the 6,000-year-old watermelon likely had bitter pulp and greenish-white flesh as today found in C. mucosospermus, given alleles in the bitterness regulators ClBT and in the red color marker LYCB; and (3) both ancient genomes showed admixture from C. mucosospermus, C. lanatus ssp. cordophanus, C. lanatus ssp. vulgaris, and even South African Citrullus amarus, and evident introgression between the Libyan seed (UMB-6) and populations of C. lanatus. An unexpected new insight is that Citrullus appears to have initially been collected or cultivated for its seeds, not its flesh, consistent with seed damage patterns induced by human teeth in the oldest Libyan material.

Genome-Wide Association Analysis of Resistance to Pseudoperonospora cubensis in Citron Watermelon

Cultivated sweet watermelon (Citrullus lanatus) is an important vegetable crop for millions of people around the world. There are limited sources of resistance to economically important diseases within C. lanatus, whereas C. amarus has a reservoir of traits that can be exploited to improve C. lanatus for resistance to biotic and abiotic stresses. Cucurbit downy mildew (CDM), caused by Pseudoperonospora cubensis, is an emerging threat to watermelon production. We screened 122 C. amarus accessions for resistance to CDM over two tests (environments). The accessions were genotyped by whole-genome resequencing to generate 2,126,759 single nucleotide polymorphic (SNP) markers. A genome-wide association study was deployed to uncover marker-trait associations and identify candidate genes underlying resistance to CDM. Our results indicate the presence of wide phenotypic variability (1.1 to 57.8%) for leaf area infection, representing a 50.7-fold variation for CDM resistance across the C. amarus germplasm collection. Broad-sense heritability estimate was 0.55, implying the presence of moderate genetic effects for resistance to CDM. The peak SNP markers associated with resistance to P. cubensis were located on chromosomes Ca03, Ca05, Ca07, and Ca11. The significant SNP markers accounted for up to 30% of the phenotypic variation and were associated with promising candidate genes encoding leucine-rich repeat receptor-like protein kinase and the WRKY transcription factor. This information will be useful in understanding the genetic architecture of the P. cubensis-Citrullus spp. patho-system as well as development of resources for genomics-assisted breeding for resistance to CDM in watermelon.

Multiple Stressors in Vegetable Production: Insights for Trait-Based Crop Improvement in Cucurbits

Vegetable production is a key determinant of contribution from the agricultural sector toward national Gross Domestic Product in a country like India, the second largest producer of fresh vegetables in the world. This calls for a careful scrutiny of the threats to vegetable farming in the event of climate extremes, environmental degradation and incidence of plant pests/diseases. Cucurbits are a vast group of vegetables grown almost throughout the world, which contribute to the daily diet on a global scale. Increasing food supply to cater to the ever-increasing world population, calls for intensive, off-season and year-round cultivation of cucurbits. Current situation predisposes these crops to a multitude of stressors, often simultaneously, under field conditions. This scenario warrants a systematic understanding of the different stress specific traits/mechanisms/pathways and their crosstalk that have been examined in cucurbits and identification of gaps and formulation of perspectives on prospective research directions. The careful dissection of plant responses under specific production environments will help in trait identification for genotype selection, germplasm screens to identify superior donors or for direct genetic manipulation by modern tools for crop improvement. Cucurbits exhibit a wide range of acclimatory responses to both biotic and abiotic stresses, among which a few like morphological characters like waxiness of cuticle; primary and secondary metabolic adjustments; membrane thermostability, osmoregulation and, protein and reactive oxygen species homeostasis and turnover contributing to cellular tolerance, appear to be common and involved in cross talk under combinatorial stress exposures. This is assumed to have profound influence in triggering system level acclimation responses that safeguard growth and metabolism. The possible strategies attempted such as grafting initiatives, molecular breeding, novel genetic manipulation avenues like gene editing and ameliorative stress mitigation approaches, have paved way to unravel the prospects for combined stress tolerance. The advent of next generation sequencing technologies and big data management of the omics output generated have added to the mettle of such emanated concepts and ideas. In this review, we attempt to compile the progress made in deciphering the biotic and abiotic stress responses of cucurbits and their associated traits, both individually and in combination.

Genetic Analysis of Fruit Quality Traits in Sweet Watermelon (Citrullus lanatus var. lanatus): A Review

Fruit quality traits of sweet watermelon (Citrullus lanatus var. lanatus) are crucial for new product development and commercialization. Sweet watermelon fruit quality traits are affected by the compositions of phytochemical compounds, phytohormones, and fruit flesh firmness which are affected by genes, the growing environment and their interaction. These compositions determine fruit ripening, eating quality, and postharvest shelf-life. Knowledge of the genetic profile and analyses of quality traits in watermelon is vital to develop improved cultivars with enhanced nutritional compositions, consumer-preferred traits, and extended storage life. This review aims to present the opportunities and progress made on the genetic analysis of fruit quality traits in watermelon as a guide for quality breeding based on economic and end-user attributes. The first section of the review highlights the genetic mechanisms involved in the biosynthesis of phytochemical compounds (i.e., sugars, carotenoids, amino acids, organic acids, and volatile compounds), phytohormones (i.e., ethylene and abscisic acid) and fruit flesh structural components (i.e., cellulose, hemicellulose, and pectin) elicited during watermelon fruit development and ripening. The second section pinpoints the progress on the development of molecular markers and quantitative trait loci (QTL) analysis for phytochemical compounds, phytohormones and fruit quality attributes. The review presents gene-editing technology and innovations associated with fruit quality traits for selection and accelerated cultivar development. Finally, the paper discussed gene actions conditioning fruit ripening in citron watermelon (C. lanatus var. citroides [L. H. Bailey] Mansf. ex Greb.) as reference genetic resources to guide current and future breeding. Information presented in this review is useful for watermelon variety design, product profiling and development to serve the diverse value chains of the crop.

Genome-wide SSR markers in bottle gourd: development, characterization, utilization in assessment of genetic diversity of National Genebank of India and synteny with other related cucurbits

Lagenaria siceraria (Molina) Standley is an important cultivated crop with its immense importance in pharmaceutical industry and as vegetable. Its seed, root, stem, leaves, flower, and fruit are used as an ointment for ailment of various diseases throughout Asia. Despite its worldwide importance, informative co-dominant microsatellite markers in the bottle gourd crop are very restricted, impeding genetic improvement, cultivar identification, and phylogenetic studies. Next-generation sequencing has revolutionized the approaches for discovery, assessment, and validation of molecular markers. We conducted a genome-wide analysis, for developing SSR markers by utilizing restriction site-associated DNA sequencing (RAD-Seq) data obtained from NCBI. By performing in silico mining of microsatellite repeat motifs, we developed 45,066 perfect SSR markers. Of which 207 markers were successfully validated and 120 (57.97%) polymorphic primer pairs were utilized for an in-depth genetic diversity and population structure analysis of 96 accessions from the National Genebank of India. Tetranucleotide repeats (∼34.3%) were the most prevalent followed by trinucleotide repeats (∼30.73%), further 21.03%, 9.6%, and 4.3% of di-, penta-, and hexa-nucleotide repeats in the bottle gourd genome, respectively. Synteny of SSR markers on 11 bottle gourd linkage groups was correlated with the 7 chromosomes of cucumber (93.2%), 12 chromosomes of melon (87.4%), and 11 of watermelon (90.8%). The generated SSR markers provide a valuable tool for germplasm characterization, genetic linkage map construction, studying synteny, gene discovery, and for breeding in bottle gourd and other cucurbits species. KEY MESSAGE: Development of 45,066 perfect microsatellite markers as a valuable tool for marker assisted selection (MAS) in plant breeding.

Cantaloupe melon genome reveals 3D chromatin features and structural relationship with the ancestral cucurbitaceae karyotype

Cucumis melo displays a large diversity of horticultural groups with cantaloupe melon the most cultivated type. Using a combination of single-molecule sequencing, 10X Genomics link-reads, high-density optical and genetic maps, and chromosome conformation capture (Hi-C), we assembled a chromosome scale C. melo var. cantalupensis Charentais mono genome. Integration of RNA-seq, MeDip-seq, ChIP-seq, and Hi-C data revealed a widespread compartmentalization of the melon genome, segregating constitutive heterochromatin and euchromatin. Genome-wide comparative and evolutionary analysis between melon botanical groups identified Charentais mono genome increasingly more divergent from Harukei-3 (reticulatus), Payzawat (inodorus), and HS (ssp. agrestis) genomes. To assess the paleohistory of the Cucurbitaceae, we reconstructed the ancestral Cucurbitaceae karyotype and compared it to sequenced cucurbit genomes. In contrast to other species that experienced massive chromosome shuffling, melon has retained the ancestral genome structure. We provide comprehensive genomic resources and new insights in the diversity of melon horticultural groups and evolution of cucurbits.

•

We provide a chromosome scale C. melo var. cantalupensis Charentais mono genome

•

Epigenomic analysis revealed a widespread compartmentalization of the melon genome

•

We reconstructed the ancestral Cucurbitaceae karyotype

•

Melon has retained the ancestral Cucurbitaceae genome structure

Cucurbitaceae genome evolution, gene function and molecular breeding

The 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 relationship, 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.

Genome-Wide Identification and Characterization of Melon bHLH Transcription Factors in Regulation of Fruit Development

The basic helix-loop-helix (bHLH) transcription factor family is one of the largest transcription factor families in plants and plays crucial roles in plant development. Melon is an important horticultural plant as well as an attractive model plant for studying fruit ripening. However, the bHLH gene family of melon has not yet been identified, and its functions in fruit growth and ripening are seldom researched. In this study, 118 bHLH genes were identified in the melon genome. These CmbHLH genes were unevenly distributed on chromosomes 1 to 12, and five CmbHLHs were tandem repeat on chromosomes 4 and 8. There were 13 intron distribution patterns among the CmbHLH genes. Phylogenetic analysis illustrated that these CmbHLHs could be classified into 16 subfamilies. Expression patterns of the CmbHLH genes were studied using transcriptome data. Tissue specific expression of the CmbHLH32 gene was analysed by quantitative RT-PCR. The results showed that the CmbHLH32 gene was highly expressed in female flower and early developmental stage fruit. Transgenic melon lines overexpressing CmbHLH32 were generated, and overexpression of CmbHLH32 resulted in early fruit ripening compared to wild type. The CmbHLH transcription factor family was identified and analysed for the first time in melon, and overexpression of CmbHLH32 affected the ripening time of melon fruit. These findings laid a foundation for further study on the role of bHLH family members in the growth and development of melon.

Rapid PCR-based method for herbivore dietary evaluation using plant-specific primers

Polyphagous pests cause significant economic loss worldwide through feeding damage on various cash crops. However, their diets in agricultural landscapes remain largely unexplored. Pest dietary evaluation in agricultural fields is a challenging task currently approached through visual observation of plant feeding and microscopic identification of semi-digested plant material in pest’s guts. While molecular gut content analysis using metabarcoding approaches using universal primers (e.g., rbcl and trnL) have been successful in evaluating polyphagous pest diet, this method is relatively costly and time-consuming. Hence, there is a need for a rapid, specific, sensitive, and cost-effective method to screen for crops in the gut of pests. This is the first study to develop plant-specific primers that target various regions of their genomes, designed using a whole plant genome sequence. We selected Verticillium wilt disease resistance protein (VE-1) and pathogenesis related protein-coding genes 1–5 (PR-1-5) as our targets and designed species-specific primers for 14 important crops in the agroecosystems. Using amplicon sizes ranging from 115 to 407 bp, we developed two multiplex primer mixes that can separate nine and five plant species per PCR reaction, respectively. These two designed primer mixes provide a rapid, sensitive and specific route for polyphagous pest dietary evaluation in agroecosystems. This work will enable future research to rapidly expand our knowledge on the diet preference and range of crops that pests consume in various agroecosystems, which will help in the redesign and development of new crop rotation regimes to minimize polyphagous pest pressure and damage on crops.

Deciphering the Genetic Architecture of Plant Virus Resistance by GWAS, State of the Art and Potential Advances

Growing virus resistant varieties is a highly effective means to avoid yield loss due to infection by many types of virus. The challenge is to be able to detect resistance donors within plant species diversity and then quickly introduce alleles conferring resistance into elite genetic backgrounds. Until now, mainly monogenic forms of resistance with major effects have been introduced in crops. Polygenic resistance is harder to map and introduce in susceptible genetic backgrounds, but it is likely more durable. Genome wide association studies (GWAS) offer an opportunity to accelerate mapping of both monogenic and polygenic resistance, but have seldom been implemented and described in the plant-virus interaction context. Yet, all of the 48 plant-virus GWAS published so far have successfully mapped QTLs involved in plant virus resistance. In this review, we analyzed general and specific GWAS issues regarding plant virus resistance. We have identified and described several key steps throughout the GWAS pipeline, from diversity panel assembly to GWAS result analyses. Based on the 48 published articles, we analyzed the impact of each key step on the GWAS power and showcase several GWAS methods tailored to all types of viruses.

PLAZA 5.0: extending the scope and power of comparative and functional genomics in plants

PLAZA is a platform for comparative, evolutionary, and functional plant genomics. It makes a broad set of genomes, data types and analysis tools available to researchers through a user-friendly website, an API, and bulk downloads. In this latest release of the PLAZA platform, we are integrating a record number of 134 high-quality plant genomes, split up over two instances: PLAZA Dicots 5.0 and PLAZA Monocots 5.0. This number of genomes corresponds with a massive expansion in the number of available species when compared to PLAZA 4.0, which offered access to 71 species, a 89% overall increase. The PLAZA 5.0 release contains information for 5 882 730 genes, and offers pre-computed gene families and phylogenetic trees for 5 274 684 protein-coding genes. This latest release also comes with a set of new and updated features: a new BED import functionality for the workbench, improved interactive visualizations for functional enrichments and genome-wide mapping of gene sets, and a fully redesigned and extended API. Taken together, this new version offers extended support for plant biologists working on different families within the green plant lineage and provides an efficient and versatile toolbox for plant genomics. All PLAZA releases are accessible from the portal website: https://bioinformatics.psb.ugent.be/plaza/.

Genome-wide association study and population structure analysis of seed-bound amino acids and total protein in watermelon

BACKGROUND

Watermelon seeds are a powerhouse of value-added traits such as proteins, free amino acids, vitamins, and essential minerals, offering a paleo-friendly dietary option. Despite the availability of substantial genetic variation, there is no sufficient information on the natural variation in seed-bound amino acids or proteins across the watermelon germplasm. This study aimed to analyze the natural variation in watermelon seed amino acids and total protein and explore underpinning genetic loci by genome-wide association study (GWAS).

METHODS

The study evaluated the distribution of seed-bound free amino acids and total protein in 211 watermelon accessions of Citrullus spp, including 154 of Citrullus lanatus, 54 of Citrullus mucosospermus (egusi) and three of Citrullus amarus. We used the GWAS approach to associate seed phenotypes with 11,456 single nucleotide polymorphisms (SNPs) generated by genotyping-by-sequencing (GBS).

RESULTS

Our results demonstrate a significant natural variation in different free amino acids and total protein content across accessions and geographic regions. The accessions with high protein content and proportion of essential amino acids warrant its use for value-added benefits in the food and feed industries via biofortification. The GWAS analysis identified 188 SNPs coinciding with 167 candidate genes associated with watermelon seed-bound amino acids and total protein. Clustering of SNPs associated with individual amino acids found by principal component analysis was independent of the speciation or cultivar groups and was not selected during the domestication of sweet watermelon. The identified candidate genes were involved in metabolic pathways associated with amino acid metabolism, such as Argininosuccinate synthase, explaining 7% of the variation in arginine content, which validate their functional relevance and potential for marker-assisted analysis selection. This study provides a platform for exploring potential gene loci involved in seed-bound amino acids metabolism, useful in genetic analysis and development of watermelon varieties with superior seed nutritional values.

Candidate gene analysis of watermelon stripe pattern locus ClSP ongoing recombination suppression

Using two segregating population, watermelon stripe pattern underlying gene ClSP was delimited to a 611.78 Kb region, consisting of four discrete haploblocks and ongoing recombination suppression. Stripe pattern is an important commodity trait in watermelon, displaying diverse types. In this study, two segregating populations were generated for genetic mapping the single dominant locus ClSP, which was finally delimited to a 611.78 Kb interval with suppression of recombination. According to polymorphism sites detected among genotypes, four discrete haploblocks were characterized in this target region. Based on reference genomes, 81 predicted genes were annotated in the ClSP interval, including seven transcription factors namely as candidate No1-No7. Meanwhile, the ortholog gene of cucumber ist responsible for the irregular stripes was considered as candidate No8. Strikingly, gene structures of No1-No5 completely varied from their reference descriptions and subsequently re-annotated. For instance, the original adjacent distribution candidates No2 and No3 were re-annotated as No2_3, while No4 and No5 were integrated as No4_5. Sequence analysis demonstrated the third polymorphism in CDS of re-annotated No4_5 resulting in truncated proteins in non-stripe plants. Furthermore, only No4_5 was down-regulated in light green stripes relative to dark green stripes. Transcriptome analysis identified 356 DEGs between dark green striped and light green striped peels, with genes involved in photosynthesis and chloroplast development down-regulated in light green stripes but calcium ion binding related genes up-regulated. Additionally, 38 DEGs were annotated as transcription factors, with the majority up-regulated in light green stripes, such as ERFs and WRKYs. This study not only contributes to a better understanding of the molecular mechanisms underlying watermelon stripe development, but also provides new insights into the genomic structure of ClSP locus and valuable candidates.

A Catalog of Natural Products Occurring in Watermelon-Citrullus lanatus

Sweet dessert watermelon (Citrullus lanatus) is one of the most important vegetable crops consumed throughout the world. The chemical composition of watermelon provides both high nutritional value and various health benefits. The present manuscript introduces a catalog of 1,679 small molecules occurring in the watermelon and their cheminformatics analysis for diverse features. In this catalog, the phytochemicals are associated with the literature describing their presence in the watermelon plant, and when possible, concentration values in various plant parts (flesh, seeds, leaves, roots, rind). Also cataloged are the chemical classes, molecular weight and formula, chemical structure, and certain physical and chemical properties for each phytochemical. In our view, knowing precisely what is in what we eat, as this catalog does for watermelon, supports both the rationale for certain controlled feeding studies in the field of precision nutrition, and plant breeding efforts for the development of new varieties with enhanced concentrations of specific phytochemicals. Additionally, improved and comprehensive collections of natural products accessible to the public will be especially useful to researchers in nutrition, cheminformatics, bioinformatics, and drug development, among other disciplines.

Genetic Resources and Vulnerabilities of Major Cucurbit Crops

The Cucurbitaceae family provides numerous important crops including watermelons (Citrullus lanatus), melons (Cucumis melo), cucumbers (Cucumis sativus), and pumpkins and squashes (Cucurbita spp.). Centers of domestication in Africa, Asia, and the Americas were followed by distribution throughout the world and the evolution of secondary centers of diversity. Each of these crops is challenged by multiple fungal, oomycete, bacterial, and viral diseases and insects that vector disease and cause feeding damage. Cultivated varieties are constrained by market demands, the necessity for climatic adaptations, domestication bottlenecks, and in most cases, limited capacity for interspecific hybridization, creating narrow genetic bases for crop improvement. This analysis of crop vulnerabilities examines the four major cucurbit crops, their uses, challenges, and genetic resources. ex situ germplasm banks, the primary strategy to preserve genetic diversity, have been extensively utilized by cucurbit breeders, especially for resistances to biotic and abiotic stresses. Recent genomic efforts have documented genetic diversity, population structure, and genetic relationships among accessions within collections. Collection size and accessibility are impacted by historical collections, current ability to collect, and ability to store and maintain collections. The biology of cucurbits, with insect-pollinated, outcrossing plants, and large, spreading vines, pose additional challenges for regeneration and maintenance. Our ability to address ongoing and future cucurbit crop vulnerabilities will require a combination of investment, agricultural, and conservation policies, and technological advances to facilitate collection, preservation, and access to critical Cucurbitaceae diversity.

Incidence of Alternaria Species Associated with Watermelon Leaf Blight in Korea

Alternaria leaf blight is one of the most common diseases in watermelon worldwide. In Korea, however, the Alternaria species causing the watermelon leaf blight have not been investigated thoroughly. A total of 16 Alternaria isolates was recovered from diseased watermelon leaves with leaf blight symptoms, which were collected from 14 fields in Korea. Analysis of internal transcribed spacer (ITS) region, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and RNA polymerase II second largest subunit (RPB2) were not competent to differentiate the Alternaria isolates. On the contrary, analysis of amplicon size of the histone H3 (HIS3) gene successfully differentiated the isolates into three Alternaria subgroups, and further sequence analysis of them identified three Alternaria spp. Alternaria tenuissima, A. gaisen, and A. alternata. Representative Alternaria isolates from three species induced dark brown leaf spot lesions on detached watermelon leaves, indicating that A. tenuissima, A. gaisen, and A. alternata are all causal agents of Alternaria leaf blight. Our results indicate that the Alternaria species associated watermelon leaf blight in Korea is more complex than reported previously. This is the first report regarding the population structure of Alternaria species causing watermelon leaf blight in Korea.

The Transcriptomic Profile of Watermelon Is Affected by Zinc in the Presence of Fusarium oxysporum f. sp. niveum and Meloidogyne incognita

Zinc (Zn) accumulation and deficiency affect plant response to pests and diseases differently in varying pathosystems. The concentrations of Zn in plants aid in priming defense signaling pathways and help in enhanced structural defenses against plant pathogens. Studies are lacking on how concentrations of Zn in watermelon plants influence defense against two important soil-borne pathogens: Fusarium oxysporum f. sp. niveum (FON) and southern root-knot nematode (RKN, Meloidogyne incognita). In this study a comparative transcriptomics evaluation of watermelon plants in response to high (1.2 ppm) and low (0.2 ppm) levels of Zn were determined. Differential transcript-level responses differed in watermelon plants when infected with FON or RKN or both under high- and low-Zn treatment regimes in a controlled hydroponics system. Higher numbers of differentially expressed genes (DEGs) were observed in high-Zn-treated than in low-Zn-treated non-inoculated plants, in plants inoculated with FON alone and in plants inoculated with RKN alone. However, in the co-inoculated system, low-Zn treatment had higher DEGs as compared to high-Zn treatment. In addition, most DEGs were significantly enriched in hormone signal transduction and MAPK signaling pathway, suggesting an induction of systemic resistance with high-Zn concentrations. Taken together, this study substantially expands transcriptome data resources and suggests a molecular potential framework for watermelon-Zn interaction in FON and RKN.

A tobacco ringspot virus-based vector system for gene and microRNA function studies in cucurbits

Cucurbits are economically important crops worldwide. The genomic data of many cucurbits are now available. However, functional analyses of cucurbit genes and noncoding RNAs have been impeded because genetic transformation is difficult for many cucurbitaceous plants. Here, we developed a set of tobacco ringspot virus (TRSV)-based vectors for gene and microRNA (miRNA) function studies in cucurbits. A TRSV-based expression vector could simultaneously express GREEN FLUORESCENT PROTEIN (GFP) and heterologous viral suppressors of RNA silencing in TRSV-infected plants, while a TRSV-based gene silencing vector could knock down endogenous genes exemplified by PHYTOENE DESATURASE (PDS) in Cucumis melo, Citrullus lanatus, Cucumis sativus, and Nicotiana benthamiana plants. We also developed a TRSV-based miRNA silencing vector to dissect the functions of endogenous miRNAs. Four representative miRNAs, namely, miR159, miR166, miR172, and miR319, from different cucurbits were inserted into the TRSV vector using a short tandem target mimic strategy and induced characteristic phenotypes in TRSV-miRNA-infected plants. This TRSV-based vector system will facilitate functional genomic studies in cucurbits.

A chromosome-level genome of a Kordofan melon illuminates the origin of domesticated watermelons

Wild relatives or progenitors of crops are important resources for breeding and for understanding domestication. Identifying them, however, is difficult because of extinction, hybridization, and the challenge of distinguishing them from feral forms. Here, we use collection-based systematics, iconography, and resequenced accessions of Citrullus lanatus and other species of Citrullus to search for the potential progenitor of the domesticated watermelon. A Sudanese form with nonbitter whitish pulp, known as the Kordofan melon (C. lanatus subsp. cordophanus), appears to be the closest relative of domesticated watermelons and a possible progenitor, consistent with newly interpreted Egyptian tomb paintings that suggest that the watermelon may have been consumed in the Nile Valley as a dessert by 4360 BP. To gain insights into the genetic changes that occurred from the progenitor to the domesticated watermelon, we assembled and annotated the genome of a Kordofan melon at the chromosome level, using a combination of Pacific Biosciences and Illumina sequencing as well as Hi-C mapping technologies. The genetic signature of bitterness loss is present in the Kordofan melon genome, but the red fruit flesh color only became fixed in the domesticated watermelon. We detected 15,824 genome structural variants (SVs) between the Kordofan melon and a typical modern cultivar, "97103," and mapping the SVs in over 400 Citrullus accessions revealed shifts in allelic frequencies, suggesting that fruit sweetness has gradually increased over the course of watermelon domestication. That a likely progenitor of the watermelon still exists in Sudan has implications for targeted modern breeding efforts.

Altered chromatin conformation and transcriptional regulation in watermelon following genome doubling

Polyploidy has played a crucial role in plant evolution, development and function. Synthetic autopolyploid represents an ideal system to investigate the effects of polyploidization on transcriptional regulation. In this study, we deciphered the impact of genome duplication at phenotypic and molecular levels in watermelon. Overall, 88% of the genes in tetraploid watermelon followed a >1:1 dosage effect, and accordingly, differentially expressed genes were largely upregulated. In addition, a great number of hypomethylated regions (1688) were identified in an isogenic tetraploid watermelon. These differentially methylated regions were localized in promoters and intergenic regions and near transcriptional start sites of the identified upregulated genes, which enhances the importance of methylation in gene regulation. These changes were reflected in sophisticated higher-order chromatin structures. The genome doubling caused switching of 108 A and 626 B compartments that harbored genes associated with growth, development and stress responses.

QTL associated with gummy stem blight resistance in watermelon

We identified QTLs associated with gummy stem blight resistance in an interspecific F_2:3 Citrullus population and developed marker assays for selection of the loci in watermelon. Gummy stem blight (GSB), caused by three Stagonosporopsis spp., is a devastating fungal disease of watermelon (Citrullus lanatus) and other cucurbits that can lead to severe yield losses. Currently, no commercial cultivars with genetic resistance to GSB in the field have been reported. Utilizing GSB-resistant cultivars would reduce yield losses, decrease the high cost of disease control, and diminish hazards resulting from frequent fungicide application. The objective of this study was to identify quantitative trait loci (QTLs) associated with GSB resistance in an F_2:3 interspecific Citrullus mapping population (N = 178), derived from a cross between Crimson Sweet (C. lanatus) and GSB-resistant PI 482276 (C. amarus). The population was phenotyped by inoculating seedlings with Stagonosporopsis citrulli 12178A in the greenhouse in two separate experiments, each with three replications. We identified three QTLs (ClGSB3.1, ClGSB5.1 and ClGSB7.1) associated with GSB resistance, explaining between 6.4 and 21.1% of the phenotypic variation. The genes underlying ClGSB5.1 includes an NBS-LRR gene (ClCG05G019540) previously identified as a candidate gene for GSB resistance in watermelon. Locus ClGSB7.1 accounted for the highest phenotypic variation and harbors twenty-two candidate genes associated with disease resistance. Among them is ClCG07G013230, encoding an Avr9/Cf-9 rapidly elicited disease resistance protein, which contains a non-synonymous point mutation in the DUF761 domain that was significantly associated with GSB resistance. High throughput markers were developed for selection of ClGSB5.1 and ClGSB7.1. Our findings will facilitate the use of molecular markers for efficient introgression of the resistance loci and development of GSB-resistant watermelon cultivars.

Combined genomic, transcriptomic, and metabolomic analyses provide insights into chayote (Sechium edule) evolution and fruit development

Chayote (Sechium edule) is an agricultural crop in the Cucurbitaceae family that is rich in bioactive components. To enhance genetic research on chayote, we used Nanopore third-generation sequencing combined with Hi-C data to assemble a draft chayote genome. A chromosome-level assembly anchored on 14 chromosomes (N50 contig and scaffold sizes of 8.40 and 46.56 Mb, respectively) estimated the genome size as 606.42 Mb, which is large for the Cucurbitaceae, with 65.94% (401.08 Mb) of the genome comprising repetitive sequences; 28,237 protein-coding genes were predicted. Comparative genome analysis indicated that chayote and snake gourd diverged from sponge gourd and that a whole-genome duplication (WGD) event occurred in chayote at 25 ± 4 Mya. Transcriptional and metabolic analysis revealed genes involved in fruit texture, pigment, flavor, flavonoids, antioxidants, and plant hormones during chayote fruit development. The analysis of the genome, transcriptome, and metabolome provides insights into chayote evolution and lays the groundwork for future research on fruit and tuber development and genetic improvements in chayote.

The genome and transcriptome analysis of snake gourd provide insights into its evolution and fruit development and ripening

Snake gourd (Trichosanthes anguina L.), which belongs to the Cucurbitaceae family, is a popular ornamental and food crop species with medicinal value and is grown in many parts of the world. Although progress has been made in its genetic improvement, the organization, composition, and evolution of the snake gourd genome remain largely unknown. Here, we report a high-quality genome assembly for snake gourd, comprising 202 contigs, with a total size of 919.8 Mb and an N50 size of 20.1 Mb. These findings indicate that snake gourd has one of the largest genomes of Cucurbitaceae species sequenced to date. The snake gourd genome assembly harbors 22,874 protein-coding genes and 80.0% of the genome consists of repetitive sequences. Phylogenetic analysis reveals that snake gourd is closely related to sponge gourd but diverged from their common ancestor ~33-47 million years ago. The genome sequence reported here serves as a valuable resource for snake gourd genetic research and comparative genomic studies in Cucurbitaceae and other plant species. In addition, fruit transcriptome analysis reveals the candidate genes related to quality traits during snake gourd fruit development and provides a basis for future research on snake gourd fruit development and ripening at the transcript level.

Elucidation of resistance signaling and identification of powdery mildew resistant mapping loci (ClaPMR2) during watermelon-Podosphaera xanthii interaction using RNA-Seq and whole-genome resequencing approach

Watermelon is an important vegetable crop and is widely cultivated in USA with an approximate global production of > 100 million tons. Powdery mildew (PM) caused by Podosphaera xanthii is a major production-limiting factor on watermelon and other cucurbits. Numerous PM and multiple disease resistant (MDR) watermelon germplasm lines have been developed by the USDA in Charleston, SC. To gain a better understanding of the innate and activated molecular defense mechanisms involved during compatible and incompatible PM-watermelon interactions, we inoculated PM susceptible (USVL677-PMS) and resistant (USVL531-MDR) watermelon plants with 10⁵ conidia ml⁻¹ of P. xanthii. RNA-seq profiling was done on leaf samples collected at 0, 1, 3, and 8 days post inoculation (DPI). A total of 2,566 unique differentially expressed genes (DEGs) were identified between compatible and incompatible interactions with P. xanthii. The compatible interactions resulted in distinct plant gene activation (> twofold unique transcripts, 335:191:1762 :: 1:3:8 DPI) as compared to incompatible interaction (> twofold unique transcripts, 314:681:487 :: 1:3:8 DPI). Further, comparative whole-genome resequencing analysis of USVL531-PMR, USVL677-PMS and four introgressed PM resistant recombinant inbred lines (RIL, USVL531-PMR × USVL677-PMS) were performed to identify the region of PM resistance introgressed break points along with other traits inherent by USVL531-PMR by comparing the SNPs and InDels. Based on SNPs identification and CAPS markers, the resistance gene was identified as ClaPMR2, Citrullus lanatus PM Resistance gene 2 {Chr2 : 26750001 .. 26753327 (−)}, a NBS-LRR resistance protein (R) with homology to the Arabidopsis thaliana PM resistance protein, RPW8. The transcriptome data also revealed a complex regulatory network associated with the introgressed junctions mediated by PM resistance R proteins (R genes) that may involve multiple signal regulators and transducers, carbohydrate metabolism, cell wall modifications and the hormone-signaling pathway.

Comparative genomics of muskmelon reveals a potential role for retrotransposons in the modification of gene expression

Melon exhibits substantial natural variation especially in fruit ripening physiology, including both climacteric (ethylene-producing) and non-climacteric types. However, genomic mechanisms underlying such variation are not yet fully understood. Here, we report an Oxford Nanopore-based high-grade genome reference in the semi-climacteric cultivar Harukei-3 (378 Mb + 33,829 protein-coding genes), with an update of tissue-wide RNA-seq atlas in the Melonet-DB database. Comparison between Harukei-3 and DHL92, the first published melon genome, enabled identification of 24,758 one-to-one orthologue gene pairs, whereas others were candidates of copy number variation or presence/absence polymorphisms (PAPs). Further comparison based on 10 melon genome assemblies identified genome-wide PAPs of 415 retrotransposon Gag-like sequences. Of these, 160 showed fruit ripening-inducible expression, with 59.4% of the neighboring genes showing similar expression patterns (r > 0.8). Our results suggest that retrotransposons contributed to the modification of gene expression during diversification of melon genomes, and may affect fruit ripening-inducible gene expression.