Resequencing of 414 cultivated and wild watermelon accessions identifies selection for fruit quality traits

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.

Watermelon domestication was shaped by stepwise selection and regulation of the metabolome

Although crop domestication has greatly aided human civilization, the sequential domestication and regulation of most quality traits remain poorly understood. Here, we report the stepwise selection and regulation of major fruit quality traits that occurred during watermelon evolution. The levels of fruit cucurbitacins and flavonoids were negatively selected during speciation, whereas sugar and carotenoid contents were positively selected during domestication. Interestingly, fruit malic acid and citric acid showed the opposite selection trends during the improvement. We identified a novel gene cluster (CGC1, cucurbitacin gene cluster on chromosome 1) containing both regulatory and structural genes involved in cucurbitacin biosynthesis, which revealed a cascade of transcriptional regulation operating mechanisms. In the CGC1, an allele caused a single nucleotide change in ClERF1 binding sites (GCC-box) in the promoter of ClBh1, which resulted in reduced expression of ClBh1 and inhibition of cucurbitacin synthesis in cultivated watermelon. Functional analysis revealed that a rare insertion of 244 amino acids, which arose in C. amarus and became fixed in sweet watermelon, in ClOSC (oxidosqualene cyclase) was critical for the negative selection of cucurbitacins during watermelon evolution. This research provides an important resource for metabolomics-assisted breeding in watermelon and for exploring metabolic pathway regulation mechanisms.

Genetic mapping of a single nuclear locus determines the white flesh color in watermelon (Citrullus lanatus L.)

INTRODUCTION

Flesh color is an important trait in watermelon (Citrullus lanatus L.). Several flesh color genes have been identified in watermelon; however, the inheritance of and the molecular basis underlying the white flesh trait remain largely unknown.

METHODS

In this study, segregation populations were constructed by crossing the canary yellow flesh line HSH-F with the white flesh line Sanbai to fine-map the white flesh gene in watermelon.

RESULTS

Genetic analysis indicated that the white flesh trait is controlled by a single recessive locus, termed Clwf2. Map-based cloning delimited the Clwf2 locus to a 132.3-kb region on chromosome 6. The candidate region contains 13 putative genes, and four of them-Cla97C06G121860, Cla97C06G121880, Cla97C06G121890, and Cla97C06G121900-were significantly downregulated in the white flesh compared to the canary yellow flesh watermelon fruits. The Cla97C06G121890 gene, which encodes a tetratricopeptide repeat protein, showed almost no expression in the white flesh fruit before maturity, whereas it had a very high expression in the canary yellow flesh fruit at 18 days after pollination. Transmission electron microscopy revealed rounded and regularly shaped chromoplasts in both the canary yellow and white flesh fruits. Further quantitative real-time PCR analysis showed that the expression levels of several key plastid division genes and almost the entire carotenoid biosynthesis pathway genes were downregulated in the white flesh compared to the canary yellow flesh fruits.

DISCUSSION

This study suggests that the proliferation inhibition of chromoplasts and downregulation of the CBP genes block the accumulation of carotenoids in watermelon and lead to white flesh. These findings advance and extend the understanding of the molecular mechanisms underlying white flesh trait formation and carotenoid biosynthesis in watermelon.

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.

Identification of the egusi seed trait locus (eg) and its suppressor gene associated with the thin seed coat trait in watermelon

Egusi watermelon has a unique egusi seed type, which could be useful for breeding both edible seeds and edible flesh in watermelon. However, the genetic basis of the unique egusi seed type is not clear. In the present study, we first reported that at least two genes with inhibitory epistasis were responsible for the thin seed coat (unique egusi seed type) in watermelon. Inheritance analysis of five populations, including F₂, BC, and BCF₂, suggested that the thin seed coat trait was controlled by a suppressor gene together with the egusi seed locus (eg) in egusi watermelon. Based on high-throughput sequencing technology, two quantitative trait loci located on chromosome 1 and chromosome 6 were identified for the thin seed coat trait in watermelon. One of the loci, the eg locus on chromosome 6, was finely mapped to a genomic region of 15.7 kb, which contained only one candidate gene. Comparative transcriptome analysis highlighted differentially expressed genes involved in cellulose and lignin synthesis between watermelon genotypes varying in the thickness of the seed coat and provided several potential candidate genes for the thin seed coat trait. Taken together, our data suggest that at least two genes are complementarily involved in the thin seed coat trait and will be useful for cloning novel genes. The results presented here provide a new reference for uncovering egusi seed genetic mechanisms and valuable information for marker-assisted selection in seed coat breeding.

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.

Comparative network analysis reveals the dynamics of organic acid diversity during fruit ripening in peach (Prunus persica L. Batsch)

BACKGROUND

Organic acids are important components that determine the fruit flavor of peach (Prunus persica L. Batsch). However, the dynamics of organic acid diversity during fruit ripening and the key genes that modulate the organic acids metabolism remain largely unknown in this kind of fruit tree which yield ranks sixth in the world.

RESULTS

In this study, we used 3D transcriptome data containing three dimensions of information, namely time, phenotype and gene expression, from 5 different varieties of peach to construct gene co-expression networks throughout fruit ripening of peach. With the network inferred, the time-ordered network comparative analysis was performed to select high-acid specific gene co-expression network and then clarify the regulatory factors controlling organic acid accumulation. As a result, network modules related to organic acid synthesis and metabolism under high-acid and low-acid comparison conditions were identified for our following research. In addition, we obtained 20 candidate genes as regulatory factors related to organic acid metabolism in peach.

CONCLUSIONS

The study provides new insights into the dynamics of organic acid accumulation during fruit ripening, complements the results of classical co-expression network analysis and establishes a foundation for key genes discovery from time-series multiple species transcriptome data.

The coiled-coil protein gene WPRb confers recessive resistance to Cucumber green mottle mosaic virus

