Next Article in Journal
Microbiome of Nodules and Roots of Soybean and Common Bean: Searching for Differences Associated with Contrasting Performances in Symbiotic Nitrogen Fixation
Next Article in Special Issue
Proteomic Analysis Revealed Different Molecular Mechanisms of Response to PEG Stress in Drought-Sensitive and Drought-Resistant Sorghums
Previous Article in Journal
Enhanced Bioactivity of Tailor-Made Glycolipid Enriched Manuka Honey
Previous Article in Special Issue
Genome–Wide Identification of the GRAS Family Genes in Melilotus albus and Expression Analysis under Various Tissues and Abiotic Stresses
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 19
10.3390/ijms231912023
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Identification of ERF Transcription Factor Family and Functional Analysis of the Drought Stress-Responsive Genes in Melilotus albus

by
Na Wei
,
Qingyan Zhai
,
Hang Li
,
Shuwen Zheng
,
Jiyu Zhang
* and
Wenxian Liu
*
State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Western China Technology Innovation Center for Grassland Industry, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou 730000, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(19), 12023; https://doi.org/10.3390/ijms231912023
Submission received: 22 July 2022 / Revised: 29 September 2022 / Accepted: 6 October 2022 / Published: 10 October 2022
(This article belongs to the Special Issue Drought-Stress Induced Physiological and Molecular Changes in Plants 2.0)
Browse Figures
Review Reports Versions Notes

Abstract

:
As an important forage legume with high values in feed and medicine, Melilotus albus has been widely cultivated. The AP2/ERF transcription factor has been shown to play an important regulatory role in plant drought resistance, but it has not been reported in the legume forage crop M. albus. To digger the genes of M. albus in response to drought stress, we identified and analyzed the ERF gene family of M. albus at the genome-wide level. A total of 100 MaERF genes containing a single AP2 domain sequence were identified in this study, named MaERF001 to MaERF100, and bioinformatics analysis was performed. Collinearity analysis indicated that segmental duplication may play a key role in the expansion of the M. albus ERF gene family. Cis-acting element predictions suggest that MaERF genes are involved in various hormonal responses and abiotic stresses. The expression patterns indicated that MaERFs responded to drought stress to varying degrees. Furthermore, four up-regulated ERFs (MaERF008, MaERF037, MaERF054 and MaERF058) under drought stress were overexpressed in yeast and indicated their biological functions to confer the tolerance to drought. This work will advance the understanding of the molecular mechanisms underlying the drought response in M. albus. Further study of the promising potential candidate genes identified in this study will provide a valuable resource as the next step in functional genomics studies and improve the possibility of improving drought tolerance in M. albus by transgenic approaches.

1. Introduction

Periodic drought is the primary limitation on plant growth and yield of crops in agricultural systems [1]. Drought-induced loss in crop yield probably exceeds losses from all other causes, since both the severity and duration of the drought stress are critical [2]. Under drought conditions, plants can resist the negative effects of harsh environments to a certain extent by mobilizing stress response genes and developing a variety of defense mechanisms at the molecular and physiological levels [3]. Transcription factors play important roles as master regulators in various biological processes and are considered to be excellent candidate genes for improving crops by genetic engineering [4]. The AP2/ERF (APETA-LA2, Ethylene Response Factor) superfamily is one of the largest transcription factor families in plants, accounting for approximately 9% of the total number of known plant transcription factor genes [5]. The AP2/ERF family of genes contains the AP2/ERF domain, which consists of 60–70 amino acids [6]. According to the sequence similarity and the number of domains, AP2/ERF family transcription factors can be divided into four subfamilies: ABI3/VP1 (RAV), AP2, ERF and Soloist [7]. Among them, the AP2 subfamily contains two AP2 conserved domains that are involved in the regulation of flower, ovule and seed development and the ability to maintain meristems [8]. The ERF subfamily contains an AP2/ERF domain, which functions in response to plant biotic and abiotic stresses and includes the ERF subfamily and the DREB subfamily. The RAV subfamily includes an AP2/ERF domain and a B3 domain with conserved DNA-binding domains present in other plant-specific transcription factors (TFs) [9]. In the past few decades, ERF family genes have attracted attention because overexpressing ERF genes in transgenic plants can improve abiotic stress tolerance [10].
The ERF (Ethylene Response Factor) transcription factor was first known as the ethylene response element binding protein, and according to the type of amino acids, 14 and 19 in the AP2/ERF domain, the ERF transcription factor family can be further divided into the ERF subfamily and DREB subfamily [11]. The ERF subfamily members usually bind to the GCC-box with AGCGCCC as the core sequence in ethylene-responsive genes, while DREB binds to the DRE cis-acting element with A/GCCGAC as the core sequence [12,13]. However, in recent years, more and more research has shown that ERF and DREB could simultaneously bind to GCC-box and DRE cis-acting elements [14,15,16], suggesting that the ERF and DREB subfamilies may have similar functions in different developmental stages of plants. In the past few decades, ERF family genes have attracted attention because overexpression of ERF genes could improve various abiotic stress tolerance in transgenic plants [10,17]. In Arabidopsis thaliana, it was reported that the AtERF1 gene responds to drought and high salt stress by combining with GCC-box and DRE elements, and overexpression of AtERF1 enhanced the tolerance to drought and salt in Arabidopsis [18]. In rice (Oryza sativa), overexpression of OsERF48 could enhance root growth and drought tolerance by regulating OsCMLI6, a calmodulin-like protein [19]. Similar reports also found that in legume, GmDREB2A and MtDREB2A could activate the expression of downstream stress-related genes and significantly enhance the tolerance to drought in soybean (Glycine max) and Medicago truncatula [20,21]. In addition, it has recently been reported that ERF transcription factors can also induce the accumulation of proline. For example, the OsERF71 gene in rice regulates the expression of proline synthesis genes under drought stress, resulting in the accumulation of proline. Under salt stress, overexpression VaERF3 of Vigna angularis increases the level of proline accumulation and improves the salt resistance in Arabidopsis [22].
M. albus is an important forage crop worldwide, and it has high economic and utilization value for the promotion of green agriculture, rural areas and sustainable agricultural production. Mainly distributed in northwest, north, and northeast China, it has the characteristics of strong vitality, wide adaptability, fast reproduction, high seed yield, and strong nitrogen fixation ability [23]. In particular, the M. albus (2n = 16) genome and transcriptome were sequenced, and the first high-quality genome was recently published [24]. The publication of this genome provides information on M. albus genotypes and permits genome-wide research on this species. However, there have not been reports of ERF genes involved in drought stress resistance in M. albus. Based on these data, we identified the MaERF gene family at the genome-wide level, analyzed their expression patterns under drought stress, and verified the function of MaERF genes under drought stress through heterologous expression in yeast. The genome-wide study of the ERF gene family in M. albus can help us understand the molecular mechanisms of its stress resistance and provide valuable clues for the functional characterization of MaERF genes in response to drought stress, which will provide genetic resources with greater resistance to stress through transgenic technology.

2. Results

2.1. Gene Identification and Chromosomal Localization of MaERFs

To identify ERF transcription factors in M. albus, previously identified ERF proteins in model plants were used as a query dataset for MaERF proteins in the M. albus genome. A total of 115 potential MaERF protein sequences were screened in the local database for M. albus. After deredundancy and domain identification, 100 MaERF protein sequences with one AP2 domain were retained for further phylogenetic and functional analysis (Table S1). All MaERF genes were named MaERF001–MaERF100 according to their locations on the eight chromosomes (Figure 1). The results showed that the 100 MaERF genes were unevenly distributed across chromosomes, with the largest number of genes on Chr1 and the smallest number of genes on Chr8.

