The Roles of GmERF135 in Improving Salt Tolerance and Decreasing ABA Sensitivity in Soybean.

Abscisic acid (ABA) mediates various abiotic stress responses, and ethylene responsive factors (ERFs) play vital role in resisting stresses, but the interaction of these molecular mechanisms remains elusive. In this study, we identified an ABA-induced soybean ERF gene GmERF135 that was highly up-regulated by ethylene (ET), drought, salt, and low temperature treatments. Subcellular localization assay showed that the GmERF135 protein was targeted to the nucleus. Promoter cis-acting elements analysis suggested that numerous potential stress responsive cis-elements were distributed in the promoter region of GmERF135, including ABA-, light-, ET-, gibberellin (GA)-, and methyl jasmonate (MeJA)-responsive elements. Overexpression of GmERF135 in Arabidopsis enhanced tolerance to drought and salt conditions. In addition, GmERF135 promoted the growth of transgenic hairy roots under salt and exogenous ABA conditions. These results suggest that soybean GmERF135 may participate in both ABA and ET signaling pathways to regulate the responses to multiple stresses.

Introduction

As a vital hormone in plants, abscisic acid (ABA) is essential to a myriad of aspects of plant growth and developmental processes, including plant embryo development, seed maturation, fruit maturity, and stomatal movement (Finkelstein et al., 2002). ABA signal transduction has been studied for many years (Pandey et al., 2006, 2009; Shen et al., 2006; Ma et al., 2009; Park et al., 2009; Yang and Guo, 2018), and the widely accepted molecular mechanism is that the pyrabactin resistant/PYR-like/regulatory component of ABA receptor (PYR/PYL/RCAR) acts as an ABA receptor which can bind to ABA, and then binds to and inhibits protein phosphatases type 2C (PP2Cs) (Santiago et al., 2009). In addition, the activity of SNF1-related protein kinases 2 (SnRK2s) is enhanced and can phosphorylate ABRE binding factors (AREB/ABFs) to induce physiological and biochemical changes in the process of ABA response (Hubbard et al., 2010; Yoshida et al., 2015). A recent study showed that ABFs can directly bind to the promoters of group A PP2C genes, and rapidly induce their expression on exogenous ABA treatments (Wang et al., 2018). However, the downstream molecular mechanism is yet not clearly understood.

Many transcription factors, such as MYC/MYB, bZIP/ABRE, AP2/ERF, and NAC, are regulated by ABA. The AP2/ERF family, the largest plant transcription factor family (Okamuro et al., 1997), can be divided into three subfamilies: the AP2, ERF, and RAV subfamilies. Among them, the ERF subfamily can specifically bind to GCC-box and/or the dehydration-responsive element/C-repeat (DRE/CRT) cis-acting elements (Allen et al., 1998).

The AP2/ERF family is involved in responses to various abiotic stresses and exogenous hormones (Xu et al., 2008, 2011; Klay et al., 2018). For example, overexpressing TaERF1 enhances tolerance of physiological and environmental stresses, such as salt, drought, low temperature, exogenous ABA, ethylene (ET), salicylic acid (SA), and disease (Xu et al., 2007). Transgenic rice plants expressing OsERF922 displayed higher susceptibility to Magnaporthe oryzae and NaCl compared to the wild type (WT), while the knockout mutant and RNAi lines enhanced resistance to these stresses (Liu et al., 2012). A recent study described how ABA modulates the expression level of ERF family, showing that the expression levels of ABA-responsive genes such as RD22, LEA3, and PODs were up-regulated after overexpressing OsERF101 in rice, which enhanced its susceptibility to ABA (Jin et al., 2018).

Transcriptomic analysis of grapevine organs treated with or without ABA showed that ERF members were involved in differently expressed genes (DEGs), and the ERF subfamily had the most significant change compared to other transcription factors (Rattanakon et al., 2016). Although the ABA signal pathway has been extensively studied, it is yet unclear how regulation of the expression of downstream genes via ERF subfamily could enhance abiotic stress tolerances in soybean. To investigate whether the ERF subfamily is modulated by the ABA signal in soybean, we studied the response to exogenous ABA and identified the functions of GmERF135 in transgenic Arabidopsis and soybean hairy roots.

