Next Article in Journal
Tea Bud Detection and 3D Pose Estimation in the Field with a Depth Camera Based on Improved YOLOv5 and the Optimal Pose-Vertices Search Method
Next Article in Special Issue
Genome-Wide Identification and Expression Analysis of the CesA/Csl Gene Superfamily in Alfalfa (Medicago sativa L.)
Previous Article in Journal
Return on Equity in Dairy Farms from Selected EU Countries: Assessment Based on the DuPont Model in Years 2004–2020
Previous Article in Special Issue
Study on Phyllosphere Microbial Community of Nettle Leaf during Different Seasons
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Agriculture
Volume 13
Issue 7
10.3390/agriculture13071404
agriculture-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

Leaf Transcription Factor Family Analysis of Halophyte Glaux maritima under Salt Stress

by
Rui Gu
1,
Zhiqiang Wan
2,
Fang Tang
1,
Fengling Shi
1,* and
Mengjiao Yan
3
1
Key Laboratory of Grassland Resources of Ministry of Education, College of Grassland, Resources and Environment, Inner Mongolia Agricultural University, Hohhot 010019, China
2
College of Geographical Science, Inner Mongolia Normal University, Hohhot 010020, China
3
Plant Protection Institute, Inner Mongolia Academy of Agricultural & Animal Husbandry Sciences, Hohhot 010031, China
*
Author to whom correspondence should be addressed.
Agriculture 2023, 13(7), 1404; https://doi.org/10.3390/agriculture13071404
Submission received: 25 May 2023 / Revised: 11 July 2023 / Accepted: 11 July 2023 / Published: 14 July 2023
(This article belongs to the Special Issue Forage Breeding and Cultivation)
Browse Figures
Review Reports Versions Notes

Abstract

:
The reduction of crop yield caused by soil salinization has become a global problem. Halophytes improve saline alkali soil, and the halophyte transcription factors that regulate salt stress are crucial for improving salt tolerance. In this study, 1466 transcription factors were identified by transcriptome sequencing analysis of Glaux maritima leaves after salt stress (0, 600, and 800 mM/L NaCl). Their genes were distributed across 57 transcription factor families. KEGG and GO analyses showed significant enrichment in 14 pathways, with a total of 54 functions annotated. Gene expression analysis showed 820 differentially expressed genes distributed in 11 transcription factor families, including ERF, bHLH, WRKY, and NAC, and 8 expression modules. KEGG analysis revealed four genes with significant positive regulation: ABF2 (Unigene0078257) in the ABA signaling pathway, EIN3 (Unigene0000457 and Unigene0012139), and EIL1 (Unigene0042139) involved in ethylene signal transduction, and two with negative regulation, MYC1/2 (Unigene0009899 and Unigene0027167) in the main regulator of Jasmonic acid signal transduction. Protein–protein interaction networks suggested ABF2 and MYC1/2 as important transcription factors regulating G. maritima salt tolerance. Overall, the salt-tolerant transcription factors discovered in this study provide genetic resources for plant salt tolerance inheritance, and lay a theoretical foundation for the study of the salt-tolerant molecular mechanism of the halophyte Glaux maritima.

1. Introduction

As one of the main abiotic stress factors, salinity plays an important role in limiting crop growth, development, and yield [1,2]. In order to adapt to adverse environmental stress, a series of genes in plants are synchronously expressed, and various mechanisms are evolved to regulate plant growth and development [3,4]. Searching for candidate genes and analyzing signal mechanisms in the abiotic stress response are crucial for cultivating stress-resistant and tolerant crops [5]. Transcription factors (TFs), as important regulators of plant growth, development, and stress responses, often participate in different and complex signal networks [6,7].
With recent in-depth and extensive research on the mechanism of the plant response to salt stress, a large number of involved TF families, including NAC (NAM, ATAF1/2, CUC2), ERF/AP2 (Ethylene Response factor/APETALA2), bZIP (basic leucine-zipper), bHLH (Basic Helix-Loop-Helix), MYB (v-myb avian mye1oblastosis viral oncogene homolog), and WRKY, have been discovered [8]. Although the main transcriptional regulator of salt stress is not yet clear, a large number of studies have shown that TFs related to the regulation of ABA, jasmonate, ethylene, and other hormones play important roles in the salt stress response. ABA mainly mediates transcriptional regulation through the AREB/ABF subfamily of bZIP TFs [9,10]. For example, ABF2 overexpression increases the resistance of plants to salt stress [11]. Overexpression of the ABA-responsive TFs dynamic influencer of gene expression (DIG) and DIG-like (DIL) leads to hypersensitivity to high-salt conditions or ABA [12]. MYC2, as the master TF of jasmonate signaling, plays a positive role in regulating plant salt tolerance [13,14]. EIN3, as a TF that mediates the core ethylene signaling pathway, improves plant salt tolerance through the expression of the DELLA protein after salt treatment [15,16]. EIN3 directly binds to a salt response gene and activates its expression through two downstream TFs, ethylene response factor 1 (ERF1) and ethylene and salt inducible 1 (ESE1), to improve plant salt tolerance [17,18]. In Arabidopsis thaliana, the nuclear calcium-binding gene AtbHLH28 (ATNIG1) is believed to play an important role in the regulation of plant salt stress signals by binding to the E-box motif (CANNTG) of salt-stress-related gene promoters [19]. Overexpression of bHLH106 in different plants increases tolerance to NaCl [4], and co-expression of bHLH17 and WRKY28 leads to resistance in A. thaliana to various abiotic stresses [20]. Salt has a significant inductive effect on wheat TaWRKY93, whose overexpression improves salt tolerance in A. thaliana. This indicates that a single WRKY gene can influence the complex process of salt tolerance in plants [21]. Although the specific mechanism of this TF is still unclear, its role in salt tolerance has generally been verified in Arabidopsis or the original species through transgenic methods [7]. At the same time, many other TFs are also involved in regulating gene expression induced by salt. High-throughput sequencing technology is helpful in identifying such TFs involved in the salt stress response [8].
Halophytes are important plant resources. Researchers mainly study their phenotype, physiology, and molecular biology to explore their salt tolerance characteristics and mechanisms and to provide theoretical bases for improving plant salt tolerance [22,23,24,25,26]. After analyzing the whole genome of Limonium bicolor, Yuan et al. [27] found the important genes affecting the formation of salt glands and the salt tolerance adaptation mechanism. Zhao et al. [28] analyzed transcriptome-level data of Arachis paniculata after salt stress and found that ethylene response factor (ERF), MYB, WRKY, and bZIP family TFs played important roles in plant salt tolerance.
As an important agricultural and animal husbandry area in northern China, Inner Mongolia has been experiencing an annual increase in saline alkali land area due to its arid climate and the fact that irrigation is the main planting method. Therefore, there is an urgent need to study salt-tolerant mesophytes that are suitable for growing in the ecological environment of the region in order to provide genetic resources for cultivating salt-tolerant and highly resistant agricultural and pastoral crops. Based on the previous analysis of the morphology, physiology, and transcriptome of Glaux maritima, this study analyzed the salt-tolerance-related TFs of the salt-secreting halophyte Glaux maritima at the transcription level and identified six important genes distributed in the EIN, bHLH, and bZIP TF families. Clustering and expression analyses of these genes revealed that they play key roles in the response to high-salt conditions and could provide valuable genetic resources for cultivating salt-tolerant crops.

2. Materials and Methods

Wild G. maritima seedlings were cultivated with 1/2 Hoagland nutrient solution for 15 days in the Key Laboratory of the Ministry of Education, College of Grassland, Resources, and Environment, Inner Mongolia Agricultural University (wild G. maritima seedlings were collected from the Hailiutu base of Inner Mongolia Agricultural University Science and Technology Park, Bikeqi Town, Tumet Left Banner, Hohhot, Inner Mongolia, N: 40°38′, E: 111°28′, Height: 1060 m.). The optimal growing conditions included a steady 25 °C temperature, 60% relative humidity, and 1600 lx of light. Selected G. maritima seedlings were subjected to NaCl concentrations of 0 (N0), 600 (N3), and 800 (N4) mM for 24 h [29,30]. (Attached Figure S1) Moreover, leaf samples of G. maritima were collected in triplicate for molecular analysis, directly frozen in liquid nitrogen, and kept at −80 °C for RNA isolation. [31].