Cucumber green mottle mosaic virus (CGMMV) is one of the major global quarantine viruses and causes severe symptoms in Cucurbit crops, particularly with regard to fruit decay. However, the genetic mechanisms that control plant resistance to CGMMV have yet to be elucidated. Here, we found that WPRb, a weak chloroplast movement under blue light 1 and plastid movement impaired 2-related protein family gene, is recessively associated with CGMMV resistance in watermelon (Citrullus lanatus). We developed a reproducible marker based on a single non-synonymous substitution (G1282A) in WPRb, which can be used for marker-assisted selection for CGMMV resistance in watermelon. Editing of WPRb conferred greater tolerance to CGMMV. We found WPRb targets to the plasmodesmata (PD) and biochemically interacts with the CGMMV movement protein, facilitating viral intercellular movement by affecting the permeability of PD. Our findings enable us to genetically control CGMMV resistance in planta by using precise genome editing techniques targeted to WPRb.

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.

Telomere-to-telomere genome assembly of bitter melon (Momordica charantia L. var. abbreviata Ser.) reveals fruit development, composition and ripening genetic characteristics

Momordica charantia L. var. abbreviata Ser. (Mca), known as bitter gourd or bitter melon, is a Momordica variety with medicinal value and belongs to the Cucurbitaceae family. In view of the lack of genomic information on bitter gourd and other Momordica species and to promote Mca genomic research, we assembled a 295.6-Mb telomere-to-telomere (T2T) high-quality Mca genome with six gap-free chromosomes after Hi-C correction. This genome is anchored to 11 chromosomes, which is consistent with the karyotype information, and comprises 98 contigs (N50 of 25.4 Mb) and 95 scaffolds (N50 of 25.4 Mb). The Mca genome harbors 19 895 protein-coding genes, of which 45.59% constitute predicted repeat sequences. Synteny analysis revealed variations involved in fruit quality during the divergence of bitter gourd. In addition, assay for transposase-accessible chromatin by high-throughput sequencing and metabolic analysis showed that momordicosides and other substances are characteristic of Mca fruit pulp. A combined transcriptomic and metabolomic analysis revealed the mechanisms of pigment accumulation and cucurbitacin biosynthesis in Mca fruit peels, providing fundamental molecular information for further research on Mca fruit ripening. This report provides a new genetic resource for Momordica genomic studies and contributes additional insights into Cucurbitaceae phylogeny.

Genome-wide association study provides genetic insights into natural variation in watermelon rind thickness and single fruit weight

Rind thickness and fruit weight are agronomic traits closely related to quality and yield, which have attracted much attention from consumers and breeders. However, the genetic mechanism of these two traits is still not well understood in natural populations. In this study, rind thickness and single fruit weight in 151 watermelon accessions were determined in 2019 and 2020, and genome-wide association analysis was performed by integrating phenotypic and genotype data. Abundant phenotypic variation was found in the test population, and the watermelon with thinner rind thickness tended to have smaller fruit weights. Five significant SNPs were closely associated with rind thickness on chromosome 2 by Genome-wide association study (GWAS), i.e., 32344170, 32321308, 32304738, 32328501, and 32311192. And there were 21 genes were annotated in the candidate interval, most notably, Cla97C02G044160 belonged to the MADS family, and part of the genes in this family played an important role in regulating organ size, further analysis of gene structure, gene expression level, and phylogenetic tree showed that Cla97C02G044160 was a candidate gene for regulating target traits. In conclusion, our study provides molecular insights into the natural variation of watermelon rind thickness and single fruit weight, meanwhile, providing data support for molecular marker-assisted breeding.

Chromosomal fragment deletion in APRR2-repeated locus modulates the dark stem color in Cucurbita pepo

Cp4.1LG15g03420 (CpDsc-1), which encodes a two-component response regulator-like protein (APRR2) in the nucleus, influences dark green stem formation in Cucurbita pepo by regulating the chlorophyll content. Stem color is an important agronomic trait in zucchini (Cucurbita pepo) for robust seeding and high yield. However, the gene controlling the stem color has not been characterized. In this study, we identified a single locus accounting for the dark green stem color of C. pepo (CpDsc-1). Genetic analysis of this trait in segregated populations derived from two parental lines (line 296 with dark green stems and line 274 with light green stems) revealed that stem color was controlled by a single dominant gene (dark green vs. light green). In bulked segregant analysis, CpDsc-1 was mapped to a 2.09-Mb interval on chromosome 15. This region was further narrowed to 65.2 kb using linkage analysis of the F₂ population. Sequencing analysis revealed a 14 kb deletion between Cp4.1LG15g03420 and Cp4.1LG15g03360; these two genes both encoded a two-component response regulator-like protein (APRR2). The incomplete structures of the two APRR2 genes and abnormal chloroplasts in line 274 might be the main cause of the light green phenotype. Gene expression pattern analysis showed that only Cp4.1LG15g03420 was upregulated in line 296. Subcellular localization analysis indicated that Cp4.1LG15g03420 was a nuclear gene. Furthermore, a co-dominant marker, G4563 (93% accuracy rate), and a co-segregation marker, Fra3, were established in 111 diverse germplasms; both of these markers were tightly linked with the color trait. This study provided insights into chlorophyll regulation mechanisms and revealed the markers valuable for marker-assisted selection in future zucchini breeding.

Architecture design of cucurbit crops for enhanced productivity by a natural allele

Increasing production efficiency is a top priority in agriculture. Optimal plant architecture is the biological basis of dense planting, high crop yield and labour cost savings, and is thus critical for improving agricultural productivity. In cucurbit crops, most species have elongated internodes, but the path to architecture improvement is still not clear. Here we identified a pumpkin accession with a dominant bushy trait, and found that the associated Bush locus harbours a cucurbit-conserved cis-regulatory element in the 5' untranslated region of a transcription factor gene YABBY1. In cucurbit crops, various B-region deletions enhance the translation of YABBY1, with consequent proportional suppression of stem length in a dose-dependent manner. Depending on different cultivation patterns, the precise deployment of these alleles has significant effects on yield improvement or labour cost saving. Our findings demonstrate that the engineering of the YABBY1 B-region is an efficient strategy to customize plant architecture in cucurbit crops.

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.

QTL associated with resistance to Stagonosporopsis citrulli in Citrullus amarus