2.2. Phylogenetic Evolution and Physicochemical Properties Analysis of MaERFs

To classify and study the evolutionary relationship of the ERF gene family in M. albus, the ERF gene classification of the model plant Arabidopsis was selected as a reference. A total of 122 AtERF genes and 100 MaERF genes with a single AP2/ERF domain were assigned to the ERF family, and based on the similarity of their encoded amino acid sequences, these genes were further classified into two subfamilies, including DREB classes (subgroups I, II, III, and IV) and ERF classes (subgroups V, VI, VI-L, VII, VIII, IX and X), based on phylogenetic analysis (Figure 2 and Figure S1).
In addition, the physicochemical properties of the 100 MaERF genes were also analyzed in this experiment, including protein length, MW (Molecular weight), pI (isoelectric point), GRAVY (Grand average of hydropathicity), and subcellular localization prediction (Table 1). The length of the 100 MaERF proteins identified ranged from 141 (MaERF037) to 506 (MaERF032) amino acids (aa). The molecular weights ranged from 14, 834 (MaERF002) to 56,148 (MaERF032) Da, and the pI ranged from 4.62 (MaERF039) to 10.19 (MaERF029). The grand average hydropathicity of all MaERF proteins was negative (<0), which shows that these MaERF proteins have good hydrophilicity. Moreover, the results of subcellular localization prediction indicated that most of the genes in this family are predicted to be located in the nucleus, followed by the cytoplasm (Table S2).

2.3. Multiple Sequence Alignment and Motif Distribution and Gene Structure Analysis of MaERFs

To study the sequence characteristics of ERF family genes, the MaERF protein sequences of M. albus were used for alignment. The results of multiple sequence alignment showed that the ERF gene family contains at least one highly conserved AP2/ERF DNA-binding domain, which consists of 57–70 amino acid residues. AP2 is characterized by two conserved regions called YRG and RAYD. Regarding the AP2 domain, the structure prediction revealed three β-sheets (β1, β2, and β3) and one α-helix region that shared significant amino acid similarity within two YRG and RAYD elements (Figure S2). In addition, the AP2 domains of ERF members were highly consistent among the groups; these domains are marked with red rectangles in the figure and include the β1-sheet G5 (glycine), β3-sheet G31 (glycine) and A39 (alanine) site in the α-helix. The conserved motifs in AP2/ERF family proteins in M. albus were investigated using MEME, revealing a total of 25 conserved motifs (designated motifs 1–25), as shown in Figure 3. Proteins in the same group contain similar motifs, while the motifs are divergent among different groups. For example, motif 6 is only present in the III group, motif 22 is only present in the VI group, and motif 16 is only present in the V and VII groups of the ERF subfamily. Similarly, motifs 1 and 2 are shared by members of the AP2 family. These results indicate that most motifs are distributed among specific groups, which is correlated with their functional divergence. Subsequently, we analyzed the gene structural characteristics of these 100 MaERF genes, and the results showed that 78% of genes in this family had no introns, and only 2% of them had UTR structural characteristics (Figure 4).

2.4. Duplication Events and Ka/Ks Analysis of the MaERFs

To further examine the evolutionary relationship of MaERF genes, segmented and tandemly repeated genome duplication events were investigated by MCScanX collinearity analysis (Figure 5). The results showed that these 19 pairs of genes were all located on different chromosomes, indicating that chromosomal segment duplication was the main method for the expansion of the ERF gene family in M. albus.
In addition, the Ka/Ks ratio is widely used to measure the genetic evolution and selection pressure of genes. The MaERF057/MaERF084 gene pair was too divergent due to gene sequence differences; thus, a numerical value could not be calculated. The Ka/Ks values of the remaining 18 gene pairs were all < 1, which indicated that purifying selection was the main force driving the evolution of the ERF gene family in M. albus (Table 2).

2.5. Expression Pattern Analysis in Response to Drought Stress of MaERFs

To verify whether MaERF genes are involved in response to drought stress, we investigated the expression patterns of the 100 MaERF genes identified in this experiment under drought stress conditions. The expression profiles of 100 MaERF genes from two plant parts, namely, drought stress of shoot (DSS) and drought stress of root (DSR) were determined. According to the hierarchical clustering result, these MaERF genes had different transcription levels under various drought treatment time points, and could be clustered into nine clades, named A to I (Figure 6). Most genes showed significant up-regulation profiles under drought treatment in roots and peaked at 3 h, such as those in A, D, G and H groups, suggesting that these genes can respond rapidly to drought stress. Genes in Group C showed significant up-regulation after drought stress and peaked at 24 h in shoots, indicating that the genes in this group were positively regulated. The expression levels of group E and group F were suppressed after drought stress both in roots and shoots, and the expression levels of group B genes were basically kept unchanged in two tissues throughout the treatment.

2.6. qRT-PCR Validation of MaERF Genes Expression

In order to verify the response of MaERF genes to drought stress, 11 genes (MaERF004, MaERF008, MaERF010, MaERF012, MaERF016, MaERF017, MaERF034, MaERF037, MaERF054, MaERF058, and MaERF085) which significantly up-regulated (more than 10-fold) under drought stress were selected and further tested by qRT-PCR (Figure 7). These 11 genes showed different expression profiles under drought stress in different tissues. In roots, almost all MaERF genes induced rapidly under drought stress within 3 h, especially for MaERF037 and MaERF054, whose expression reached peaks at 1 h and expression levels were 10 times higher than that of control. Compared to root, most MaERF genes peaked at 6 h or 24 h under drought stress in shoot.

2.7. Analysis of cis-Acting Elements and Protein-Protein Interaction Networks of MaERFs

The promoter is a DNA sequence that determines the correct and efficient transcription of genes in plants. To study the function of MaERF promoters, the promoters of the 11 genes were extracted from 2000 bp upstream of the start codons using TBtools. Subsequently, potential functional elements in the promoter regions of these genes were analyzed by querying the PlantCARE database. The promoters are mainly grouped into four categories: abiotic, biotic, light responsiveness, and plant growth and development (Figure 8). We found that all 11 MaERF genes contain cis-acting elements related to abiotic stress, such as ABRE (abscisic acid response element); AuxRR-core (auxin response element); LTR (low-temperature response element); MBS (drought stress-related elements) and other abiotic stress-related elements. Except for MaERF004 and MaERF016, the remaining genes contain an ABRE, and only MaERF016 and MaERF058 contain an MBS element. These binding sites are located upstream of coding sequences and can provide binding sites for transcription factors that respond to abiotic stresses. In addition, except for MaERF016, all other genes had biological stress-related elements, and the MaERF017 gene contained the largest number of light-response-related cis-elements. At the same time, not all 11 genes had elements related to growth and development.
Furthermore, the STRING softer was used to determine the functional and physical relationships of 11 MaERF proteins predicted to be stress related through an M. truncatula association model (Figure S3). Using the model crop M. truncatula as the reference genome, the protein interaction prediction of these 11 orthologs of MaERF proteins was performed. The predicted result shows that an NAC-like transcription factor (AET04823) and a gene of unknown function (AES64073) interacted with ERF1A in M. truncatula. A previous study found that overexpression of the NAC-like transcription factor confers tolerance to drought and high salt stress through an ABA-independent signaling pathway [25]. Whether the interaction between ERF1A and NAC-like transcription factors can improve plant stress resistance remains to be verified, which provides a reference for predicting the potential regulatory roles of MaERF proteins in M. albus.

2.8. Transgenic Yeast to Analysis Drought Stresses

