The draft genome of watermelon (Citrullus lanatus) and resequencing of 20 diverse accessions

Triumphs of genomic-assisted breeding in crop improvement

Conventional breeding approaches have played a significant role in meeting the food demand remarkably well until now. However, the increasing population, yield plateaus in certain crops, and limited recombination necessitate using genomic resources for genomics-assisted crop improvement programs. As a result of advancements in the next-generation sequence technology, GABs have developed dramatically to characterize allelic variants and facilitate their rapid and efficient incorporation in crop improvement programs. Genomics-assisted breeding (GAB) has played an important role in harnessing the potential of modern genomic tools, exploiting allelic variation from genetic resources and developing cultivars over the past decade. The availability of pangenomes for major crops has been a significant development, albeit with varying degrees of completeness. Even though adopting these technologies is essentially determined on economic grounds and cost-effective assays, which create a wealth of information that can be successfully used to exploit the latent potential of crops. GAB has been instrumental in harnessing the potential of modern genomic resources and exploiting allelic variation for genetic enhancement and cultivar development. GAB strategies will be indispensable for designing future crops and are expected to play a crucial role in breeding climate-smart crop cultivars with higher nutritional value.

Chromosome-level assembly of Lindenbergia philippensis and comparative genomic analyses shed light on genome evolution in Lamiales

Lamiales, comprising over 23,755 species across 24 families, stands as a highly diverse and prolific plant group, playing a significant role in the cultivation of horticultural, ornamental, and medicinal plant varieties. Whole-genome duplication (WGD) and its subsequent post-polyploid diploidization (PPD) process represent the most drastic type of karyotype evolution, injecting significant potential for promoting the diversity of this lineage. However, polyploidization histories, as well as genome and subgenome fractionation following WGD events in Lamiales species, are still not well investigated. In this study, we constructed a chromosome-level genome assembly of Lindenbergia philippensis (Orobanchaceae) and conducted comparative genomic analyses with 14 other Lamiales species. L. philippensis is positioned closest to the parasitic lineage within Orobanchaceae and has a conserved karyotype. Through a combination of Ks analysis and syntenic depth analysis, we reconstructed and validated polyploidization histories of Lamiales species. Our results indicated that Primulina huaijiensis underwent three rounds of diploidization events following the γ-WGT event, rather than two rounds as reported. Besides, we reconfirmed that most Lamiales species shared a common diploidization event (L-WGD). Subsequently, we constructed the Lamiales Ancestral Karyotype (LAK), comprising 11 proto-chromosomes, and elucidated its evolutionary trajectory, highlighting the highly flexible reshuffling of the Lamiales paleogenome. We identified biased fractionation of subgenomes following the L-WGD event across eight species, and highlighted the positive impacts of non-WGD genes on gene family expansion. This study provides novel genomic resources and insights into polyploidy and karyotype remodeling of Lamiales species, essential for advancing our understanding of species diversification and genome evolution.

Telomere-to-telomere Citrullus super-pangenome provides direction for watermelon breeding

To decipher the genetic diversity within the cucurbit genus Citrullus, we generated telomere-to-telomere (T2T) assemblies of 27 distinct genotypes, encompassing all seven Citrullus species. This T2T super-pangenome has expanded the previously published reference genome, T2T-G42, by adding 399.2 Mb and 11,225 genes. Comparative analysis has unveiled gene variants and structural variations (SVs), shedding light on watermelon evolution and domestication processes that enhanced attributes such as bitterness and sugar content while compromising disease resistance. Multidisease-resistant loci from Citrullus amarus and Citrullus mucosospermus were successfully introduced into cultivated Citrullus lanatus. The SVs identified in C. lanatus have not only been inherited from cordophanus but also from C. mucosospermus, suggesting additional ancestors beyond cordophanus in the lineage of cultivated watermelon. Our investigation substantially improves the comprehension of watermelon genome diversity, furnishing comprehensive reference genomes for all Citrullus species. This advancement aids in the exploration and genetic enhancement of watermelon using its wild relatives.

Genome-Wide Analysis and Expression Profiling of Lectin Receptor-like Kinase Genes in Watermelon (Citrullus lanatus)

Watermelon is one of the most important edible plants worldwide. Owing to its special cultivation conditions, watermelon is exposed to many biological and abiotic stresses during its development. Lectin receptor-like kinases (LecRLKs) are plant-specific membrane proteins that play important roles in sensing and responding to environmental stimuli. Although the LecRLK gene family has been identified in a variety of plants, a comprehensive analysis has not yet been undertaken in watermelon. In this study, 61 putative LecRLK genes were identified in watermelon, consisting of 36 G-type, 24 L-type, and 1 C-type LecRLK genes. They were distributed in clusters on chromosomes, and members from the same subfamily were mostly clustered together. The analysis of the phylogenetic tree and conserved motif indicated that there were obvious differences among three ClaLecRLK subfamilies, and there was also rich diversity in the C-terminal within subfamilies. A collinear analysis revealed that the evolution of the ClaLecRLK gene family in different Cucurbitaceae crops was asynchronous. Furthermore, the analysis of the ClaLecRLK protein structure showed that not all proteins contained signal peptides and a single transmembrane domain. A subcellular localization assay confirmed that the number and position of transmembrane domains did not affect ClaLecRLK protein localization in cells. Transcriptome data revealed distinct expression patterns of LecRLK genes of watermelon in various tissues, and their responses to different fungi infection were also significantly different. Finally, the potential binding sites of the ClaLecRLK genes targeted by miRNA were predicted. This study enhances the understanding of the characteristics and functions of the LecRLK gene family in watermelon and opens up the possibility of exploring the roles that LecRLK genes may play in the life cycle of Cucurbitaceae plants.

Quantitative Trait Loci Mapping and Comparative Transcriptome Analysis of Fruit Weight (FW) in Watermelon (Citrullus lanatus L.)

The watermelon (Citrullus lanatus L.) holds substantial economic value as a globally cultivated horticultural crop. However, the genetic architecture of watermelon fruit weight (FW) remains poorly understood. In this study, we used sh14-11 with small fruit and N14 with big fruit to construct 100 recombinant inbred lines (RILs). Based on whole-genome resequencing (WGR), 218,127 single nucleotide polymorphisms (SNPs) were detected to construct a high-quality genetic map. After quantitative trait loci (QTL) mapping, a candidate interval of 31-38 Mb on chromosome 2 was identified for FW. Simultaneously, the bulked segregant analysis (BSA) in the F2 population corroborated the identification of the same interval, encompassing the homologous gene linked to the known FW-related gene fas. Additionally, RNA-seq was carried out across 11 tissues from sh14-11 and N14, revealing expression profiles that identified 1695 new genes and corrected the annotation of 2941 genes. Subsequent differential expression analysis unveiled 8969 differentially expressed genes (DEGs), with 354 of these genes exhibiting significant differences across four key developmental stages. The integration of QTL mapping and differential expression analysis facilitated the identification of 14 FW-related genes, including annotated TGA and NAC transcription factors implicated in fruit development. This combined approach offers valuable insights into the genetic basis of FW, providing crucial resources for enhancing watermelon cultivation.

Deciphering recursive polyploidization in Lamiales and reconstructing their chromosome evolutionary trajectories