2.1. RNA Isolation

Total RNA was extracted using a TRIzol reagent kit (Invitrogen, Carlsbad, CA, USA), as per the manufacturer’s instructions [31]. To verify the integrity of the RNA, we processed it through RNase-free agarose gel electrophoresis and then analyzed the results using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). Following total RNA isolation, DynabeadsTM Oligo(dT)25 beads (Thermo Fisher Scientific, Waltham, MA, USA) were used to remove rRNA from eukaryotic mRNA, and a Ribo-ZeroTM Magnetic Kit was used to do the same for bacterial mRNA (Illumina, San Diego, CA, USA) (Attached Figure S2, Supplementary Table S1). The enriched mRNA was fragmented using a fragmentation buffer before being used for reverse transcription into cDNA. Second-strand cDNA was synthesized using DNA polymerase I, deoxynucleoside triphosphates, RNase H, and a buffer. Furthermore, cDNA extraction was carried out using a QiaQuick PCR extraction kit (Qiagen, Venlo, The Netherlands). The cDNA products were end-repaired, had an extra base added, and ligated using Illumina sequencing adapters. Finally, agarose gel electrophoresis was used to separate the ligation products by size, which were then amplified by PCR and sequenced on an Illumina Novaseq 6000 (Illumina).

2.2. qRT-PCR Verification

TransScript One-Step gDNA Removal and cDNA Synthesis SuperMix (CAS:AT311; Transgen Biotech Co., Ltd., Beijing, China) was used to reverse transcribe RNA. PerfectStartTM Green qPCR SuperMix (CAS:AQ601, Transgen) containing 50 ng cDNA template and 50× passive reference dye for a final volume of 20 μL was used for qRT-PCR on each sample in a LightCycler 480 (Roche, Basel, Switzerland). The amplification protocol was as follows: 94 °C for 30 s, 45 cycles at 94 °C for 5 s, 60 °C for 34 s, and 72 °C for 10 s, and then melting curve analysis was carried out at 95 °C 15 s, 60 °C for 1 min, and 95 °C for 15 s. Supplementary Table S4 lists the sequences of the specific primers used for qRT-PCR. β-Actin was used as an internal control. All qRT-PCR analyses were performed in three replicates.

2.3. Data Analysis