Considering both yeast and plants belong to eukaryotes, and yeast has been widely used as a model research system to study the functions of transcription factors under multiple abiotic stresses [10], in this study, we chose the yeast system to test the potential functions of four MaERFs in response to drought stress. The results showed that the transgenic yeast cells and untransformed empty yeast cells grew well under the control treatment. After simulating drought stress with 30% PEG-6000, the growth of transgenic yeast cells was significantly better than that of untransformed empty yeast cells (Figure 9). This indicated that heterologous expression of MaERF017, MaERF037, MaERF054, and MaERF058 improved the tolerance of yeast to osmotic stress. It was preliminarily proven that these four drought-tolerant candidate genes were possibly involved in the drought-tolerant response of M. albus.

3. Discussion

AP2/ERF transcription factors are an important family of transcription factors in all plant species that can participate in the regulation of various growth and development processes in plants and the responses to various stresses [26]. At present, many studies have confirmed that the ERF family plays an important role in the response to plant abiotic stresses [11]. However, there have been no studies on the molecular function of ERF TFs in response to abiotic stresses in M. albus. Therefore, this study focused on systematically exploring the potential molecular functions of MaERF genes under drought stress.
First, we identified a total of 100 MaERF genes in M. albus (Figure S1). This result is similar to the number of genes in M. truncatula (n = 107) [27] and soybean (n=98) [28] but slightly less than that in rice (n = 139) and Arabidopsis (n = 122). In total, 100 ERF proteins were grouped into 11 different subfamilies through phylogenetic analysis. The phylogenetic analysis results of most of the ERF proteins were consistent with a previous report in Arabidopsis [29], which showed that the results of the phylogenetic analysis in the present study have a high degree of credibility (Figure S2). The nucleotide sequence lengths, molecular weights and theoretical isoelectric points of the MaERF genes varied. It is worth noting that all MaERF genes have good hydrophilicity, and most genes are predicted to be localized to the nucleus (Table S2), which is consistent with the fact that transcription factors regulate the transcription of downstream genes in the nucleus to adapt to changes in the external environment. In domain research, it was found that the MaERF domain contains two important regions, namely, YRG and RAYD (Figure S2). The YRG region contains 20 amino acids in the β1 sheet of the AP2 domain and has been shown to play a key role in establishing direct contact with DNA molecules. In contrast, the RAYD region contains approximately 40 amino acids in the α-helical region and is involved in protein‒protein interactions [30]. In addition, it has also been reported that the hydrophobic surface of the α-helix of the RAYD region interacts with the major groove of DNA [31]. Moreover, we found that gene members belonging to the same subfamily had similar motifs, and found that the motif structures of MaERF proteins in each group were highly conserved and all amino acid sequences of MaERF proteins contained a conserved AP2 motif (Figure 3). The above results are consistent with the reports of Arabidopsis, rice and soybean [29,32].
Gene structure analysis showed that the MaERF gene family had fewer introns (Figure 4). The number and distribution of introns may be linked to plant evolution, and the intron number of ERF subfamily genes has likely decreased during plant evolution [33,34], which was also confirmed in Arabidopsis and M. truncatula [27,29]. In addition, we found that gene duplication occurred in the MaERF gene family, and a total of 19 duplication events were identified. All repetitive genes show segmental duplication, indicating that segmental duplication is important to the process of gene duplication in the MaERF gene family. Meanwhile, ERF gene duplication has been found in many species, including Rosa chinensis [35] and Gossypium barbadense [36] (Figure 5). Interestingly, 18 of the MaERF genes had Ka/Ks ratios < 1, indicating that the MaERF gene family has undergone purifying rather than positive selection (Figure 6) [37].
A cis-regulatory element is a specific motif located in the promoter region of a gene that acts as a binding site for a gene and can play an important role in the stress response by regulating the transcription of downstream genes [38]. In addition, some phytohormones (abscisic acid, salicylic acid, jasmonic acid and ethylene) are also involved in regulating the adaptive response of plants to abiotic stresses [39]. After the analysis of cis-elements in the promoter region of the candidate genes screened in the previous stage [40], it was found that the promoter elements involved in stress response and hormone response, including ABRE, LTR, MBS and TC-rich elements, exist in the promoter region of MaERF genes. These elements are involved in the ABA response and low temperature, drought and other biotic or abiotic stress defense processes (Figure 7). The results indicate that the MaERF gene promoter may play an important role in the transcriptional activation of stress response-related genes and the process of the plant stress response [41].
Studies have shown that ERF transcription factors can be potential candidates for crop improvement as they are key regulators in different plant developmental processes and responses to various stresses [42]. Nevertheless, the function of the MaERF gene in M. albus is still unclear. Therefore, it is necessary to analyze the transcriptional regulation of ERF transcription factors in M. albus, in order to use them to improve the quality and abiotic stress tolerance of M. albus. Here, we first systematically analyzed the expression patterns of these genes under drought stress and determined the potential functions of MaERF genes under drought stress. We observed a significant up-regulated expression trend of MaERF genes under drought stress, implying that they may be broadly related to the regulation of M. albus in response to drought (Figure 8). In addition, the expression levels of MaERF017, MaERF037, MaERF054, and MaERF058 genes were more fold-up-regulated, suggesting that these genes may play roles in key processes. In particular, MaERF058 showed a continuous up-regulation with a 20-fold increase. Interestingly, MtERF1A (the homologous gene of MaERF058 in M. truncatula) mediates resistance to a subset of root pathogens, playing an important function in M. truncatula responding to biotic stress [43]. As previously reported, overexpression of the BrERF4 gene in Arabidopsis can improve drought and salt tolerance and overexpression of SodERF3 in tobacco significantly improved the drought and salt tolerance of transgenic tobacco [41,44]. In this study, MaERF017, MaERF037, MaERF054, and MaERF058 were validated in transgenic yeast cells, and the results showed that transgenic yeast cells grew better than empty vector yeast cells after drought resistance treatment (Figure 9). It is speculated that the drought resistance of transgenic yeast is related to the combination with the DRE/CRT (dehydration response element/C-repeat) cis-acting element (TAC-CGACAT) that can regulate the expression of a series of downstream stress response genes and functional proteins [42]. However, the specific reasons for this up-regulation of these genes are still to be verified, and whether these MaERF genes can function as negative regulators in drought stress in M. albus still needs to be further investigated. The results of this study will lay a foundation for analyzing the drought resistance function of the MaERF gene family in M. albus and creating new germplasm of M. albus with high drought resistance through genetic engineering technology.

4. Materials and Methods

4.1. Identification of the ERF Gene Family in M. albus