Lamiales is an order of core eudicots with abundant diversity, and many Lamiales plants have important medicinal and ornamental values. Here, we comparatively reanalyzed 11 Lamiales species with well-assembled genome sequences and found evidence that Lamiales plants, in addition to a hexaploidization or whole-genome triplication (WGT) shared by core eudicots, experienced further polyploidization events, establishing new groups in the order. Notably, we identified a whole-genome duplication (WGD) occurred just before the split of Scrophulariaceae from the other Lamiales families, such as Acanthaceae, Bignoniaceae, and Lamiaceae, suggesting its likely being the causal reason for the establishment and fast divergence of these families. We also found that a WGT occurred ∼68 to 78 million years ago (Mya), near the split of Oleaceae from the other Lamiales families, implying that it may have caused their fast divergence and the establishment of the Oleaceae family. Then, by exploring and distinguishing intra- and intergenomic chromosomal homology due to recursive polyploidization and speciation, respectively, we inferred that the Lamiales ancestral cell karyotype had 11 proto-chromosomes. We reconstructed the evolutionary trajectories from these proto-chromosomes to form the extant chromosomes in each Lamiales plant under study. We must note that most of the inferred 11 proto-chromosomes, duplicated during a WGD thereafter, have been well preserved in jacaranda (Jacaranda mimosifolia) genome, showing the credibility of the present inference implementing a telomere-centric chromosome repatterning model. These efforts are important to understand genome repatterning after recursive polyploidization, especially shedding light on the origin of new plant groups and angiosperm cell karyotype evolution.

Arbuscular mycorrhizal fungi by inducing watermelon roots secretion phthalates, altering soil enzyme activity and bacterial community composition to alleviate the watermelon wilt

BACKGROUND

Long-term continuous cropping has resulted in the frequent occurrence of fusarium wilt of watermelon (Citrullus lanatus). AMF inoculation can alleviate the continuous cropping barrier and reduce the incidence of fusarium wilt of watermelon. Our previous study found that the root exudates of mycorrhizal watermelon can enhance watermelon resistance to this disorder. It is necessary to further isolate and identify the specific compounds in root exudates of mycorrhizal watermelon and explore their control effects on fusarium wilt of continuous cropping watermelon.

RESULT

The results of this study showed that the root system of watermelon seedlings inoculated with AMF (Funneliformis mosseae or Glomus versiforme) secreted diisooctyl phthalate (A) and dibutyl phthalate (B). Compared with water treatment, treatment with 0.1 ml/L (A1, B1), 0.5 ml/L (A2, B2) and 1 ml/L (A3, B3) of A or B significantly increased soil enzyme activities, the numbers of bacteria and actinomycetes, and the bacteria/fungi ratio in the rhizosphere. Furthermore, the Disease indexes (DI) of A1 and B3 were 25% and 20%, respectively, while the prevention and control effects (PCE) were 68.8% and 75%, respectively. In addition, diisooctyl phthalate or dibutyl phthalate increased the proportions of Gemmatimonadetes, Chloroflexi, and Acidobacteria in the rhizosphere of continuous cropping watermelon, and decreased the proportions of Proteobacteria and Firmicutes, with Novosphingobium, Kaistobacter, Bacillus, and Acinetobacter as the predominant bacteria. Compared with the water treatment, the abundance of Neosphingosaceae, Kateybacterium and Bacillus in the A1 group was increased by 7.33, 2.14 and 2.18 times, respectively, while that in the B2 group was increased by 60.05%, 80.24% and 1 time, respectively. In addition, exogenous diisooctyl phthalate and dibutyl phthalate were shown to promote growth parameters (vine length, stem diameter, fresh weight and dry weight) and antioxidant enzyme system activities (SOD, POD and CAT) of continuous cropping watermelon.

CONCLUSION

Lower watermelon fusarium wilt incidence in mycorrhizal watermelons was associated with phthalate secretion in watermelons after AMF inoculation. Exogenous diisooctyl phthalate and dibutyl phthalate could alleviate the continuous cropping disorder of watermelon, reduce the incidence of fusarium wilt, and promote the growth of watermelon by increasing the enzyme activities and the proportion of beneficial bacteria in rhizosphere soil. In addition, the low concentration of phthalate diisooctyl and high concentration of phthalic acid dibutyl works best. Therefore, a certain concentration of phthalates in the soil can help alleviate continuous cropping obstacles.

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.

Identification of Key Genes of Fruit Shape Variation in Jujube with Integrating Elliptic Fourier Descriptors and Transcriptome

Jujube (Ziziphus jujuba) exhibits a rich diversity in fruit shape, with natural occurrences of gourd-like, flattened, and other special shapes. Despite the ongoing research into fruit shape, studies integrating elliptical Fourier descriptors (EFDs) with both Short Time-series Expression Miner (STEM) and weighted gene co-expression network analysis (WGCNA) for gene discovery remain scarce. In this study, six cultivars of jujube fruits with distinct shapes were selected, and samples were collected from the fruit set period to the white mature stage across five time points for shape analysis and transcriptome studies. By combining EFDs with WGCNA and STEM, the study aimed to identify the critical periods and key genes involved in the formation of jujube fruit shape. The findings indicated that the D25 (25 days after flowering) is crucial for the development of jujube fruit shape. Moreover, ZjAGL80, ZjABI3, and eight other genes have been implicated to regulate the shape development of jujubes at different periods of fruit development, through seed development and fruit development pathway. In this research, EFDs were employed to precisely delineate the shape of jujube fruits. This approach, in conjunction with transcriptome, enhanced the precision of gene identification, and offered an innovative methodology for fruit shape analysis. This integration facilitates the advancement of research into the morphological characteristics of plant fruits, underpinning the development of a refined framework for the genetic underpinnings of fruit shape variation.

Identification of allelic relationship and translocation region among chromosomal translocation lines that leads to less-seed watermelon

Less-seed and seedless traits are desirable characteristics in watermelon (Citrullus lanatus). Hybridization between watermelon chromosomal translocated lines and wild lines significantly reduced seed counts in the hybrid fruits, approaching even seedless. However, the allelic relationships and the chromosomal translocation breakpoints from different sources are unclear, which limits their utility in breeding practices. This study focused on three groups of chromosomal translocation materials from different sources and conducted inheritance and allelic relationship analysis of translocation points. The results from third-generation genome sequencing and fluorescence in situ hybridization (FISH) revealed that the specific translocations in the naturally mutated material MT-a involved reciprocal translocations between Chr6 and Chr10. The Co60γ radiation-induced mutant material MT-b involved reciprocal translocations between Chr1 and Chr5, Chr4 and Chr8. The Co60γ radiation-induced mutant material MT-c involved complex translocations among Chr1, Chr5, and Chr11. Cytological observation showed that heterozygous translocation hybrids showed chromosomal synapsis abnormalities during meiotic diakinesis. Further, dominant and codominant molecular markers were developed on both sides of the translocation breakpoints, which could facilitate rapid and efficient identification of chromosome translocation lines. This study provides technical guidance for utilizing chromosomal translocation materials in the development of less-seed watermelon varieties.

Promoter variations of ClERF1 gene determines flesh firmness in watermelon

BACKGROUND

