Next Article in Journal
COVID-19 Inflammatory Markers and Vitamin D Relationship in Pediatric Patients
Next Article in Special Issue
Role of Sodium Nitroprusside on Potential Mitigation of Salt Stress in Centaury (Centaurium erythraea Rafn) Shoots Grown In Vitro
Previous Article in Journal
Cellular Behaviors of Periodontal Ligament Stem Cells in the Presence of Bone Grafting Biomaterials, In-Vitro Study
Previous Article in Special Issue
Applications of Molecular Markers for Developing Abiotic-Stress-Resilient Oilseed Crops
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

DsDBF1, a Type A-5 DREB Gene, Identified and Characterized in the Moss Dicranum scoparium

by
Alfred O. Onele
,
Anastasia B. Mazina
,
Ilya Y. Leksin
and
Farida V. Minibayeva
*
Kazan Institute of Biochemistry and Biophysics, FRC Kazan Scientific Center, P.O. Box 261, 420111 Kazan, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Life 2023, 13(1), 90; https://doi.org/10.3390/life13010090
Submission received: 5 December 2022 / Revised: 23 December 2022 / Accepted: 25 December 2022 / Published: 28 December 2022
(This article belongs to the Special Issue Plant Biotic and Abiotic Stresses)

Abstract

:
Plant dehydration-responsive element binding (DREB) transcription factors (TFs) play important roles during stress tolerance by regulating the expression of numerous genes involved in stresses. DREB TFs have been extensively studied in a variety of angiosperms and bryophytes. To date, no information on the identification and characterization of DREB TFs in Dicranum scoparium has been reported. In this study, a new DBF1 gene from D. scoparium was identified by cloning and sequencing. Analysis of the conserved domain and physicochemical properties revealed that DsDBF1 protein has a classic AP2 domain encoding a 238 amino acid polypeptide with a molecular mass of 26 kDa and a pI of 5.98. Subcellular prediction suggested that DsDBF1 is a nuclear and cytoplasmic protein. Phylogenetic analysis showed that DsDBF1 belongs to group A-5 DREBs. Expression analysis by reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) revealed that DsDBF1 was significantly upregulated in response to abiotic stresses such as desiccation/rehydration, exposure to paraquat, CdCl2, high and freezing temperatures. Taken together, our data suggest that DsDBF1 could be a promising gene candidate to improve stress tolerance in crop plants, and the characterization of TFs of a stress tolerant moss such as D. scoparium provides a better understanding of plant adaptation mechanisms.

1. Introduction

Abiotic stresses, such as drought, salinity, heavy metals, high and low temperatures, can disrupt cellular homeostasis, resulting in redox imbalances and the accumulation of reactive oxygen species (ROS), which can limit plant development and crop productivity [1,2]. To react and adapt to these environmental challenges, plants have developed complex mechanisms including physiological, biochemical, and molecular processes [3,4,5]. Significant progress has been made over the past two decades regarding the identification and characterization of stress-responsive genes and proteins that directly protect plants against stresses [5,6]. Numerous genes are regulated by transcription factors (TF) in response to various abiotic stimuli [7]. Transcription factors play important roles in controlling the expression of genes in various signaling pathways due to their DNA-binding specificity [8,9].
Dehydration-responsive element binding (DREB) TFs play critical roles in plant tolerance by regulating the expression of stress-inducible genes during abiotic stresses [10]. They have been extensively identified in a wide variety of higher plants (angiosperms), such as Arabidopsis thaliana [11], soybean (Glycine max) [12], rice (Oryza sativa) [13,14], maize (Zea mays) [15], barley (Hordeum vulgare) [16], and others. The dehydration-responsive element (DRE)-binding factor, also known as DBF1, belongs to the APETALA2/Ethylene-Responsive Factor (AP2/ERF) TF family, which has been demonstrated to be involved in various biological processes in plants, including metabolism, development, and stress response [12,17,18]. So far, AP2/ERF genes annotated in the mosses Physcomitrium patens and Sphagnum fallax are the largest TF families found in the plant TF databases (TFDB), although the AP2/ERF gene family has been rarely studied in the moss species [19,20,21]. Additionally, it has been shown that AP2/ERFs are regulated in response to numerous stresses, such as salinity and UV in P. patens [22] and the gene PpDBF1 played a role in drought, salt, and cold tolerance in transgenic tobacco [23]. Furthermore, in the desiccation tolerant moss Syntrichia caninervis, AP2/ERFs were found to be the most abundant TFs [24].
Dicranum scoparium is a category “A” moss, one of the most desiccation tolerant moss species [25]. It is a widely distributed Holarctic moss that grows in various habitats and is one of the most polymorphic species in their genus [26]. Our preliminary analysis of class III peroxidase activity in three feather mosses such as D. scoparium, Hylocomium splendens and Pleurozium schreberi growing together in Aisha forest, Tatarstan, Russia revealed that they have high peroxidase activity and diverse peroxidase isoforms [27]. However, D. scoparium had the highest activity, approximately double that of H. splendens and P. schreberi, and this activity was stimulated by the desiccation/rehydration cycle. Therefore, D. scoparium was chosen for further investigation of desiccation tolerance mechanisms. Although to date the full genome of D. scoparium has not been sequenced and no reports of DREB families are available for this species, we have previously described in D. scoparium how temperature and desiccation/rehydration stresses change the expression of genes encoding Class I ascorbate peroxidase (DsAPX) and Class III peroxidases (DsPODs) [28]. We hypothesized that in Dicranum, abiotic stress will also influence the expression of DBF1. In this study, we report the isolation of a cDNA from D. scoparium that encodes a new DNA-binding TF, designated as DsDBF1. Furthermore, we analyzed the physico-chemical properties and subcellular localization of protein and gene expression patterns after desiccation/rehydration, high and low temperature, paraquat, DCMU, and CdCl2 stresses. Our study showed that DsDBF1 was significantly upregulated after exposure of D. scoparium to abiotic stresses, especially desiccation/rehydration, freezing temperature, paraquat, and CdCl2, suggesting that this TF plays multiple roles in the tolerance of the moss to abiotic stresses.

2. Materials and Methods

2.1. Plant Material

Dicranum scoparium Hedw. was collected in the Aisha Forest in Tatarstan, Russia (55°53 21.3 N 48°38 14.3 E). Plant material was placed between sheets of paper and left to dry slowly in the open air for 2 days before being stored in the refrigerator at +4 °C in the dark until usage [28].

2.2. Identification and Retrieval of DsDBF1protein

