The weeping forsythia genome was obtained from National Center for Biotechnology Information (accession no. JAHHPY000000000; Li et al., 2022). The keyword “WRKY” was used to search WRKY genes in the annotation file, and blasted the candidate genes in NCBI (https://blast.ncbi.nlm.nih.gov/; Altschul et al. 1990 (link)) to identify the WRKY domain. The genes with conserved WRKY domain were considered the true WRKY genes. ExPASy online tool (https://web.expasy.org/cgi-bin/protparam/protparam; Artimo et al., 2012 (link)) was used to predict and analyze the physicochemical properties of WRKY protein in weeping forsythia, including molecular weight, isoelectric point, amino acid number, fat index, instability index, and hydrophobicity. The subcellular localization of weeping forsythia WRKY genes was predicted by Plant-mPLoc online software (http://www.csbio.sjtu.edu.cn/bioinf/plant-multi/, Chou and Shen, 2010 (link)), and the online software SOPMA (https://npsa-prabi.ibcp.fr/cgi-bin/secpred_sopma.pl; Geourjon and Deleage, 1995 (link)) was used to predict the secondary structure of WRKY protein.
According to the identified WRKY gene ID of weeping forsythia and the genome sequence of weeping forsythia, the chromosome position information of WRKY gene family members was obtained. The position of WRKY gene of weeping forsythia on chromosome was visualized by TBtools software (Chen et al. 2020a (link)). The genome and protein sequence data of Arabidopsis thaliana is from the Arabidopsis information resource (https://www.arabidopsis.org/; Eulgem et al., 2000 (link)).