Gummy stem blight (GSB) is a fungal disease affecting cucurbit crops, including watermelon (Citrullus lanatus), leading to significant yield losses. The disease is caused by three Stagonosporopsis species, of which Stagonosporopsis citrulli is the most common in the southeastern United States. Currently no gummy stem blight-resistant watermelon cultivars are available to growers. In this study, QTL-seq in an interspecific population developed from Sugar Baby × PI 189225 (Citrullus amarus) identified QTL on chromosomes 2, 5, 9 and 11. A novel QTL on chromosome 5 (Qgsb5.2) associated with resistance to S. citrulli (PVE = 13.3%) was confirmed by genetic mapping. KASP marker assays were developed for selection of Qgsb5.2 to allow breeders to track the allele contributing resistance to GSB, reducing the need for laborious phenotyping. Pyramiding different GSB resistance QTL could be a useful strategy to develop GSB resistant watermelon cultivars.

Chromosome-level genome assembly of bunching onion illuminates genome evolution and flavor formation in Allium crops

The Allium genus is cultivated globally as vegetables, condiments, or medicinal plants and is characterized by large genomes and strong pungency. However, the genome evolution and genomic basis underlying their unique flavor formation remain poorly understood. Herein, we report an 11.27-Gb chromosome-scale genome assembly for bunching onion (A. fistulosum). The uneven bursts of long-terminal repeats contribute to diversity in genome constituents, and dispersed duplication events largely account for gene expansion in Allium genomes. The extensive duplication and differentiation of alliinase and lachrymatory factor synthase manifest as important evolutionary events during flavor formation in Allium crops. Furthermore, differential selective preference for flavor-related genes likely lead to the variations in isoalliin content in bunching onions. Moreover, we reveal that China is the origin and domestication center for bunching onions. Our findings provide insights into Allium genome evolution, flavor formation and domestication history and enable future genome-assisted breeding of important traits in these crops.

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.

An allelic variant in the ACS7 gene promotes primary root growth in watermelon

Gene mining in a C. lanatus × C. amarus population revealed one gene, ACS7, linked to primary root elongation in watermelon. Watermelon is a xerophytic crop characterized by a long primary root and robust lateral roots. Therefore, watermelon serves as an excellent model for studying root elongation and development. However, the genetic mechanism underlying the primary root elongation in watermelon remains unknown. Herein, through bulk segregant analysis we identified a genetic locus, qPRL.Chr03, controlling primary root length (PRL) using two different watermelon species (Citrullus lanatus and Citrullus amarus) that differ in their root architecture. Fine mapping revealed that xaa-Pro dipeptidase and 1-aminocyclopropane-1-carboxylate synthase 7 (ACS7) are candidate regulators of the primary root growth. Allelic variation in the delimited region among 193 watermelon accessions indicated that the long-root alleles might only exist in C. amarus. Interestingly, the discrepancy in PRL among the C. amarus accessions was clearly associated with a nonsynonymous single nucleotide polymorphism variant within the ACS7 gene. The ACS7 expression and ethylene levels in the primary root tips suggested that ethylene is a negative regulator of root elongation in watermelon, as supported by the application of 1-aminocyclopropane-1-carboxylate (ACC, the ethylene precursor) or 2-aminoethoxyvinyl glycine (AVG, an ACS inhibitor). To the best of our knowledge, these findings provide the first description of the genetic basis of root elongation in watermelon. The detected markers of the ACS7 gene will facilitate marker-assisted selection for the PRL trait to improve water and nutrient use efficacy in watermelon and beyond.

Genomic analyses of rice bean landraces reveal adaptation and yield related loci to accelerate breeding

Rice bean (Vigna umbellata) is an underexploited domesticated legume crop consumed for dietary protein in Asia, yet little is known about the genetic diversity of this species. Here, we present a high-quality reference genome for a rice bean landrace (FF25) built using PacBio long-read data and a Hi-C chromatin interaction map, and assess the phylogenetic position and speciation time of rice bean within the Vigna genus. We sequence 440 landraces (two core collections), and GWAS based on data for growth sites at three widely divergent latitudes reveal loci associated with flowering and yield. Loci harboring orthologs of FUL (FRUITFULL), FT (FLOWERING LOCUS T), and PRR3 (PSEUDO-RESPONSE REGULATOR 3) contribute to the adaptation of rice bean from its low latitude center of origin towards higher latitudes, and the landraces which pyramid early-flowering alleles for these loci display maximally short flowering times. We also demonstrate that copy-number-variation for VumCYP78A6 can regulate seed-yield traits. Intriguingly, 32 landraces collected from a mountainous region in South-Central China harbor a recently acquired InDel in TFL1 (TERMINAL FLOWER1) affecting stem determinacy; these materials also have exceptionally high values for multiple human-desired traits and could therefore substantially advance breeding efforts to improve rice bean.

Rice bean is an underexploited legume crop that has many desirable properties against bio and abiotic stresses. Here, the authors report the genome assembly of this species, conduct population genetics studies and reveal the genetic variations associated with adaptation and yield traits.

Genome-wide association study of traits in sacred lotus uncovers MITE-associated variants underlying stamen petaloid and petal number variations

Understanding the genetic variants responsible for floral trait diversity is important for the molecular breeding of ornamental flowers. Widely used in water gardening for thousands of years, the sacred lotus exhibits a wide range of diversity in floral organs. Nevertheless, the genetic variations underlying various morphological characteristics in lotus remain largely unclear. Here, we performed a genome-wide association study of sacred lotus for 12 well-recorded ornamental traits. Given a moderate linkage disequilibrium level of 32.9 kb, we successfully identified 149 candidate genes responsible for seven flower traits and plant size variations, including many pleiotropic genes affecting multiple floral-organ-related traits, such as NnKUP2. Notably, we found a 2.75-kb presence-and-absence genomic fragment significantly associated with stamen petaloid and petal number variations, which was further confirmed by re-examining another independent population dataset with petal number records. Intriguingly, this fragment carries MITE transposons bound by siRNAs and is related to the expression differentiation of a nearby candidate gene between few-petalled and double-petalled lotuses. Overall, these genetic variations and candidate genes responsible for diverse lotus traits revealed by our GWAS highlight the role of transposon variations, particularly MITEs, in shaping floral trait diversity.