Excel 2019 was used for preliminary data processing. Several variables were analyzed using univariate analysis in SPSS (version 19.0; IBM, Armonk, NY, USA) to demonstrate their importance across the range of salt-stress treatments, and the differentially expressed gene (DEG) statistics were based on the results of differential analysis. The screening of significantly DEGs was based on |log2 fold change| > 1 and a false discovery rate < 0.05. The expression patterns of the genes were clustered using trend analysis, followed by analysis using the Short Time-series Expression Miner (STEM v1.3.13; http://www.cs.cmu.edu/~jernst/stem accessed on 22 July 2022), and enrichment analysis was carried out using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) resources. To obtain enriched GO keywords and pathways, we used a threshold of Q ≤ 0.05. Relative gene expression levels were calculated using the 2−ΔΔCt method (Livak, 2001) and analyzed via analysis of variance (ANOVA) with a significance threshold of p < 0.05). The DEGs of the protein–protein interaction (PPI) network were predicted by STRING [32] (https://string-db.org/ accessed on 22 October 2022), and the drawn network diagram and network node degree were analyzed by Cytoscape 3.5.1. TFs were identified using Plant-TFDB 4.0 (http://planttfdb.cbi.pku.edu.cn/ accessed on 1 October 2022).

3. Results

3.1. Overall Analysis of Transcription Factors of Glaux Maritima under Salt Stress

Through assembly analysis and comparison of transcriptome data from 9 cDNA libraries of Glaux maritima, 1466 TF genes were counted and distributed across 57 TF families, mainly the ERF, bHLH, MYB-related, NAC, and WRKY families. The TF families with the largest numbers of genes were ERF (133, 9.07%), bHLH (123, 8.39%), MYB-related (95, 6.48%), NAC (94, 6.41%), WRKY (84, 5.73%), C2H2 (81, 5.53%), bZIP (70, 4.77%), and MYB (66, 4.50%) (Figure 1).

3.2. KEGG and GO Enrichment Analysis of All Transcription Factors

In order to better understand the main biochemical and signal transduction pathways related to the obtained TFs, KEGG pathway analysis was conducted on the TF family (Figure 2). A total of 95 TF genes were significantly enriched in 14 pathways, including plant hormone signal transduction (66 unigenes), the MAPK (Mitogen-activated protein kinases) signaling pathway (34 unigene), plant–pathogen interaction (18 unigene), circadian rhythm (5 unigene), base exception repair (1 unigene), and DNA replication (1 unigene). These annotations provide important information for studying the specific functions, processes, and mechanisms involved in the salt tolerance of halophytes.
To understand the salt tolerance characteristics of the TFs, GO annotation and enrichment analysis were carried out on all TFs after salt treatment (Figure 3). A total of 54 functions were annotated, and 1191 genes were enriched in these processes. The highly enriched molecular functions were binding (GO:0005488) and nucleic acid binding transcription factor activity (GO:0001071), and the highly enriched cell components were cell (GO:0005623), cell part (GO:0044464), and organelle (GO:0043226). In addition, biological processes such as the cellular process (GO:0009987), the metabolic process (GO:008152), the regulation of biological process (GO:0050789), and the biological regulation (GO:0065007) were highly enriched at all stages.

3.3. Trends in Transcription Factor Gene Expression under Different Levels of Salt Stress

All TF genes were subjected to trend analysis, and 820 genes with significant differences in expression were identified. Differential gene expression analysis showed that these genes were distributed in eight expression modules (Figure 4). The genes with significantly different expression levels after different salt treatments (p < 0.05) were all distributed in the four modules of Profile 0, 1, 6, and 7, with 83, 239, 222, and 64 genes, respectively. These DEGs were mainly distributed in the ERF (96, 11.71%), bHLH (75, 9.15%), WRKY (67, 8.17%), NAC (65, 7.93%), MYB-related (49, 5.98%), C2H2 (47, 5.73%), MYB (45, 5.49%), bZIP (37, 4.51%), C3H (25, 3.05%), GRAS (25, 3.05%), and HD-ZIP (22, 2.68%) families (Figure 5). Comparing the numbers of DEGs in all TF families showed that most DEGs were concentrated in the top ten TF families mentioned above (Figure 6).

3.4. Differentially Expressed Transcription Factors under Salt Stress

The differentially expressed TF genes were identified. Under N3 and N4 salt concentration stress, 33 up-regulated and 35 down-regulated DEGs were identified (Figure 7a,c). Differentially expressed TFs included ERF, WRKY, EIN, MYC, MYB, ABF, ARF, and NAC family members. To further understand the expression trends of these single genes, a trend change diagram (Figure 7b,d) was constructed from the log2 FPKM values of these 68 genes. This comprehensive analysis showed that under N3 and N4 conditions, the expression of all single genes increased or decreased in a consistent trend.

3.5. Analysis of ABF, EIN, and MYC Transcription Factors Involved in the Salt Stress Response of Glaux Maritima

The plant hormone ABA activates a series of signaling pathways to activate the resistance mechanisms of plants under abiotic stress. ABF TFs play an important role in the regulation of the ABA signal pathway. They specifically recognize the AREB (ABA responsive element binding protein) class of basic leucine zipper proteins and are mainly involved in the regulation of ABA and stress. Positive regulatory ABF TFs related to salt stress have been found in Arabidopsis, tobacco, sweet potato, rape, and other crops [9,27,33,34]. In this study, five DEGs related to ABF TFs (Unigene0003302, Unigene0007479, Unigene0007480, Unigene0007482, and Unigene0078257) were found after salt stress in Glaux maritima. Unigene0078257 expression levels significantly increased with the increase in NaCl concentration, while the expression levels of the other four genes showed no significant changes. In order to better understand the relationship between these TFs and salt tolerance, a phylogenetic tree of ABF TFs was constructed using ABF protein sequences downloaded from the Plant TF Database (Figure 8a). Phylogenetic analysis showed that the five genes were listed in two branches. Unigene0007479, Unigene0007480, and Unigene0007482 were in the same branch as ABF TF proteins in model plants A. thaliana (AT4G34000.1 and AT4G34000.2) and Nicotiana tabacum (XP016458525.1 and XP016451400.1). Meanwhile, Unigene0003302 and Unigene0078257 had high homology with the ABF TF protein in Actinidia chinensis Planch (Achn026271 and Achn131571). These results provide a good theoretical basis for the cloning of and functional research on the ABF TFs of Glaux maritima.
bHLH, as the second largest eukaryotic TF family after MYB, participates in the defensive responses of plants and helps resist damage when plants are subjected to abiotic stress. After MYC2 was first found and its function identified in A. thaliana, a large number of studies showed that this TF participated in various stress responses and that it was also the main regulator of jasmonate signaling [35]. In this study, seven regulatory DEGs related to the TFs MYC1 and MYC2 (Unigene0006978, Unigene0009899, Unigene0027154, Unigene0027167, Unigene0046064, Unigene0069618, and Unigene0077023) were also found in G. maritima after salt stress. Other than Unigene0077023 and Unigene0006978, whose expression levels did not change under the different concentrations of salt stress, the genes were all negatively regulated. Among them, Unigene0009899 and Unigene0027167 expression levels decreased significantly as salt concentration increased, while the expression levels of other genes did not significantly change under salt stress. A phylogenetic tree analysis (Figure 8b) showed that the seven genes belonged to two different branches. Unigene0009899 and Unigene0027154, both related to MYC2 TFs, were in the same branch as MYC TFs in N. tabacum (XP016462596.1, XP016508165.1, and XP016449248.1) and A. chinensis (Achn31034.1). The other five genes were associated with Gossypium hirsutum (Gh_D12G2074), N. tabacum (XP016441449.1 and XP016499262.1), Capsella grandiflora (Cagra. 3356s0080.1. p), A. thaliana (AT5G46830.1, AT4G17880.1, and AT5G46760.1), and Arabidopsis lyrate (493087). Additionally, the MYC TF family genes in Brassica napus (GSBRNA2T00131278001 and GSBRNA2T00080269001) and Brassica oleracea (XP_013605597.1) are in the same branch. These results expand our understanding of the function of MYC TF family genes in salt tolerance and lay the foundation for further research.
In this study, nine regulatory genes related to EIN3/EIL TFs were found in the TF analysis of G. maritima. Three genes (Unigene0000457, Unigene0012139, and Unigene0036810) were related to EIL1/3 TFs, and six (Unigene0036956, Unigene0036957, Unigene0042139, Unigene0071273, Unigene0071274, and Unigene0071276) to EIN3 TFs. Phylogenetic tree analysis showed that Unigene0036957, Unigene0000457, and the other seven genes belong to three different branches. Unigene0000457 and genes found in Citrus clementina (Ciclev10000607m), A. thaliana (AT2G27050.1) and Nicotiana sylvestris (XP_009801714.1) belong to the same branch, while Unigene0036810, Unigene0036956, Unigene0042139, Unigene0071273, Unigene0071274, Unigene0071276, and genes found in A. chinensis (Achn347321 and Achn018431), A. thaliana (AT3G20770.1), Carica papaya (evm. model. superconting_103.60), and Citrus sinensis (orange1.1g007186m) were in adjacent branches. Unigene0000457, Unigene0012139, and Unigene0042139 expression levels were significantly up-regulated in response to salt stress, but no other genes of the EIL TF family showed significant expression changes.

3.6. Protein–Protein Interaction Network and Topological Analysis in Glaux Maritima

To better understand how DEGs in these TF families participate in protein regulation and mediation, we constructed a protein–protein interaction (PPI) network of the DEGs in the ABF, MYC, and EIN TF families. The results showed that the DEG related to the ABF family (Unigene0078257) interacted with 106 genes, forming a total of 472 interactions (Figure 9a). The DEGs related to the MYC family (Unigene0009899 and Unigene0027167) interacted with 349 and 348 genes, respectively (Figure 9b,c), to form 4531 and 3526 interactions, respectively, but there were no interactions with the DEGs related to the EIL family.

3.7. KEGG Pathway Analysis of Differentially Expressed Genes

KEGG analysis showed that salt stress significantly increased the expression of the ABF TF family gene Unigene0078257 downstream of SnRK2 in the ABA signaling pathway, and Unigene0078257 expression more than doubled under N3 and N4 conditions (Figure 10). In the MAPK signaling pathway, the three EIN3/EIL TF family DEGs related to ethylene signal transduction were all up-regulated. Compared with levels after control treatment, expression levels increased more than 1.1 times after N3 and N4 treatments. In the Jasmonate signaling pathway, the two DEGs related to the MYC2 family were down-regulated, with expression levels less than 0.4 times those observed after control treatment (Figure 11).

4. Discussion

Transcriptional regulation plays an important role in plant resistance to stress. Most regulation will be realized through interactions between specific TFs and corresponding cis-acting elements. TFs, as DNA-binding proteins that specifically interact with cis-acting elements in the gene promoter region, will activate or inhibit transcription on their own or with other proteins so that genes are expressed in specific places at specific times. Many different TF families related to salt stress, including WRKY, avian myeloblastosis viral oncogene homolog (MYB), NAM, ATAF1/2, and CUC2 (NAC), basic leucine zipper (bZIP), APETALA2/ethylene response factor (AP2/ERF), and basic helix loop helix (bHLH), among others, have been identified by analyzing the whole genomes of different plants [36]. With the increasing global saline-alkali land area, halophytes are receiving increasing attention from researchers. Previous studies on the salt tolerance characteristics and molecular mechanisms of halophytes have been conducted at phenotypic, physiological, and molecular levels [22,23,24,25,26]. Halophytes could provide better genetic resources for improving plants salt tolerance. At present, the salt-tolerant TF families found in halophytes mainly include AP2/EREBP, MYB, WRKY, NAC, HD ZIP, zinc finger domain (ZnF), G-box binding proteins, 14-3-3 protein analogs, MADS box proteins, and bHLH. These TFs show different expression patterns under salt stress and have specific or universal response mechanisms for certain times and locations. Additionally, they participate in complex regulatory networks in plants [37,38,39,40,41,42]. This is similar to the results of the study on the salt-tolerant TF family in the halophytic plant G. maritima. This study also found that TFs in the ERF, bHLH, and bZIP TF families play a crucial role in the adaptation of G. maritima to high-salt environments. In this study, by combining KEGG and GO enrichment analysis and trend analysis, we analyzed the DEGs in these TF families and found that six salt-tolerant key genes in the EIN, bZIP, and bHLH TF families were expressed in the Plant hormone signal pathway. These were similar to the research results on the gene expression regulation of ERF, MYB, bZIP, and bHLH TF families in Arabidopsis, tobacco, soybeans, Medicago truncatula, and Zea mays to enhance plant salt tolerance [43,44,45,46]. The molecular mechanism of salt tolerance in the halophyte G. maritima may be directly related to the regulation of these TFs.
Phytohormones are important endogenous molecules that regulate plant development and salt tolerance and play an important role in plants subjected to salt, drought, temperature, and other stressors. Therefore, the genes involved in the regulation of plant salt tolerance and in plant hormone signal transduction have attracted the attention of researchers. In attempts to understand and improve plant salt tolerance, researchers have revealed and verified a large number of positive regulatory genes in the ABA and ethylene-related pathways [47,48]. Using an immunoblotting experiment, Zhu et al. [49] found that AtCDPK4 and AtCDPK11 in Arabidopsis can phosphorylate AtABF1 and AtABF4 and improve tolerance to high-salt and drought stress by regulating stomatal movement. Liu et al. [50] found that under salt stress, the chlorophyll content of A. thaliana abf3 and abf4 mutants decreased, and the yellowing rate increased after salt stress, indicating that the AtABF3 and AtABF4 genes played key roles in the salt-stress response. When subjected to salt stress, transgenic potatoes overexpressing AtABF4 exhibited higher relative water content and lower stomatal conductance and transpiration rates than levels in wild-type potatoes, hallmarks of stronger salt resistance [51]. Transgenic A. thaliana overexpressing the Ipomoea batatas (L.) IbABF4 gene showed a salt-resistant phenotype with higher photosynthetic efficiency and lower malondialdehyde and hydrogen peroxide content than levels in wild-type plants [52]. BnaABF2 overexpression in Arabidopsis showed that transgenic plant leaves can improve salt tolerance by reducing stomatal aperture size and inhibiting water loss after salt stress [53]. In this study, we found that five ABF TF family genes involved in the ABA signal pathway may be key regulators of the G. maritima salt stress response. Phylogenetic tree analysis of the genes showed similarity to ABF TF protein sequences found in plants such as A. thaliana, N. tabacum, and A. chinensis. Additionally, PPI network analysis of Unigene0078257, a DEG that was up-regulated by salt stress, showed connections with 106 genes to form 472 interactions (Figure 9a). This result suggests that this gene plays a key role in regulating the salt tolerance of G. maritima and represents an excellent genetic resource for the breeding of salt-resistant crops.
As the first discovered gaseous plant hormone, ethylene is involved in the complete processes of plant growth and development and in various stress responses. The ethylene signaling pathway has also been clearly described. Ethylene insensitive 3/ethylene insensitive 3-like (EIN3/EIL) is a family of small TFs in plants mainly involved in ethylene signal transduction [54]. Some studies have shown that overexpression of EIN3/EIL TFs in Arabidopsis and grapes can enhance salt resistance [17,55,56]. As an important family of nuclear TFs in the ethylene signaling pathway, EIN3/EIL has attracted extensive research attention. The genes in this small TF family play a crucial role in the salt tolerance and stress resistance of plants [15,17,56,57,58,59]. In Arabidopsis, as the core component of ethylene signal transduction, EIN2 mutants are very sensitive to salt stress. CEND of EIN2 attenuates the salt hypersensitivity phenotype of the mutant. The EIN3 mutant responds normally to salt stress, while the EIN3/EIL1 double mutant has a phenotype similar to that of the EIN2 mutant and is extremely sensitive to high-salt conditions [60]. EIN3 and EIL1 regulate the ethylene response in seedlings, and seedlings have higher salt tolerance than do mature plants [15,16]. In this study, nine genes of the EIN3/EIL TF family were found to be enriched in the ethylene signaling pathway and may be related to salt tolerance in G. maritima. According to changes in gene expression after salt stress, three significantly up-regulated DEGs were discovered. These results extend our understanding of ethylene signaling pathway regulation in response to salt stress and help explain the function of EIN3/EIL genes in salt tolerance. However, PPI network analysis showed no interacting proteins, indicating that these three genes may be highly conserved. These results lay a theoretical foundation for the functional verification of regulatory elements of the ethylene signaling pathway.
In this study, in addition to finding the above DEGs coding for positive regulatory TFs related to the ABA and ethylene signaling pathways, we also found DEGs coding for negative regulators in the MYC TF family related to the jasmonate signaling pathway. MYC TFs are widely expressed in animals and plants and have multiple regulatory functions [61]. When the GhMYC4 TF in upland cotton was overexpressed in Arabidopsis, its salt tolerance and cold resistance were significantly improved, enhancing resistance of the plants to high-salt and drought conditions [62]. Identification and expression analysis of MYC TFs in tea plants showed that the TFs were related to amino acid anabolism under jasmonate stress and that CsMYC1, 3, 5, and 8 were significantly related to photosynthesis-related pathways. These results showed that under stress, jasmonate signaling may change the amino acid content and photosynthetic efficiency of tea plants through MYC family expression, thus improving adaptability under stress [63]. MYC2, a negative regulatory TF found in Arabidopsis, mediates salt sensitivity by inhibiting seed germination and delaying root growth; its knockout improves seed germination and plant root growth to better cope with salt stress [64]. These results are similar to the expression results of the two MYC family DEGs found in G. maritima. We found that these two DEGs interact with many proteins (Figure 9b,c), indicating that MYC TFs may play a key role in salt tolerance and stress resistance in G. maritima. The above results indicate that DEGs related to the MYC family in G. maritima participate most in protein regulation and mediation, which lays a theoretical foundation for explaining the function of salt resistance genes in G. maritima. These results have improved the theoretical basis for our subsequent functional verification experiments.

5. Conclusions

In this study, 1466 TFs genes were counted by analyzing the TF families in the TF data of the leaves of G. maritima under different salt concentrations; these genes were distributed across 57 TF families. After trend analysis, a total of 820 DEGs were counted, mainly distributed across 8 expression modules. After qRT-PCR expression verification, it was found that the TFs ABF2 (Unigene0078257) were related to ABA, EIN3 (Unigene0000457 and Unigene0012139) and EIL1 (Unigene0042139) were related to ethylene regulation, and MYC1/2 (Unigene0009899 and Unigene0027167) were related to the bHLH TF family. These TF genes may play an important role in regulating salt tolerance in G. maritima. These results provided a theoretical basis for improving the salt tolerance and stress resistance of plants and exploring the genetic resources for salt tolerance in halophytes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agriculture13071404/s1, Figure S1: Physiological and biochemical indexes of G. maritima leaves under salt stress.; Figure S2: Agilent 2100 RNA 6000 Nano kit detection Table S1: the quality of the isolated RNA samples; Table S2: Characteristics of G. maritima transcriptome sequencing; Table S3: Sequencing data of G. maritima leaves under different salt concentrations; Table S4: Primer sequences used for real-time fluorescence quantitative analysis.

Author Contributions

F.S. designed the experiment, R.G., Z.W. and M.Y. performed the experiments, R.G. and F.T. analyzed the data, R.G. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Science and Technology Project of the Inner Mongolia Autonomous Region (2020GG0063) and the Key Projects in Science and Technology of Inner Mongolia (2021ZD0031).

Data Availability Statement

Datasets are available at the NCBI project PRJNA845001. The transcriptome datasets are available in the NCBI Short Read Archive (SRA) database with the accession numbers SAMN28853762, SAMN28853763, SAMN28853764, SAMN28853765, SAMN28853766, SAMN28853767, SAMN28853768, SAMN28853769, and SAMN28853770.

Acknowledgments

Thanks to the Key Laboratory of Grassland Resources, Ministry of Education of P.R. of China for providing facilities.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lindemose, S.; O’Shea, C.; Jensen, M.K.; Skriver, K. Structure, function and networks of transcription factors involved in abiotic stress responses. Int. J. Mol. Sci. 2013, 14, 5842–5878. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Zandalinas, S.I.; Fichman, Y.; Devireddy, A.R.; Sengupta, S.; Azad, R.K.; Mittler, R. Systemic signaling during abiotic stress combination in plants. Proc. Natl. Acad. Sci. USA 2020, 117, 13810–13820. [Google Scholar] [CrossRef] [PubMed]
  3. Lamers, J.; van der Meer, T.; Testerink, C. How plants sense and respond to stressful environments. Plant Physiol. 2020, 182, 1624–1635. [Google Scholar] [CrossRef] [Green Version]
  4. Ahmad, A.; Niwa, Y.; Goto, S.; Ogawa, T.; Shimizu, M.; Suzuki, A.; Kobayashi, K.; Kobayashi, H. bHLH106 integrates functions of multiple genes through their G-Box to confer salt tolerance on Arabidopsis. PLoS ONE 2015, 10, e0126872. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Anwar, A.; Kim, J.K. Transgenic breeding approaches for improving abiotic stress tolerance: Recent progress and future perspectives. Int. J. Mol. Sci. 2020, 21, 2695. [Google Scholar] [CrossRef] [Green Version]
  6. Shamimuzzaman, M.; Vodkin, L. Genome-wide identification of binding sites for NAC and YABBY transcription factors and co-regulated genes during soybean seedling development by ChIP-seq and RNA-seq. BMC Genom. 2013, 14, 477. [Google Scholar] [CrossRef] [Green Version]
  7. Baillo, E.H.; Kimotho, R.N.; Zhang, Z.; Xu, P. Transcription factors associated with abiotic and biotic stress tolerance and their potential for crops improvement. Genes 2019, 10, 771. [Google Scholar] [CrossRef] [Green Version]
  8. Zhao, C.; Zhang, H.; Song, C.; Zhu, J.; Shabala, S. Mechanisms of Plant Responses and Adaptation to Soil Salinity. Innovation 2020, 1, 100017. [Google Scholar] [CrossRef]
  9. Choi, H.I.; Hong, J.H.; Ha, J.O.; Kang, J.Y.; Kim, S.Y. ABFs, a family of ABA-responsive element binding factors. J. Biol. Chem. 2000, 275, 1723–1730. [Google Scholar] [CrossRef] [Green Version]
  10. Uno, Y.; Furihata, T.; Abe, H.; Yoshida, R.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Arabidopsis basic leucine zipper transcription factors involved in an abscisic acid-dependent signal transduction pathway under drought and high-salinity conditions. Proc. Natl. Acad. Sci. USA 2000, 97, 11632–11637. [Google Scholar] [CrossRef]
  11. Kim, S.; Kang, J.Y.; Cho, D.I.; Park, J.H.; Soo, Y.K. ABF2, an ABRE-binding bZIP factor, is an essential component of glucose signaling and its overexpression affects multiple stress tolerance. Plant 2004, 40, 75–87. [Google Scholar] [CrossRef] [PubMed]
  12. Song, L.; Huang, S.C.; Wise, A.; Castanoz, R.; Nery, J.R.; Chen, H.; Watanabe, M.; Thomas, J.; Bar-Joseph, Z.; Ecker, J.R. A transcription factor hierarchy defines an environmental stress response network. Science 2016, 354, aag1550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Wei, D.; Wang, M.; Xu, F.; Quan, T.; Peng, K.; Xiao, L.; Xia, G. Wheat oxophytodienoate reductase gene TaOPR1 confers salinity tolerance via enhancement of abscisic acid signaling and reactive oxygen species scavenging. Plant Physiol. 2013, 161, 1217–1228. [Google Scholar]
  14. Zhao, Y.; Dong, W.; Zhang, N.; Ai, X.; Wang, M.; Huang, Z.; Xiao, L.; Xia, G. A wheat allene oxide cyclase gene enhances salinity tolerance via jasmonate signaling. Plant Physiol. 2014, 164, 1068–1076. [Google Scholar] [CrossRef] [Green Version]
  15. Achard, P.; Cheng, H.; De Grauwe, L.; Decat, J.; Schoutteten, H.; Moritz, T.; Van, D.; Straeten, D.; Peng, J.; Harberd, N.P. Integration of plant responses to environmentally activated phytohormonal signals. Science 2006, 311, 91–94. [Google Scholar] [CrossRef]
  16. Peng, J.; Li, Z.; Wen, X.; Li, W.; Shi, H.; Yang, L.; Zhu, H.; Guo, H. Salt-induced stabilization of EIN3/EIL1 confers salinity tolerance by deterring ROS accumulation in Arabidopsis. PLoS Genet. 2014, 10, e1004664. [Google Scholar] [CrossRef] [Green Version]
  17. Zhang, L.; Li, Z.; Quan, R.; Li, G.; Wang, R.; Huang, R. An AP2 domain-containing gene, ESE1, targeted by the ethylene signaling component EIN3 is important for the salt response in Arabidopsis. Plant Physiol. 2011, 157, 854–865. [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] [Green Version]
  19. Kim, J.; Kim, H.Y. Functional analysis of a calcium-binding transcription factor involved in plant salt stress signaling. FEBS Lett. 2006, 580, 5251–5256. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Babitha, K.C.; Ramu, S.V.; Pruthvi, V.; Mahesh, P.; Nataraja, K.; Udayakumar, M. Co-expression of AtbHLH17 and AtWRKY28 confers resistance to abiotic stress in Arabidopsis. Transgenic Res. 2013, 22, 327–341. [Google Scholar] [CrossRef]
  21. Qin, Y.; Tian, Y.; Liu, X. A wheat salinity-induced WRKY transcription factor TaWRKY93 confers multiple abiotic stress tolerance in Arabidopsis thaliana. Biochem. Biophys. Res. Commun. 2015, 464, 428–433. [Google Scholar] [CrossRef] [PubMed]
  22. Wilson, H.; Mycock, D.; Weiersbye, I.M. The salt glands of Tamarix usneoides E. Mey. Ex-Bunge (South African salt cedar). Int. J. Phytoremediat. 2017, 19, 587–595. [Google Scholar] [CrossRef]
  23. Yuan, F.; LYu, M.; Leng, B.Y.; Zheng, G.Y.; Feng, Z.T.; Li, P.H.; Zhu, X.G.; Wang, B.S. Comparative transcriptome analysis of developmental stages of the limonium bicolor leaf generates insights into salt gland differentiation. Plant Cell Environ. 2015, 38, 1637–1657. [Google Scholar] [CrossRef] [PubMed]
  24. Shabala, S.; Bose, J.; Hedrich, R. Salt bladders: Do they matter? Trends Plant Sci. 2014, 19, 687–691. [Google Scholar] [CrossRef] [PubMed]
  25. Yuan, Y.; Chen, M.; Chen, Q. Effects of NaCl stress on physiological and biochemical characteristics of pitaya (Hylocereus undatus) seedlings. Chin. J. Trop. Crops 2012, 33, 1256–1260. [Google Scholar]
  26. Ghanem, E.; Mohamed, E.; Kasem, A.; El-Ghamery, A.A. Differential salt tolerance strategies in three halophytes from the same ecological habitat: Augmentation of antioxidant enzymes and compounds. Plants 2021, 10, 1100. [Google Scholar] [CrossRef]
  27. Yuan, F.; Wang, X.; Zhao, B.; Xu, X.; Shi, M.; Leng, B.; Dong, X.; Lu, C.; Feng, Z.; Guo, J. The genome of the recreto halophyte Limonium bicolor provides insights into salt gland development and salinity adaptation during terrestrial evolution. Mol. Plants 2022, 15, 21. [Google Scholar]
  28. Zhao, H.; Jia, F.; Zhang, F.; Wang, Y. The transcriptome information analysis of differentially expressed genes of Halostachys caspica under salt stress. Chin. J. Bioinform. 2014, 12, 90–98. [Google Scholar]
  29. Zhang, H.H.; Li, X.; Guan, Y.P.; Li, M.; Wang, Y.; An, M.J.; Zhang, Y.H.; Liu, G.J.; Xu, N.; Sun, G.Y. Physiological and proteomic responses of reactive oxygen species metabolism and antioxidant machinery in mulberry (Morus alba L.) seedling leaves to NaCl and NaHCO3 stress. Ecotoxicol. Environ. Saf. 2020, 193, 110259. [Google Scholar]
  30. Ben, R.; Ben, R.W.; Bouteraa, M.T.; Olfa, J.; Anis, B.H. Lobularia maritima thioredoxin-h2 gene mitigates salt and osmotic stress damage in tobacco by modeling plant antioxidant system. Plant Growth Regul. 2022, 97, 101–115. [Google Scholar]
  31. Pandey, A.; Khan, M.K.; Hamurcu, M.; Brestic, M.; Topal, A.; Gezgin, S. Insight into the Root Transcriptome of a Boron-Tolerant Triticum zhukovskyi Genotype Grownunder Boron Toxicity. Agronomy 2022, 12, 2421. [Google Scholar] [CrossRef]
  32. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T) (Delta Delta C) method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  33. Joo, H.; Lim, C.W.; Lee, S.C. Roles of pepper bZIP transcription factor CaATBZ1 and its interacting partner RING-type E3 ligase CaASRF1 in modulation of ABA signaling and drought tolerance. Plant J. 2019, 100, 399–410. [Google Scholar] [CrossRef]
  34. Drogelaser, W.; Snoek, B.L.; Snel, B.; Weiste, C. The Arabidopsis bZIP transcription factor family-An update. Curr. Opin. Plant Biol. 2018, 45, 36–49. [Google Scholar] [CrossRef] [PubMed]
  35. Chen, X.; Huang, H.; Qi, T.; Liu, B.; Song, S. New perspective of the bHLH-MYB complex in jasmonate-regulated plant fertility in Arabidopsis. Plant Signal. Behav. 2016, 11, e1135280. [Google Scholar] [CrossRef] [Green Version]
  36. Khan, S.A.; Li, M.Z.; Wang, S.M.; Yin, H.J. Revisiting the role of plant transcription factors in the battle against abiotic stress. Int. J. Mol. Sci. 2018, 19, 1634. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Gruber, V.; Blanchet, S.; Diet, A.; Zahaf, O.; Boualem, A.; Kakar, K.; Alunni, B.; Udvardi, M.; Frugier, F.; Crespi, M. Identification of transcription factors involved in root apex responses to salt stress in Medicago truncatula. Mol. Genet. Genom. 2009, 281, 55–66. [Google Scholar] [CrossRef] [Green Version]
  38. Miyama, M.; Tada, Y. Transcriptional and physiological study of the response of Burma mangrove (Bruguiera gymnorhiza) to salt and osmotic stress. Plant Mol. Biol. 2008, 68, 119–129. [Google Scholar] [CrossRef]
  39. Wang, Y.; Chu, Y.; Liu, G.; Wang, M.H.; Jing, J.; Hou, Y.; Qu, G.; Yang, C. Identification of expressed sequence tags in an alkali grass (Puccinellia tenuiflora) cDNA library. Plant Physiol. 2007, 164, 78–89. [Google Scholar] [CrossRef]
  40. Wang, Y.; Yang, C.; Liu, G.; Liu, G.F.; Zhang, G.D.; Ban, Q.Y. Microarray and suppression subtractive hybridization analyses of gene expression in Puccnellia tenuiflora after exposure to NaHCO3. Plant Sci. 2007, 173, 309–320. [Google Scholar] [CrossRef]
  41. Wang, Y.; Yang, C.; Liu, G.; Zhang, G.; Ban, Q. Development of a cDNA microarray to identify gene expression of Puccinellia tenuiflora under saline-alkali stress. Plant Physiol. Biochem. 2007, 45, 567–576. [Google Scholar] [CrossRef] [PubMed]
  42. Jiang, Y.; Deyholos, M.K. Comprehensive transcriptional profiling of NaCl-stressed Arabidopsis roots reveals novel classes of responsive genes. BMC Plant Biol. 2006, 6, 25. [Google Scholar] [CrossRef] [Green Version]
  43. Chandler, J.W. Class VIIIb APETALA2 ethylene response factors in plant development. Trends Plant Sci. 2018, 23, 151–162. [Google Scholar] [CrossRef]
  44. Wu, L.J.; Chen, X.L.; Ren, H.Y.; Zhang, Z.J.; Zhang, H.W.; Wang, J.Y.; Wang, X.C.; Huang, R.F. ERF protein JERF1 that transcriptionally modulates the expression of abscisic acid biosynthesis-related gene enhances the tolerance under salinity and cold in tobacco. Planta 2007, 226, 815–825. [Google Scholar] [CrossRef] [PubMed]
  45. Mizoi, J.; Ohori, T.; Moriwaki, T.; Kidokoro, S.; Todaka, D.; Maruyama, K.; Kusakabe, K.; Osakabe, Y.; Shinozaki, K.; Yamaguchi-Shinozaki, K. GmDREB2A a canonical dehydration-responsive element-binding protein 2-type transcription factor in soybean, is post translationally regulated and mediates dehydration-responsive element-dependent gene expression. Plant Physiol. 2013, 161, 346–361. [Google Scholar] [CrossRef] [Green Version]
  46. Qin, F.; Kakimoto, M.; Sakuma, Y.; Maruyama, K.; Osakabe, Y.; Tran, L.S.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Regulation and functional analysis of ZmDREB2A in response to drought and heat stresses in Zea mays L. Plant 2007, 50, 54–69. [Google Scholar] [CrossRef]
  47. Chen, N.; Zhang, G.; Zhang, L.; An, Y.; Du, J.; Wang, D.; Guo, C. The Role of ABF Transcription Factors in Response to Abiotic Stress in Plant. J. Plant Genet. Resour. 2021, 22, 930–938. [Google Scholar]
  48. Zheng, Y.; Zhu, Z. Relaying the ethylene signal: New roles for EIN2. Trends Plant Sci. 2016, 21, 2–4. [Google Scholar] [CrossRef]
  49. Zhu, S.Y.; Yu, X.C.; Wang, X.J.; Zhao, R.; Li, Y.; Fan, R.C.; Shang, Y.; Du, S.Y.; Wang, X.F.; Wu, F.Q.; et al. Two calcium-dependent protein kinases, CPK4 and CPK11, regulate abscisic acid signal transduction in Arabidopsis. Plant Cell 2007, 19, 3019–3036. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Liu, S.; Zhu, M.K.; Liu, X.; Li, L. Studies on the responses to ABA and salt stress in ABF3 and ABF4-related mutants of Arabidopsis. J. South China Norm. Univ. 2012, 44, 101–104. [Google Scholar]
  51. García, M.N.M.; Cortelezzi, J.I.; Fumagalli, M.; Capiati, D.A. Expression of the Arabidopsis ABF4 gene in potato increases tuber yield improves tuber quality and enhances salt and drought tolerance. Plant Mol. Biol. 2018, 98, 137–152. [Google Scholar] [CrossRef]
  52. Wang, W.B.; Qiu, X.P.; Yang, Y.X.; Kim, H.S.; Jia, X.Y.; Yu, H.; Kwak, S.S. Sweet potato bZIP transcription factor IbABF4 confers tolerance to multiple abiotic stresses. Front. Plant Sci. 2019, 10, 630–645. [Google Scholar] [CrossRef] [Green Version]
  53. Zhao, B.Y.; Hu, Y.F.; Li, J.J.; Yao, X.; Liu, K.D. BnaABF2, a bZIP transcription factor from rapeseed (Brassica napus L.), enhances drought and salt tolerance in transgenic Arabidopsis. Bot. Stud. 2016, 57, 12–24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Song, J.; Zhu, C.; Zhang, X.; Zhang, X.; Wen, X.; Liu, L.; Peng, J.; Guo, H.; Yi, C. Biochemical and structural insights into the mechanism of DNA recognition by Arabidopsis ETHYLENE INSENSITIVE3. PLoS ONE 2015, 10, e0137439. [Google Scholar] [CrossRef]
  55. Wang, Y.; Liu, C.; Li, K.; Sun, F.; Hu, H.; Li, X.; Zhao, Y.; Han, C.; Zhang, W.; Duan, Y.; et al. Arabidopsis EIN2 modulates stress response through abscisic acid response pathway. Plant Mol. Biol. 2007, 64, 633–644. [Google Scholar] [CrossRef]
  56. Nai, G.; Lu, S.; Ma, W.; Li, Y.; Chen, B.; Mao, J. Genome-wide identification and expression characteristic analysis of EIN3/EIL transcription factor family in grape. J. Fruit Sci. 2021, 38, 856–870. [Google Scholar]
  57. Cao, W.H.; Liu, J.; He, X.J.; Mu, R.L.; Zhou, H.L.; Chen, S.Y.; Zhang, J.S. Modulation of ethylene responses affects plant salt-stress responses. Plant Physiol. 2007, 143, 707–719. [Google Scholar] [CrossRef] [Green Version]
  58. Guo, H.; Ecker, J.R. The ethylene signaling pathway: New insights. Curr. Opin. Plant Biol. 2004, 7, 40–49. [Google Scholar] [CrossRef] [PubMed]
  59. Lin, Z.; Zhong, S.; Grierson, D. Recent advances in ethylene research. J. Exp. Bot. 2009, 60, 3311–3336. [Google Scholar] [CrossRef] [Green Version]
  60. Lei, G.; Shen, M.; Li, Z.G.; Li, Z.; Zhang, B.; Duan, K.; Wang, N.; Cao, Y.; Zhang, W.; Ma, B.; et al. EIN2 regulates salt stress response and interacts with a MA3 domain containing protein ECIP1 in Arabidopsis. Plant Cell Environ. 2011, 34, 1678–1692. [Google Scholar] [CrossRef]
  61. Pires, N.; Dolan, L. Origin and diversification of basic helic-loop-helix proteins in plants. Mol. Biol. Evol. 2010, 27, 862–874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Gao, L.; Liu, B.; Li, J.; Wu, Y.; Tang, Y. Cloning and Function Analysis of bHLH Transcription Factor Gene GhMYC4 from Gossypium hirsutism L. J. Agric. Sci. Technol. 2016, 18, 33–41. [Google Scholar]
  63. Zheng, Y.C.; Gu, M.Y.; Bi, W.J.; Hu, Q.; Wang, P.; Ye, N.; Sun, Y. Genome-wide Analysis and Expression Pattern of MYC Family in Camellia sinensis. Fujian J. Agric. Sci. 2021, 36, 1007–1016. [Google Scholar]
  64. Verma, D.; Jalmi, S.K.; Bhagat, P.K.; Verma, N.; Sinha, A.K. A bHLH transcription factor, MYC2, imparts salt intolerance by regulating proline biosynthesis in Arabidopsis. FEBS J. 2020, 287, 2560–2576. [Google Scholar] [CrossRef] [PubMed]
