Skip to main content

Transcriptome-wide identification and expression profiles of the WRKY transcription factor family in Broomcorn millet (Panicum miliaceum L.)



WRKY genes, as the most pivotal transcription factors in plants, play the indispensable roles in regulating various physiological processes, including plant growth and development as well as in response to stresses. Broomcorn millet is one of the most important crops in drought areas worldwide. However, the WRKY gene family in broomcorn millet remains unknown.


A total of 32 PmWRKY genes were identified in this study using computational prediction method. Structural analysis found that PmWRKY proteins contained a highly conserved motif WRKYGQK and two common variant motifs, namely WRKYGKK and WRKYGEK. Phylogenetic analysis of PmWRKYs together with the homologous genes from the representative species could classify them into three groups, with the number of 1, 15, and 16, respectively. Finally, the transcriptional profiles of these 32 PmWRKY genes in various tissues or under different abiotic stresses were systematically investigated using qRT-PCR analysis. Results showed that the expression level of 22 PmWRKY genes varied significantly under one or more abiotic stress treatments, which could be defined as abiotic stress-responsive genes.


This was the first study to identify the organization and transcriptional profiles of PmWRKY genes, which not only facilitates the functional analysis of the PmWRKY genes, and also lays the foundation to reveal the molecular mechanism of stress tolerance in this important crop.


Broomcorn millet (Panicum miliaceum L.) is one of the world’s oldest cultivated crops and also is one of the most important crops in drought areas. Before the domestication of rice and wheat, broomcorn millet was the main food staple in many semiarid regions of Asia, including China, Korea, India, Russia, Japan, and even in the entire Eurasian continent [13]. It’s supposed that broomcorn millet had been cultivated in East Asia since 10,000 years ago and played an essential impact on human civilisation [4]. Moreover, broomcorn millet has many favored agronomic traits, such as short growing season, high productivity, little nutrient, and water requirements as well as excellent tolerance to salt, drought, high temperature, and other extreme conditions.

With the advent of global climate changes, abiotic stresses such as salt, drought, freezing, and heat have become the hazardous threats to world’s agricultural production. Growth and development of crops were suppressed under abiotic stresses and then result in huge yield loss of production. It has been discovered that a large number of genes were induced when plants suffered to abiotic stresses, which could be were considered as stress-responsive genes. In general, the stress-responsive genes could be divided into functional proteins and regulatory factors [5]. Transcription factors, which are the dominant groups of regulatory factors to control the expression of other target genes, play the core roles in gene expression regulation and signal transduction.

WRKY transcription factors contain a highly conserved 60 amino acid DNA-binding regions, which designated as WRKY domain consisting of the conserved peptide sequence WRKYGQK (N-terminus) and a zinc-finger motif (C-terminus) [6]. According to the number of WRKY domains and structure features of zinc-finger motif, WRKY gene family can be classified into three groups. For the number of WRKY domain, group I of WRKY members contain two WRKY domains, while group II and group III have only one WRKY domain, respectively. For the structure features of zinc-finger motif, group I and group II have the C2H2 (C–X4–C–X22–23–HXH) zinc finger motif, while group III has the C2CH (C-X7-C-X23-HXC) motif. Furthermore, the Group II can be divided into five subgroups, including from IIa to IIe [7].

WRKY gene family has been revealed to be one of the most important transcription factor families in plant, which played the important roles in regulating various development and physiological processes [8]. In Arabidopsis, it was found that AtWRKY75 has a negative effect on root hair development [9]. AtWRKY18 and AtWRKY53 have a positive effect on leaf senescence. In addition, over-expression of AtWRKY18 and AtWRKY53 could results in delaying senescence [10], and AtWRKY34 is required for male gametogenesis [11]. It has also proven that WRKY proteins were involved in regulating pathogen infection and various abiotic stresses defenses [1214]. Previous report has found that over-expression of GhWRKY44 not only enhanced resistance to fungal pathogen R. solani in cotton, but also improved tolerance to bacterial pathogen R. solanacearum [15]. In terms of GhWRKY39-1, its expression was induced by pathogen infection and salt stresses [16]. Additionally, TaWRKY10 gene was found to be played a crucial role in wheat under salinity, cold, and drought stress responses [17]. Furthermore, co-expression of WRKYs had capacity to form the network to resist the stress. For instance, over-expression of OsWRKY11 and OsWRKY45 genes improved drought tolerance in rice [18]. Over-expression of GmWRKY54 and GmWRKY13 improved salt tolerance in Arabidopsis [19]. At present, extensive studies have been performed to identify and dissect the function of WRKY gene family in different plant species, included Arabidopsis [20], rice [21], wheat [22], grape [23], maize [24], and Brachypodium distachyon [25]. However, the organization and function of WRKY genes in broomcorn millet is completely unknown.

Given the above arguments, this study was performed with the following four objectives: (1) to systematically explore the WRKY genes in Broomcorn millet using the available EST and unigene sequences; (2) to investigate the gene composition, protein structure and orthologs construction of these identified WRKY genes; (3) to elucidate the evolutionary relationship and classification among the WRKY genes in broomcorn millet; (4) to comprehensively investigate expression profiles of these PmWRKY genes in various tissues or abiotic stress response to identify appropriate candidates for further functional studies. Taken together, the study will provide the novel insight into protein structures, evolutionary relationships, and expression pattern of WRKYs in Broomcorn millet, which will facilitate for further investigation of the functions of PmWRKY genes.

