Mechanisms governing melon fruit skin pigmentation: Insights from transcriptome sequencing and whole-genome bisulfite sequencing analyses

Fruit skin color is a critical agronomic trait determining fruit quality, but the mechanisms governing melon (Cucumis melo L.) fruit skin pigmentation remain inadequately characterized. We performed cytological, transcriptomic, and whole-genome bisulfite sequencing (WGBS) comparative analyses on two near-isogenic lines (NILs) at different developmental stages: the green-skinned (GS) and white-skinned (WS) melon. The enriched functions among the 536 identified differentially expressed genes (DEGs) were related to chloroplast development, pigment metabolism, and photosynthetic pathway. The genes in these pathways were downregulated in the WS line, potentially contributing to the decreased chlorophyll content and abnormal chloroplast development in WS. WGBS analysis revealed that the levels of genomic DNA methylation in WS progressively exceeded those in GS as the fruits developed. The integrated analysis of the transcriptome and methylome identified DMR-associated DEGs (CHLI, CRD1, PORA, and HCAR) involved in chlorophyll metabolism. Notably, PORA exhibited downregulated expression and was hypermethylated in WS fruits during two developmental stages, showing a coordinated expression with the chlorophyll contents in developing WS and GS fruits. Furthermore, we identified a transcription factor gene, MELO3C011576, that was closely associated with the expression of PORA and chlorophyll levels across three developmental stages in both NILs. We speculate that PORA positively regulates chlorophyll biosynthesis in the two NILs, with its expression may be co-controlled by DNA methylation and transcription factors.