- Research
- Open access
- Published:
Phytoalexin deficient 4 is associated with the lesion mimic trait in watermelon clalm mutant (Citrullus lanatus)
BMC Plant Biology volume 25, Article number: 92 (2025)
Abstract
In watermelon (Citrullus lanatus), lesion mimic is a rare, valuable trait that can be used by breeders for selection at early growth stages. In this study, we tested a seven-generation family to determine the inheritance and genetic basis of this trait. As revealed by analysis of the lesion mimic mutant clalm, this trait is controlled by a single dominant gene. Whole genome resequencing–bulked segregant analysis demonstrated that this gene is located on chromosome 4 from 3,760,000 bp to 7,440,000 bp, a region corresponding to a physical distance of 3.68 Mb encompassing approximately 72 annotated genes. There are 6 genes with non-synonymous mutation SNP sites. The predicted target gene, ClCG04G001930, encodes a Phytoalexin deficient 4 (PAD4), a protein that plays an important regulatory role in leaf senescence in many plant species. According to quantitative real-time PCR analysis, the expression level of ClCG04G001930 was significantly higher in the clalm mutant than in normal watermelon. Twenty-five SNPs were identified in the ClCG04G001930 gene of F2 individuals of the clalm mutant. Overexpression the ClCG04G001930 gene, designated as ClPAD4, yielded transgenic lines whose leaves gradually developed chlorotic lesions over 3 weeks. RNA interference of the ClPAD4 yielded transgenic lines whose cotyledon prone to diseased over 2 weeks. Our results suggest that ClPAD4 might be the candidate gene responsible for lesion mimic in the clalm mutant. Our findings may serve as a foundation for elucidating the mechanism underlying the molecular metabolism of programmed cell death and should be useful for marker-assisted selection breeding in watermelon.
Introduction
Plant lesion mimic mutants (LMMs) exhibit spontaneous cell death in the absence of pathogen attack [1]. The LMMs of rice, such as spl7, spl11, spl28, spl30 and lmr, exhibit small, reddish-brown lesions or spots that develop at leaf tips or bases and gradually spread on fully expanded leaves [2,3,4,5]. In Arabidopsis, overexpression of UGT76D1, which plays an important role in SA homeostasis, results in a hypersensitive response (HR)-like lesion mimic phenotype [6]. In birch, BpGH3.5, an early auxin-response factor regulating root elongation, gives rise to typical LMM characteristics and accelerates leaf senescence [7]. In switchgrass, co-silencing of methylenetetrahydrofolate reductase and caffeic acid O-methyltransferase results in a novel lesion-mimic leaf phenotype [8]. In contrast to these findings, the genetic mechanism underlying the lesion-mimic leaf phenotype of watermelon (Citrullus lanatus) remains unknown.
The genome of watermelon (2n = 2x = 22) comprises 11 chromosomes and has an approximate size of 425 Mb. The annotated watermelon genome sequence, which was released in 2013, includes 23,440 predicted protein-coding genes and is thus suitable for next-generation sequencing to accelerate the identification of candidate genes controlling important agronomic traits in this species. Watermelon fruit shape, skin color, and rind stripe patterning have been previously investigated by bulked segregant analysis and next-generation sequencing technology (BSA-seq) [9,10,11]. Many dwarfing genes have been identified through BSA-seq. Two single-nucleotide polymorphism (SNP) mutation of GA20ox, one SNP mutation of GA3ox and a SNP deletion in an ABC transporter gene lead to dwarfism phenotypes in different watermelon [12,13,14,15]. The SNP-based genetic linkage map was constructed using a recombinant inbred lines derived from a cross between the 97103 (Citrullus lanatus var lanatus) and the 296341-FR (Citrullus lanatus var citroides) [16,17,18]. The objective of the present study was the identification of the LMM gene in the watermelon clalm mutant. We re-sequenced the whole genome of two DNA bulks (i.e., spotted leaf and normal pools) developed from plants in an F2 population. Our study results provide preliminary evidence that ClCG04G001930 encoding phytoalexin deficient 4 (PAD4) might be the candidate gene responsible for. To our knowledge, our study is the first reported gene mapping and elucidation of the inheritance mechanism controlling lesion mimic in watermelon. To improve the quality of watermelon cultivars, breeders have focused on the introduction of novel traits into existing germplasm. The application of the spotted leaf trait as a selection tool in breeding programs will thus help improve the ability of breeders to make selections at early growth stages, thus accelerating the watermelon breeding program. In practical production and application, the seedlings with spots should be retained, and the seedlings without spots should be removed. The purity of watermelon can be guaranteed by observation, which is simple and fast.
Methods
Plant materials
The clalm mutant (lesion mimic mutants of Citrullus lanatus var lanatus) was detected in an inbred watermelon line derived from ‘V603’ (Citrullus lanatus var lanatus) in July 2012, after years of screening and cultivating. The clalm mutant were analyzed by conducting testcross experiments. The clalm mutant were hand-pollinated with pollen from four inbred lines ‘Kexi’, ‘Hei2’, ‘Cai10’, and ‘K1’ (Citrullus lanatus var lanatus) to generate the F1 generation in the spring of 2013. The wild type watermelon cultivars ‘Zd’ (Citrullus lanatus var lanatus) is male parent with no spotting phenotypes. The clalm mutant is female parent with yellow spots spreading over leaves and fruits. Watermelon cultivars ‘Zd’ and clalm mutant were self-pollinated for seven generations to obtain stable phenotypes. F2 population obtained by clalm and ‘Zd’ crosses. For use as BSA-seq research materials, 60 F2 segregation populations (30 + 30 mixed pools with extreme characteristics), 10 parent plants clalm and ‘Zd’ were grown and evaluated at the Henan University Genetics and Breeding Base in the spring of 2018.
Trypan blue staining
To analyze the state of lesion mimic leaf cells, leaves were incubated in fixative solution, immersed and rinsed five times in ultrapure water, and placed on filter paper. After absorption of excess water, the leaves were stained at room temperature for 2 to 6Â h in darkness. Next, the leaves were immersed and rinsed five time in ultrapure water and then immersed and stored in plant trypan blue dye solution B for 3 to 16Â h. After incubation for 30Â min at room temperature in 20Â ml of plant trypan blue dye solution, the leaves were photographed.
Microscopic observations
For light microscopic examination, leaves were fixed in fixative solution (4% paraformaldehyde) for 24 h and then dehydrated, embedded in paraffin, and sliced. The slices were stained with toluidine blue and sealed in neutral gum heated to 38 °C. The prepared samples were observed with an Eclipse E100 microscope (Nikon, Tokyo, Japan) at 100 × magnification. For staining prior to transmission electron microscopy (TEM) examination, targeted fresh tissues were selected to minimize mechanical damage. Fresh tissue blocks (≤ 1 mm2) were quickly cut with a sharp blade and harvested within 1 to 3 min. The washed tissue blocks were immediately fixed with electron microscopy fixative for 2 h at room temperature and then transferred to 4 °C for preservation and transportation.
Samples were fixed in a solution of 1% OsO4 in PB (0.1 M, pH 7.4) for 7 h at room temperature. After removal of OsO4, the tissues were rinsed three times in 0.1 M PB (pH 7.4), dehydrated in a graded ethanol series (30%, 50%, 70%, 80%, 90%, 95%, and 100%) for 1 h at each concentration, resin penetrated, and embedded at 37 °C in a graded reagent series (3:1 acetone:EMBed 812 for 2–4 h; 1:1 acetone:EMBed 812 overnight; 1:3 acetone:EMBed 812 for 2–4 h; and pure EMBed 812 for 5–8 h). Tissues were inserted into embedding molds filled with pure EMBed 812 resin and then incubated overnight in a 37 °C oven. The embedding molds containing resin and samples were moved to a 65 °C oven and allowed to polymerize for more than 48 h. The resin blocks were cut into 60- to 80-nm sections using an ultramicrotome, and the tissues were fished out onto 150-mesh cuprum grids with formvar film. The grid-attached tissues were stained with 2% uranium acetate-saturated alcohol solution in darkness for 8 min, rinsed three times each with 70% ethanol and ultrapure water, stained with 2.6% lead citrate in the absence of CO2 for 8 min, and rinsed three times with ultrapure water. After blotting with filter paper, the cuprum grids were placed on a grid board and dried overnight at room temperature. The samples were examined under an HT7700 transmission electron microscope (Hitachi, Japan).
Sample collection and BSA library preparation
Young leaves of the parents and F2 plants were collected and immediately placed in liquid nitrogen, and total genomic DNA was isolated using the CTAB method. Genomic DNA samples of 60 F2 segregation populations and 10 parent plants were subjected to whole-genome resequencing performed by the Beijing Biomarker Technologies Corporation (Beijing, China). To generate bulked samples, equal amounts of DNA from each plant per group were mixed to form yellow-spot (Y-pool) and normal (N-pool) sets at a final concentration of 40 ng/μL. Prior to high-throughput sequencing, DNA samples were sonicated to produce 350-bp fragments. After trimming of barcodes, clean, high-quality reads from identical samples were mapped onto the C. lanatus subsp. vulgaris ‘Charleston Gray’ genome sequence (http://cucurbitgenomics.org/organism/4). The reason for selecting C. lanatus subsp. vulgaris ‘Charleston Gray’ genome sequence as the reference genome is that the character of clalm mutant is similar to that. All identified SNPs shared across the bulk were considered polymorphic in association studies. To represent the difference between the SNP index of the two pools, we calculated ΔSNP and ΔInDel indexes, which are association analysis metrics used to find significant differences in genotype frequency between two pools [19, 20]. Candidate regions over the threshold (99th percentile) were extracted from each linkage group.
Expression analysis of candidate genes
We investigated the expression patterns of ClCG04G000300, ClCG04G000420, ClCG04G001450, ClCG04G001740, ClCG04G001900, and ClCG04G001930 using qRT-PCR. 'Zd' and clalm mutant seedlings were grown for 10, 30, and 60 days. Spotted leaves and normal leaves were collected from three plants each, and each sample collection was repeated three times. After addition of SYBR Green I to each reaction mixture, qRT-PCR amplifications were performed on a Roche LightCycler 480 II instrument using the following thermal profile: a 30-s hot start at 95 °C followed by 40 cycles of 95 °C for 5 s and 60 °C for 30 s. The ClYLS8 gene was included as a control for normalizing gene expression data [12]. The primers used to amplify genes are listed in Table 1.
Bioinformatic analysis
The homology of nucleotide sequences of candidate gene was analyzed online (http://www.ncbi.nlm.nih.gov). A phylogenetic tree comprising the target gene and its homologous sequences was constructed using the MEGA 7.0 software [21]. A bootstrap analysis (1,000 replicates) was completed to assess the reliability of the tree [22]. The high quality tree figure was created online (https://itol.embl.de/). The protein tertiary structure of ClPAD4 was predicted online ( https://swissmodel.expasy.org/) [23].
Watermelon transformation
The binary vectors pCambia 1300 were constructed. The candidate gene was inserted between SalI and SpeI sites in the vector. The expression vector pCAM1300 were constructed. The candidate gene was inserted between SalI and SpeI sites in the vector. The binary RNAi interference vectors DHpart27RNAi-FADP1P4 were constructed as described by Dong et al. (2021) [24]. In brief, 300-bp target fragments of the candidate gene were inserted between HindIII and XbalI sites in the vector, and the reverse sequence was inserted between XhoI and KpnI sites. The 300-bp target fragments contain the non-synonymous coding SNP site. The watermelon explants were transformed according to a modified method [17, 25, 26]. In particular, surface-sterilized watermelon seeds of watermelon cultivars Kexi (the wild type) were sown on basic Murashige-Skoog solid medium supplemented with 3% Suc and maintained for 2Â days in darkness and for 2 to 3Â days under a 16-h/8-h light/dark photoperiod. The cotyledons were excised, and hypocotyls and the apical portion were discarded. Agrobacterium tumefaciens strain GV3101 harboring the binary vector was used for the transformation. The cotyledon explants were co-cultivated in darkness for 4Â days. After transfer and maintenance on selective induction medium for 4Â weeks, the explants were transferred onto selective elongation medium containing 0.2Â mg/L KT, 60Â mg/L Kan, and 400Â mg/L Cef and incubated for 2Â weeks. Plantlets with well-developed roots were taken from the rooting medium and placed in plastic cups containing vermiculite. The traits of T0 plants were detected and the transgenic plants following self-pollination. The transgenic plants were self-pollinated for five generations to obtain stable phenotypes.
Results
Morphology of clalm mutant
To evaluate the genetic stability, clalm mutant were hand-pollinated using pollen from four inbred lines (‘Kexi’, ‘Hei2’, ‘Cai10’, and ‘K1’) to produce F1 and F2 generations. All plants of the F1 generation developed spotted leaves (Fig. 1 a–d), whereas the ratio of spotted- to non-spotted-leaf plants in the F2 population conformed to a Mendelian segregation ratio of 3:1 (χ2 ≤ 0.5, P < 0.05) (Table 2). The parents, F1 and F2 phenotypes of the constructed population were shown in fig. S3.
Necrotic tissue of clalm mutant spotted-leaf
In contrast to living tissue, necrotic tissue can be stained with trypan blue. Plants of clalm mutant were self-pollinated for seven generations to yield stable lesion mimic on leaves and fruits. Yellow spots appeared on leaves 10 days after sowing and then gradually spread over entire leaves and fruits during the life cycle of clalm mutant, the female parent (Fig. 2 a-d). Dead leaf-membrane cells will appear dark blue, whereas other parts should be almost colorless. We therefore stained leaves of the wild type and clalm mutant with trypan blue, which revealed that leaves of clalm mutant, but not those of the wild type, harbored small, spontaneous lesions (Fig. 2 e, f).
Comparative morphology of clalm mutantand wild-type ‘Zd’ watermelon. a-b Plants of lesion mimic mutant clalm (a) and ‘Zd’ (b) growing at the experimental breeding site. (c) Spotted leaf of lesion mimic mutant clalm. d Non-spotted leaf of ‘Zd’. e-f Leaves of lesion mimic mutant clalm (e) and ‘Zd’. f stained with trypan blue
Microscopic observations
No significant differences between speckled leaf areas and normal leaves were observed under a light microscope (Fig. 3a, b). Obvious differences were detected, however, using TEM. First, mesophyll cells of normal leaves were rich in chloroplasts, whereas those of clalm mutant spotted leaves contained very little chlorophyll (Fig. 3c, e). Second, mesophyll cell chloroplasts of normal leaves were rich in starch grains, but starch grains were absent from those of spotted leaves (Fig. 3d, f). Finally, mesophyll cell chloroplasts of normal leaves were rich in grana; in contrast, those of spotted leaves had no grana and only a few thylakoids (Fig. 3d, f).
Optical microscope imaging of wild-type ‘Zd’ (a) and lesion mimic mutant clalm (b) leaf cross sections. Electron microscope imaging of wild-type ‘Zd’ (c, d) lesion mimic mutant clalm (e, f) leaf cross sections. E, epidermis; PM, palisade mesophyll cells; SM, spongy mesophyll cells; S, stoma; V, vein; C, chloroplasts; A, amyloplasts; P, peroxisome; SG, starch grain
Whole-genome resequencing analysis
Filtering of the raw data yielded 43.89 G of clean data for further analysis. The Q30 ratio of the clean data was 85%, and the GC content was 41.67%, with an average genome resequencing depth of 22.75 × and a genome coverage of 98.81%. The data have been submitted to the NCBI database. Compared with the reference genome, 165,590 SNPs and 71,589 insertions–deletions (InDels) were identified between spotted- and non-spotted-leaf parents (Fig. S1). A Circos plot of the chromosomal distribution of candidate regions among samples is shown in Fig. 4. In addition, 12,784 SNPs and 19,212 InDels were found in the two bulk segregant populations. Association analyses between the two bulks based on SNP and InDel indexes narrowed the candidate region to approximately 4.04 Mb (Fig. 5 a) and 4.00 Mb (Fig. 5 b), respectively. Collectively, these two indexes indicated the candidate gene was located in an approximately 3.68 Mb region on chromosome 4 containing roughly 72 annotated genes. There are 6 genes with non-synonymous mutation SNP sites, ClCG04G000300, ClCG04G000420, ClCG04G001450, ClCG04G001740, ClCG04G001900, and ClCG04G001930, which were likely directly associated with leaf spotting.
Chromosomal distributions of SNP, ΔSNP, InDel, and ΔInDel index values. a Chromosomal distributions of SNP index values in the non-spotted-leaf plant bulk from the F2 population (top), SNP index values in the spotted-leaf bulk from the F2 population (middle), and the difference in SNP index values (ΔSNP) between the two bulks (bottom). b Distributions of InDel index values in the non-spotted-leaf plant bulk from the F2 population (top),InDel index values of the spotted-leaf bulk from the F2 population (middle), andΔIndel index values between the two bulks (bottom). The 99% threshold in each figure is indicated by a red dashed line, and Citrullus lanatus chromosome numbers are shown at the bottom
Identification of the LMM gene
We examined expression levels of six candidate genes (ClCG04G000300, ClCG04G000420, ClCG04G001450, ClCG04G001740, ClCG04G001900, and ClCG04G001930) in the two parental watermelon lines by qRT-PCR to assess whether their expression levels were correlated with the spotted leaf phenotype (Fig. 6). The expression level of ClCG04G001930 was significantly higher in the clalm mutant than in normal watermelon (P < 0.05), thus suggesting that this gene might be the candidate gene responsible for leaf spotting in clalm mutant. The SNP of ClCG04G001930 gene in the F2 population was examined by sequencing. Twenty-five SNPs were identified in ClCG04G001930, including a non-synonymous coding SNP, five SNPs in the ClCG04G001930 promoter and seventeen SNPs in the intron (Table 3). The SNP analysis resulted in one non-synonymous coding SNP locating in the ClCG04G001930 gene exon. Mutation of amino acid codon from GTA to ATA (Fig. S2). Amino acid changes from Ile to Val.
ClCG04G001930 was predicted to encode a PAD4, a protein that plays an important regulatory role in leaf senescence in many plant species. We therefore hypothesized that ClCG04G001930 is a PAD4 homolog in watermelon and accordingly named this gene ClPAD4. According to an amino acid sequence multiple alignment (Fig. 7 a), the ClPAD4 protein is 100% similar to PAD4 proteins of Cucumis sativus (accession number XP_011653897.1), Cucumis melo (XP_008442139.1), Cucurbita maxima (XP_022966244.1), and Momordica charantia (XP_022149409.1). The constructed phylogenetic tree clarified the molecular evolutionary relationship between ClPAD4 and its homologs (Fig. 7 b).
Bioinformatics characteristics of ClPAD4. a Results of Blast comparison of ClPAD4 proteins from various angiosperms. Query_10001 = Citrullus lanatus PAD4 (Gene ID:ClCG04G001930); Query_10002 = Benincasa hispida PAD4 (gene ID:LOC120072712); Query_10003 = Cucumis sativus PAD4 (gene ID:LOC101222098); Query_10004 = Cucumis melo PAD4 (gene ID:LOC103486086); Query_10005 = Cucurbita maxima PAD4 (gene ID:LOC111465963); Query_10006 = Momordica charantia PAD4 (gene ID:LOC111017842); Query_10007 = Quercus suber PAD4 (gene ID:LOC112032821); Query_10008 = Ziziphus jujuba PAD4 (gene ID:LOC107404877); Query_10009 = Tripterygium wilfordii PAD4 (gene ID:LOC119990818); Query_10010 = Arabidopsis thaliana PAD4 (gene ID:At3g52430). b Phylogenetic tree based on ClPAD4 in different organisms. The GenBank accession numbers are provided. The sequences are in Supplementary Material
Over-expression of ClPAD4 in watermelon
To characterize the function of ClPAD4, we generated transgenic watermelon plants for this gene. We obtained 35 transgenic plants and screened them for positive clones (Fig. 8a, b, c). The yellow spots gradually appeared on leaves of T0 ClPAD4 over-expression after Agrobacterium infection (Fig. 8c). The rooted shoots were transferred to plastic cups containing vermiculite (Fig. 8d). The chlorotic lesions on the transgenic leaves are marked with red boxes.
RNAi of ClPAD4 in watermelon
To characterize the function of ClPAD4, we generated transgenic knockdown watermelon plants for this gene. We obtained 56 transgenic plants and screened them for positive clones (Fig. 9a). The rooted shoots were transferred to plastic cups containing vermiculite (Fig. 9b). And the transgenic plants do not exhibit specific traits. The homozygous fruits were produced by the transgenic plants following self-pollination (Fig. 9c). It's very interesting that the cotyledons of T5 ClPAD4 RNAi watermelon seedlings susceptible to mosaic virus infection in a normal production environment which was not sterile environment (Fig. 9d). After transplantation, the leaves of the transgenic seedlings grew normally (Fig. 9e).
Discussion
The strategy used by watermelon breeders to improve the quality of watermelon cultivars is the introduction of novel traits, especially ones useful for seedling screening, into existing germplasm. The underlying gene in LMM of watermelon and other cucurbitaceous crops has not yet been well studied.
The timing of spot appearance can be used to circumscribe vegetative, reproductive, and complete growth periods [1, 27]. In our study, lesion mimic emerged on watermelon cotyledons 10Â days after sowing, which makes their presence an obvious marker for identifying impure and aberrant varieties at the seedling stage. The seedlings grow very slowly and flowered late, which was consistent with previous reports. In rice LMMs, leaf spots appear during the seedling stage and are relatively weakly correlated with the duration of the entire growth period [28]. We found that watermelon lesion mimicry is controlled by a single dominant gene. LMMs are relatively rare in cucurbitaceous species, with the spotting trait mainly manifested in fruits and seeds. In addition, a single dominant gene is responsible for fruit spotting in non-spotted varieties of Cucumis melo, Cucurbita pepo, and watermelon [6, 29,30,31,32]. In rice, a single dominant gene controls the LMM phenotype of NH1, spl12, spl13, spl15, and spl24 mutants, whose leaves display small, reddish-brown lesions or spots [33,34,35]. Many LMMs exhibit altered disease resistance and are thus considered ideal for studying signaling pathways in species such as Arabidopsis, maize, and rice [36,37,38,39,40].
Watermelon LMM genes are valuable experimental materials for studying molecular mechanisms and creating new germplasm resources. The advent of BSA-seq has accelerated the identification of candidate genes controlling important agronomic traits [12, 41]. This method has been used to map major QTLs for powdery mildew resistance to chromosome 12 in melon [42]. In addition, the candidate genes Csa2M435460.1 and Csa5M579560.1 conferring resistance to cucumber powdery mildew have been identified using BSA-seq [43]. Furthermore, a genome-wide analysis of SNPs resulted in the detection of a genomic region harboring the candidate dwarfism gene Cla010726 [12]. In contrast, few researchers have investigated the gene responsible for spontaneous lesion mimicry in watermelon and other cucurbitaceous crops. In this study, we used BSA-seq to map ClCG04G001930 of watermelon for the first time, which might be the candidate gene responsible for lesion mimic. We localized this gene to chromosome 4 in a region between 3,760,000 bp to 7,440,000 bp, corresponding to a physical distance of 3.68 Mb. ClCG04G001930 is a PAD4 homolog in watermelon. ClPAD4 and AtPAD4 have a protein sequence similarity of 62.72%, which suggests that they have similar functions (Fig. 8). The PAD4 orthologs are present in many plant species and essential for systemic resistance against biotic stress in angiosperms [44,45,46,47,48]. AtPAD4 is a member of a small family of sequence-related immunity regulator [49, 50]. Heterologous expression of the AtPAD4 gene in soybean roots inhibits the development of plant parasitic nematodes. GbPAD4 is up-regulated in cotton during pathogen infection [51] and LePAD4 expression is elevated in response to green peach aphid infestation in tomato [52]. The PAD4 protein in grape supports the response of the SA defense pathway to biotic stress [53]. ICS1, NPR1-3, PRs, EDS1, PAD4, and FMO1 signaling is strongly elicited during rust disease infection in rice [54].
The PAD4 is involved in the regulation of programmed cell death and acclimation to biotic and abiotic stresses [47, 48]. We detected an obvious difference between LMM and normal watermelon leaves under an electron microscope: the cells of watermelon LMM leaves contained very little chlorophyll, whereas those of normal leaves were rich in chloroplasts and grana (Fig. 3). Chloroplasts play an important role in regulating PAD4-modulated stress responses [55,56,57]. Our experimental result is thus very interesting and requires detailed study. The ClPAD4 over-expression transgenic watermelon plants exhibited spontaneous lesions (Fig. 8b). Overexpression of TaPAD4 produce necrotic spots to prevent the spread of powdery mildew, which validate the function of TaPAD4 in wheat powdery mildew resistance [58]. The ClPAD4 RNAi transgenic watermelon seedlings susceptible to mosaic virus infection (Fig. 9d). Silencing of OsPAD4 increases sensitivity to biotrophic pathogens in rice [46], and silencing of GmPAD4 reduces SA accumulation and enhances soybean susceptibility to virulent pathogens [59]. Previous studies have shown that the LMM phenotype of the mutant cpr5 is controlled by PAD4-dependent SA accumulation [37, 60]. Our experimental results, which are consistent with previous studies, indicate that the mutation of ClPAD4 in watermelon most likely led to the spontaneous lesions.
Genes used as screening markers for resistance breeding, biomass production, and productivity enhancement are important from an agricultural point of view. Watermelon breeders have always focused on the introduction of novel traits into existing germplasm as a means of improving cultivar quality. The clalm mutant exhibited obvious chlorotic lesions, which makes this application useful for seedling screening. Future elucidation of the mechanism of the PAD4 gene with respect to the spotted leaf trait should be useful in molecular marker-assisted selection breeding of watermelon.
Data availability
The datasets generated during the current study are available in the NCBI (BioSample accessions: SAMN45962846, SAMN45962847, SAMN45962848, SAMN45962849, SAMN45962850, SAMN45962851, SAMN45962852, SAMN45962853).
References
Qiao Y, Jiang W, Lee JH, Park B, Choi MS, Piao R, Woo MO, Roh JH, Han L, Paek NC, Seo HS, Koh HJ. SPL28 encodes a clathrin-associated adaptor protein complex 1, medium subunit μ1 (AP1M1) and is responsible for spooted leaf and early senescence in rice (Oryza sativa). New Phytol. 2010;185(1):258–74. https://doi.org/10.1111/j.1469-8137.2009.03047.x.
Yamanouchi U, Yano M, Lin HX, Ashikari M, Yamada K. A rice spotted leaf gene, spl7, encodes a heat shock transcription factor protein. Proc Natl Acad Sci U S A. 2002;99(11):7531–5. https://doi.org/10.1073/pnas.112209199.
Zeng LR, Qu S, Bordeos A, Yang C, Baraoidan M, Yan H, Xie Q, Nahm BH, Leung H, Wang GL. Spotted leaf 11, a negative regulator of plant cell death and defense, encodes a U-box/armadillo repeat protein endowed with E3 ubiquitin ligase activity. Plant Cell. 2004;16(10):2795–808. https://doi.org/10.1105/tpc.104.025171.
Fekih R, Tamiru M, Kanzaki H, Abe A, Yoshida K, Kanzaki E, Saitoh H, Takagi H, Natsume S, Undan JR, Undan J, Terauchi R. The rice (Oryza sativa L.) Lesion mimic esembling, which encodes an AAA-type ATPase, is implicated in defense response. Mol Genet Genemics. 2015;290(2):611–22. https://doi.org/10.1007/s00438-014-0944-z.
Ruan B, Hua Z, Zhao J, Zhang B, Ren D, Liu C, Yang S, Zhang A, Jiang H, Yu H, Hu J, Zhu L, Chen G, Shen L, Dong G, Zhang G, Zeng D, Guo L, Qian Q, Gao Z. OsACL-A2 negatively regulates cell death and disease resistance in rice. Plant Biotechnol J. 2019;17(7):1344–56. https://doi.org/10.1111/pbi.13058.
Lv R, Li Z, Li M, Dogra V, Lv S, Liu R, Lee KP, Kim C. Uncoupled expression of nuclear and plastid photosynthesis-associated genes contributes to cell death in a lesion mimic mutant. Plant Cell. 2019;31(1):210–30. https://doi.org/10.1105/tpc.18.00813.
Li R, Chen S, Liu G, Han R, Jiang J. Characterization and Identification of a woody lesion mimic mutant lmd, showing defence response and resistance to Alternaria alternate in birch. Sci Rep. 2017;7(1):11308. https://doi.org/10.1038/s41598-017-11748-2.
Liu S, Fu C, Gou J, Sun L, Huhman D, Zhang Y, Wang ZY. Simultaneous downregulation of MTHFR and COMT in switchgrass affects plant performance and induces lesion-mimic cell death. Front Plant Sci. 2017;8:982. https://doi.org/10.3389/fpls.2017.00982.
Kim H, Han D, Kang J, Choi Y, Levi A, Lee GP, Park Y. Sequence-characterized amplified polymorphism markers for selecting rind stripe pattern in watermelon (Citrullus lanatus L.). Hortic Environ Biote. 2015;56(3):341–9. https://doi.org/10.1007/s13580-015-0017-1.
Dou J, Lu X, Ali A, Zhao S, Zhang L, He N, Liu W. Genetic mapping reveals a marker for yellow skin in watermelon (Citrullus lanatus L.). PLoS One. 2018;13(9):e0200617. https://doi.org/10.1371/journal.pone.0200617.
Dou J, Zhao S, Lu X, He N, Zhang L, Ali A, Kuang H, Liu W. Genetic mapping reveals a candidate gene (ClFS1) for fruit shape in watermelon (Citrullus lanatus L.). Theor Appl Genet. 2018;131(4):947–58. https://doi.org/10.1007/s00122-018-3050-5.
Dong W, Wu D, Li G, Wu D, Wang Z. Next-generation sequencing from bulked segregant analysis identifies a dwarfism gene in watermelon. Sci Rep. 2018;8(1):2908. https://doi.org/10.1038/s41598-018-21293-1.
Gebremeskel H, Dou J, Li B, Zhao S, Muhammad U, Lu X, He N, Liu W. Molecular mapping and candidate gene analysis for GA3 responsive short internode in watermelon (Citrullus lanatus). Int J Mol Sci. 2019;21(1):290. https://doi.org/10.3390/ijms21010290.
Wei C, Zhu C, Yang L, Zhao W, Ma R, Li H, Zhang Y, Ma J, Yang J, Zhang X. A point mutation resulting in a 13 bp deletion in the coding sequence of Cldf leads to a GA-deficient dwarf phenotype in watermelon. Hortic Res. 2019;6:132. https://doi.org/10.1038/s41438-019-0213-8.
Zhu H, Zhang M, Sun S, Yang S, Li J, Li H, Yang H, Zhang K, Hu J, Liu D, Yang L. A single nucleotide deletion in an ABC transporter geneleads to a dwarf phenotype in watermelon. Front Plant Sci. 2019;10:1399. https://doi.org/10.3389/fpls.2019.01399.
Sandlin K, Prothro J, Heesacker A, Khalilian N, Okashah R, Xiang W, Bachlava E, Caldwell DG, Taylor CA, Seymour DK, et al. Comparative mapping in watermelon [Citrullus lanatus (Thunb.) Matsum. et Nakai]. Theor Appl Genet. 2012;125:1603–18. https://doi.org/10.1007/s00122-012-1938-z.
Ren Y, Guo S, Zhang J, He H, Sun H, Tian S, Gong G, Zhang H, Levi A, Tadmor Y, Xu Y. A tonoplast sugar transporter underlies a sugar accumulation QTL in watermelon. Plant Physiol. 2018;176(1):836–50. https://doi.org/10.1104/pp.17.01290.
Ren Y, McGregor C, Zhang Y, Gong G, Zhang H, Guo S, Sun H, Cai W, Zhang J, Xu Y. An integrated genetic map based on four mapping populations and quantitative trait loci associated with economically important traits in watermelon (Citrullus lanatus). BMC Plant Biol. 2014;14:33. https://doi.org/10.1186/1471-2229-14-33.
Fekih R, Takagi H, Tamiru M, Abe A, Natsume S, Yaegashi H, Sharma S, Sharma S, Kanzaki H, Matsumura H, Saitoh H, Mitsuoka C, Utsushi H, Uemura A, Kanzaki E, Kosugi S, Yoshida K, Cano L, Kamoun S, Terauchi R. MutMap+: genetic mapping and mutant identification without crossing in rice. PLoS One. 2013;8(7): e68529. https://doi.org/10.1371/journal.pone.0068529.
Hill JT, Demarest BL, Bisgrove BW, Gorsi B, Su YC, Yost HJ. MMAPPR: mutation mapping analysis pipeline for pooled RNA-seq. Genome Res. 2013;23(4):687–97. https://doi.org/10.1101/gr.146936.112.
Kumar S, Tamura K, Jakobsen IB, Nei M. MEGA2: molecular evolutionary genetics analysis software. Bioinformatics. 2001;17(12):1244–5. https://doi.org/10.1093/bioinformatics/17.12.1244.
Felsenstein,. Phylogenies from restriction sites: a maximum-likelihood approach. Evolution. 1992;46(1):159–73. https://doi.org/10.1111/j.1558-5646.1992.tb01991.x.
Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, Gumienny R, Heer FT, de Beer TAP, Rempfer C, Bordoli L, Lepore R, Schwede T. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 2018;46(W1):W296–303. https://doi.org/10.1093/nar/gky427.
Dong W, Wu DW, Wang CH, Liu Y, Wu DF. Characterization of the molecular mechanism underlying the dwarfism of dsh mutant watermelon plants. Plant SCI. 2021;313(2021):111074. https://doi.org/10.1016/j.plantsci.2021.111074.
Yu TA, Chiang CH, Wu HW, Li CM, Yang CF, Chen JH, Chen YW, Yeh SD. Generation of transgenic watermelon resistant to Zucchini yellow mosaicvirus and papaya ringspot virus type W. Plant Cell Rep. 2011;30(3):359–71. https://doi.org/10.1007/s00299-010-0951-4.
Liu L, Gu Q, Ijaz R, Zhang J, Ye Z. Generation of transgenic watermelon resistance to cucumber mosaicvirus facilitated by an effective Agrobacterium-mediatedtransformation method. Sci Hortic. 2016;205:32–8. https://doi.org/10.1016/j.scienta.2016.04.013.
Ma J, Wang Y, Ma X, Meng L, Jing R, Wang F, Wang S, Cheng Z, Zhang X, Jiang L, Wang J, Wang J, Zhao Z, Guo X, Lin Q, Wu F, Zhu S, Wu C, Ren Y, Lei C, Zhai H, Wan J. Affiliations expand disruption of gene SPL35, encoding a novel CUE domain-containing protein, leads to cell death and enhanced disease response in rice. Plant Biotechnol J. 2019;17(8):1679–93. https://doi.org/10.1111/pbi.13093.
Zhang XB, Feng BH, Wang HM, Xu X, Shi YF, He Y, Chen Z, Sathe AP, Shi L, Wu JL. A substitution mutation in Os PELOTA confers bacterial blight resistance by activating the salicylic acid pathway. J Integr Plant Biol. 2018;60(2):160–72. https://doi.org/10.1111/jipb.12613.
Paris HS. Gene for broad, contiguous dark stripes in cocozelle squash (Cucurbita pepo). Euphytica. 2000;115(3):191–6. https://doi.org/10.1023/A:1004074618717.
Paris HS. Genetic control of irregular striping, a new phenotype in Cucurbita pepo. Euphytica. 2002;129(1):119–26. https://doi.org/10.1023/A:1021564610029.
Ntui V, Uyoh E. Inheritance of stripe pattern on fruits and seed colour in Begusi melon, Colocynthis citrullus L. Glob J Agric Sci. 2005;4(1):29–32. https://doi.org/10.4314/gjass.v4i1.2251.
Pitrat M. 2002 gene list for melon. Cucurb Genet Coop Rep. 2002;25:1–18.
Mizobuchi R, Hirabayashi H, Kaji R, Nishizawa Y, Yoshimura A, Satoh H, Ogawa T, Okamoto M. Isolation and characterization of rice lesion mimic mutants with enhanced resistance to rice blast and bacterial blight. Plant Sci. 2002;163:345–53. https://doi.org/10.1016/S0168-9452(02)00134-6.
Wu C, Bordeos A, Madamba MR, Baraoidan M, Ramos M, Wang GL, Leach JE, Leung H. Rice lesion mimic mutants with enhanced resistance to diseases. Mol Genet Genomics. 2008;279(6):605–19. https://doi.org/10.1007/s00438-008-0337-2.
Yuan Y, Zhong S, Li Q, Zhu Z, Lou Y, Wang L, Wang J, Wang M, Li Q, Yang D, He Z. Functional analysis of rice NPR1-like genes reveals that Os NPR1/NHI is the rice orthologue conferring disease resistance with enhanced herbivore susceptibility. Plant Biotechno J. 2007;5(2):313–24. https://doi.org/10.1111/j.1467-7652.2007.00243.x.
Johal GS, Hulbert S, Briggs SP. Disease lesion mimics of maize: a model for cell death in plants. BioEssays. 1995;17(8):685–92. https://doi.org/10.1002/bies.950170805.
Lorrain S, Lin B, Auriac MC, Kroj T, Saindrenan P, Nicole M, Balagué C, Roby D. Vascular associated death1, a novel GRAM domain-containing protein, is a regulator of cell death and defense responses in vascular tissues. Plant Cell. 2004;16(8):2217–32. https://doi.org/10.1105/tpc.104.022038.
Shirsekar GS, Vega-Sanchez ME, Bordeos A, Baraoidan M, Swisshelm A, Fan J, Park CH, Leung H, Wang GL. Identification and characterization of suppressor mutants of spl11-mediated cell death in rice. Mol Plant Microbe Interact. 2014;27(6):528–36. https://doi.org/10.1094/MPMI-08-13-0259-R.
Wang C, Zhang X, Fan Y, Gao Y, Zhu Q, Zheng C, Qin T, Li Y, Che J, Zhang M, Yang B, Liu Y, Zhao K. XA23 is an executor R protein and confers broad-spectrum disease resistance in rice. Mol Plant. 2015;8:290–302. https://doi.org/10.1016/j.molp.2014.10.010.
Wang S, Lei C, Wang J, Ma J, Tang S, Wang C, Zhao K, Tian P, Zhang H, Qi C, Cheng Z, Zhang X, Guo X, Liu L, Wu C, Wan J. SPL33, encoding an eEF1A-like protein, negatively regulates cell death and defense responses in rice. J Exp Bot. 2017;68(5):899–913. https://doi.org/10.1093/jxb/erx001.
Lv J, Fu Q, Lai Y, Zhou M, Wang H. Inheritance and gene mapping of spotted to non-spotted trait gene CmSp-1 in melon (Cucumis melo L. var. chinensis Pangalo). Mol Breeding. 2018;38:105. https://doi.org/10.1007/s11032-018-0860-8.
Li B, Zhao Y, Zhu Q, Zhang Z, Fan C, Luan F, Gao P. Melon powdery mildew resistance gene location and candidate gene analysis. China vegetables. 2017;2017(6):17–24.
Xu Q, Yang S, Yu T, Xu X, Yan Y, Qi X, Chen X. Wholegenome resequencing of a cucumber chromosome segment substitution line and its recurrent parent to identify candidate genes governing powdery mildew resistance. PLoS One. 2016;11(10): e0164469. https://doi.org/10.1371/journal.pone.0164469.
Rusterucci C, Aviv DH, Holt BF 3rd, Dangl JL, Parker JE. The disease resistance signalling components EDS1 and PAD4 are essential regulators of the cell death pathway controlled by LSD1 in Arabidopsis. Plant Cell. 2001;13(10):2211–24. https://doi.org/10.1105/tpc.010085.
Gao F, Dai R, Pike SM, Qiu W, Gassmann W. Functions of EDS1-like and PAD4 genes in grapevine defenses against powdery mildew. Plant Mol Biol. 2014;86(4–5):381–93. https://doi.org/10.1007/s11103-014-0235-4.
Ke Y, Liu H, Li X, Xiao J, Wang S. Rice OsPAD4 functions differently from Arabidopsis AtPAD4 in host-pathogen interactions. Plant J. 2014;78(4):619–31. https://doi.org/10.1111/tpj.12500.
Yan Z, Xingfen W, Wei R, Jun Y, Zhiying M. Island cotton enhanced disease susceptibility 1 gene encoding a lipase-like protein plays a crucial role in response to Verticillium dahliae by regulating the SA level and H2O2 accumulation. Front Plant Sci. 2016;7:1830. https://doi.org/10.3389/fpls.2016.01830.
Chen G, Wei B, Li G, Gong C, Fan R, Zhang X. TaEDS1 genes positively regulate resistance to powdery mildew in wheat. Plant Mol Biol. 2018;96(6):607–25. https://doi.org/10.1007/s11103-018-0718-9.
Feys BJ, Wiermer M, Bhat RA, Moisan LJ, Medina-Escobar N, Neu C, Cabral A, Parker JE. Arabidopsis senescence-associated gene101 stabilizes and signals within an Enhanced disease susceptibility1 complex in plant innate immunity. Plant Cell. 2005;17(9):2601–13. https://doi.org/10.1105/tpc.105.033910.
Lapin D, Kovacova V, Sun X, Dongus JA, Bhandari D, Born P, Bautor J, Guarneri N, Rzemieniewski J, Stuttmann J, Beyer A, Parker JE. A coevolved EDS1-SAG101-NRG1 module mediates cell death signaling by TIR-domain immune receptors. Plant Cell. 2019;31(10):2430–55. https://doi.org/10.1105/tpc.19.00118.
Zhang B, Yang Y, Chen T, Yu W, Liu T, Li H, Fan X, Ren Y, Shen D, Liu L, Dou D, Chang Y. Island cotton Gbve1 gene encoding a receptor-like protein confers resistance to both defoliating and non-defoliating isolates of Verticillium dahliae. PLoS One. 2012;7(12):e51091. https://doi.org/10.1371/journal.pone.0051091.
Singh V, Shah J. Tomato responds to green peach aphid infestation with the activation of trehalose metabolism and starch accumulation. Plant Signal Behav. 2012;7(6):605–7. https://doi.org/10.4161/psb.20066.
Tandon G, Jaiswal S, Iquebal MA, Kumar S, Kaur S, Rai A, Kumar D. Evidence of salicylic acid pathway with EDS1 and PAD4 proteins by molecular dynamics simulation for grape improvement. J Biomol Struct Dyn. 2015;33(10):2180–91. https://doi.org/10.1080/07391102.2014.996187.
Sahu R, Prabhakaran N, Kundu P, Kumar A. Differential response of phytohormone signalling network determines nonhost resistance in rice during wheat stem rust (Puccinia graminis f. sp. tritici) colonization. Plant Pathol. 2021;70(6):1409–20. https://doi.org/10.1111/ppa.13376.
Mühlenbock P, Szechynska-Hebda M, Płaszczyca M, Baudo M, Mateo A, Mullineaux PM, Parker JE, Karpinska B, Karpinski S. Chloroplast signaling and LESION SIMULATING DISEASE1 regulate crosstalk between light acclimation and immunity in Arabidopsis. Plant Cell. 2008;20(9):2339–56. https://doi.org/10.1105/tpc.108.059618.
Karpi ´nski S, Szechy ´nska-Hebda M, Wituszy ´nska W, Burdiak P. Light acclimation, retrograde signalling, cell death and immune defences in plants. Plant Cell Environ. 2013;36(4):736–44. https://doi.org/10.1111/pce.12018.
Bernacki M, Czarnocka W, Szechy ´nska-Hebda M, Mittler R, Karpi ´nski S. Biotechnological potential of LSD1, EDS1, and PAD4 in the improvement of crops and industrial plants. Plants. 2019;8(8):290. https://doi.org/10.3390/plants8080290.
Song N, Lin J, Liu X, Liu Z, Liu D, Chu W, Li J, Chen Y, Chang S, Yang Q, Liu X, Guo W, Xin M, Yao Y, Peng H, Ni Z, Xie C, Sun Q, Hu Z. Histone acetyltransferase TaHAG1 interacts with TaPLATZ5 to activate TaPAD4 expression and positively contributes to powdery mildew resistance in wheat. New Phytol. 2022;236(2):590–607. https://doi.org/10.1111/nph.18372. Epub 2022 Jul 28 PMID: 35832009.
Wang J, Shine MB, Gao QM, Navarre D, Jiang W, Liu C, Chen Q, Hu G, Kachroo A. Enhanced disease susceptibility1 mediates pathogen resistance and virulence function of a bacterial effector in soybean. Plant Physiol. 2014;165(3):1269–84. https://doi.org/10.1104/pp.114.242495.
Jirage D, Zhou N, Cooper B, Clarke JD, Dong X, Glazebrook J. Constitutive salicylic acid-dependent signaling in cpr1 and cpr6 mutants requires PAD4. Plant J. 2001;26(4):395–407. https://doi.org/10.1046/j.1365-313x.2001.2641040.x.
Acknowledgements
We gratefully acknowledge the support provided by the National Natural Science Foundation of China (31801882) and Henan Province Key Science and Technology Projects (242102111085).
Funding
This work was supported by funding from the National Natural Science Foundation of China (31801882) and Henan Province Key Science and Technology Projects (242102111085).
Author information
Authors and Affiliations
Contributions
W.D. wrote the main manuscript text. W.D. and D.W. prepared figures 1-3 and figure S3. J.S., C.Y., Y.Q., Q.L. and S.H. prepared figure 4-9 and figure S1-S2. J.S. and C.Y. contributed equally. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
12870_2025_6071_MOESM1_ESM.pdf
Supplementary Material 1. Fig. S1 Venn diagrams of small InDels (a) and SNPs (b) in spotted-leaf and non-spotted-leaf watermelon. R01, clalm mutant; R02, ‘Zd’; R03, spotted-leaf plants in the F2 population; R04, non-spotted-leaf plants in the F2 population.
12870_2025_6071_MOESM2_ESM.pdf
Supplementary Material 2. Fig. S2 Predicted protein tertiary structure of ClPAD4. Mutation of amino acid codon from Ile to Val. Marked with red boxes and tables.
12870_2025_6071_MOESM3_ESM.pdf
Supplementary Material 3. Fig. S3 The photos of F1 and F2 phenotypes of the constructed population. (a) The experimental materials before collection. (b) All plants of the F1 generation developed spotted leaves. (c) The non-spotted-leaf plants in the F2. (d) The spotted-leaf plants in the F2.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Shi, J., Yang, C., Qin, Y. et al. Phytoalexin deficient 4 is associated with the lesion mimic trait in watermelon clalm mutant (Citrullus lanatus). BMC Plant Biol 25, 92 (2025). https://doi.org/10.1186/s12870-025-06071-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12870-025-06071-2