Result and discussion

Identification of PmWRKY genes

As one of the most important transcription factors in plant, WRKY gene family play the pivotal role in regulating plant growth and development as well as stress response [8]. Although the functions of several WRKY genes in Arabidopsis and other model crops have been systematically studied [1619], little is known about this gene family in the oldest food crops, broomcorn millet. At present, there is not any WRKY gene has been cloning or reported from broomcorn millet, which limited the study to reveal the molecular mechanism of stress resistance of this inimitable stress-tolerant crops.

To explore the organization and function of WRKY genes in broomcorn millet, the available ESTs and transcriptome contigs of this species were applied to identify WRKY genes through computational prediction. By blast search of the foxtail millet WRKY genes, a total of 59 unique ESTs/contigs showed high similarity (e-values range from 5.0 × 10−23 to 3.0 × 10−108). These contigs were considered as the putative WRKY genes in broomcorn millet. Then, the obtained sequences were submitted to the NCBI-CDD web server to analyze their conserved protein domain. Furthermore, the sequences contained the complete WRKY domain were further filtered to remove repeats by BLAST against foxtail millet transcription factor database ( Finally, a total of 32 unique PmWRKY genes includes complete conserved domain were identified in this study.

Phylogenetic analysis of PmWRKY genes

The WRKY gene family, which are high conserved in both monocots and eudicots, could be divided into three groups according to the protein structure and sequence similarity [21, 24]. To categorize and investigate the evolutionary relationships of PmWRKY genes, all putative 32 PmWRKY proteins in broomcorn millet, 75 WRKY proteins in rice, 50 WRKY proteins in maize, and 103 WRKY proteins in foxtail millet were selected to perform phylogenetic analysis. As shown in Fig. 1, phylogenetic results revealed that these PmWRKY proteins could be categorized into three groups, including Group I, Group II, and Group III. PmWRKY16 protein were considered as Group I, which included two WRKY domains and C2H2 (C–X4–C–X22–23–HXH) zinc-finger motif. A total of 15 PmWRKY proteins contained one WRKY domain and C2H2 (C-X4–5-C-X23-HXH) zinc-binding motif, which were classified to group II. The remaining 16 genes were assigned to group III, which had single WRKY domain and C2CH (C-X7-C-X23-HXC) zinc-binding motif. According to the WRKY subgroup classification of rice, maize, and foxtail millet, the group II of PmWRKYs was further classified into five subgroups, including group IIa (1), IIb (1) IIc (8), IId (1), and IIe (4).

Fig. 1
figure 1

Phylogenetic analysis of WRKY proteins among Broomcorn millet (32), foxtail millet (103), rice (75) and maize (50)

It was reported that variety size of WRKY transcription factor mainly depended on the variations of group III among WRKY genes, and group III of WRKY genes also played a pivotal role in evolutionary relationship [23]. Previous studies have reported that group III was the largest group of WRKY gene family in rice and wheat, which accounted for 38 and 41 %, respectively [26, 27]. while in Arabidopsis, the group II was the largest group, accounting for 24 % [7]. In this study, the group III is also found to be the largest group of WRKY gene family in broomcorn millet, accounting for 50 % of all PmWRKYs, which was consistent with the results of rice and wheat but different from Arabidopsis (Fig. 2).

Fig. 2
figure 2

The distribution of WRKY transcription factors from Broomcorn millet, Arabidopsis thaliana, Oryza sativa, Brachypodium distachyon, Triticum aestivum, Hordeum vulgare, and Zea mays

Protein structure analysis of PmWRKY

The protein structure of PmWRKY was further analyzed. The results showed that 23 PmWRKYs contained highly conserved sequence WRKYGQK. Among them, five PmWRKYs, including PmWRKY2, PmWRKY15, PmWRKY23, PmWRKY24, and PmWRKY28 proteins also have the most common variant sequence WRKYGKK. Adiditionally, four PmWRKYs, including PmWRKY5, PmWRKY6, PmWRKY8, and PmWRKY20 proteins contained the less common variant sequence WRKYGEK (Additional file 1: Figure S1). These results indicated that nine PmWRKY proteins had a single amino acid variation in their WRKY domain. It was reported that a single amino acid variation for WRKY domain had widely distributed in many species [28, 29]. Previous study demonstrated that five unusual VvWRKY domains possessed variations in grape, and seven unusual TaWRKY domains showed variations in wheat [23].

The WRKY domain is the most vital structural of WRKY proteins, and the usual WRKY domain is WRKYGQK motif, which can interact with the TTGACY core motif to activate downstream genes [7]. It has demonstrated that variance of the WRKYGQK motif could influence the normal activity of DNA binding [30]. When WRKY genes lack the WRKYGQK motif, the WRKYGKK motif may recognize binding sequences excluding the W-box element. For instance, tobacco WRKY protein NtWRKY12 and soybean WRKY protein GmWRKY6 cannot bind W-box element, while they can bind to WK-box (TTTTCCAC) [19, 31]. In addition, the variations were also found in zinc-finger motif of PmWRKYs, which was also reported in three VvWRKY proteins. However, the function of variations in zinc-finger motif remains unclear, which had only influence on the classification of the WRKY genes [23].

Conserved motifs analysis of PmWRKY

MEME program was used to predict the conserved motifs of PmWRKY proteins. A total of 20 motifs were identified to illustrate the WRKY protein structure in broomcorn millet (Fig. 3). Details of 20 conserved motifs were showed in Additional file 2: Figure S2. Results found that the number of motifs in PmWRKYs ranged from 1 to 7, and the length of motifs ranged from 7 to 50 amino acids. In addition, 3 motifs, namely motif 1, 2, and 3 were found in the WRKY domain. The others 17 motifs were found to be located outside the WRKY domain (Additional file 2: Figure S2). Meanwhile, motif 1, 2, and 3 were shared by 26 PmWRKYs, of which motif 1 was shared by PmWRKY5, PmWRKY9, PmWRKY11, PmWRKY26, and PmWRKY32, while motif 1 and motif 3 were shared by PmWRKY29. Furthermore, three group III members, viz PmWRKY18, PmWRKY31 and PmWRKY32 shared the motif 4. What is more, motif 8, constituted calmodulin-binding domain (CaMBD), were shared by group I (PmWRKY16), IIe (PmWRKY25), IId (PmWRKY26) and III (PmWRKY21) [32]. Additionally, it is generally recognized that WRKYs have a consensus coactivator domain, which was defined as LXXLL [33]. In the present study, it is found that motif 13 have LXXLL-like motif. Two PmWRKY proteins, including PmWRKY25 and PmWRKY26 were found to possess motif 8 (Fig. 3 and Additional file 2: Figure S2).

Fig. 3
figure 3

Unrooted phylogenetic tree (left) and conserved motifs (right) of 32 PmWRKY proteins. The phylogenetic tree was constructed by MEGA4.0. Different color represented various groups. MEME was used to predict motifs and these motifs represented with boxes

Homologous gene and functional annotation of PmWRKY proteins in broomcorn millet and rice

Orthologues and paralogues, two types of homologous genes, diverged through speciation in different species or duplication in single specie, respectively. Orthologous genes are widely distributed species and they are generally assumed to perform equivalent biological functions to share key properties of other species [34]. Many methods have been developed to identify the orthologous gene for studying the gene function [35]. Among them, the phylogenetic analysis is one of the most rapid, simple and relative accurate approach to evaluate the orthologs, which has been widely used in different organisms nowadays. As shown in Additional file 3: Table S1, several orthologous genes had been identified. Among these genes, one PmWRKY gene was found to be considered as the homologous genes of several OsWRKYs. For instance, PmWRKY5 is the homologous genes of OsWRKY18 and OsWRKY55. At the same time, several PmWRKY proteins were also found to have the same homologous gene in rice. For instance, PmWRKY4 and PmWRKY27 have the same homologous gene OsWRKY45.

It has been widely proved that the function of the given genes could be putatively predicted according to their identified homologous genes. Due to the function of several OsWRKYs had been well studied [18, 3648], the homologous WRKYs in rice could offer a reference for exploring the potential roles of PmWRKYs. Depending on the function of rice homologs, the function of some PmWRKYs was further putatively analyzed. As shown in Fig. 4, four groups were identified, of which seven PmWRKY genes were involved in abiotic stress responses. In addition, 14 PmWRKYs were found to involve in biotic stresses response, and the remaining 11 PmWRKYs were involved in other functions. Previous researches have reported that several WRKY genes had been well described in terms of their important roles in responses to abiotic and biotic stress [49, 50]. It was reported that about 51 % of AtWRKY genes in roots showed significantly high expression levels under salt stress [51], and 57 % of OsWRKY were differentially expressed in response to cold, drought, and salt stress [52] as well as, 70–90 % of VvWRKY genes have significantly differential expressions under various abiotic or biotic stresses [23]. Result showed that the proportion of the WRKYs involving in stresses response in broomcorn millet is similar with other plant.

Fig. 4
figure 4

Functional annotation of PmWRKY proteins based on homologous gene of rice. Unrooted phylogenetic tree were constructed for rice and Broomcorn millet proteins to identify homologous gene pairs

Tissue-specific expression patterns of PmWRKY

To explore the tissue specificity of PmWRKY genes, the expression levels of 32 PmWRKY genes were detected by RT-PCR in root, stem, leaf, and spike (Fig. 5 and Additional file 4: Figure S3). Almost one-third of PmWRKY genes represented clear tissue-specific expression. Among them, PmWRKY3 and PmWRKY17 had high expression levels in spike and leaf, respectively. The lowest expression of PmWRKY1 was found in root, while the strongest expressions of PmWRKY1 were detected in the stem, leaf, and spike. Furthermore, the expression levels of some PmWRKY genes were found to have no significant differences. Additionally, tissue-specific expression patterns of some PmWRKY genes were found to be consistents with their orthologous genes in rice (Additional file 4: Figure S3). PmWRKY2 represented relatively strongest expression in leaf tissue, and its orthologous gene OsWRKY34 had highly expression level in the same tissue. Nevertheless, some PmWRKY genes and its orthologous gene appeared to differential expression profiles. It is found that PmWRKY31 gene had a leaf-specific expression, while its orthologous gene OsWRKY47 had highly expression levels in root. Totally, more than half of PmWRKY genes showed relative consistent expression profiles with their orthologous genes, suggesting that the orthologous genes have similar biological functions with each other (Fig. 5).

Fig. 5
figure 5

Relative expression level of 32PmWRKY genes (left) and 29 homologous gene OsWRKY(right) in specific tissues. The result of qRT-PCR was analyzed by TanonGis Tools. Red represents decreased expression level and green represents increased expression level

Expression analysis of PmWRKY genes under abiotic stress

In order to verify the function of PmWRKY genes in broomcorn millet and provide the effective informations for future functional analysis, qRT-PCR analysis was used to reveal the differential expressions of each PmWRKYs under different abiotic stresses. The result indicated that 10 and 16 PmWRKY genes were significant induced by drought stress and cold stress treatment, respectively. Additionally, 11 genes were found to be induced by salt and heat stresses. Finally, a total of 22 PmWRKY genes had significantly differential expression under one or more stresses, which could be considered as stress-responsive genes. Interestingly, a diverse range of stress-specific expression profiles of these stress-responsive genes were found. The expression levels of PmWRKY32 and PmWRKY11 under stresses conditions were 26-fold and 108-fold higher than that of control condition, respectively. It is incredible that the expression level of PmWRKY7 was rapidly increased to 216-fold, 1500-fold under drought and cold treatment in comparsion with control group, respectively (Fig. 6). WRKY proteins domain have W box binding site (a four-stranded β-sheet with a zinc-binding pocket) to bind W box [53], which is the key regulatory element involving in stress response. It has reported that 11 OsWRKY genes and 8 TaWRKY genes were induced by salt, drought, cold, and heat stresses [22, 37]. OsWRKY69 was found to specifically bind ABL1 to regulates the rice stress responses [54]. BhWRKY1 were reported to bind W boxes of BhGolS1 under dehydration and ABA stresses [55]. Moreover, some WRKY proteins could interact with promoter elements, such as WK-box [56]. Therefore, these genes may be the invaluable factors in the stress metabolic network of broomcorn millet, and played the important role in abiotic stress response.

Fig. 6
figure 6

qRT-PCR was used to analyze the expression profiles of 32 PmWRKY genes under salt, drought, cold and heat stress. CK is untreated seedling. Error bars were obtained from three technical replicates. Asterisks reveal the gene significantly up-regulated or down-regulated under abiotic stresses by t-test (*P <0.05, **P <0.01)

Furthermore, it was found that orthologous genes have difference in expression patterns under stress treatment. PmWRKY25 was found to be up-regulated under salt, drought, cold and heat treatments, while its orthologous gene OsWRKY17 was down-regulated under these stresses treatments [57]. Furthermore, the expression levels of two PmWRKY genes, including PmWRKY10 and PmWRKY32 were the most significant different in all abiotic stresses. Meanwhile, four PmWRKY genes inculding PmWRKY3, PmWRKY5, PmWRKY23, and PmWRKY29 were rapidly and significantly expressed under one given stress, but these genes have no significant variance in other stresses. The expression level of PmWRKY23 have significantly up-regulated under heat stress, while it was down-regulated under salt or cold stresses and had no significant variance under drought stress. These results indicated the different PmWRKYs played the different role in regulating stress response, and further investigation of the functions for these PmWRKY genes is necessary. The results reported here provided some insight into the role of PmWRKYs playing in stress response and also provided the candidates for future functional studies.


In this study, 32 PmWRKY genes were identified, which was the first study to investigate the organization and abundance of WRKY in broomcorn millet. The protein structure, conserved motif composition, gene classification and phylogenetic relationship of these PmWRKYs were systematically analyzed and compared. Furthermore, in comparision with rice, orthologous genes analysis was performed to putatively predict the biological function of PmWRKYs. Finally, qRT-PCR analysis was used to investigate the transcriptional profiles in various tissues and under abiotic stresses. Based on the qPCR results, the tissues-specific and stress-responsive PmWRKYs were obtained. This study provide the helpful informations for further investigation of the function of PmWRKYs, and also shed light on the roles of WRKY genes playing in regulating development and physiological processes in this important stress-tolerance crops.


WRKY identification

A total of 211 original ESTs of broomcorn millet were obtained from the NCBI dbEST database ( The plant tissues of two Broomcorn millet cultivars (Yumi No.2 and Yumi No.3), including leaves, stems, roots, shoots, and spikes of different grown stages, were separately selected from three different plants for RNA extraction. The equal quantity RNA of each sample from the same variety was pooled and used for RNA-seq. Then, the public available ESTs together with the unigenes obtained from RNA-seq were used to predict WRKY genes. The WRKY cDNA sequences of Setaria italic were obtained from the Foxtail millet Transcription Factor Database (, and then these sequences were used to search the EST database through BLAST search with an e-value cut-off <10−5. The best hits were extracted as the EST representing PmWRKY. After manual curation, these ESTs were clustered and assembled by CAP3 tool and then obtained contigs and singletons were checked using NCBI-CDD search ( to identify the conserved protein domain. The sequences contained the complete WRKY domain was further used to BLAST against the FmTFDb database to remove redundant. Finally, a total of 32 putative PmWRKY genes were obtained and termed as PmWRKY1 to PmWRKY32.

Multiple sequence alignment, conserved motifs identification and Phylogenetic analysis

The Arabidopsis thaliana WRKY genes were obtained from DATF, while those of rice and maize were obtained from previous studies [24, 58]. To investigate the evolutionary relationships among these WRKYs, the predicted amino acid sequences of the WRKY genes of broomcorn millet and other species were aligned using ClutsalX1.83 program and then sequence identity and similarity were calculated. The phylogenetic tree was constructed using MEGA6.0 and neighbor-joining (NJ) method were adopted by 1,000 bootstrap replications [7]. Motifs of PmWRKYs were determined by using the MEME program ( and a schematic diagram of amino acid motifs of each PmWRKY gene was drawn accordingly.

Plant materials and abiotic stress treatments

Broomcorn millet cultivar Yumi No.3 was kindly provided by Prof. Baili Feng in this study. Seeds were planted in flowerpots and filled with a mixture of soil and sand (at a ratio of 1:1, v/v), and then it were adequately watered and grown in growth room under the condition of 22 °C, 16 h photoperiod (10,000 lux) and 20 °C, 8 h dark period. Leafs and roots were harvested at the given time after sowing: seedling stage (the 2th week), jointing stage (the 4th week), booting stage (the 6th week), filling stage (the 8th week) and mature stage (the 14th week). Stems were harvested at jointing stage (the 4th week), booting stage (the 6th week), filling stage (the 8th week) and mature stage (the 14th week). Spikes were harvested at mature stage (the 14th week). Meantime, in the fourth weeks, old whole seedlings of Yumi No.3 were treated by 200 mM NaCl, 19.2 % PEG, 4 °C or 38 °C for 24 h, which represented salt, drought, cold and heat treatment respectively. Then, these plants materials were collected and immediately frozen in liquid nitrogen until RNA extraction.

RT-PCR and qRT-PCR analysis

Total RNA from collected samples were isolated and performed with three independent experiments using the plant RNA isolation kit (Omega BioTek, USA). Concentration of the RNA was measured, and then equivalent mixed of three independent RNA for each sample to ensure that each reverse transcriptase reactions contained an equal amount of RNA from three independent seedlings. 1 μg high-quality RNA was used to synthesize the complementary DNA by the cDNA amplification kit (Vazyme, China) for downstream use. Primer Premier 5.0 was used to designed the primers for each PmWRKY gene and actin as reference gene [59] (Additional file 5: Table S2). The thermal profile for RT-PCR was as follows: 30 s at 94 °C, 35 cycles of 5 s at 94 °C, and 30 s at 58–60 °C, last 72 °C for 1 min. PCR products were checked by 1.0 % (w/v) agarose gel electrophoresis. The expression levels of specific tissues were analyzed from RT-PCR data with using the tanon Gis tool. qRT-PCR were performed using ABI 7300 real-time PCR system (Applied Biosystems, Grand Island, NY) based on SYBR Green II method. The 20 μl qRT-PCR mixtures contained 10 μl SYBR Premix buffer, 0.8 μl each of the primers (10 μM), 0.4 μl ROX Reference Dye (10 μmol/L), 2 μl cDNA and 6.8 μl PCRgrade water. The thermal profile for qRT-PCR was as follows: 30 s at 94 °C, 40 cycles of 5 s at 94 °C, and 30 s at 60 °C, last 72 °C for 1 min. Each reaction was run in triplicate to obtain the average value and Ct method was applied for the analysis gene expression profiles under abiotic stresses, and no template reactions were used as negative controls.


  1. Bellwood P, Gamble C, Le Blanc SA, Pluciennik M, Richards M, Terrell JE. First Farmers: the Origins of Agricultural Societies, by Peter Bellwood. Malden (MA): Blackwell, 2005; ISBN 0-631-20565-9 hardback£ 60; ISBN 0-631-20566-7 paperback£ 17.99, xix+ 360 pp., 59 figs., 3 tables. Camb Archaeol J. 2007;17(01):87–109.

    Article  Google Scholar 

  2. Fuller DQ. Agricultural origins and frontiers in South Asia: a working synthesis. J World Prehist. 2006;20(1):1–86.

    Article  Google Scholar 

  3. Crawford GW. East Asian Plant Domestication. In: Stark MT, editors. Archaeology of Asia. Oxford: Blackwell Publishing Ltd; 2006. p. 77-95. 

  4. Lu H, Zhang J, Liu K-b, Wu N, Li Y, Zhou K, Ye M, Zhang T, Zhang H, Yang X. Earliest domestication of common millet (Panicum miliaceum) in East Asia extended to 10,000 years ago. Proc Natl Acad Sci. 2009;106(18):7367–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Yamaguchi-Shinozaki K, Shinozaki K. Improving plant drought, salt and freezing tolerance by gene transfer of a single stress-inducible transcription factor. Rice Biotechnology: Improving Yield, Stress Tolerance and Grain Quality. 2001;236:176–89.

    CAS  Google Scholar 

  6. Rushton PJ, Somssich IE, Ringler P, Shen QJ. WRKY transcription factors. Trends Plant Sci. 2010;15(5):247–58.

    Article  CAS  PubMed  Google Scholar 

  7. Eulgem T, Rushton PJ, Robatzek S, Somssich IE. The WRKY superfamily of plant transcription factors. Trends Plant Sci. 2000;5(5):199–206.

    Article  CAS  PubMed  Google Scholar 

  8. Wu Z-J, Li X-H, Liu Z-W, Li H, Wang Y-X, Zhuang J. Transcriptome-wide identification of Camellia sinensis WRKY transcription factors in response to temperature stress. Mol Gen Genomics. 2016;291(1):255–69.

    Article  CAS  Google Scholar 

  9. Rishmawi L, Pesch M, Juengst C, Schauss AC, Schrader A, Hülskamp M. Non-cell-autonomous regulation of root hair patterning genes by WRKY75 in Arabidopsis. Plant Physiol. 2014;165(1):186–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Potschin M, Schlienger S, Bieker S, Zentgraf U. Senescence networking: WRKY18 is an upstream regulator, a downstream target gene, and a protein interaction partner of WRKY53. J Plant Growth Regul. 2014;33(1):106–18.

    Article  CAS  Google Scholar 

  11. Guan Y, Meng X, Khanna R, LaMontagne E, Liu Y, Zhang S. Phosphorylation of a WRKY transcription factor by MAPKs is required for pollen development and function in Arabidopsis. 2014.

    Google Scholar 

  12. Zhang L, Gu L, Ringler P, Smith S, Rushton PJ, Shen QJ. Three WRKY transcription factors additively repress abscisic acid and gibberellin signaling in aleurone cells. Plant Sci. 2015;236:214–22.

    Article  CAS  PubMed  Google Scholar 

  13. Hrmova M, Lopato S. Enhancing abiotic stress tolerance in plants by modulating properties of stress responsive transcription factors. In: Roberto Tuberosa R, Graner A, Frison E, editors. Genomics of Plant Genetic Resources. New York: Springer Publishing Company; 2014. p. 291–316.

  14. Banerjee A, Roychoudhury A. WRKY proteins: signaling and regulation of expression during abiotic stress responses. Sci World J. 2015;2015:807560.

    Article  Google Scholar 

  15. Li J, Wang J, Wang N, Guo X, Gao Z. GhWRKY44, a WRKY transcription factor of cotton, mediates defense responses to pathogen infection in transgenic Nicotiana benthamiana. Plant Cell Tissue Organ Cult. 2015;121(1):127–40.

    Article  CAS  Google Scholar 

  16. Shi W, Hao L, Li J, Liu D, Guo X, Li H. The Gossypium hirsutum WRKY gene GhWRKY39-1 promotes pathogen infection defense responses and mediates salt stress tolerance in transgenic Nicotiana benthamiana. Plant Cell Rep. 2014;33(3):483–98.

    Article  CAS  PubMed  Google Scholar 

  17. Bold O, Jeevan RJ, Lim YP, Vanjildorj E. Differatation of wheat WRKY transcription factor TAWRKY10 gene expression in abiotic stress resistance. Mongolian J Agric Sci. 2015;13(2):136–40.

    Article  Google Scholar 

  18. Qiu Y, Yu D. Over-expression of the stress-induced OsWRKY45 enhances disease resistance and drought tolerance in Arabidopsis. Environ Exp Bot. 2009;65(1):35–47.

    Article  CAS  Google Scholar 

  19. Zhou QY, Tian AG, Zou HF, Xie ZM, Lei G, Huang J, Wang CM, Wang HW, Zhang JS, Chen SY. Soybean WRKY‐type transcription factor genes, GmWRKY13, GmWRKY21, and GmWRKY54, confer differential tolerance to abiotic stresses in transgenic Arabidopsis plants. Plant Biotechnol J. 2008;6(5):486–503.

    Article  CAS  PubMed  Google Scholar 

  20. Eulgem T, Somssich IE. Networks of WRKY transcription factors in defense signaling. Curr Opin Plant Biol. 2007;10(4):366–71.

    Article  CAS  PubMed  Google Scholar 

  21. Wu K-L, Guo Z-J, Wang H-H, Li J. The WRKY family of transcription factors in rice and Arabidopsis and their origins. DNA Res. 2005;12(1):9–26.

    Article  CAS  PubMed  Google Scholar 

  22. Wu H, Ni Z, Yao Y, Guo G, Sun Q. Cloning and expression profiles of 15 genes encoding WRKY transcription factor in wheat (Triticum aestivem L.). Prog Nat Sci. 2008;18(6):697–705.

    Article  CAS  Google Scholar 

  23. Guo C, Guo R, Xu X, Gao M, Li X, Song J, Zheng Y, Wang X. Evolution and expression analysis of the grape (Vitis vinifera L.) WRKY gene family. J Exp Bot. 2014;65(6):1513–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Wei K-F, Chen J, Chen Y-F, Wu L-J, Xie D-X. Molecular phylogenetic and expression analysis of the complete WRKY transcription factor family in maize. DNA Res. 2012;19(2):153–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Wen F, Zhu H, Li P, Jiang M, Mao W, Ong C, Chu Z. Genome-wide evolutionary characterization and expression analyses of WRKY family genes in Brachypodium distachyon. DNA Res. 2014;21(3):327–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Xie Z, Zhang Z-L, Zou X, Huang J, Ruas P, Thompson D, Shen QJ. Annotations and functional analyses of the rice WRKY gene superfamily reveal positive and negative regulators of abscisic acid signaling in aleurone cells. Plant Physiol. 2005;137(1):176–89.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Zhu X, Liu S, Meng C, Qin L, Kong L, Xia G. WRKY transcription factors in wheat and their induction by biotic and abiotic stress. Plant Mol Biol Report. 2013;31(5):1053–67.

    Article  CAS  Google Scholar 

  28. Satapathy L, Singh D, Ranjan P, Kumar D, Kumar M, Prabhu KV, Mukhopadhyay K. Transcriptome-wide analysis of WRKY transcription factors in wheat and their leaf rust responsive expression profiling. Mol Gen Genomics. 2014;289(6):1289–306.

    Article  CAS  Google Scholar 

  29. Ding M, Chen J, Jiang Y, Lin L, Cao Y, Wang M, Zhang Y, Rong J, Ye W. Genome-wide investigation and transcriptome analysis of the WRKY gene family in Gossypium. Mol Gen Genomics. 2015;290(1):151–71.

    Article  CAS  Google Scholar 

  30. Maeo K, Hayashi S, Kojima-Suzuki H, Morikami A, Nakamura K. Role of conserved residues of the WRKY domain in the DNA-binding of tobacco WRKY family proteins. Biosci Biotechnol Biochem. 2001;65(11):2428–36.

    Article  CAS  PubMed  Google Scholar 

  31. van Verk MC, Pappaioannou D, Neeleman L, Bol JF, Linthorst HJ. A novel WRKY transcription factor is required for induction of PR-1a gene expression by salicylic acid and bacterial elicitors. Plant Physiol. 2008;146(4):1983–95.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Park CY, Lee JH, Yoo JH, Moon BC, Choi MS, Kang YH, Lee SM, Kim HS, Kang KY, Chung WS. WRKY group IId transcription factors interact with calmodulin. FEBS Lett. 2005;579(6):1545–50.

    Article  CAS  PubMed  Google Scholar 

  33. Savkur R, Burris T. The coactivator LXXLL nuclear receptor recognition motif. J Pept Res. 2004;63(3):207–12.

    Article  CAS  PubMed  Google Scholar 

  34. Cai H, Tian S, Dong H. Large scale in silico identification of MYB family genes from wheat expressed sequence tags. Mol Biotechnol. 2012;52(2):184–92.

    Article  CAS  PubMed  Google Scholar 

  35. Remm M, Storm CE, Sonnhammer EL. Automatic clustering of orthologs and in-paralogs from pairwise species comparisons. J Mol Biol. 2001;314(5):1041–52.

    Article  CAS  PubMed  Google Scholar 

  36. Swarbrick P, Huang K, Liu G, Slate J, Press M, Scholes J. Global patterns of gene expression in rice cultivars undergoing a susceptible or resistant interaction with the parasitic plant Striga hermonthica. New Phytol. 2008;179(2):515–29.

    Article  CAS  PubMed  Google Scholar 

  37. Qiu Y, Jing S, Fu J, Li L, Yu D. Cloning and analysis of expression profile of 13WRKY genes in rice. Chin Sci Bull. 2004;49(20):2159–68.

    CAS  Google Scholar 

  38. Gloria M. Analysis of the interaction transcriptome during biotrophicinvasion by the blast fungus. In: Magnaporthe oryzae [PhD dissertation]. USA: Kansas State University Kansas; 2007. p. 88–90.

    Google Scholar 

  39. Lan A, Huang J, Zhao W, Peng Y, Chen Z, Kang D. A salicylic acid‐induced rice (Oryza sativa L.) transcription factor OsWRKY77 is involved in disease resistance of Arabidopsis thaliana. Plant Biol. 2013;15(3):452–61.

    Article  CAS  PubMed  Google Scholar 

  40. Wei K, Chen J, Chen Y, Wu L, Xie D. Multiple-strategy analyses of ZmWRKY subgroups and functional exploration of ZmWRKY genes in pathogen responses. Mol BioSyst. 2012;8(7):1940–9.

    Article  CAS  PubMed  Google Scholar 

  41. Qiu D, Xiao J, Xie W, Cheng H, Li X, Wang S. Exploring transcriptional signalling mediated by OsWRKY13, a potential regulator of multiple physiological processes in rice. BMC Plant Biol. 2009;9(1):74.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Zhang J, Peng Y, Guo Z. Constitutive expression of pathogen-inducible OsWRKY31 enhances disease resistance and affects root growth and auxin response in transgenic rice plants. Cell Res. 2008;18(4):508–21.

    Article  CAS  PubMed  Google Scholar 

  43. Wang H-J, Wan A-R, Hsu C-M, Lee K-W, Yu S-M, Jauh G-Y. Transcriptomic adaptations in rice suspension cells under sucrose starvation. Plant Mol Biol. 2007;63(4):441–63.

    Article  CAS  PubMed  Google Scholar 

  44. Raineri J, Wang S, Peleg Z, Blumwald E, Chan RL. The rice transcription factor OsWRKY47 is a positive regulator of the response to water deficit stress. Plant Mol Biol. 2015;88(4-5):401–13.

    Article  CAS  PubMed  Google Scholar 

  45. Liu X, Li Y, Wang L, Liu X, Wang C, Wang L, Pan Q. The effect of the rice blast resistance gene Pi36 on the expression of disease resistance-related genes. Chin Sci Bull. 2010;55(18):1881–8.

    Article  CAS  Google Scholar 

  46. Xu Q, Feng W, Peng H, Ni Z, Sun Q. TaWRKY71, a WRKY transcription factor from wheat, enhances tolerance to abiotic stress in transgenic Arabidopsis thaliana. Cereal Res Commun. 2013;42(1):47–57.

    Article  CAS  Google Scholar 

  47. Shimono M, Sugano S, Nakayama A, Jiang C-J, Ono K, Toki S, Takatsuji H. Rice WRKY45 plays a crucial role in benzothiadiazole-inducible blast resistance. Plant Cell. 2007;19(6):2064–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Cui X-W, Zhang Y, Qi F-Y, Gao J, Chen Y-W, Zhang C-L. Overexpression of a moso bamboo (Phyllostachys edulis) transcription factor gene PheWRKY1 enhances disease resistance in transgenic Arabidopsis thaliana. Botany. 2013;91(7):486–94.

    Article  CAS  Google Scholar 

  49. Jiang Y, Duan Y, Yin J, Ye S, Zhu J, Zhang F, Lu W, Fan D, Luo K. Genome-wide identification and characterization of the Populus WRKY transcription factor family and analysis of their expression in response to biotic and abiotic stresses. J Exp Bot. 2014;65(22):6629–44.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Fan X, Guo Q, Xu P, Gong Y, Shu H, Yang Y, Ni W, Zhang X, Shen X. Transcriptome-Wide Identification of Salt-Responsive Members of the WRKY Gene Family in Gossypium aridum. 2015.

    Google Scholar 

  51. Jiang Y, Deyholos MK. Comprehensive transcriptional profiling of NaCl-stressed Arabidopsis roots reveals novel classes of responsive genes. BMC Plant Biol. 2006;6(1):25.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Ramamoorthy R, Jiang S-Y, Kumar N, Venkatesh PN, Ramachandran S. A comprehensive transcriptional profiling of the WRKY gene family in rice under various abiotic and phytohormone treatments. Plant Cell Physiol. 2008;49(6):865–79.

    Article  CAS  PubMed  Google Scholar 

  53. Tripathi P, Rabara RC, Rushton PJ. A systems biology perspective on the role of WRKY transcription factors in drought responses in plants. Planta. 2014;239(2):255–66.

    Article  CAS  PubMed  Google Scholar 

  54. Yang X, Yang Y-N, Xue L-J, Zou M-J, Liu J-Y, Chen F, Xue H-W. Rice ABI5-Like1 regulates abscisic acid and auxin responses by affecting the expression of ABRE-containing genes. Plant Physiol. 2011;156(3):1397–409.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Wang Z, Zhu Y, Wang L, Liu X, Liu Y, Phillips J, Deng X. A WRKY transcription factor participates in dehydration tolerance in Boea hygrometrica by binding to the W-box elements of the galactinol synthase (BhGolS1) promoter. Planta. 2009;230(6):1155–66.

    Article  CAS  PubMed  Google Scholar 

  56. Niu CF, Wei W, Zhou QY, Tian AG, Hao YJ, Zhang WK, Ma B, Lin Q, Zhang ZB, Zhang JS. Wheat WRKY genes TaWRKY2 and TaWRKY19 regulate abiotic stress tolerance in transgenic Arabidopsis plants. Plant Cell Environ. 2012;35(6):1156–70.

    Article  CAS  PubMed  Google Scholar 

  57. Tyagi A, Kapoor S, Khurana J, Ray S. Expression data for stress treatment in rice seedlings. NCBI; Accessed 10 Oct 2015.

  58. Ross CA, Liu Y, Shen QJ. The WRKY gene family in rice (Oryza sativa). J Integr Plant Biol. 2007;49(6):827.

    Article  CAS  Google Scholar 

  59. Van den Berg N, Crampton BG, Hein I, Birch PR, Berger DK. High-throughput screening of suppression subtractive hybridization cDNA libraries using DNA microarray analysis. 2004.

    Google Scholar 

Download references


This project was financially supported by Open project of the State Key Laboratory of Crop Stress Biology in Arid Areas, Northwest A&F University (Grant no.CSBAA2014002). We sincerely appreciate Prof. Baili Feng for kindly providing broomcorn millet seeds.

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Weining Song or Xiaojun Nie.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

XN and WS design and interpretation of this experiment. XD collected experimental materials. HY, MW and SL carried out experimental work. HY data analysis and prepared this manuscript. All authors read and approved the final manuscript.

Additional files

Additional file 1: Figure S1.

Multiple-sequence alignment of the WRKY protein domain from PmWRKYs and OsWRKYs. Conserved amino acids were indicated by blank background, conserved WRKY domains and zinc-finger motifs were indicated by red box. (DOC 1365 kb)

Additional file 2: Figure S2.

Sequence logos of PmWRKY domain. The PmWRKY proteins domain submitted to MEME server. The total height of stack was used to shows ‘information content’ of that position in the motif. Height of letters in stack suggests probability of each amino acids at that position. (DOC 558 kb)

Additional file 3: Table S1.

The identified PmWRKY proteins in Broomcorn millet and their putative orthologous gene in rice. (DOC 47 kb)

Additional file 4: Figure S3.

Expression level of 32 PmWRKY genes in different tissues. Actin was an internal reference gene. (DOC 797 kb)

Additional file 5: Table S2.

Primers of 32 PmWRKY genes. (DOC 54 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yue, H., Wang, M., Liu, S. et al. Transcriptome-wide identification and expression profiles of the WRKY transcription factor family in Broomcorn millet (Panicum miliaceum L.). BMC Genomics 17, 343 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: