Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 25 May 2022
Sec. Plant Abiotic Stress
This article is part of the Research Topic Salt Tolerance: Molecular and Physiological Mechanisms and Breeding Applications, Volume II View all 20 articles

A Cys2His2 Zinc Finger Transcription Factor BpSZA1 Positively Modulates Salt Stress in Betula platyphylla

\r\nXi ZhangXi ZhangQing GuoQing GuoLinlin QinLinlin QinLi Li
Li Li*
  • State Key Laboratory of Tree Genetics and Breeding, Northeast Forestry University, Harbin, China

Zinc finger proteins (ZFPs) are widely involved in plant growth and abiotic stress responses, however, few of these proteins have been functionally characterized in tree species. In this study, we cloned and characterized the BpSZA1 gene encoding a C2H2-type ZFP from Betula platyphylla. BpSZA1 is a transcription factor localized in the nucleus, with a transcription activation domain located at the N-terminus. BpSZA1 was predominantly expressed in stems and was induced by salt. We generated transgenic birch lines displaying overexpression (OE) or RNAi silencing (Ri) of BpSZA1 and exposed these along with wild-type birch seedlings to salinity. Phenotypic and physiological parameters such as superoxide dismutase, peroxisome, H2O2 content, proline content, water loss rate, and malondialdehyde content were examined. Overexpression of BpSZA1 in birch conferred increased salt tolerance. Chromatin immunoprecipitation-qPCR and RNA-seq showed that BpSZA1 binds to the GAGA-motif in the promoter of downstream target genes including BpAPX1, BpAPX2, BpCAT, and Bp6PGDH to activate their transcription. BpSZA1 also participates in abscisic acid (ABA) biosynthesis, proline biosynthesis, and the ABA/jasmonic acid pathway to enhance the salt stress of B. platyphylla.

Introduction

Plants typically suffer multiple abiotic stresses during their life cycle, which limit their growth and development. Climate extremes and modern agricultural practices are leading to continuous increases in soil salinity, making this a major abiotic stress (Abe et al., 2003; Franco-Zorrilla et al., 2014; Shrivastava and Kumar, 2015; Zhu, 2016). High salinity causes osmotic salt stress by reducing water uptake; it also accelerates the production of active oxygen radicals, leading to oxidative bursts that can damage or even kill plants (Munns and Tester, 2008; Franco-Zorrilla et al., 2014). To cope with salt stress, plants have evolved various defense mechanisms. Organic solutes including sucrose, proline, and glycinebetaine, accumulate in the cytosol and organelles to balance the osmotic potential of vacuolar Na+ (Hasegawa et al., 2000; Seki et al., 2003). In addition, antioxidant enzymes such as superoxide dismutase (SOD), peroxidase (POD), ascorbate peroxidase (APX), and catalase (CAT) and non-enzymatic compounds (ascorbate and reduced glutathione) are essential for reactive oxygen species (ROS) homeostasis, detoxifying ROS and increasing ROS scavenging ability (Franco-Zorrilla et al., 2014; Gupta and Huang, 2014). Those defense mechanisms are known to be regulated by specific transcription factors (TFs) including APETALA2/ethylene-responsive factor (AP2/ERF), v-myb avian myeloblastosis viral oncogene homolog (MYB), WRKY, basic region-leucine zipper (bZIP), and zinc finger proteins (ZFPs) (Riechmann et al., 2000; Seki et al., 2003; Guo et al., 2005; Iida et al., 2005; Wang et al., 2019). Thus, TFs serve as important participants in plant responses to salt stress by combining with various cis-elements in the promoter regions of downstream genes to modify their expression (Abe et al., 2003; Franco-Zorrilla et al., 2014).

Zinc finger proteins are found widely in animals, plants, and microorganisms (Miller et al., 1985). Plant ZFPs are divided into nine subfamilies, C2H2, C2HC, C4, C6, C8, C2HC5, C3HC4, C4HC3, and CCCH depending on the number and position of cysteine (Cys) and histidine (His) residues (Berg and Shi, 1996). Plant C2H2 type ZFPs from one of the largest subfamilies of ZFPs and play a crucial role in plant development (Englbrecht et al., 2004; Takatsuji, 1998). STAMENLESS 1, a C2H2 ZFP from rice, positively modulated floral organ identity by participating in transcriptional regulation of SPW1/OsMADS16 (Xiao et al., 2009). Also in rice, NON-STOP GLUMES1 regulated Spikelet development by inhibiting the expression of DROOPING LEAF (DL), LONGSTERILELEMMAS 1 (G1), and MOSAIC FLORAL ORGANS 1 (MFO1) (Zhuang et al., 2020). In wheat, Tipped1, containing an EAR motif, acts as a transcriptional repressor to regulate awn development; lines with a high expression level of Tipped1 displayed awn inhibition (Huang et al., 2020). Several studies also suggest that plant C2H2 ZFPs are involved in responses to salt stresses (Kiełbowicz-Matuk, 2012; Wang et al., 2019). For example, the heterologous expression of ZjZFN1 enhances the salt tolerance of transgenic Arabidopsis plants (Teng et al., 2018). Tomato SlZF3 interacts with CSN5B to confer greater tolerance to salt treatment in lines overexpressing SlZF3 (Li et al., 2018). OsZFP213, a C2H2 ZFP from wheat, confers salt tolerance by interacting with OsMAPK3 (Zhang et al., 2018). In a recent study, overexpression of GhZAT34 and GhZAT79 in cotton enhanced salt tolerance (Rehman et al., 2021). Taken together, these studies demonstrate that C2H2 ZFPs play an important role in plant development and the response of plants to salt stress.

To screen differentially expressed genes (DEGs) at the transcription level, RNA-seq can be used and requires neither genome annotation nor pre-synthesized nucleotides as probes. More and more studies have identified genes differentially expressed under abiotic stress using RNA-seq. Evidence suggests that there are 97 putative WRKY proteins in pearl millet, with PgWRKY33, PgWRKY62, and PgWRKY65 responding to dehydration and salinity stress (Chanwala et al., 2020). Similarly, 132 NAC proteins were identified in peanuts using RNA-seq; these were classified into eight subgroups, where the members of groups IV, VII, and VIII are induced by drought and abscisic acid (ABA) stresses (Li et al., 2021). In addition, genes encoding 197 bHLH proteins have been identified in pear, and PbrbHLH7, PbrbHLH8, PbrbHLH128, PbrbHLH160, PbrbHLH161, and PbrbHLH195 are up-regulated under cold and drought stress (Dong et al., 2021). RNA-seq is thus a useful tool for identifying TFs that respond to abiotic stress.

Chromatin immunoprecipitation (ChIP) is used to detect interactions between TFs and DNA in vivo and can also identify downstream genes regulated by TFs (Nelson et al., 2006). Recent work established that, in abiotic stress, WRKY33 and WRKY12 can form a complex under abiotic stress increasing the expression of RAP2.2 to enhance hypoxia tolerance in Arabidopsis (Tang et al., 2021). Similarly, the ChIP-seq assay combined with RNA-seq revealed that SlGRAS4 forms a homodimer to activate its transcription as well as bind directly to the SlCBF promoter to enhance cold tolerance (Liu et al., 2020). ZmNST3, a member of the NAC TF family, directly combines with the promoter of GST/GlnRS to activate its transcription and increase ROS scavenging, thereby enhancing drought tolerance in maize (Ren et al., 2020). Currently, however, there are limited reports on the application of ChIP technology to forests species.

Birch (Betula platyphylla Suk.) is widely distributed in cool-temperate and boreal forests in the northern hemisphere (Kazuhiko, 2007) and is widely used in building and wooden furniture (Li et al., 2002). However, birch is more sensitive to salt than Ulmus pumila and Fraxinus mandshurica. So it is important to study the salt-tolerance mechanisms of B. platyphylla for cultivating salt-tolerant birch varieties. In this study, we cloned a gene encoding a C1-2i subfamily member of the C2H2 ZFPs from B. platyphylla that is induced by salt. The gene was named BpSZA1, and we generated lines overexpressing (OE) or showing RNAi silencing (Ri) lines of BpSZA1 for gain- and loss- of function analysis. BpSZA1 conferred increased tolerance to salt stress in OE lines compared with wild-type (WT) and Ri lines by participating in ROS scavenging, ABA biosynthesis, proline biosynthesis, and jasmonic acid (JA) responses. Our results provide useful information on the function of C2H2 ZFPs in response to abiotic stress.

Materials and Methods

Plant Materials and Growth Conditions

