Background: Fruit skin color play important role in commercial value of cucumber, which is mainly determined by the content and composition of chlorophyll and anthocyanins. Therefore, understanding the related genes and metabolomics involved in composition of fruit skin color is essential for cucumber quality and commodity value. Results: The results showed that chlorophyll a, chlorophyll b and carotenoid content in fruit skin were higher in Lv (dark green skin) than Bai (light green skin) on fruit skin. Cytological observation showed more chloroplast existed in fruit skin cells of Lv. A total of 162 significantly different metabolites were found between the fruit skin of the two genotypes by metabolome analysis, including 40 flavones, 9 flavanones, 8 flavonols, 6 anthocyanins, and other compounds. Crucial anthocyanins and flavonols for fruit skin color, were detected significantly decreased in fruit skin of Bai compared with Lv. By RNA-seq assay, 4516 differentially expressed genes (DEGs) were identified between two cultivars. Further analyses suggested that low expression level of chlorophyll biosynthetic genes, such as chlM, por and NOL caused less chlorophylls or chloroplast in fruit skin of Bai. Meanwhile, a predicted regulatory network of anthocyanin biosynthesis was established to illustrate involving many DEGs, especially 4CL, CHS and UFGT. Conclusions: This study uncovered significant differences between two cucumber genotypes with different fruit color using metabolome and RNA-seq analysis. We lay a foundation to understand molecular regulation mechanism on formation of cucumber skin color, by exploring valuable genes, which is helpful for cucumber breeding and improvement on fruit skin color.

Drought stress is one of the most serious threats to cucumber quality and yield. To gain a good understanding of the molecular mechanism upon water deficiency, we compared and analyzed the RNA sequencing-based transcriptomic responses of two contrasting cucumber genotypes, L-9 (drought-tolerant) and A-16 (drought-sensitive). In our present study, combining the analysis of phenotype, twelve samples of cucumber were carried out a transcriptomic profile by RNA-Seq under normal and water-deficiency conditions, respectively. A total of 1008 transcripts were differentially expressed under normal conditions (466 up-regulated and 542 down-regulated) and 2265 transcripts under drought stress (979 up-regulated and 1286 down-regulated). The significant positive correlation between RNA sequencing data and a qRT-PCR analysis supported the results found. Differentially expressed genes (DEGs) involved in metabolic pathway and biosynthesis of secondary metabolism were significantly changed after drought stress. Several genes, which were related to sucrose biosynthesis (Csa3G784370 and Csa3G149890) and abscisic acid (ABA) signal transduction (Csa4M361820 and Csa6M382950), were specifically induced after 4 days of drought stress. DEGs between the two contrasting cultivars identified in our study provide a novel insight into isolating helpful candidate genes for drought tolerance in cucumber.

Temperature rising caused by global warming has imposed significant negative effects on crop qualities and yields. To get the well-known molecular mechanism upon the higher temperature, we carefully analyzed the RNA sequencing-based transcriptomic responses of two contrasting chieh-qua genotypes: A39 (heat-tolerant) and H5 (heat-sensitive). In this study, twelve cDNA libraries generated from A39 and H5 were performed with a transcriptome assay under normal and heat stress conditions, respectively. A total of 8705 differentially expressed genes (DEGs) were detected under normal conditions (3676 up-regulated and 5029 down-regulated) and 1505 genes under heat stress (914 up-regulated and 591 down-regulated), respectively. A significant positive correlation between RNA-Seq data and qRT-PCR results was identified. DEGs related to heat shock proteins (HSPs), ubiquitin-protein ligase, transcriptional factors, and pentatricopeptide repeat-containing proteins were significantly changed after heat stress. Several genes, which encoded HSPs (CL2311.Contig3 and CL6612.Contig2), cytochrome P450 (CL4517.Contig4 and CL683.Contig7), and bHLH TFs (CL914.Contig2 and CL8321.Contig1) were specifically induced after four days of heat stress. DEGs detected in our study between these two contrasting cultivars would provide a novel basis for isolating useful candidate genes of heat stress responses 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.

Bottle gourd (Lagenaria siceraria (Molina) Standl) is a widely distributed Cucurbitaceae species, but gaps and low-quality assemblies have limited its genomic study. To address this, we assembled a nearly complete, high-quality genome of the bottle gourd (Pugua) using PacBio HiFi sequencing and Hi-C correction. The genome, being 298.67 Mb long with a ContigN50 of 28.55 Mb, was identified to possess 11 chromosomes, 11 centromeres, 18 telomeres, and 24,439 predicted protein-coding genes; notably, gap-free telomere-to-telomere assembly was accomplished for seven chromosomes. Based on the Pugua genome, the transcriptomic and metabolomic combined analyses revealed that amino acids and lipids accumulate during the expansion stage, while sugars and terpenoids increase during ripening. GA4 and genes of the Aux/IAA family mediate fruit expansion and maturation, while cell wall remodeling is regulated by factors such as XTHs, EXPs, polyphenols, and alkaloids, contributing to environmental adaptation. GGAT2 was positively correlated with glutamate, a source of umami, and SUS5 and SPS4 expression aligned with sucrose accumulation. This study provides a valuable genetic resource for bottle gourd research, enhancing the understanding of Cucurbitaceae evolution and supporting further studies on bottle gourd development, quality, and genetic improvement.