Materials and Methods

Plant Materials and Treatment

Soybean cultivar “Tiefeng 8” was sown in pots containing vermiculite and grown at 25°C for 14 days and then treated with various abiotic stresses. For the various treatments, soybean plants were exposed in the air for rapid drought, dipped into 200 mM NaCl for salt stress, placed in 42°C/4°C chamber for high/low temperature stress, respectively. They were also placed in an airtight container filled with ET, or dipped in 100 μM ABA, 50 μM salicylic acid (SA), or 50 μM jasmonate (JA) for exogenous hormone stresses. Leaves of these plants were collected after 0, 0.5, 1, 2, 5, 12, and 24 h treatment and then immediately stored at −80°C for RNA extraction.

Gene Structure and Protein Domain of ERFs in Soybean

The whole genome data of the candidate ERF genes was obtained from JGI Glyma1.0 annotation[1] (Goodstein et al., 2012). Gene structure was analyzed by submitting the CDS of the candidate genes and whole genomic DNA sequences to the Gene Structure Display Server (GSDS) website[2]. Protein Fold Recognition Server (PHYRE2)[3] was used for analysis of structural homology modeling of these genes, and DOG 2.0 was used to draw the protein domains.

Quantitative Real-Time PCR (qRT-PCR)

Total RNA of the soybean plants was extracted using the TRIzol reagent (Invitrogen, Carlsbad, CA, United States). The specific primer pairs for the 16 genes were designed by Primer Premier 5.0 according to the cDNAs. qRT-PCR was conducted using the SYBR Premix Ex TaqTM kit (TAKARA, Kusatsu, Japan) and the ABI Prism 7500 real-time PCR system (Thermo Fisher Scientific, Waltham, MA, United States). The 2–ΔΔCT method was used to conduct qRT-PCR analysis (Le et al., 2011). Soybean Actin (U60506) was used as an internal control and sequence of specific primers were shown in Supplementary Table S1.

Isolation and Promoter Analysis of GmERF135

The GmERF135 gene was amplified by PCR and the primers for cloning the GmERF135 were 5′-AATCATTATGTGTGGCGGTGCC-3′ and 5′-TATTCCTCGC TAATCGAAACTCCAGAG-3′. The PCR product was then cloned into the pEASY-T1 vector (TransGen, Beijing, China). 1,886 bp promotor region of GmERF135 was cloned and submitted to PLACE[4] and PlantCARE[5] databases to analyze the putative cis-acting elements in the promoter region (Lescot et al., 2002).

Subcellular Localization Analysis

The cDNA of GmERF135 were augmented with PCR, connected to the N-terminus of humanized green fluorescent protein (hGFP) reporter gene under the control of the double Cauliflower Mosaic Virus (2 × CaMV) 35S promoter, and a recombinant plasmid was obtained (Scott et al., 1999). The recombinant plasmid was introduced into onion epidermal cells, while the onion epidermal cells with hGFP vector acted as the control. Fluorescence microscopy was used to identify hGFP expression (Xu et al., 2007; He et al., 2016).

Generation of Transgenic Arabidopsis and Stress Treatments

AT5G47230 is the orthologue of GmERF135 in Arabidopsis, which was named ERF2 and used to investigate function in responses to various stresses. The seeds of two mutants, erf2-1 and erf2-2 (SALK_126889, SALK_076967) were mutated via T-DNA insertion.

The cDNA of GmERF135 was obtained by using the specific primer pairs: 5′-TGATTACGCCAAGCTTATGTGTG GCGGTGCC-3′, 5′-CCGGGGATCCTCTAGAATCGAAACTC CAGAG-3′, and then were cloned into pBI121 under the control of the CaMV 35S promoter. The recombinant plasmid 35S::GmERF135 was sequenced and transformed into wild-type (WT) and the two mutants Arabidopsis lines using the vacuum infiltration method (Bechtold and Pelletier, 1998; Liu et al., 2013). T3 seeds of transgenic lines were selected for further analysis.

