DNA methylation is an important epigenetic mark associated with plant immunity, but little is known about its roles in viral infection of watermelon. We carried out whole-genome bisulfite sequencing of watermelon leaves at 0 h (ck), 48 h, and 25 days post-inoculation with Cucumber green mottle mosaic virus (CGMMV). The number of differentially methylated regions (DMRs) increased during CGMMV infection and 2788 DMR-associated genes (DMGs) were screened out among three libraries. Most DMRs and DMGs were obtained under the CHH context. These DMGs were significantly enriched in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of secondary biosynthesis and metabolism, plant–pathogen interactions, Toll-like receptor signaling, and ABC transporters. Additionally, DMGs encoding PR1a, CaMs, calcium-binding protein, RIN4, BAK1, WRKYs, RBOHs, STKs, and RLPs/RLKs were involved in the watermelon–CGMMV interaction and signaling. The association between DNA methylation and gene expression was analyzed by RNA-seq and no clear relationship was detected. Moreover, downregulation of genes in the RdDM pathway suggested the reduced RdDM-directed CHH methylation plays an important role in antiviral defense in watermelon. Our findings provide genome-wide DNA methylation profiles of watermelon and will aid in revealing the molecular mechanism in response to CGMMV infection at the methylation level.

Watermelon (Citrullus lanatus L.), which is an economically important cucurbit crop that is cultivated worldwide, is vulnerable to various adverse environmental conditions. Small heat shock protein 20s (HSP20s) are the most abundant plant HSPs and they play important roles in various biotic and abiotic stress responses. However, they have not been systematically investigated in watermelon. In this study, we identified 44 watermelon HSP20 genes and analyzed their gene structures, conserved domains, phylogenetic relationships, chromosomal distributions, and expression profiles. All of the watermelon HSP20 proteins have a conserved the α-crystallin (ACD) domain. Half of the ClHSP20s arose through gene duplication events. Plant HSP20s were grouped into 18 subfamiles and a new subfamily, nucleo-cytoplasmic XIII (CXIII), was identified in this study. Numerous stress- and hormone-responsive cis-elements were detected in the putative promoter regions of the watermelon HSP20 genes. Different from that in other species, half of the watermelon HSP20s were repressed by heat stress. Plant HSP20s displayed diverse responses to different virus infections and most of the ClHSP20s were generally repressed by Cucumber green mottle mosaic virus (CGMMV). Some ClHSP20s exhibited similar transcriptional responses to abscisic acid, melatonin, and CGMMV. Subcellular localization analyses of six selected HSP20- green fluorescence protein fusion proteins revealed diverse subcellular targeting. Some ClHSP20 proteins were affected by CGMMV, as reflected by changes in the size, number, and distribution of fluorescent granules. These systematic analyses provide a foundation for elucidating the physiological functions and biological roles of the watermelon HSP20 gene family.

Cucumber green mottle mosaic virus (CGMMV) is a member of the genus Tobamovirus, which cause diseases in cucurbits, especially watermelon. In watermelon, symptoms develop on the whole plant, including leaves, stems, peduncles, and fruit. To better understand the molecular mechanisms of watermelon early responses to CGMMV infection, a comparative transcriptome analysis of 24 h CGMMV-infected and mock-inoculated watermelon leaves was performed. A total of 1641 differently expressed genes (DEGs) were identified, with 886 DEGs upregulated and 755 DEGs downregulated after CGMMV infection. A functional analysis indicated that the DEGs were involved in photosynthesis, plant–pathogen interactions, secondary metabolism, and plant hormone signal transduction. In addition, a few transcription factor families, including WRKY, MYB, HLH, bZIP and NAC, were responsive to the CGMMV-induced stress. To confirm the high-throughput sequencing results, 15 DEGs were validated by qRT-PCR analysis. The results provide insights into the identification of candidate genes or pathways involved in the responses of watermelon leaves to CGMMV infection.