*4.4. Gene Isolation and Phylogenetic Analysis of GmWRKY12*

Venn2.0 (http://bioinfogp.cnb.csic.es/tools/venny/index.html) was used to screen GmWRKYs that respond to both drought and salt treatment, then qRT-PCR was used to find genes highly expressed under stresses. Full-length *GmWRKY12* was amplified by PCR with specific primers from soybean cDNA (*Williams 82*); primers of *GmWRKY12* are available in Table S4. PCR products were cloned into pLB vector (TIANGEN, Beijing, China) and sequenced for further study. The amino acid sequence of WRKY12 in different species were searched for in the NCBI database on account of the amino acid similarity between GmWRKY12 and WRKY12 in different species is more than 50%. DNAMAN was applied for multiple sequence alignment on the basis of the amino acid similarity between GmWRKY12 and WRKY12 in different species is more than 60%. Phylogenetic trees were constructed using MEGA 6.0 with the neighbor-joining method [66] and 1000 bootstrap replications. Information of WRKY12 in different species is listed in Table S3.