Flesh firmness is a critical factor that influences fruit storability, shelf-life and consumer's preference as well. However, less is known about the key genetic factors that are associated with flesh firmness in fresh fruits like watermelon.

RESULTS

In this study, through bulk segregant analysis (BSA-seq), we identified a quantitative trait locus (QTL) that influenced variations in flesh firmness among recombinant inbred lines (RIL) developed from cross between the Citrullus mucosospermus accession ZJU152 with hard-flesh and Citrullus lanatus accession ZJU163 with soft-flesh. Fine mapping and sequence variations analyses revealed that ethylene-responsive factor 1 (ClERF1) was the most likely candidate gene for watermelon flesh firmness. Furthermore, several variations existed in the promoter region between ClERF1 of two parents, and significantly higher expressions of ClERF1 were found in hard-flesh ZJU152 compared with soft-flesh ZJU163 at key developmental stages. DUAL-LUC and GUS assays suggested much stronger promoter activity in ZJU152 over ZJU163. In addition, the kompetitive allele-specific PCR (KASP) genotyping datasets of RIL populations and germplasm accessions further supported ClERF1 as a possible candidate gene for fruit flesh firmness variability and the hard-flesh genotype might only exist in wild species C. mucosospermus. Through yeast one-hybrid (Y1H) and dual luciferase assay, we found that ClERF1 could directly bind to the promoters of auxin-responsive protein (ClAux/IAA) and exostosin family protein (ClEXT) and positively regulated their expressions influencing fruit ripening and cell wall biosynthesis.

CONCLUSIONS

Our results indicate that ClERF1 encoding an ethylene-responsive factor 1 is associated with flesh firmness in watermelon and provide mechanistic insight into the regulation of flesh firmness, and the ClERF1 gene is potentially applicable to the molecular improvement of fruit-flesh firmness by design breeding.

Mapping and validation of Fusarium wilt race 2 resistance QTL from Citrullus amarus line USVL246-FR2

KEY MESSAGE

Fon race 2 resistant QTLs were identified on chromosomes 8 and 9. Families homozygous for resistance alleles at a haplotype of three KASP markers had 42% lower disease severity than those with susceptible alleles in an independent, interspecific validation population confirming their utility for introgression of Fusarium wilt resistance. Fusarium oxysporum f. sp. niveum (Fon) race 2 causes Fusarium wilt in watermelon and threatens watermelon production worldwide. Chemical management options are not effective, and no resistant edible watermelon cultivars have been released. Implementation of marker-assisted selection to develop resistant cultivars requires identifying sources of resistance and the underlying quantitative trait loci (QTL), developing molecular markers associated with the QTL, and validating marker-phenotype associations with an independent population. An intraspecific Citrullus amarus recombinant inbred line population from a cross of resistant USVL246-FR2 and susceptible USVL114 was used for mapping Fon race 2 resistance QTL. KASP markers were developed (N = 51) for the major QTL on chromosome 9 and minor QTL on chromosomes 1, 6, and 8. An interspecific F2:3 population was developed from resistance donor USVL246-FR2 (C. amarus) and a susceptible cultivar 'Sugar Baby' (Citrullus lanatus) to validate the utility of the markers for introgression of resistance from the wild crop relative into cultivated watermelon. Only 16 KASP markers segregated in the interspecific C. amarus/lanatus validation population. Four markers showed significant differences in the separation of genotypes based on family-mean disease severity, but together explained only 16% of the phenotypic variance. Genotypes that inherited homozygous resistant parental alleles at three KASP markers had 42% lower family-mean disease severity than homozygous susceptible genotypes. Thus, haplotype analysis was more effective at predicting the mean disease severity of families than single markers. The haplotype identified in this study will be valuable for developing Fon race 2 resistant watermelon cultivars.

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.

Identification of WRKY Family Members and Characterization of the Low-Temperature-Stress-Responsive WRKY Genes in Luffa (Luffa cylindrica L.)

The plant-specific WRKY transcription factor family members have diverse regulatory effects on the genes associated with many plant processes. Although the WRKY proteins in Arabidopsis thaliana and other species have been thoroughly investigated, there has been relatively little research on the WRKY family in Luffa cylindrica, which is one of the most widely grown vegetables in China. In this study, we performed a genome-wide analysis to identify L. cylindrica WRKY genes, which were subsequently classified and examined in terms of their gene structures, chromosomal locations, promoter cis-acting elements, and responses to abiotic stress. A total of 62 LcWRKY genes (471-2238 bp) were identified and divided into three phylogenetic groups (I, II, and III), with group II further divided into five subgroups (IIa, IIb, IIc, IId, and IIe) in accordance with the classification in other plants. The LcWRKY genes were unevenly distributed across 13 chromosomes. The gene structure analysis indicated that the LcWRKY genes contained 0-11 introns (average of 4.4). Moreover, 20 motifs were detected in the LcWRKY proteins with conserved motifs among the different phylogenetic groups. Two subgroup IIc members (LcWRKY16 and LcWRKY31) contained the WRKY sequence variant WRKYGKK. Additionally, nine cis-acting elements related to diverse responses to environmental stimuli were identified in the LcWRKY promoters. The subcellular localization analysis indicated that three LcWRKY proteins (LcWRKY43, LcWRKY7, and LcWRKY23) are localized in the nucleus. The tissue-specific LcWRKY expression profiles reflected the diversity in LcWRKY expression. The RNA-seq data revealed the effects of low-temperature stress on LcWRKY expression. The cold-induced changes in expression were verified via a qRT-PCR analysis of 24 differentially expressed WRKY genes. Both LcWRKY7 and LcWRKY12 were highly responsive to the low-temperature treatment (approximately 110-fold increase in expression). Furthermore, the LcWRKY8, LcWRKY12, and LcWRKY59 expression levels increased by more than 25-fold under cold conditions. Our findings will help clarify the evolution of the luffa WRKY family while also providing valuable insights for future studies on WRKY functions.

Genome-Wide Identification of Vascular Plant One-Zinc-Finger Gene Family in Six Cucurbitaceae Species and the Role of CmoVOZ2 in Salt and Drought Stress Tolerance

As a plant-specific transcription factor, the vascular plant one-zinc-finger (VOZ) plays a crucial role in regulating various biological processes. In this study, a total of 17 VOZ genes in the Cucurbitaceae family were investigated using various bioinformatics methods. The 17 VOZ genes in Cucurbitaceae are distributed across 16 chromosomes. Based on the affinity of VOZ proteins to AtVOZ proteins, these 17 proteins were categorized into two groups: group I encompassed eight VOZ members, while group II comprised nine VOZ members. The expression profiles of CmoVOZs under various hormonal and abiotic stresses indicated that these genes were induced differentially by JA, ABA, GA, salt, and drought stress. Subsequently, CmoVOZ1 and CmoVOZ2 were found to be transcriptionally active, with the CmoVOZ2 protein being located mainly in the nucleus. Further experiments revealed that yeast cells expressing CmoVOZ2 gene showed increased tolerance to salt stress and drought stress. These results suggest that the VOZ gene family is not only important for plant growth and development but also that this mechanism may be universal across yeast and plants.

Association Analysis Provides Insights into Plant Mitonuclear Interactions