For root length assay, 30 seeds of each line were sown on 1/2-strength MS growth medium with or without 6% PEG, 75 mM NaCl, 1-aminocyclopropane-1-carboxylic acid (ACC), or in dark treatment for growth, respectively. At least 30 seedlings per line were randomly selected to measure root length. The cotyledon pieces of each line were recorded every 12 h. Each treatment contained three independent replicates.

Soybean Hairy Root Induction and Stress Treatments

The Superroot of Lotus corniculatus and Cucumopine-type Agrobacterium rhizogene strain K599 with pGFPGUSPlus were provided by Professor Tian-Fu Han (CAAS, China). Seedling growth, rooting, hairy root induction, and hairy root transformation were performed as described by Chen et al. (2004, 2014)

GFP positive (GFP+) hairy roots were induced by K599 carrying the pGFPGUSPlus-GmERF135 binary vector. These hairy roots were cultured on 1/2 MS medium supplemented with 0, 50, 85, 120, or 150 mM NaCl for salt treatment, and 50, 100, or 150 μM ABA for hormone treatment, respectively. They were then incubated at 24°C under a 16/8 h light/dark cycle condition for 2 weeks. After 24 h incubation at 105°C, the increase in dry weight (30 roots per unit) was measured and recorded. Each treatment contained three independent replicates.

Results

Molecular Characterization of Soybean Targeted ERF Genes

In a previous study, 160 non-redundant soybean ERFs were identified using the Pfam database and SMART program, and these soybean ERFs were clustered into eight groups (Supplementary Table S2; Zhao et al., unpublished). To comprehensively understand the responses of soybean ERFs to ABA, 16 ERF genes were selected for further investigation according to the phylogenetic tree. Gene structure analysis showed that six genes had no introns, including GmERF106, GmERF132, GmERF135, GmERF41, GmERF49, and GmERF84 (Figure 1). The remaining 10 ERF genes contained one intron, which was distributed in each cDNA region except for GmERF103, GmERF111, and GmERF15. The protein domains of each soybean ERF gene were drawn by DOG2.0 (Figure 2). All the proteins contained a conservative AP2/ERF domain, which was distributed in different positions of each protein sequence. The AP2/ERF domain of the proteins in each group displayed similar locations such as Group I, III, IV, V, VII, and VIII.

FIGURE 1

Intron-exon structure of candidate soybean ERF genes. Genome data for candidate ERF genes was obtained from JGI Glyma1.0 annotation. Intron-exon structure was produced using the GSDS online tool by submitting CDSs and genomic sequences of candidate genes. Yellow boxes represent exons, blue boxes represent introns, and black lines represent untranslated regions (UTRs).

FIGURE 2

Protein domains of candidate soybean ERFs. Conserved motifs were analyzed using MEME, and DOG 2.0 was used to draw the domains. The conserved AP2/ERF domain is indicated by the green boxes. The number under green boxes shows the location of each domain.

An expression pattern map of soybean ERFs was drawn based on the gene-chip data downloaded from the soybean genome database[6]. As shown in Supplementary Figure S1A, the soybean ERFs were expressed at different levels in various tissues and organs. Among them, GmERF111 and GmERF135 showed high expression levels in almost all tissues and organs. GmERF75 showed high expression levels in root and nodule, and GmERF49 was highly expressed in nodule and 10 DAF pod shell. GmERF15 mainly expressed in young leaf. All the candidate genes had different expression patterns.

Isolation and Subcellular Localization Characterization of GmERF135

