Cloning and bioinformatics analysis of HcERF4
The ERF transcriptional factor, HcERF4, was cloned from kenaf according to the transcriptome data (Chen et al. 2020). The open reading framer (ORF) is 702bp in length, encoding 233 amino acids, containing a typical AP2 motif from 25 to 57. (Fig. 1A). The molecular weight of the protein is 25.53 KDa and the theoretical pI value is 8.89. The average hydrophilic coefficient (GRAVY) of HcERF4 is -0.640, indicating it is a hydrophilic protein. The instability coefficient of this protein is 65.71, which means it is an unstable protein. To understand the phylogenetic relationships between HcERF4 and its homologous proteins, multiple sequence alignment was carried out using the Jalview program, the result showed that they shared high sequence similarity, especially the AP2/ERF domain (Fig. 1B). The three-dimensional structure predicated by Phyre2 showed that HcERF4 has one trans-membrane α-helices (Fig. 1C). Phylogenetic analysis based on amino acid sequence of HcERF4 and homologs from other species revealed that the HcERF4 gene has the closest homology relationship with two ERF4 proteins in Hibiscus syriacus and Gossypium hirsutum which belong to the same family of Malvaceae (Fig. 1D).
HcERF4 is a nuclear-localized protein
Firstly, the online software WoLF PSORT was used to predict subcellular localization, and the results showed that HcERF4 is localized in the nucleus. In order to further determine whether HcERF4, as an ERF-type transcription factor, is localized in the nucleus, we conducted transient expression experiments in tobacco leaf cells using an HcERF4 (green fluorescent protein) fusion vector. Compared with the epidermal cells transformed with an empty GFP vector alone, the GFP fluorescence signal was only observed in the nucleus with the HcERF4-GFP fusion construct (Fig. 2). These results determined that HcERF4 is a nucleus localized protein.
HcERF4 acts as a transcriptional activator in kenaf
The yeast test was used to study the transcriptional activity of the HcERF4 transcription factor (Fig. 3). The yeast transformants containing either pGBKT7 (negative control) or pGBKT7-HcERF4 grew normally on SD/Trp (single dropout supplements) medium. However, only the yeast transformed with pGBKT7-HcERF4 grew normally on SD/-Trp-His-Leu (triple dropout supplements) medium. When added with X-α-gal on the medium, the yeast transformed with pGBKT7-HcERF4 showed blue, the yeast transformed with pGBKT7 could not. These results indicate that the HcERF4 transcription factor has transcriptional activation activity.
Expression profile analysis of HcERF4 response to stress
qRT-PCR was employed to analyze the HcERF4 expression profile response to stress. Firstly, the HcERF4 expression level in different tissues was investigated, and the result indicated that HcERF4 was expressed extensively in all tested tissues including leaves, roots and stems under normal condition. HcERF4 has the highest expression level in leaves, followed by in roots, and the lowest in stems (Fig. 4A).
To further reveal its detailed expression model to salt or drought stress, kenaf leaves samples from different stress condition and different time were selected for qRT-PCR analysis. Under salt stress conditions, the expression of HcERF4 in kenaf leaves increased with the increase of salt concentration (Fig. 4B). With the duration of high salt concentration, there is a trend of first increasing and then decreasing. The expression level of HcERF4 under 100mM NaCl treatment reached the highest at 3h, which was 4.77 folds that of the control, and under 200mM NaCl stress, the expression level reached the highest at 1h, which was 17.05 folds that of the control, and the expression level gradually recovered after 12h. Under drought stress, with the increase of PEG concentration, the gene expression level showed a trend of first increasing and then decreasing (Fig. 4C). With the increase of stress time, under 10% PEG stress, the expression level reached its peak at 1 hour, which was 19.55 folds that of the control, and then increased with time. It showed a gradual downward trend, and the expression level showed an upward trend in 24h. Under 20%PEG treatment, the expression level reached the highest at 1h, and was 6.43 folds that of the control, and then gradually decreased with the increase of time. Our results show that the expression pattern of HcERF4 is different under different stress treatments. Under the same stress treatment, the expression patterns of HcERF4 at different concentrations are also different. These results indicate that HcERF4 may be involved in salt and drought signal response pathways.
VIGS (Virus-induced gene silencing) analysis of HcERF4
To reveal the role of HcERF4 under salt and drought stress in kenaf, VIGS experiment analysis were performed. The pTRV2-HcERF4 recombinant vectors were constructed and the positive plasmid was transferred into Agrobacterium GV3101, and the primers of HcERF4 fragment were used for PCR verification. A fragment of about 250 bp was obtained, indicating that the expression vector has been successfully transformed into Agrobacterium GV3101, and were injected into kenaf seedlings with the syringe infiltration method. After 14 days, infected leaves (sampled from the third true leaf) were selected for qRT-PCR analysis. As a represent example shown in Fig. 5, the result showed that the expression of HcERF4 in two of the five VIGS seedlings (2# and 5#) was significantly lower than that of the control, indicating that the HcERF4 was successfully silenced in these two kenaf seedlings (Fig. 5). The authentic silenced plant seedlings were used in subsequent phenotype and physiological analyses.
HcERF4 VIGS silenced kenaf reduced the tolerance to salt stress
In order to study the role of HcERF4 in response to salt stress, the HcERF4 silenced kenaf seedlings were subjected to stress treatments under 150mM NaCl stress for 7 days. As shown in Fig. 6, there was no significant difference in appearance of CK and pTRV2, while HcERF4 silenced plants showed significantly lower plant height, smaller leaf area, and less developed root system compare with CK and pTRV2 (Fig. 6A). The agronomic characteristics of VIGS plants including plant height, stem diameter, fresh weight and relative water content (RWC) were measured under salt stress condition. As shown in Fig. 6, the plant height (Fig. 6B), stem diameter (Fig. 6C), fresh weight (Fig. 6D) and RWC (Fig. 6E) of HcERF4 silenced plants were significantly lower than those of CK and pTRV2 plants. The results showed that HcERF4 VIGS silenced inhibited the growth of kenaf seedlings and reduced its salt tolerance.
The physiological indexes were analyzed and the results showed that the contents of MDA, O2− and H2O2 in pTRV2-HcERF4 plants were significantly increased by 64.3%, 43.9% and 21.6% under salt stress, while the contents of MDA, O2− and H2O2 in pTRV2 plants have no major changes (Fig. 7A, B, C). Silenced plants under salt stress The CAT and SOD were significantly reduced by 65.6% and 91.3% (Fig. 7D, F), and POD activity was significantly increased by 69.8% (Fig. 7E). Salt stress also significantly increased the proline content in silent plants (Fig. 7G). These results indicate that HcERF4 silent plants are sensitive to salt stress.
Silencing of HcERF4 in kenaf reduced the tolerance to drought stress
In order to study the role on the drought tolerance, HcERF4 VIGS silenced seedlings were treated with 15% PEG6000 for 7 days. By observing the phenotypes, it was found that there was no significant difference between CK and pTRV2 plants. However, HcERF4 silenced plants were significantly lower than CK and pTRV2 plants, and the leaves are curled and wilted shapes, the root system is under developed (Fig. 8A). This phenomenon indicates that the silence of VIGS reduces the drought tolerance of kenaf.
As shown in Fig. 8, after the drought stress treatment, the agronomic indicators of the treated plants were determined. Among them, CK and pTRV2 had no significant differences in plant height (Fig. 8B), stem diameter (Fig. 8C), fresh weight (Fig. 8D) and relative water content (Fig. 8E). However, after the VIGS silencing, the indicators decreased with exception of stem diameter, indicating that the HcERF4 may positively regulates the drought resistance in kenaf.
To study the effect of HcERF4 on drought tolerance, the content of MDA content and antioxidant activity were determined. The content of MDA can reflect the degree of damage to cell membranes caused by reactive oxygen species. Under drought stress, the MDA content of silenced plants was significantly reduced, but the MDA content of pTRV2 plants did not change significantly (Fig. 9A). After drought stress, the O2− and H2O2 content of silent plants increased by 1.28 and 1.31 folds, respectively, compared with the control (Fig. 9B, C), and there was no significant change in pTRV2 plants. Further study the oxidative stress under drought stress by analyzing the activity of antioxidant enzymes. The results showed that drought stress had a significant effect on the antioxidant enzyme activity of silenced plants. SOD is the first line of defense against ROS-mediated oxidative stress. After drought stress, the SOD activity of silent plants was significantly reduced, about 5.98 folds lower than that of the control, and the difference of pTRV2 under 15% PEG stress was not significant (Fig. 9D). Under drought stress, POD activity decreased significantly, about 1.48 folds lower than the control (Fig. 9E). CAT activity and POD activity showed similar trends. Compared with the control, the CAT activity of silent plants under drought stress was significantly reduced by 3.08 folds (Fig. 9F). However, the chlorophyll and proline contents of silent plants under drought stress were not significantly different from those of control plants (Fig. 9G). These results indicated that the decline of POD and SOD activity after HcERF4 silencing led to the increase in ROS levels, thereby aggravating membrane damage and reducing the drought tolerance in kenaf.
Histochemical detection of O2− and H2O2
Additionally, in order to analyze the accumulation of ROS in the HcERF4-silenced and control plants after salt stress, NBT and DAB staining were used to detect O2− and H2O2 levels. As shown in Fig. 10, after 7 days of salt stress or drought stress, the NBT and DAB stained area of silent plants was significantly larger than that of wild-type plants. In addition, the O2−and H2O2 content in the silent plants was significantly higher than that of wild-type plants (Fig. 7B-C and Fig. 9B-C). These results indicated that HcERF4 silenced plants aggravated its susceptibility to salt and drought stress.