Zhuomi Xie
Chuyun Yang
Siyi Liu
Mingjie Li
Li Gu
Xin Peng
Zhongyi Zhang- 1College of Agriculture, Fujian Agriculture and Forestry University, Fuzhou, China
- 2Key Laboratory of Ministry of Education for Genetics, Breeding and Multiple Utilization of Crops, Fujian Agriculture and Forestry University, Fuzhou, China
- 3Ningbo Municipal Hospital of TCM, Affiliated Hospital of Zhejiang Chinese Medical University, Ningbo, China
- 4Medicinal Plant Resource Center, Ningbo Research Institute of Traditional Chinese Medicine, Ningbo, China
Tetrastigma hemsleyanum (T. hemsleyanum) is a traditional medicinal plant that is widely used in China. Cultivated T. hemsleyanum usually encounters cold stress, limiting its growth and quality at key developmental stages. APETALA2 (AP2)/ethylene-responsive factor (ERF) transcription factors (TFs) comprise one of the largest gene superfamilies in plants and are widely involved in biotic and abiotic stresses. To reveal the roles of AP2/ERF TFs during T. hemsleyanum development, 70 AP2/ERF TFs were identified in T. hemsleyanum. Among them, 18 and 2 TFs were classified into the AP2 and RAV families, respectively. The other 50 TFs belonged to the ERF family and were further divided into the ERF and (dehydration reaction element binding factor) DREB subfamilies. The ERF subfamily contained 46 TFs, while the DREB subfamily contained 4 TFs. Phylogenetic analysis indicated that AP2/ERF TFs could be classified into five groups, in which 10 conserved motifs were confirmed. Several motifs were group- or subgroup-specific, implying that they were significant for the functions of the AP2/ERF TFs of these clades. In addition, 70 AP2/ERF TFs from the five groups were used for an expression pattern analysis under three low-temperature levels, namely, –4, 0, and 4°C. The majority of these AP2/ERF TFs exhibited a positive response to cold stress conditions. Specifically, ThERF5, ThERF31, ThERF46, and ThERF55 demonstrated a more sensitive response to cold stress. Moreover, AP2/ERF TFs exhibited specific expression patterns under cold stress. Transient overexpression and RNA interference indicated that ThERF46 has a specific tolerance to cold stress. These new insights provide the basis for further studies on the roles of AP2/ERF TFs in cold stress tolerance in T. hemsleyanum.
Introduction
Tetrastigma hemsleyanum (T. hemsleyanum) is one of the traditional medicinal plants that is widely used in China. The root tuber of T. hemsleyanum has unique antiviral, antitumor, and anti-inflammatory effects. The growth and development of cultivated T. hemsleyanum are mainly affected by abiotic stress, with cold stress being the key environmental factor limiting its growth and quality. Transcription factors (TFs), which can specifically regulate the gene expression of genes by binding to them (Dimova et al., 2003), play critical roles in related regulatory networks or signaling pathways responding to abiotic stress (Sharma et al., 2010).
The APETALA2/ethylene-responsive factor (AP2/ERF) is one of the largest TF families in plants. The protein encoded by the AP2/ERF TF is composed of 60–70 amino acids, mainly including a DNA binding domain, a transcriptional activation or inhibition domain, an oligomerization site, and a nuclear localization signal (Nakano et al., 2006). The biological functions of the AP2/ERF gene family have been extensively studied in model plants. The AP2/ERF family has been shown to play an important role in plant growth, development (Kubo and Kakimoto, 2000; Rashotte et al., 2006; Dietz et al., 2010), and abiotic stress response (Liu et al., 1998; Xu et al., 2011). Moreover, AP2/ERF TFs have been implicated in the plant defense response, signal transmission, stress response, gene expression regulation, plant growth, and developmental regulation. In these processes, different domains of the AP2/ERF family play different roles. The AP2/ERF family is divided into four subfamilies according to the characteristics of their conserved domains, namely, ethylene reaction element binding factor (ERF), dehydration reaction element binding factor (DREB), AP2 related to ABI3/VP1 (RAV), and soloist. Both ERF and DREB subfamily members contain an AP2 domain (Sharma et al., 2010). The Ap2 domain can regulate the binding activity of cis-acting elements, such as the dehydration response element/C repeat (DRE/CRT) and GCC box in the promoter region of AP2/ERF TFs’ target genes (Ohme-Takagi and Shinshi, 1995; Jofuku et al., 2005). Ap2 subfamily members contain two adjacent AP2 domains, whereas RAV subfamily members have both an AP2 domain and a B3 domain (Nakano et al., 2006). The whole-genome identification and analysis of the AP2/ERF family have been carried out in several plant species, and AP2/ERF family has been deeply studied in many plants, including Arabidopsis (Nakano et al., 2006), rice (Sharoni et al., 2011), alfalfa (Shu et al., 2016), and other plants (Zhuang et al., 2008, 2011; Licausi et al., 2010; Pirrello et al., 2012; Song et al., 2013).
The AP2 subfamily is involved in regulating plant flowering time (Jofuku et al., 1994), flower organ growth and development (Aukerman and Sakai, 2003), ovule development (Klucher et al., 1996), determining spikelet meristem formation (Guillaumot et al., 2008), and leaf epidermal cell formation (Moose and Sisco, 1996). However, the AP2/ERF gene family is regulated by hormones such as ethylene and brassinolide, and its DNA-binding domain can directly interact with the GCC box to regulate the expression of downstream target genes (Lata et al., 2014; Xu et al., 2020). Overexpression of the AP2/ERF gene can improve the resistance of plants to salt stress and drought stress (Sohn et al., 2006; Li et al., 2015). AP2/ERF TFs are also involved in salt stress (Wu et al., 2008) and high- and low-temperature stress (Ito et al., 2006; Qin et al., 2007). These findings have greatly promoted the understanding of the biological function of AP2/ERF family TFs. Moreover, numerous studies have identified the roles of AP2/ERF TFs in plants, but the role of AP2/ERF TFs in T. hemsleyanum is still unclear.
This study was the first to identify AP2/ERF TFs in T. hemsleyanum using a bioinformatics approach. Moreover, it analyzed in detail the expression pattern of the AP2/ERF gene family under cold stress. The functions of the AP2/ERF gene were identified in detail by transient overexpression and RNA-interference expression. The results provide new insights for further investigation of the molecular mechanism of AP2/ERF TFs in T. hemsleyanum under cold stress.
Materials and methods
Plant materials and stress treatments
Healthy cutting seedlings of T. hemsleyanum were grown under standard greenhouse conditions set as follows: a 16-h day/8-h night cycle, a 23 ± 2°C temperature, and 50% relative humidity. Then, 6-month-old seedlings of T. hemsleyanum with identical growth potentials were selected for a transient expression experiment of purpose genes. Agrobacterium tumefaciens strain GV3101 containing a 35S promoter was used to transiently overexpress or interfere with AP2/ERFs. The density of the GV3101 bacterium solution was amplified to OD600 = 1.0 and set to 5,000/rpm for 5 min. The bacteria were resuspended with the infection solution (containing 100 μM acetosyringone, 10 mM 2-morpholinoethanesulfonic acid, 10 mM MgCl2, and pH = 5.7), and the concentration of the solution was adjusted to OD600 = 1.0. The infection solution was injected into the back of leaves by syringe, and a total of 25 plants were infected as repeats. After 72 h of infection, the expression levels of ThERF46 were detected to determine whether the gene was expressed or silenced by qRT-PCR. When the transient expression was successful, T. hemsleyanum was treated at 4°C for 6 h. Phenotypic characteristics and electrolyte leakage were carefully analyzed. Simultaneously, the leaves were sampled, immediately frozen in liquid nitrogen, and stored at –80°C for further analysis on gene expression and biochemical index. The gene expression and biochemical index analysis were repeated three times, and each biological replicate consisted of a sample pool of 25 seedlings.
Identification of AP2/ERF transcription factors in Tetrastigma hemsleyanum
Transcriptome data were obtained in our previous study (accession number: PRJNA797653) (Xie et al., 2022). AP2/ERF TF family proteins of Arabidopsis thaliana were compared with the T. hemsleyanum transcriptome data, and the T. hemsleyanum transcriptome was searched using the HMMER software domain. A hidden Markov model (HMM) of the AP2 (PF00847) was downloaded from the Pfam database1 (El-Gebali et al., 2019), and the AP2/ERF family TFs were identified using HMMER3 (version 3.0) software (Eddy, 2011) with a defined threshold of E < 1e–5. The NCBI Conserved Domain Search Service (CD Search) (Marchler-Bauer et al., 2017) was used to confirm manually the predicted AP2/ERF family TFs.
Conserved motif and evolutionary relationship analysis
The conserved motifs in the T. hemsleyanum AP2/ERF TFs were identified using the online motif finding tool MEME 4.11.22 (Bailey et al., 2009). Arabidopsis and Vitis vinifera ERF-related proteins were downloaded from the TAIR website3 and NCBI, respectively. The evolutionary relationship of three kinds of ERF-related proteins, namely, AtERF, VvERF, and ThERF, was analyzed using the MEGA 7.0 software.
Subcellular localization prediction and gene structure analysis
The subcellular localization and physicochemical properties of ThERF proteins were analyzed using the Wolf PSORT software4 and ProtParam software5, respectively. Three-dimensional models of the ThERF proteins were analyzed using the SWISS-MODEL online software6.
Expression pattern analysis
Total RNA was isolated from 100 mg of leaves (fresh-weight) using a plant RNA extraction kit (Nanjing Vazyme Biotech Co., Ltd.). cDNA synthesis was performed with the Evo M-MLV Mix Kit with gDNA Clean for qPCR AG11728 [Accurate Biotechnology (Hunan) Co., Ltd.]. Each reaction contained 10 μl of 2*SYBR Green Pro Taq HS Premix AG11701 [Accurate Biotechnology (Hunan) Co., Ltd.], 2 μl of template cDNA, and 0.4 μl of each forward and reverse primers (10 μM). The GAPDH gene was used as an internal reference. All primers used in this assay are shown in Supplementary Table S1. The PCR reaction procedure was performed as follows: incubation at 95°C for 2 min followed by 40 cycles at 95°C for 5 s and 60°C for 30 s. Each gene was tested in biological triplicates with three technical repeats. The expression level for each sample was expressed as 2–ΔΔCT. The data were exhibited as the mean ± SD of three independent experiments.
Construction of expression vectors and subcellular localization
To construct a transient overexpression vector for ThERF46, its full-length of ORF was cloned into the entry vector pBI121-GFP digested by KpnI using the following primers: ThERF-forward, 5′ggtaccATGGCGGTCGAGCCCCTC3′ and ThERF-reverse, and 5′ggtaccTCAGAAGAGCGGGGCCGG3′. To build an RNA-interference vector for ThERF46, a fragment from ORF of ThERF46 was amplified using the following primers: ThERF-forward, 5′accag gtctcaggagATCTCATCCTTCCCGTTCTG3′ and ThERF-reverse, and 5′accaggtctcatcgtGGCAGCATCGTCGTAAGC3′. This fragment was then inserted into the entry vector pRNAiGG. To detect the subcellular localization of ThERF46, its ORF was further inserted into the pBI121-GFP vector to form fused proteins with GFP. The expression of ThERF46-GFP fusion was driven by double 35S promoters. Agrobacterium tumefaciens strain GV3101 containing 35S: ThERF46: GFP and 35S: GFP were grown overnight in LB culture solution, and then resuspended to OD600 = 1.0. Agrobacterium tumefaciens strains harboring the GFP fusion constructs were filtrated into N. benthamiana leaves.
Biochemical index analysis
The transient expression of T. hemsleyanum was used to determine superoxide dismutase (SOD, EC 1.15.1.1), peroxidase (POD, EC 1.11.1.7) activity, and malondialdehyde (MDA) content. POD activity was determined as guaiacol oxidation by H2O2. SOD activity was analyzed based on the inhibiting rates of the reduction of nitro blue tetrazolium (NBT), whereas MDA content was extracted and determined using the thiobarbituric acid reaction method according to our previous report (Peng et al., 2019). Proline content was determined using a colorimetric method (Bates et al., 1973). Leaf electrolyte leakage (EL) was calculated based on the formula (%) = Cinitial/Cmax × 100, where Cinitial indicates the initial conductivity and Cmax represents the max conductivity. Fresh leaves (0.2 g) were soaked instantly in 50 ml of distilled water for 24 h to detect the C−initial. These samples were autoclaved at 120°C for 20 min and cooled down to room temperature to detect the C−max using a conductivity meter (YSI Model 32, Yellow Spring, OH, United States) (Blum and Ebercon, 1981). Evans blue was used to detect the activity of leaf cells in different treatments according to our previous study (Xie et al., 2022). In brief, the leaves were immersed in an aqueous solution of 1 mg⋅ml–1 Evans blue and incubated for 12 h in the dark. The leaves were set to the ethanol:lactic acid:glycerol (3:1:1) mixture to boil for 5 min and stored in 60% glycerol before photography.
Results
Identification, physicochemical properties, and subcellular localization prediction of Tetrastigma hemsleyanum’s AP2/ERF transcription factors
A total of 104 ThERF sequences were identified in the T. hemsleyanum transcriptome. These sequences were first annotated using the NR and Swiss-Prot databases, then they were identified and analyzed using the NCBI-CDD online website. Finally, 70 ThERF sequences were obtained. The physicochemical properties of the proteins encoded by the ThERF TF family were predicted. The results indicated that the protein length, molecular weight, and theoretical isoelectric point of ThERF TFs were 128–684 aa, 13.93–76.07 kD, and 4.34–9.85, respectively. The instability coefficients were in the range of 32.72–81.44, with 94.28% of the proteins being unstable proteins with an instability coefficient of over 40. The average hydrophilic coefficients were lower than 0, and all ThERFs were hydrophilic proteins. The results of the subcellular localization prediction suggested that 52 ThERF TFs were located in the nucleus, accounting for 74.28%; 10 ThERF TFs were located in the cytoplasm, accounting for 14.28%; and 7 ThERF TFs were located in the cytoplasm and nucleus, accounting for 10%. Notably, ThERF5 was located in the cytoplasm and mitochondria (Table 1). The ThERF proteins were mainly composed of acidic and unstable hydrophilic proteins, which were primarily expressed in the nucleus.
Table 1. The AP2/ERF protein information of T. hemsleyanum.
Identification of the conserved motifs of Tetrastigma hemsleyanum’s AP2/ERF transcription factors
Multiple sequence alignments and conserved domains of AP2/ERF proteins were identified using the DNAMAN software and NCBI CD search online tools, respectively (Supplementary Figure S1). All AP2/ERF proteins were divided into cl00033 (95.71%), cl15242 (2.85%), and cl11268 (1.42%). Cl00033 was located at the C-terminal of the ThERF8 protein domain (Supplementary Table S2). The conserved motifs of 70 AP2/ERF proteins were analyzed using the MEME online tool. A total of 10 conserved motifs were detected (Figure 1). Motif 1, motif 2, and motif 3 located all AP2/ERF proteins (100%). Motif 4 to motif 10 existed mainly in ThERF17, ThERF30, ThERF52, ThERF64, and ThERF66. Motif 4 was also identified in ThERF6, ThERF7, ThERF10, ThERF16, and ThERF37. Motif 4 to motif 10 were unavailable in AP2/ERF proteins, suggesting that they may have been caused by the long evolution of T. hemsleyanum.
Figure 1. Distribution of conserved motifs of the AP2/ERF transcription factors in T. hemsleyanum. Different motifs are highlighted in different colors.
Phylogeny of Tetrastigma hemsleyanum’s AP2/ERF transcription factors
To confirm the classification and evolutionary relationships of AP2/ERF TFs in T. hemsleyanum, the full-length sequences of the putative proteins were compared, and a phylogenetic tree analysis was conducted. All AP2/ERF proteins could be classified into five groups (Figure 2). Group I included 34 TFs, containing motif 1, motif 6, and motif 9. Group II included 22 TFs, almost all of which contained motif 1 and motif 6. Motif 9 was not found in ThERF35, ThERF45, and ThERF68. Groups I and II were considered members of the DREB subfamily. A large proportion of the TFs in these groups containing two AP2 domains was classified into the AP2 family. Group III, similar to the second group, contained 3 TFs. ThERF5 and ThERF67 included motif 1 and motif 6 but had lost motif 9. Most TFs in this group were classified into the RAV family. Group IV included 11 TFs, and almost all ThERF TFs contained motif 4 and motif 8. The TFs, including ThERF7, ThERF30, ThERF64, and ThERF66 in other clades, contained almost all identified motifs.
Figure 2. Phylogenetic tree of AP2/ERF transcription factors in T. hemsleyanum. Deduced full-length amino acid sequences were used to construct the phylogenetic tree with the MEGA 7.0 software using a neighbor-joining method with 1,000 bootstrap replicates. Five groups are highlighted in different colors.
Three-dimensional structure models of Tetrastigma hemsleyanum’s AP2/ERF transcription factors
Swiss-model software was used to predict three-dimensional (3D) structure models of key ThERF TFs (Supplementary Figure S2). As a result, ThERFs could be divided into four types of structures, and the identified similarity of all sequences ranged from 36.99 to 98.53. Group I consisted of 10 ThERF proteins, including ThERF31, ThERF32, ThERF40, ThERF44, ThERF45, ThERF46, ThERF48, ThERF49, ThERF55, and ThERF56, which had been identified as GCC-box containing, ethylene-responsive TFs. Group II and group III included ThERF5 and ThERF67, and ThERF35 and ThERF37 with similar 3D models, which were identified as the DNA-binding protein RAV1 solution structures of the B3 DNA-binding domain of RAV1 and a B3 domain-containing transcription reporter, respectively. Group IV contained ThERF14, which was described as a NAD (P) H-quinone oxidoreductase subunit M chloroplast NDH complex special structure (Supplementary Table S3). The results indicated that these ThERF proteins had different functions in the growth and development of T. hemsleyanum.
Expression profiling of AP2/ERF genes under cold stress
Numerous studies suggested that the AP2/ERF family plays an important role in plants’ resistance to cold stress. To prove this conjecture, the expression profiles of AP2/ERF TFs were explored under different cold stress conditions for 6 h. All 70 ThERFs were shown to have different expression levels under cold stress (Supplementary Figures S3–S5). To facilitate the observation, the expression levels of all ThERFs in T. hemsleyanum were described as a heat map (Figure 3A). In addition, the 15 key ThERFs, including ThERF5, ThERF14, ThERF31, ThERF32, ThERF35, ThERF37, ThERF40, ThERF44, ThERF45, ThERF46, ThERF48, ThERF49, ThERF55, ThERF56, and ThERF67 were selected for detailed analysis (Figure 3B). A total of 13 ThERFs increased 1.89-fold to 68.4-fold compared with the control (P < 0.01), while the relative expression levels of ThERF48 and ThERF56 decreased 0.39-fold and 0.94-fold compared with the control, respectively. The expression levels of 15 key ThERFs increased significantly at 0 and -4°C, increasing 1.28-fold to 174.57-fold and 1.37-fold to 90.57-fold compared with the 25°C treatment, respectively (P < 0.01). The findings suggested that these key ThERFs were significantly influenced by cold stress and played an important role in the resistance to cold stress.
Figure 3. Expression profiles of AP2/ERF TFs at (–4, 0, and 4°C) cold treatments. (A) Heat map for gene expression of 70 ThERFs in T. hemsleyanum under cold treatments. (B) Gene expression of 15 key ThERFs in T. hemsleyanum under cold treatments. The heat map is constructed based on the Log2 (relative expression level of each transcription factor) detected by qRT-PCR. The color scale represents the Log2 (relative expression level of each transcription factor), with white denoting low expression and red denoting high expression. Data are means ± standard error of the three biological replicates for each treatment. *P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001.
Identification of the functions of ThERF46 in Tetrastigma hemsleyanum under cold stress
To further determine the biological function of ThERFs, the ORFs of 15 key ThERFs were predicted, and the complete CDS sequence of ThERF46 was cloned in overexpression and RNA-interference expression vectors. The transient ThERF46-OE and ThERF46-Ri leaves were subjected to cold stress for 6 h. After 6 h of 4°C cold treatment, it can be clearly found that the leaves of the empty vector (EV) and ThERF46-Ri treatment had water stains and wilted, and could not extend completely under cold stress by phenotypic observation. In addition, after transient overexpression of ThERF46, the leaves did not show obvious damage from cold stress. The ThERF46-Ri and EV leaves were damaged more severely than ThERF46-OE leaves (Figure 4A). Subcellular localization analysis found that ThERF46 was located in the nucleus and cell membrane of N. benthamiana (Figure 4B). The expression levels of ThERF46-OE and ThERF46-Ri were significantly increased and decreased 46-fold and 0.3-fold compared with the EV leaves, respectively. Under cold stress conditions, the MDA content and relative electrolyte leakage of ThERF46-OE and ThERF46-Ri leaves were quantified. The resulting MDA content and relative electrolyte leakage leaves were significantly decreased to 1.6 ± 3.28% and 31.2 ± 6.32% in ThERF46-OE compared with EV, respectively. The SOD enzyme activity, POD enzyme activity, and proline content showed a dramatic increase to 421.4 ± 12.3%, 223.8 ± 11.2%, and 27.3 ± 9.1%, respectively (Figure 5A). Moreover, the results also showed there was a higher content of proline and less damage area in ThERF46-OE leaves compared with EV and ThERF46-Ri treatments (Figure 5B). These results showed that the transient overexpression of ThERF46 reduced the cell membrane damage and positively regulated cold tolerance by enhancing the activity of the plant cell and antioxidant enzyme system of T. hemsleyanum under cold stress.
Figure 4. Phenotype and subcellular localization analysis of ThERF46 in T. hemsleyanum. (A) The phenotype of ThERF46-OE, ThERF46-Ri, and EV treatment under cold stress. EV, empty vector (the EV without an exogenous gene was used as the control); R1, Repeat 1; R2, Repeat 2. (B) Subcellular localization of ThERF46 in N. benthamiana. GFP, green fluorescence protein; RFP, red fluorescence protein (used as a marker of nucleus localization signals). Three biological replicates for each treatment. Bar, 0.5 cm.
Figure 5. Analysis of the biochemical indexes of ThERF46 in T. hemsleyanum. (A) Analysis of the physiological indexes of T. hemsleyanum. (B) Analysis of cell activity of T. hemsleyanum. EV, empty vector (the EV without an exogenous gene was used as the control); R1, Repeat 1; R2, Repeat 2. Data are means ± standard error of the three biological replicates for each treatment. Bar, 1 cm. *P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001.
Discussion
In contrast to animals, plants must be subjected to different biological and abiotic stresses during their growth and development, such as temperature, light, precipitation, humidity, and soil conditions. Among them, cold stress is the key factor affecting plant growth and development (Chinnusamy et al., 2010). The special growth requirements of T. hemsleyanum are an important problem challenge for the artificial cultivation industry. The quality and quantity of T. hemsleyanum have always suffered from freezing damage. Response to cold stress is a very complex process in plants, involving various pathways, metabolic pathways, gene expression, intercellular changes, cold signal transduction, gene transcription, and post-translational alterations in plants (Chinnusamy et al., 2010). Several mechanisms would minimize the potential damage caused by cold stress in plants, involving a series of physiological and biochemical modifications (Guo et al., 2021). During this period, the expression of many genes, proteins, and metabolites in plants would also change in response to cold stress (Guo et al., 2021). Therefore, it is important to further study the response mechanism of T. hemsleyanum under cold stress.
AP2/ERF TFs are one of the most specific and largest TF families in plants (Rashid et al., 2012; Licausi et al., 2013). They not only participate in the regulation of plant growth and development but also play an important role in plant stress. In recent years, with the development of sequencing technology in plant genomes, the distribution of AP2/ERF TFs in different plants has been gradually revealed. The AP2/ERF TF subfamily has also been studied extensively in many plants. More than 100 AP2/ERF TFs have been identified and investigated in many plants, but there has not been a relevant investigation of the AP2/ERF family in T. hemsleyanum yet (Zhuang et al., 2008, 2011; Licausi et al., 2010; Sharoni et al., 2011; Pirrello et al., 2012; Song et al., 2013; Shu et al., 2016). In this study, 70 AP2/ERF TFs were identified in T. hemsleyanum, which were fewer than those identified in Vitis vinifera (149). This may be due to long-term evolutionary relationships or maybe an inevitable factor caused by transcriptome data. A further phylogenetic tree analysis of AP2/ERF among T. hemsleyanum, Vitis vinifera, and Arabidopsis indicated that the number of TFs in the AP2 family and RAV and DERB subfamilies in T. hemsleyanum was significantly lower than that in Vitis vinifera and Arabidopsis. The findings suggested that the difference in the number of AP2/ERF TFs in T. hemsleyanum was mainly caused by the contractions of AP2 and DREB TFs and the amplification of some ERF TFs (Li et al., 2017). VvERF4 in Vitis vinifera was classified into a single group, indicating a long evolutionary distance of genetic relationships between T. hemsleyanum and Vitis vinifera TFs, which may be the result of a long evolution.
In the phylogenetic analysis, a few ThERF TFs were identified in group I, which included most AtERF TFs. ThERF37, ThERF45, ThERF56, and ThERF67 in T. hemsleyanum differed at the evolutionary level but also showed similar expression changes, suggesting that they had the same functions under cold stress (Licausi et al., 2010; Zhang et al., 2022). In general, the domains or amino acid motifs of TFs are often involved in nuclear localization, DNA-binding, protein-protein interaction, and transcriptional activity (Nakano et al., 2006). The B3 domain in RAV TFs played a vital role in abiotic stress and pathogen resistance (Kim et al., 2005; Cao et al., 2020; Cui et al., 2021; Zhang et al., 2021). The core sequence of the B3 special domain was identified in ThERF5, ThERF35, and ThERF67, which belonged to the RAV TF family. The RAV (B3) special domains were confirmed to be a transcriptional repressor in plants (Hu et al., 2004; Mittal et al., 2015). The qRT-PCR results showed that cold stress activated the expression of ThERF5, ThERF35, and ThERF67, which contained a B3 domain. These findings indicate that a function of the B3 domain in ThERFs under cold stress should be further explored. In addition, most AP2/ERF TFs with an AP2 domain in T. hemsleyanum were specifically expressed at high levels during cold stress resistance, suggesting that they played the same functions as the B3 domains. Some domains or motifs of ThERFs were highly conserved, and the newly evolved motifs had been generated in the long evolutionary process of T. hemsleyanum, which might play an important role in the subfunctionalization or new functions of AP2/ERF TFs in specific plant species. Notably, the Ndhm domain, which was annotated as a medium for electron transfer during the process of PS I in vascular plants, was identified in ThERF14 (Yamamoto et al., 2016). The Ndhm domains were not common in the AP2/ERF TF family. This may be a specific domain in the AP2/ERF TF family in T. hemsleyanum. The function and regulation of these newly evolved motifs in AP2/ERF TFs should be further explored, and this study provides important insights for further explaining the evolution and functions of AP2/ERF TFs.
To further study the potential mechanisms of AP2/ERF TFs in response to cold stress, the expression patterns of 70 AP2/ERF TFs were analyzed. Almost all 70 AP2/ERF TFs were activated under cold stress, and the 15 key ThERFs rapidly responded to the beginning cold stress. These findings are in line with previous reports on the involvement of AP2/ERF TFs in abiotic stress (Mizoi et al., 2012; Velivelli et al., 2015; Shu et al., 2016). Then, the relevant resistance genes were expressed at high levels to resist cold stress during the intermediate stages of cold stress processes. Eleven key AP2/ERF proteins, including ThERF31, ThERF32, ThERF37, ThERF40, ThERF44, ThERF45, ThERF46, ThERF48, ThERF49, ThERF55, and ThERF56, showed the same expression profiles under these conditions. These findings indicate that ThERFs play an important role in the tolerance to cold stress in T. hemsleyanum, and their role should be explored further. Then, with stress exacerbation, plants could not resist the stress only via their own regulation, resulting in metabolic disorder, cell activity, and decreased gene expression (Fowler and Thomashow, 2002; Dong and Liu, 2010). During this period, the expression levels of ThERFs also showed a downward trend compared to the initial stage of cold stress (Figure 3B).
Plant cell tissue is the primary system to respond to cold stress, and cold stress affects the membrane of plant cells, which is considered to be the main reaction to trigger the cold stress response in plants (Orvar et al., 2000). The membrane system changes mainly in membrane lipid transformation and membrane lipid peroxidation under cold stress conditions (Ye et al., 2019; Zhuang et al., 2019). Membrane permeability is also an important physiological index to evaluate cold resistance during membrane lipid transformation. The electrolyte leakage was negatively correlated with plant resistance under cold stress. In the process of membrane lipid peroxidation, the activity of SOD and POD enzymes is rapidly induced under cold stress and decreases the accumulation of H2O2 and ROS in plant cells (Theocharis et al., 2012). In addition, MDA is the final product of lipid peroxidation, and over-accumulation of MDA can change the structure of the cell membrane and aggravate membrane lipid phase transformation and membrane damage (Huang et al., 2017; Yu et al., 2018; Peng et al., 2019). Proline is very hydrophilic. It can stabilize the protoplast colloid and the metabolic process and reduce the chilling injury to prevent cell dehydration (Guo et al., 2021). Therefore, the content of these molecules and enzyme activity are widely used to evaluate the resistance of plants under cold stress. Gene function analysis was considered an important approach to understanding the molecular mechanisms of plant response to abiotic stress. In this study, ThERF46 was identified and cloned, and subcellular localization determined its functional location. The detection of physiological functions under cold stress confirmed the role of ThERF46 in the resistance of T. hemsleyanum. Transient overexpression of ThERF46 increased the tolerance by enhancing the activity of cell and antioxidant enzyme systems of T. hemsleyanum, while transient RNA-interference expression of ThERF46 decreased the tolerance of cold stress in T. hemsleyanum. These results indicate that ThERF46 plays an important role in T. hemsleyanum during the cold stress process. Nevertheless, although several AP2/ERF TFs have been proven to play an important role in the tolerance of plants to various biological and abiotic stresses, the function of AP2/ERF TFs in T. hemsleyanum is still largely unknown. This study identified potential candidate genes for the AP2/ERF TFs, especially ThERF46, to further explore their role in T. hemsleyanum under cold stress.
Conclusion
In conclusion, 70 AP2/ERF TFs were identified in T. hemsleyanum. Clustering and phylogenetic analyses were conducted to divide these TFs into five groups. Ten conserved motifs were identified in the 70 AP2/ERF TFs. AP2/ERF TFs exhibiting specific expression patterns under cold stress were also confirmed. Transient overexpression and RNA-interference expression of ThERF46 increased and decreased the tolerance to cold stress, respectively. These findings provide new insights into the AP2/ERF TFs in T. hemsleyanum and identify candidate AP2/ERF TFs to further elucidate their roles in cold stress tolerance.
Data availability statement
The transcriptome datasets presented in this study can be found in the NCBI database (PRJNA797653).
Author contributions
ZZ and XP conceived the study. ML designed the experiment. ZX performed the physiological indicator determination and expression analysis. CY and LG performed the RNA extraction and quality determination. SL carried out the analysis. All authors have read and approved the manuscript.
Funding
This study was supported by the National Natural Science Foundation of China (81903743), China Postdoctoral Science Special Foundation (2019T120552), Ningbo City Science and Technology Innovation 2025 Major Research Project (2019B10008), Major Science and Technology Projects of Breeding New Varieties of Agriculture in Zhejiang Province (2021C02074), and Zhejiang Provincial Natural Science Foundation (LQ19H280005).
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.936602/full#supplementary-material
Supplementary Figure 1 Multiple sequence alignments of AP2/ERF TFs in T. hemsleyanum.
Supplementary Figure 2 Protein 3D model of AP2/ERF transcription factors in T. hemsleyanum. Full-length amino acid sequences were used to construct the 3D model (similarity model > 30%).
Supplementary Figure 3 Gene expression of ThERF1–ThERF25 in T. hemsleyanum under cold stress.
Supplementary Figure 4 Gene expression of ThERF26–ThERF50 in T. hemsleyanum under cold stress.
Supplementary Figure 5 Gene expression of ThERF51–ThERF70 in T. hemsleyanum under cold stress.
Footnotes
- ^ http://pfam.Xfam.org/
- ^ http://meme-suite.org/tools/meme
- ^ http://arabidopsis.org
- ^ https://wolfpsort.hgc.jp/
- ^ https://web.expasy.org/protparam/
- ^ https://swissmodel.expasy.org/
References
Aukerman, M. J., and Sakai, H. (2003). Regulation of flowering time and floral organ identity by a MicroRNA and its APETALA2-like target genes. Plant Cell. 15, 2730–2741. doi: 10.1105/tpc.016238
Bailey, T. L., Boden, M., Buske, F. A., Frith, M., Grant, C. E., Clementi, L., et al. (2009). MEME SUITE: Tools for motif discovery and searching. Nucleic Acids Res. 37, W202–W208. doi: 10.1093/nar/gkp335
Bates, L. S., Waldren, R. P., and Teare, I. D. (1973). Rapid determination of free proline for water-stress studies. Plant Soil. 39, 205–207. doi: 10.1007/BF00018060
Blum, A., and Ebercon, A. (1981). Cell membrane stability as a measure of drought and heat tolerance in wheat. Crop Sci. 21, 43–47.
Cao, S. L., Wang, Y., Li, X. T., Gao, F., Feng, J. C., and Zhou, Y. J. (2020). Characterization of the AP2/ERF transcription factor family and expression profiling of DREB subfamily under cold and osmotic stresses in Ammopiptanthus nanus. Plants (Basel) 9:455. doi: 10.3390/plants9040455
Chinnusamy, V., Zhu, J. K., and Sunkar, R. (2010). Gene regulation during cold stress acclimation in plants. Methods Mol. Biol. 639, 39–55. doi: 10.1007/978-1-60761-702-0_3
Cui, M. J., Haider, M. S., Chai, P. P., Guo, J. J., Du, P., Li, H. Y., et al. (2021). Arachis hypogaea genome-wide identification and expression analysis of AP2/ERF transcription factor related to drought stress in cultivated peanut (L.). Front. Genet. 12:750761. doi: 10.3389/fgene.2021.750761
Dietz, K. J., Vogel, M. O., and Viehhauser, A. (2010). AP2/EREBP transcription factors are part of gene regulatory networks and integrate metabolic, hormonal and environmental signals in stress acclimation and retrograde signalling. Protoplasma 245, 3–14. doi: 10.1007/s00709-010-0142-8
Dimova, D. K., Stevaux, O., Frolov, M. V., and Dyson, N. J. (2003). Cell cycle-dependent and cell cycle-independent control of transcription by the drosophila E2F/RB pathway. Genes Dev. 17, 2308–2320. doi: 10.1101/gad.1116703
Dong, C. J., and Liu, J. Y. (2010). The Arabidopsis EAR-motif-containing protein RAP2.1 functions as an active transcriptional repressor to keep stress responses under tight control. BMC Plant Biol. 10:47. doi: 10.1186/1471-2229-10-47
Eddy, S. R. (2011). Accelerated profile HMM searches. PLoS Comput. Biol. 7:e1002195. doi: 10.1371/journal.pcbi.1002195
El-Gebali, S., Mistry, J., Bateman, A., Eddy, S. R., Luciani, A., Potter, S. C., et al. (2019). The Pfam protein families database in 2019. Nucleic Acids Res. 47, D427–D432. doi: 10.1093/nar/gky995
Fowler, S., and Thomashow, M. F. (2002). Arabidopsis transcriptome profiling indicates that multiple regulatory pathways are activated during cold acclimation in addition to the CBF cold response pathway. Plant Cell. 14, 1675–1690. doi: 10.1105/tpc.003483
Guillaumot, D., Lelu-Walter, M., Germot, A., Meytraud, F., Gastinel, L., and Riou-Khamlichi, C. (2008). Expression patterns of LmAP2L1 and LmAP2L2 encoding two-APETALA2 domain proteins during somatic embryogenesis and germination of hybrid larch (larix x marschlinsii). J. Plant. Physiol. 165, 1003–1010. doi: 10.1016/j.jplph.2007.08.009
Guo, Q. Q., Li, X., Niu, L., Jameson, P. E., and Zhou, W. B. (2021). Transcription-associated metabolomic adjustments in maize occur during combined drought and cold stress. Plant. Physiol. 186, 677–695. doi: 10.1093/plphys/kiab050
Hu, Y. X., Wang, Y. X., Liu, X. F., and Li, J. Y. (2004). Arabidopsis RAV1 is down-regulated by brassinosteroid and may act as a negative regulator during plant development. Cell Res. 14, 8–15. doi: 10.1038/sj.cr.7290197
Huang, X. B., Shi, H. Y., Hu, Z. G., Liu, A., Amombo, E., Chen, L., et al. (2017). ABA is involved in regulation of cold stress response in bermudagrass. Front. Plant Sci. 8:1613. doi: 10.3389/fpls.2017.01613
Ito, Y., Katsura, K., Maruyama, K., Taji, T., and Yamaguchi-Shinozaki, K. (2006). Functional analysis of rice dreb1/cbf-type transcription factors involved in cold-responsive gene expression in transgenic rice. Plant Cell Physiol. 47, 141–153. doi: 10.1093/pcp/pci230
Jofuku, K. D., den Boer, B. G., Van, M. M., and Okamuro, J. K. (1994). Control of Arabidopsis flower and seed development by the homeotic gene APETALA2. Plant Cell. 6, 1211–1225. doi: 10.1105/tpc.6.9.1211
Jofuku, K. D., Omidyar, P. K., Gee, Z., and Okamuro, J. K. (2005). Control of seed mass and seed yield by the floral homeotic gene APETALA2. Proc. Natl. Acad. Sci. U.S.A. 102, 3117–3122. doi: 10.1073/pnas.0409893102
Kim, S. Y., Kim, Y. C., Lee, J. H., Oh, S. K., Chung, E., Lee, S., et al. (2005). Identification of a CaRAV1 possessing an AP2/ERF and B3 DNA-binding domain from pepper leaves infected with Xanthomonas axonopodis pv. glycines 8ra by differential display. Biochim. Biophys. Acta. 1729, 141–146. doi: 10.1016/j.bbaexp.2005.04.009
Klucher, K. M., Chow, H., and Fischer, R. (1996). The aintegumenta gene of Arabidopsis required for ovule and female gametophyte development is related to the floral homeotic gene apetala2. Plant Cell Online 8, 137–153. doi: 10.1105/tpc.8.2.137
Kubo, M., and Kakimoto, T. (2000). The cytokinin-hypersensitive genes of Arabidopsis negatively regulate the cytokinin-signaling pathway for cell division and chloroplast development. Plant J. 23, 385–394. doi: 10.1046/j.1365-313x.2000.00796.x
Lata, C., Mishra, A. K., Muthamilarasan, M., Bonthala, V. S., Khan, Y., and Prasad, M. (2014). Genome-wide investigation and expression profiling of AP2/ERF transcription factor superfamily in foxtail millet (Setaria italica L.). PLoS One 9:e113092. doi: 10.1371/journal.pone.0113092
Li, H., Wang, Y., Wu, M., Li, L. H., Li, C., Han, Z. P., et al. (2017). Genome-wide identification of AP2/ERF transcription factors in cauliflower and expression profiling of the ERF family under salt and drought stresses. Front. Plant Sci. 8:946. doi: 10.3389/fpls.2017.00946
Li, X. J., Li, M., Zhou, Y., Hu, S., Hu, R., Chen, Y., et al. (2015). Overexpression of cotton RAV1 gene in Arabidopsis confers transgenic plants high salinity and drought sensitivity. PLoS One 10:e0118056. doi: 10.1371/journal.pone.0118056
Licausi, F., Giorgi, F. M., Zenoni, S., Osti, F., Pezzotti, M., and Perata, P. (2010). Genomic and transcriptomic analysis of the AP2/ERF superfamily in vitis vinifera. BMC Geno. 11:719. doi: 10.1186/1471-2164-11-719
Licausi, F., Ohme-Takagi, M., and Perata, P. (2013). APETALA2/ethylene responsive factor (AP2/ERF) transcription factors: mediators of stress responses and developmental programs. New Phytol. 199, 639–649. doi: 10.1111/nph.12291
Liu, Q., Kasuga, M., Sakuma, Y., Abe, H., Miura, S., Yamaguchi-Shinozaki, K., et al. (1998). Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low-temperature-responsive gene expression, respectively, in Arabidopsis. Plant Cell. 10, 1391–1406. doi: 10.1105/tpc.10.8.1391
Marchler-Bauer, A., Bo, Y., Han, L. Y., He, J., Lanczycki, C. J., Lu, S. N., et al. (2017). CDD/SPARCLE: Functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 45, D200–D203. doi: 10.1093/nar/gkw1129
Mittal, A., Jiang, Y., Ritchie, G. L., Burke, J. J., and Rock, C. D. (2015). AtRA V1 and AtRA V2 overexpression in cotton increases fiber length differentially under drought stress and delays flowering. Plant Sci. 241, 78–95. doi: 10.1016/j.plantsci.2015.09.013
Mizoi, J., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2012). AP2/ERF family transcription factors in plant abiotic stress responses. Biochim. Biophys. Acta 1819, 86–96. doi: 10.1016/j.bbagrm.2011.08.004
Moose, S. P., and Sisco, P. H. (1996). Glossy15, an APETALA2-like gene from maize that regulates leaf epidermal cell identity. Genes Dev. 10, 3018–3027. doi: 10.1101/gad.10.23.3018
Nakano, T., Suzuki, K., Fujimura, T., and Shinshi, H. (2006). Genome-wideanalysis of the ERF gene family in Arabidopsis and rice. Plant Physiol. 140, 411–432. doi: 10.1104/pp.105.073783
Ohme-Takagi, M., and Shinshi, H. (1995). Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. Plant Cell. 7, 173–182. doi: 10.1105/tpc.7.2.173
Orvar, B. L., Sangwan, V., Omann, F., and Dhindsa, R. S. (2000). Early steps in cold sensing by plant cells: the role of actin cytoskeleton andmembrane fluidity. Plant J. 23, 785–794. doi: 10.1046/j.1365-313x.2000.00845.x
Peng, X., Wu, H., Chen, H., Zhang, Y., Qiu, D., and Zhang, Z. (2019). Transcriptome profiling reveals candidate flavonol-relatedgenes of Tetrastigma hemsleyanum under cold stress. BMC Geno. 20:687. doi: 10.1186/s12864-019-6045-y
Pirrello, J., Prasad, B. C. N., Zhang, W. S., Chen, K. S., Mila, I., Zouine, M., et al. (2012). Functional analysis and binding affinity of tomato ethylene response factors provide insight on the molecular bases of plant differential responses to ethylene. BMC Plant Biol. 12:190. doi: 10.1186/1471-2229-12-190
Qin, F., Kakimoto, M., Sakuma, Y., Maruyama, K., Osakabe, Y., Tran, L. S. P., et al. (2007). Regulation and functional analysis of ZmDREB2A in response to drought and heat stresses in Zea mays L. Plant J. 50, 54–69. doi: 10.1111/j.1365-313X.2007.03034.x
Rashid, M., Guangyuan, H., Guangxiao, Y., Hussain, J., and Xu, Y. (2012). AP2/ERF transcription factor in rice: Genome-wide canvas and syntenic relationships between monocots and eudicots. Evol. Bioinform. Online 8, 321–355. doi: 10.4137/EBO.S9369
Rashotte, A. M., Mason, M. G., Hutchison, C. E., Ferreira, F. J., Schaller, G. E., and Kieber, J. J. (2006). A subset of Arabidopsis AP2 transcription factors mediates cytokinin responses in concert with a two-component pathway. Proc. Natl. Acad. Sci. U.S.A. 103, 11081–11085. doi: 10.1073/pnas.0602038103
Sharma, M. K., Kumar, R., Solanke, A. U., Sharma, R., Tyagi, A. K., and Sharma, A. K. (2010). Identification, phylogeny, and transcript profiling of ERF family genes during development and abiotic stress treatments in tomato. Mol. Genet. Geno. 284, 455–475. doi: 10.1007/s00438-010-0580-1
Sharoni, A. M., Nuruzzaman, M., Satoh, K., Shimizu, T., Kondoh, H., Sasaya, T., et al. (2011). Gene structures, classification and expression models of the AP2/EREBP transcription factor family in rice. Plant Cell Physiol. 52, 344–360. doi: 10.1093/pcp/pcq196
Shu, Y., Liu, Y., Zhang, J., Song, L., and Guo, C. (2016). Genome-wide analysis of the AP2/ERF superfamily genes and their responses to abiotic stress in Medicago truncatula. Front. Plant Sci. 6:1247. doi: 10.3389/fpls.2015.01247
Sohn, K. H., Lee, S. C., Jung, H. W., Hong, J. K., and Hwang, B. K. (2006). Expression and functional roles of the pepper pathogen-induced transcription factor rav1 in bacterial disease resistance, and drought and salt stress tolerance. Plant Mol. Biol. 61, 897–915. doi: 10.1007/s11103-006-0057-0
Song, X. M., Li, Y., and Hou, X. L. (2013). Genome-wide analysis of the AP2/ERF transcription factor superfamily in Chinese cabbage (Brassica rapa ssp. pekinensis). BMC Geno. 14:573. doi: 10.1186/1471-2164-14-573
Theocharis, A., Christophe, C., and Essaïd, A. B. (2012). Physiological and molecular changes in plants grown at low temperatures. Planta 235, 1091–1105.
Velivelli, S. L., Lojan, P., Cranenbrouck, S., de Boulois, H. D., Suarez, J. P., Declerck, S., et al. (2015). The induction of ethylene response factor 3 (ERF3) in potato as a result of co-inoculation with Pseudomonas sp. R41805 and rhizophagus irregularis MUCL 41833-a possible role in plant defense. Plant Signal. Behav. 10:e988076. doi: 10.4161/15592324.2014.988076
Wu, L. J., Zhang, Z. J., Zhang, H. W., Wang, X. C., and Huang, R. F. (2008). Transcriptional modulation of ethylene response factor protein JERF3 in the oxidative stress response enhances tolerance of tobacco seedlings to salt, drought, and freezing. Plant Physiol. 148, 1953–1963. doi: 10.1104/pp.108.126813
Xie, Z., Yang, C., Li, M., Zhang, Z., Wu, Y., Gu, L., et al. (2022). Nitric oxide crosstalk with phytohormone is involved in enhancing photosynthesis of Tetrastigma hemsleyanum for photovoltaic adaptation. Front. Plant Sci. 13:852956. doi: 10.3389/fpls.2022.852956
Xu, L., Feng, G. Y., Yang, Z. F., Xu, X. H., Huang, L. K., Yang, Q. C., et al. (2020). Genome-wide AP2/ERF gene family analysis reveals the classification, structure, expression profiles and potential function in orchardgrass (Dactylis glomerata). Mol. Biol. Rep. 47, 5225–5241. doi: 10.1007/s11033-020-05598-x
Xu, Z. S., Chen, M., Li, L. C., and Ma, Y. Z. (2011). Functions and application of the AP2/ERF transcription factor family in crop improvement. J. Int. Plant Biol. 53, 570–585. doi: 10.1111/j.1744-7909
Yamamoto, H., Fan, X. Y., Sugimoto, K., Fukao, Y., Peng, L. W., and Shikanai, T. (2016). Chlororespiratory reduction 9 is a novel factor required for formation of subcomplex a of the chloroplast NADH dehydrogenase-like complex. Plant Cell Physiol. 57, 2122–2132. doi: 10.1093/pcp/pcw130
Ye, K. Y., Li, H., Ding, Y. L., Shi, Y. T., Song, C. P., Gong, Z. Z., et al. (2019). brassinosteroid-insensitive2 negatively regulates the stability of transcription factor ICE1 in response to cold stress in Arabidopsis. Plant Cell 31, 2682–2696. doi: 10.1105/tpc.19.00058
Yu, C. W., Qiao, G. R., Qiu, W. M., Yu, D. B., Zhou, S. R., Shen, Y., et al. (2018). Molecular breeding of water lily: Engineering cold stress tolerance into tropical water lily. Hortic. Res. 5:73. doi: 10.1038/s41438-018-0086-2
Zhang, H. N., Pan, X. L., Liu, S. H., Lin, W. Q., Li, Y. H., and Zhang, X. M. (2021). Genome-wide analysis of AP2/ERF transcription factors in pineapple reveals functional divergence during flowering induction mediated by ethylene and floral organ development. Genomics 113, 474–489. doi: 10.1016/j.ygeno.2020.10.040
Zhang, J., Liao, J. Y., Ling, Q. Q., Xi, Y., and Qian, Y. X. (2022). Genome-wide identification and expression profiling analysis of maize AP2/ERF superfamily genes reveal essential roles in abiotic stress tolerance. BMC Geno. 23:125. doi: 10.1186/s12864-022-08345-7
Zhuang, J., Anyia, A., Vidmar, J., Xiong, A. S., and Jian, Z. (2011). Discovery and expression assessment of the ap2-like genes in hordeum vulgare. Acta Physiol. Plant. 33, 1639–1649.
Zhuang, J., Cai, B., Peng, R. H., Zhu, B., Jin, X. F., Xue, Y., et al. (2008). Genome-wide analysis of the AP2/ERF gene family in Populus trichocarpa. Biochem. Biophys. Res. Commun. 371, 468–474. doi: 10.1016/j.bbrc.2008.04.087
Keywords: Tetrastigma hemsleyanum Diels at Gilg, cold stress, expression profiling, AP2/ERF transcription factor (TFs), functional analysis (FA)
Citation: Xie Z, Yang C, Liu S, Li M, Gu L, Peng X and Zhang Z (2022) Identification of AP2/ERF transcription factors in Tetrastigma hemsleyanum revealed the specific roles of ERF46 under cold stress. Front. Plant Sci. 13:936602. doi: 10.3389/fpls.2022.936602
Received: 05 May 2022; Accepted: 18 July 2022;
Published: 09 August 2022.
Edited by:
Yongfeng Guo, Tobacco Research Institute (CAAS), ChinaReviewed by:
Charu Lata, CSIR-National Institute of Science Communication and Policy Research, IndiaLiang Chen, University of Chinese Academy of Sciences, China
Copyright © 2022 Xie, Yang, Liu, Li, Gu, Peng and Zhang. 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: Xin Peng, cGVuZ3hAbml0LnpqdS5lZHUuY24=; Zhongyi Zhang, enl6aGFuZ0BmYWZ1LmVkdS5jbg==