Birch (B. platyphylla Suk.) plantlets were cultured in vitro on solid Lloyd and McCown’s Woody Plant Basal Salts (WPM) medium as previously described (Wang et al., 2016). After regeneration, plantlets were transplanted into individual pots containing a mixture of perlite/soil/vermiculite (1:4:1, by vol) for cultivation for 8 weeks in a greenhouse under controlled conditions (16/8 h light/dark, 25 ± 1°C, and 70–75% relative humidity) for tissue-specific expression and salt stress. For tissue-specific expression, the following tissues were harvested: including YL: young leaves, leaves from the first to second internodes; AL: adult leaves, leaves from the third to fourth internodes, OL: old leaves, leaves from the fifth to sixth internodes; AS: adult stems, stems with the third to fourth internodes; OS: old stems, stems with the fifth to sixth internodes. For salt stress, the root of 8-week-old birch plants prepared as described above was watered directly with 200 mM NaCl solution, and samples were harvested after treatment for 0, 6, 12, 24, 48, and 72 h. Three independent biological replicates were performed. Leaves from the birch seedlings were harvested for RNA isolation.

RNA Isolation and cDNA Synthesis, qRT-PCR

Total RNA was isolated from each sample using a Universal Plant Total RNA Extraction Kit according to the manufacturer’s instructions (BioTeke Corporation, China). qRT-PCR analysis, assays were performed using an SYBR Premix Ex Taq II (TaKaRa, Dalian, China) kit according to the manufacturer’s specifications. Quantification was performed using the 2–ΔΔt method, where the birch Tubulin gene BpTubulin is a housekeeping gene for normalization (Regier and Frey, 2010). All primers used in this study are detailed in Supplementary Table 1.

Bioinformatics Analysis of BpSZA1

A previously released dataset of B. platyphylla (Wang et al., 2018) was analyzed to identify C2H2 ZFP genes exhibiting induced by salt stress. Birch_GENE_10017081 was selected for further analysis because its expression was upregulated at all-time points. Protein information was obtained from the amino acid sequence encoded by the gene and was analyzed using ExPASY.1 The ZFP sequences in Arabidopsis, Populus pilosa, and Betula pendula subsp. were obtained from TAIR,2 the phytozome database,3 and the CoGe database,4 respectively. A phylogenetic tree was reconstructed in MEGA7 using the maximum likelihood method with 1,000 bootstraps replicates. Conserved motifs in BpSZA1 were predicted using Multiple Expectation maximizations for Motif Elicitation (MEME5).

Binary Vector Generation and Plant Transformation

Overexpression and RNAi constructs were generated using the BP and LR Gateway cloning system following the manufacturer’s instructions (Thermo Fisher Scientific, United States). The full-length of BpSZA1 CDS (open reading frame) without termination codon was amplified using gene-specific primers for the overexpression construct. To generate the RNAi (silencing) vector, a pair of primers flanking the CDS sequence from position 532–741 bp was used. Amplified fragments were first cloned into the pENTR™/SD/D-TOPO™ vector using BP-cloning and then into the binary vectors pGWB5 (overexpression) or pB7GWIWG2 (I) through LR reactions (Karimi et al., 2002). After sequence verification, the binary plasmids were transformed into Agrobacterium tumefaciens strain EHA105 using the freeze/thaw method. Agrobacterium-mediated transformation was performed as previously described (Holsters et al., 1978). Transgene presence was verified using the PCR primers LR-35S-F(CCTCGGATTCCATTGCCCAGCTA) and LR-GFP R(GTCGATGCCCTTCAGCTCGAT). Sequences of all primers used for the cloning are shown in Supplementary Table 2.

Subcellular Localization Assay

The full-length BpSZA1 coding sequence (without the stop codon) was amplified and cloned into a pBI121-eGFP vector upstream of the eGFP sequence and behind the CaMV 35S promoter, according to the reading frame, to generate 35S:BpSZA1-GFP. The method of tobacco subcellular localization is based on the method published by Sparkes et al. (2006). The 35S:BpSZA1-GFP plasmid was transferred to A. tumefaciens GV3101 competent cells and transformed into Nicotiana benthamiana; 35S-GFP was used as a control. After 2–3 days, fluorescence in the tobacco cells was observed at 488 nm under a confocal laser scanning microscope LSM700 (Zeiss, Jena, Germany). Sequences of all primers used for the cloning are shown in Supplementary Table 3.

Transactivation Assay

To examine the transactivation activity of BpSZA1 and determine the minimal domain required for activation, the complete BpSZA1 CDS along with seven truncated BpSZA1 fragments were individually fused to the GAL4 DNA binding domain in the pGBKT7 vector using the In-Fusion HD Cloning System CE (TaKaRa, Dalian, China). The recombinant plasmids and the pGBKT7 vector (negative control) were individually transformed into yeast strain Y2HGold. Transformants were successfully dropped onto SD/-Trp and SD/-Trp/-His/-Ade/X-α-gal selective medium and incubated at 30°C for 3–5 days. The transactivation activity was evaluated according to the growth status of colonies and the blue color manifested by X-α-gal activity. Sequences of all primers used for cloning are shown in Supplementary Table 4.

Salt Stress Treatments

For salt stress treatments, 8-week-old WT and transgenic B. platyphylla seedlings, three OE lines (OE2, OE6, and OE9) and three Ri lines (Ri5, Ri7, and Ri8), were treated with 200 mM NaCl for 7 days and the phenotypes were photographed.

For plant height, fresh weight, and root length measurements, 7-day-old birch seedlings (WT, OE2, OE6, OE9, Ri5, Ri7, and Ri8) were cultured in half MS minimal medium containing 40 mM NaCl. After 4 weeks of growth, plant height, fresh weight, and root length were measured. Birch seedlings grown in half MS minimal medium were used as a control.

Analysis of Physiological Parameters

Similar-sized birch seedlings were treated with 200 mM NaCl for 3 days, and adult leaves (WT, OE6, OE9, Ri5, and Ri7) were harvested. SOD activity was measured using a SOD test kit (Nanjing Jiancheng Bioengineering Institute, A001-1), POD activity was measured using a POD test kit (Nanjing Jiancheng Bioengineering Institute, A084-3), malondialdehyde (MDA) content was measured as described previously (Heath and Packer, 1968), proline content was determined according to a previously published protocol (Bates, 1973), and total chlorophyll content was assayed using a method previously published (Minotti et al., 1994). For electrolyte leakage analysis, leaf samples of B. platyphylla were collected using a 1-cm punch with 10 pieces for each sample, as previously described (Ben-Amor et al., 1999). Leaves of B. platyphylla were treated with 200 mM NaCl for 6 h, and leaves were then detached and stained with 3,3-diaminobenzidine (DAB) or nitro blue tetrazolium (NBT) solution to detect H2O2 or O22–, respectively, according to a previously published method (Fryer et al., 2002).

Transcriptome Determination and Analysis

Leaves of 8-week-old B. platyphylla OE5, Ri7, and WT seedlings treated with 200 mM NaCl were rapidly collected after 3 days of treatment. Six samples with three repeats were used for RNA sequence analysis—Ri before/OE before, Ri before/WT before, WT before/OE before, Ri after/OE after, Ri after/WT after, and WT after/OE after—where “before” indicates samples collected before salt stress and “after” indicates samples collected after salt stress. A total of 1 μg RNA per sample was used as input material for RNA sample preparations. Sequencing libraries were generated using a NEBNext Ultra™ RNA Library Prep Kit for Illumina (NEB, United States) following the manufacturer’s recommendations. Clustering of index-coded samples was performed on a cBot Cluster Generation System using a TruSeq PE Cluster Kit v4-cBot-HS (Illumina) according to the manufacturer’s instructions. After cluster generation, library preparations were sequenced on an Illumina platform and paired-end reads were generated. Adaptor sequences and low-quality sequence reads were removed from the data sets. Raw sequences were transformed into clean reads after data processing. Hisat26 tools software was used to map reads to the reference birch genome.7 Transcripts were assembled using StringTie.8 Differential expression analysis of pairs of samples was performed using EBseq.9 The false discover rate (FDR) < 0.01 and |log2(fold change)| ≥ 2 were set as the threshold for significantly differential expression. Gene Ontology (GO) enrichment analysis of DEGs was implemented using the GOseq R packages based on Wallenius non-central hypergeometric distribution (Young et al., 2010), which can adjust for gene length bias in DEGs.

Chromatin Immunoprecipitation-PCR and Chromatin Immunoprecipitation-qPCR Analysis

To identify putative target genes of BpSZA1, ChIP-PCR, and ChIP-qPCR were used. Sonication conditions of the samples were optimized according to a previously described method (Li et al., 2014). Briefly, approximately 2 g pGWB5-BpSZA1-OE birch seedlings were immersed in 1% (w/v) formaldehyde for protein cross-linking and protein binding to DNA. Next, ChIP was followed by extraction and shearing of the chromatin precipitated using a ChIP analysis kit (P2078, Beyotime) according to the manufacturer’s instructions. An anti-GFP antibody (AG281, Beyotime) (ChIP+) was then used and chromatin was immunoprecipitated; an anti-IgG2a antibody was used as a negative control (Mock). Chromatin without immunoprecipitation was used as the input control (Input). Three biological replicates were carried out. All primers used for ChIP-PCR and ChIP-qPCR analysis are detailed in Supplementary Table 6.

Statistical Analysis