Metatranscriptome data for the moss D. scoparium deposited to the Sequence Read Archive in the NCBI under accession numbers: PRJEB21674, ID: 393814 and PRJNA499105, ID: 499105 were extracted from the database [29,30]. The files were downloaded using the SRA Toolkit [31] and then converted to fastq format. FastQC software [32] was used to evaluate library quality control (QC). Adapter removal and trimming was done using Trimmomatic software version 0.39 [33]. After trimming, the reads were reassessed using the FastQC software [32]. Library assembly was performed using Trinity software [34]. All contaminants and foreign fungal and bacterial sequences were removed from the original data [29,30].
BLASTX [35] of representative sequences from mosses and other predicted taxa was used to determine the taxonomic classification of the identified transcripts (E-value < 1 × 10−6). To reduce transcript redundancy, moss transcripts containing the top hits were isolated to a separate file and filtered using the EvidentialGene package (https://sourceforge.net/projects/evidentialgene/, accessed on 10 October 2022) with the default parameters. The DBF1 amino acid sequence from P. patens was used as a query in a TBLASTN search [36] for a similar protein among the filtered transcripts from D. scoparium.
A DBF1 transcript was found after using TBLASTN and open reading frames (ORFs) were detected using the Augustus gene prediction and AssemblyPostProcessor tools in Galaxy version 1.0.3.0 (https://usegalaxy.org/, accessed on 10 October 2022). To confirm the domain-identifying members of the DBF1 gene family, the predicted sequence was submitted against PFAM [37], NCBI Conserved Domains Database (CDD) [38], InterProScan [39], and HMMER [40].

2.3. Cloning and Sequencing of DBF1 Gene

Total RNA was extracted using the GeneJET Plant RNA Purification Mini Kit (Thermo Scientific, Vilnius, Lithuania). RNA concentration and purity were assessed using a NanoDrop® ND-1000 spectrophotometer (Thermo Scientific, Waltham, MA, USA) and integrity was confirmed using 1% (m/v) agarose gel electrophoresis. First strand and double strand cDNA were synthesized using Evrogen Mint 2 synthesis kit according to manufacturer’s protocols.
To verify the DBF1 from D. scoparium identified in silico, the DBF1 sequence with the highest homology to DBF1 from P. patens was cloned into the pAL2-T vector (Evrogen, Moscow, Russia) using primers: F TGGGTTCACACGATGCGGA; R ACGCTTTGAATCCACTGACGG and then sequenced.

2.4. Sequence Analysis

BLASTN software available online at (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 12 October 2022) was used to perform a homology search to compare our sequenced DsDBF1 with other genes in the database. Files in Fasta format were downloaded from the NCBI database after BLAST search and then subjected to multiple sequence alignments using Clustal Omega [41] and ClustalW [42] in MEGA X [43]. The Expasy ProtParam tool [44] was used to predict the physico-chemical properties of the DsDBF1 protein, including molecular weight, isoelectric point, instability index, and grand average of hydropathicity (GRAVY). Subcellular localization was predicted by MULocDeep [45].
The homologous sequences of DsDBF1 proteins obtained after BLASTX and other known DREB proteins from the NCBI database were aligned by ClustalW [42] in MEGA X [43]. A phylogenetic tree was constructed in MEGA X [43] using the neighbor joining method for 1500 bootstraps [46]. Evolutionary distances were calculated using the Poisson correction method [47] and all ambiguous positions were removed by pairwise deletion.
The MEME suite (http://meme-suite.org/index.html, accessed on 12 October 2022) was used to analyze DREB protein sequences to find conserved motifs with the following parameters: zero or one site per sequence, number of motifs (1–10), motif width (6–50) [48]. After MEME, the motif map was rebuilt using the TBtools software [49].

2.5. Stress Treatments

For stress treatment in this study, we followed the protocol developed in our early studies [28]. Before the experiment, 2 cm apical stem segments of dry mosses were pre-hydrated at +4 °C for 24 h on wet filter paper. For stress treatments, 0.2 g moss segments were incubated in 20 mL of 100 µM paraquat (1,1-dimethyl-4,4-bipyridylium dichloride), 100 µM DCMU (3-(3,4-dichlorophenyl)-1,1-dimethylurea), or 100 µM CdCl2 for 1 and 12 h. Hydrated apical stem segments were also thermally stressed by their exposure to −20 °C or +30 °C for 1 or 12 h in a dark temperature-controlled chamber (Thermostat LOIP, St. Petersburg, Russia). In all treatments, hydrated mosses kept at room temperature served as controls.
For desiccation stress, three biological replicates per treatment were used, each containing 0.17 g dry mass from 2 cm apical stem segments. Initially, air-dry mosses were fully hydrated by immersing them in a 20 mL volume of distilled water for 1 h while slowly shaking them on an orbital shaker. Then, the hydrated moss was gently blotted with filter paper and placed in the desiccator above silica gel. Here, moss samples were taken at time 0 (after 1 h of hydration), and after 2, 24, and 72 h of desiccation. After 72 h of desiccation, moss samples were rehydrated for 0.5 and 2 h. The change in relative water content (RWC) was monitored according to the protocol previously described in [28].

2.6. RNA Extraction, cDNA Synthesis and RT-qPCR

Samples exposed to stresses were immersed in liquid nitrogen, then, each sample was ground into a fine powder. For RT-qPCR, 0.1 g of material from each replicate was immediately frozen in liquid nitrogen and stored at −80 °C until use. Extraction of total RNA from D. scoparium thalli was performed using the RNeasy Plant Mini kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. RNA concentration and purity were measured with NanoDrop® ND-1000 spectrophotometer (Thermo Scientific, Waltham, MA, USA), and the integrity was further evaluated by gel electrophoresis in a 1% (w/v) agarose gel. First strand cDNA was synthesized using protocols from the Evrogen Mint 2 synthesis kit.
The vector NTI Suite 9 software was used to design RT-qPCR primers with the following parameters: amplicon length from 60 to 300 bp and a Tm range of 55 to 65 °C. RT-qPCR was performed on CFX Connect™ Real-Time System (Bio-Rad Laboratories, Singapore) with qPCRmix-HS SYBR (Evrogen). The templates were amplified three times at 95 °C for 3 min followed by 40 cycles of amplification (94 °C for 10 s and 55/60 °C for 40 s). Melting curve analysis after RT-qPCR and gel electrophoresis examination of the amplified products were used to assess the specificity of the primers. The gene-specific primers used for RT-qPCR are listed in Table S1. Ribosomal RNA (18S), glyceraldehyde-3-phosphate dehydrogenase (GAPDH2) and α-tubulin (α-TUB1 and α-TUB2) were used as internal controls for RT-qPCR normalization [28].

2.7. Statistical Analysis

Three biological and six analytical replicates were used to run all reactions. Gene expression differences were assessed using normalized expression (Cq) in the Bio-Rad CFX MaestroTM/Software version 2.3 and were found to be significant for p ≤ 0.05 (*), p ≤ 0.01 (**), p ≤ 0.001 (***) after ANOVA and Shapiro–Wilk Normality tests. The standard errors of the mean are shown as vertical bars (n = 6).

3. Results

3.1. Characterization and Phylogenetic Analysis of DsDBF1

A DBF1 gene was identified based on the metatranscriptome data for the moss D. scoparium downloaded from the Sequence Read Archive in the NCBI. To verify the DBF1 gene identified from D. scoparium, specific primers were designed and the PCR product (717 bp) was cloned into the pAL2-T vector (Evrogen, Moscow, Russia) and then sequenced. Blasting the DBF1 sequence from D. scoparium after cloning and sequencing revealed high homology with ERF/DREBs of other mosses and vascular plants in the NCBI database. Further analyses of the protein sequence using PFAM [37], NCBI CDD [38], InterProScan [39] and HMMER [40] databases revealed that this protein had a classic AP2 domain structure (Figure 1) and was named DsDBF1. In addition, coding domain sequence (CDS) length (bp), subcellular localization, and physico-chemical properties such as protein length (aa), molecular weight (MW, kDa), isoelectric point (pI), instability index, and GRAVY were predicted (Table 1). The results showed that the cloned DsDBF1 encoded a 238 amino acid polypeptide (Figure S1) with a predicted molecular weight of 26 kDa and isoelectric point of 5.98. Calculation of the instability index classified the protein as unstable with a value greater than 40. A negative value of GRAVY indicated that DsDBF1 was hydrophilic and subcellular prediction showed that the protein was localized within the nucleus and cytoplasm (Table 1).
Sequence alignment analysis indicated that DsDBF1 shared high homology and a conserved AP2/ERF domain with other DREBs (Figure 1), but with low similarity in their overall amino acid sequences (Figure S2). Additionally, two conserved elements (YRG and RAYD) were found in the AP2/ERF domain after sequence analysis, although arginine (R) is replaced by lysine (K) in both the first YRG and second RAYD elements in DsDBF1 and some DREBs of other mosses (Figure 1). The homologous protein sequences obtained after BLASTP search of the DsDBF1 sequence and other known ERF/DREB proteins from GenBank were used to construct a phylogenetic tree demonstrating the evolutionary relationship between DsDBF1 and other similar sequences from mosses and vascular plants. The evolutionary tree showed that DsDBF1 belongs to the A-5 group of the DREB subfamily as it shared a common ancestry and homology with other known A-5 DREBs from mosses such as S. caninervis, P. patens, Bryum argenteum, Pohlia nuntans, and vascular plants such as Selaginella moellendorffii, G. max, O. sativa, Citrus sinensis, Theobroma cacao, and Gossypium hirsutum (Figure 2). Furthermore, it was found that group A-5 was divided into seven subgroups, with S. moellendorffii positioning between the protein subgroups of mosses and the vascular plants. As shown in Figure 2, all other known DREBs from vascular plants were clustered into different DREB subfamilies such as A-1, A-2, A-4, and A-6.
Additionally, the results of MEME analyses showed that DsDBF1 contained a total of six motifs, among them, motifs 1–3 represented the basic conserved motifs that made up the AP2 domain (Figure 3). Motif 4 was absent only in the DREB protein of Triticum aestivum (TaDREBP1_AAL01124.1), while motif 5 was absent in the DREB proteins of Gossypium hirsutum (GhDREBP1 AAO43165.1 and GhDBP RAP2-4-like_NP 001314591.1), G. max (GmDREBP3 ABB36646.1 and GmDREBP_NP 001345278.1), O. sativa (OsRAP_XP 468111.1 and DREBP1A_XP 015610912.1), Z. mays (ZmDBF1_AAM80486.1) including T. aestivum (TaDREBP1_AAL01124.1). However, an additional motif 9 was detected in DsDBF1, which was only conserved in S. caninervis (DREBP5_AMT92109.1) and two in P. patens (ERF RAP2-1-like_XP 024372564.1 and PpDBF1_ABA43687.2) (Figure 3).

3.2. Expression Patterns of DsDBF1 in Response to Abiotic Stress Treatments

The expression pattern of DsDBF1 was studied after application of abiotic stresses such as desiccation/rehydration, exposure to DCMU, CdCl2, paraquat, high and freezing temperatures to moss apical segments. Desiccation of the hydrated moss for 2, 24, and 72 h over silica gel resulted in almost up to 94% loss of RWC in the moss samples, accompanied by downregulation of DsDBF1, with the lowest expression observed after 24 h of dehydration (Figure 4A). Rehydration of the mosses after 72 h of desiccation showed a gradual increase in DsDBF1 expression after 0.5 h and further 2 h with the expression of DsDBF1 2-fold higher compared to that in the hydrated mosses before desiccation (Figure 4A). Treatment of moss segments with an inhibitor of photosynthesis DCMU downregulated DsDBF1 expression after 1 and 12 h (Figure 4B). Subjecting the mosses to heavy metal CdCl2 and prooxidant paraquat significantly increased the expression of DsDBF1 after 1 h (Figure 4B); however, further treatment for 12 h downregulated DsDBF1 expression. No significant changes in DsDBF1 expression were observed after exposing moss to +30 °C for 1 and 12 h. Exposure of mosses to a freezing temperature of −20 °C reduced the level of DsDBF1 expression after 1 h (Figure 4B); however, further exposure for 12 h at −20 °C upregulated gene expression almost 10-fold compared to a 1 h treatment (Figure 4B).

4. Discussion

Members of the AP2/ERF family of TFs are among the most important key regulators of genes responsible for stress tolerance and developmental transitions of plants. These TFs regulate transcriptional networks to activate or repress gene expression in response to biotic and abiotic factors through the modulation of several signaling pathways [7,8,17,50]. In the last few years, many DREBs have been identified and characterized in several angiosperms, including A. thaliana [11], rice (O. sativa) [13,51], soybean (G. max) [52], maize (Z. mays) [15,53], cotton (G. hirsutum) [54,55], barley (H. vulgare) [16,56,57], wheat (T. aestivum) [58], Populus euphratica [59], Caragana korshinskii [60], and others. Surprisingly, the AP2/ERF gene family has been rarely studied in stress-tolerant moss species [19,20,21]. Several recent studies have shown that AP2/ERF TFs play an important role in the developmental processes and stress responses in some moss species, such P. patens [22,23,61], S. caninervis [24,62,63,64,65], B. argenteum [21], and P. nutans [66]. Dicranum scoparium is a desiccation-tolerant moss [25] whose genome has not been fully sequenced, and no TF families of this species have been reported to date. In this present study, we first identified in silico a cDNA of the DBF1 gene in the moss D. scoparium. Then, the identified gene was verified by cloning and sequencing. In addition, we performed molecular characterization of the protein, including analysis of the conserved domain, physico-chemical properties, subcellular localization, phylogenetic relationship, and motif analyses of identified DsDBF1 and DREBs of other plants, and finally, we examined the expression patterns of this gene in response to abiotic stresses. Our results demonstrate that the expression of DsDBF1 is strongly induced by rehydration after desiccation, and treatments with CdCl2, paraquat, and freezing temperature, providing insights into the roles of DBF1 in response of D. scoparium to abiotic stresses.
Analyses of the physico-chemical properties and the subcellular localization showed that DsDBF1 encodes a 238-amino acid polypeptide with a molecular weight of 25 kDa and a pI of 5.98 and the protein is localized within the nucleus and cytoplasm (Table 1). While the majority of TFs are nuclear localized, some are not when initially synthesized [67]. Some of these TFs are kept inactive in the cytoplasm when synthesized or expressed as membrane proteins, but when stimulated, they are activated by proteolytic cleavage, releasing the active form, which enters the nucleus and activates target genes [67,68].
Furthermore, the BLASTP search of the NCBI database revealed that DsDBF1 shares high sequence similarities with some DREBs from angiosperms and mosses. In addition, some uncharacterized proteins from mosses such as Ceratodon purpureus and S. fallax also show very high similarities to DsDBF1. The amino acid composition of the AP2 domain of DsDBF1 revealed that it contains 65 amino acid residues (Figure 1), which approximately corresponds to the conserved 60 amino acids of the AP2/ERF domain found in all DREBs [11]. Amino acid alignments of DREB proteins from different plants show high sequence similarity in the middle of AP2/ERF domain of these proteins (Figure 1), which is a significant feature of plant DREBs [64,69]. However, in general, outside the domain box, low similarity is observed in their overall amino acid sequences (Figure S2).
Analysis of the AP2/ERF domain after multiple sequence alignments revealed the presence of two conserved YRG and RYAD elements (Figure 1), although only glycine (G) is conserved in the YRG element among all the DREBs, while alanine (A) and aspartic acid (D) are conserved in the RYAD elements (Figure 1). Furthermore, in the first YRG element, tyrosine (Y) and arginine (R) are replaced by phenylalanine (F) and lysine (K), respectively, whereas in the second RYAD element, R is substituted by K, leucine (L), and histidine (H), and the Y is substituted by F and H. The AP2/ERF domain is a type of DNA-binding module that contains two known conserved elements (YRG and RAYD), and these two elements can bind with the promoter sequence or some other interacting proteins [69,70]. Studies have shown that YRG is involved in DNA binding activity and is the basic hydrophilic N-terminal side of the AP2/ERF domain. The N-terminal region is approximately 19 to 22 amino acids in length [50,69,70]. In addition, the second element, RAYD, is located in the acidic C-terminal region of the AP2/EREBP domain with a length of 42 to 43 amino acids. It is suggested that the RAYD element plays a crucial role in mediating protein–protein interactions [50,70]. However, in this study, the substitution of amino acids observed at various positions within the conserved elements in the AP2/ERF domain after multiple sequence alignments (Figure 1) of DREB proteins, may imply their functional divergence within DREB subfamilies.
To understand the evolutionary relationship between DsDBF1 and other well-known DREBs from other plants, a neighbor-joining tree was constructed using the deduced amino acid residues of these DREB proteins (Figure 2). In this analysis, DsDBF1 was found to belong to the A-5 group of the DREB subfamily as it shares a common ancestor with other known A-5 DREBs from mosses such as S. caninervis, P. patens, B. argenteum, P. nuntans, a lycophyte, for example, S. moellendorffii, and the angiosperms, such as G. max, O. sativa, C. sinensis, T. cacao, and G. hirsutum (Figure 2). In the A-5 subgroup, S. moellendorffii branches from the moss subgroup, positioning itself between the mosses and the angiosperms. This supports the report of early divergence of vascular plants from the ancient non-vascular plants [71]. It has been proposed that PpDBF1, an A-5 type DREB from P. patens, is an ancestor of DREB proteins and plays a general role in various stresses in non-vascular moss, which has diverged into different subclasses with different functions in the higher plants [23]. Consequently, the grouping of DsDBF1 and some other A-5 DREB proteins from mosses, lycophytes, and angiosperms in one clade suggests that they were established in the early stages of land plant evolution. Additionally, it was found that all other known DREBs from vascular plants diverged into different DREB subfamilies such as A-1, A-2, A-4, and A-6 (Figure 2). The DREB gene subfamily may have evolved and assumed new roles as a result of the divergence of the AP2 genes. The functional diversity and divergence of DREB genes during the adaptive evolution of stress signaling pathways in plants is most likely the result of subsequent duplication and transposition events [23].
Moreover, an investigation of the conserved motifs in DsDBF1 and other selected DREB proteins was carried out using MEME software. From the results, DsDBF1 contains a total of six motifs. Motifs 1–3 represent the conserved motifs of the AP2 domain (Figure 3). Furthermore, an additional motif 9 was detected in DsDBF1. This motif is present in S. caninervis (DREBP5_AMT92109.1) and P. patens (ERF RAP2-1-like_XP 024372564.1 and PpDBF1_ABA43687.2) (Figure 3), suggesting their common origin. Genome-wide sequence analysis of AP2/ERF family TFs in numerous plants revealed conserved regions and motifs on both sides of the AP2/ERF domain with important roles in transcriptional activity, protein–protein interactions, and nuclear localization. These conserved motifs can serve as an evidence for further classification of subgroups [50,72].
Plant DREB TFs play critical roles in the response to dehydration, salinity, and cold stresses [73,74]. To further understand the role of DsDBF1 in response to stresses, we examined the expression profile of DBF1 gene by RT-qPCR in the D. scoparium subjected to desiccation/rehydration, exposure to DCMU, CdCl2, paraquat, heat and freezing temperature. Our results indicate that DsDBF1 gene is upregulated by most of these stresses, suggesting that this gene is involved in D. scoparium response to abiotic stresses (Figure 4A,B). Surprisingly, DsDBF1 gene is downregulated following exposure of the moss to DCMU (Figure 4B). It has been reported that genes assigned to different groups within the same gene family show diverse stress response patterns and stress tolerance [64,75]. To date, most reports on DREB and Cold binding factors (CBFs) have mainly focused on DREBA1 and DREBA2, the largest among the subgroups [69,76]. The A-1 type DREBs (DREB1) are induced by cold and improve plant stress tolerance to low temperatures [73,77], whereas A-2 type DREBs (DREB2) play a major role in response to dehydration and heat stress, and improve drought and salt tolerance in plants [78]. A-5 DREBs have rarely been studied, and the functional and stress response mechanisms are still unclear [63].
Interestingly, desiccation of the hydrated mosses for 2, 24, and 72 h decreased DsDBF1 expression (Figure 4A). However, rehydrating moss thalli after 72 h of desiccation progressively increased DsDBF1 expression to 2-fold higher compared to the hydrated mosses before desiccation (Figure 4A). PpDBF1, a homolog of DsDBF1, was weakly induced by dehydration stress but strongly induced by ABA [23]. Out of ten A-5 type DREBs from S. caninervis, ScDREB5 was downregulated under rapid desiccation stress over silica gel, while ScDREB3, ScDREB9, and ScDREB10 were poorly induced by desiccation [63]. However, these four DREBs were significantly induced by cold stress, while ScDREB3 and ScDREB5 were upregulated during heat stress [63]. The rapid desiccation used in our experiment, in which moss thalli were dried over silica gel and reached an RWC of 6% [28] after 24 h, never occurs in boreal forests, where drying rates of mosses are much slower [27,79]. Other A-5-type DREBs such as GhDBP1 and GmDREB3 have been shown to improve plant stress tolerance [23,55,77]. Moreover, photosynthesis inhibitor DCMU decreased DsDBF1 expression (Figure 4B). Meanwhile, CdCl2 and paraquat significantly altered the expression of DsDBF1, as short-term exposure increased DsDBF1 expression after 1 h (Figure 4B). Long-term treatment of moss samples with paraquat and CdCl2 resulted in downregulation of DsDBF1 expression after 12 h. Furthermore, exposure of moss to +30 °C had little effect on DsDBF1 expression, although a freezing temperature of −20 °C for 12 h upregulated gene expression almost 10-fold compared to a 1 h cold treatment (Figure 4B). It has been reported that PpDBF1, GmDREB2, StDREB2, ScDREB1, ScDREB2, ScDREB4, ScDREB6, ScDREB7, and ScDREB8 responded to drought, salt, and cold treatment among members of the A-5 subgroups [23,63,80,81]. Taken together, the upregulation of DsDBF1 during rehydration after desiccation, exposure to CdCl2, paraquat, and freezing-temperature stress suggests that DsDBF1, like other A-5 DREBs, plays important roles in D. scoparium stress tolerance.

5. Conclusions

An A-5 type gene, DsDBF1, encoding DRE-binding transcription factor TF was identified and cloned in the moss D. scoparium. DsDBF1 protein was predicted to be localized within the nucleus and cytoplasm. Furthermore, RT-qPCR analysis showed that DsDBF1 expression was significantly induced in response to abiotic stresses such as desiccation/rehydration, exposure to paraquat, CdCl2, high and freezing temperatures. D. scoparium is a desiccation tolerant moss species. Based on our results, we believe that DsDBF1 could be a promising gene candidate to improve stress tolerance in various crop plants, and characterization of transcription factors of a stress-tolerant moss such as D. scoparium provides a better understanding of plant response and adaptation mechanisms.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/life13010090/s1, Figure S1: DsDBF1 coding domain sequence and protein sequence; Figure S2: Sequence alignments of DsDBF1 and other known ERF/DREB proteins from mosses and vascular plants such as Syntrichia caninervis (DREBP5_AMT92109.1), Physcomitrium patens (ERF RAP2-1-like_XP 024372564.1; PpDBF1 ABA43687.2; ERF_TINY-like_XP 024390306.1), Gossypium hirsutum (GhDREBP1_AAO43165.1; GhDBP_RAP2-4-like_NP 001314591.1), Glycine max (GmDREBP3_ABB36646.1; GmDREBP_NP 001345276.1; GmDREBP_NP 001345278.1), Oryza sativa (OsRAP_XP 468111.1; OsDREBP2A_XP 025878770.1; DREBP1A_XP 015610912.1; OsDREB1F_NP 001359120.1; ERF038_XP 015614793.1), Triticum aestivum (TaDREBP1_AAL01124.1; CRT/DREBP_XP 044398325.1), Zea mays (ZmDBF1_AAM80486.1; ZmDBF2_AAM80485.1), Capsicum annuum (CaCBF1B_AAQ88400.1), Capsella bursa-pastoris (CbCBF_AAR26658.1), Brassica napus (BnCBF_AAL38243.1), Prunus avium (PaDREB1F_XP 021803652.1), Arabidopsis thaliana (TINY2_NP 196720.1), Bryum argenteum (ERF_QDB64575.1), Pohlia nutans (ERF_QCF46602.1), Selaginella moellendorffii (ERF039-like_XP 024530345.1), Citrus sinensis (ERF016_XP 006474696.1), Apostasia shenzhenica (ERF034_PKA61103.1) and Theobroma cacao (ERF016_XP 007012585.2). Multiple alignment was performed using Clustal Omega. Amino acid sequences are highlighted with different colors. Sequences marked by (*) show conserved amino acid residues, Table S1: Primers of RT-qPCR.

Author Contributions

Conceptualization, A.B.M. and F.V.M.; methodology, A.B.M. and A.O.O.; software, A.O.O., A.B.M. and I.Y.L.; formal analysis, A.O.O. and I.Y.L.; investigation, A.B.M.; writing—original draft preparation, A.O.O.; writing—review and editing, F.V.M.; supervision, F.V.M.; project administration, F.V.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the state assignment of the FRC KazSC RAS. F.V.M. thanks Russian Foundation for Basic Research [grant number 20-04-00721] for partial financial support (gene cloning and sequencing).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Sachdev, S.; Ansari, S.A.; Ansari, M.I.; Fujita, M.; Hasanuzzaman, M. Abiotic stress and reactive oxygen species: Generation, signaling, and defense mechanisms. Antioxidants 2021, 10, 277. [Google Scholar] [CrossRef] [PubMed]
  2. Bhuyan, M.B.; Hasanuzzaman, M.; Parvin, K.; Mohsin, S.M.; Al Mahmud, J.; Nahar, K.; Fujita, M. Nitric oxide and hydrogen sulfide: Two intimate collaborators regulating plant defense against abiotic stress. Plant Growth Regul. 2020, 90, 409–424. [Google Scholar] [CrossRef]
  3. Garg, A.K.; Kim, J.-K.; Owens, T.G.; Ranwala, A.P.; Choi, Y.D.; Kochian, L.V.; Wu, R.J. Trehalose accumulation in rice plants confers high tolerance levels to different abiotic stresses. Proc. Natl. Acad. Sci. USA 2002, 99, 15898–15903. [Google Scholar] [CrossRef] [Green Version]
  4. Munns, R.; Tester, M. Mechanisms of salinity tolerance. Annu. Rev. Plant Biol. 2008, 59, 651–681. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Yamaguchi-Shinozaki, K.; Shinozaki, K. Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses. Annu. Rev. Plant Biol. 2006, 57, 781–803. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Chinnusamy, V.; Zhu, J.; Zhu, J.-K. Salt stress signaling and mechanisms of plant salt tolerance. Genet. Eng. 2006, 27, 141–177. [Google Scholar] [CrossRef] [Green Version]
  7. Kavar, T.; Maras, M.; Kidrič, M.; Šuštar-Vozlič, J.; Meglič, V. Identification of genes involved in the response of leaves of Phaseolus Vulgaris to drought stress. Mol. Breed. 2008, 21, 159–172. [Google Scholar] [CrossRef]
  8. Shinozaki, K.; Yamaguchi-Shinozaki, K. Gene networks involved in drought stress response and tolerance. J. Exp. Bot. 2007, 58, 221–227. [Google Scholar] [CrossRef] [Green Version]
  9. Shinozaki, K.; Yamaguchi-Shinozaki, K.; Seki, M. Regulatory network of gene expression in the drought and cold stress responses. Curr. Opin. Plant Biol. 2003, 6, 410–417. [Google Scholar] [CrossRef]
  10. Seki, M.; Kamei, A.; Yamaguchi-Shinozaki, K.; Shinozaki, K. Molecular responses to drought, salinity and frost: Common and different paths for plant protection. Curr. Opin. Biotechnol. 2003, 14, 194–199. [Google Scholar] [CrossRef]
  11. Stockinger, E.J.; Gilmour, S.J.; Thomashow, M.F. Arabidopsis thaliana CBF1 encodes an AP2 domain-containing transcriptional activator that binds to the C-repeat/DRE, a cis-acting DNA regulatory element that stimulates transcription in response to low temperature and water deficit. Proc. Natl. Acad. Sci. USA 1997, 94, 1035–1040. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Mizoi, J.; Ohori, T.; Moriwaki, T.; Kidokoro, S.; Todaka, D.; Maruyama, K.; Kusakabe, K.; Osakabe, Y.; Shinozaki, K.; Yamaguchi-Shinozaki, K. GmDREB2A;2, a canonical dehydration-responsive element-binding protein2-type transcription factor in soybean, is posttranslationally regulated and mediates dehydration-responsive element-dependent gene expression. Plant Physiol. 2013, 161, 346–361. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Dubouzet, J.G.; Sakuma, Y.; Ito, Y.; Kasuga, M.; Dubouzet, E.G.; Miura, S.; Seki, M.; Shinozaki, K.; Yamaguchi-Shinozaki, K. OsDREB genes in rice, Oryza sativa L., encode transcription activators that function in drought-, high-salt- and cold-responsive gene expression. Plant J. 2003, 33, 751–763. [Google Scholar] [CrossRef] [PubMed]
  14. Nayak, S.N.; Balaji, J.; Upadhyaya, H.D.; Hash, C.T.; Kishor, P.B.K.; Chattopadhyay, D.; Rodriquez, L.M.; Blair, M.W.; Baum, M.; McNally, K.; et al. Isolation and sequence analysis of DREB2A homologues in three cereal and two legume species. Plant Sci. 2009, 177, 460–467. [Google Scholar] [CrossRef] [Green Version]
  15. Qin, F.; Kakimoto, M.; Sakuma, Y.; Maruyama, K.; Osakabe, Y.; Tran, L.-S.P.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Regulation and functional analysis of ZmDREB2A in response to drought and heat stresses in Zea mays L. Plant J. 2007, 50, 54–69. [Google Scholar] [CrossRef]
  16. Xue, G.-P. An AP2 domain transcription factor HvCBF1 activates expression of cold-responsive genes in barley through interaction with a (G/a)(C/t)CGAC motif. Biochim. Biophys. Acta 2002, 1577, 63–72. [Google Scholar] [CrossRef]
  17. Matsukura, S.; Mizoi, J.; Yoshida, T.; Todaka, D.; Ito, Y.; Maruyama, K.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Comprehensive analysis of rice DREB2-type genes that encode transcription factors involved in the expression of abiotic stress-responsive genes. Mol. Genet. Genom. 2010, 283, 185–196. [Google Scholar] [CrossRef]
  18. Licausi, F.; Ohme-Takagi, M.; Perata, P. APETALA2/Ethylene Responsive Factor (AP2/ERF) transcription factors: Mediators of stress responses and developmental programs. New Phytol. 2013, 199, 639–649. [Google Scholar] [CrossRef]
  19. Jin, J.; Tian, F.; Yang, D.-C.; Meng, Y.-Q.; Kong, L.; Luo, J.; Gao, G. PlantTFDB 4.0: Toward a central hub for transcription factors and regulatory interactions in plants. Nucleic Acids Res. 2017, 45, D1040–D1045. [Google Scholar] [CrossRef] [Green Version]
  20. Jin, J.; Zhang, H.; Kong, L.; Gao, G.; Luo, J. PlantTFDB 3.0: A portal for the functional and evolutionary study of plant transcription factors. Nucleic Acids Res. 2014, 42, D1182–D1187. [Google Scholar] [CrossRef]
  21. Li, X.; Gao, B.; Zhang, D.; Liang, Y.; Liu, X.; Zhao, J.; Zhang, J.; Wood, A.J. Identification, classification, and functional analysis of AP2/ERF family genes in the desert moss Bryum argenteum. Int. J. Mol. Sci. 2018, 19, 3637. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Hiss, M.; Laule, O.; Meskauskiene, R.M.; Arif, M.A.; Decker, E.L.; Erxleben, A.; Frank, W.; Hanke, S.T.; Lang, D.; Martin, A.; et al. Large-scale gene expression profiling data for the model moss Physcomitrella patens aid understanding of developmental progression, culture and stress conditions. Plant J. 2014, 79, 530–539. [Google Scholar] [CrossRef] [PubMed]
  23. Liu, N.; Zhong, N.-Q.; Wang, G.-L.; Li, L.-J.; Liu, X.-L.; He, Y.-K.; Xia, G.-X. Cloning and functional characterization of PpDBF1 gene encoding a DRE-Binding transcription factor from Physcomitrella patens. Planta 2007, 226, 827–838. [Google Scholar] [CrossRef] [PubMed]
  24. Gao, B.; Zhang, D.; Li, X.; Yang, H.; Wood, A.J. De novo assembly and characterization of the transcriptome in the desiccation-tolerant moss Syntrichia caninervis. BMC Res. Notes 2014, 7, 490. [Google Scholar] [CrossRef] [Green Version]
  25. Wood, A.J. The nature and distribution of vegetative desiccation-tolerance in hornworts, liverworts and mosses. Bryologist 2007, 110, 163–177. [Google Scholar] [CrossRef]
  26. Lang, A.; Naciri, Y. New chloroplast primers for intraspecific variation in Dicranum scoparium Hedw. (Dicranaceae) and amplification success in other bryophyte species. Mol. Ecol. Resour. 2010, 10, 735–737. [Google Scholar] [CrossRef]
  27. Onele, A.O.; Chasov, A.; Viktorova, L.; Beckett, R.P.; Trifonova, T.; Minibayeva, F. Biochemical characterization of peroxidases from the moss Dicranum scoparium. S. Afr. J. Bot. 2018, 119, 132–141. [Google Scholar] [CrossRef]
  28. Onele, A.O.; Chasov, A.V.; Viktorova, L.V.; Minibayeva, F.V.; Beckett, R.P. Characterization and expression analysis of ascorbate peroxidase from the moss Dicranum scoparium during abiotic stresses. Bryologist 2021, 124, 68–84. [Google Scholar] [CrossRef]
  29. Chen, K.-H.; Liao, H.-L.; Arnold, A.E.; Bonito, G.; Lutzoni, F. RNA-based analyses reveal fungal communities structured by a senescence gradient in the moss Dicranum scoparium and the presence of putative multi-trophic fungi. New Phytol. 2018, 218, 1597–1611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Chen, K.-H.; Liao, H.-L.; Bellenger, J.-P.; Lutzoni, F. Differential gene expression associated with fungal trophic shifts along the senescence gradient of the moss Dicranum scoparium. Environ. Microbiol. 2019, 21, 2273–2289. [Google Scholar] [CrossRef]
  31. Leinonen, R.; Sugawara, H.; Shumway, M.; International nucleotide sequence database collaboration. The sequence read archive. Nucleic Acids Res. 2011, 39, D19–D21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Andrews, S. FASTQC. A Quality Control Tool for High throughput Sequence Data. 2010. Available online: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 7 November 2022).
  33. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Haas, B.J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.D.; Bowden, J.; Couger, M.B.; Eccles, D.; Li, B.; Lieber, M.; et al. De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nat. Protoc. 2013, 8, 1494–1512. [Google Scholar] [CrossRef] [PubMed]
  35. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  36. Altschul, S.F.; Madden, T.L.; Schäffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [Green Version]
  37. El-Gebali, S.; Mistry, J.; Bateman, A.; Eddy, S.R.; Luciani, A.; Potter, S.C.; Qureshi, M.; Richardson, L.J.; Salazar, G.A.; Smart, A.; et al. The Pfam protein families database in 2019. Nucleic Acids Res. 2019, 47, D427–D432. [Google Scholar] [CrossRef]
  38. Marchler-Bauer, A.; Bo, Y.; Han, L.; He, J.; Lanczycki, C.J.; Lu, S.; Chitsaz, F.; Derbyshire, M.K.; Geer, R.C.; Gonzales, N.R.; et al. CDD/SPARCLE: Functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 2017, 45, D200–D203. [Google Scholar] [CrossRef] [Green Version]
  39. Mitchell, A.L.; Attwood, T.K.; Babbitt, P.C.; Blum, M.; Bork, P.; Bridge, A.; Brown, S.D.; Chang, H.-Y.; El-Gebali, S.; Fraser, M.I.; et al. InterPro in 2019: Improving coverage, classification and access to protein sequence annotations. Nucleic Acids Res. 2019, 47, D351–D360. [Google Scholar] [CrossRef] [Green Version]
  40. Potter, S.C.; Luciani, A.; Eddy, S.R.; Park, Y.; Lopez, R.; Finn, R.D. HMMER web server: 2018 update. Nucleic Acids Res. 2018, 46, W200–W204. [Google Scholar] [CrossRef] [Green Version]
  41. Sievers, F.; Higgins, D.G. Clustal Omega for making accurate alignments of many protein sequences. Protein Sci. 2018, 27, 135–145. [Google Scholar] [CrossRef]
  42. Thompson, J.D.; Higgins, D.G.; Gibson, T.J. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22, 4673–4680. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. 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–1549. [Google Scholar] [CrossRef] [PubMed]
  44. Gasteiger, E.; Hoogland, C.; Gattiker, A.; Duvaud, S.; Wilkins, M.R.; Appel, R.D.; Bairoch, A. Protein Identification and Analysis Tools on the ExPASy Server. In The Proteomics Protocols Handbook; Walker, J.M., Ed.; Springer Protocols Handbooks; Humana Press: Totowa, NJ, USA, 2005; pp. 571–607. [Google Scholar] [CrossRef]
  45. Jiang, Y.; Wang, D.; Yao, Y.; Eubel, H.; Künzler, P.; Møller, I.M.; Xu, D. MULocDeep: A deep-learning framework for protein subcellular and suborganellar localization prediction with residue-level interpretation. Comput. Struct. Biotechnol. J. 2021, 19, 4825–4839. [Google Scholar] [CrossRef] [PubMed]
  46. Saitou, N.; Nei, M. The Neighbor-joining Method: A New Method for Reconstructing Phylogenetic Trees. Mol. Biol. Evol. 1987, 4, 406–425. [Google Scholar] [CrossRef]
  47. Zuckerkandl, E.; Pauling, L. Evolutionary Divergence and Convergence in Proteins. In Evolving Genes and Proteins; Bryson, V., Vogel, H.J., Eds.; Academic Press: Cambridge, MA, USA, 1965; pp. 97–166. [Google Scholar] [CrossRef]
  48. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME SUITE: Tools for motif discovery and searching. Nucleic Acids Res. 2009, 37, W202–W208. [Google Scholar] [CrossRef]
  49. 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]
  50. Rehman, S.; Mahmood, T. Functional role of DREB and ERF transcription factors: Regulating stress-responsive network in plants. Acta Physiol. Plant. 2015, 37, 178. [Google Scholar] [CrossRef]
  51. Tian, X.-H.; Li, X.-P.; Zhou, H.-L.; Zhang, J.-S.; Gong, Z.-Z.; Chen, S.-Y. OsDREB4 genes in rice encode AP2-containing proteins that bind specifically to the dehydration-responsive element. J. Integr. Plant Biol. 2005, 47, 467–476. [Google Scholar] [CrossRef]
  52. Li, X.-P.; Tian, A.-G.; Luo, G.-Z.; Gong, Z.-Z.; Zhang, J.-S.; Chen, S.-Y. Soybean DRE-binding transcription factors that are responsive to abiotic stresses. Theor. Appl. Genet. 2005, 110, 1355–1362. [Google Scholar] [CrossRef]
  53. Kizis, D.; Pagès, M. Maize DRE-binding proteins DBF1 and DBF2 are involved in rab17 regulation through the drought-responsive element in an ABA-dependent pathway. Plant J. 2002, 30, 679–689. [Google Scholar] [CrossRef] [Green Version]
  54. Huang, B.; Liu, J.-Y. Cloning and functional analysis of the novel gene GhDBP3 encoding a DRE-Binding transcription factor from Gossypium hirsutum. Biochim. Biophys. Acta 2006, 1759, 263–269. [Google Scholar] [CrossRef] [PubMed]
  55. Huang, B.; Liu, J.-Y. A cotton dehydration responsive element binding protein functions as a transcriptional repressor of DRE-mediated gene expression. Biochem. Biophys. Res. Commun. 2006, 343, 1023–1031. [Google Scholar] [CrossRef] [PubMed]
  56. Xu, Z.-S.; Ni, Z.-Y.; Li, Z.-Y.; Li, L.-C.; Chen, M.; Gao, D.-Y.; Yu, X.-D.; Liu, P.; Ma, Y.-Z. Isolation and functional characterization of HvDREB1-a gene encoding a dehydration-responsive element binding protein in Hordeum Vulgare. J. Plant Res. 2009, 122, 121–130. [Google Scholar] [CrossRef]
  57. Choi, D.-W.; Rodriguez, E.M.; Close, T.J. Barley Cbf3 gene identification, expression pattern, and map location. Plant Physiol. 2002, 129, 1781–1787. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Xu, Z.-S.; Ni, Z.-Y.; Liu, L.; Nie, L.-N.; Li, L.-C.; Chen, M.; Ma, Y.-Z. Characterization of the TaAIDFa gene encoding a CRT/DRE-binding factor responsive to drought, high-salt, and cold stress in wheat. Mol. Genet. Genom. 2008, 280, 497–508. [Google Scholar] [CrossRef]
  59. Chen, J.; Xia, X.; Yin, W. Expression profiling and functional characterization of a DREB2-type gene from Populus euphratica. Biochem. Biophys. Res. Commun. 2009, 378, 483–487. [Google Scholar] [CrossRef] [PubMed]
  60. Wang, X.; Chen, X.; Liu, Y.; Gao, H.; Wang, Z.; Sun, G. CkDREB gene in Caragana korshinskii is involved in the regulation of stress response to multiple abiotic stresses as an AP2/EREBP transcription factor. Mol. Biol. Rep. 2011, 38, 2801–2811. [Google Scholar] [CrossRef]
  61. Do, T.H.; Pongthai, P.; Ariyarathne, M.; Teh, O.-K.; Fujita, T. AP2/ERF Transcription factors regulate salt-induced chloroplast division in the moss Physcomitrella patens. J. Plant Res. 2020, 133, 537–548. [Google Scholar] [CrossRef]
  62. Li, X.; Zhang, D.; Gao, B.; Liang, Y.; Yang, H.; Wang, Y.; Wood, A.J. Transcriptome-wide identification, classification, and characterization of AP2/ERF family genes in the desert moss Syntrichia caninervis. Front. Plant Sci. 2017, 8, 262. [Google Scholar] [CrossRef] [Green Version]
  63. Li, H.; Zhang, D.; Li, X.; Guan, K.; Yang, H. Novel DREB A-5 subgroup transcription factors from desert moss (Syntrichia caninervis) confers multiple abiotic stress tolerance to yeast. J. Plant Physiol. 2016, 194, 45–53. [Google Scholar] [CrossRef]
  64. Liang, Y.; Li, X.; Zhang, D.; Gao, B.; Yang, H.; Wang, Y.; Guan, K.; Wood, A.J. ScDREB8, a novel A-5 type of DREB gene in the desert moss Syntrichia caninervis, confers salt tolerance to Arabidopsis. Plant Physiol. Biochem. 2017, 120, 242–251. [Google Scholar] [CrossRef] [PubMed]
  65. Li, X.; Liang, Y.; Gao, B.; Mijiti, M.; Bozorov, T.A.; Yang, H.; Zhang, D.; Wood, A.J. ScDREB10, an A-5c type of DREB gene of the desert moss Syntrichia caninervis, confers osmotic and salt tolerances to Arabidopsis. Genes 2019, 10, 146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Li, C.; Liu, S.; Zhang, W.; Chen, K.; Zhang, P. Transcriptional profiling and physiological analysis reveal the critical roles of ROS-scavenging system in the Antarctic moss Pohlia nutans under Ultraviolet-B radiation. Plant Physiol. Biochem. 2019, 134, 113–122. [Google Scholar] [CrossRef] [PubMed]
  67. Pu, X.; Yang, L.; Liu, L.; Dong, X.; Chen, S.; Chen, Z.; Liu, G.; Jia, Y.; Yuan, W.; Liu, L. Genome-wide analysis of the MYB transcription factor superfamily in Physcomitrella Patens. Int. J. Mol. Sci. 2020, 21, 975. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Seo, P.J.; Kim, S.-G.; Park, C.-M. Membrane-bound transcription factors in plants. Trends Plant Sci. 2008, 13, 550–556. [Google Scholar] [CrossRef]
  69. Zhou, M.; Ma, J.-T.; Pang, J.-F.; Zhang, Z.-L.; Tang, Y.-X.; Sun, Z.-M. Regulation of plant stress response by dehydration responsive element binding (DREB) transcription factors. Afr. J. Biotechnol. 2010, 9, 9255–9269. [Google Scholar] [CrossRef]
  70. 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]
  71. Rensing, S.A.; Lang, D.; Zimmer, A.D.; Terry, A.; Salamov, A.; Shapiro, H.; Nishiyama, T.; Perroud, P.-F.; Lindquist, E.A.; Kamisugi, Y.; et al. The Physcomitrella genome reveals evolutionary insights into the conquest of land by plants. Science 2008, 319, 64–69. [Google Scholar] [CrossRef]
  72. Yan, H.W.; Hong, L.; Zhou, Y.Q.; Jiang, H.Y.; Zhu, S.W.; Fan, J.; Cheng, B.J. A genome-wide analysis of the ERF gene family in sorghum. Genet. Mol. Res. 2013, 12, 2038–2055. [Google Scholar] [CrossRef]
  73. Liu, Q.; Kasuga, M.; Sakuma, Y.; Abe, H.; Miura, S.; Yamaguchi-Shinozaki, K.; Shinozaki, K. Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low-temperature-responsive gene expression, respectively, in Arabidopsis. Plant Cell 1998, 10, 1391–1406. [Google Scholar] [CrossRef] [Green Version]
  74. Agarwal, P.K.; Agarwal, P.; Reddy, M.K.; Sopory, S.K. Role of DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Rep. 2006, 25, 1263–1274. [Google Scholar] [CrossRef] [PubMed]
  75. Sakuma, Y.; Liu, Q.; Dubouzet, J.G.; Abe, H.; Shinozaki, K.; Yamaguchi-Shinozaki, K. DNA-Binding Specificity of the ERF/AP2 Domain of Arabidopsis DREBs, Transcription Factors Involved in Dehydration- and Cold-Inducible Gene Expression. Biochem. Biophys. Res. Commun. 2002, 290, 998–1009. [Google Scholar] [CrossRef] [PubMed]
  76. Bouaziz, D.; Pirrello, J.; Charfeddine, M.; Hammami, A.; Jbir, R.; Dhieb, A.; Bouzayen, M.; Gargouri-Bouzid, R. Overexpression of StDREB1 transcription factor increases tolerance to salt in transgenic potato plants. Mol. Biotechnol. 2013, 54, 803–817. [Google Scholar] [CrossRef] [Green Version]
  77. Chen, M.; Xu, Z.; Xia, L.; Li, L.; Cheng, X.; Dong, J.; Wang, Q.; Ma, Y. Cold-induced modulation and functional analyses of the DRE-binding transcription factor gene, GmDREB3, in soybean (Glycine max L.). J. Exp. Bot. 2009, 60, 121–135. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Agarwal, P.; Agarwal, P.K.; Joshi, A.J.; Sopory, S.K.; Reddy, M.K. Overexpression of PgDREB2A transcription factor enhances abiotic stress tolerance and activates downstream stress-responsive genes. Mol. Biol. Rep. 2010, 37, 1125–1135. [Google Scholar] [CrossRef] [PubMed]
  79. Zotz, G.; Schweikert, A.; Jetz, W.; Westerman, H. Water relations and carbon gain are closely related to cushion size in the moss Grimmia pulvinata. New Phytol. 2000, 148, 59–67. [Google Scholar] [CrossRef] [Green Version]
  80. Chen, M.; Wang, Q.-Y.; Cheng, X.-G.; Xu, Z.-S.; Li, L.-C.; Ye, X.-G.; Xia, L.-Q.; Ma, Y.-Z. GmDREB2, a soybean DRE-binding transcription factor, conferred drought and high-salt tolerance in transgenic plants. Biochem. Biophys. Res. Commun. 2007, 353, 299–305. [Google Scholar] [CrossRef]
  81. Bouaziz, D.; Pirrello, J.; Ben Amor, H.; Hammami, A.; Charfeddine, M.; Dhieb, A.; Bouzayen, M.; Gargouri-Bouzid, R. Ectopic expression of dehydration responsive element binding proteins (StDREB2) confers higher tolerance to salt stress in potato. Plant Physiol. Biochem. 2012, 60, 98–108. [Google Scholar] [CrossRef]
