Background {#Sec1}
==========

Sesame (*Sesamum indicum* L.) is an important, and probably the most ancient, oil crop and is grown widely in tropical and subtropical regions of the world \[[@CR1]\]. Recently, the demand for sesame has increased, but sesame yields have been poor compared with those of other oil crops (e.g., rapeseed: 1939.3 kg/ha; soybean: 2498.5 kg/ha; and peanut: 1657.6 kg/ha). Average sesame yields were alarmingly low between 2010 and 2014, with only 576.1 kg/ha being produced in 73 countries (<http://faostat.fao.org>). This low yield may be attributable to a variety of factors, although abiotic stresses are certainly one of the most significant.

For sesame, the most important abiotic stresses that limit plant growth, development, and yield are drought and waterlogging. In central China, which is the major sesame production area, sesame is generally planted during the rainy season, when waterlogging is the most significant problem and can decrease sesame yields by more than 80% \[[@CR2]\]. Sesame is also grown extensively in the tropical regions of Africa and South America. Here, drought presents a major challenge and can limit the yield from sesame by affecting the number of capsules produced by each plant \[[@CR3]\]. Therefore, there is an urgent requirement to understand the molecular mechanisms that underlie the ability of sesame plants to tolerate both drought and waterlogging stresses.

Abiotic stress responses and gene regulation have been studied in a number of plant species, including *Arabidopsis*, rice, maize, and tomato. Several families of genes are particularly associated with significant improvements in abiotic stress tolerance, including the *WRKY*, *NAC*, and *ERF* gene families \[[@CR4]--[@CR6]\]. Numerous studies have demonstrated that *WRKY* genes are expressed strongly and rapidly in response to particular abiotic stresses, including wounding, waterlogging, drought, and salt stress \[[@CR7]--[@CR9]\]. In *Arabidopsis*, *AtWRKY30* is induced by methyl viologen, hydrogen peroxide, arsenic, drought, sodium chloride, and mannitol \[[@CR10]\]. Nuruzzaman et al. \[[@CR11]\] identified five *OsWRKY* genes expressed at higher levels in drought-tolerant rice compared with those in drought-susceptible rice under experimental water-deficit conditions. Overexpression of *OsWRKY47* increased both the yield and drought tolerance compared with wild-type plants \[[@CR12]\]. Meng et al. \[[@CR13]\] discovered that 10 selected *WRKY* genes showed differential expression patterns under waterlogging and drought stress in an apple rootstock. These observations suggest that studying *WRKY* gene families may provide valuable insights into the mechanism underlying abiotic stress tolerance in plants. Furthermore, although drought and waterlogging may primarily affect plants grown in different parts of the world, very little is known about the identity and functions of *WRKY* genes in sesame.

Over the last decade, WRKY transcription factors have become one of the most extensively studied gene families involved in regulating plant abiotic stress tolerance \[[@CR14]\]. WRKY proteins have one or two unique DNA-binding domains that are approximately 60 amino acids (aa) in length and contain the WRKYGQK sequence followed by a C~2~H~2~ zinc-finger-like motif \[[@CR15]\]. The DNA-binding region is designated a WRKY domain because the WRKYGQK aa sequence is completely conserved. The WRKY proteins are classified into three major groups (I--III) based on the number of WRKY domains and the pattern of zinc-finger-like motifs. Group II is further divided into five distinct subgroups (IIa--IIe) \[[@CR15]\]. The first identified *WRKY* gene, *SPF1*, was cloned from sweet potato (*Ipomoea batatas*) 20 years ago \[[@CR16]\]. Since then, a large number of *WRKY* genes have been identified, including 74 from *Arabidopsis thaliana* \[[@CR17]\], 103 from *Oryza sativa* \[[@CR18]\], 45 from *Hordeum vulgare* \[[@CR19]\], 119 from *Zea mays* \[[@CR20]\], and many more from other plant species \[[@CR21]--[@CR23]\].

Sequencing of the entire sesame genome, and annotation of its 24,148 putative genes, provides an opportunity to identify all the sesame *WRKY* genes \[[@CR24]--[@CR26]\]. In this study, 71 *WRKY* genes were identified from the sesame genome and analyses of their structure, phylogeny, chromosomal distribution and duplication, conserved motifs, and stimulation in response to waterlogging and drought were performed. The results provide insights into the evolution of the sesame *WRKYs* and their functions in abiotic stress responses. The identification and characterization of these *WRKY* genes may provide opportunities to improve the stress tolerance of sesame.

Results {#Sec2}
=======

Identification of *WRKY* family genes in sesame {#Sec3}
-----------------------------------------------

All *Arabidopsis* WRKY protein sequences were used as queries for the Basic Local Alignment Search Tool (BLAST) to identify sesame WRKY proteins. In total, 61 putative *WRKY* genes were identified and predicted protein sequences without a WRKY domain were excluded. A Hidden Markov Model (HMM) search was also performed against the sesame protein database using the WRKY-domain PF03106. An additional 10 protein sequences containing the complete WRKY domain were identified. In total, 71 WRKY proteins were identified in the sesame genome (Table [1](#Tab1){ref-type="table"}).Table 1Informations of *SiWRKY* genesGene symbolGene locusLinkage groupPeptide lengthpIMWGroup*SiWRKY1*SIN_1000785scaffold002331867.2620.85NG*SiWRKY2*SIN_1001523scaffold001643025.6333.94IIe*SiWRKY3*SIN_1001786LG123006.3033.82IIc*SiWRKY4*SIN_1001880LG047606.2884.34I*SiWRKY5*SIN_1002759scaffold001243555.9540.01III*SiWRKY6*SIN_1002960scaffold001201229.7414.44IIc*SiWRKY7*SIN_1003153LG147285.9678.03I*SiWRKY8*SIN_1003599scaffold00109367NoNoIIc*SiWRKY9*SIN_1003920LG055406.2760.10I*SiWRKY10*SIN_1003975LG131769.2620.25IIc*SiWRKY11*SIN_1004161LG151879.4721.25IIc*SiWRKY12*SIN_1004874LG152919.7231.60IId*SiWRKY13*SIN_1005422LG022979.2132.27IId*SiWRKY14*SIN_1005676LG113154.8634.41IIe*SiWRKY15*SIN_1005706LG112015.1322.75III*SiWRKY16*SIN_1006024LG085616.1460.65IIb*SiWRKY17*SIN_1006129LG074975.7754.61NG*SiWRKY18*SIN_1006550LG085138.2856.11I*SiWRKY19*SIN_1006749LG124625.3849.64IIe*SiWRKY20*SIN_1006978LG125736.4461.81IIb*SiWRKY21*SIN_1007987LG154796.0652.05IIb*SiWRKY22*SIN_1008040LG155118.0455.40I*SiWRKY23*SIN_1009399scaffold000576268.2869.68I*SiWRKY24*SIN_1009643LG063459.5538.85IId*SiWRKY25*SIN_1009858LG111516.2317.23NG*SiWRKY26*SIN_1010783LG012968.6732.83IIa*SiWRKY27*SIN_1010982LG113276.4336.72III*SiWRKY28*SIN_1011023LG115265.7958.18I*SiWRKY29*SIN_1011192LG1111418.49125.94IIc*SiWRKY30*SIN_1011284LG114529.1049.38I*SiWRKY31*SIN_1011416LG033469.7037.75IId*SiWRKY32*SIN_1012054LG043165.2735.60III*SiWRKY33*SIN_1012055LG043375.7037.18III*SiWRKY34*SIN_1012623LG062814.8331.76NG*SiWRKY35*SIN_1012631LG062935.2632.48IIc*SiWRKY36*SIN_1012891LG063369.7037.73IId*SiWRKY37*SIN_1014111LG013326.3136.70IIe*SiWRKY38*SIN_1014143LG013505.0539.25III*SiWRKY39*SIN_1014268LG127236.1378.07I*SiWRKY40*SIN_1014366LG125846.5263.04I*SiWRKY41*SIN_1014422LG123977.2143.06IIb*SiWRKY42*SIN_1015494LG061878.5521.30IIc*SiWRKY43*SIN_1015496LG063385.5137.95IIe*SiWRKY44*SIN_1016166LG033726.1741.48IIc*SiWRKY45*SIN_1016382LG036025.6264.69IIb*SiWRKY46*SIN_1016491LG043654.8139.49IIe*SiWRKY47*SIN_1016829LG165916.0862.98IIb*SiWRKY48*SIN_1017975LG013305.3037.22III*SiWRKY49*SIN_1017989LG021529.3017.70IIc*SiWRKY50*SIN_1018215LG025648.0661.82I*SiWRKY51*SIN_1018227LG023345.5935.28IIc*SiWRKY52*SIN_1018859LG043166.0034.30IIc*SiWRKY53*SIN_1019334LG143167.7435.63IIc*SiWRKY54*SIN_1019555LG08444NoNoIIe*SiWRKY55*SIN_1019627LG083126.1134.81IIc*SiWRKY56*SIN_1019661LG082854.8132.22IIe*SiWRKY57*SIN_1019937LG055046.1654.40IIb*SiWRKY58*SIN_1020605LG063469.6536.99IId*SiWRKY59*SIN_1020883LG061607.7618.49IIc*SiWRKY60*SIN_1021497LG011686.5218.90IIc*SiWRKY61*SIN_1021618LG012688.9229.50IIa*SiWRKY62*SIN_1021622LG012558.8128.19IIa*SiWRKY63*SIN_1021665LG013349.5436.05IId*SiWRKY64*SIN_1021953LG033089.3034.42IIa*SiWRKY65*SIN_1022426LG063416.8338.16IIc*SiWRKY66*SIN_1022971LG084938.3754.34IIb*SiWRKY67*SIN_1023226LG065428.0859.73I*SiWRKY68*SIN_1026059LG104936.5854.73IIb*SiWRKY69*SIN_1026464LG081628.8918.95IIc*SiWRKY70*SIN_1026809LG083658.8539.90IIb*SiWRKY71*SIN_1026948LG085646.1960.99IIb*pI* proteins' isoelectric point, *MW* molecular weight

The 71 sesame WRKY proteins ranged from 122 (*SiWRKY6*) to 1141 (*SiWRKY29*) aa in length, with an average length of approximately 390 aa. The molecular weights (MWs) ranged from 14.44 kDa (*SiWRKY6*) to 125.94 kDa (*SiWRKY29*). The isoelectric points (pIs) of the WKRY proteins ranged from 4.81 (*SiWRKY46* and *SiWRKY56*) to 9.74 (*SiWRKY6*), with 39 pIs \<7 and the remaining pIs \>7 (Table [1](#Tab1){ref-type="table"}). Similar observations were made in Chinese cabbage \[[@CR27]\], which has 56 *BcWRKYs* with pIs ranging from 4.69 to 10.45 and MWs ranging from 20.44 kDa to 119.84 kDa.

Chromosomal locations of and duplication events of the *SiWRKY* genes {#Sec4}
---------------------------------------------------------------------

Of the 71 *SiWRKY* genes, 65 mapped to 15 sesame linkage groups (LGs), with the exception of LG09. Six genes (*SiWRKY1*, *2*, *5*, *6*, *8*, and *23*) mapped to unanchored scaffolds (Fig. [1](#Fig1){ref-type="fig"}, Table [1](#Tab1){ref-type="table"}). LG06 contained the greatest number of sesame *WRKY* genes (10, 15.38%), whereas LG07, LG10, LG13, and LG16 contained only one gene each.Fig. 1Distribution of *SiWRKY* genes within the sesame linkage group (LG). Vertical bars represent the LGs within the sesame genome. The LG number is indicated at the top of each LG. The scale on the right is in 1 million bases (Mb)

Syntenic analysis and comparison with the grapevine genome revealed that the sesame genome was duplicated in its entirety approximately 71 million years ago, creating two syntenic subgenomes \[[@CR26]\]. Based on the synteny of these subgenomes, we identified 10 pairs of duplicated sesame *WRKY* genes (Additional file [1](#MOESM1){ref-type="media"}, Table [2](#Tab2){ref-type="table"}). We were unable to identify any sesame *WRKY* genes using datasets for tandemly duplicated genes obtained from PTGBase \[[@CR28]\], indicating that the *WRKY* gene family did not undergo tandem gene duplication; this finding is consistent with a previous report \[[@CR29]\]. These results indicate that the *WRKY* gene family underwent whole genome duplication, without tandem gene duplication events (Table [2](#Tab2){ref-type="table"}).Table 2Genome-wide duplication of *SiWRKY* genes*Grapevine*Subgenome1Subgenome2GSVIVT01008046001*SiWRKY16SiWRKY45*GSVIVT01010525001*SiWRKY69SiWRKY49*GSVIVT01014854001*SiWRKY7SiWRKY39*GSVIVT01021397001*SiWRKY55SiWRKY3*GSVIVT01026965001*SiWRKY37SiWRKY14*GSVIVT01027069001*SiWRKY38SiWRKY15*GSVIVT01030258001*SiWRKY28SiWRKY9*GSVIVT01033063001*SiWRKY11SiWRKY10*GSVIVT01033188001*SiWRKY63SiWRKY58*GSVIVT01035426001*SiWRKY60SiWRKY59*

Classification and phylogenetic analysis of the *SiWRKY* genes {#Sec5}
--------------------------------------------------------------

We performed multiple sequence alignments to examine the structural features of each *SiWRKY* protein (Additional file [2](#MOESM2){ref-type="media"}). The results showed that 69 *SiWRKY* proteins contained one or two identical WRKYGQK domains. Although the WRKYGQK domain is highly conserved in WRKY proteins, *SiWRKY59* and *SiWRKY60* differed at one residue, with a glutamine being replaced by a lysine residue; this change is also found in *WRKYs* from tomato, *Arabidopsis*, and other plant species \[[@CR13], [@CR15], [@CR19], [@CR22]\]. Additionally, most of the *SiWRKY* proteins contained the C-X~4--7~-C-X~23~-H motif that forms the C~2~H~2~/C~2~HC-type zinc-finger structure.

A phylogenetic tree was constructed using the neighbor-joining (NJ) method and based on multiple alignments of sesame and *Arabidopsis* WRKY domain aa sequences \[[@CR15]\]. As shown in Fig. [2](#Fig2){ref-type="fig"}, the 71 *SiWRKYs* were classified into three groups (I, II, and III), and the *WRKYs*in Group II were further subdivided into five subgroups (IIa--e). Groups I, II, and III consisted of 12, 48, and seven *SiWRKY* proteins, respectively. A total of four, 11, 18, seven, and eight proteins were assigned to subgroups IIa, IIb, IIc, IId, and IIe, respectively.Fig. 2Phylogenetic analysis of the WRKY proteins in sesame and *Arabidopsis.* Multiple sequence alignments of WRKY amino-acid sequences were performed using ClustalX, and the phylogenetic tree was constructed using MEGA5 by the neighbor-joining (NJ) method and 1000 bootstrap replicates. The tree was divided into seven phylogenetic subgroups, designated I, IIa--e, and III. The bootstrap values were ≥50%

Conserved motifs and structure of the *SiWRKY* family genes {#Sec6}
-----------------------------------------------------------

Using the *SiWRKY* phylogenetic relationships data, we identified structural features of the sesame *WRKYs*, including conserved motifs and the locations of exons and introns. Using Multiple Em for Motif Elicitation (MEME) and InterPro Scan 5, we identified 10 conserved motifs in the sesame *WRKYs* (Fig. [3](#Fig3){ref-type="fig"}, Additional file [3](#MOESM3){ref-type="media"}) \[[@CR30]\]. Motifs 1 and 4 were annotated as WRKY DNA-binding motifs, which is the fundamental characteristic of WRKY proteins. The motif 4 region sequence is conserved in N-terminal WRKY domains. All *SiWRKYs* contained at least one of these motifs, indicating the existence of features conserved in the *WRKY* gene family among the sesame WRKYs identified in this study. Group I proteins had two WRKY domains, each consisting of the conserved aa sequence WRKYGQK and a novel zinc-finger-like motif \[[@CR15]\]. Group I might include the original genes from the other groups \[[@CR30]\]. The gene structure predictions (Fig. [4](#Fig4){ref-type="fig"}) revealed that the *SiWRKY* genes had between one (*SiWRKY14*, *15*, *10*, *11*, *49*, *69*, *42*, and *6*) and 11 (*SiWRKY29*) introns.Fig. 3Conserved motifs of the *SiWRKY* proteins arranged according to their phylogenetic relationships. The NJ tree was constructed from the amino-acid sequences of sesame *WRKYs* using ClustalX and MEGA5 with 1000 bootstrap replicates. The conserved motifs in the *SiWRKY* proteins were identified using Multiple Em for Motif Elicitation (MEME). In total, 10 motifs were identified and are shown in different colors. Motif locations are also indicated Fig. 4Structures of the 71 *SiWRKY* genes arranged in families. The NJ tree was constructed from the amino-acid sequences of the sesame *WRKYs* using ClustalX and MEGA5 with 1000 bootstrap replicates. Structural analyses of the *SiWRKY* genes were performed using the gene structure display server. The exons and introns are represented by colored boxes and black lines, respectively

Tissue-specific expression profiling of the *SiWRKY* genes {#Sec7}
----------------------------------------------------------

To generate expression profiles of the *SiWRKY* genes under normal conditions, RNA sequence transcriptome data were collected and analyzed. The expression levels of the 71 *SiWRKY* genes were obtained on the basis of the reads per kilobase of transcript per million mapped reads (RPKM) values from six tissue samples (roots, stem, flowers, leaves, capsules, and seeds). The RPKM values of the transcripts were clustered hierarchically and displayed in a heat map (Fig. [5](#Fig5){ref-type="fig"}).Fig. 5Expression profile analysis of *SiWRKY* genes in different sesame tissues. Transcriptome data (Reads Per Kilobase per Million mapped reads; RPKM) were used to measure the expression levels of *SiWRKY* genes in roots, stem, flowers, leaves, capsules, and seeds. The colored scale for the different expression levels is shown

Quantification of transcript levels expressed in different tissues can be useful in determining gene function. The *SiWRKY* genes displayed diverse expression patterns, possibly reflecting the distinct roles of the different gene family members. A total of 66.20% (47/71), 30.99% (22/71), 32.39% (23/71), 21.13% (15/71), 43.66% (31/71), and 21.13% (15/71) of the *SiWRKY* genes were highly expressed (values \>1) in roots, stem, flowers, leaves, capsules, and seeds, respectively. Most of the *SiWRKY* genes were expressed in all tissues, although *SiWRKY4*, *SiWRKY32*, *SiWRKY61*, *SiWRKY62*, and *SiWRKY70* were only expressed at low levels. Additionally, six *SiWRKY* genes (*SiWRKY22* and *SiWRKY39* in Group I, *SiWRKY16* and *SiWRKY21* in Group IIb, *SiWRKY58* in Group IId, and *SiWRKY29* in Group IIc) were continuously expressed at high levels (values \>1) in all six organs, suggesting that these genes may be important for plant growth and organ development.

Expression patterns of *SiWRKY*s in response to waterlogging and drought stresses {#Sec8}
---------------------------------------------------------------------------------

The expression of *WRKY* genes has been examined under different stress conditions, including high salinity, drought, and high temperature; however, plant gene expression in response to waterlogging stress has not been studied extensively \[[@CR13]\]. In this study, we investigated the expression of *SiWRKY* genes in the roots of sesame cultivars that were tolerant or sensitive to waterlogging stress using quantitative real-time polymerase chain reaction (qRT-PCR). As shown in Fig. [6](#Fig6){ref-type="fig"}, the majority of the *SiWRKY* genes (42 in the tolerant and 40 in the sensitive cultivar) were upregulated in both tolerant and sensitive waterlogged cultivars. Among these upregulated genes, \>2-fold increases in expression (*P* \< 0.01) were observed in 26 of the waterlogging-tolerant cultivars and 22 of the waterlogging-sensitive cultivars. Moreover, the same 18 *SiWRKY* genes (*SiWRKY8*, *SiWRKY13*, *SiWRKY16*, *SiWRKY19*, *SiWRKY30*, *SiWRKY35*, *SiWRKY41*, *SiWRKY43*, *SiWRKY46*, *SiWRKY49*, *SiWRKY51*, *SiWRKY54*, SiWRKY55, *SiWRKY56*, *SiWRKY64*, *SiWRKY66*, *SiWRKY68*, and *SiWRKY71*) displayed \>2-fold increases in expression level in both waterlogging-tolerant and -sensitive cultivars. In addition, the *SiWRKY68* gene exhibited the highest expression level, with \>10-fold increases in both waterlogging-tolerant and -sensitive cultivars. Waterlogging also decreased the transcript abundance of a large number of *SiWRKY* genes in roots. In total, 30 (42.3%) and 29 (40.8%) *SiWRKY* genes exhibited \>2-fold downregulation (*P* \< 0.01) in waterlogging-sensitive and -tolerant cultivars, respectively. In particular, the same 15 *SiWRKY* genes (*SiWRKY1*, *SiWRKY6*, *SiWRKY7*, *SiWRKY12*, *SiWRKY17*, *SiWRKY27*, *SiWRKY39*, *SiWRKY42*, *SiWRKY47*, *SiWRKY57*, *SiWRKY59*, *SiWRKY60*, *SiWRKY62*, *SiWRKY63*, and *SiWRKY70*) displayed exhibited \>2-fold downregulation in both waterlogging-sensitive and -tolerant cultivars. Therefore, our results show differential expression (\>2-fold upregulation or downregulation by a factor of two or more) of 33 genes in both waterlogging-sensitive and -tolerant cultivars, suggesting that these sesame genes play important roles in the response to waterlogging. Additionally, 27 of these 33 genes belong to *SiWRKY* gene Group II, while only three (*SiWRKY7*, *SiWRKY30*, and *SiWRKY39*), one (*SiWRKY27*), and two (*SiWRKY1* and *SiWRKY17*) belong to Group I, Group III, and the unknown group, respectively. Of the 18 upregulated genes, only *SiWRKY30* did not belong to Group II. In addition, 29 significant different expression *SiWRKYs* were found between waterlogging-tolerant and -sensitive cultivars, and 18 of them had high different expression level (\>2 or \<−2) (Additional file [4](#MOESM4){ref-type="media"}).Fig. 6*SiWRKY* gene expression in sesame roots treated for 8 h with waterlogging stress compared with untreated controls. Transcript abundance was quantified using quantitative real-time polymerase chain reaction (qRT-PCR) and expression levels were normalized using sesame *β*-actin (SIN_1009011) as a reference gene. The mean expression levels from three independent biological replicates were analyzed for significance using *t*-tests (*p* \< 0.01). The histograms represent the relative expression levels and rates of gene induction (stress/control). An asterisk (\*) indicates a significant (2-fold) increase in gene expression in plants treated with waterlogging stress compared with untreated controls

As shown in Fig. [7](#Fig7){ref-type="fig"}, drought stress decreased *SiWRKY* gene expression in sesame roots. Most of the *SiWRKY* genes were downregulated in both types of cultivars (53 drought-tolerant and 51 drought-sensitive cultivars). More genes were downregulated by \>2-fold among the drought-tolerant (32 genes) than in the drought-sensitive (19 genes) cultivars. In contrast, 20 and 18 *SiWRKY* genes were upregulated in the drought-sensitive and drought-resistant sesame cultivars, respectively. The expression of five genes (*SiWRKY11*, *SiWRKY33*, *SiWRKY49*, *SiWRKY55*, and *SiWRKY59*) was upregulated by \>2-fold and the expression of 19 genes (*SiWRKY5*, *SiWRKY6*, *SiWRKY8*, *SiWRKY16*, *SiWRKY17*, *SiWRKY21*, *SiWRKY24*, *SiWRKY26*, *SiWRKY27*, *SiWRKY32*, *SiWRKY38*, *SiWRKY40*, *SiWRKY43*, *SiWRKY47*, *SiWRKY48*, *SiWRKY56*, *SiWRKY57*, *SiWRKY62*, and *SiWRKY70*) was downregulated by \>2-fold. Two *SiWRKY* genes (*SiWRKY42* and *SiWRKY61*) displayed increases in expression by \>2-fold in the drought-sensitive cultivar, while decreases in expression by \>2-fold were detected in the drought-tolerant cultivar. These 26 sesame genes displaying marked changes in expression might play important roles in drought stress responses. Eighteen of these genes belonged to Group II, one to Group I, six to Group III, and one to the unknown gene group. Additionaly, 33 significant different expression *SiWRKYs* were found between drought-tolerant and -sensitive cultivars, and 26 of them had high different expression level (\>2 or \<−2) (Additional file [5](#MOESM5){ref-type="media"}).Fig. 7*SiWRKY* gene expression in sesame roots treated for 11 d with drought stress compared with untreated controls. Transcripts abundance was quantified using qRT-PCR and expression levels were normalized using sesame *β*-actin (SIN_1009011) as a reference gene. The mean expression levels from three independent biological replicates were analyzed for significance using *t*-tests (*p* \< 0.01). The histograms represent the relative expression levels and rates of gene induction (stress/control). An asterisk (\*) indicates a significant (2-fold) increase in gene expression in plants treated with drought stress compared with untreated controls

Expression of selected *SiWRKY* genes in response to waterlogging and drought stresses {#Sec9}
--------------------------------------------------------------------------------------

To confirm the identities of some of the genes important for waterlogging- and drought tolerance, 12 differential expression *SiWRKY* genes were selected and their expression levels quantified by qRT-PCR at different time-points after the onset of each abiotic stress. As shown in Fig. [8](#Fig8){ref-type="fig"}, six *SiWRKY* genes were expressed at different times following the start of the waterlogging treatment in both the waterlogging-tolerant and -sensitive cultivars (*P* \< 0.05). The expression levels of *SiWRKY13*, *SiWRKY35*, and *SiWRKY43* increased during the waterlogging treatment, although the expression of each gene peaked at a different time. The peak expression of *SiWRKY35* occurred before that of *SiWRKY13* and *SiWRKY43*. In contrast, the expression of *SiWRKY17*, *SiWRKY59*, and *SiWRKY63* was downregulated by waterlogging. For the waterlogging -tolerant and --sensitive cultivars, the difference of the expression of these six genes mainly appeared at 36 h. We noticed that *SiWRKY35*, *SiWRKY43*, *SiWRKY59* and *SiWRKY63* expressed in a higher level in waterlogging sensitive cultivar than that in waterlogging tolerant cultivar, especially for the *SiWRKY63*. This result suggested that the different of the expression of *SiWRKY* genes might be one of the reasons for the tolerance of waterlogging for sesame varieties.Fig. 8Expression profiles of six highly-expressed *SiWRKY* genes in response to waterlogging stress. The relative expression levels of six *SiWRKY* genes were measured from plants treated with waterlogging for 0, 8, 16, and 36 h, and also at 12 h after water was withdrawn (WD12h) from plants waterlogged for 36 h. The different columns represent different cultivars: a waterlogging-tolerant cultivar (WT) and a waterlogging-sensitive cultivar (WS). Three independent replicates were used to generate each expression value. The error bars represent standard deviations. Values with the same letter were not significantly different when assessed using Duncan's multiple range test (*p* \< 0.05, *n* = 3)

As shown in Fig. [9](#Fig9){ref-type="fig"}, the expression of *SiWRKY* genes in response to drought differed significantly between drought-tolerant and -sensitive cultivars (*P* \< 0.05). The expression of *SiWRKY6*, *SiWRKY11*, *SiWRKY42*, *SiWRKY55*, and *SiWRKY59* was considerably increased by drought stress. However, the expression of *SiWRKY6* was suppressed by severe drought conditions (5% soil water content) and recovered following re-watering (REW), indicating that severe drought decreases the expression of some sesame *WRKY* genes. Gene expression of *SiWRKY* under drought stress also showed significant different between drought- tolerant and - sensitive cultivars. Similar to the waterlogging stress, gene expression of most *SiWRKY* genes in the sensitive cultivar were higher than that in tolerant cultivar. The time-point that the largest difference of the gene expression appeared varied in each gene. For example, *SiWRKY11* had a much higher gene expression level in 5% soil water content of drought-tolerant cultivars, while *SiWRKY42* expressed highest in 15% soil water content.Fig. 9Expression profiles of six highly-expressed *SiWRKY* genes in response to drought stress. The relative expression levels of six *SiWRKY* genes were measured from plants treated with drought stress and harvested at a soil water content of 35%, 15%, 10%, 5%, and 35% after re-watering (REW35%). The different columns represent different cultivars: a drought-tolerant cultivar (DT) and a drought-sensitive cultivar (DS). Three independent replicates were used to generate each expression value. The error bars represent standard deviations. Values with the same letter were not significantly different when assessed using Duncan's multiple range test (*p* \< 0.05, *n* = 3)

Discussion {#Sec10}
==========

Number and type of sesame *WRKY* genes {#Sec11}
--------------------------------------

The WRKY transcription factor family is one of the most important gene families involved in plant development and stress responses, and *WRKY* genes have been identified in many species, including *Arabidopsis*, rice, grape, maize, and cucumber \[[@CR15], [@CR20], [@CR31]--[@CR33]\]. Table [3](#Tab3){ref-type="table"} summarizes the numbers and types of *WRKY* genes found in higher plants, and illustrates their diversity among species that have had their genomes sequenced; the number of genes ranges from 55 in cucumber to 343 in rapeseed \[[@CR34]\]. In this study, we identified 71 *WRKY* genes from a total of 27,148 annotated genes in the sesame genome. Relative to the genome size, the sesame *WRKY* gene family (350 Mb, 71 *WRKY* genes) is large compared with that of grape, cucumber, and castor bean. However, it is small compared with the *WRKY* gene families of *Arabidopsis* (107 Mb, 72 *WRKY* genes) and rice (440 Mb, 103 *WRKY* genes). Table [3](#Tab3){ref-type="table"} shows that one key difference between the sesame, *Arabidopsis*, and rice genomes is the number of Group III *WRKY* genes in each. The much more numerous Group III *WRKY* genes in *Arabidopsis* and rice is explained by tandem duplication and recent duplication events, which has led to a large-scale expansion of the gene families in these genomes \[[@CR31], [@CR35]\]. Recent gene duplication and tandem duplication events are the most important factors in the rapid expansion and evolution of gene families \[[@CR35]\]. Previous research has demonstrated that the *Arabidopsis* Group III *WRKY* gene family expanded rapidly as a result of recent segmental- and tandem duplication events. Additionally, all of the tandemly duplicated *WRKY* genes in *Arabidopsis* belong to Group III, whereas we identified no segmentally or tandemly duplicated Group III *WRKY* genes in sesame. In this study, we identified 10 pairs of segmentally duplicated *SiWRKY* genes, but none of these belonged to Group III. Additionally, no *SiWRKY* genes had been generated by tandem duplication events in the sesame genome. Therefore, the small size of Group III in sesame is probably due to the absence of *WRKY* gene tandem duplication events.Table 3Numbers and types of *WRKY* genes in higher plantsPlantSpeciesGenome sizeNameTotalGroupNGIIIaIIbIIcIIdIIeIIIsesame*Sesamum indicum*350 Mb*SiWRKY*7112411187874cucumber*Cucumis sativus*367 Mb*CsWRKY*551044168760*ArabidopsisArabidopsis thaliana*107 Mb*AtWRKY*7213471879140grape*Vitis vinifera*490 Mb*VvWRKY*591238157662rice*Oryza sativa*440 Mb*OsWRKY*103154815711360tomato*Solanum lycopersicum*900 Mb*SlWRKY*78155816617113flax*Linum usitatissimum*373 Mb*LuWRKY*9724413161112152soybean*Glycine max*1.1 Gb*GmWRKY*188321433422120260Castor bean*Ricinus communis*350 Mb*CbWRKY*479310123550*Brachypodium distachyonBrachypodium distachyon*272 Mb*BdWRKY*86153621610232maize*Zea mays*2.3 Gb*ZmWRKY*13627711291417310cotton*Gossypium raimondii*761 Mb*GrWRKY*11622616331512120rapeseed*Brassica napus*630 Mb*BnWRKY*34312111345528305113barley*Hordeum vulgare*5.1 Gb*HvWRKY*458411153130pear*Pyrus bretschneideri*527 Mb*PbWRKY*10317610241516150

Conserved motifs and structures of the sesame *WRKY* genes {#Sec12}
----------------------------------------------------------

Almost all of the *SiWRKYs* contained the WRKYGQK domain, although two Group IIc WRKYs (*SiWRKY59* and *SiWRKY60*) contained the WRKYGKK domain. This variant of the WRKY domain has also been found in pepper \[[@CR9]\], tea \[[@CR36]\], and apple proteins \[[@CR13]\]. Waterlogging induced the expression of *SiWRKY59* and *SiWRKY60* (Fig. [6](#Fig6){ref-type="fig"}), indicating that these genes may be involved in sesame abiotic stress responses.

The conserved motifs and structural features of the sesame *WRKYs* were identified using MEME and InterPro Scan 5. The InterPro Scan 5 analysis suggested that *SiWRKY8* and *SiWRKY35* are PWRKY transcription factors, which represented a subset of Group IIc probable WRKY transcription factors from plants. PWRKY transcription factors, with the InterPro number IPR017396, were known to regulate various abiotic stress responses \[[@CR37]\]. Thus, *SiWRKY8* and *SiWRKY35* might regulate the biotic and abiotic stress responses. This is consistent with the high expression of these genes in response to waterlogging. The *SiWRKY12*, *SiWRKY13*, *SiWRKY24*, *SiWRKY31*, *SiWRKY36*, *SiWRKY58*, and *SiWRKY63* genes contained Zn-cluster-domain sequences (IPR018872) and encoded WRKY-GCM1 zinc-finger-domain proteins, indicating that these genes acquired their functional diversity as developmental and regulatory genes \[[@CR38]\]. *SiWRKY29* encodes an ATP-dependent metallopeptidase that belongs to the FtsH (IPR005936) protein family \[[@CR39]\].

The exon/intron structural diversity found among the *SiWRKY* genes is related to their evolution \[[@CR40]\]. An exon/intron distribution analysis demonstrated that *WRKY* genes from sesame had greater structural diversity than those from *Populus trichocarpa* or cassava. Most of the *SiWRKY* genes (33/71) had two introns, which is common in other plants, including *Pyrus bretschneideri* (59/103), *Populus trichocarpa* (49/104), cassava (42/85), and physic nut (30/58) \[[@CR4], [@CR21], [@CR30], [@CR41]\]. Most of the sesame MADS-box genes (63.2%, 36/57) have 0 or 1 intron, but most of *SiWRKY* genes (81.7%, 58/71) have 2 to 4 introns \[[@CR25]\]. The sesame *WRKY* genes had significantly more introns than the sesame MADS-box genes, indicating that the WRKY gene structure in sesame is more complex.

Diverse expression patterns of *SiWRKY* genes in different tissues {#Sec13}
------------------------------------------------------------------

We analyzed the expression of *SiWRKY* genes in six different tissues. The results demonstrated variation in the expression patterns of *SiWRKY* genes. Most *SiWRKY* genes were highly expressed in roots, whereas a few *SiWRKY* genes were expressed in developing seeds. This is consistent with observations made in other plants, including rice \[[@CR18]\], cucumber \[[@CR33]\], grape \[[@CR42]\], apple \[[@CR13]\], cassava \[[@CR30]\], cotton \[[@CR43]\], physic nut \[[@CR41]\], and cabbage \[[@CR44]\]. Our results revealed that *SiWRKY* genes are expressed tissue-specifically and the high expression levels observed in roots might reflect their roles in responses to abiotic and biotic stresses that first affect plants below ground.

In total, 15 *SiWRKY* genes were highly expressed in at least five sesame tissues. Six of these genes (*SiWRKY16*, *SiWRKY21*, *SiWRKY22*, *SiWRKY29*, *SiWRKY39*, and *SiWRKY58*) were highly expressed in all sesame tissues. Highly expressed *WRKY* genes usually play important roles in plant development \[[@CR45]\]. Therefore, we concluded that the 15 highly expressed *SiWRKY* genes might be important regulatory factors in sesame development, although further studies are required to verify the function of these genes. Most of these highly expressed *SiWRKY* genes belong to Groups I and IId. Previous research has demonstrated that Group I *WRKY* genes are ancestral to other *WRKY* genes in plants and are more likely to be constitutively expressed in different tissues \[[@CR46]\]. For example, the Group I genes *SiWRKY28*, *SiWRKY29*, and *SiWRKY67* are expressed in most sesame tissues and highly expressed in response to waterlogging and drought stresses.

In contrast, 17 *SiWRKY* genes were expressed at low levels in all sesame tissues, and 19 *SiWRKY* genes were specifically expressed in only one tissue. Among the specifically expressed *SiWRKY* genes, *SiWRKY51* was expressed only in capsules and the remaining *SiWRKY* genes were expressed only in roots. These specifically or minimally expressed *SiWRKY* genes were found in all the *WRKY* gene subgroups, but many were found in Groups IIc and IIe. A number of Group IIc *WRKY* genes in *Arabidopsis* (e.g., *AtWRKY8*, *AtWRKY48*, *AtWRKY50*, and *AtWRKY57*) are involved in responses to bacterial and fungal pathogens, and in the jasmonic acid- and salicylic acid-mediated signaling pathways \[[@CR27]\]. Therefore, although some Group IIc *SiWRKY* genes were expressed at low levels in most sesame tissues, they may play key roles in responses to biotic and abiotic stresses. In this study, *SiWRKY51* and *SiWRKY65* were highly expressed in the roots of waterlogged plants, whereas *SiWRKY10* and *SiWRKY53* were highly expressed in response to drought stress. These results indicate that some *SiWRKY* genes might only be expressed in response to particular abiotic stresses.

Identification of *SiWRKY* genes involved in responses to abiotic stresses {#Sec14}
--------------------------------------------------------------------------

Waterlogging and drought are the most serious abiotic stresses for sesame and result in significant losses (20%--50%) in sesame production within China \[[@CR24], [@CR47]\]. However, few abiotic stress tolerance genes have been identified in sesame. Recent research has demonstrated that *WRKY* genes are involved in responses to various stresses and there is now compelling evidence that *WRKYs* are plant transcription factors that regulate tolerance to abiotic stresses \[[@CR48]\]. Gene expression studies have shown that 20 *AtWRKY* genes in *Arabidopsis*, 41 *OsWRKY* genes in rice, 66 *GmWRKY* genes in soybean, 41 *BrWRKY* genes in *Brassica rapa*, and 74 *BnWRKY* genes in rapeseed are involved in responses to abiotic stresses \[[@CR14], [@CR18], [@CR34], [@CR49]--[@CR51]\]. In this study, 44 *SiWRKY* genes were expressed differentially in response to waterlogging and drought stresses, indicating that these genes may also be involved in responses to abiotic stresses. To identify the *WRKY* genes that regulate tolerance to abiotic stresses in sesame, waterlogging- and drought-tolerant and sensitive cultivars were investigated. As shown in Figs. [6](#Fig6){ref-type="fig"} and [7](#Fig7){ref-type="fig"}, the expression of some *SiWRKY* genes differed significantly between the tolerant and sensitive sesame cultivars. For example, *SiWRKY10* was highly expressed in tolerant cultivars in response to waterlogging. Further analysis showed that responses to abiotic stresses occurred at different time-points. *SiWRKY17* and *SiWRKY59* were highly expressed after 8 h of waterlogging, whereas the expression of *SiWRKY13* and *SiWRKY43* peaked at 36 h after waterlogging began. This suggests that these *SiWRKY* genes might play important regulatory roles in sesame abiotic stress tolerance and may act at different stages of the stress response.

Compared with the *WRKY* genes that involved in the response of drought, cold and heat stresses, few *WRKY* genes that responded to waterlogging stress have been identified in plant. In addition, the expression pattern of *WRKY* genes under waterlogging stress was also unclear. Therefore, the expression of *SiWRKY* genes under waterlogging stress and the expression level of six highly expressed *SiWRKY* genes during the waterlogging treatment were detected in the present study. The qPCR results showed that 33 *SiWRKY* genes either increase or decrease their expression by a factor of two or more in both waterlogging tolerant and sensitive cultivars. With increasing treatment duration, both up-regulation and down-regulation of the *SiWRKY* genes were found. Previous studies have revealed that one *WRKY*gene could function in several disparate signaling pathways. For example, *AtWRKY70* functions in plant growth, drought response and cell death \[[@CR52], [@CR53]\]. Interestingly, we found that the *SiWRKY17* and *SiWRKY43* was induced to highly express under waterlogging stress, while the orthologous genes of them, *AtWRKY32* and *AtWRKY29*, was reported to respond to UV irradiation and heat stress, respectively. The result indicated that the orthologous *WRKY* genes might mediate different pathways and play different roles under the abiotic stress response in different species.

Both phylogeny-based and BLAST-based methods were used to identify *WRKY* gene orthologs in comparisons of sesame and *Arabidopsis*. A phylogenetic tree based on WRKY protein sequences from sesame and *Arabidopsis* was constructed and nodes with bootstrap values \>50 were used to identify possible orthologs. In addition, standard BLASTP searches were applied to verify possible orthologs of *WRKY* genes in sesame and *Arabidopsis*, with relatively strict criteria (Additional file [6](#MOESM6){ref-type="media"}). In total, we found 10 orthologous pairs shared between sesame and *Arabidopsis* (Table [4](#Tab4){ref-type="table"}). The functions of the 10 *AtWRKY* gene orthologs have been determined and all are involved in responses to abiotic stresses in *Arabidopsis*. *AtWRKY13* and *AtWRKY20* regulate tolerance to drought stress \[[@CR54], [@CR55]\] and their sesame orthologs, *SiWRKY6* and *SiWRKY40*, are also highly expressed in response to drought. Therefore, we conclude that the *SiWRKY* orthologs of *AtWRKY* genes also play key roles in the tolerance of abiotic stresses in sesame.Table 4Orthologous *WRKY* genes in sesame and *Arabidopsis*Sesame*Arabidopsis*FunctionReference*SiWRKY6AtWRKY13*Drought\[[@CR54]\]*SiWRKY17AtWRKY32*UV irradiation, heavy metals\[[@CR62]\]*SiWRKY29AtWRKY49*H~2~O~2~\[[@CR63]\]*SiWRKY30AtWRKY44*drought\[[@CR64]\]*SiWRKY33AtWRKY55*Oxidative stress\[[@CR65]\]*SiWRKY35AtWRKY23*H~2~O~2~, ABA, mannitol\[[@CR66]\]*SiWRKY40AtWRKY20*Drought\[[@CR55]\]*SiWRKY42AtWRKY43*Nitrogen\[[@CR67]\]*SiWRKY43AtWRKY29*Heat\[[@CR68]\]*SiWRKY67AtWRKY26*Cold\[[@CR69]\]

Functional divergence and segmental duplication of *WRKY* genes {#Sec15}
---------------------------------------------------------------

Recent segmental duplication has occurred frequently in plant genomes because most plants are diploidized polyploids, and many duplicated chromosomal blocks have been retained \[[@CR28]\]. Comparisons with grape suggest that the sesame genome underwent a recent genome duplication event, approximately 71 ± 19 million years ago \[[@CR26]\]. In total, 10 pairs of sesame *WRKY* genes were identified as segmentally duplicated. In the *Arabidopsis* genome, the Group III *WRKY* genes are located in a recently segmentally duplicated region of the genome and are highly expressed in response to abiotic stresses. However, in this study, there was little correlation among the expression patterns of the duplicated sesame *WRKY* genes. For example, *SiWRKY16* was highly expressed in all sesame tissues, whereas the duplicated *SiWRKY45* gene was only expressed in roots. Additionally, the expression of *SiWRKY55* in response to drought was much higher than that of *SiWRKY3*.

Research using *Arabidopsis*, rice, and soybean has focused on identifying the gene targets of *WRKYs* and understanding the associated regulatory networks. One study showed that co-regulated networks involving *WRKY* genes were important in regulating the responses of pak-choi to a variety of abiotic stresses \[[@CR27]\]. Wheat *TaWRKY19* regulates the expression of *DREB2A*, which encodes a key transcription factor that controls the expression of drought-related genes \[[@CR56]\]. Therefore, the regulatory roles of *WRKY* genes in response to abiotic stresses are complex and further studies are required to understand their functions in sesame.

Conclusions {#Sec16}
===========

In this study, we identified a total of 71 sesame *WRKY* genes and focused on those involved in responses to waterlogging and drought stresses. The distribution, classification, gene structure, and evolutionary characteristics of the sesame *WRKY* genes were investigated. The differential expression patterns of *SiWRKY* genes in the tissues of selected cultivars showed that these genes play different roles in sesame development and many exhibit tissue-specific expression patterns. Additionally, *SiWRKY* gene expression analyses revealed that some were markedly upregulated or downregulated in response to waterlogging and drought stresses. Our results also revealed significant differences in the abiotic-stress-induced expression of *WRKYs* in stress-sensitive and -tolerant sesame cultivars, indicating the involvement of these *WRKY* genes in abiotic stress tolerance in sesame. In conclusion, our study establishes a structural and functional framework to investigate sesame WRKY proteins. Although the sesame genome was sequenced several years ago, the identification of sesame abiotic-stress-related genes and investigations into their functions are still at an early stage. Our results will facilitate further studies into the functions of *WRKY* genes important in responses to abiotic stresses and the development of molecular breeding programs to enhance abiotic stress tolerance in sesame.

Methods {#Sec17}
=======

Identification of the *WRKY* gene family in sesame {#Sec18}
--------------------------------------------------

All sesame protein sequences were obtained from the sesame genome database ([http://ocri-genomics.org/Sinbase/](http://ocri-genomics.org/Sinbase)) \[[@CR26]\]. The *Arabidopsis thaliana AtWRKY* gene sequences were downloaded from UniProt (<http://www.uniprot.org>/). The HMM profile for the WRKY DNA-binding domain (PF03106) was downloaded from the PFAM protein families database (<http://pfam.xfam.org>) and used to identify *WRKY* genes from the sesame genome with HMMER 3.0 (<http://hmmer.janelia.org>/). BLAST analyses with all the *Arabidopsis WRKYs* were used to check the predicted *WRKYs* from the sesame database. The CDD ([http://www.ncbi.nlm.nih.gov/cdd/](http://www.ncbi.nlm.nih.gov/cdd)) and PFAM databases (<http://pfam.xfam.org>/) were used to validate all the potential sesame *WRKY* genes identified by HMM and BLAST if they contained a WRKY domain. Multiple sequence alignments were used to confirm the conserved domains from the predicted WRKY sequences.

Chromosomal location and phylogenetic analysis of the *WRKY* gene family in sesame {#Sec19}
----------------------------------------------------------------------------------

The physical positions of the *SiWRKY* genes were established using Sinbase ([http://ocri-genomics.org/Sinbase/](http://ocri-genomics.org/Sinbase)) and mapped to 16 LGs in the sesame genome using MapChart 2.3 \[[@CR57]\]. Additionally, Clustal X 2.1 and MEGA 5.2 were used to construct a NJ phylogenetic tree based on the aa sequences of the sesame WRKY domains and selected *Arabidopsis WRKYs*, with 1000 bootstrap replicates. An alignment of sesame WRKY domains is shown in Additional file [1](#MOESM1){ref-type="media"}.

Protein properties and sequence analysis {#Sec20}
----------------------------------------

Protein MWs and isoelectric points (pIs) were predicted using the ProtParam program (ExPASy tools) based on their deduced aa sequences. The conserved motifs in the full-length WRKY proteins were identified using the MEME program (<http://alternate.meme-suite.org/tools/meme>). The parameters employed in the analysis were as follows: maximum number of motifs = 10; optimum width of motifs = 15--50 \[[@CR30]\]. Additionally, all of the identified motifs were annotated using InterProScan (<http://www.ebi.ac.uk/interpro/search/sequence-search>). The exon/intron structures of the *SiWRKY* genes were determined by comparing their predicted coding sequence (CDS) with genomic sequences using the gene structure display server web-based bioinformatics tool (<http://gsds.cbi.pku.edu.cn>/) \[[@CR58]\].

Analysis of *SiWRKY* gene expression in different organs using transcriptomic data {#Sec21}
----------------------------------------------------------------------------------

Total RNA was extracted from roots, shoots, leaves, seed capsules, and seeds of Zhongzhi No. 13 grown under normal conditions. RNA pools were constructed using 3 μg of RNA per sample according to the manufacturer's instructions and sequenced on a Gene Analyzer II system (Illumina, Inc., San Diego, CA, USA) according to the Illumina RNA-seq protocol. Gene expression levels were calculated in RPKM by taking into account the length of each gene and the read counts mapped. The sesame *WRKY* gene expression pattern analyses were performed using Gene Cluster 3.0, and the RPKM values for each gene in all tissue samples were log10 transformed. Finally, a heat map was generated using TreeView 1.0.4 \[[@CR59]\].

Plant materials and treatments {#Sec22}
------------------------------

Waterlogging-tolerant (WT) cultivar 2541, waterlogging-sensitive (WS) cultivar 4508, drought-tolerant (DT) cultivar 0635, and drought-sensitive (DS) cultivar 4728 were all selected from sesame germplasm provided by the Oil Crops Research Institute, Chinese Academy of Agricultural Sciences, Wuhan, China.

For the waterlogging treatment, sesame plants at anthesis were irrigated until the soil surface was covered by a thin layer of water and this was maintained for 36 h. The plants were harvested 8 h later and their roots were immediately frozen in liquid nitrogen and stored at −80 °C prior to further analysis. Control plants were harvested 15 h before the waterlogging treatment. To further investigate waterlogging resistance-related genes, we harvested plants that had been waterlogged for 0, 8, 16, and 36 h, and also 12 h after water was withdrawn (WD12h) from plants waterlogged for 36 h. The roots were harvested as described previously.

For the drought stress treatment, water was withheld for 11 d from sesame plants at anthesis. Root samples were collected immediately thereafter and frozen in liquid nitrogen prior to analysis \[[@CR24]\]. To further investigate genes important for drought resistance, we harvested plants when the soil water content was reached 35% (0 d controls), 15% (3 d), 10% (7 d), 5% (11 d), and 35% (14 d) after REW. The roots were harvested as described previously.

qRT-PCR analyses of *SiWRKY* gene expression in response to waterlogging and drought stresses {#Sec23}
---------------------------------------------------------------------------------------------

RNA was extracted from the roots of each of the four cultivars using the EASYspin Plus Plant RNA Kit (Aidlab Biotechnologies, Beijing, China) \[[@CR60]\] according to the manufacturer's instructions. The RNA was quantified using a BIOMATE 3 spectrophotometer (Thermo Scientific, Worcester, MA, USA) and its integrity was confirmed using 1% agarose gel electrophoresis. A total of 1 mg of RNA was reverse-transcribed into cDNA using the iScript cDNA Synthesis kit (Bio-Rad, Hercules, CA, USA). A control amplicon was generated using the following primers for amplification of *β*-actin (SIN_1009011): forward primer, 5′-TTTGAGCAGGAACTGGACACT-3′, and reverse primer, 5′-ACAACACTTCTGGACAACGGA-3′. Gene expression levels were determined by performing qRT-PCR in triplicate on an Icycler iQ5 (Bio-Rad) using the SYBR Green Supermix kit (Bio-Rad), all according to the manufacturer's instructions. Data were analyzed using iQ5 2.1 software (Bio-Rad) and the 2^--ΔΔCT^ method \[[@CR61]\].

Additional files
================

 {#Sec24}

Additional file 1:Synteny of subgenomes in the sesame. The green bars represent the sesame chromosomes. The numbers 01--16 represent LGs within the sesame genome. Black lines on the green bars indicate the locations of sesame genes within the LGs. Colored lines indicate subgenomes in sesame (PDF 1332 kb) Additional file 2:Alignment of *SiWRKY* domain aa sequences. The alignment was performed using Jalview. The conserved WRKY aa and zinc-finger motifs are highlighted in blue. Gaps are indicated by dashes. (PDF 506 kb) Additional file 3:Conserved motifs of WRKY proteins in sesame. Significant motifs of more than 10 aa in length were predicted using MEME analysis. The motif IDs, consensus sequence lengths in aa, and *e*-value of each predicted motif are shown. (PDF 234 kb) Additional file 4:The cultivar-specific *SiWRKY* gene expression in sesame roots treated for 8 h with waterlogging stress compared with untreated controls in cultivars. Transcript abundance was quantified using quantitative real-time polymerase chain reaction (qRT-PCR) and expression levels were normalized using sesame *β*-actin (SIN_1009011) as a reference gene. The mean expression levels from three independent biological replicates were analyzed for significance using *t*-tests (*p* \< 0.01). The histograms represent the relative expression levels and rates of gene induction (stress/control). An asterisk (\*) indicates a significant (2-fold) increase in gene expression in plants treated with waterlogging stress compared with untreated controls. (PDF 341 kb) Additional file 5:The cultivar-specific *SiWRKY* gene expression in sesame roots treated for 11 d with drought stress compared with untreated controls in cultivars. Transcripts abundance was quantified using qRT-PCR and expression levels were normalized using sesame *β*-actin (SIN_1009011) as a reference gene. The mean expression levels from three independent biological replicates were analyzed for significance using *t*-tests (*p* \< 0.01). The histograms represent the relative expression levels and rates of gene induction (stress/control). An asterisk (\*) indicates a significant (2-fold) increase in gene expression in plants treated with drought stress compared with untreated controls. (PDF 544 kb) Additional file 6:Synteny between *SiWRKY* genes in the sesame and *Arabidopsis* genomes. The green bars represent the chromosomes of the two species. The numbers 01--16 represent LGs within the sesame genome and the five *Arabidopsis* chromosomes are labeled Chr1--Chr5. Black lines on the green bars indicate the locations of *SiWRKY* genes on the chromosomes/within the LGs. Colored lines indicate orthologous genes in sesame and *Arabidopsis*. (PDF 425 kb)

Aa

:   amino acid(s)

ATP

:   Adenosine triphosphate

BLAST

:   Basic Local Alignment Search Tool

BLASTP

:   BLAST for protein sequences

CDD

:   Conserved domain database

CDS

:   Coding DNA sequence

DS

:   Drought-sensitive

DT

:   Drought-tolerant

FtsH

:   Filamentation temperature-sensitive H

GSDS

:   Gene structure display server

HMM

:   Hidden Markov model

LGs

:   Linkage groups

MEME

:   Multiple Em for Motif Elicitation

MW

:   Molecular weight

NJ

:   Neighbor-joining

pI

:   isoelectric point

qRT-PCR

:   quantitative real-time polymerase chain reaction

REW

:   After re-watering

RPKM

:   Reads per kilobase of transcript per million mapped reads

TD

:   Tandem duplication

WGD

:   Whole genome duplication

WS

:   Waterlogging-sensitive

WT

:   Waterlogging-tolerant

**Electronic supplementary material**

The online version of this article (10.1186/s12870-017-1099-y) contains supplementary material, which is available to authorized users.

We sincerely thank Ms. Yuan Gao and Ms. Wenjuan Yang for laboratory assistance and Doc. Jun You for the language editing on the manuscript.

Funding {#FPar1}
=======

This work was funded by the National Natural Science Foundation of China (grant number: 31401412), the China Agriculture Research System (grant number: CARS-14), the Agricultural Science and Technology Innovation Project of the Chinese Academy of Agricultural Sciences (grant number: CAAS-ASTIP-2013-OCRI), and by using the Fundamental Research Fund from the Central Non-profit Scientific Institution (grant number: 1610172014003).

Availability of data and materials {#FPar2}
==================================

The datasets supporting the conclusions of this article are included within the article and its additional files; Sesame sequences in this article can be found from the sesame genome database ([http://ocri-genomics.org/Sinbase/](http://ocri-genomics.org/Sinbase)); The *Arabidopsis thaliana* gene sequences in this article were downloaded from UniProt (<http://www.uniprot.org>/). The raw RNA-seq reads and WRKY sequences are available at SesameFG (<http://ncgr.ac.cn/SesameFG>). All plant materials were selected from sesame germplasm provided by the Oil Crops Research Institute, Chinese Academy of Agricultural Sciences, Wuhan, China.

DL, XW, and ZX conceived and designed the experiments. DL and PL performed the experiments. DL, PL, XW, KD, JY and RZ analyzed the data. LW and YZ provided transcriptome data. DL and XW wrote the manuscript. All authors have read and approved the final manuscript.

Ethics approval and consent to participate {#FPar3}
==========================================

Not applicable.

Consent for publication {#FPar4}
=======================

Not applicable.

Competing interests {#FPar5}
===================

The authors declare that they have no competing interest.

Publisher's Note {#FPar6}
================

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