Abscisic acid plays an important role in plant growth and development which is closely related to quality and yield in plants. To investigate how ABA affected the expression pattern of the 16 soybean ERF genes, qRT-PCR was conducted. Almost all the ERFs were up-regulated by ABA, except for GmERF49 (Figure 3). The ERF gene GmERF135 was the highest expressed gene after exogenous ABA treatment except for GmERF75 which has been studied (Zhao et al., unpublished). The transcript level of GmERF135 was activated and showed a 17-fold increase within 2 h after treatment and then declined to normal. Considering the high expression level and ABA-responsive increase in plants, GmERF135 was selected for further study.

FIGURE 3

Relative expression levels of ERF genes after ABA treatment. Total RNA was extracted from soybean seedlings after ABA treatment for 0, 0.5, 1, 2, 5, 12, and 24 h and used for qRT-PCR. The data were shown as the means ± SD obtaining from three biological replicates.

The subcellular localization assay was conducted to provide clues to understand intrinsic characteristics of cell activities. hGFP reporter gene fused to C-terminus of GmERF135, which could fluoresce when laser irradiated. The GmERF135::hGFP fused protein only fluoresced in nucleus, while the control fluoresced throughout the entire cell (Figure 4). This result showed that GmERF135 functions in nucleus.

FIGURE 4

Subcellular localization of GmERF135. GmERF135::hGFP and control vector were introduced to onion epidermal cells and GmERF135::hGFP fusion protein and the control hGFP protein were expressed in the cells. Photographs were taken by confocal microscope (Leica) in the dark field for green fluorescence (A and D), in bright light for morphology of the cells (B and E) and in combination (C and F).

GmERF135 Promoter Region Comprises Diverse Stress-Responsive Elements

The promoter region is an important part of the gene which could regulate gene expression and control gene action. To investigate the potential regulation mechanism, the promoter region of GmERF135, which has 1886 bp length upstream of the start codon, was isolated. The PLACE and PlantCARE databases were used to analyze the putative regulatory elements in the promoter region. Several regulatory elements were identified to be involved in responses to abiotic stresses and plant hormones (Table 1).

TABLE 1

Analysis of putative cis-acting elements of GmERF135 promoter.

Cis-acting element	Position	Strand	Matrix score	Sequences	Function
ABRE	747	+	6	CACGTG	ABA responsive element
ABRE	1221	−	5	ACGTG	ABA responsive element
ABRE	748	+	5	ACGTG	ABA responsive element
Box 4	175	+	6	ATTAAT	Light responsive element
CGTCA-motif	1022	+	5	CGTCA	MeJA-responsive element
G-Box	747	+	6	CACGTG	Light responsive element
G-Box	1221	+	6	CACGTT	Light responsive element
MBS	1295	+	6	CAACTG	ABA and stress responsive element
TCT-motif	201	+	6	TCTTAC	light responsive element
TGACG-motif	1022	−	5	TGACG	MeJA responsive element
Core of MYBST1	57	−	5	GGATA	ABA and stress responsive element
Core of MYBST1	532	−	5	GGATA	ABA and stress responsive element
Core of MYBST1	1725	+	5	GGATA	ABA and stress responsive element
Core of MYBST1	1818	−	5	GGATA	ABA and stress responsive element
DPBF binding site	108	+	7	ACACNNG	ABA responsive element
DPBF binding site	134	+	7	ACACNNG	ABA responsive element
GARE	1,207	+	7	TAACAAR	Gibberellin responsive element
GARE	654	+	7	TAACAAR	Gibberellin responsive element
ERE	151	+	8	ATTTCAAA	Ethylene responsive element

The GmERF135 promoter region contains many ABA and stress responsive elements, including ABA-responsive elements (ABREs, 3 hits), MYBST1 core sequences (4 hits), MYB binding sites (MBS, 1 hit), and DPBF binding sites (2 hits) (Table 1). Except for ABREs, three hormone-responsive elements were predicted, including an ethylene- responsive element (ERE, 1 hit), MeJA-responsive element TGACG-motif (1 hit) and CGTCA-motif (1 hit), and GA-responsive element (GARE, 2 hits). In addition, TCT-motif (1 hit), G-box (2 hits), and Box 4 (1 hit) were found, which could act as light-responsive elements in the promoter region of GmERF135. These elements suggested that GmERF135 may be involved in responses to multiple abiotic stresses and exogenous hormones.

