16The WRKY proteins constitute a large family of transcription factors that have been known to 17 play a wide range of regulatory roles in multiple biological processes. Over the past few years, 18 many reports have focused on analysis of evolution and biological function of WRKY genes at the 19 whole genome level in different plant species. However, little information is known about WRKY 20 genes in melon (Cucumis melo L.). In the present study, a total of 56 putative WRKY genes were 21 identified in melon, which were randomly distributed on their respective chromosomes. A 22 multiple sequence alignment and phylogenetic analysis using melon, cucumber and watermelon 2 23 predicted WRKY domains indicated that melon WRKY proteins could be classified into three 24 main groups (I-III). Our analysis indicated that no recent duplication events of WRKY genes were 25 detected in melon, and strong purifying selection was observed among the 85 orthologous pairs of 26 Cucurbitaceae species. Expression profiles of CmWRKY derived from RNA-seq data and 27 quantitative RT-PCR (qRT-PCR) analyses showed distinct expression patterns in various tissues, 28 and the expression of 16 CmWRKY were altered following powdery mildew infection in melon.
29Besides, we also found that a total of 24 WRKY genes were co-expressed with 11 VQ family genes 30 in melon. Our comparative genomic analysis provides a foundation for future functional dissection 31 and understanding the evolution of WRKY genes in cucurbitaceae species, and will promote 32 powdery mildew resistance study in melon. 33 Keywords: cucurbitaceae, WRKY, phylogenetic relationship, evolution, powdery mildew 34 disease 35 36 39 number of functions in a wide range of physiological and biochemical processes [1-5]. The 40 WRKY family proteins are defined by a highly conserved domain with approximately 60 amino 41 acid residues in length, which contains one or two highly conserved short peptide WRKYGQK as 42 well as a conserved C2H2-or C2HC-type zinc finger motif [6-8]. The conserved short peptide also 43 consists of various forms, such as WRKYGKK, WRKYDQK, and WRKYDHK. Based on the 44 number of WRKY domains and the type of zinc finger motif, the WRKY family proteins can be 3 45 classified into three groups [7,8]. Group I contains two WRKY motifs and a C2H2 type zinc 46 finger motif (C-X4-5-C-X22-23-H-X1-H). Group II contains one WRKY motif and a C2H2 type 47 59 correlated with different seed size [19]. Besides, VvWRKY26, a TTG2-like homolog protein, 60 controls vacuolar acidification, transport and flavonoid biosynthesis in grapevine, and plays 61 important regulatory roles in fleshy fruit development [17].
62Recent studies revealed that a large part of WRKY proteins were also involved in response to 63 biotic (bacterial, fungal and viral pathogens) defense and abiotic stress (salinity, drought, heat, 64 cold and osmotic stress) [21]. In many plants, expression of WRKY family genes was significantly 100