All experiments were conducted with at least three biological replicates unless otherwise mentioned, and the standard error of the mean was computed in each case. The Student’s t-test was performed for the estimation of statistical significances (SPSS 18; IBM Corp., Armonk, NY, United States). Data points representing statistically significant differences between WT and transgenic lines or between control and stress conditions are indicated.

Results

Identification and Structural Analysis of BpSZA1

To identify valuable salt-induced genes from B. platyphylla and to cultivate salt-tolerant birch varieties, we used a next-generation sequencing technique were used to construct a cDNA library of B. platyphylla (200 mmol L–1 NaCl) (Wang et al., 2018). Analysis of RNA-seq data revealed 30 DEGs annotated as C2H2 ZFPs, of which 12 were down-regulated and 20 were up-regulated by salt, respectively. Among these salt-induced C2H2 ZFPs, birch_GENE_10017081, named BpSZA1 (GenBank ID: MZ544846), was steadily up-regulated at all-time points under salt stress compared with the control condition (Figure 2A). We used PCR to amplify the cDNA of BpSZA1, and sequencing results indicated that the full-length BpSZA1 is 744 bp without introns, encoding a 247 aa protein with a predicted mass of 26.24 kDa and a theoretical pI of 8.66. A phylogenetic tree was reconstructed using nucleotide sequences of 20 members of the ZFPs family from Arabidopsis thaliana and 13 from Populus trichocarpa available in the Phytozome database and 12 genes from B. pendula subsp. available in the CoGe database. BpSZA1 was found to cluster with AtAZF2 (AT3G19580), AtAZF3 (AT5G43170), AtZAT6 (AT5G04340), AtZAT10 (AT1G27730), and AtZAT13 (AT3G49930) from Arabidopsis (Figure 1A). Furthermore, we used these six sequences to predict conserved motifs using MEME, identifying the presence of five conserved motifs (Figure 1B). There were two zinc finger domains of QALGGH among the five motifs (motif 1 and motif 3). Specific motif sequences are shown in Supplementary Table 7.

FIGURE 1
www.frontiersin.org

Figure 1. Phylogenetic and bioinformatics analyses of BpSZA1 from Betula platyphylla. (A) Phylogenetic analysis of BpSZA1 and 45 members of the C2H2-type transcription factor from Arabidopsis thaliana (20 members, prefixed with At), Populus trichocarpa (13 members, prefixed with Potri), and Betula pendula Roth (12 members, prefixed with Bpev) by using the maximum likelihood method. (B) Conserved protein motifs of BpSZA1 and five C2H2-type transcription factor C1-2i subfamily members in Arabidopsis thaliana.

FIGURE 2
www.frontiersin.org

Figure 2. Expression of genes encoding zinc finger proteins in B. platyphylla under 200 mM NaCl salt stress, and BpSZA1 expression patterns in different tissues and under different treatments. (A) Bioinformatic analysis of the transcriptome data set showing expression levels for 30 genes encoding zinc finger proteins that were upregulated and downregulated after exposure of birch to 200 mM NaCl for 6, 12, 24, and 48 h. The color scale represents change. (B) Transcript levels of BpSZA1 in different tissues of B. platyphylla. YL, young leaf; AL, adult leaf; OL, old leaf; AS, adult stem; OS, old stem; R, root. (C) Relative expression levels of BpSZA1 were measured by RT-qPCR in response to NaCl. Values are mean ± SE (n = 15). Asterisks denote significant differences: *P ≤ 0.05; **P ≤ 0.01. Error bars represent standard error for three replicates.

Expression Analysis of BpSZA1

To verify the accuracy of RNA-seq, we treated birch seedlings with NaCl and examined BpSZA1 transcript levels. Under salt treatment (NaCl) (Figure 2C), BpSZA1 expression increased over the first 6 h, showing a 3.48-fold increase compared with expression at 0 h, and then gradually declined. This result verified the accuracy of RNA-seq and indicated that the expression of BpSZA1 is responsive to salt stress.

To investigate the spatial expression of BpSZA1, we used RT-PCR in different tissues (Figure 2B). BpSZA1 was ubiquitously expressed in all tissues, with the highest expression levels in adult stems, followed by old stems, roots, old leaves, and adult leaves. Expression levels of BpSZA1 were lowest among all sampled tissues in young leaves, which were used as a control. These results indicated that BpSZA1 is expressed specifically in different tissues.

Subcellular Localization of BpSZA1 Protein

To determine whether the C2H2 TF BpSZA1, is localized to the nucleus, we constructed the fusion vector 35s:BpSZA1: GFP and transfected this into tobacco epidermal cells, with the 35: GFP vector used as a control. Fluorescence signals in tobacco epidermal cells were observed under an LSM700 confocal laser microscope. The 35s:BpSZA1: GFP fusion protein was mainly visible in the nucleus, while the control (35s: GFP) was observed throughout the cells (Figure 3A). Hence, BpSZA1 is localized to the nucleus.

FIGURE 3
www.frontiersin.org

Figure 3. Subcellular localization and transactivation of BpSZA1. (A) Subcellular localization of BpSZA1 in transiently expressed tobacco leaves shows BpSZA1 is targeted to the nucleus. Bars = 50 μm. (B) Truncated sequences of the CDS of BpSZA1 were fused in-frame, respectively, to the GAL4 DNA-binding domain in pGBKT7 and transformed into Y2H yeast cells. The transformed cells were plated onto SD/-Trp (growth control) or SD/-Trp/-His/-Ade/X-a-Gal medium. pGBKT7 empty vector was used as a control.

Transcriptional Activity of BpSZA1

Transactivation activity of BpSZA1 was examined using a series of deletions of the BpSZA1 CDS fused with the GAL4 DNA-binding domain sequence in the pGBKT7 vector and transformed into Y2H Gold cells. Cells containing the full-length CDS of BpSZA1 grew equally well in both SD/-Trp (control) and SD/-Trp/-Leu/-Ade/X-a-Gal medium but transformed yeast cells containing the full-length BpSZA1 CDS turned blue in SD/-Trp/-Leu/-Ade/X-a-Gal medium, indicating the BpSZA1 is a TF. Meanwhile, assays using the deletion mutants showed that the transcription activation domain of the BpSZA1 protein is located between amino acids 212aa and 247aa at the C-terminus of the protein (Figure 3B).

Overexpression of BpSZA1 Enhanced Salt Tolerance

To compare their respective stress tolerance, 8-week-old WT and OE2, OE6, OE9, Ri5, Ri7, and Ri8 transgenic birch plants grown in soil, were exposed to 200 mM NaCl. Under control conditions, there were no significant differences in the appearance of OE lines, Ri lines, and WT plants (Figure 4A). However, after 7 days of salt stress, leaves of the Ri lines showed obvious withering and yellowing compared with their appearance before stress. Obvious leaf curling and wilting also appeared in the WT seedlings, but the leaves were greener and more stretched than those of the Ri lines. By contrast, the OE lines showed only a slight curl of the leaves under salt stress, and the leaves became pale (Figure 4A). Meanwhile, there were no significant differences in plant height, fresh weight, or root length in the different lines under control conditions (Figures 4B–D), indicating that BpSZA1 does not affect plant growth or development. Under salt stress, the plant height, fresh weight, and root length of Ri lines decreased by an average of 59, 56, and 63%, respectively, compared with measurement before stress with 48, 50, and 49% average decrease, respectively, in the WT. The reduction in plant height, fresh weight, and root length of OE lines were minimal under salt stress, being 28, 26, and 29%, respectively (Figures 4B–D). To compare growth under salt stress, we analyzed chlorophyll contents among all salt-treated lines. Under control conditions, all salt-treated lines had approximately equal chlorophyll contents (Figure 4E). However, after salt treatment, OE lines had higher chlorophyll contents compared with WT plants, but Ri lines had lower chlorophyll contents (Figure 4E). Therefore, overexpression of BpSZA1 improved the tolerance of birch to salt stress, while inhibition of BpSZA1 expression reduced the tolerance of birch to salt stress.

FIGURE 4
www.frontiersin.org

Figure 4. BpSZA1 confers salt stress tolerance. (A) Phenotypic differences following salt treatment. (B) Plant height. (C) Fresh weight. (D) Relative root length. (E) Chlorophyll content. Asterisks denote significant differences: *P ≤ 0.05. Error bars represent standard error for three replicates.

BpSZA1 Decreases Membrane Lipid Peroxidation and Increases Proline Content Under Salt Stress Conditions

To further determine differences in salt tolerance among the WT and OE2, OE6, OE9, Ri5, Ri7, and Ri9 transgenic birch plants, we measured electrolyte leakage, MDA content, and proline content. Under control conditions, there was no significant difference in electrolyte leakage rates among WT and transgenic lines. Nonetheless, after salt stress, however, electrolyte leakage rates were increased in all salt-treated lines. Electrolyte leakage rates of WT plants were significantly lower than those of Ri lines but were significantly higher than those of OE lines (Figure 5B). Meanwhile, there was no difference in MDA content among birch lines under control conditions. When exposed to salt stress conditions, MDA contents were increased in all lines; however, OE lines had the lowest MDA contents, followed by WT, with Ri lines having the highest MDA contents (Figure 5D). Similarly, the proline content of OE lines was greater than that of WT plants under salt stress conditions; however, the proline content of Ri lines was less than that of WT plants (Figure 5C). These results suggested that compared with WT plants, OE lines showed increased salt tolerance, and Ri lines demonstrated decreased salt tolerance.