Cytonuclear interaction refers to the complex and ongoing process of coevolution between nuclear and organelle genomes, which are responsible for cellular respiration, photosynthesis, lipid metabolism, etc. and play a significant role in adaptation and speciation. There have been a large number of studies to detect signatures of cytonuclear interactions. However, identification of the specific nuclear and organelle genetic polymorphisms that are involved in these interactions within a species remains relatively rare. The recent surge in whole genome sequencing has provided us an opportunity to explore cytonuclear interaction from a population perspective. In this study, we analyzed a total of 3,439 genomes from 7 species to identify signals of cytonuclear interactions by association (linkage disequilibrium) analysis of variants in both the mitochondrial and nuclear genomes across flowering plants. We also investigated examples of nuclear loci identified based on these association signals using subcellular localization assays, gene editing, and transcriptome sequencing. Our study provides a novel perspective on the investigation of cytonuclear coevolution, thereby enriching our understanding of plant fitness and offspring sterility.

High-quality Momordica balsamina genome elucidates its potential use in improving stress resilience and therapeutic properties of bitter gourd

Introduction

Momordica balsamina is the closest wild species that can be crossed with an important fruit vegetable crop, Momordica charantia, has immense medicinal value, and placed under II subclass of primary gene pool of bitter gourd. M. balsamina is tolerant to major biotic and abiotic stresses. Genome characterization of Momordica balsamina as a wild relative of bitter gourd will contribute to the knowledge of the gene pool available for improvement in bitter gourd. There is potential to transfer gene/s related to biotic resistance and medicinal importance from M. balsamina to M. charantia to produce high-quality, better yielding and stress tolerant bitter gourd genotypes.

Methods

The present study provides the first and high-quality chromosome-level genome assembly of M. balsamina with size 384.90 Mb and N50 30.96 Mb using sequence data from 10x Genomics, Nanopore, and Hi-C platforms.

Results

A total of 6,32,098 transposons elements; 2,15,379 simple sequence repeats; 5,67,483 transcription factor binding sites; 3,376 noncoding RNA genes; and 41,652 protein-coding genes were identified, and 4,347 disease resistance, 67 heat stress-related, 05 carotenoid-related, 15 salt stress-related, 229 cucurbitacin-related, 19 terpenes-related, 37 antioxidant activity, and 06 sex determination-related genes were characterized.

Conclusion

Genome sequencing of M. balsamina will facilitate interspecific introgression of desirable traits. This information is cataloged in the form of webgenomic resource available at http://webtom.cabgrid.res.in/mbger/. Our finding of comparative genome analysis will be useful to get insights into the patterns and processes associated with genome evolution and to uncover functional regions of cucurbit genomes.

Analyzing genomic variation in cultivated pumpkins and identification of candidate genes controlling seed traits

Pumpkins are important vegetable crops widely grown worldwide, and seeds are considered a popular nutraceutical food and an excellent source of protein, oil, and vitamins. Seed size is one of the most important targets for commercial breeding in Cucurbita species; studies have shown that pumpkin seed size variation has a similar trend with fruit size, shape, and seed yield. However, few studies have been conducted to identify genetic loci controlling seed-related traits in cultivated pumpkins. This study analyzed the genomic characteristics of pumpkin breeding materials of 321 Cucurbita accessions collected worldwide, including Cucurbita moschata, Cucurbita maxima, and Cucurbita pepo, using extensive single nucleotide polymorphisms obtained from the genotyping-by-sequencing method, significant genetic variations were identified within and between Cucurbita species. Four major cultivar fruit types were further revealed in C. moschata species, and significant differentiation patterns were detected in several chromosomal regions. A total of 15 significant loci associated with pumpkin seed traits were mapped through a genome-wide association approach; 32 genes previously reported to be associated with seed size regulation in Arabidopsis and Oryza sativa were located in the intervals defined by linkage disequilibrium. Through this study, we gained a deep understanding of the genomic variation distribution across Cucurbita species. The available genetic resources and the associated genetic contents could be used in commercial pumpkin breeding and will facilitate molecular marker-assisted selection in pumpkin seed trait improvement.

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.

Comprehensive Profiling of Alternative Splicing and Alternative Polyadenylation during Fruit Ripening in Watermelon (Citrullus lanatus)

Fruit ripening is a highly complicated process that is accompanied by the formation of fruit quality. In recent years, a series of studies have demonstrated post-transcriptional control play important roles in fruit ripening and fruit quality formation. Till now, the post-transcriptional mechanisms for watermelon fruit ripening have not been comprehensively studied. In this study, we conducted PacBio single-molecule long-read sequencing to identify genome-wide alternative splicing (AS), alternative polyadenylation (APA) and long non-coding RNAs (lncRNAs) in watermelon fruit. In total, 6,921,295 error-corrected and mapped full-length non-chimeric (FLNC) reads were obtained. Notably, more than 42,285 distinct splicing isoforms were derived from 5,891,183 intron-containing full-length FLNC reads, including a large number of AS events associated with fruit ripening. In addition, we characterized 21,506 polyadenylation sites from 11,611 genes, 8703 of which have APA sites. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that fructose and mannose metabolism, starch and sucrose metabolism and carotenoid biosynthesis were both enriched in genes undergoing AS and APA. These results suggest that post-transcriptional regulation might potentially have a key role in regulation of fruit ripening in watermelon. Taken together, our comprehensive PacBio long-read sequencing results offer a valuable resource for watermelon research, and provide new insights into the molecular mechanisms underlying the complex regulatory networks of watermelon fruit ripening.

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 watermelon H3K4 and H3K27 genes and their expression profiles during watermelon fruit development

BACKGROUND

The ClaH3K4s and ClaH3K27s gene families are subfamilies of the SET family, each with a highly conserved SET structure domain and a PHD structural domain. Both participate in histone protein methylation, which affects the chromosome structure and gene expression, and is essential for fruit growth and development.

METHODS AND RESULTS

In order to demonstrate the structure and expression characteristics of ClaH3K4s and ClaH3K27s in watermelon, members of the watermelon H3K4 and H3K27 gene families were identified, and their chromosomal localization, gene structure, and protein structural domains were analyzed. The phylogeny and covariance of the gene families with other species were subsequently determined, and the expression profiles were obtained by performing RNA-Seq and qRT-PCR. The watermelon genome had five H3K4 genes with 3207-8043 bp nucleotide sequence lengths and four H3K27 genes with a 1107-5499 bp nucleotide sequence. Synteny analysis revealed the close relationship between watermelon and cucumber, with the majority of members displaying a one-to-one covariance. Approximately half of the 'Hua-Jing 13 watermelon' ClaH3K4s and ClaH3K27s genes were expressed more in the late fruit development stages, while the changes were minimal for the remaining half. H3K4-2 expression was observed to be slightly greater on day 21 compared to other periods. Moreover, ClaH3K27-1 and ClaH3K27-2 were hardly expressed throughout the developing period, and ClaH3K27-4 exhibited the highest expression.

CONCLUSION

These results serve as a basis for further functional characterization of the H3K4 and H3K27 genes in the fruit development of watermelon.

A Genome-Wide Identification and Comparative Analysis of the Heavy-Metal-Associated Gene Family in Cucurbitaceae Species and Their Role in Cucurbita pepo under Arsenic Stress