The Roles of GmERF135 in Responding to Multiple Stimuli

Quantitative Real-Time PCR was conducted to investigate the expression level of GmERF135 under abiotic stresses, such as drought, salt, high/low temperature, and exogenous hormones including ET, SA, and JA. GmERF135 was induced by almost all the abiotic stresses and exogenous hormones. Remarkably, the expression level of the gene was extremely up-regulated by drought, low temperature, and ET treatment. The transcription of GmERF135 peaked at 5 h under drought or low-temperature treatment which had 18-fold / 22-fold increases, respectively (Figure 5). The peak of GmERF135 transcription appeared at 12 h after treatment with ET, which showed a 28-fold increase. In contrast, GmERF135 rapidly reached the maximum transcript level within 0.5 h under the NaCl and JA treatments, and then declined to baseline after 1 h. Under the SA treatment, GmERF135 was induced and reached the peak within 2 h (about 7-fold) and then declined to normal level in 24 h. These results showed that GmERF135 may be involved in responses to multiple stimuli.

FIGURE 5

Expression patterns of GmERF135 under multiple stimuli. Two-week soybean seedlings were given various abiotic stresses and exogenous hormones for 0, 0.5, 1, 2, 5, 12, and 24 h, which were used to extract RNA and obtained cDNA. The transcript levels of GmERF135 after multiple stresses were quantified by qRT-PCR, including drought, NaCl, high/low temperature, ET (ethylene), JA (jasmonate), and SA (salicylic acid). Data were shown referring to three biological replicates.

GmERF135 Rescued Two erf2 Mutants and Affected Growth of the Root and Cotyledon

Total RNA was extracted for semi-quantitative PCR from hypocotyls, root, stem, and leaf tissues of normal growth soybean seedlings. Actin primers were used as a parallel reaction to normalize the added template amounts. GmERF135 was predominantly expressed in the leaf, less in the hypocotyl and stem, and very little in the root (Supplementary Figure S1B).

To assess whether GmERF135 could rescue the two mutants, the seeds of transgenic GmERF135::erf2 lines, WT, and two erf2 mutants were sown on the 1/2 MS medium for growth and the plant growth rate was recorded each 12 h. Six day later, the phenotypes were imaged (Figures 6A,B).

FIGURE 6

GmERF135 rescued the speed of cotyledon growth in the two mutants. (A) Phenotype of two erf2 mutants and WT Arabidopsis grown on the MS solid medium. (B) Phenotype of two GmERF135::erf2 lines and WT Arabidopsis grown on 1/2 MS solid medium. (C) Days needed for growing the cotyledon. Statistic data were recorded every 12 h. Three independent experiments (20 plants each) were performed. Error bars indicate SD of three biological replicates (n = 3). (D) Average root lengths of WT, erf2, and GmERF135::erf2 after grown on 1/2 MS medium for 6 days. Data were shown as the means ± SD of three biological replicates (at least 30 individuals per line).

Growth rate assay showed GmERF135::erf2 lines displayed faster growth rate compared to the mutants (Figure 6C). WT plants required 2 d of growth to sprout the two pieces of the cotyledon while the two mutants required 4 days. After introduction of GmERF135 in erf2 lines, 3 days was required for the same process. The two mutants needed 6 days to reach 4-leaf stage, while the transgenic GmERF135::erf2 lines had similar growth rate with WT, which only needed 5 d. Root length assay showed there is an approximately 25% decrease in erf2 mutants compared to WT lines (Figure 6D). After overexpressing GmERF135 in the erf2 mutants, the lines had similar root lengths to WT and the phenotype of which was rescued. These results showed that overexpression of GmERF135 in the two erf2 mutants partly rescued two erf2 mutants and affected growth of root and cotyledon in Arabidopsis.

