Abstract
(1) Rhododendron is one of the top ten traditional flowers in China, with both high ornamental and economic values. However, with the change of the environment, Rhododendron suffers from various biological stresses. The WRKY transcription factor is a member of the most crucial transcription factor families, which plays an essential regulatory role in a variety of physiological processes and developmental stresses. (2) In this study, 57 RsWRKYs were identified using genome data and found to be randomly distributed on 13 chromosomes. Based on gene structure and phylogenetic relationships, 57 proteins were divided into three groups: I, II, and III. Multiple alignments of RsWRKYs with Arabidopsis thaliana homologous genes revealed that WRKY domains in different groups had different conserved sites. RsWRKYs have a highly conserved domain, WRKYGQK, with three variants, WRKYGKK, WRKYGEK, and WRKYGRK. Furthermore, cis-acting elements analysis revealed that all of the RsWRKYs had stress and plant hormone cis-elements, with figures varying by group. Finally, the expression patterns of nine WRKY genes treated with gibberellin acid (GA), methyl jasmonate (MeJA), heat, and drought in Rhododendron were also measured using quantitative real-time PCR (qRT-PCR). The results showed that the expression levels of the majority of RsWRKY genes changed in response to multiple phytohormones and abiotic stressors. (3) This current study establishes a theoretical basis for future studies on the response of RsWRKY transcription factors to various hormone and abiotic stresses as well as a significant foundation for the breeding of new stress-tolerant Rhododendron varieties.