QTL associated with Gummy Stem Blight (GSB) resistance in watermelon

Background

Gummy stem blight (GSB), caused by Didymella bryoniae (syn. Stagonosporopsis cucurbitacearum), produces devastating symptoms on whole plants of watermelon (Citrullus lanatus) and other cucurbits, significantly reducing yield and quality. Identification of genetic determinants and sources of resistance to this devastating GSB disease in watermelon is essential for developing resistant varieties.

Results

In this study, we aimed at identifying quantitative trait loci (QTLs) linked to GSB resistance in melon. We identified the genome-wide single nucleotide polymorphisms (SNPs) by genotyping by sequencing (GBS) of an F₂ population developed from C. lanatus lines, ‘PI 279461’ (resistant) ✕ ‘PI 223764’ (susceptible). Inheritance analysis indicated that resistance to GSB is a multi-genic trait in this population. Three QTLs namely, ClGSB1.1, ClGSB10.1, and ClGSB11.1 associated with GSB resistance, explaining approximately 10% of the phenotypic variation, were identified. Among these, the QTL ClGSB1.1 on chromosome 1 is identified as a major QTL harboring five candidate genes associated with GSB resistance including two RLKs (ClC01G014900 and ClC01G015010), two WRKY transcription factors (ClC01G014910 and ClC01G014990), and one AvrRpt-cleavage domain protein (ClC01G015130).

Conclusion

Two high resolution melting (HRM) markers, WmGSB1.1–2 and WmGSB1.1–7 having a high positive correlation with the phenotypic variations, were developed. Five potential candidate genes were predicted to be associated with GSB resistance. These findings will help breeders to develop watermelon cultivars resistant to GSB.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08849-2.

A novel strategy for improving watermelon resistance to cucumber green mottle mosaic virus by exogenous boron application

The molecular mode controlling cucumber green mottle mosaic virus (CGMMV)-induced watermelon blood flesh disease (WBFD) is largely unknown. In this study, we have found that application of exogenous boron suppressed CGMMV infection in watermelon fruit and alleviated WBFD symptoms. Our transcriptome analysis showed that the most up-regulated differentially expressed genes (DEGs) were associated with polyamine and auxin biosynthesis, abscisic acid catabolism, defence-related pathways, cell wall modification, and energy and secondary metabolism, while the down-regulated DEGs were mostly involved in ethylene biosynthesis, cell wall catabolism, and plasma membrane functions. Our virus-induced gene silencing results showed that silencing of SPDS expression in watermelon resulted in a higher putrescine content and an inhibited CGMMV infection correlating with no WBFD symptoms. SBT and TUBB1 were also required for CGMMV infection. In contrast, silencing of XTH23 and PE/PEI7 (low-level lignin, cellulose and pectin) and ATPS1 (low-level glutathione) promoted CGMMV accumulation. Furthermore, RAP2-3, MYB6, WRKY12, H2A, and DnaJ11 are likely to participate in host antiviral resistance. In addition, a higher (spermidine + spermine):putrescine ratio, malondialdehyde content, and lactic acid content were responsible for fruit decay and acidification. Our results provide new knowledge on the roles of boron in watermelon resistance to CGMMV-induced WBFD. This new knowledge can be used to design better control methods for CGMMV in the field and to breed CGMMV resistant watermelon and other cucurbit crops.

Transcriptomic and metabonomic insights into the biocontrol mechanism of Trichoderma asperellum M45a against watermelon Fusarium wilt

Watermelon (Citrullus lanatus) is one of the most popular fruit crops. However, Fusarium wilt (FW) is a serious soil-borne disease caused by Fusarium oxysporum f. sp. niveum (FON) that severely limits the development of the watermelon industry. Trichoderma spp. is an important plant anti-pathogen biocontrol agent. The results of our previous study indicated that Trichoderma asperellum M45a (T. asperellum M45a) could control FW by enhancing the relative abundance of plant growth-promoting rhizobacteria (PGPR) in the rhizosphere of watermelon. However, there are few studies on its mechanism in the pathogen resistance of watermelon. Therefore, transcriptome sequencing of T. asperellum M45a-treated watermelon roots combined with metabolome sequencing of the rhizosphere soil was performed with greenhouse pot experiments. The results demonstrated that T. asperellum M45a could stably colonize roots and significantly increase the resistance-related enzymatic activities (e.g., lignin, cinnamic acid, peroxidase and peroxidase) of watermelon. Moreover, the expression of defense-related genes such as MYB and PAL in watermelon roots significantly improved with the inoculation of T. asperellum M45a. In addition, KEGG pathway analysis showed that a large number of differentially expressed genes were significantly enriched in phenylpropane metabolic pathways, which may be related to lignin and cinnamic acid synthesis, thus further inducing the immune response to resist FON. Furthermore, metabolic analysis indicated that four differential metabolic pathways were enriched in M45a-treated soil, including six upregulated compounds and one down-regulated compound. Among them, galactinol and urea were significantly positively correlated with Trichoderma. Hence, this study provides insight into the biocontrol mechanism of T. asperellum M45a to resist soil-borne diseases, which can guide its industrial application.

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.

Comparative Transcriptome Analysis Reveals Differential Gene Expression in Resistant and Susceptible Watermelon Varieties in Response to Meloidogyne incognita