The heavy-metal-associated (HMA) proteins are a class of PB1-type ATPases related to the intracellular transport and detoxification of metals. However, due to a lack of information regarding the HMA gene family in the Cucurbitaceae family, a comprehensive genome-wide analysis of the HMA family was performed in ten Cucurbitaceae species: Citrullus amarus, Citrullus colocynthis, Citrullus lanatus, Citrullus mucosospermus, Cucumis melo, Cucumis sativus, Cucurbita maxima, Cucurbita moschata, Cucurbita pepo, and Legenaria siceraria. We identified 103 Cucurbit HMA proteins with various members, ranging from 8 (Legenaria siceraria) to 14 (Cucurbita pepo) across species. The phylogenetic and structural analysis confirmed that the Cucurbitaceae HMA protein family could be further classified into two major clades: Zn/Co/Cd/Pb and Cu/Ag. The GO-annotation-based subcellular localization analysis predicted that all HMA gene family members were localized on membranes. Moreover, the analysis of conserved motifs and gene structure (intron/exon) revealed the functional divergence between clades. The interspecies microsynteny analysis demonstrated that maximum orthologous genes were found between species of the Citrullus genera. Finally, nine candidate HMA genes were selected, and their expression analysis was carried out via qRT-PCR in root, leaf, flower, and fruit tissues of C. pepo under arsenic stress. The expression pattern of the CpeHMA genes showed a distinct pattern of expression in root and shoot tissues, with a remarkable expression of CpeHMA6 and CpeHMA3 genes from the Cu/Ag clade. Overall, this study provides insights into the functional analysis of the HMA gene family in Cucurbitaceae species and lays down the basic knowledge to explore the role and mechanism of the HMA gene family to cope with arsenic stress conditions.

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.

Fine-Scale analysis of both wild and cultivated horned galls provides insight into their quality differentiation

BACKGROUND

Galla chinensis is a traditional Chinese medicine (TCM) produced due to the interaction between the Fordinae aphids and the Rhus plant species. Horned galls with high tannin content are the most widely cultivated gall type, and Wufeng county of Hubei province in China is the center of cultivation. However, long-term artificial cultivation and domestication of horned galls to meet the increasing production demand have led to quality degradation. Understanding the reasons underlying quality degradation is urgent for horned gall production and application. The present study used a combination of metabolic, genetic, and ecological analyses to investigate the quality and genetic differentiation of the horned galls under long-term domestication as well as the potential relationships between them.

RESULTS

Analysis of gallic acid content and other three phenotypic traits (fresh weight, gall size, and wall thickness) revealed quality differentiation of horned galls collected from five locations in Wufeng, in which the cultivated samples from Wang Jiaping (WJP) showed the highest degradation. Genetic differentiation between the cultivated and wild Rhus chinensis trees in WJP, and between WJP and the other populations was detected based on SSR molecular markers, however, no significant difference in genetic structure was seen for the aphid populations. Among the various ecological factors examined, temperature was identified as the primary one affecting the quality of horned galls.

CONCLUSIONS

Both genetic and ecological factors caused quality differentiation of horned galls. The collection of diverse germplasm of host trees and aphids will help reduce the quality degradation of horned galls in Wufeng.

Comprehensive analyses of microtubule-associated protein MAP65 family genes in Cucurbitaceae and CsaMAP65s expression profiles in cucumber

MAP65 is a microtubule-binding protein family in plants and plays crucial roles in regulating cell growth and development, intercellular communication, and plant responses to various environmental stresses. However, MAP65s in Cucurbitaceae are still less understood. In this study, a total of 40 MAP65s were identified from six Cucurbitaceae species (Cucumis sativus L., Citrullus lanatus, Cucumis melo L., Cucurbita moschata, Lagenaria siceraria, and Benincasa hispida) and classified into five groups by phylogenetic analysis according to gene structures and conserved domains. A conserved domain (MAP65_ASE1) was found in all MAP65 proteins. In cucumber, we isolated six CsaMAP65s with different expression patterns in tissues including root, stem, leaf, female flower, male flower, and fruit. Subcellular localizations of CsaMAP65s verified that all CsaMAP65s were localized in microtubule and microfilament. Analyses of the promoter regions of CsaMAP65s have screened different cis-acting regulatory elements involved in growth and development and responses to hormone and stresses. In addition, CsaMAP65-5 in leaves was significantly upregulated by salt stress, and this promotion effect was higher in cucumber cultivars with salt tolerant than that without salt tolerant. CsaMAP65-1 in leaves was significantly upregulated by cold stress, and this promotion was higher in cold-tolerant cultivar than intolerant cultivar. With the genome-wide characterization and phylogenetic analysis of Cucurbitaceae MAP65s, and the expression profile of CsaMAP65s in cucumber, this study laid a foundation for further study on MAP65 functions in developmental processes and responses to abiotic stress in Cucurbitaceae species.

Allelic variations of ClACO gene improve nitrogen uptake via ethylene-mediated root architecture in watermelon

KEY MESSAGE

The ClACO gene encoding 1-aminocyclopropane-1-carboxylate oxidase enabled highly efficient 15N uptake in watermelon. Nitrogen is one of the most essential nutrient elements that play a pivotal role in regulating plant growth and development for crop productivity. Elucidating the genetic basis of high nitrogen uptake is the key to improve nitrogen use efficiency for sustainable agricultural productivity. Whereas previous researches on nitrogen absorption process are mainly focused on a few model plants or crops. To date, the causal genes that determine the efficient nitrogen uptake of watermelon have not been mapped and remains largely unknown. Here, we fine-mapped the 1-aminocyclopropane-1-carboxylate oxidase (ClACO) gene associated with nitrogen uptake efficiency in watermelon via bulked segregant analysis (BSA). The variations in the ClACO gene led to the changes of gene expression levels between two watermelon accessions with different nitrogen uptake efficiencies. Intriguingly, in terms of the transcript abundance of ClACO, it was concomitant with significant differences in ethylene evolutions in roots and root architectures between the two accessions and among the different genotypic offsprings of the recombinant BC2F1(ZJU132)-18. These findings suggest that ethylene as a negative regulator altered nitrogen uptake efficiency in watermelon by controlling root development. In conclusion, our current study will provide valuable target gene for precise breeding of 'green' watermelon varieties with high-nitrogen uptake efficiencies.

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.

The Progression in Developing Genomic Resources for Crop Improvement

Sequencing technologies have rapidly evolved over the past two decades, and new technologies are being continually developed and commercialized. The emerging sequencing technologies target generating more data with fewer inputs and at lower costs. This has also translated to an increase in the number and type of corresponding applications in genomics besides enhanced computational capacities (both hardware and software). Alongside the evolving DNA sequencing landscape, bioinformatics research teams have also evolved to accommodate the increasingly demanding techniques used to combine and interpret data, leading to many researchers moving from the lab to the computer. The rich history of DNA sequencing has paved the way for new insights and the development of new analysis methods. Understanding and learning from past technologies can help with the progress of future applications. This review focuses on the evolution of sequencing technologies, their significant enabling role in generating plant genome assemblies and downstream applications, and the parallel development of bioinformatics tools and skills, filling the gap in data analysis techniques.