FIGURE 5
www.frontiersin.org

Figure 5. Determination of malondialdehyde (MDA) content, electrolyte leakage, proline content, MDA content, ROS accumulation (Histochemical staining and H2O2 content), superoxide dismutase (SOD) activity, peroxidase (POD) activity. Birch plantlets were treated with 200 mM NaCl for 0 h (control) and 6 h, and then their leaves were detached and stained with NBT (A) and DAB (B) to detect H2O2 and O2. (B) Electrolyte leakage, (C) proline content, (D) MDA content, (E) H2O2 content, (F) SOD activity, and (G) POD activity in the leaves of WT, BpSZA1 OE, and Ri plants under normal, salt stress. Asterisks denote significant differences: *P ≤ 0.05. Error bars represent standard error for three replicates.

BpSZA1 Increases Reactive Oxygen Species Scavenging and Proline Content Under Salt Stress

We used DAB and NBT in situ staining to measure the accumulation of H2O2 and O22–. Staining results revealed that under control conditions, OE, WT, and Ri lines had similar H2O2 and O22– levels (Figure 5A). However, when exposed to salt stress conditions, we visualized stronger staining by NBT (for O2) or DAB (for H2O2) in Ri lines compared with WT, whereas OE lines showed less accumulation of H2O2 and O22– than WT plants (Figure 5A). To further clarify whether the histochemical staining differences were caused by ROS scavenging ability, we compared H2O2 content and the activities of some ROS scavenging enzymes, including POD and SOD, among the WT, OE, and Ri lines. After salt treatment, Ri lines showed higher H2O2 accumulation compared with WT plants, while OE lines had lower H2O2 accumulation (Figure 5E). However, there was no remarkable difference in H2O2 accumulation among OE lines, Ri lines, or the WT under control conditions (Figure 5E). At the same time, under control conditions, there was no difference in SOD or POD activities among any of the lines. However, when OE, Ri, and WT plants were treated with 200 mM NaCl, both SOD and POD activities were significantly lower in Ri lines compared with the WT but were significantly higher in the OE lines (Figures 5F,G). Taken together, these results indicated that overexpression of BpSZA1 increases ROS scavenging capability under salt stress.

Determination of the Genes Regulated by BpSZA1 on the Genome Level

To identify genes regulated by BpSZA1, we performed global transcriptomic profiling by RNA-seq on OE6, WT, and Ri7 lines. Plants were untreated plants (0 day) or salt-treated for 3 days (ST 3 days). The complete list of DEGs identified by pairwise comparison is given in Supplementary Table 8. Under normal growth conditions, 884 genes were up-regulated and 716 genes were down-regulated by BpSZA1 expression (with a P-value < 0.01 adjusted by the FDR) (Supplementary Figure 1A), whereas 1,119 genes were up-regulated and 1,841 genes were down-regulated by BpSZA1 expression under salt stress condition (with a P-value < 0.01 adjusted by the FDR) (Supplementary Figure 1B).

Gene Ontology term analysis revealed that, before salt stress, oxidation-reduction process, protein phosphorylation, and metabolic process were enriched among biological process GO terms. Among cellular components GO terms, integral components of membrane, nucleus, and plasma membrane terms were enriched. Among molecular function GO terms, ATP binding, metal ion binding, and protein serine/threonine kinase activity were enriched. After salt stress, oxidation-reduction process, protein phosphorylation, and metabolic process were also mainly enriched among biological process GO terms. Meanwhile, response to oxidative stress and salt stress terms were also enriched among biological process GO terms. Among cellular components GO terms, integral components of membrane, plasma membrane, and membrane terms were enriched. For molecular function GO terms, ATP binding and heme binding, and metal ion binding terms were enriched, meanwhile, oxidoreductase activity and peroxidase activity were also showed genes enrichment (Supplementary Figure 2). These results indicated that genes associated with these GO terms might be regulated by BpSZA1 in response to salt stress.

To identify the genes modulated by BpSZA1, we analyzed the expression levels of putative downstream genes in OE5, Ri7, and WT using RT-PCR. We selected the following genes according to physiological parameters and analysis of transcriptome data: BpSOD2, BpAPX1, BpAPX2, pyrroline-5-carboxylate synthase (BpP5CS), delta-1-pyrroline-5-carboxylate dehydrogenase (BpP5CDH) and 6-phosphogluconate dehydrogenase (Bp6PGDH), BpTIFY5A, Zeaxanthin epoxidase gene1 (BpZEP1), and BpCAT. All genes were up-regulated in the OE line except for TIFY5A and P5CDH, which were up-regulated in the Ri line (Figure 6). Genes related to ROS scavenging BpSOD2, BpAPX1, BpAPX2, BpCAT, and Bp6PGDH were significantly up-regulated after salt treatment, with the expression level of BpSOD2 after salt stress reaching as much as 58-fold of that under control conditions in the OE line (Figures 6A–E). Bp6PGDH expression in the OE lines was up-regulated by 41.4-fold under salt treatment (Figure 6C), while the expression of BpAPX1, BpAPX2, and BpCAT was up-regulated by 4.32-, 13.02-, and 7.05-fold, receptively, under salt treatment in the OE line (Figures 6A,B,D). Similarly, BpP5CS, which participated in proline biosynthesis was up-regulated by 34.18-fold under salt conditions in the OE line (Figure 6H), whereas expression of BpP5CDH, a participant in proline degradation, was up-regulated by 3.7-fold in the Ri line and remarkably down-regulated in OE the line after salt treatment (Figure 6G). Meanwhile, expression of BpTIFY5A, a suppressor of JA responses, was also up-regulated in the Ri line, reaching as much as 9.95-fold that under control conditions after salt treatment (Figure 6I). Notably, BpZEP1 is an important participant in ABA biosynthesis and was significantly down-regulated in the Ri line after salt stress, however, expression of BpZEP1 showed no difference under control or salt-treated conditions in the OE line and the WT (Figure 6F). These results showed that BpSZA1 has a positive function in response to salt stress. All primers used for RT-PCR analysis are detailed in Supplementary Table 5.

FIGURE 6
www.frontiersin.org

Figure 6. RT-PCR analysis of the expression levels of salt-responsive genes in Ri lines, WT, and OE lines before and after salt stress. (A) BpAPX1, (B) BpAPX2, (C) Bp6PGDH, (D) BpCAT, (E) BpSOD2, (F) BpZEP1, (G) BpP5CDH, (H) BpP5CS, and (I) BpTIFY5A. Asterisks denote significant differences: **P ≤ 0.01. Error bars represent standard error for three replicates.

Identification of Target Genes Directly Regulated by BpSZA1

To identify the motif that binds BpSZA1 for activating gene transcription following exposure to abiotic stress, we applied the MEME motif discovery tool (see text footnote 5). The promoters of 20 genes based on the qRT-PCR and RNA-seq data, that were up-regulated by BpSZA1 were used for further study. MEME results revealed that there was an 11-base conserved sequence present in most of the promoters studied (Figure 7B). The fourth to seventh bases of the conserved sequences appeared with the highest frequency and might be the core sequences of this motif; therefore, they were used for further experiments (Figure 7A). The 11-base conserved sequences were named the GAGA-motif (AGAGAGAGGGA, occupy the highest proportio). ChIP-PCR and ChIP-qPCR demonstrated that BpSZA1 protein bound to the GAGA-motif in the promoters of BpAPX1, BpAPX2, BpCAT, and Bp6PGDH (Figures 7C,D), indicating that the GAGA-motif combines with BpSZA1 to regulate the expression of BpAPX1, BpAPX2, BpCAT, and Bp6PGDH under salt stress. In addition, ChIP results indicated that genes associated with ROS scavenging ability, such as BpAPX1, BpAPX2, BpCAT, and Bp6PGDH, were predominantly regulated directly by BpSZA1.

FIGURE 7
www.frontiersin.org

Figure 7. (A) Multiple Expectation maximizations for Motif Elicitation (MEME) analysis of the conserved sequence GAGA-box present in the promoters of genes regulated by BpSZA1. (B) The schematic diagram shows the positions of the GAGA-box in the promoters. (C) Enrichment of the promoters of genes containing GAGA-box. (D) ChIP-qPCR analysis shows enrichment of the promoter sequences of genes after ChIP. ChIP+: the sonicated chromatin was immunoprecipitated with GFP antibody; mock: the sonicated chromatin was immunoprecipitated with IgG2a; input: the sonicated chromatin was used as a positive control. Three biological replicates were performed. Asterisks denote significant differences: **P ≤ 0.01. Error bars represent standard error for three replicates.

Discussion

BpSZA1 Belongs to C2H2 Zinc Finger Proteins and Confers Salt Tolerance

