*2.1. Identification of GmWRKYs Up-Regulated under Drought/Salt Treatment*

The GmWRKYs are distributed in different tissues or located upstream of soybean genes to bind the W-box consensus (TTGACY) in the promoters of target genes, initiating functions such as plant development, pathogen defense, insect resistance, response to biotic and abiotic stress and participating in signal transduction mediated by plant hormones [59,60]. In order to identify the function of genes or to explore whether GmWRKY mRNA expression goes up under biotic and abiotic stress, we conducted RNA-Seq (Tables S5 and S6). RNA-Seq data were used to screen GmWRKYs that are responsive to drought and salt. There were 105 GmWRKYs upregulated after drought treatment and fifty-three GmWRKYs were selected based on the rule that log<sup>2</sup> (GH\_treat/CK1\_treat) >1 (Table 1). Nine GmWRKYs were selected from salt treatment RNA-Seq data based on the rule that log<sup>2</sup> (NaCl\_treat/CK2\_treat) >1 (Table 2).


**Table 1.** Annotation of *Glycine max* WRKY transcription factors responding to drought stress (up-regulation).


**Table 1.** *Cont.*

<sup>a</sup>—The annotated GmWRKYs according to NCBI (https://www.ncbi.nlm.nih.gov/pubmed) and. PlantTFDB (http://planttfdb.cbi.pku.edu.cn/); <sup>b</sup>—The names of GmWRKYs are given according to SoyDB (http://soykb.org/); <sup>c</sup>—The grouping is according to [30,61].

**Table 2.** Annotation of *Glycine max* WRKY transcription factors responding to salt stress (up-regulation).


<sup>a</sup>—The annotated GmWRKYs according to NCBI (https://www.ncbi.nlm.nih.gov/pubmed) and PlantTFDB (http: //planttfdb.cbi.pku.edu.cn/). <sup>b</sup>—The names of GmWRKYs are given according to SoyDB (http://soykb.org/). <sup>c</sup>—The grouping is according to [30,61].
