Abstract
The multigene family WRKY is plant-specific transcription factors, which is involved in numerous aspects of plant growth and development. Here, we report a genome-wide search for the rubber tree WRKY (HbWRKY) genes and their expression profiles at various tissues including latexes. We also present a comprehensive overview of whole HbWRKY gene family, including gene structure, chromosomal locations, conserved protein domains, gene duplications and phylogenetic inferences. We identify a total of 108 full-length HbWRKY genes, which were classified into the three major subgroups. The genome-scale syntenic analysis showed that WGDs might have played a major role in the evolution of the HbWRKY gene family. Our further real-time reverse transcription-PCR (qRT-PCR) experiments validated that the HbWRKY genes are divergently expressed and gene expression divergence seemingly existed among the duplicated WRKY genes, suggesting a fundamentally functional divergence of the duplicated WRKY paralogs in rubber tree. We also investigated the gene expression profiling of the HbWRKYs under salt stress, to which the majority of genes differentially responded. This study acquires a new insight into the evolution of the HbWRKYs, which will help further functional characterization of theses candidate genes in rubber tree.
Similar content being viewed by others
References
Ahmad A, Niwa Y, Goto S, Ogawa T, Shimizu M, Suzuki A, Kobayashi K, Kobayashi H (2015) bHLH106 integrates functions of multiple genes through their G-box to confer salt tolerance on Arabidopsis. PLoS One 10(5):e0126872. https://doi.org/10.1371/journal.pone.0126872
Backhaus RA (1985) Rubber formation in plants-a mini-review. ISR J Plant Sci 34(2–4):283–293. https://doi.org/10.1080/0021213X.1985.10677027
Bailey TL, Boden M, Buske FA et al (2009) MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res 37(suppl_2):W202–W208. https://doi.org/10.1093/nar/gkp335
Bi C, Xu Y, Ye Q, Yin T, Ye N (2016) Genome-wide identification and characterization of WRKY gene family in Salix suchowensis. Peerj 4:e2437. https://doi.org/10.7717/peerj.2437
Bowers JE (1990) Natural rubber-producing plants for the United States. Proc Natl Acad Sci USA 100(3):1352–1357. https://doi.org/10.1073/pnas.242738899
Cannon SB, Mitra A, Baumgarten A, Young ND, May G (2004) The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol 4:10. https://doi.org/10.1186/1471-2229-4-10
Ciolkowski I, Wanke D, Birkenbihl RP, Somssich IE (2008) Studies on DNA-binding selectivity of WRKY transcription factors lend structural clues into WRKY-domain function. Plant Mol Biol 68(1–2):81–92
Eulgem T, Somssich IE (2007) Networks of WRKY transcription factors in defense signaling. Curr Opin Plant Biol 10(4):366–371. https://doi.org/10.1016/j.pbi.2007.04.020
Eulgem T, Rushton PJ, Schmelzer E, Hahlbrock K, Somssich IE (1999) Early nuclear events in plant defense: rapid gene activation by WRKY transcription factors. EMBO J 18(17):4689–4699
Eulgem T, Rushton PJ, Robatzek S, Somssich IE (2000) The WRKY superfamily of plant transcription factors. Trends Plant Sci 5(5):199–206. https://doi.org/10.1016/S1360-1385(00)01600-9
Guo C, Guo R, Xu X, Gao M, Li X, Song J, Zheng Y, Wang X (2014) Evolution and expression analysis of the grape (Vitis vinifera L.) WRKY gene family. J Exp Bot 65(6):1513–1528. https://doi.org/10.1093/jxb/eru007
Hellsten U, Khokha MK, Grammer TC, Harland RM, Richardson P, Rokhsar DS (2007) Accelerated gene evolution and subfunctionalization in the pseudotetraploid frog Xenopus laevis. BMC Biol 5:31. https://doi.org/10.1186/1741-7007-5-31
Huang X, Li K, Xu X, Yao Z, Jin C, Zhang S (2015) Genome-wide analysis of WRKY transcription factors in white pear (Pyrus bretschneideri) reveals evolution and patterns under drought stress. BMC Genomics 16:1104. https://doi.org/10.1186/s12864-015-2233-6
International Rice Genome Sequencing P (2005) The map-based sequence of the rice genome. Nature 436(7052):793–800. https://doi.org/10.1038/nature03895
Ishiguro S, Nakamura K (1994) Characterization of a cDNA encoding a novel DNA-binding protein, SPF1, that recognizes SP8 sequences in the the 5′upstream regions of genes coding for sporamin and β-amylase from sweet potato. Mol Gen Genet 244(6):563–571
Johnson CS (2002) TRANSPARENT TESTA GLABRA2, a trichome and seed coat development gene of Arabidopsis, encodes a WRKY transcription factor. Plant Cell 14(6):1359–1375. https://doi.org/10.1105/tpc.001404
Kim K-C, Lai Z, Fan B, Chen Z (2008) Arabidopsis WRKY38 and WRKY62 transcription factors interact with histone deacetylase 19 in basal defense. Plant Cell 20(9):2357–2371
Kumar S, Stecher G, Tamura K (2016) MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol 33(7):1870–1874. https://doi.org/10.1093/molbev/msw054
Lagacé M, Matton DP (2004) Characterization of a WRKY transcription factor expressed in late torpedo-stage embryos of Solanum chacoense. Planta 219(1):185–189
Larkin MA, Blackshields G, Brown N et al (2007) Clustal W and Clustal X version 2.0. Bioinformatics 23(21):2947–2948. https://doi.org/10.1093/bioinformatics/btm404
Li WH (1993) Unbiased estimation of the rates of synonymous and nonsynonymous substitution. J Mol Evol 36:96–99
Li WH, Yang J, Gu X (2005) Expression divergence between duplicate genes. Trends Genet 21(11):602–607. https://doi.org/10.1016/j.tig.2005.08.006
Li HL, Guo D, Yang ZP, Tang X, Peng SQ (2014) Genome-wide identification and characterization of WRKY gene family in Hevea brasiliensis. Genomics 104(1):14–23. https://doi.org/10.1016/j.ygeno.2014.04.004
Li P, Song A, Gao C, Wang L, Wang Y, Sun J, Jiang J, Chen F, Chen S (2015) Chrysanthemum WRKY gene CmWRKY17 negatively regulates salt stress tolerance in transgenic chrysanthemum and Arabidopsis plants. Plant Cell Rep 34(8):1365–1378. https://doi.org/10.1007/s00299-015-1793-x
Lieberei R (2007) South American leaf blight of the rubber tree (Hevea spp.): new steps in plant domestication using physiological features and molecular markers. Ann Bot-London 100(6):1125–1142. https://doi.org/10.1093/aob/mcm133
Liu JJ, Ekramoddoullah AK (2009) Identification and characterization of the WRKY transcription factor family in Pinus monticola. Genome 52(1):77–88. https://doi.org/10.1139/G08-106
Liu J et al (2020) The chromosome-based rubber tree genome provides new insights into spurge genome evolution and rubber biosynthesis. Mol Plant 13(2):336–350. https://doi.org/10.1016/j.molp.2019.10.017
Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2-△△CT method. Methods 25(4):402–408. https://doi.org/10.1006/meth.2001.1262
Makova KD, Li WH (2003) Divergence in the spatial pattern of gene expression between human duplicate genes. Genome Res 13(7):1638–1645. https://doi.org/10.1101/gr.1133803
Mishra P, Singh N, Jain A et al (2018) Identification of cis-regulatory elements associated with salinity and drought stress tolerance in rice from co-expressed gene interaction networks. Bioinformation 14(3):123–131. https://doi.org/10.6026/97320630014123
Mzid R, Marchive C, Blancard D, Deluc L, Barrieu F, Corio-Costet MF, Drira N, Hamdi S, Lauvergeat V (2007) Overexpression of VvWRKY2 in tobacco enhances broad resistance to necrotrophic fungal pathogens. Physiol Plant 131(3):434–447. https://doi.org/10.1111/j.1399-3054.2007.00975.x
Nan H, Gao LZ (2019) Genome-wide analysis of WRKY genes and their response to hormone and mechanic stresses in carrot. Front Genet 10:363. https://doi.org/10.3389/fgene.2019.00363
Nan H, Li W, Lin YL, Gao LZ (2020) Genome-wide analysis of WRKY genes and their response to salt stress in the wild progenitor of Asian cultivated rice, Oryza rufipogon. Front Genet 11:359
Narusaka Y, Nakashima K, Shinwari ZK, Sakuma Y, Furihata T, Abe H, Narusaka M, Shinozaki K, Yamaguchi-Shinozaki K (2003) Interaction between two cis-acting elements, ABRE and DRE, in ABA-dependent expression of Arabidopsis rd29A gene in response to dehydration and high-salinity stresses. Plant J 34:137–148
Oku E, Iwara A, Ekukinam E (2012) Effects of age of rubber (Hevea brasiliensis Muell Arg.) plantation on pH, organic carbon, organic matter, nitrogen and micronutrient status of Ultisols in the humid forest zone of Nigeria. Kasetsart Journal-Natural Science 46(5):684–693
Pandey SP, Somssich IE (2009) The role of WRKY transcription factors in plant immunity. Plant Physiol 150(4):1648–1655. https://doi.org/10.1104/pp.109.138990
Rahman AYA, Usharraj AO, Misra BB, Thottathil GP, Jayasekaran K, Feng Y et al (2013) Draft genome sequence of the rubber tree (Hevea brasiliensis). BMC Genomics:14–75
Rushton PJ, Somssich IE, Ringler P, Shen QJ (2010) WRKY transcription factors. Trends Plant Sci 15(5):247–258. https://doi.org/10.1016/j.tplants.2010.02.006
Song H, Wang P, Lin JY, Zhao C, Bi Y, Wang X (2016) Genome-wide identification and characterization of WRKY gene family in peanut. Front Plant Sci 7:534. https://doi.org/10.3389/fpls.2016.00534
Sun C (2003) A novel WRKY transcription factor, SUSIBA2, participates in sugar signaling in barley by binding to the sugar-responsive elements of the iso1 promoter. Plant Cell 15(9):2076–2092. https://doi.org/10.1105/tpc.014597
Tang J, Wang F, Hou XL, Wang Z, Huang ZN (2013) Genome-wide fractionation and identification of WRKY transcription factors in Chinese cabbage (Brassica rapa ssp. pekinensis) reveals collinearity and their expression patterns under abiotic and biotic stresses. Plant Mol Biol Rep 32(4):781–795. https://doi.org/10.1007/s11105-013-0672-2
Tang C et al (2016) The rubber tree genome reveals new insights into rubber production and species adaptation. Nat Plants 2(6):16073. https://doi.org/10.1038/nplants.2016.73
Tao Z, Kou Y, Liu H, Li X, Xiao J, Wang S (2011) OsWRKY45 alleles play different roles in abscisic acid signalling and salt stress tolerance but similar roles in drought and cold tolerance in rice. J Exp Bot 62(14):4863–4874. https://doi.org/10.1093/jxb/err144
Todaka D, Nakashima K, Shinozaki K, Yamaguchi-Shinozaki K (2012) Toward understanding transcriptional regulatory networks in abiotic stress responses and tolerance in rice. Rice (N Y) 5(1):6–6. https://doi.org/10.1186/1939-8433-5-6
van Verk MC, Pappaioannou D, Neeleman L, Bol JF, Linthorst HJM (2008) A novel WRKY transcription factor is required for induction of PR-1a gene expression by salicylic acid and bacterial elicitors. Plant Physiol 146(4):1983–1995. https://doi.org/10.1104/pp.107.112789
Wang Y, Guo D, Li H-L, Peng SQ (2013) Characterization of HbWRKY1, a WRKY transcription factor from Hevea brasiliensis that negatively regulates HbSRPP. Plant Physiol Bioch 71:283–289. https://doi.org/10.1016/j.plaphy.2013.07.020
Wang M, Vannozzi A, Wang G et al (2014) Genome and transcriptome analysis of the grapevine (Vitis vinifera L.) WRKY gene family. Hortic Res 1:14016. https://doi.org/10.1038/hortres.2014.16
Wu KL, Guo ZJ, Wang HH, Li J (2005) The WRKY family of transcription factors in rice and Arabidopsis and their origins. DNA Res 12(1):9–26. https://doi.org/10.1093/dnares/12.1.9
Wu Z, Xu X, Xiong W, Wu P, Chen Y, Li M, Wu G, Jiang H (2015) Genome-wide analysis of the NAC gene family in physic nut (Jatropha curcas L.). PLoS One 10(6):e0131890. https://doi.org/10.1371/journal.pone.0131890
Xiao Y, Zhou L, Lei X, Cao H, Wang Y, Dou Y, Tang W, Xia W (2017) Genome-wide identification of WRKY genes and their expression profiles under different abiotic stresses in Elaeis guineensis. PLoS One 12(12):e0189224. https://doi.org/10.1371/journal.pone.0189224
Xu H, Watanabe KA, Zhang L, Shen QJ (2016) WRKY transcription factor genes in wild rice Oryza nivara. DNA Res 23(4):311–323. https://doi.org/10.1093/dnares/dsw025
Yang Z (2007) PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol 24(8):1586–1591. https://doi.org/10.1093/molbev/msm088
Yang ZH, Bielawski JP (2000) Statistical methods for detecting molecular adaptation. Trends Ecol Evol 15:496–503
Yang Y, Zhou Y, Chi Y, Fan B, Chen Z (2017) Characterization of soybean WRKY gene family and identification of soybean WRKY genes that promote resistance to soybean cyst nematode. Sci Rep 7(1):17804. https://doi.org/10.1038/s41598-017-18235-8
Yao QY, Xia EH, Liu FH, Gao LZ (2015) Genome-wide identification and comparative expression analysis reveal a rapid expansion and functional divergence of duplicated genes in the WRKY gene family of cabbage, Brassica oleracea var capitata. Gene 557(1):35–42. https://doi.org/10.1016/j.gene.2014.12.005
Yu J et al (2005) The genomes of Oryza sativa: a history of duplications. PLoS Biol 3(2):e38. https://doi.org/10.1371/journal.pbio.0030038
Zhang J (2003) Evolution by gene duplication: an update. Trends Ecol Evol 18(6):292–298. https://doi.org/10.1016/s0169-5347(03)00033-8
Zhang Y, Wang L (2005) The WRKY transcription factor superfamily: its origin in eukaryotes and expansion in plants. BMC Evol Biol 5(1):1
Zhang CQ, Xu Y, Lu Y, Yu HX, Gu MH, Liu QQ (2011) The WRKY transcription factor OsWRKY78 regulates stem elongation and seed development in rice. Planta 234(3):541–554. https://doi.org/10.1007/s00425-011-1423-y
Zhou QY, Tian AG, Zou HF, Xie ZM, Lei G, Huang J, Wang CM, Wang HW, Zhang JS, Chen SY (2008) Soybean WRKY-type transcription factor genes, GmWRKY13, GmWRKY21, and GmWRKY54, confer differential tolerance to abiotic stresses in transgenic Arabidopsis plants. Plant Biotechnol J 6(5):486–503. https://doi.org/10.1111/j.1467-7652.2008.00336.x
Acknowledgments
We appreciate the anonymous reviewers for their comments on this manuscript. This study was supported by MOST (2018YFD1000502) (to H. Huang) and Yunnan Innovation Team Project (to L. Z. Gao).
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of Interest
The authors declare no conflicts of interest.
Additional information
Communicated by: Ray Ming
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
ESM 1
Multiple sequence alignment of domains among the HbWRKY proteins. The alignment was performed using Clustal W. Dashes indicate gaps. ‘N’ and ‘C’ denote the N-terminal and C-terminal WRKY domains of a specific WRKY protein. (PNG 5171 kb)
ESM 2
Phylogenetic analysis of the 324 WRKY proteins in Arabidopsis, rice and rubber tree. A. The phylogenetic tree was constructed using Neighbor-Joining method based on WRKY domains. B. The phylogenetic tree was constructed using maximum-likelihood method based on WRKY domains. Reliability of the predicted trees was tested using bootstrapping with 1000 replicates. (PNG 1845 kb)
ESM 3
The schematic distribution of motifs within the HbWRKY proteins. (PNG 771 kb)
ESM 4
Phylogenetic relationships and conserved motif compositions of the WRKY genes in the rubber tree genome. Phylogenetic tree was constructed using the maximum-likelihood method based on alignments of the predicted complete proteins. The conserved motifs were detected using MEME software and represented by colored boxes. Reliability of the predicted tree was tested using bootstrapping with 1000 replicates. The length of WRKY proteins can be estimated using the scale at the bottom and conserved motif were shown in Table S1. (PNG 626 kb)
ESM 5
The exon-intron structures of the HbWRKY genes. The exons and introns are indicated by blue rectangles and black lines, respectively. (PNG 388 kb)
ESM 6
(XLS 47 kb)
ESM 7
(XLS 39 kb)
ESM 8
(XLS 13 kb)
ESM 9
(XLS 25 kb)
ESM 10
(XLSX 5732 kb)
ESM 11
(XLS 30 kb)
Rights and permissions
About this article
Cite this article
Nan, H., Lin, YL., Liu, J. et al. Genome-Wide Analysis of the WRKY Transcription Factor Gene Family and Their Response to Salt Stress in Rubber Tree. Tropical Plant Biol. 14, 22–33 (2021). https://doi.org/10.1007/s12042-020-09268-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s12042-020-09268-x