Multiple sequence alignment revealed that BpSZA1 clustered with AtAZF2, AtAZF3, AtZAT6, AtZAT10, and AtZAT13 from Arabidopsis (Figure 1B). Many studies have shown that important information on these genes plays an important role in the response to abiotic stress. According to reports in 2011, AZF1-OE and AZF2-OE lines in Arabidopsis show significant growth inhibition and sensitivity to abiotic stress (Kodaira et al., 2011). AtZAT6 not only participates in root development but also increases tolerance to cadmium in transgenic Arabidopsis (Chen et al., 2016). Among poplars, ZAT10 OE lines display increased cold and salt tolerance compared with WT plants (He et al., 2019, 2020). In the current study, we found that BpSZA1 is induced by salt. In addition, BpSZA1 plays a positive role in the response of birch plants to salt stress. Many physiological parameters, including POD activity (Figure 5F), SOD activity (Figure 5G), and proline content (Figure 4C), are modulated through the expression of BpSZA1. Overall, these results indicate that BpSZA1 belongs to C2H2 ZFPs and positively modulates salt stress in birch.

BpSZA1 Binds to GAGA-Motif to Regulate Gene Expression

Transcription factors usually regulate downstream gene expression by combining with cis-acting elements, so it is important to identify which cis-element combines with each TF when studying the function of TFs. Utilizing MEME, we found that BpAPX1, BpAPX2, BpCAT, and Bp6PGDH all contain the GAGA-rich motif sequence (Figure 7A). Further, ChIP-PCR and ChIP-qPCR analysis showed that BpSZA1 can directly regulate the expression of these genes by binding to GAGA-rich motifs in their promoters (Figures 7C,D). This result is consistent with those of previous studies. In Drosophila, the Tr1 gene encodes a protein with a C2H2 domain that can bind to a GAGA element to regulate the expression of downstream genes (Pedone et al., 1996). A similar mechanism may exist in plants. RAMOSA1, with a single C2H2 ZFP domain, can bind to the GAGA element to regulate the inflorescence structure of maize (Eveland et al., 2014). Furthermore, it was recently predicted that MaC2H2-4 and MaC2H2-5 in bananas can regulate the ripening of banana fruits by combining with different TF binding sites, including two GA-rich motifs (Kuang et al., 2021). In summary, these results show that BpSZA1 may regulate downstream gene expression by binding to the GAGA-repeat sequence.

BpSZA1 Enhanced Salt Tolerance by Directly Regulating Genes Related to Reactive Oxygen Species Scavenging

Abiotic stress including high salt and drought, causes apoplast oxidative stress in plants, accelerating the production of ROS including H2O2, ⋅O2, 1O2, and ⋅OH, which can damage and kill plants (Mittler et al., 2004; Suzuki et al., 2012; Baxter et al., 2014). However, the ROS generated can be scavenged by numerous enzymes. Major ROS-scavenging enzymes include SOD, APX, CAT, glutathione peroxidase (GPX), and peroxiredoxin (PrxR) in plants (Nakano and Asada, 1981; Rizhsky et al., 2004; Yang and Guo, 2018). In previous studies, ZAT10 OE lines displayed enhanced salt tolerance by directly regulating the expression of APX2 (Mittler et al., 2006; He et al., 2019). This result is consistent with our findings. BpSZA1 bound to BpAPX1 and BpAPX2 promoter regions containing the GAGA-motif, and expression of BpAPX1, BpAPX2 in the OE line was significantly increased after salt stress compared with that in WT and Ri line (Figures 6A,B), which may improve the salt tolerance of plants. Another important finding was that BpSZA1 could directly bind to the GAGA-motif of the BpCAT promoter region (Figures 7C,D). At the same time, expression levels of BpCAT in OE lines and the WT were highly up-regulated after salt treatment compared with those under control conditions, but the expression level in the Ri lines was lower (Figure 6D). This is also in accordance with earlier observations, which showed that NtERF172, a TF of the AP2 family, also binds to the NtCAT promoter region to activate NtCAT transcription, and overexpression lines of NtERF172 displayed greater ROS scavenging and better drought tolerance (Zhao et al., 2020). Furthermore, we found no significant difference in the expression of BpSOD2 among WT, OE, and Ri lines. However, after salt treatment, expression levels were 1.5-fold higher than that in the WT and 58-fold than in the Ri line (Figure 6E). This indicates that BpSZA1 could enhance salt tolerance by regulating the expression of BpSOD2.

Besides ROS-scavenging enzymes, the antioxidants ascorbic acid and glutathione can scavenge ROS produced in the stroma (Yang and Guo, 2018). The reduced coenzyme NADPH is an electron donor in the ascorbate-glutathione biosynthesis phase and is important for maintaining ascorbic acid and glutathione in a reduced state (Corpas et al., 1998). 6PGDH, G6DPH, malic enzyme, and isocitrate dehydrogenases are the main cellular sources of NADPH (Mittler et al., 2004; Marino et al., 2012; Suzuki et al., 2012; Baxter et al., 2014). We found that BpSZA1 could bind to the Bp6PGDH promoter region containing a GAGA-motif. The expression level of Bp6PGDH in OE lines after salt stress was 40-fold higher than that without salt stress. However, there was no significant difference in Bp6PGDH expression in WT or Ri lines between salt treatment and control conditions (Figure 6C). This finding has also been reported in previous research. Salinity-tolerant lines of barley have a higher expression level of 6PGDH (Witzel et al., 2010). Overexpression of BpSZA1 conferring tolerance to salt might increase the expression of Bp6PGDH. BpSZA1 might therefore enhance salt tolerance by regulating genes related to ROS scavenging, including BpAPX1, BpAPX2, BpCAT, BpSOD2, and Bp6PGDH.

BpSZA1 Improves Salt Tolerance by Participating in Abscisic Acid Synthesis, Proline Synthesis, and Activating the Abscisic Acid/Jasmonic Acid Pathways

The gene ZEP1 is a key gene for plant biosynthesis of ABA (Zhu, 2002). In 1998, the ZEP gene was reported to be involved in the biosynthesis of ABA precursors (Audran et al., 1998). Meanwhile, in a recent study, a ZEP gene MsZEP from alfalfa (Medicago sativa), overexpression of MsZEP enhances drought and salt tolerance by increasing ABA biosynthesis (Zhang et al., 2016). This is consistent with our study. After salt stress, the expression level of ZEP1 in the Ri lines was lower than that in WT plants; however, the expression level of ZEP1 in the OE lines and WT plants did not vary significantly before and after the stress. This shows that the absence of BpSZA1 might inhibit the activity of the BpZEP1 enzyme and affect the ABA-dependent pathway response of plants to salt stress (Figure 6F).

Proline is a component of plant proteins and has a critical role in regulating cell redox potential (Kavi Kishor et al., 1995; Zhu, 2002). In this study, expression levels of some genes related to proline biosyntheses, such as BpP5CS (Figure 6H) and BpP5CDH (Figure 6G), were significantly increased in OE lines after salt stress. At the same time, proline content was also increased after salt stress in OE lines, balancing the osmotic pressure inside and outside the cell and protecting the cell from damage. This indicates that BpSZA1 positively regulates proline biosynthesis to improve the salt tolerance of plants.

It is well documented that because ZFPs contain an EAR motif, they may act as transcriptional repressors in regulating the tolerance of plants to abiotic stress (Ciftci-Yilmaz et al., 2007; Xie et al., 2019; Wang et al., 2020). We found that only the expression level of TIFY5A/JAZ8 in Ri lines increased instantaneously after salt stress, indicating that BpSZA1 might inhibit the expression of TIFY5A/JAZ8. TIFY5A belongs to the JAZ subfamily of the TIFY TF family and is also known as JAZ8 (Shyu et al., 2012). It has previously been observed that ZAT10 can regulate the expression of TIFY10a/JAZ1 in the JA signaling cascade (Pauwels et al., 2008). When not responding to the JA signal, JAZ proteins combine with and inhibit downstream TFs expression to suppress JA responses. After being induced by JA signals, JAZ proteins are ubiquitinated and degraded through the 26S proteasome, thereby activating downstream TF transcription and JA responses (Fernandez-Calvo et al., 2011; Ruan et al., 2019; Yang et al., 2019). JA participates in plant salt responses. Extensive research has shown that the exogenous application of JA enhances salt tolerance by maintaining ROS or ion homeostasis (Qiu et al., 2014; Jiang et al., 2016; Farhangi-Abriz and Ghassemi-Golezani, 2018; Tavallali and Karimi, 2019). Overall, this indicates that inhibiting BpSZA1 expression may increase the expression level of BpTIFY5A and further suppress JA responses, thereby destroying ROS homeostasis and ion homeostasis and resulting in Ri lines showing more sensitivity to salt treatment (Figure 6I). Consequently, BpSZA1 may enhance the salt tolerance of plants by participating in JA responses. However, how BpSZA1 responds to salt stress through involvement in the JA signaling pathway requires further study.

Conclusion