M. incognita is a major parasitic plant disease in watermelon production, causing serious economic losses. Although there are many studies on root-knot nematode, the resistance mechanism is still unclear. In this study, in order to fully understand the mechanism of watermelon resistance to root-knot nematode, the relatively strongly resistant ‘Hongzi watermelon’ variety and the susceptible ‘M16’ watermelon variety were used as materials, combined with RNA sequencing (RNA-seq), to analyze the expression abundance of resistant and susceptible varieties at 0, 2, 8 and 15 days post-infection (DPI) by M. incognita. The number of differentially expressed genes (DEGs) in the four comparison groups (A0_B0, A1_B1, A2_B2 and A3_B3) was 3645, 2306, 4449 and 2362, respectively, and there were 835 shared DEGs among them. GO annotation and KEGG pathway enrichment analysis showed that 835 DEGs were mainly involved in phenylpropane biosynthesis and carbon metabolism. Furthermore, lignin-biosynthesis-related genes (4CL (4-coumaric acid-CoA ligase), C3H (coumaric acid 3-hydroxylase), CSE (caffeoyl shikimate esterase), COMT (caffeic acid-O-methyltransferase), CCR (cinnamyl CoA reductase) and PRX (peroxidase)), defense-related proteins (UDP-glucoronosyl/UDP-glucosyl transferase, UGT84A13; salicylic acid binding protein, SABP2) and some transcription factors (TFs) were highlighted, which may be potential candidate genes for further analysis in the infection process of M. incognita. These results suggest that watermelon can achieve resistance to M. incognita by increasing the content of lignin and phenols in root or improving ROS level. These RNA-seq data provide new knowledge for future functional studies and will be helpful to further elucidate the molecular mechanism of resistance to M. incognita in watermelon.

Combined nature and human selections reshaped peach fruit metabolome

Background

Plant metabolites reshaped by nature and human beings are crucial for both their lives and human health. However, which metabolites respond most strongly to selection pressure at different evolutionary stages and what roles they undertake on perennial fruit crops such as peach remain unclear.

Results

Here, we report 18,052 significant locus-trait associations, 12,691 expression-metabolite correlations, and 294,676 expression quantitative trait loci (eQTLs) for peach. Our results indicate that amino acids accumulated in landraces may be involved in the environmental adaptation of peaches by responding to low temperature and drought. Moreover, the contents of flavonoids, the major nutrients in fruits, have kept decreasing accompanied by the reduced bitter flavor during both domestication and improvement stages. However, citric acid, under the selection of breeders’ and consumers’ preference for flavor, shows significantly different levels between eastern and western varieties. This correlates with differences in activity against cancer cells in vitro in fruit from these two regions. Based on the identified key genes regulating flavonoid and acid contents, we propose that more precise and targeted breeding technologies should be designed to improve peach varieties with rich functional contents because of the linkage of genes related to bitterness and acid taste, antioxidant and potential anti-cancer activity that are all located at the top of chromosome 5.

Conclusions

This study provides powerful data for future improvement of peach flavor, nutrition, and resistance in future and expands our understanding of the effects of natural and artificial selection on metabolites.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02719-6.

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.

A natural mutation of the NST1 gene arrests secondary cell wall biosynthesis in the seed coat of a hull-less pumpkin accession

Hull-less pumpkins (Cucurbita pepo L.) are naturally occurring novel variants known as oilseed or naked-seeded pumpkins, and are characterized by the absence of a normal lignified seed coat. Due to a specialized seed coat structure, these variants serve as a good model for studying seed coat formation and simplify the processing of pumpkin seeds. However, causal genes for this hull-less trait still remain unknown. Here, by bulked segregant analysis and fine mapping, we found that mutation of a single gene, NAC SECONDARY WALL THICKENING PROMOTING FACTOR 1 (NST1), accounts for the hull-less trait. A 14-bp sequence insertion in the CpNST1 gene causes premature termination of CpNST1 translation, leading to lack of secondary cell wall (SCW) biosynthesis in hull-less seed coats. In situ hybridization analysis provided further evidence for the role of CpNST1 in pumpkin seed coat SCW biosynthesis. Interestingly, through secondary cell wall compositional analysis, we found that the main SCW components differed among cell layers in the seed coat. RNA-seq analysis indicated an upstream role of CpNST1 in the SCW biosynthesis network. Collectively, our findings provide mechanistic insight into seed coat SCW biosynthesis, and a target gene for breeders to introduce this hull-less trait for commercial exploitation.

Melatonin-mediated postharvest quality and antioxidant properties of fresh fruits: A comprehensive meta-analysis

At postharvest, fruits have a short shelf life. Recently, there has been much literature on the effects of melatonin on the postharvest quality of horticultural crops. However, reports of various findings comprise mixed claims and product-specific conclusions. Therefore, a meta-analysis systematically dissects the comprehensive effect on several fruits. In this meta-analysis, standard mean difference (SMD) was adopted using a random-effect model. The study used 36 articles and isolated 24 indicator parameters of postharvest quality and antioxidant properties based on the inclusion criteria. As exhibited in the forest plot, melatonin reduced chilling injury, weight loss, respiration rate, and ethylene content (SMD -0.90, 95% CI [-1.14, -0.65]; I²  = 81%; p < .00001). Similarly, the application of melatonin significantly suppressed electrolyte leakage, malondialdehyde (MDA), hydrogen peroxide, superoxide anion, lipoxygenase, and polyphenol oxidase (SMD -0.89, 95% CI [-1.09, -0.69]; I²  = 70%; p < .00001). In addition, exogenous melatonin application induced endogenous melatonin content, phenolic content, and flavonoid and anthocyanin contents (SMD 1.15, 95% CI [0.91, 1.39]; I²  = 71%; p = .01). Moreover, melatonin treatment enhanced antioxidant activities (catalase, superoxide dismutase, peroxidase, ascorbate peroxidase, and phenylalanine ammonia-lyse) (SMD 1.37, 95% CI [1.03, 1.71]; I²  = 86%; p < .00001). Thus, in the whole study, the overall effect was significantly high in treated fruit (p < .0001), and the overall heterogeneity was above (I² ) > 70%. In addition, the funnel plot showed symmetry in the most selected studies. To sum up, the result gives a further understanding of melatonin's capabilities in reducing postharvest losses and maintaining the quality of fresh fruits.

The integrated genomics of crop domestication and breeding

As a major event in human civilization, wild plants were successfully domesticated to be crops, largely owing to continuing artificial selection. Here, we summarize new discoveries made during the past decade in crop domestication and breeding. The construction of crop genome maps and the functional characterization of numerous trait genes provide foundational information. Approaches to read, interpret, and write complex genetic information are being leveraged in many plants for highly efficient de novo or re-domestication. Understanding the underlying mechanisms of crop microevolution and applying the knowledge to agricultural productions will give possible solutions for future challenges in food security.