Agriculture 13 01404 g001 550
Figure 1. Pie chart showing the percent share of identified transcription factors in each Transcription factor family. Note: The figure shows the distribution of 1466 TF genes in the top 8 TF families in the Glaux maritima. The 8 TF families were ERF (133, 9.07%), bHLH (123, 8.39%), MYB related (95, 6.48%), NAC (94, 6.41%), WRKY (84, 5.73%), C2H2 (81, 5.53%), bZIP (70, 4.77%), and MYB (66, 4.50%).
Figure 1. Pie chart showing the percent share of identified transcription factors in each Transcription factor family. Note: The figure shows the distribution of 1466 TF genes in the top 8 TF families in the Glaux maritima. The 8 TF families were ERF (133, 9.07%), bHLH (123, 8.39%), MYB related (95, 6.48%), NAC (94, 6.41%), WRKY (84, 5.73%), C2H2 (81, 5.53%), bZIP (70, 4.77%), and MYB (66, 4.50%).
Agriculture 13 01404 g001
Agriculture 13 01404 g002 550
Figure 2. KEGG pathway enrichment analysis of the transcription factor genes. Note: under salt stress, a total of 95 TF genes were successfully annotated in the KEGG database for Glaux maritima TF genes. The KEGG pathways with the highest number of annotated genes were Plant hormone signaling pathway (66), MAPK signaling pathway-plant (34) and plant pathogen interaction (18).
Figure 2. KEGG pathway enrichment analysis of the transcription factor genes. Note: under salt stress, a total of 95 TF genes were successfully annotated in the KEGG database for Glaux maritima TF genes. The KEGG pathways with the highest number of annotated genes were Plant hormone signaling pathway (66), MAPK signaling pathway-plant (34) and plant pathogen interaction (18).
Agriculture 13 01404 g002
Agriculture 13 01404 g003 550
Figure 3. GO enrichment analysis of the transcription factor genes. Note: under salt stress, a total of 54 functions were anno-tated, and 1191 genes were enriched in GO annotation and enrichment analysis. Among the three ontological processes, binding (GO:0005488) and nuclear acid binding transcription factor activity (GO:0001071) had the highest distribution of Molecular Func-tion; (GO:0005623), cell part (GO:0044464), and organelle (GO:0043226) had the highest distribution of Cell Component; cellular process (GO:0009987), metabolic process (GO:008152), regulation of biological process (GO:0050789), and biological regulation (GO:0065007) had the highest distribution of Biological Process.
Figure 3. GO enrichment analysis of the transcription factor genes. Note: under salt stress, a total of 54 functions were anno-tated, and 1191 genes were enriched in GO annotation and enrichment analysis. Among the three ontological processes, binding (GO:0005488) and nuclear acid binding transcription factor activity (GO:0001071) had the highest distribution of Molecular Func-tion; (GO:0005623), cell part (GO:0044464), and organelle (GO:0043226) had the highest distribution of Cell Component; cellular process (GO:0009987), metabolic process (GO:008152), regulation of biological process (GO:0050789), and biological regulation (GO:0065007) had the highest distribution of Biological Process.
Agriculture 13 01404 g003
Agriculture 13 01404 g004 550
Figure 4. Gene expression module diagram for trend analysis of all transcription factors. Note: After trend analysis, there were a total of 820 differentially expressed genes distributed in 8 expression modules. They were profile 0 (83 genes), profile 1 (239 genes), profile 2 (60 genes), profile 3 (62 genes), profile 4 (28 genes), profile 5 (62 genes), profile 6 (222 genes), and profile 7 (64 genes). In addition, the gene expression in the four modules of profile 0, 1, 6, and 7 showed significant differences after different concentrations of salt stress (p < 0.05). N0: 0 mM/L NaCl stress; N3: 600 mM/L NaCl stress; N4: 800 mM/L NaCl stress.
Figure 4. Gene expression module diagram for trend analysis of all transcription factors. Note: After trend analysis, there were a total of 820 differentially expressed genes distributed in 8 expression modules. They were profile 0 (83 genes), profile 1 (239 genes), profile 2 (60 genes), profile 3 (62 genes), profile 4 (28 genes), profile 5 (62 genes), profile 6 (222 genes), and profile 7 (64 genes). In addition, the gene expression in the four modules of profile 0, 1, 6, and 7 showed significant differences after different concentrations of salt stress (p < 0.05). N0: 0 mM/L NaCl stress; N3: 600 mM/L NaCl stress; N4: 800 mM/L NaCl stress.
Agriculture 13 01404 g004
Agriculture 13 01404 g005 550
Figure 5. Differentially expressed transcription factor gene statistics. Note: The figure shows the distribution of 820 TF genes in the top 11 TF families in the Glaux maritima. The 11 TF families were ERF (96, 11.71%), bHLH (75, 9.15%), WRKY (67, 8.17%), NAC (65, 7.93%), MYB-related (49, 5.98%), C2H2 (47, 5.73%), MYB (45, 5.49%), bZIP (37, 4.51%), C3H (25, 3.05%), GRAS (25, 3.05%), and HD-ZIP (22, 2.68%).
Figure 5. Differentially expressed transcription factor gene statistics. Note: The figure shows the distribution of 820 TF genes in the top 11 TF families in the Glaux maritima. The 11 TF families were ERF (96, 11.71%), bHLH (75, 9.15%), WRKY (67, 8.17%), NAC (65, 7.93%), MYB-related (49, 5.98%), C2H2 (47, 5.73%), MYB (45, 5.49%), bZIP (37, 4.51%), C3H (25, 3.05%), GRAS (25, 3.05%), and HD-ZIP (22, 2.68%).
Agriculture 13 01404 g005
Agriculture 13 01404 g006 550
Figure 6. Numbers of genes in the top 10 transcription factor families. Note: Statistics on the total number of TF in the top 10 TF families of Glaux maritima and the number of differentially expressed genes in trend analysis. The 10 TF families were ERF (TF: 133 genes, DEGs: 96 genes), bHLH (TF: 123 genes, DEGs: 75 genes), MYB-related (TF: 95 genes, DEGs: 49 genes), NAC (TF: 94 genes, DEGs: 65 genes), WRKY (TF: 84 genes, DEGs: 67 genes), C2H2 (TF: 81 genes, DEGs: 47 genes), bZIP (TF: 70 genes, DEGs: 37 genes), MYB (TF: 66 genes, DEGs: 45 genes), C3H (TF: 53 genes, DEGs: 25 genes) and GRAS (TF:48 genes, DEGs: 25 genes). TF: transcription factor, DEGs: differentially expressed genes.
Figure 6. Numbers of genes in the top 10 transcription factor families. Note: Statistics on the total number of TF in the top 10 TF families of Glaux maritima and the number of differentially expressed genes in trend analysis. The 10 TF families were ERF (TF: 133 genes, DEGs: 96 genes), bHLH (TF: 123 genes, DEGs: 75 genes), MYB-related (TF: 95 genes, DEGs: 49 genes), NAC (TF: 94 genes, DEGs: 65 genes), WRKY (TF: 84 genes, DEGs: 67 genes), C2H2 (TF: 81 genes, DEGs: 47 genes), bZIP (TF: 70 genes, DEGs: 37 genes), MYB (TF: 66 genes, DEGs: 45 genes), C3H (TF: 53 genes, DEGs: 25 genes) and GRAS (TF:48 genes, DEGs: 25 genes). TF: transcription factor, DEGs: differentially expressed genes.
Agriculture 13 01404 g006
Agriculture 13 01404 g007 550
Figure 7. Heat map of differential transcription factor gene expression and log2 RPKM expression trends. Note: (a): Heat map of 33 up-regulated differentially expressed genes under different salt concentrations; (b): log2 RPKM expression trends of 33 upregulated differentially expressed genes under different salt concentrations; (c): Heat map of 35 down-regulated differentially expressed genes under different salt concentrations; (d): log2 RPKM expression trends of 35 down-regulated differentially expressed genes under different salt concentrations. N0: 0 mM/L NaCl stress; N3: 600 mM/L NaCl stress; N4: 800 mM/L NaCl stress.3.5. Analysis of ABF, EIN, and MYC transcription factors involved in the salt stress response of Glaux maritima.
Figure 7. Heat map of differential transcription factor gene expression and log2 RPKM expression trends. Note: (a): Heat map of 33 up-regulated differentially expressed genes under different salt concentrations; (b): log2 RPKM expression trends of 33 upregulated differentially expressed genes under different salt concentrations; (c): Heat map of 35 down-regulated differentially expressed genes under different salt concentrations; (d): log2 RPKM expression trends of 35 down-regulated differentially expressed genes under different salt concentrations. N0: 0 mM/L NaCl stress; N3: 600 mM/L NaCl stress; N4: 800 mM/L NaCl stress.3.5. Analysis of ABF, EIN, and MYC transcription factors involved in the salt stress response of Glaux maritima.
Agriculture 13 01404 g007
Agriculture 13 01404 g008a 550Agriculture 13 01404 g008b 550
Figure 8. Phylogenetic tree analysis of genes in the ABF, MYC, and EIN3/EIL transcription factor families. Note: The genes were identified via trend analysis of TFs in Glaux maritima and are related to plant hormone signal transduction. (a): Phylogenetic tree analysis of ABF2 TFs. (b): Phylogenetic tree analysis of MYC1 and MYC2 TFs. (c): Phylogenetic tree analysis of EIN3 TFs.
Figure 8. Phylogenetic tree analysis of genes in the ABF, MYC, and EIN3/EIL transcription factor families. Note: The genes were identified via trend analysis of TFs in Glaux maritima and are related to plant hormone signal transduction. (a): Phylogenetic tree analysis of ABF2 TFs. (b): Phylogenetic tree analysis of MYC1 and MYC2 TFs. (c): Phylogenetic tree analysis of EIN3 TFs.
Agriculture 13 01404 g008aAgriculture 13 01404 g008b
Agriculture 13 01404 g009a 550Agriculture 13 01404 g009b 550Agriculture 13 01404 g009c 550
Figure 9. PPI network analysis on DEGs found in ABFs and MYC transcription factor families. Note: (a): The ABFs TF related differential gene (Unigene0078257) interacted with 106 genes, forming a total of 472 interactions; (b): The MYC TF related differential gene (Unigene0009899) interacted with 349 genes, forming a total of 4531 interactions; (c): The MYC TF related differential gene (Unigene0027167) interacted with 348 genes, forming a total of 3526 interactions.
Figure 9. PPI network analysis on DEGs found in ABFs and MYC transcription factor families. Note: (a): The ABFs TF related differential gene (Unigene0078257) interacted with 106 genes, forming a total of 472 interactions; (b): The MYC TF related differential gene (Unigene0009899) interacted with 349 genes, forming a total of 4531 interactions; (c): The MYC TF related differential gene (Unigene0027167) interacted with 348 genes, forming a total of 3526 interactions.
Agriculture 13 01404 g009aAgriculture 13 01404 g009bAgriculture 13 01404 g009c
Agriculture 13 01404 g010 550
Figure 10. Differentially expressed genes involved in plant hormone signal transduction. Note: N0: 0 mM/L NaCl stress; N3: 600 mM/L NaCl stress; N4: 800 mM/L NaCl stress.
Figure 10. Differentially expressed genes involved in plant hormone signal transduction. Note: N0: 0 mM/L NaCl stress; N3: 600 mM/L NaCl stress; N4: 800 mM/L NaCl stress.
Agriculture 13 01404 g010
Agriculture 13 01404 g011 550
Figure 11. Differentially expressed genes enriched in the MAPK signaling pathway-plant. Note: N0: 0 mM/L NaCl stress; N3: 600 mM/L NaCl stress; N4: 800 mM/L NaCl stress.
Figure 11. Differentially expressed genes enriched in the MAPK signaling pathway-plant. Note: N0: 0 mM/L NaCl stress; N3: 600 mM/L NaCl stress; N4: 800 mM/L NaCl stress.
Agriculture 13 01404 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gu, R.; Wan, Z.; Tang, F.; Shi, F.; Yan, M. Leaf Transcription Factor Family Analysis of Halophyte Glaux maritima under Salt Stress. Agriculture 2023, 13, 1404. https://doi.org/10.3390/agriculture13071404

AMA Style

Gu R, Wan Z, Tang F, Shi F, Yan M. Leaf Transcription Factor Family Analysis of Halophyte Glaux maritima under Salt Stress. Agriculture. 2023; 13(7):1404. https://doi.org/10.3390/agriculture13071404

Chicago/Turabian Style

Gu, Rui, Zhiqiang Wan, Fang Tang, Fengling Shi, and Mengjiao Yan. 2023. "Leaf Transcription Factor Family Analysis of Halophyte Glaux maritima under Salt Stress" Agriculture 13, no. 7: 1404. https://doi.org/10.3390/agriculture13071404

APA Style

Gu, R., Wan, Z., Tang, F., Shi, F., & Yan, M. (2023). Leaf Transcription Factor Family Analysis of Halophyte Glaux maritima under Salt Stress. Agriculture, 13(7), 1404. https://doi.org/10.3390/agriculture13071404

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