In summary, we identified a novel C2H2 ZFP, BpSZA1, which plays a positive role in the plant salt stress response. We analyzed the function of BpSZA1 by generating OE and Ri in transgenic birch plants. We further showed that BpSZA1 can bind directly to the GAGA-motif in the promoter regions of various genes related to ROS scavenging to regulate their expression and maintain ROS homeostasis (Figure 8). Thus, our findings provide new insights into the function of C2H2 ZFPs in abiotic stress and further efforts to cultivate salt-tolerant birch varieties.

FIGURE 8
www.frontiersin.org

Figure 8. The model of BpSZA1 regulation network in response to Salt stress. Salt stress signals activate the expression of upstream genes of BpSZA1, which in turn activates BpSZA1. BpSZA1 binds to the GAGA-box in the promoter region, thereby activating the expression of downstream genes such as BpAPX1, BpAPX2, BpCAT, and Bp6PGDH to increase ROS scavenging. The result is tolerance to salt stress is increased in the plants.

Data Availability Statement

The original contributions presented in the study are publicly available. This data can be found here: NCBI, SRA, PRJNA808444.

Author Contributions

LL designed the research. XZ, QG, and LQ conducted the experiments and data analysis. XZ wrote the manuscript and performed the transcriptome data analysis. All authors read and approved the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 31470676).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.823547/full#supplementary-material

Footnotes

  1. ^ http://web.expasy.org/cgi-bin/protparam/
  2. ^ https://www.arabidopsis.org
  3. ^ https://phytozome.jgi.doe.gov/pz/portal.html
  4. ^ https://genomevolution.org/coge/GenomeInfo.pl?gid=35079
  5. ^ http://meme-suite.org/
  6. ^ http://ccb.jhu.edu/software/hisat2/index.shtml
  7. ^ http://tophat.cbcb.umd.edu/
  8. ^ https://ccb.jhu.edu/software/stringtie/index.shtml
  9. ^ https://www.biostat.wisc.edu/∼kendzior/EBSEQ/

References

Abe, H., Urao, T., Ito, T., Seki, M., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2003). Arabidopsis AtMYC2 (bHLH) and AtMYB2 (MYB) function as transcriptional activators in abscisic acid signaling. Plant Cell 15, 63–78. doi: 10.1105/tpc.006130

PubMed Abstract | CrossRef Full Text | Google Scholar

Audran, C., Borel, C., Frey, A., Sotta, B., Meyer, C., Simonneau, T., et al. (1998). Expression studies of the zeaxanthin epoxidase gene in Nicotiana plumbaginifolia. Plant Physiol. 118, 1021–1028. doi: 10.1104/pp.118.3.1021

PubMed Abstract | CrossRef Full Text | Google Scholar

Bates, L. S. (1973). Rapid determination of free proline for water-stress studies Summary. Plant Soil 39, 205–207. doi: 10.1007/bf00018060

CrossRef Full Text | Google Scholar

Baxter, A., Mittler, R., and Suzuki, N. (2014). ROS as key players in plant stress signalling. J. Exp. Bot. 65, 1229–1240. doi: 10.1093/jxb/ert375

PubMed Abstract | CrossRef Full Text | Google Scholar

Ben-Amor, M., Flores, B., Latché, A., Bouzayen, M., Pech, J. C., and Fomojaro, F. (1999). Inhibition of ethylene biosynthesis by antisense ACC oxidase RNA prevents chilling injury in Charentais cantaloupe melons. Plant Cell Environ. 22, 1579–1586. doi: 10.1046/j.1365-3040.1999.00509.x

CrossRef Full Text | Google Scholar

Berg, J. M., and Shi, Y. (1996). The galvanization of biology: a growing appreciation for the roles of zinc. Science 271, 1081–1085. doi: 10.1126/science.271.5252.1081

PubMed Abstract | CrossRef Full Text | Google Scholar

Chanwala, J., Satpati, S., Dixit, A., Parida, A., Giri, M. K., and Dey, N. (2020). Genome-wide identification and expression analysis of WRKY transcription factors in pearl millet (Pennisetum glaucum) under dehydration and salinity stress. BMC Genomics 21:231. doi: 10.1186/s12864-020-6622-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J., Yang, L., Yan, X., Liu, Y., Wang, R., Fan, T., et al. (2016). Zinc-Finger transcription factor ZAT6 positively regulates cadmium Tolerance through the glutathione-dependent pathway in Arabidopsis. Plant Physiol. 171, 707–719. doi: 10.1104/pp.15.01882

PubMed Abstract | CrossRef Full Text | Google Scholar

Ciftci-Yilmaz, S., Morsy, M. R., Song, L., Coutu, A., Krizek, B. A., Lewis, M. W., et al. (2007). The EAR-motif of the Cys2/His2-type zinc finger protein Zat7 plays a key role in the defense response of Arabidopsis to salinity stress. J. Biol. Chem. 282:9260. doi: 10.1074/jbc.M611093200

PubMed Abstract | CrossRef Full Text | Google Scholar

Corpas, F. J., Barroso, J. B., Sandalio, L. M., Distefano, S., Palma, J. M., Lupianez, J. A., et al. (1998). A dehydrogenase-mediated recycling system of NADPH in plant peroxisomes. Biochem. J. 330(Pt 2), 777–784. doi: 10.1042/bj3300777

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, H., Chen, Q., Dai, Y., Hu, W., Zhang, S., and Huang, X. (2021). Genome-wide Identification of PbrbHLH Family Genes, and Expression Analysis in Response to Drought and Cold Stresses in Pear (Pyrus Bretschneideri). BMC Plant Biol. 21:86. doi: 10.1186/s12870-021-02862-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Englbrecht, C., Schoof, H., and Böhm, S. (2004). Conservation, diversification and expansion of C2H2 zinc finger proteins in the Arabidopsis thaliana genome. BMC Genomics 5:39. doi: 10.1186/1471-2164-5-39

PubMed Abstract | CrossRef Full Text | Google Scholar

Eveland, A. L., Goldshmidt, A., Pautler, M., Morohashi, K., Liseron-Monfils, C., Lewis, M. W., et al. (2014). Regulatory modules controlling maize inflorescence architecture. Genome Res. 24, 431–443. doi: 10.1101/gr.166397.113

PubMed Abstract | CrossRef Full Text | Google Scholar

Farhangi-Abriz, S., and Ghassemi-Golezani, K. (2018). How can salicylic acid and jasmonic acid mitigate salt toxicity in soybean plants? Ecotoxicol. Environ. Saf. 147, 1010–1016. doi: 10.1016/j.ecoenv.2017.09.070

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernandez-Calvo, P., Chini, A., Fernandez-Barbero, G., Chico, J. M., Gimenez-Ibanez, S., Geerinck, J., et al. (2011). The Arabidopsis bHLH transcription factors MYC3 and MYC4 Are targets of JAZ repressors and act additively with MYC2 in the activation of jasmonate responses. Plant Cell 23, 701–715. doi: 10.1105/tpc.110.080788

PubMed Abstract | CrossRef Full Text | Google Scholar

Franco-Zorrilla, J. M., López-Vidriero, I., Carrasco, J. L., Godoy, M., Vera, P., and Solano, R. (2014). DNA-binding specificities of plant transcription factors and their potential to define target genes. Proc. Natl. Acad. Sci. U.S.A. 111, 2367–2372. doi: 10.1073/pnas.1316278111

PubMed Abstract | CrossRef Full Text | Google Scholar

Fryer, M. J., Oxborough, K., Mullineaux, P. M., and Baker, N. R. (2002). Imaging of photo-oxidative stress responses in leaves. J. Exp. Bot. 53, 1249–1254. doi: 10.1093/jxb/53.372.1249

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, A., He, K., Liu, D., Bai, S., Gu, X., Wei, L., et al. (2005). DATF: a database of Arabidopsis transcription factors. Bioinformatics 21, 2568–2569. doi: 10.1093/bioinformatics/bti334

PubMed Abstract | CrossRef Full Text | Google Scholar

Gupta, B., and Huang, B. (2014). Mechanism of salinity tolerance in plants: physiological, biochemical, and molecular characterization. Int. J. Genomics 2014:701596. doi: 10.1155/2014/701596

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasegawa, P. M., Bressan, R. A., Zhu, J. K., and Bohnert, H. J. (2000). Plant cellular and molecular responses to high salinity. Annu. Rev. Plant Biol. 51, 463–499. doi: 10.1146/annurev.arplant.51.1.463

PubMed Abstract | CrossRef Full Text | Google Scholar

He, F., Li, H. G., Wang, J. J., Su, Y., Wang, H. L., Feng, C. H., et al. (2019). PeSTZ1, a C2H2-type zinc finger transcription factor from Populus euphratica, enhances freezing tolerance through modulation of ROS scavenging by directly regulating PeAPX2. Plant Biotechnol. J. 17, 2169–2183. doi: 10.1111/pbi.13130

PubMed Abstract | CrossRef Full Text | Google Scholar