Development of Whole Genome SNP-CAPS Markers and Preliminary QTL Mapping of Fruit Pedicel Traits in Watermelon

Fruit pedicel (FP) is an important determinant of premium fruit quality that directly affects commercial market value. However, in-depth molecular and genetic basis of pedicel-related traits has not been identified in watermelon. Herein, a quantitative trait locus (QTL) mapping strategy was used to identify the potential genetic regions controlling FP traits based on newly derived whole-genome single nucleotide polymorphism based cleaved amplified polymorphism sequence (SNP-CAPS) markers. Next-generation sequencing based whole-genome re-sequencing of two watermelon parent lines revealed 98.30 and 98.40% of average coverage, 4,989,869 SNP variants, and 182,949 CAPS loci pairs across the reference genome, respectively. A total of 221 sets of codominant markers exhibited 46.42% polymorphism rate and were effectively genotyped within 100-F_2:3 derived mapping population. The developed linkage map covered a total of 2,630.49 cM genetic length with averaged 11.90 cM, and depicted a valid marker-trait association. In total, 6 QTLs (qFPL4.1, qFPW4.1, qFPD2.1, qFPD2.2, qFPD8.1, qFPD10.1) were mapped with five major effects and one minor effect between the whole genome adjacent markers positioned over distinct chromosomes (02, 04, 08, 10), based on the ICIM-ADD mapping approach. These significant QTLs were similarly mapped in delimited flanking regions of 675.10, 751.38, 859.24, 948.39, and 947.51 kb, which collectively explained 8.64-13.60% PVE, respectively. A highly significant and positive correlation was found among the observed variables. To our knowledge, we first time reported the mapped QTLs/genes affecting FP traits of watermelon, and our illustrated outcomes will deliver the potential insights for fine genetic mapping as well as functional gene analysis through MAS-based breeding approaches.

Identification, Characterization and Function of Orphan Genes Among the Current Cucurbitaceae Genomes

Orphan genes (OGs) that are missing identifiable homologs in other lineages may potentially make contributions to a variety of biological functions. The Cucurbitaceae family consists of a wide range of fruit crops of worldwide or local economic significance. To date, very few functional mechanisms of OGs in Cucurbitaceae are known. In this study, we systematically identified the OGs of eight Cucurbitaceae species using a comparative genomics approach. The content of OGs varied widely among the eight Cucurbitaceae species, ranging from 1.63% in chayote to 16.55% in wax gourd. Genetic structure analysis showed that OGs have significantly shorter protein lengths and fewer exons in Cucurbitaceae. The subcellular localizations of OGs were basically the same, with only subtle differences. Except for aggregation in some chromosomal regions, the distribution density of OGs was higher near the telomeres and relatively evenly distributed on the chromosomes. Gene expression analysis revealed that OGs had less abundantly and highly tissue-specific expression. Interestingly, the largest proportion of these OGs was significantly more tissue-specific expressed in the flower than in other tissues, and more detectable expression was found in the male flower. Functional prediction of OGs showed that (1) 18 OGs associated with male sterility in watermelon; (2) 182 OGs associated with flower development in cucumber; (3) 51 OGs associated with environmental adaptation in watermelon; (4) 520 OGs may help with the large fruit size in wax gourd. Our results provide the molecular basis and research direction for some important mechanisms in Cucurbitaceae species and domesticated crops.

The vegetable SNP database: An integrated resource for plant breeders and scientists