Figure 1. Sequence alignments of DsDBF1 and other known ERF/DREB proteins showing the classical AP2 domains from mosses and vascular plants such as Syntrichia caninervis (DREBP5_AMT92109.1), Physcomitrium patens (ERF RAP2-1-like_XP 024372564.1; PpDBF1 ABA43687.2; ERF_TINY-like_XP 024390306.1), Gossypium hirsutum (GhDREBP1_AAO43165.1; GhDBP_RAP2-4-like_NP 001314591.1), Glycine max (GmDREBP3_ABB36646.1; GmDREBP_NP 001345276.1; GmDREBP_NP 001345278.1), Oryza sativa (OsRAP_XP 468111.1; OsDREBP2A_XP 025878770.1; DREBP1A_XP 015610912.1; OsDREB1F_NP 001359120.1; ERF038_XP 015614793.1), Triticum aestivum (TaDREBP1_AAL01124.1; CRT/DREBP_XP 044398325.1), Zea mays (ZmDBF1_AAM80486.1; ZmDBF2_AAM80485.1), Capsicum annuum (CaCBF1B_AAQ88400.1), Capsella bursa-pastoris (CbCBF_AAR26658.1), Brassica napus (BnCBF_AAL38243.1), Prunus avium (PaDREB1F_XP 021803652.1), Arabidopsis thaliana (TINY2_NP 196720.1), Bryum argenteum (ERF_QDB64575.1), Pohlia nutans (ERF_QCF46602.1), Selaginella moellendorffii (ERF039-like_XP 024530345.1), Citrus sinensis (ERF016_XP 006474696.1), Apostasia shenzhenica (ERF034_PKA61103.1), and Theobroma cacao (ERF016_XP 007012585.2). Amino acid sequences are highlighted with different colors. Sequences marked by (*) show conserved amino acid residues. Two conserved elements (YRG and RAYD) are marked by black horizontal lines. The differences in the moss conserved elements are shown in black frames. The red line shows the classical AP2 domain.
Figure 1. Sequence alignments of DsDBF1 and other known ERF/DREB proteins showing the classical AP2 domains from mosses and vascular plants such as Syntrichia caninervis (DREBP5_AMT92109.1), Physcomitrium patens (ERF RAP2-1-like_XP 024372564.1; PpDBF1 ABA43687.2; ERF_TINY-like_XP 024390306.1), Gossypium hirsutum (GhDREBP1_AAO43165.1; GhDBP_RAP2-4-like_NP 001314591.1), Glycine max (GmDREBP3_ABB36646.1; GmDREBP_NP 001345276.1; GmDREBP_NP 001345278.1), Oryza sativa (OsRAP_XP 468111.1; OsDREBP2A_XP 025878770.1; DREBP1A_XP 015610912.1; OsDREB1F_NP 001359120.1; ERF038_XP 015614793.1), Triticum aestivum (TaDREBP1_AAL01124.1; CRT/DREBP_XP 044398325.1), Zea mays (ZmDBF1_AAM80486.1; ZmDBF2_AAM80485.1), Capsicum annuum (CaCBF1B_AAQ88400.1), Capsella bursa-pastoris (CbCBF_AAR26658.1), Brassica napus (BnCBF_AAL38243.1), Prunus avium (PaDREB1F_XP 021803652.1), Arabidopsis thaliana (TINY2_NP 196720.1), Bryum argenteum (ERF_QDB64575.1), Pohlia nutans (ERF_QCF46602.1), Selaginella moellendorffii (ERF039-like_XP 024530345.1), Citrus sinensis (ERF016_XP 006474696.1), Apostasia shenzhenica (ERF034_PKA61103.1), and Theobroma cacao (ERF016_XP 007012585.2). Amino acid sequences are highlighted with different colors. Sequences marked by (*) show conserved amino acid residues. Two conserved elements (YRG and RAYD) are marked by black horizontal lines. The differences in the moss conserved elements are shown in black frames. The red line shows the classical AP2 domain.
Life 13 00090 g001
Figure 2. Phylogenetic analyses of DsDBF1 and other ERF/DREB proteins from mosses and vascular plant constructed using the neighbor-joining method with 1500 bootstrap test showing the relationship between the amino acid sequences. Evolutionary distances were calculated using the Poisson correction method and all ambiguous positions were removed by pairwise deletion. Amino acid sequences used for phylogenetic tree construction were retrieved, in part, from GenBank and after blast analysis from GenBank: Syntrichia caninervis (DREBP5_AMT92109.1), Physcomitrium patens (ERF RAP2-1-like_XP 024372564.1; PpDBF1 ABA43687.2; ERF_TINY-like_XP 024390306.1), Gossypium hirsutum (GhDREBP1_AAO43165.1; GhDBP_RAP2-4-like_NP 001314591.1), Glycine max (GmDREBP3_ABB36646.1; GmDREBP_NP 001345276.1; GmDREBP_NP 001345278.1), Oryza sativa (OsRAP_XP 468111.1; OsDREBP2A_XP 025878770.1; DREBP1A_XP 015610912.1; OsDREB1F_NP 001359120.1; ERF038_XP 015614793.1), Triticum aestivum (TaDREBP1_AAL01124.1; CRT/DREBP_XP 044398325.1), Zea mays (ZmDBF1_AAM80486.1; ZmDBF2_AAM80485.1), Capsicum annuum (CaCBF1B_AAQ88400.1), Capsella bursa-pastoris (CbCBF_AAR26658.1), Brassica napus (BnCBF_AAL38243.1), Prunus avium (PaDREB1F_XP 021803652.1), Arabidopsis thaliana (TINY2_NP 196720.1), Bryum argenteum (ERF_QDB64575.1), Pohlia nutans (ERF_QCF46602.1), Selaginella moellendorffii (ERF039-like_XP 024530345.1), Citrus sinensis (ERF016_XP 006474696.1), Apostasia shenzhenica (ERF034_PKA61103.1), and Theobroma cacao (ERF016_XP 007012585.2).
Figure 2. Phylogenetic analyses of DsDBF1 and other ERF/DREB proteins from mosses and vascular plant constructed using the neighbor-joining method with 1500 bootstrap test showing the relationship between the amino acid sequences. Evolutionary distances were calculated using the Poisson correction method and all ambiguous positions were removed by pairwise deletion. Amino acid sequences used for phylogenetic tree construction were retrieved, in part, from GenBank and after blast analysis from GenBank: Syntrichia caninervis (DREBP5_AMT92109.1), Physcomitrium patens (ERF RAP2-1-like_XP 024372564.1; PpDBF1 ABA43687.2; ERF_TINY-like_XP 024390306.1), Gossypium hirsutum (GhDREBP1_AAO43165.1; GhDBP_RAP2-4-like_NP 001314591.1), Glycine max (GmDREBP3_ABB36646.1; GmDREBP_NP 001345276.1; GmDREBP_NP 001345278.1), Oryza sativa (OsRAP_XP 468111.1; OsDREBP2A_XP 025878770.1; DREBP1A_XP 015610912.1; OsDREB1F_NP 001359120.1; ERF038_XP 015614793.1), Triticum aestivum (TaDREBP1_AAL01124.1; CRT/DREBP_XP 044398325.1), Zea mays (ZmDBF1_AAM80486.1; ZmDBF2_AAM80485.1), Capsicum annuum (CaCBF1B_AAQ88400.1), Capsella bursa-pastoris (CbCBF_AAR26658.1), Brassica napus (BnCBF_AAL38243.1), Prunus avium (PaDREB1F_XP 021803652.1), Arabidopsis thaliana (TINY2_NP 196720.1), Bryum argenteum (ERF_QDB64575.1), Pohlia nutans (ERF_QCF46602.1), Selaginella moellendorffii (ERF039-like_XP 024530345.1), Citrus sinensis (ERF016_XP 006474696.1), Apostasia shenzhenica (ERF034_PKA61103.1), and Theobroma cacao (ERF016_XP 007012585.2).
Life 13 00090 g002
Figure 3. Motif analysis of DsDBF1 and other known classic ERF/DREB proteins from mosses and vascular plants: Syntrichia caninervis (DREBP5_AMT92109.1), Physcomitrium patens (ERF RAP2-1-like_XP 024372564.1; PpDBF1 ABA43687.2; ERF_TINY-like_XP 024390306.1), Gossypium hirsutum (GhDREBP1_AAO43165.1; GhDBP_RAP2-4-like_NP 001314591.1), Glycine max (GmDREBP3_ABB36646.1; GmDREBP_NP 001345276.1; GmDREBP_NP 001345278.1), Oryza sativa (OsRAP_XP 468111.1; OsDREBP2A_XP 025878770.1; DREBP1A_XP 015610912.1; OsDREB1F_NP 001359120.1; ERF038_XP 015614793.1), Triticum aestivum (TaDREBP1_AAL01124.1; CRT/DREBP_XP 044398325.1), Zea mays (ZmDBF1_AAM80486.1; ZmDBF2_AAM80485.1), Capsicum annuum (CaCBF1B_AAQ88400.1), Capsella bursa-pastoris (CbCBF_AAR26658.1), Brassica napus (BnCBF_AAL38243.1), Prunus avium (PaDREB1F_XP 021803652.1), Arabidopsis thaliana (TINY2_NP 196720.1), Bryum argenteum (ERF_QDB64575.1), Pohlia nutans (ERF_QCF46602.1), Selaginella moellendorffii (ERF039-like_XP 024530345.1), Citrus sinensis (ERF016_XP 006474696.1), Apostasia shenzhenica (ERF034_PKA61103.1), and Theobroma cacao (ERF016_XP 007012585.2). Distribution of 10 putative conserved motifs in DREB proteins is shown. Conserved motifs are represented by different colored boxes numbered 1–10.
Figure 3. Motif analysis of DsDBF1 and other known classic ERF/DREB proteins from mosses and vascular plants: Syntrichia caninervis (DREBP5_AMT92109.1), Physcomitrium patens (ERF RAP2-1-like_XP 024372564.1; PpDBF1 ABA43687.2; ERF_TINY-like_XP 024390306.1), Gossypium hirsutum (GhDREBP1_AAO43165.1; GhDBP_RAP2-4-like_NP 001314591.1), Glycine max (GmDREBP3_ABB36646.1; GmDREBP_NP 001345276.1; GmDREBP_NP 001345278.1), Oryza sativa (OsRAP_XP 468111.1; OsDREBP2A_XP 025878770.1; DREBP1A_XP 015610912.1; OsDREB1F_NP 001359120.1; ERF038_XP 015614793.1), Triticum aestivum (TaDREBP1_AAL01124.1; CRT/DREBP_XP 044398325.1), Zea mays (ZmDBF1_AAM80486.1; ZmDBF2_AAM80485.1), Capsicum annuum (CaCBF1B_AAQ88400.1), Capsella bursa-pastoris (CbCBF_AAR26658.1), Brassica napus (BnCBF_AAL38243.1), Prunus avium (PaDREB1F_XP 021803652.1), Arabidopsis thaliana (TINY2_NP 196720.1), Bryum argenteum (ERF_QDB64575.1), Pohlia nutans (ERF_QCF46602.1), Selaginella moellendorffii (ERF039-like_XP 024530345.1), Citrus sinensis (ERF016_XP 006474696.1), Apostasia shenzhenica (ERF034_PKA61103.1), and Theobroma cacao (ERF016_XP 007012585.2). Distribution of 10 putative conserved motifs in DREB proteins is shown. Conserved motifs are represented by different colored boxes numbered 1–10.
Life 13 00090 g003
Figure 4. Expression patterns of DsDBF1 under abiotic stress treatments analyzed using RT-qPCR. (A) Relative expression of DsBF1 during desiccation over silica gel and rehydration. Shaded bars represent the hydrated and rehydrated moss, and solid bars represent the desiccated moss. (B) Relative expression of DsDBF1 exposed to DCMU, CdCl2, paraquat and high/low temperature for 1 and 12 h. Open bars correspond to control samples of mosses kept at room temperature. Red bars with white dots represent mosses treated with 100 µM DCMU, bars with orange horizontal stripes show moss treated with 100 µM CdCl2, bars with green stripes correspond to samples subjected to 100 µM paraquat, yellow and blue bars correspond to mosses exposed to +30 °C and −20 °C, respectively. p ≤ 0.05 (*), p ≤ 0.01 (**), p ≤ 0.001 (***). The vertical bars indicate the standard errors of the mean (n = 6).
Figure 4. Expression patterns of DsDBF1 under abiotic stress treatments analyzed using RT-qPCR. (A) Relative expression of DsBF1 during desiccation over silica gel and rehydration. Shaded bars represent the hydrated and rehydrated moss, and solid bars represent the desiccated moss. (B) Relative expression of DsDBF1 exposed to DCMU, CdCl2, paraquat and high/low temperature for 1 and 12 h. Open bars correspond to control samples of mosses kept at room temperature. Red bars with white dots represent mosses treated with 100 µM DCMU, bars with orange horizontal stripes show moss treated with 100 µM CdCl2, bars with green stripes correspond to samples subjected to 100 µM paraquat, yellow and blue bars correspond to mosses exposed to +30 °C and −20 °C, respectively. p ≤ 0.05 (*), p ≤ 0.01 (**), p ≤ 0.001 (***). The vertical bars indicate the standard errors of the mean (n = 6).
Life 13 00090 g004
Table 1. Physico-chemical properties and subcellular localization of DsDBF1.
Table 1. Physico-chemical properties and subcellular localization of DsDBF1.
ParametersDsDBF1
CDS length, bp717
Number of amino acids238
Molecular weight (kDa)26
Theoretical pI5.98
Instability index64.99
Grand average of hydropathicity (GRAVY)−0.224
Subcellular localization predictionNucleus/cytoplasm
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

Onele, A.O.; Mazina, A.B.; Leksin, I.Y.; Minibayeva, F.V. DsDBF1, a Type A-5 DREB Gene, Identified and Characterized in the Moss Dicranum scoparium. Life 2023, 13, 90. https://doi.org/10.3390/life13010090

AMA Style

Onele AO, Mazina AB, Leksin IY, Minibayeva FV. DsDBF1, a Type A-5 DREB Gene, Identified and Characterized in the Moss Dicranum scoparium. Life. 2023; 13(1):90. https://doi.org/10.3390/life13010090

Chicago/Turabian Style

Onele, Alfred O., Anastasia B. Mazina, Ilya Y. Leksin, and Farida V. Minibayeva. 2023. "DsDBF1, a Type A-5 DREB Gene, Identified and Characterized in the Moss Dicranum scoparium" Life 13, no. 1: 90. https://doi.org/10.3390/life13010090

APA Style

Onele, A. O., Mazina, A. B., Leksin, I. Y., & Minibayeva, F. V. (2023). DsDBF1, a Type A-5 DREB Gene, Identified and Characterized in the Moss Dicranum scoparium. Life, 13(1), 90. https://doi.org/10.3390/life13010090

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