He, F., Niu, M. X., Feng, C. H., Li, H. G., Su, Y., Su, W. L., et al. (2020). PeSTZ1 confers salt stress tolerance by scavenging the accumulation of ROS through regulating the expression of PeZAT12 and PeAPX2 in Populus. Tree Physiol. 40, 1292–1311. doi: 10.1093/treephys/tpaa050

PubMed Abstract | CrossRef Full Text | Google Scholar

Heath, R. L., and Packer, L. (1968). Photoperoxidation in isolated chloroplasts. Arch. Biochem. Biophys. 125, 189–198. doi: 10.1016/0003-9861(68)90654-1

CrossRef Full Text | Google Scholar

Holsters, M., de Waele, D., Depicker, A., Messens, E., van Montagu, M., and Schell, J. (1978). Transfection and transformation of Agrobacterium tumefaciens. Mol. Gen. Genet. 163, 181–187. doi: 10.1007/BF00267408

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, D., Zheng, Q., Melchkart, T., Bekkaoui, Y., Konkin, D. J. F., Kagale, S., et al. (2020). Dominant inhibition of awn development by a putative zinc-finger transcriptional repressor expressed at the B1 locus in wheat. New Phytol. 225, 340–355. doi: 10.1111/nph.16154

PubMed Abstract | CrossRef Full Text | Google Scholar

Iida, K., Seki, M., Sakurai, T., Satou, M., Akiyama, K., Toyoda, T., et al. (2005). RARTF: database and tools for complete sets of Arabidopsis transcription factors. DNA Res. 12, 247–256. doi: 10.1093/dnares/dsi011

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, M., Xu, F., Peng, M., Huang, F., and Meng, F. (2016). Methyl jasmonate regulated diploid and tetraploid black locust (Robinia pseudoacacia L.) tolerance to salt stress. Acta Physiol. Plant. 38, 1–13.

Google Scholar

Karimi, M., Inzé, D., and Depicker, A. (2002). GATEWAY vectors for Agrobacterium-mediated plant transformation. Trends Plant Sci. 7, 193–195. doi: 10.1016/S1360-1385(02)02251-3

CrossRef Full Text | Google Scholar

Kavi Kishor, P. B., Hong, Z., Miao, G.-H., Hu, C.-A. A., and Verma, D. P. S. (1995). Overexpression of Δ1-pyrroline-5-carboxylate synthetase increases proline production and confers osmotolerance in transgenic plants. Plant Physiol. 108, 1387–1394. doi: 10.1104/pp.108.4.1387

PubMed Abstract | CrossRef Full Text | Google Scholar

Kazuhiko, M. (2007). Floral sex allocation at individual and branch levels in Betula platyphylla var. japonica (Betulaceae), a tall, wind-pollinated monoecious tree species. Am. J. Bot. 94, 1450–1458. doi: 10.3732/ajb.94.9.1450

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiełbowicz-Matuk, A. (2012). Involvement of plant C2H2-type zinc finger transcription factors in stress responses. Plant Sci. 18, 78–85. doi: 10.1016/j.plantsci.2011.11.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Kodaira, K. S., Qin, F., Tran, L. S. P., Maruyama, K., Kidokoro, S., Fujita, Y., et al. (2011). Arabidopsis Cys2/His2 zinc-finger proteins AZF1 and AZF2 negatively regulate abscisic acid-repressive and auxin-inducible genes under abiotic stress conditions. Plant Physiol. 157, 742–756. doi: 10.1104/pp.111.182683

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuang, J. F., Wu, C. J., Guo, Y. F., Walther, D., Shan, W., Chen, J., et al. (2021). Deciphering transcriptional regulators of banana fruit ripening by regulatory network analysis. Plant Biotechnol. J. 19, 477–489. doi: 10.1111/pbi.13477

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C. Y., Puhakainen, T., Welling, A., Viherae-Aarnio, A., Ernstsen, A., Junttila, O., et al. (2002). Cold acclimation in silver birch (Betula pendula): development of freezing tolerance in different tissues and climatic ecotypes. Physiol. Plant. 116, 478–488. doi: 10.1034/j.1399-3054.2002.1160406.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, P., Peng, Z., Xu, P., Tang, G., Ma, C., Zhu, J., et al. (2021). Genome-wide identification of NAC transcription factors and their functional prediction of abiotic stress response in peanut. Front. Genet. 12:630292. doi: 10.3389/fgene.2021.630292

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, W., Lin, Y. C., Li, Q., Shi, R., Lin, C. Y., Chen, H., et al. (2014). A robust chromatin immunoprecipitation protocol for studying transcription factor-DNA interactions and histone modifications in wood-forming tissue. Nat. Protoc. 9, 2180–2193. doi: 10.1038/nprot.2014.146

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Chu, Z., Luo, J., Zhou, Y., Cai, Y., Lu, Y., et al. (2018). The C2H2 zinc-finger protein Sl ZF 3 regulates AsA synthesis and salt tolerance by interacting with CSN5B. Plant Biotechnol. J. 16, 1201–1213. doi: 10.1111/pbi.12863

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y. D., Shi, Y., Zhu, N., Zhong, S. L., Bouzayen, M., and Li, Z. G. (2020). SlGRAS4 mediates a novel regulatory pathway promoting chilling tolerance in tomato. Plant Biotechnol. J. 18, 1620–1633. doi: 10.1111/pbi.13328

PubMed Abstract | CrossRef Full Text | Google Scholar

Marino, D., Dunand, C., Puppo, A., and Pauly, N. (2012). A burst of plant NADPH oxidases. Trends Plant Sci. 17, 9–15. doi: 10.1016/j.tplants.2011.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, J., Mclachlan, A. D., and Klug, A. (1985). Repetitive zinc-binding domains in the protein transcription factor IIIA from Xenopus oocytes. EMBO J. 4, 1609–1614. doi: 10.1002/j.1460-2075.1985.tb03825.x

CrossRef Full Text | Google Scholar

Minotti, P. L., Halseth, D. E., and Sieczka, J. B. (1994). Field chlorophyll measurements to assess the nitrogen status of potato varieties. Hortscience 29, 1497–1500. doi: 10.21273/hortsci.29.12.1497

CrossRef Full Text | Google Scholar

Mittler, R., Kim, Y., Song, L., Coutu, J., Coutu, A., Ciftci-Yilmaz, S., et al. (2006). Gain- and loss-of-function mutations in Zat10 enhance the tolerance of plants to abiotic stress. FEBS Lett. 580, 6537–6542. doi: 10.1016/j.febslet.2006.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Mittler, R., Vanderauwera, S., Gollery, M., and Breusegem, F. V. (2004). Reactive oxygen gene network of plants. Trends Plant Sci. 9:490. doi: 10.1016/j.tplants.2004.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Munns, R., and Tester, M. (2008). Mechanisms of salinity tolerance. Annu. Rev. Plant Biol. 59, 651–681. doi: 10.1146/annurev.arplant.59.032607.092911

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakano, Y., and Asada, K. (1981). Hydrogen peroxide is scavenged by ascorbate-specific peroxidase in spinach chloroplasts. Plant Cell Physiol. 22, 867–880.

Google Scholar

Nelson, J. D., Denisenko, O., and Bomsztyk, K. (2006). Protocol for the fast chromatin immunoprecipitation (ChIP) method. Nat. Protoc. 1, 179–185. doi: 10.1038/nprot.2006.27

PubMed Abstract | CrossRef Full Text | Google Scholar

Pauwels, L., Morreel, K., De Witte, E., Lammertyn, F., Van Montagu, M., Boerjan, W., et al. (2008). Mapping methyl jasmonate-mediated transcriptional reprogramming of metabolism and cell cycle progression in cultured Arabidopsis cells. Proc. Natl. Acad. Sci. U.S.A. 105, 1380–1385. doi: 10.1073/pnas.0711203105

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedone, P. V., Ghirlando, R., Clore, G. M., Gronenborn, A. M., Felsenfeld, G., and Omichinski, J. G. (1996). The single Cys(2)-His(2) zinc finger domain of the GAGA protein flanked by basic residues is sufficient for high-affinity specific DNA binding. Proc. Natl. Acad. Sci. U.S.A. 93, 2822–2826. doi: 10.1073/pnas.93.7.2822

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, Z., Guo, J., Zhu, A., Zhang, L., and Zhang, M. (2014). Exogenous jasmonic acid can enhance tolerance of wheat seedlings to salt stress. Ecotoxicol. Environ. Saf. 104, 202–208. doi: 10.1016/j.ecoenv.2014.03.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Regier, N., and Frey, B. (2010). Experimental comparison of relative RT-qPCR quantification approaches for gene expression studies in poplar. BMC Mol. Biol. 11:57. doi: 10.1186/1471-2199-11-57

PubMed Abstract | CrossRef Full Text | Google Scholar

