S1). By comparison of the amino acid sequences in the WRKY domain regions from Gossypium and Arabidopsis, 120 cotton WRKY candidate genes were classified

into three groups (groups I, II, and III), and group II genes were further classified into five subgroups (groups IIa–e; Fig. 2), based on the classification rules employed for the WRKY family genes in Arabidopsis [4]. Among the three groups, there were 20 members in group I, 88 in group II, and 12 in group III. Furthermore, in group II, subgroups IIa–e contained 7, 16, 37, 15, and 13 members, respectively. The types and chromosome distribution of these members are described in Table 1. It is noteworthy that WRKY108 in group I contained three WRKY domains (WRKY108N1, WRKY108N2, and WRKY108C). However, the three WRKY Belinostat mw domains were not clustered in the N-terminal WRKY domain (NTWD) and the C-terminal WRKY domain (CTWD). The phylogenetic results showed that WRKY108N1, WRKY108N2, and WRKY108C were clustered into group IIc, group III, and group IId, respectively ( Fig. 2). According to D5 genomic sequence information, there was at least one intron insert in the WRKY candidate genes, with WRKY108 and WRKY109 having the most complex structures. The intron splices

in the conserved WRKY domain could be classified into GW 572016 two major types, the R type and the V type. V-type introns were observed only in groups IIa and IIb ( Fig. S2). In addition to the WRKY domain, the WRKY family members were also predicted by MEME to contain other conserved motifs. However, six WRKY proteins,

encoded by WRKY14, WRKY21, WRKY35, WRKY46, WRKY77, and WRKY90, contained only a WRKY domain ( Fig. S3). WRKYGQK residues are considered to be important regions of the WRKY transcription factor family. However, we found some genes with diverse amino acid residues PLEK2 in this region. Among the seven amino acid residues (WRKYGQK), mutations at the W and K sites were not observed; most variations involved Q to T, H, or K substitutions. For WRKY109 in group I, there were large variations in this seven residue regions in both NTWD and CTWD, with variations in three and four amino acid residues, respectively. In total, ten members showed divergence in the WRKY domain, of which seven belonged to group IIc (Table S3). In addition to the variations in amino acid residues in the WRKY DNA binding domain, some mutations were discovered in the zinc finger motif regions. Four members, including WRKY35 and WRKY114 in group I and WRKY108 and WRKY109 in group III, exhibited variations in amino acid residues in this motif (Table S4). By designing gene-specific primers (Table S5), we performed PCR cloning of WRKY genes and amplified the transcripts in given tissues of G. hirsutum acc. TM-1.