GmERF135 Promotes Plant Growth Under Drought and Salt Stresses in Arabidopsis

To determine whether GmERF135 confers abiotic stress tolerance to Arabidopsis plants, 3-day-old WT, erf2 mutants, and GmERF135 overexpression seedlings were transferred to 1/2 MS medium with or without 6% PEG, 75 mM NaCl, ACC, and dark environment grown for 3 days. The growth rate and the root length of 30 seedlings in each line were recorded (Figures 7B,C). Statistics showed that the average root length of erf2 was shorter than the WT and transgenic Arabidopsis plants under all conditions except for ACC treatment (Figures 7A,B). The transgenic plants displayed larger cotyledon than the WT and mutants under NaCl and PEG treatments, which suggested that GmERF135 enhanced the growth rate of Arabidopsis plants. Overexpression of GmERF135 in Arabidopsis resulted in a longer root compared to the WT and two erf2 mutants under 6% PEG, 75 mM NaCl, and dark conditions. No detectable difference was observed between the transgenic lines and WT under ACC conditions.

FIGURE 7

GmERF135 transgenic Arabidopsis lines improved cotyledon growth under the NaCl and PEG conditions. (A) Phenotype of GmERF135 transgenic Arabidopsis lines under various stresses. Three-day-old GmERF135 transgenic Arabidopsis lines were transferred to 1/2 MS medium with 6% PEG, 75 mM NaCl, ACC, or grown in darkness. (B) Average root lengths of WT, erf2, and 35S::GmERF135 lines. Three-day-old seedlings were treated with various stresses for 3 days and the data were shown as the means ± SD of three biological replicates (at least 30 individuals per line). (C) Days needed for growing cotyledon. Statistic data were recorded every 12 h. Three biological replicates were performed.

To further understood the response to various stresses in Arabidopsis, some marker genes, such as NCED3 (Ruggiero et al., 2004), RD29A (Yamaguchi-Shinozaki and Shinozaki, 1994), COR15A (Xu et al., 2011), DREB2A (Liu et al., 1998), ABA1, ABA2 (He et al., 2016), ACO4, and ACS2 (Kim et al., 2013) were selected for qRT-PCR (Supplementary Figure S2).

GmERF135 Improves Tolerance of Salt and ABA in Soybean

To investigate the roles of GmERF135 in soybean, the pGUS-GmERF135 expression vector was constructed and transformed into soybean hairy roots, which were studied on MS basal medium supplemented with different concentrations of ABA and NaCl. GmERF135 transgenic soybean hairy roots experienced greater growth under the treatments compared to the vector control, especially under the 100 μM ABA and 85 mM NaCl conditions (Figure 8A).

FIGURE 8

GmERF135 transgenic soybean hairy roots enhanced tolerance to ABA and salt. (A) Phenotypes of WT and GmERF135 transgenic soybean hairy roots grown on MS medium with different concentrations of ABA and NaCl. Vc for Vector control. (B) Dry weight of WT and transgenic soybean hairy roots after 2 weeks of treatment. Data are shown as means ± SD of three biological experiments. Student’s t test was used for statistical analysis and asterisks indicate significant differences from the vector control at *P < 0.05 and ∗∗P < 0.01.

Dry weight measurement results showed that GmERF135 hairy roots increased by about 10% compared to the vector under the MS0 condition, while transgenic hairy roots displayed a 6.24–75% improvement compared to WT under NaCl and exogenous ABA at the different concentrations (Figure 8B). For NaCl treatment, there has a 6.24–75.03% increase in the average dry weight of hairy roots overexpressing GmERF135. For exogenous ABA, the transgenic hairy roots displayed a 36–65% improvement compared to vector control under exogenous ABA treatment at different concentrations (Figure 8B). It is noteworthy that a total of 64.46% improvement in GmERF135 hairy roots was observed under the 100 μM ABA condition, and an improvement of 75.03% with 85 mM NaCl (Figure 8B). Results of dry weight measurements confirmed the conclusion that GmERF135 could promote plant growth under ABA and NaCl conditions.