To identify ERF gene family members in M. albus, this study first obtained the ERF protein sequences of 122 Arabidopsis, 139 O. sativa and 98 M. truncatula proteins as a query [27,29]. The MaERF protein and gene sequences were obtained by searching the genome sequence of M. albus with the local BLASTP program, and the E value was 10−5 [45]. The Pfam program online tool (http://pfam.xfam.org/search#tabview=tab1, accessed on 9 November 2021) and the NCBI-CD Search (https://www. ncbi. nlm. nih. gov/Structure/cdd/wrpsb.Cgi, accessed on 12 November 2021) were used to determine that these sequences contained only one AP2 domain [46]. The CD-Hit tool (http://weizhongli-lab.org/cd-hit/, accessed on 12 November 2021) was used to remove redundant sequences in the remaining ERF protein sequences [40]. The molecular weight (MW), grand average hydropathicity (GRAVY), and isoelectric point (pI) of M. albus ERF gene family members were detected using ExPASy (http://web.expasy.org/protparam/), accessed on 18 November 2021 [47].

4.2. Chromosomal Location and Phylogenetic Analysis of MaERF Genes

Through the M. albus genomic gff3 data, the corresponding position of the MaERF gene on the chromosome was obtained. The chromosome location map was drawn using R 3.6.3 software. Then, we named the genes according to their position on the chromosome. Phylogenetic analysis of protein sequences was performed using MEGA-X software [48]. The phylogenetic tree was constructed based on the neighbor-joining (NJ) method, the bootstrapping multiple calculation method was set to 1000, and other parameters were set to the default [49].

4.3. AP2/ERF Domain Analysis of MaERF Genes

To analyze the sequence features of the AP2 domain in the MaERF protein, we performed a multiple sequence alignment of the ERF protein sequence in M. albus using DNAMAN [50]. The sequence logos for AP2 repeats were generated using WebLogo online tool (http://weblogo.threeplusone.com/, accessed on 12 December 2021) with the default settings [51].

4.4. Conserved Motif and Gene Structure Analysis of MaERF Genes

We performed conserved motif analysis using the online tool MEME Suit 5.3.3 (https://meme-suite.org/meme/tools/meme, accessed on 18 December 2021) [52]. The conserved motifs in full-length ERF proteins were identified using MEME, with the maximum number of motifs as 25 [32]. Gene structure analysis of the ERF gene family members in M. albus was performed using the GSDS 2.0 (gene structure display server) (http://gsds.gao-lab.Org/, accessed on 22 December 2021) online website [53].

4.5. Gene Duplication Events and Ka/Ks Analysis of MaERF Genes

We used the MCScanX tool to assess the complete genome sequences and genome annotation files of M. albus [54]. TBtools software was used to visualize the obtained results [55]. Advanced Circos was used to visualize the collinearity of homologous genes based on homology and their positions in the genome. In addition, nonsynonymous (Ka) and synonymous (Ks) substitution rates were calculated to explore the mechanism of gene divergence after duplication by using TBtools.

4.6. Expression Analysis of MaERF Genes under Drought Stresses

To analyze the expression profiles of the MaERFs in shoot stress (DSS) and the root stress (DRS) under drought stress, the expression data were obtained as reported previously [56], and used for heatmap creation to display the gene expression profiles with TBtools software.

4.7. Analysis of Cis-Acting Elements and Protein-Protein Interaction Networks of MaERF Genes

According to the results of the expression pattern analysis, the cis-acting elements of 11 MaERF genes with high differential expression were analyzed. Subsequently, TBtools was used to extract the promoter sequences 2000 bp upstream of the start codons of these genes. Next, potential cis-acting elements in the promoter regions of these genes were analyzed using the PlantCARE database (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 29 December 2021) and graphs were drawn using R 3.6.3 software [57]. Using the model plant M. truncatula as a reference, the network structure interactions of 11 MaERF proteins were predicted, and the MaERF proteins network interactions were analyzed based on STRING software, where the high confidence was set to 0.4 and the max number was set to 10 [58].

4.8. Plant Cultivation and Drought Treatments

The variety “Ma46” of M. albus was chosen as the study material. First, we selected seeds with full grains and good shape for sterilization and germination. When the seeds germinated and revealed white buds, we selected seedlings with good growth and uniformity, transferred them to a hydroponic box containing 1/2 MS (pH = 5.8), and regularly replaced the culture medium in the hydroponic box [59]. The guaranteed incubation conditions were as follows: 16 h light/8 h dark cycle, 80% relative humidity, and 22 °C. Previous studies have indicated that H2O2, proline, malondialdehyde (MDA) and soluble sugar can be used as indicators to evaluate drought resistance in plants [60]. Our previous study found that under drought stress, the content of H2O2 peaked up to 24 h and the content of malondialdehyde (MDA), proline, and soluble sugar increased continuously up to 48 h, while the soluble protein content remained basically unchanged for 24 h and then decreased significantly [24]. So, in this study, we used 20% PEG6000 to simulate the drought stress and set different time points within 24 h [56]. To avoid the influence of circadian rhythm on plant growth and plant gene expression, the root and shoot tissue were uniformly sampled after 24 h, flash-frozen in liquid nitrogen and stored at −80 °C.

4.9. Plant RNA Extraction and Quantitative Real-Time PCR Analysis

Total RNA was extracted from M. albus under different treatment times. A cDNA synthesis kit (FastKing gDNA Dispelling RT SuperMix, Tiangen Biotech, Bei Jing) was used for reverse transcription, and first-strand cDNA was obtained after genomic DNA was removed. The concentration of cDNA was detected by a NanoDrop 2000 UV spectrophotometer (ND-8000, Xi’an) and then uniformly diluted to 100 ng·µL−1 for qRT‒PCR. Gene-specific primers for qRT‒PCR analysis were designed with NCBI Primer-BLAST (Table S3). qRT‒PCR analysis was performed using a 7500 fast real-time PCR system (Applied Biosystems, USA), and the experiment was performed in 3 technical replicates. The qRT‒PCR system included 5 µL 2SG Fast qPCR Master Mix (B639271-0005, Shanghai), 0.5 µL forward and reverse primers, 1 µL DNF Buffer, 1 µL cDNA and 2 µL ddH2O. The qRT‒PCR conditions were as follows: reaction denaturation (95 °C for 30 s), and 40 cycles of 95 °C for 5 s and 60 °C for 30 s. β-Tubulin was used as an internal reference gene for the ΔCt method to calculate the relative fold difference, and the gene expression level was calculated using the FC = 2−ΔΔCt method [61].

4.10. Yeast Heterologous Expression and Functional Validation of MaERF Genes

Total RNA was extracted from the M. albus seedlings, and the first-strand cDNA was obtained after genomic DNA was removed. After the target gene was cloned and sequenced, the amplified fragment was ligated into the pYES2 (Invitrogen, Carlsbad, CA, USA) expression vector by seamless cloning; the primer sequences are shown in Table S4. An empty pYES2 control plasmid and five expression vectors were transformed into the yeast strain INVSc1 by using a lithium acetate procedure according to the pYES2 vector kit instructions [37]. Following a previous method [62], we performed the transformation of yeast strains, and the transformants were incubated in a synthetic complete (SC) medium devoid of uracil with 2% (w/v) glucose at 30 °C for 36 h at 220 rpm. The drought stress evaluation experiment was carried out in the SC-Ura medium. We added 30% PEG-6000 stress solution to the yeast cells and added water to another group as a blank control. Afterwards, we serially diluted the 36 h-stressed culture and spotted it on SC-Ura agar plates. The spotted SC-Ura medium was placed upside down at 30 °C and observed after 48 h.

5. Conclusions

This study explored the potential molecular function of MaERF genes in M. albus. We identified the MaERF gene family at the genome-wide level and studied the basic gene characteristics of the family. In addition, the gene structure, motif composition, phylogenetic relationship, chromosomal location, collinearity analysis, cis-regulatory elements and other characteristics of the ERF gene family were studied. At the same time, we studied the expression changes of this gene family under drought stress at the molecular level and verified the function of MaERF genes by heterologous expression in transgenic yeast. These results will enrich our knowledge of the MaERF gene family and lay a foundation for further exploration of the function of MaERF genes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms231912023/s1.

Author Contributions

The experiments were conceived and designed by W.L and J Z, Q.Z. performed seedling maintenance. H.L. and S.Z. performed qRT-PCR experiments. N.W. performed data analysis and data visualization. N.W. wrote the first draft. W.L. and N.W. revised the manuscript. W.L. and J.Z. received financial support for the project leading to this publication. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (32061143035), the Natural Science Foundation of Gansu Province (20JR10RA622), and the Fundamental Research Funds for the Central Universities (lzujbky-2021-ct21, lzujbky-2021-it05, lzujbky-2021-ct14).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Boyer, J.S. Plant productivity and environment. Science 1982, 218, 443–448. [Google Scholar] [CrossRef] [PubMed]
  2. Farooq, M.; Wahid, A.; Kobayashi, N.; Fujita, D.; Basra, S. Plant drought stress: Effects, mechanisms and management. In Sustainable Agriculture; Springer: Berlin, Germany, 2009; pp. 153–188. [Google Scholar]
  3. Nakashima, K.; Ito, Y.; Yamaguchi-Shinozaki, K. Transcriptional regulatory networks in response to abiotic stresses in Arabidopsis and grasses. Plant Physiol. 2009, 149, 88–95. [Google Scholar] [CrossRef] [Green Version]
  4. Almeida, D.M.; Gregorio, G.B.; Oliveira, M.M.; Saibo, N.J. Five novel transcription factors as potential regulators of OsNHX1 gene expression in a salt tolerant rice genotype. Plant Mol. Biol. 2017, 93, 61–77. [Google Scholar] [CrossRef]
  5. Mizoi, J.; Shinozaki, K.; Yamaguchi-Shinozaki, K. AP2/ERF family transcription factors in plant abiotic stress responses. Biochim. Et Biophys. Acta (BBA)-Gene Regul. Mech. 2012, 1819, 86–96. [Google Scholar] [CrossRef]
  6. Riechmann, J.L.; Heard, J.; Martin, G.; Reuber, L.; Jiang, C.-Z.; Keddie, J.; Adam, L.; Pineda, O.; Ratcliffe, O.; Samaha, R. Arabidopsis transcription factors: Genome-wide comparative analysis among eukaryotes. Science 2000, 290, 2105–2110. [Google Scholar] [CrossRef] [PubMed]
  7. Xu, Z.S.; Chen, M.; Li, L.C.; Ma, Y.Z. Functions and application of the AP2/ERF transcription factor family in crop improvement F. J. Integr. Plant Biol. 2011, 53, 570–585. [Google Scholar] [CrossRef] [PubMed]
  8. Trupiano, D.; Yordanov, Y.; Regan, S.; Meilan, R.; Tschaplinski, T.; Scippa, G.S.; Busov, V. Identification, characterization of an AP2/ERF transcription factor that promotes adventitious, lateral root formation in Populus. Planta 2013, 238, 271–282. [Google Scholar] [CrossRef] [PubMed]
  9. Jin, X.; Yin, X.; Ndayambaza, B.; Zhang, Z.; Min, X.; Lin, X.; Wang, Y.; Liu, W. Genome-wide identification and expression profiling of the ERF gene family in Medicago sativa L. under various abiotic stresses. DNA Cell Biol. 2019, 38, 1056–1068. [Google Scholar] [CrossRef] [PubMed]
  10. Li, X.; Zhang, D.; Li, H.; Wang, Y.; Zhang, Y.; Wood, A.J. EsDREB2B, a novel truncated DREB2-type transcription factor in the desert legume Eremosparton songoricum, enhances tolerance to multiple abiotic stresses in yeast and transgenic tobacco. BMC Plant Biol. 2014, 14, 44. [Google Scholar] [CrossRef] [Green Version]
  11. Ohme-Takagi, M.; Shinshi, H. Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. Plant Cell 1995, 7, 173–182. [Google Scholar]
  12. Hao, D.; Ohme-Takagi, M.; Sarai, A. Unique mode of GCC box recognition by the DNA-binding domain of ethylene-responsive element-binding factor (ERF domain) in plant. J. Biol. Chem. 1998, 273, 26857–26861. [Google Scholar] [CrossRef]
  13. Hao, D.; Yamasaki, K.; Sarai, A.; Ohme-Takagi, M. Determinants in the sequence specific binding of two plant transcription factors, CBF1 and NtERF2, to the DRE and GCC motifs. Biochemistry 2002, 41, 4202–4208. [Google Scholar] [CrossRef] [PubMed]
  14. Lin, Z.; Hong, Y.; Yin, M.; Li, C.; Zhang, K.; Grierson, D. A tomato HD-Zip homeobox protein, LeHB-1, plays an important role in floral organogenesis and ripening. Plant J. 2008, 55, 301–310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Lee, S.-Y.; Hwang, E.Y.; Seok, H.-Y.; Tarte, V.N.; Jeong, M.S.; Jang, S.B.; Moon, Y.-H. Arabidopsis AtERF71/HRE2 functions as transcriptional activator via cis-acting GCC box or DRE/CRT element and is involved in root development through regulation of root cell expansion. Plant Cell Rep. 2015, 34, 223–231. [Google Scholar] [CrossRef] [PubMed]
  16. Klay, I.; Gouia, S.; Liu, M.; Mila, I.; Khoudi, H.; Bernadac, A.; Bouzayen, M.; Pirrello, J. Ethylene Response Factors (ERF) are differentially regulated by different abiotic stress types in tomato plants. Plant Sci. 2018, 274, 137–145. [Google Scholar] [CrossRef]
  17. Figueroa-Yañez, L.; Pereira-Santana, A.; Arroyo-Herrera, A.; Rodriguez-Corona, U.; Sanchez-Teyer, F.; Espadas-Alcocer, J.; Espadas-Gil, F.; Barredo-Pool, F.; Castaño, E.; Rodriguez-Zapata, L.C. RAP2. 4a is transported through the phloem to regulate cold and heat tolerance in papaya tree (Carica papaya cv. Maradol): Implications for protection against abiotic stress. PLoS ONE 2016, 11, e0165030. [Google Scholar] [CrossRef] [Green Version]
  18. Cheng, M.-C.; Liao, P.-M.; Kuo, W.-W.; Lin, T.-P. The Arabidopsis Ethylene Response Factor1 regulates abiotic stress-responsive gene expression by binding to different cis-acting elements in response to different stress signals. Plant Physiol. 2013, 162, 1566–1582. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Jung, H.; Chung, P.J.; Park, S.H.; Redillas, M.C.F.R.; Kim, Y.S.; Suh, J.W.; Kim, J.K. Overexpression of Os ERF 48 causes regulation of Os CML 16, a calmodulin-like protein gene that enhances root growth and drought tolerance. Plant Biotechnol. J. 2017, 15, 1295–1308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Sakuma, Y.; Maruyama, K.; Osakabe, Y.; Qin, F.; Seki, M.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Functional analysis of an Arabidopsis transcription factor, DREB2A, involved in drought-responsive gene expression. Plant Cell 2006, 18, 1292–1309. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Chen, J.-R.; Lü, J.-J.; Wang, T.-X.; Chen, S.-Y.; Wang, H.-F. Activation of a DRE-binding transcription factor from Medicago truncatula by deleting a Ser/Thr-rich region. Vitr. Cell. Dev. Biol. -Plant 2009, 45, 1–11. [Google Scholar] [CrossRef]
  22. Li, W.-Y.; Wang, C.; Shi, H.-H.; Wang, B.; Wang, J.-X.; Liu, Y.-S.; Ma, J.-Y.; Tian, S.-Y.; Zhang, Y.-W. Genome-wide analysis of ethylene-response factor family in adzuki bean and functional determination of VaERF3 under saline-alkaline stress. Plant Physiol. Biochem. 2020, 147, 215–222. [Google Scholar] [CrossRef] [PubMed]
  23. Luo, K.; Jahufer, M.; Wu, F.; Di, H.; Zhang, D.; Meng, X.; Zhang, J.; Wang, Y. Genotypic variation in a breeding population of yellow sweet clover (Melilotus officinalis). Front. Plant Sci. 2016, 7, 972. [Google Scholar] [CrossRef] [PubMed]
  24. Wu, F.; Duan, Z.; Xu, P.; Yan, Q.; Meng, M.; Cao, M.; Jones, C.S.; Zong, X.; Zhou, P.; Wang, Y. Genome and systems biology of Melilotus albus provides insights into coumarins biosynthesis. Plant Biotechnol. J. 2022, 20, 592. [Google Scholar] [CrossRef]
  25. Xie, L.-N.; Chen, M.; Min, D.-H.; Feng, L.; Xu, Z.-S.; Zhou, Y.-B.; Xu, D.-B.; Li, L.-C.; Ma, Y.-Z.; Zhang, X.-H. The NAC-like transcription factor SiNAC110 in foxtail millet (Setaria italica L.) confers tolerance to drought and high salt stress through an ABA independent signaling pathway. J. Integr. Agric. 2017, 16, 559–571. [Google Scholar] [CrossRef]
  26. Li, M.-Y.; Xu, Z.-S.; Huang, Y.; Tian, C.; Wang, F.; Xiong, A.-S. Genome-wide analysis of AP2/ERF transcription factors in carrot (Daucus carota L.) reveals evolution and expression profiles under abiotic stress. Mol. Genet. Genom. 2015, 290, 2049–2061. [Google Scholar] [CrossRef] [PubMed]
  27. Shu, Y.; Liu, Y.; Zhang, J.; Song, L.; Guo, C. Genome-wide analysis of the AP2/ERF superfamily genes and their responses to abiotic stress in Medicago truncatula. Front. Plant Sci. 2016, 6, 1247. [Google Scholar] [CrossRef] [Green Version]
  28. Yu, Y.; Duan, X.; Ding, X.; Chen, C.; Zhu, D.; Yin, K.; Cao, L.; Song, X.; Zhu, P.; Li, Q. A novel AP2/ERF family transcription factor from Glycine soja, GsERF71, is a DNA binding protein that positively regulates alkaline stress tolerance in Arabidopsis. Plant Mol. Biol. 2017, 94, 509–530. [Google Scholar] [CrossRef]
  29. Nakano, T.; Suzuki, K.; Fujimura, T.; Shinshi, H. Genome-wide analysis of the ERF gene family in Arabidopsis and rice. Plant Physiol. 2006, 140, 411–432. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Faraji, S.; Filiz, E.; Kazemitabar, S.K.; Vannozzi, A.; Palumbo, F.; Barcaccia, G.; Heidari, P. The AP2/ERF gene family in Triticum durum: Genome-wide identification and expression analysis under drought and salinity stresses. Genes 2020, 11, 1464. [Google Scholar] [CrossRef]
  31. Okamuro, J.K.; Caster, B.; Villarroel, R.; Van Montagu, M.; Jofuku, K.D. The AP2 domain of APETALA2 defines a large new family of DNA binding proteins in Arabidopsis. Proc. Natl. Acad. Sci. USA 1997, 94, 7076–7081. [Google Scholar] [CrossRef] [Green Version]
  32. Zhang, G.; Chen, M.; Chen, X.; Xu, Z.; Guan, S.; Li, L.-C.; Li, A.; Guo, J.; Mao, L.; Ma, Y. Phylogeny, gene structures, and expression patterns of the ERF gene family in soybean (Glycine max L.). J. Exp. Bot. 2008, 59, 4095–4107. [Google Scholar] [CrossRef]
  33. Tang, Y.; Qin, S.; Guo, Y.; Chen, Y.; Wu, P.; Chen, Y.; Li, M.; Jiang, H.; Wu, G. Genome-wide analysis of the AP2/ERF gene family in physic nut and overexpression of the JcERF011 gene in rice increased its sensitivity to salinity stress. PLoS ONE 2016, 11, e0150879. [Google Scholar] [CrossRef]
  34. Zhang, Z.; Li, X. Genome-wide identification of AP2/ERF superfamily genes and their expression during fruit ripening of Chinese jujube. Sci. Rep. 2018, 8, 15612. [Google Scholar] [CrossRef] [Green Version]
  35. Li, D.; Liu, X.; Shu, L.; Zhang, H.; Zhang, S.; Song, Y.; Zhang, Z. Global analysis of the AP2/ERF gene family in rose (Rosa chinensis) genome unveils the role of RcERF099 in Botrytis resistance. BMC Plant Biol. 2020, 20, 533. [Google Scholar] [CrossRef]
  36. Zafar, M.M.; Rehman, A.; Razzaq, A.; Parvaiz, A.; Mustafa, G.; Sharif, F.; Mo, H.; Youlu, Y.; Shakeel, A.; Ren, M. Genome-wide characterization and expression analysis of Erf gene family in cotton. BMC Plant Biol. 2022, 22, 134. [Google Scholar] [CrossRef]
  37. Wang, Q.; Wei, N.; Jin, X.; Min, X.; Ma, Y.; Liu, W. Molecular characterization of the COPT/Ctr-type copper transporter family under heavy metal stress in alfalfa. Int. J. Biol. Macromol. 2021, 181, 644–652. [Google Scholar] [CrossRef]
  38. Nakashima, K.; Yamaguchi-Shinozaki, K.; Shinozaki, K. The transcriptional regulatory network in the drought response and its crosstalk in abiotic stress responses including drought, cold, and heat. Front. Plant Sci. 2014, 5, 170. [Google Scholar] [CrossRef] [Green Version]
  39. Santner, A.; Estelle, M. Recent advances and emerging trends in plant hormone signalling. Nature 2009, 459, 1071–1078. [Google Scholar] [CrossRef]
  40. Zhou, Q.; Jia, C.; Ma, W.; Cui, Y.; Jin, X.; Luo, D.; Min, X.; Liu, Z. MYB transcription factors in alfalfa (Medicago sativa): Genome-wide identification and expression analysis under abiotic stresses. PeerJ 2019, 7, e7714. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Trujillo, L.; Sotolongo, M.; Menendez, C.; Ochogavia, M.; Coll, Y.; Hernandez, I.; Borras-Hidalgo, O.; Thomma, B.; Vera, P.; Hernandez, L. SodERF3, a novel sugarcane ethylene responsive factor (ERF), enhances salt and drought tolerance when overexpressed in tobacco plants. Plant Cell Physiol. 2008, 49, 512–525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Xie, Z.; Nolan, T.M.; Jiang, H.; Yin, Y. AP2/ERF transcription factor regulatory networks in hormone and abiotic stress responses in Arabidopsis. Front. Plant Sci. 2019, 10, 228. [Google Scholar] [CrossRef]
  43. Anderson, J.P.; Lichtenzveig, J.; Gleason, C.; Oliver, R.P.; Singh, K.B. The B-3 ethylene response factor MtERF1-1 mediates resistance to a subset of root pathogens in Medicago truncatula without adversely affecting symbiosis with rhizobia. Plant Physiol. 2010, 154, 861–873. [Google Scholar] [CrossRef] [Green Version]
  44. Seo, Y.J.; Park, J.-B.; Cho, Y.-J.; Jung, C.; Seo, H.S.; Park, S.-K.; Nahm, B.H.; Song, J.T. Overexpression of the ethylene-responsive factor gene BrERF4 from Brassica rapa increases tolerance to salt and drought in Arabidopsis plants. Mol. Cells 2010, 30, 271–277. [Google Scholar] [CrossRef]
  45. Min, X.; Jin, X.; Zhang, Z.; Wei, X.; Ndayambaza, B.; Wang, Y.; Liu, W. Genome-Wide Identification of NAC Transcription Factor Family and Functional Analysis of the Abiotic Stress-Responsive Genes in Medicago sativa L. J. Plant Growth Regul. 2020, 39, 324–337. [Google Scholar] [CrossRef]
  46. Zhang, Z.; Jin, X.; Liu, Z.; Zhang, J.; Liu, W. Genome-wide identification of FAD gene family and functional analysis of MsFAD3. 1 involved in the accumulation of α-linolenic acid in alfalfa. Crop Sci. 2021, 61, 566–579. [Google Scholar] [CrossRef]
  47. Gasteiger, E.; Gattiker, A.; Hoogland, C.; Ivanyi, I.; Appel, R.D.; Bairoch, A. ExPASy: The proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003, 31, 3784–3788. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547. [Google Scholar] [CrossRef] [PubMed]
  49. Zhao, S.-P.; Xu, Z.-S.; Zheng, W.-J.; Zhao, W.; Wang, Y.-X.; Yu, T.-F.; Chen, M.; Zhou, Y.-B.; Min, D.-H.; Ma, Y.-Z. Genome-wide analysis of the RAV family in soybean and functional identification of GmRAV-03 involvement in salt and drought stresses and exogenous ABA treatment. Front. Plant Sci. 2017, 8, 905. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Edwards, D.J.; Grodevant, N.W.; Lee, P.J.; Peralta, J.B. DNA-MAN: Dynamic Natural Attributes for Synthetic Military Forces. In Proceedings of the 2007 IEEE Systems and Information Engineering Design Symposium, Charlottesville, VA, USA, 27 April 2007; pp. 1–5. [Google Scholar]
  51. Crooks, G.E.; Hon, G.; Chandonia, J.-M.; Brenner, S.E. WebLogo: A sequence logo generator. Genome Res. 2004, 14, 1188–1190. [Google Scholar] [CrossRef] [Green Version]
  52. Bailey, T.L.; Johnson, J.; Grant, C.E.; Noble, W.S. The MEME suite. Nucleic Acids Res. 2015, 43, W39–W49. [Google Scholar] [CrossRef] [Green Version]
  53. Hu, B.; Jin, J.; Guo, A.-Y.; Zhang, H.; Luo, J.; Gao, G. GSDS 2.0: An upgraded gene feature visualization server. Bioinformatics 2015, 31, 1296–1297. [Google Scholar] [CrossRef] [PubMed]
  54. Wang, Y.; Tang, H.; DeBarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.-H.; Jin, H.; Marler, B.; Guo, H. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef] [PubMed]
  56. Duan, Z.; Yan, Q.; Wu, F.; Wang, Y.; Wang, S.; Zong, X.; Zhou, P.; Zhang, J. Genome-Wide Analysis of the UDP-Glycosyltransferase Family Reveals Its Roles in Coumarin Biosynthesis and Abiotic Stress in Melilotus albus. Int. J. Mol. Sci. 2021, 22, 10826. [Google Scholar] [CrossRef]
  57. Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van de Peer, Y.; Rouzé, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30, 325–327. [Google Scholar] [CrossRef]
  58. Franceschini, A.; Szklarczyk, D.; Frankild, S.; Kuhn, M.; Simonovic, M.; Roth, A.; Lin, J.; Minguez, P.; Bork, P.; von Mering, C.; et al. STRING v9.1: Protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013, 41, D808–D815. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Ma, Y.; Wei, N.; Wang, Q.; Liu, Z.; Liu, W. Genome-wide identification and characterization of the heavy metal ATPase (HMA) gene family in Medicago truncatula under copper stress. Int. J. Biol. Macromol. 2021, 193, 893–902. [Google Scholar] [CrossRef] [PubMed]
  60. Luo, D.; Liu, J.; Wu, Y.; Zhang, X.; Zhou, Q.; Fang, L.; Liu, Z. Nuclear transport factor 2-LIKE improves drought tolerance by modulating leaf water loss in alfalfa (Medicago sativa L.). Plant J. 2022. [Google Scholar] [CrossRef]
  61. Dong, R.; Dong, D.; Luo, D.; Zhou, Q.; Chai, X.; Zhang, J.; Xie, W.; Liu, W.; Dong, Y.; Wang, Y. Transcriptome analyses reveal candidate pod shattering-associated genes involved in the pod ventral sutures of common vetch (Vicia sativa L.). Front. Plant Sci. 2017, 8, 649. [Google Scholar] [CrossRef]
  62. Kawai, S.; Hashimoto, W.; Murata, K. Transformation of Saccharomyces cerevisiae and other fungi: Methods and possible underlying mechanism. Bioeng. Bugs 2010, 1, 395–403. [Google Scholar] [CrossRef]
Ijms 23 12023 g001 550
Figure 1. Chromosomal location distribution of 100 MaERF genes. The position of each gene shown on the graph is the average of the two ends of the gene, and the chromosome numbers are arranged in order. The scale on the left is in megabases (Mb). Chromosome colors from blue to red represent gene density from low to high gradually. The chromosome location map was generated using the R 3.6.3 software.
Figure 1. Chromosomal location distribution of 100 MaERF genes. The position of each gene shown on the graph is the average of the two ends of the gene, and the chromosome numbers are arranged in order. The scale on the left is in megabases (Mb). Chromosome colors from blue to red represent gene density from low to high gradually. The chromosome location map was generated using the R 3.6.3 software.
Ijms 23 12023 g001
Ijms 23 12023 g002 550
Figure 2. Phylogenetic relationship of 100 MaERF proteins in M. albus. An unrooted phylogenetic tree was constructed based on the multiple alignments of 100 MaERFs amino acid sequences with 1000 directed repeats. The classification of each subgroup is represented by a colored box. According to the classification method of Arabidopsis, 100 MaERF genes are divided into the ERF subfamily and DREB subfamily with 10 groups. Sequence-conserved motifs that are genetically similar in each group are identified by MEME.
Figure 2. Phylogenetic relationship of 100 MaERF proteins in M. albus. An unrooted phylogenetic tree was constructed based on the multiple alignments of 100 MaERFs amino acid sequences with 1000 directed repeats. The classification of each subgroup is represented by a colored box. According to the classification method of Arabidopsis, 100 MaERF genes are divided into the ERF subfamily and DREB subfamily with 10 groups. Sequence-conserved motifs that are genetically similar in each group are identified by MEME.
Ijms 23 12023 g002
Ijms 23 12023 g003 550
Figure 3. MEME was used to analyze protein motifs of MaERF proteins, the phylogenetic tree was constructed using MEGAX by the neighbor-joining method with 1000 bootstrap replicates.
Figure 3. MEME was used to analyze protein motifs of MaERF proteins, the phylogenetic tree was constructed using MEGAX by the neighbor-joining method with 1000 bootstrap replicates.
Ijms 23 12023 g003
Ijms 23 12023 g004 550
Figure 4. Structural analysis of the MaERF genes of M. albus. We clustered genes according to their phylogenetic relationships, with different colors representing different groups.
Figure 4. Structural analysis of the MaERF genes of M. albus. We clustered genes according to their phylogenetic relationships, with different colors representing different groups.
Ijms 23 12023 g004
Ijms 23 12023 g005 550
Figure 5. The collinearity analysis was performed on the MaERF genes, and only 18 pairs of gene names with collinearity were marked on the chromosome, and the green lines in the circle indicated the collinear relationship between genes.
Figure 5. The collinearity analysis was performed on the MaERF genes, and only 18 pairs of gene names with collinearity were marked on the chromosome, and the green lines in the circle indicated the collinear relationship between genes.
Ijms 23 12023 g005
Ijms 23 12023 g006 550
Figure 6. Heatmap representation of MaERF gene expression patterns under drought stress. Groups A to I exhibited nine expression patterns of MaERF genes. The bars at the right of the heatmap represent relative expression values. Root and shoot of drought stress expression levels at 0 h, 3 h and 24 h are shown.
Figure 6. Heatmap representation of MaERF gene expression patterns under drought stress. Groups A to I exhibited nine expression patterns of MaERF genes. The bars at the right of the heatmap represent relative expression values. Root and shoot of drought stress expression levels at 0 h, 3 h and 24 h are shown.
Ijms 23 12023 g006
Ijms 23 12023 g007 550
Figure 7. Expression patterns of MaERF genes in response to drought stress. The relative expression levels of the MaERF gene in M. albus root (A) and shoot (B) treated with PEG6000 (20%) were determined by qRT-PCR. Different letters indicate significant differences between different treatment time (p < 0.05). The name of the gene is written on the top of each bar diagram (error bars indicate the standard deviation from three replicates).
Figure 7. Expression patterns of MaERF genes in response to drought stress. The relative expression levels of the MaERF gene in M. albus root (A) and shoot (B) treated with PEG6000 (20%) were determined by qRT-PCR. Different letters indicate significant differences between different treatment time (p < 0.05). The name of the gene is written on the top of each bar diagram (error bars indicate the standard deviation from three replicates).
Ijms 23 12023 g007
Ijms 23 12023 g008 550
Figure 8. Analysis of cis-acting elements in the promoter in MaERF genes of the M. albus. A total of 20 cis-acting elements, including biotic, abiotic stress and phylogenetically related elements were counted, and their numbers were counted with different colors. The picture was generated using the R 3.6.3 software.
Figure 8. Analysis of cis-acting elements in the promoter in MaERF genes of the M. albus. A total of 20 cis-acting elements, including biotic, abiotic stress and phylogenetically related elements were counted, and their numbers were counted with different colors. The picture was generated using the R 3.6.3 software.
Ijms 23 12023 g008
Ijms 23 12023 g009 550
Figure 9. Expression of MsERF017, MsERF037, MsERF054 and MsERF058 in yeast-induced sensitivity to drought stress. The left panel represents the growth of yeast cells under normal conditions, and the right panel represents the growth of yeast cells under drought stress.
Figure 9. Expression of MsERF017, MsERF037, MsERF054 and MsERF058 in yeast-induced sensitivity to drought stress. The left panel represents the growth of yeast cells under normal conditions, and the right panel represents the growth of yeast cells under drought stress.
Ijms 23 12023 g009
Table
Table 1. Analysis of physicochemical properties of each group in MaERF genes of M. albus.
Table 1. Analysis of physicochemical properties of each group in MaERF genes of M. albus.
ERF GroupsNo of GenesLength (aa)Molecular Weight (Da)pIGRAVY
Max.Min.Avg.Max.Min.Avg.Max.Min.Avg.Max.Min.Avg.
ERF-I643229636248,78226,88840,6308.75.186.70−0.50−1.02−0.67
ERF-II1024515519926,88817,82222,1348.734.676.26−0.36−1.05−0.59
ERF-III2030817322434,90618,88624,9687.914.625.83−0.31−0.94−0.62
ERF-IV950619529556,14820,92732,7309.214.866.82−0.44−0.99−0.75
ERF-V830917622539,74734,57325,2999.325.77.22−0.61−0.81−0.70
ERF-VI635521627139,74824,19630,3279.194.775.72−0.62−0.73−0.68
ERF-VI-L238526533841,53532,02636,7806.685.156.005−0.66−0.79−0.72
ERF-VII437128732942,40229,22736,7425.824.755.09−0.47−0.96−0.73
ERF-VIII1335214424838,84616,46527,30910.194.887.84−0.21−1.07−0.68
ERF-IX1741714124246,53214,83427,1329.334.737.05−0.41−0.87−0.66
ERF-X541717130445,28819,03233,4638.696.086.93−0.71−1.08−0.81
All10050614127656,14814,83430,68310.194.626.49−0.21−1.08−0.69
Note: pI, isoelectric point, GRAVY, Grand average of hydropathicity.
Table
Table 2. Ka/ks analysis in the MaERF genes of M. albus.
Table 2. Ka/ks analysis in the MaERF genes of M. albus.
Sequance1Sequance2KaKsKa_KsEffectiveLenAverageS-SitesAverageN-Sites
MaERF10MaERF890.6421150.9660720.664665852200.8333333651.1666667
MaERF10MaERF1000.3604330.6763640.5328981017250.5833333766.4166667
MaERF13MaERF190.206721.4329350.144263519120.5398.5
MaERF14MaERF230.2337320.8631670.27078546598.16666667366.8333333
MaERF15MaERF180.3780741.8421930.20523720154.8333333565.1666667
MaERF20MaERF640.3388631.7417910.194549606133.1666667472.8333333
MaERF21MaERF290.3446931.5214770.226552561125.25435.75
MaERF26MaERF880.4225221.1917330.354545873193.5679.5
MaERF26MaERF970.2807050.8124290.345514990225.9166667764.0833333
MaERF27MaERF990.2357740.4066620.579779642144.4166667497.5833333
MaERF34MaERF790.2040080.8398050.242923780167.4166667612.5833333
MaERF42MaERF800.2466290.8716250.282953579135.3333333443.6666667
MaERF43MaERF480.4247682.1267460.199727558136.3333333421.6666667
MaERF48MaERF830.3728311.5837920.235404594138.1666667455.8333333
MaERF50MaERF650.2787260.6725960.414404846184.8333333661.1666667
MaERF56MaERF660.5047231.1585540.435649873194.5678.5
MaERF57MaERF840.570869NaNNaN567130.5833333436.4166667
MaERF61MaERF740.3137031.2640340.248176576130.6666667445.3333333
MaERF90MaERF190.3512292.1451540.163731537125.1666667411.8333333
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wei, N.; Zhai, Q.; Li, H.; Zheng, S.; Zhang, J.; Liu, W. Genome-Wide Identification of ERF Transcription Factor Family and Functional Analysis of the Drought Stress-Responsive Genes in Melilotus albus. Int. J. Mol. Sci. 2022, 23, 12023. https://doi.org/10.3390/ijms231912023

AMA Style

Wei N, Zhai Q, Li H, Zheng S, Zhang J, Liu W. Genome-Wide Identification of ERF Transcription Factor Family and Functional Analysis of the Drought Stress-Responsive Genes in Melilotus albus. International Journal of Molecular Sciences. 2022; 23(19):12023. https://doi.org/10.3390/ijms231912023

Chicago/Turabian Style

Wei, Na, Qingyan Zhai, Hang Li, Shuwen Zheng, Jiyu Zhang, and Wenxian Liu. 2022. "Genome-Wide Identification of ERF Transcription Factor Family and Functional Analysis of the Drought Stress-Responsive Genes in Melilotus albus" International Journal of Molecular Sciences 23, no. 19: 12023. https://doi.org/10.3390/ijms231912023

APA Style

Wei, N., Zhai, Q., Li, H., Zheng, S., Zhang, J., & Liu, W. (2022). Genome-Wide Identification of ERF Transcription Factor Family and Functional Analysis of the Drought Stress-Responsive Genes in Melilotus albus. International Journal of Molecular Sciences, 23(19), 12023. https://doi.org/10.3390/ijms231912023

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop