1. Introduction
A key event in the life cycle of angiosperms is the formation of flowers, which marks the transitions of plants from vegetative growth to reproductive growth [
1]. It not only directly affects the growth, development, and yield of crop plants, but also has a great impact on horticultural plants as the flower is the most important ornamental part in landscaping [
2]. Flowering is controlled by a very large and fine-tuned regulatory network, which integrates various exogenous and endogenous factors to determine the most appropriate flowering time [
3,
4,
5].
On the other hand, plants are also challenged with multiple types of abiotic stress and must respond to survive. Among them, high salinity and drought are the main abiotic stresses that affect plant growth and development and even cause death [
6]. Under high salinity, ions from the soil accumulate in the cytoplasm and thus disrupt the osmotic regulation. It further induces secondary oxidative stress, leading to the inhibition of photosynthesis and metabolism and ultimately a deduction in plant yield [
7]. Similar to high salinity, drought causes osmotic stress and secondary stress including oxidative stress, damage to cellular components (e.g., membrane lipids, proteins, and nucleic acids), and metabolic dysfunction [
8]. To cope with these abiotic stresses, evolution has equipped plants with a series of complex regulatory networks for survival [
9].
The two seemingly unrelated processes of flowering and response to abiotic stress have one thing in common: both are under tight control of transcriptional factors (TFs). Indeed, the critical roles of TFs in flowering have been extensively studied [
10]. In fact, some of the most important proteins in flowering pathways are TFs or transcriptional regulators including Constans (CO) in the photoperiod pathway and flowering locus C (FLC), a MADS-box TF that represses flowering. Similarly, TFs play a critical role in the responses to abiotic stresses [
11]. By activating or inhibiting the expression of target genes, they regulate plant metabolism and adaptability to abiotic stress, leading to an increased interest in using TFs to combat environmental stresses in the face of global warming [
12,
13].
Interestingly, TFs may regulate multiple independent biological processes because they can interact with distinct sets of substrates. The WRKY family is a good example that controls a diverse range of physiological events via a complex gene network [
14,
15]. For instance, WRKY regulates plant tolerance to salt, drought, low temperature, high temperature, and other stresses as well as flowering by interacting with abscisic acid (ABA), jasmonic acid (JA), and ROS. WRKYs are named after the highly conserved WRKYGQK domain at the N-terminus [
16]. The C terminus contains a CX4-7CX22 or 23HXH/C zinc finger that is important for protein interactions and DNA binding [
17,
18]. Based on the number of WRKY domains and the type of zinc finger structure, WRKYs are grouped into group Ⅰ (two WRKY domains and a C2H2 type zinc finger structure), group II (one WRKY domain and a C2H2 type zinc finger structure), and group III (one WRKY domain and a C2HC type zinc finger structure) [
19]. Advances in sequencing have led to the genome-wide identification of WRKYs in many plant species including 74 members in
Arabidopsis, 59 in grape, 127 in apple, 102 in rice, and 119 in corn [
20].
Iris laevigata Fisch is native in Northeast China with strong cold resistance and beautiful flowers of unique shape and bright color[
21]. Such desirable traits make
I. laevigata a good candidate for both landscaping and industrial use as a cutting flower. However, the potential is dampened because of the short flowering period with only 2 days for a single flower and 20 days for the population[
21].
I. laevigata is also sensitive to abiotic stress[
22]. Thus, extending the flowering period and increasing the resistance to abiotic stress are two central goals in the field of breeding
I. laevigata of novel properties[
23,
24]. One prerequisite to achieve these goals is to identify key genes controlling the target traits, which has been propelled tremendously by the advances in genome sequencing [
25,
26,
27].
However, mechanistic studies at the molecular levels in
Iris lag far behind compared with other plant species due to the lack of genome information. Our group pioneered the first transcriptomics studies in
I. laevigata[
21]. The dataset also allowed us to perform genome-wide identification of TFs. In this study, we report the identification of 68 WRKYs in
I. laevigata, which are distributed to the three groups that regulate multiple biological processes. We further experimentally studied the involvement of
IlWRKY22 in the content of both flowering and abiotic stress response.
WRKY22 has been shown to play an important role in the defense against aphids[
28], bacteria[
29], low temperature[
30], flowering[
31] and senescence[
32]. Our data demonstrated that
IlWRKY22 is a negative regulator in flowering and a positive regulator in resistance to salt and drought. Thus,
IlWRKY22 is a multifaceted TF that can be tweaked to improve multiple traits.
3. Discussion
Using our
I. laevigata transcriptome data from the flower, we identified a total of 68 WRKYs and classified them into three groups: 22 in Group I, 35 in Group II, and 11 in Group III. These three groups of WRKYs may participate in various biological processes or pathways based on previous functional studies in multiple plant species (
Table S2). Group I WRKYs are involved in plant disease resistance, abiotic stress (e.g., salt, cold, and heat stress) responses, senescence, and developmental processes (e.g., pollen development and seed coat development). Group II genes are mainly involved in defense responses to bacterial and fungal infection, abiotic stresses (salt, cadmium, drought stress) response, senescence, and developmental processes (e.g., lignin synthesis, pollen development, and flowering). They also participate in auxin transport and jasmonic acid signaling pathways. The biological function of Group III is also diverse, including defense against bacteria and/or fungi, abiotic stress (e.g., ozone, temperature, and drought), senescence, development of lateral roots and leaves, and interaction with plant hormone pathways (salicylic acid, abscisic acid, brassinosteroids, ethylene hormone, and jasmonic acid). The functional redundancy of these WRKYs indicates their key roles in responding to both external and internal stimuli [
35]. It can also be reflected by the presence of redundant sequences in many plant species including
Xanthoceras sorbifolium[
36] and tomato[
37], possibly due to duplication during plant evolution. Indeed, several sub-genome duplication events may have led to the formation of a large WRKY family in most modern eukaryotes[
38].
Due to the large number of WRKYs in
I. laevigata, we selected
IlWRKY22 for functional studies because it is less explored compared to other
WRKYs. We first focused on its role in flowering and found that overexpressing
IlWRKY22 delays flowering in
Arabidopsis. This provides another line of evidence to support the key roles of WRKYs in controlling flowering time. For example, heterologous expression of the cotton
GbWRKY1 gene in
Arabidopsis promotes flowering by regulating the transcription of
SOC1[
39]. Similarly, an up-regulation of
BcWRKY22 in Chinese cabbage by low temperature promotes bolting and flowering via increasing
BcSOC1[
40]. In addition, a
WRKY11 mutant in rice causes a late flowering phenotype [
41]. Because flowering is a complex process that is sensitive to multiple environmental factors[
42], such data collectively indicates that WRKYs can sense and integrate perturbations in the environment and initiate dynamic regulation of flowering [
43].
Corresponding to the multiple environmental aspects affecting flowering, plants evolved multiple flowering transition pathways including the de-repression of CDF on
CO and
FT by the photoperiod pathway, the de-repression of
FLC on
FT in the vernalization, temperature, and autonomous pathways, as well as the de-repression of
SOC1 by degrading DELLA proteins in the GA pathway[
44]. To delineate the exact mechanisms in which
IlWRKY22 was involved, we quantified the expression of 10 genes in multiple pathways (Figure 4). Our data showed a significant up-regulation of
CO and
GA20OX, indicating that the primary networks downstream of
IlWRKY22 were the photoperiod pathway and the GA pathway. This in turn resulted in altered metabolism of other players for flowering control including
SPL3. A lower-level expression of
SPL3 in the transgenic plants may further lead to a down-regulation of
SOC1, which ultimately causes a delay in flowering. Interestingly,
SPL3 is a major effector in the aging pathway as overexpressing or knocking-out darktreated
AtWRKY22 led to accelerated and delayed senescence phenotypes, respectively, in
Arabidopsis[
45]. Thus, it raises the possibility that
IlWRKY22 regulates flowering transition via the aging pathway and that
IlWRKY22 is the converging point of multiple pathways including flowering and senescence.
In addition to flowering control, our study further revealed a role of
IlWRKY22 in responding to abiotic stress. Compared to controls, plants overexpressing
IlWRKY22 showed much-improved growth under both salt and drought stress (Figure 5). This is in line with the finding that
WRKY22 is significantly up-regulated by drought stress[
46]. The enhanced resistance to abiotic stress can be attributed to up-regulation in photosynthesis (higher chlorophyll content, net photosynthesis, and Fv/Fm) and water retention capacity (lower stomatal conductance, intercellular CO
2, and transpiration rate, Figure 6). Photosynthesis is an important indicator of the ability of plants to respond to abiotic stress[
47]. Under high-salt and drought conditions, plants lose water due to high osmotic pressure, leading to stomatal closing, disrupted chloroplast structures, decreased activity of chlorophyll enzymes, and eventually inhibition of photosynthesis[
48,
49,
50]. Our data suggested that the protective role of
IlWRKY22 involved multiple processes that can alleviate the adverse impact of abiotic stress.
Although the exact molecular mechanisms by which
IlWRKY22 protects plants against abiotic stress are not clear, one possible mechanism is via ROS regulation. This is evidenced by the observation that the transgenic plants overexpressing
IlWRKY22 accumulate less
, H
2O
2, and MDA compared to that of the control plants (Figure 7). The generation of ROS and secondary damage to the cell membrane including membrane lipid peroxidation under environmental challenges have been well-documented in plants[
51,
52]. By lowering the level of oxidative stress,
IlWRKY22 may protect plants from excessive damage to the cellular membrane under stress. In addition, numerous studies have shown the key role of antioxidant enzymes including SOD, PDA, and CAT in counterbalancing ROS in plants[
52,
53]. Consistent with these studies, we found that the enzyme activities of SOD, CAT, and POD were increased in
IlWRKY22-overexpressing plants under stress (Figure 7). In summary, both a decrease in ROS accumulation and an increase in antioxidant enzyme activity led to a lower level of oxidative stress under salt and/or drought treatment in transgenic plants, which contributed to a better photosynthetic capacity and enhanced resistance under unfavorable conditions.
Figure 1.
Phylogenetic analysis of WRKYs from I. laevigata (in red) and Arabidopsis (in black).
Figure 1.
Phylogenetic analysis of WRKYs from I. laevigata (in red) and Arabidopsis (in black).
Figure 2.
Sequence phylogenetic analysis and subcellular localization of IlWRKY22. A IlWRKY22 was aligned with its homologs from Zingiber officinale (XP 042466882.1), Zingiber officinale (XP 042374628.1), Carex littledalei (KAF3333198.1), Lolium rigidum (XP 047076847.1), Lolium perenne (XP 0512272831), Dioscorea cavenensis subsp.rotundata (XP 039124149.1), Dendrobium catenatum (XP 020684026.1), Apostasia shenzhenica (PKA59051.1), Gastrodia elata (UUG47650.1), Nymphaea colorata (XP 031504405.1), Nvmphaea thermarum (KAF3792846.1), Cinnamomum micranthum f. kanehirae (RWR80315.1), Cocos nucifera (KAG1367745.1), Elaeis guineensis (XP 010924759.1), Phoenix dactvlifera (XP 0087924072), Gossypium raimondii (XP 0124900221), Durio zibethinus (XP 022757461.1), Cajanus cajan (XP 020220494), and Sesamum indicum (XP 0110929901). The red box was the highly conserved WRKYGQK domain. B A phylogenetic tree of IlWRKY22. C Subcellular localization of IlWRKY22. From left to right: Green fluorescence of GFP, cells under bright field, and cells under superposition of GFP and bright field. Scale bar was 20 μm.
Figure 2.
Sequence phylogenetic analysis and subcellular localization of IlWRKY22. A IlWRKY22 was aligned with its homologs from Zingiber officinale (XP 042466882.1), Zingiber officinale (XP 042374628.1), Carex littledalei (KAF3333198.1), Lolium rigidum (XP 047076847.1), Lolium perenne (XP 0512272831), Dioscorea cavenensis subsp.rotundata (XP 039124149.1), Dendrobium catenatum (XP 020684026.1), Apostasia shenzhenica (PKA59051.1), Gastrodia elata (UUG47650.1), Nymphaea colorata (XP 031504405.1), Nvmphaea thermarum (KAF3792846.1), Cinnamomum micranthum f. kanehirae (RWR80315.1), Cocos nucifera (KAG1367745.1), Elaeis guineensis (XP 010924759.1), Phoenix dactvlifera (XP 0087924072), Gossypium raimondii (XP 0124900221), Durio zibethinus (XP 022757461.1), Cajanus cajan (XP 020220494), and Sesamum indicum (XP 0110929901). The red box was the highly conserved WRKYGQK domain. B A phylogenetic tree of IlWRKY22. C Subcellular localization of IlWRKY22. From left to right: Green fluorescence of GFP, cells under bright field, and cells under superposition of GFP and bright field. Scale bar was 20 μm.
Figure 5.
Overall plant phenotypes under abiotic stress at 0, 7 and 14 days post treatment. A Top view and elevation of the growth state of WT, EV, and OE lines of tobacco at 0, 7, and 14 days in the normal growth state. B Top view and elevation of the growth status of WT, EV, and OE lines of tobacco at 0, 7, and 14 days under natural drought stress. C Top view and elevation of the growth status of WT, EV, and OE lines of tobacco at 0, 7, and 14 days under NaCl stress at 300mM. The scale bar is 20 cm.
Figure 5.
Overall plant phenotypes under abiotic stress at 0, 7 and 14 days post treatment. A Top view and elevation of the growth state of WT, EV, and OE lines of tobacco at 0, 7, and 14 days in the normal growth state. B Top view and elevation of the growth status of WT, EV, and OE lines of tobacco at 0, 7, and 14 days under natural drought stress. C Top view and elevation of the growth status of WT, EV, and OE lines of tobacco at 0, 7, and 14 days under NaCl stress at 300mM. The scale bar is 20 cm.
Figure 6.
Impact of overexpressing IlWRKY22 on photosynthetic and transpiration rate indexes under abiotic stress. Net photosynthetic rate, total chlorophyll content, Fv/Fm, stomatal conductance, intercellular CO2 content, and transpiration rate at 0, 7, and 14 days under 300 mM NaCl (A-F) or drought (G-L) were shown. Different letters above the bars indicated significant differences (p < 0.05).
Figure 6.
Impact of overexpressing IlWRKY22 on photosynthetic and transpiration rate indexes under abiotic stress. Net photosynthetic rate, total chlorophyll content, Fv/Fm, stomatal conductance, intercellular CO2 content, and transpiration rate at 0, 7, and 14 days under 300 mM NaCl (A-F) or drought (G-L) were shown. Different letters above the bars indicated significant differences (p < 0.05).
Figure 7.
Involvement of ROS regulation. MDA, O2-, H2O2 and the activities of SOD, CAT and POD, as well as leaf disc staining at 0, 7, and 14 days under 300 mM NaCl (A-F, and M) or drought (G-L, and N) were shown. Different letters above the bars indicated significant differences (p < 0.05).
Figure 7.
Involvement of ROS regulation. MDA, O2-, H2O2 and the activities of SOD, CAT and POD, as well as leaf disc staining at 0, 7, and 14 days under 300 mM NaCl (A-F, and M) or drought (G-L, and N) were shown. Different letters above the bars indicated significant differences (p < 0.05).
Figure 8.
Expression of abiotic stress-responsive genes. The expression levels of NtCAT, NtHAK1, NtPMA4, NtPOD, NtSOD, and NtSOS1 at 0, 7, and 14 days under 300 mM NaCl (A-F) or drought (G-L) were shown. Different letters above the bars indicated significant differences (p < 0.05).
Figure 8.
Expression of abiotic stress-responsive genes. The expression levels of NtCAT, NtHAK1, NtPMA4, NtPOD, NtSOD, and NtSOS1 at 0, 7, and 14 days under 300 mM NaCl (A-F) or drought (G-L) were shown. Different letters above the bars indicated significant differences (p < 0.05).