Rehman, A., Wang, N., Peng, Z., He, S., Zhao, Z., Gao, Q., et al. (2021). Identification of C2H2 subfamily ZAT genes in Gossypium species reveals GhZAT34 and GhZAT79 enhanced salt tolerance in Arabidopsis and cotton. Int. J. Biol. Macromol. 184, 967–980. doi: 10.1016/j.ijbiomac.2021.06.166

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, Z., Zhang, D., Cao, L., Zhang, W., Zheng, H., Liu, Z., et al. (2020). Functions and regulatory framework of ZmNST3 in maize under lodging and drought stress. Plant Cell Environ. 43, 2272–2286. doi: 10.1111/pce.13829

PubMed Abstract | CrossRef Full Text | Google Scholar

Riechmann, J. L., Heard, J., Martin, G., Reuber, L., Jiang, C. Z., Keddie, J., et al. (2000). Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science 290, 2105–2110. doi: 10.1126/science.290.5499.2105

PubMed Abstract | CrossRef Full Text | Google Scholar

Rizhsky, L., Davletova, S., Liang, H., and Mittler, R. (2004). The zinc finger protein Zat12 is required for cytosolic ascorbate peroxidase 1 expression during oxidative stress in Arabidopsis. J. Biol. Chem. 279, 11736–11743. doi: 10.1074/jbc.m313350200

PubMed Abstract | CrossRef Full Text | Google Scholar

Ruan, J. J., Zhou, Y. X., Zhou, M. L., Yan, J., Khurshid, M., Weng, W. F., et al. (2019). Jasmonic acid signaling pathway in plants. Int. J. Mol. Sci. 20:2479. doi: 10.3390/ijms20102479

PubMed Abstract | CrossRef Full Text | Google Scholar

Seki, M., Ayako, K., Kazuko, Y. S., and Kazuo, S. (2003). Molecular responses to drought, salinity and frost: common and different paths for plant protection. Curr. Opin. Biotechnol. 14, 194–199. doi: 10.1016/S0958-1669(03)00030-2

CrossRef Full Text | Google Scholar

Shrivastava, P., and Kumar, R. (2015). Soil salinity: a serious environmental issue and plant growth promoting bacteria as one of the tools for its alleviation. Saudi J. Biol. Sci. 22, 123–131. doi: 10.1016/j.sjbs.2014.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Shyu, C., Figueroa, P., DePew, C. L., Cooke, T. F., Sheard, L. B., Moreno, J. E., et al. (2012). JAZ8 lacks a canonical degron and has an EAR motif that mediates transcriptional repression of jasmonate responses in Arabidopsis. Plant Cell 24, 536–550. doi: 10.1105/tpc.111.093005

PubMed Abstract | CrossRef Full Text | Google Scholar

Sparkes, I. A., Runions, J., Kearns, A., and Hawes, C. (2006). Rapid, transient expression of fluorescent fusion proteins in tobacco plants and generation of stably transformed plants. Nat. Protoc. 1, 2019–2025. doi: 10.1038/nprot.2006.286

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, N., Koussevitzky, S., Mittler, R. O. N., and Miller, G. A. D. (2012). ROS and redox signalling in the response of plants to abiotic stress. Plant Cell Environ. 35, 259–270. doi: 10.1111/j.1365-3040.2011.02336.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Takatsuji, H. (1998). Zinc-finger transcription factors in plants. Cell. Mol. Life Sci. 54, 582–596. doi: 10.1007/s000180050186

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, H., Bi, H., Liu, B., Lou, S., Song, Y., Tong, S., et al. (2021). WRKY33 interacts with WRKY12 protein to up-regulate RAP2.2 during submergence induced hypoxia response in Arabidopsis thaliana. New Phytol. 229, 106–125. doi: 10.1111/nph.17020

PubMed Abstract | CrossRef Full Text | Google Scholar

Tavallali, V., and Karimi, S. (2019). Methyl jasmonate enhances salt tolerance of almond rootstocks by regulating endogenous phytohormones, antioxidant activity and gas-exchange. J. Plant Physiol. 23, 98–105. doi: 10.1016/j.jplph.2019.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Teng, K., Tan, P., Guo, W., Yue, Y., Fan, X., and Wu, J. (2018). Heterologous expression of a novel Zoysia japonica C2H2 zinc finger gene. ZjZFN1, improved salt tolerance in Arabidopsis. Front. Plant Sci. 9:1159. doi: 10.3389/fpls.2018.01159

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J. Q., Zhang, X., and Li, L. (2018). Bioinformatic and Expression Analysis of HD-Zip Family Gene in Betula platyphylla. Bull. Bot. Res. 38, 931–938.

Google Scholar

Wang, K., Ding, Y., Cai, C., Chen, Z., and Zhu, C. (2019). The role of C2H2 zinc finger proteins in plant responses to abiotic stresses. Physiol. Plant. 165, 690–700. doi: 10.1111/ppl.12728

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W., Wang, X., Wang, Y., Zhou, G., Wang, C., Hussain, S., et al. (2020). SlEAD1, an EAR motif-containing ABA down-regulated novel transcription repressor regulates ABA response in tomato. GM Crops Food 11:275. doi: 10.1080/21645698.2020.1790287

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y. H., Qin, L. L., and Li, L. (2016). Identification and expression analysis of genes encoding zinc finger proteins in Betula platyphylla. Bull. Bot. Res. 36, 395–400.

Google Scholar

Witzel, K., Weidner, A., Surabhi, G. K., Varshney, R. K., Kunze, G., Buck-Sorlin, G. H., et al. (2010). Comparative analysis of the grain proteome fraction in barley genotypes with contrasting salinity tolerance during germination. Plant Cell Environ. 33, 211–222. doi: 10.1111/j.1365-3040.2009.02071.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, H., Tang, J., Li, Y., Wang, W., Li, X., Jin, L., et al. (2009). STAMENLESS 1, encoding a single C2H2 zinc finger protein, regulates floral organ identity in rice. Plant J. 59, 789–801. doi: 10.1111/j.1365-313X.2009.03913.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, M., Sun, J., Gong, D., and Kong, Y. (2019). The roles of Arabidopsis C1-2i subclass of C2H2-type zinc-finger transcription factors. Genes 10:653. doi: 10.3390/genes10090653

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Ahammed, G. J., Wan, C., Liu, H., Chen, R., and Zhou, Y. (2019). Comprehensive analysis of TIFY transcription factors and their expression profiles under jasmonic acid and abiotic stresses in watermelon. Int. J. Genomics 2019:6813086. doi: 10.1155/2019/6813086

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., and Guo, Y. (2018). Unraveling salt stress signaling in plants. J. Integr. Plant Biol. 60, 796–804. doi: 10.1111/jipb.12689

PubMed Abstract | CrossRef Full Text | Google Scholar

Young, M. D., Wakefield, M. J., Smyth, G. K., and Oshlack, A. (2010). Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 11:R14. doi: 10.1186/gb-2010-11-2-r14

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Liu, H., Sun, C., Ma, Q., Bu, Hs., Chong, K., et al. (2018). A C2H2 zinc-finger protein OsZFP213 interacts with OsMAPK3 to enhance salt tolerance in rice. J. Plant Physiol. 229, 100–110. doi: 10.1016/j.jplph.2018.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Wang, Y., Chang, L., Zhang, T., An, J., Liu, Y., et al. (2016). MsZEP, a novel zeaxanthin epoxidase gene from alfalfa (Medicago sativa), confers drought and salt tolerance in transgenic tobacco. Plant Cell Rep. 35, 439–453. doi: 10.1007/s00299-015-1895-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Q., Hu, R. S., Liu, D., Liu, X., Wang, J., Xiang, X. H., et al. (2020). The AP2 transcription factor NtERF172 confers drought resistance by modifying NtCAT. Plant Biotechnol. J. 18, 2444–2455. doi: 10.1111/pbi.13419

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhuang, H., Wang, H. L., Zhang, T., Zeng, X. Q., Chen, H., Wang, Z. W., et al. (2020). NONSTOP GLUMES1 encodes a C2H2 zinc finger protein that regulates spikelet development in rice. Plant Cell 32, 392–413. doi: 10.1105/tpc.19.00682

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, J. K. (2002). Salt and drought stress signal transduction in plants. Annu. Rev. Plant Biol. 53, 247–273. doi: 10.1146/annurev.arplant.53.091401.143329

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, J. K. (2016). Abiotic stress signaling and responses in plants. Cell 167, 313–324. doi: 10.1016/j.cell.2016.08.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Betula platyphylla, BpSZA1, salt stress, ChIP, ROS scavenging

Citation: Zhang X, Guo Q, Qin L and Li L (2022) A Cys2His2 Zinc Finger Transcription Factor BpSZA1 Positively Modulates Salt Stress in Betula platyphylla. Front. Plant Sci. 13:823547. doi: 10.3389/fpls.2022.823547

Received: 27 November 2021; Accepted: 29 March 2022;
Published: 25 May 2022.

Edited by:

Loredana F. Ciarmiello, University of Campania Luigi Vanvitelli, Italy

Reviewed by:

Necla Pehlivan, Recep Tayyip Erdoğan University, Turkey
Mingjun Li, Northwest A&F University, China

Copyright © 2022 Zhang, Guo, Qin and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Li Li, bGlsaUBuZWZ1LmVkdS5jbg==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.