Single nucleotide polymorphisms (SNPs) are widely used in genetic research and molecular breeding. To date, the genomes of many vegetable crops have been assembled, and hundreds of core germplasms for each vegetable have been sequenced. However, these data are not currently easily accessible because they are stored on different public databases. Therefore, a vegetable crop SNP database should be developed that hosts SNPs demonstrated to have a high success rate in genotyping for genetic research (herein, "alpha SNPs"). We constructed a database (VegSNPDB, http://www.vegsnpdb.cn/) containing the sequence data of 2032 germplasms from 16 vegetable crop species. VegSNPDB hosts 118,725,944 SNPs of which 4,877,305 were alpha SNPs. SNPs can be searched by chromosome number, position, SNP type, genetic population, or specific individuals, as well as the values of MAF, PIC, and heterozygosity. We hope that VegSNPDB will become an important SNP database for the vegetable research community.

ClZISO mutation leads to photosensitive flesh in watermelon

The mutation of ClZISO identified in EMS-induced watermelon leads to photosensitive flesh in watermelon. Watermelon (Citrullus lanatus) has a colorful flesh that attracts consumers and benefits human health. We developed an ethyl-methanesulfonate mutation library in red-fleshed line '302' to create new flesh color lines and found a yellow-fleshed mutant which accumulated ζ-carotene. The initial yellow color of this mutant can be photobleached within 10 min under intense sunlight. A long-term light-emitting diode (LED) light treatment turned flesh color from yellow to pink. We identified this unique variation as photosensitive flesh mutant ('psf'). Using bulked segregant analysis, we fine-mapped an EMS-induced G-A transversion in 'psf' which leads to a premature stop codon in 15-cis-ζ-carotene isomerase (ClZISO) gene. We detected that wild-type ClZISO is expressed in chromoplasts to catalyze the conversion of 9,15,9'-tri-cis-ζ-carotene to 9,9'-di-cis-ζ-carotene. The truncated ClZISO_mu protein in psf lost this catalytic function. Light treatment can partially compensate ClZISO_mu isomerase activity via photoisomerization in vitro and in vivo. Transcriptome analysis showed that most carotenoid biosynthesis genes in psf were downregulated. The dramatic increase of ABA content in flesh with fruit development was blocked in psf. This study explores the molecular mechanism of carotenoid biosynthesis in watermelon and provides a theoretical and technical basis for breeding different flesh color lines in watermelon.

Efficient genetic transformation and CRISPR/Cas9-mediated genome editing of watermelon assisted by genes encoding developmental regulators

Cucurbitaceae is an important family of flowering plants containing multiple species of important food plants, such as melons, cucumbers, squashes, and pumpkins. However, a highly efficient genetic transformation system has not been established for most of these species (Nanasato and Tabei, 2020). Watermelon (Citrullus lanatus), an economically important and globally cultivated fruit crop, is a model species for fruit quality research due to its rich diversity of fruit size, shape, flavor, aroma, texture, peel and flesh color, and nutritional composition (Guo et al., 2019). Through pan-genome sequencing, many candidate loci associated with fruit quality traits have been identified (Guo et al., 2019). However, few of these loci have been validated. The major barrier is the low transformation efficiency of the species, with only few successful cases of genetic transformation reported so far (Tian et al., 2017; Feng et al., 2021; Wang JF et al., 2021; Wang YP et al., 2021). For example, Tian et al. (2017) obtained only 16 transgenic lines from about 960 cotyledon fragments, yielding a transformation efficiency of 1.67%. Therefore, efficient genetic transformation could not only facilitate the functional genomic studies in watermelon as well as other horticultural species, but also speed up the transgenic and genome-editing breeding.

Unravelling the Genetic Architecture of Rust Resistance in the Common Bean (Phaseolus vulgaris L.) by Combining QTL-Seq and GWAS Analysis

The common bean (Phaseolus vulgaris L.) is the most important legume crop directly used for human consumption worldwide. Bean rust, caused by Uromyces appendiculatus, is a devastating disease and usually causes severe loss of seed yield and pod quality. Deployment of resistant cultivars is the best strategy to combat this disease. However, despite being the largest snap bean-producing country, the genetic basis research of rust resistance has largely lagged in China. In this study, an RIL population and a diversity panel were evaluated for rust resistance against a purified rust isolate Cua-LS using a detached leaf assay. Deploying a QTL-Seq analysis in the RIL population, a 1.81 Mb interval on chromosome 4, a 2.73 Mb interval on chromosome 5 and a 1.26 Mb interval on chromosome 6 were identified as major QTLs for rust resistance, designated as Qur-1, Qur-2 and Qur-3, respectively. Through a GWAS diversity panel, 64 significant SNPs associated with rust resistance were detected, distributed in all 11 chromosomes and explaining 19–49% of the phenotypic variation. Synteny analysis showed that Qur-2 was validated in GWAS, but the rust QTL/SNPs detected in our study were different from the known genes, except Ur-11. A total of 114 candidate genes, including the typical NBS-LRR genes, protein kinase superfamily proteins and ABC transporter family proteins, were identified and proposed as the likely candidates. The identified 17 resistant accessions will enrich the resistant germplasm resources, and the detected QTLs/SNPs will facilitate the molecular breeding of rust resistance in the common bean.

Quantitative Transcriptomic and Proteomic Analysis of Fruit Development and Ripening in Watermelon (Citrullus lanatus)

Fruit ripening is a highly complicated process, which is modulated by phytohormones, signal regulators and environmental factors playing in an intricate network that regulates ripening-related genes expression. Although transcriptomics is an effective tool to predict protein levels, protein abundances are also extensively affected by post-transcriptional and post-translational regulations. Here, we used RNA sequencing (RNA-seq) and tandem mass tag (TMT)-based quantitative proteomics to study the comprehensive mRNA and protein expression changes during fruit development and ripening in watermelon, a non-climacteric fruit. A total of 6,226 proteins were quantified, and the large number of quantitative proteins is comparable to proteomic studies in model organisms such as Oryza sativa L. and Arabidopsis. Base on our proteome methodology, integrative analysis of the transcriptome and proteome showed that the mRNA and protein levels were poorly correlated, and the correlation coefficients decreased during fruit ripening. Proteomic results showed that proteins involved in alternative splicing and the ubiquitin proteasome pathway were dynamically expressed during ripening. Furthermore, the spliceosome and proteasome were significantly enriched by Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, suggesting that post-transcriptional and post-translational mechanisms might play important roles in regulation of fruit ripening-associated genes expression, which might account for the poor correlation between mRNAs and proteins during fruit ripening. Our comprehensive transcriptomic and proteomic data offer a valuable resource for watermelon research, and provide new insights into the molecular mechanisms underlying the complex regulatory networks of fruit ripening.

Integrative analysis of transcriptomic and metabolomic profiles reveals new insights into the molecular foundation of fruit quality formation in Citrullus lanatus (Thunb.) Matsum. & Nakai

In this study, an integrated transcriptome and metabolome analysis was used to explore the molecular foundation of fruit quality in two parent lines of Citrullus lanatus with distinct flesh characteristics, including ‘14-1’ (sweet, red, and soft) and ‘W600’ (bitter, light yellow, and firm), as well as the corresponding F1 population (bitter, light yellow, and firm). Numerous differentially expressed genes (DEGs) were identified in the fruit samples: 3,766 DEGs for ‘14-1’ vs. ‘W600’, 2,767 for ‘14-1’ vs. F1, and 1,178 for F1 vs. ‘W600’ at the transition stage; and 4,221 for ‘14-1’ vs. ‘W600’, 2,447 for ‘14-1’ vs. F1, and 446 for F1 vs. ‘W600’ at the maturity stage. Weighted gene co-expression network analysis (WGCNA) revealed that a gene module including 1,111 DEGs was closely associated with flesh taste and color, and another gene module including 1,575 DEGs contributed significantly to flesh texture. The metabolomic results showed that there were 447 differential metabolites (DMs) for ‘14-1’ vs. ‘W600’ fruits, 394 for ‘14-1’ vs. F1, and 298 for F1 vs. ‘W600’ at the maturity stage. Combining WGNCA and metabolomic results, several DEGs and DMs were further identified as hub players in fruit quality formation: six DEGs with four DMs for flesh sweetness; six DEGs with 13 DMs for bitterness; nine DEGs with 10 DMs for flesh color; and nine DEGs with four DMs for flesh texture. Altogether, these observations not only expand our knowledge of the molecular basis of fruit quality in watermelon, but also provide potential targets for future watermelon improvement.

Pangenomes as a Resource to Accelerate Breeding of Under-Utilised Crop Species

Pangenomes are a rich resource to examine the genomic variation observed within a species or genera, supporting population genetics studies, with applications for the improvement of crop traits. Major crop species such as maize (Zea mays), rice (Oryza sativa), Brassica (Brassica spp.), and soybean (Glycine max) have had pangenomes constructed and released, and this has led to the discovery of valuable genes associated with disease resistance and yield components. However, pangenome data are not available for many less prominent crop species that are currently under-utilised. Despite many under-utilised species being important food sources in regional populations, the scarcity of genomic data for these species hinders their improvement. Here, we assess several under-utilised crops and review the pangenome approaches that could be used to build resources for their improvement. Many of these under-utilised crops are cultivated in arid or semi-arid environments, suggesting that novel genes related to drought tolerance may be identified and used for introgression into related major crop species. In addition, we discuss how previously collected data could be used to enrich pangenome functional analysis in genome-wide association studies (GWAS) based on studies in major crops. Considering the technological advances in genome sequencing, pangenome references for under-utilised species are becoming more obtainable, offering the opportunity to identify novel genes related to agro-morphological traits in these species.

Comparative genomics analysis of bHLH genes in cucurbits identifies a novel gene regulating cucurbitacin biosynthesis

The basic helix-loop-helix (bHLH) family of transcription factors (TFs) participate in a variety of biological regulatory processes in plants, and have undergone significant expansion during land plant evolution by gene duplications. In cucurbit crops, several bHLH genes have been found to be responsible for the agronomic traits such as bitterness. However, the characterization of bHLH genes across the genomes of cucurbit species has not been reported, and how they have evolved and diverged remains largely unanswered. Here we identified 1160 bHLH genes in seven cucurbit crops and performed a comprehensive comparative genomics analysis. We determined orthologous and paralogous bHLH genes across cucurbit crops by syntenic analysis between or within species. Orthology and phylogenetic analysis of the tandem-duplicated bHLH genes in the Bt cluster which regulate the biosynthesis of cucurbitacins suggest that this cluster is derived from three ancestral genes after the cucurbit-common tetraploidization event. Interestingly, we identified a new conserved cluster paralogous to the Bt cluster that includes two tandem bHLH genes, and the evolutionary history and expression profiles of these two genes in the new cluster suggest the involvement of one gene (Brp) in the regulation of cucurbitacin biosynthesis in roots. Further biochemical and transgenic assays in melon hairy roots support the function of Brp. This study provides useful information for further investigating the functions of bHLH TFs and novel insights into the regulation of cucurbitacin biosynthesis in cucurbit crops and other plants.

Chromosomal genome and population genetic analyses to reveal genetic architecture, breeding history and genes related to cadmium accumulation in Lentinula edodes

BACKGROUND

Lentinula edodes (Berk.) is the second most productive mushroom in the world. It contains compounds effective for antiviral, antitumor, antioxidant and immune regulation. Although genomes have previously been reported for this species, a high-quality chromosome-level reference for L. edodes is unavailable. This hinders detailed investigation of population genetics, breeding history of strains and genes related to environmental stress responses.

RESULTS

A high-quality chromosome-level genome was constructed. We separated a monokaryon from protoplasts of the commercial L. edodes strain L808 and assembled the genome of L. edodes using PacBio long-read and Illumina short-read sequencing, along with the high-throughput chromatin conformation capture (Hi-C) technique. We assembled a 45.87 Mb genome, and 99% of the sequences were anchored onto 10 chromosomes. The contig and scaffold N50 length were 2.17 and 4.94 Mb, respectively. Over 96% of the complete Benchmarking Universal Single-Copy Orthologs (BUSCO) were identified, and 9853 protein-coding genes were predicted. We performed population genome resequencing using 34 wild strains and 65 commercial cultivars of L. edodes originating from China, Japan, the United States and Australia. Based on whole-genome variants, we showed substantial differences in the Chinese wild population, which divided into different branches according to the main areas of their geographical distribution. We also determined the breeding history of L. edodes at the molecular level, and demonstrated that the cultivated strains in China mainly originated from wild strains from China and Northeast Asia. Phenotypic analysis showed that 99 strains exhibited differences on the Cd accumulation. Three significant loci in the of L. edodes genome were identified using the genome-wide association study (GWAS) of Cd accumulation traits. Functional genes associated with Cd accumulation traits were related to DNA ligase and aminoacyl tRNA synthetase, indicating that DNA damage repair and in vivo protein translation may be responses to Cd stress.

CONCLUSIONS

A high-quality chromosome-level genome and population genetic data of L. edodes provide genetic resources for functional genomic, evolutionary and artificial breeding studies for L. edodes.

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 association analysis provides molecular insights into the natural variation of watermelon seed size

Seed-consumption watermelon tend to have larger-sized seeds, while flesh-consumed watermelon often require relatively smaller seed. Therefore, the seed size of watermelon has received extensive attention from consumers and breeders. However, the study on the natural variation and genetic mechanism of watermelon seed size is not clear enough. In the present study, 100 seed weight, seed hilum length, seed length, seed width, and seed thickness in 197 watermelon accessions were examined. Furthermore, association analysis was conducted between seed size traits and high-quality SNP data. The results revealed that there was a strong correlation between the five seed traits. And seed enlargement was an important feature during watermelon seed size domestication. Meanwhile, the seed consumption biological species C. mucosospermu and C. lanatus edible seed watermelon had a significantly bigger seed size than other species's. Eleven non-repeating significant SNPs above the threshold line were obtained by GWAS analysis. Four of them on chromosome 5 were considered to be closely associated with seed size traits, i.e. S5: 32250307, S5: 32250454, S5: 32256177, S5: 32260870, which could be used as potential molecular markers for the breeding of watermelon cultivars with target seed size. In addition, combined with gene annotation information and previous reports, five genes near the four significant SNPs may regulate seed size. And qRT-PCR analysis showed that two genes Cla97C05G104360 and Cla97C05G104380, which may be involved in abscisic acid metabolism, may play an important role in regulating the seed size of watermelon. Our findings provide molecular insights into natural variation in watermelon seed size, and gives valuable information of molecular marker-assisted breeding.