Discussion

Abscisic acid is well-known to be involved in responses to multiple stresses, such as drought, salt, and cold, and to induce expression of stress-related genes in plants (Leung et al., 1997; Rohde et al., 2000; Rook et al., 2001; Xu et al., 2007). The interaction of the ABA and ET signal pathway is extremely complex and intricate. AREB/ABFs, members of the bZIP family of transcription factors, participate in the ABA signal pathway which can specifically bind to ABRE cis-element to modulate the expression of downstream target genes (Bonetta and McCourt, 1998; Leung and Giraudat, 1998; Li et al., 2012). Fujita et al. (2005) observed that AREB1 can act as a transcription activator in ABRE-dependent ABA signaling, which enhances drought tolerance in Arabidopsis. Overexpression of the AREB1 gene in Arachis hypogaea could enhance drought tolerance by modulating ROS scavenging and maintaining content (Li et al., 2013). Phosphorylated AREB/ABFs may bind to the ABRE cis-element of GmERF135 promotor region to activate its transcription level (Fujita et al., 2013), which could activate or repress the transcription of targeted genes. In our study, analysis of the GmERF135 promoter region showed that three ABREs were located at the promoter region, which may be binding sites of AREB/ABFs to affect the transcript of GmERF135. GmERF135 has a conversed AP2/ERF domain, which could specifically bind to GCC-box and/or dehydration-responsive element/C-repeat (DRE/CRT) cis-acting elements to modulate the expression of pathogenesis- and abiotic stresses-related genes (Hao et al., 1998; Chakravarthy et al., 2003; Zarei et al., 2011).

Previous studies showed there is an antagonistic interaction between ABA and JA-ET signaling pathways ABA for disease response (Beaudoin et al., 2000; Anderson et al., 2004; Martin-Rodriguez et al., 2011). A recent study showed that SlAREB1 could enhance the expression of ethylene biosynthetic genes ACS and ACO in tomato fruits, which are the two key genes of ET synthesis (Mou et al., 2018). In our study, qRT-PCR analysis showed that the transcription of both AtACO4 and AtACS2 increased (Supplementary Figure S2) in transgenic Arabidopsis lines, which suggested acceleration of ET production. ET acceleration could trigger a series of reactions of ET and activate the transcript of ERF responsive genes (Guo and Ecker, 2004). At the same time, the transcript levels of ABA1 and ABA2, key factors of ABA synthesis, were also upregulated after overexpressing GmERF135 in Arabidopsis lines. These results suggested that GmERF135 may participate in both ET and ABA signaling pathways, and the regulation between the two signaling pathways needs further research.

Except for ABREs, other stress-related elements were also distributed in the promotor region of GmERF135, such as the MYBST1 core sequence, MBS, and ERE (Table 1). The MYBST1 core sequence and MBS elements have been demonstrated to be involved in drought, low temperature, salt, and ABA stress responses (Abe et al., 2003). qRT-PCR showed that GmERF135 is induced by multiple stresses, including drought, salt, low temperature, exogenous ET, SA, and JA. These changes may be caused by some corresponding cis-elements in promoter regions which can be bound by specific transcription factors (Huang et al., 2012; Buscaill and Rivas, 2014). These results suggest that soybean ERF gene GmERF135 is a key factor which participates in multiple signaling pathways to regulate expression levels of stress-related genes.

Author Contributions

Z-SX and Y-ZM coordinated the project, conceived and designed experiments, and edited the manuscript. M-JZ and L-JY conducted the bioinformatic work, performed experiments and wrote the first draft. JM revised the manuscript and figures. J-CZ conducted the bioinformatic work. Y-XW, J-HL, and J-DF contributed with valuable discussions. MC provided analytical tools and managed reagents. All authors have read and approved the final manuscript.

Conflict of Interest Statement

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.