Comprehensive Characterization of B-Box Zinc Finger Genes in Citrullus lanatus and Their Response to Hormone and Abiotic Stresses

Plant B-BOX (BBX) zinc finger transcription factors play crucial roles in growth and development and the stress response. Although the BBX family has been characterized in various plants, systematic analysis in watermelon is still lacking. In this study, 25 watermelon ClBBX genes were identified. ClBBXs were grouped into five clades (Clade I, II, III, IV, and V) based on their conserved domains and phylogenetic relationships. Most of the ClBBXs (84%) might be localized in the nuclei or cytoplasm. The classification of ClBBXs was consistent with their gene structures. They were unevenly distributed in nine chromosomes except for Chr4 and Chr10, with the largest number of six members in Chr2. Segmental duplications were the major factor in ClBBX family expansion. Some BBXs of watermelon and Arabidopsis evolved from a common ancestor. In total, 254 hormonal and stress-responsive cis elements were discovered in ClBBX promoters. ClBBXs were differentially expressed in tissues, and the expression levels of ClBBX15 and 16 were higher in aboveground tissues than in roots, while the patterns of ClBBX21a, 21b, 21c, 28 and 30b were the opposite. With salicylic acid, methyl jasmonate and salt stress conditions, 17, 18 and 18 ClBBXs exhibited significant expression changes, respectively. In addition, many ClBBXs, including ClBBX29b, 30a and 30b, were also responsive to cold and osmotic stress. In summary, the simultaneous response of multiple ClBBXs to hormonal or abiotic stress suggests that they may have functional interactions in the stress hormone network. Clarifying the roles of key ClBBXs in transcriptional regulation and mediating protein interactions will be an important task. Our comprehensive characterization of the watermelon ClBBX family provides vital clues for the in-depth analysis of their biological functions in stress and hormone signaling pathways.

Genome-wide characterization of long terminal repeat retrotransposons provides insights into trait evolution of four cucurbit species

Cucurbits are a diverse plant family that includes economically important crops, such as cucumber, watermelon, melon, and pumpkin. Knowledge of the roles that long terminal repeat retrotransposons (LTR-RTs) have played in diversification of cucurbit species is limited; to add to understanding of the roles of LTR-RTs, we assessed their distributions in four cucurbit species. We identified 381, 578, 1086, and 623 intact LTR-RTs in cucumber (Cucumis sativus L. var. sativus cv. Chinese Long), watermelon (Citrullus lanatus subsp. vulgaris cv. 97103), melon (Cucumis melo cv. DHL92), and Cucurbita (Cucurbita moschata var. Rifu), respectively. Among these LTR-RTs, the Ale clade of the Copia superfamily was the most abundant in all the four cucurbit species. Insertion time and copy number analysis revealed that an LTR-RT burst occurred approximately 2 million years ago in cucumber, watermelon, melon, and Cucurbita, and may have contributed to their genome size variation. Phylogenetic and nucleotide polymorphism analyses suggested that most LTR-RTs were formed after species diversification. Analysis of gene insertions by LTR-RTs revealed that the most frequent insertions were of Ale and Tekay and that genes related to dietary fiber synthesis were the most commonly affected by LTR-RTs in Cucurbita. These results increase our understanding of LTR-RTs and their roles in genome evolution and trait characterization in cucurbits.

Mutagenesis-based plant breeding approaches and genome engineering: A review focused on tomato

Breeding is the most important and efficient method for crop improvement involving repeated modification of the genetic makeup of a plant population over many generations. In this review, various accessible breeding approaches, such as conventional breeding and mutation breeding (physical and chemical mutagenesis and insertional mutagenesis), are discussed with respect to the actual impact of research on the economic improvement of tomato agriculture. Tomatoes are among the most economically important fruit crops consumed worldwide because of their high nutritional content and health-related benefits. Additionally, we summarize mutation-based mapping approaches, including Mutmap and MutChromeSeq, for the efficient mapping of several genes identified by random indel mutations that are beneficial for crop improvement. Difficulties and challenges in the adaptation of new genome editing techniques that provide opportunities to demonstrate precise mutations are also addressed. Lastly, this review focuses on various effective and convenient genome editing tools, such as RNA interference (RNAi), zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and clustered regularly interspaced short palindromic repeats (CRISPR/Cas9), and their potential for the improvement of numerous desirable traits to allow the development of better varieties of tomato and other horticultural crops.

Pan-genome of Citrullus genus highlights the extent of presence/absence variation during domestication and selection

The rich genetic diversity in Citrullus lanatus and the other six species in the Citrullus genus provides important sources in watermelon breeding. Here, we present the Citrullus genus pan-genome based on the 400 Citrullus genus resequencing data, showing that 477 Mb contigs and 6249 protein-coding genes were absent in the Citrullus lanatus reference genome. In the Citrullus genus pan-genome, there are a total of 8795 (30.5%) genes that exhibit presence/absence variations (PAVs). Presence/absence variation (PAV) analysis showed that a lot of gene PAV were selected during the domestication and improvement, such as 53 favorable genes and 40 unfavorable genes were identified during the C. mucosospermus to C. lanatus landrace domestication. We also identified 661 resistance gene analogs (RGAs) in the Citrullus genus pan-genome, which contains 90 RGAs (89 variable and 1 core gene) located on the pangenome additional contigs. By gene PAV-based GWAS, 8 gene presence/absence variations were found associated with flesh color. Finally, based on the results of gene PAV selection analysis between watermelon populations with different fruit colors, we identified four non-reference candidate genes associated with carotenoid accumulation, which had a significantly higher frequency in the white flesh. These results will provide an important source for watermelon breeding.

Allele mining through TILLING and EcoTILLING approaches in vegetable crops

MAIN CONCLUSION

The present review illustrates a comprehensive overview of the allele mining for genetic improvement in vegetable crops, and allele exploration methods and their utilization in various applications related to pre-breeding of economically important traits in vegetable crops. Vegetable crops have numerous wild descendants, ancestors and terrestrial races that could be exploited to develop high-yielding and climate-resilient varieties resistant/tolerant to biotic and abiotic stresses. To further boost the genetic potential of economic traits, the available genomic tools must be targeted and re-opened for exploitation of novel alleles from genetic stocks by the discovery of beneficial alleles from wild relatives and their introgression to cultivated types. This capability would be useful for giving plant breeders direct access to critical alleles that confer higher production, improve bioactive compounds, increase water and nutrient productivity as well as biotic and abiotic stress resilience. Allele mining is a new sophisticated technique for dissecting naturally occurring allelic variants in candidate genes that influence important traits which could be used for genetic improvement of vegetable crops. Target-induced local lesions in genomes (TILLINGs) is a sensitive mutation detection avenue in functional genomics, particularly wherein genome sequence information is limited or not available. Population exposure to chemical mutagens and the absence of selectivity lead to TILLING and EcoTILLING. EcoTILLING may lead to natural induction of SNPs and InDels. It is anticipated that as TILLING is used for vegetable crops improvement in the near future, indirect benefits will become apparent. Therefore, in this review we have highlighted the up-to-date information on allele mining for genetic enhancement in vegetable crops and methods of allele exploration and their use in pre-breeding for improvement of economic traits.

Creation of a watermelon haploid inducer line via ClDMP3-mediated single fertilization of the central cell

The use of doubled haploids is one of the most efficient breeding methods in modern agriculture. Irradiation of pollen grains has been shown to induce haploids in cucurbit crops, possibly because it causes preferential fertilization of the central cell over the egg cell. Disruption of the DMP gene is known to induce single fertilization of the central cell, which can lead to the formation of haploids. In the present study, a detailed method of creating a watermelon haploid inducer line via ClDMP3 mutation is described. The cldmp3 mutant induced haploids in multiple watermelon genotypes at rates of up to 1.12%. These haploids were confirmed via fluorescent markers, flow cytometry, molecular markers, and immuno-staining. The haploid inducer created by this method has the potential to greatly advance watermelon breeding in the future.

Characterization of HSP70 family in watermelon (Citrullus lanatus): identification, structure, evolution, and potential function in response to ABA, cold and drought stress

Watermelon (Citrullus lanatus) as a crop with important economic value, is widely cultivated around the world. The heat shock protein 70 (HSP70) family in plant is indispensable under stress conditions. However, no comprehensive analysis of watermelon HSP70 family is reported to date. In this study, 12 ClHSP70 genes were identified from watermelon, which were unevenly located in 7 out of 11 chromosomes and divided into three subfamilies. ClHSP70 proteins were predicted to be localized primarily in cytoplasm, chloroplast, and endoplasmic reticulum. Two pairs of segmental repeats and 1 pair of tandem repeats existed in ClHSP70 genes, and ClHSP70s underwent strong purification selection. There were many abscisic acid (ABA) and abiotic stress response elements in ClHSP70 promoters. Additionally, the transcriptional levels of ClHSP70s in roots, stems, true leaves, and cotyledons were also analyzed. Some of ClHSP70 genes were also strongly induced by ABA. Furthermore, ClHSP70s also had different degrees of response to drought and cold stress. The above data indicate that ClHSP70s may be participated in growth and development, signal transduction and abiotic stress response, laying a foundation for further analysis of the function of ClHSP70s in biological processes.

Lateral Root Initiation in Cucumber (Cucumis sativus): What Does the Expression Pattern of Rapid Alkalinization Factor 34 (RALF34) Tell Us?

In Arabidopsis, the small signaling peptide (peptide hormone) RALF34 is involved in the gene regulatory network of lateral root initiation. In this study, we aimed to understand the nature of the signals induced by RALF34 in the non-model plant cucumber (Cucumis sativus), where lateral root primordia are induced in the apical meristem of the parental root. The RALF family members of cucumber were identified using phylogenetic analysis. The sequence of events involved in the initiation and development of lateral root primordia in cucumber was examined in detail. To elucidate the role of the small signaling peptide CsRALF34 and its receptor CsTHESEUS1 in the initial stages of lateral root formation in the parental root meristem in cucumber, we studied the expression patterns of both genes, as well as the localization and transport of the CsRALF34 peptide. CsRALF34 is expressed in all plant organs. CsRALF34 seems to differ from AtRALF34 in that its expression is not regulated by auxin. The expression of AtRALF34, as well as CsRALF34, is regulated in part by ethylene. CsTHESEUS1 is expressed constitutively in cucumber root tissues. Our data suggest that CsRALF34 acts in a non-cell-autonomous manner and is not involved in lateral root initiation in cucumber.

Comparative chloroplast genome analysis of seven extant Citrullus species insight into genetic variation, phylogenetic relationships, and selective pressure

Citrullus ecirrhosus, Citrullus rehmii, and Citrullus naudinianus are three important related wild species of watermelon in the genus Citrullus, and their morphological differences are clear, however, their chloroplast genome differences remain unknown. This study is the first to assemble, analyze, and publish the complete chloroplast genomes of C. ecirrhosus, C. rehmii, and C. naudinianus. A comparative analysis was then conducted among the complete chloroplast genomes of seven extant Citrullus species, and the results demonstrated that the average genome sizes of Citrullus is 157,005 bp, a total of 130-133 annotated genes were identified, including 8 rRNA, 37 tRNA and 85-88 protein-encoding genes. Their gene content, order, and genome structure were similar. However, noncoding regions were more divergent than coding regions, and rps16-trnQ was a hypervariable fragment. Thirty-four polymorphic SSRs, 1,271 SNPs and 234 INDELs were identified. Phylogenetic trees revealed a clear phylogenetic relationship of Citrullus species, and the developed molecular markers (SNPs and rps16-trnQ) could be used for taxonomy in Citrullus. Three genes (atpB, clpP1, and rpoC2) were identified to undergo selection and would promote the environmental adaptation of Citrullus.

DNA methylation signatures of duplicate gene evolution in angiosperms

Gene duplication is a source of evolutionary novelty. DNA methylation may play a role in the evolution of duplicate genes (paralogs) through its association with gene expression. While this relationship has been examined to varying extents in a few individual species, the generalizability of these results at either a broad phylogenetic scale with species of differing duplication histories or across a population remains unknown. We applied a comparative epigenomics approach to 43 angiosperm species across the phylogeny and a population of 928 Arabidopsis (Arabidopsis thaliana) accessions, examining the association of DNA methylation with paralog evolution. Genic DNA methylation was differentially associated with duplication type, the age of duplication, sequence evolution, and gene expression. Whole genome duplicates were typically enriched for CG-only gene-body methylated or unmethylated genes, while single-gene duplications were typically enriched for non-CG methylated or unmethylated genes. Non-CG methylation, in particular, was characteristic of more recent single-gene duplicates. Core angiosperm gene families differentiated into those which preferentially retain paralogs and 'duplication-resistant' families, which convergently reverted to singletons following duplication. Duplication-resistant families that still have paralogous copies were, uncharacteristically for core angiosperm genes, enriched for non-CG methylation. Non-CG methylated paralogs had higher rates of sequence evolution, higher frequency of presence-absence variation, and more limited expression. This suggests that silencing by non-CG methylation may be important to maintaining dosage following duplication and be a precursor to fractionation. Our results indicate that genic methylation marks differing evolutionary trajectories and fates between paralogous genes and have a role in maintaining dosage following duplication.

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.

Comparative bioinformatics analysis and abiotic stress responses of expansin proteins in Cucurbitaceae members: watermelon and melon

Watermelon and melon are members of the Cucurbitaceae family including economically significant crops in the world. The expansin protein family, which is one of the members of the cell wall, breaks down the non-covalent bonds between cell wall polysaccharides, causing pressure-dependent cell expansion. Comparative bioinformatics and molecular characterization analysis of the expansin protein family were carried out in the watermelon (Citrullus lanatus) and melon (Cucumis melo) plants in the study. Gene expression levels of expansin family members were analyzed in leaf and root tissues of watermelon and melon under ABA, drought, heat, cold, and salt stress conditions by quantitative real-time PCR analysis. After comprehensive searches, 40 expansin proteins (22 ClaEXPA, 14 ClaEXPLA, and 4 ClaEXPB) in watermelon and 43 expansin proteins (19 CmEXPA, 15 CmEXPLA, 3 CmEXPB, and 6 CmEXPLB) in melon were identified. The greatest orthologous genes were identified with soybean expansin genes for watermelon and melon. However, the latest divergence time between orthologous genes was determined with poplar expansin genes for watermelon and melon expansin genes. ClaEXPA-04, ClaEXPA-09, ClaEXPB-01, ClaEXPB-03, and ClaEXPLA-13 genes in watermelon and CmEXPA-12, CmEXPA-10, and CmEXPLA-01 genes in melon can be involved in tissue development and abiotic stress response of the plant. The current study combining bioinformatics and experimental analysis can provide a detailed characterization of the expansin superfamily which has roles in growth and reaction to the stress of the plant. The study ensures detailed data for future studies examining gene functions including the roles in plant growth and stress conditions.

Comparative karyotype analysis of eight Cucurbitaceae crops using fluorochrome banding and 45S rDNA-FISH

To have an insight into the karyotype variation of eight Cucurbitaceae crops including Cucumissativus Linnaeus, 1753, Cucumismelo Linnaeus, 1753, Citrulluslanatus (Thunberg, 1794) Matsumura et Nakai, 1916, Benincasahispida (Thunberg, 1784) Cogniaux, 1881, Momordicacharantia Linnaeus, 1753, Luffacylindrica (Linnaeus, 1753) Roemer, 1846, Lagenariasicerariavar.hispida (Thunberg, 1783) Hara, 1948 and Cucurbitamoschata Duchesne ex Poiret, 1819, well morphologically differentiated mitotic metaphase chromosomes were prepared using the enzymatic maceration and flame-drying method, and the chromosomal distribution of heterochromatin and 18S-5.8S-26S rRNA genes (45S rDNA) was investigated using sequential combined PI and DAPI (CPD) staining and fluorescence in situ hybridization (FISH) with 45S rDNA probe. Detailed karyotypes were established using the dataset of chromosome measurements, fluorochrome bands and rDNA FISH signals. Four karyotype asymmetry indices, CVCI, CVCL, MCA and Stebbins' category, were measured to elucidate the karyological relationships among species. All the species studied had symmetrical karyotypes composed of metacentric and submetacentric or only metacentric chromosomes, but their karyotype structure can be discriminated by the scatter plot of MCA vs. CVCL. The karyological relationships among these species revealed by PCoA based on x, 2n, TCL, MCA, CVCL and CVCI was basically in agreement with the phylogenetic relationships revealed by DNA sequences. CPD staining revealed all 45S rDNA sites in all species, (peri)centromeric GC-rich heterochromatin in C.sativus, C.melo, C.lanatus, M.charantia and L.cylindrica, terminal GC-rich heterochromatin in C.sativus. DAPI counterstaining after FISH revealed pericentromeric DAPI+ heterochromatin in C.moschata. rDNA FISH detected two 45S loci in five species and five 45S loci in three species. Among these 45S loci, most were located at the terminals of chromosome arms, and a few in the proximal regions. In C.sativus, individual chromosomes can be precisely distinguished by the CPD band and 45S rDNA signal patterns, providing an easy method for chromosome identification of cucumber. The genome differentiation among these species was discussed in terms of genome size, heterochromatin, 45S rDNA site, and karyotype asymmetry based on the data of this study and previous reports.

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.

Drought stress tolerance mechanisms and their potential common indicators to salinity, insights from the wild watermelon (Citrullus lanatus): A review

Climate change has escalated the effect of drought on crop production as it has negatively altered the environmental condition. Wild watermelon grows abundantly in the Kgalagadi desert even though the environment is characterized by minimal rainfall, high temperatures and intense sunshine during growing season. This area is also characterized by sandy soils with low water holding capacity, thus bringing about drought stress. Drought stress affects crop productivity through its effects on development and physiological functions as dictated by molecular responses. Not only one or two physiological process or genes are responsible for drought tolerance, but a combination of various factors do work together to aid crop tolerance mechanism. Various studies have shown that wild watermelon possess superior qualities that aid its survival in unfavorable conditions. These mechanisms include resilient root growth, timely stomatal closure, chlorophyll fluorescence quenching under water deficit as key physiological responses. At biochemical and molecular level, the crop responds through citrulline accumulation and expression of genes associated with drought tolerance in this species and other plants. Previous salinity stress studies involving other plants have identified citrulline accumulation and expression of some of these genes (chloroplast APX, Type-2 metallothionein), to be associated with tolerance. Emerging evidence indicates that the upstream of functional genes are the transcription factor that regulates drought and salinity stress responses as well as adaptation. In this review we discuss the drought tolerance mechanisms in watermelons and some of its common indicators to salinity at physiological, biochemical and molecular level.

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.

Chromosome-Level Genome Assembly of Herpetospermum pedunculosum (Cucurbitaceae)

This study presents a chromosome-level reference genome assembly of a traditional Tibetan medicinal plant, Herpetospermum pedunculosum belonging to the Cucurbitaceae family. Following a combined PacBio high-fidelity sequencing and Hi-C analysis, a final H. pedunculosum genome assembly, 804.11 Mb in length was obtained, 90.45% of which was anchored into ten pseudochromosomes with a contig N50 of 24.39 Mb. In addition, 579.55 Mb repetitive sequences and 23,924 high-confidence protein-coding genes were annotated. Phylogenetic analysis revealed that H. pedunculosum was sister to a clade formed by cucumber, zucchini, and wax gourd. Further whole-genome duplication analysis revealed no recent polyploidization event in the H. pedunculosum genome. The high-quality H. pedunculosum genome presented here will be highly useful in investigating the molecular mechanisms underlying the biosynthesis of its active compounds and adaptation strategies to the extreme environment. It will also provide great insights into comparative genomic studies of Cucurbitaceae and flowering plants.

Synergistic effect of Bacillus subtilis and Paecilomyces lilacinus in alleviating soil degradation and improving watermelon yield

Continuous cropping of watermelon (Citrullus lanatus) may lead to soil degradation. As a soil conditioner, microbial agent has great potential in improving soil function and enhancing plant growth. In this study, we aimed to explore how microbial agent relieves the soil sickness of watermelon by analyzing watermelon performance, soil physicochemical properties and microbial community structures. Results suggested that microbial agent treatments significantly changed the photosynthetic efficiency of upper and lower leaves, which helped improve the growth of watermelon. The single fruit weight, fruit sugar degree and total phosphorus of soil following treatment with a mixture of Paecilomyces lilacinus DZ910 and Bacillus subtilis KC1723 (treatment D_K) were higher than those in single biofertilizer treatments and control. The soil microbial community under microbial agent treatments also changed significantly, indicating the feasibility of using microbial agents as soil remediations. The proportions of Pseudomonas and Flavobacterium, changed significantly after using microbial agents. Pseudomonas increased significantly after B. subtilis KC1723 and D_K treatments, while Flavobacterium increased significantly after using all three kinds of microbial agents compared to control. Increases in these bacteria were positively correlated with agronomic variables of watermelon. The fungi Aspergillus and Neocosmospora in the soil, which create an soil sickness of watermelon, decreased after KC1723 and D_K treatments. Meanwhile, Aspergillus and Neocosmospora were positively related to Myceliophthora incidence and negatively correlated with watermelon growth (single fruit weight and photosynthetic efficiency of upper leaves). Our microbial agent, especially D_K, represents a useful technique for alleviating soil sickness in watermelon.

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.

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.