Next Article in Journal
Detrimental Roles of Hypoxia-Inducible Factor-1α in Severe Hypoxic Brain Diseases
Next Article in Special Issue
Unraveling a Small Secreted Peptide SUBPEP3 That Positively Regulates Salt-Stress Tolerance in Pyrus betulifolia
Previous Article in Journal
Anti-Neuroinflammatory Effects of a Macrocyclic Peptide-Peptoid Hybrid in Lipopolysaccharide-Stimulated BV2 Microglial Cells
Previous Article in Special Issue
Genome-Wide Identification and Analysis of the EIN3/EIL Transcription Factor Gene Family in Doubled Haploid (DH) Poplar
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 25
Issue 8
10.3390/ijms25084464
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Systematic Analysis of DNA Demethylase Gene Families in Foxtail Millet (Setaria italica L.) and Their Expression Variations after Abiotic Stresses

by
Yingying Sun
1,
Xin Wang
1,
Yunfei Di
1,
Jinxiu Li
1,
Keyu Li
1,
Huanhuan Wei
1,
Fan Zhang
1 and
Zhenxia Su
1,2,*
1
College of Life Sciences, Shanxi University, Taiyuan 030006, China
2
Xinghuacun College (Shanxi Institute of Brewing Technology and Industry), Shanxi University, Taiyuan 030006, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(8), 4464; https://doi.org/10.3390/ijms25084464
Submission received: 12 March 2024 / Revised: 15 April 2024 / Accepted: 17 April 2024 / Published: 18 April 2024
(This article belongs to the Special Issue Latest Research on Plant Genomics and Genome Editing)
Browse Figures
Versions Notes

Abstract

:
DNA methylation is a highly conserved epigenetic modification involved in many biological processes, including growth and development, stress response, and secondary metabolism. DNA demethylase (DNA-deMTase) genes have been identified in some plant species; however, there are no reports on the identification and analysis of DNA-deMTase genes in Foxtail millet (Setaria italica L.). In this study, seven DNA-deMTases were identified in S. italica. These DNA-deMTase genes were divided into four subfamilies (DML5, DML4, DML3, and ROS1) by phylogenetic and gene structure analysis. Further analysis shows that the physical and chemical properties of these DNA-deMTases proteins are similar, contain the typical conserved domains of ENCO3c and are located in the nucleus. Furthermore, multiple cis-acting elements were observed in DNA-deMTases, including light responsiveness, phytohormone responsiveness, stress responsiveness, and elements related to plant growth and development. The DNA-deMTase genes are expressed in all tissues detected with certain tissue specificity. Then, we investigated the abundance of DNA-deMTase transcripts under abiotic stresses (cold, drought, salt, ABA, and MeJA). The results showed that different genes of DNA-deMTases were involved in the regulation of different abiotic stresses. In total, our findings will provide a basis for the roles of DNA-deMTase in response to abiotic stress.

1. Introduction

DNA methylation is a prominent epigenetic modification in many eukaryotes [1,2,3,4]. Research has shown that epigenetics is the study of heritable phenotypic changes that do not involve changes in the DNA sequence [5], and epigenetic mechanisms regulate critical agronomic traits in crops via DNA methylation, histone modifications, and small RNAs that affect gene expression and impact growth, seeding, germination, and fruit development [6]. DNA methylation is involved in transposons [7,8,9,10], the suppression of gene silencing, genomic imprinting [11], cell differentiation [12], etc., and is a necessary set up for the normal growth and development of organisms. DNA methylation has been shown to be a dynamic process that can be regulated according to different development periods or environmental conditions [13,14], thereby reducing the stress of plant survival in harsh environments and improving the ability of plants to cope with stress [15,16]. The reversibility of DNA methylation allows for the rapid and reversible modification of plant genomic DNA, avoiding excessive gene recombination and population diversity [17], thereby regulating gene expression and maintaining growth and development [18,19,20,21].
The precise status of DNA methylation depends on the function of DNA methyltransferase and DNA demethylase [22]. DNA methylation is catalysed by DNA methyltransferases and occurs by the addition of a methyl group to the C-5 site of cytosine, the N-6 site of adenine, and the N-7 site of guanine [23]. DNA demethylases contain a conserved DNA glycosilase motif and play a critical role in DNA demethylation by excision of the 5-methylcytosine in CG, CHG and CHH (where H stands for A, C or T) sequence contexts [24]. Increasing evidence shows that plant DNA methylation and demethylation are closely associated with various environmental stresses, including cold, drought, heat, heavy metal, salt, and ultraviolet stresses [25,26,27,28,29,30]. Under conditions of stress, the plant genome can overcome the limitation of genome instability through DNA methylation by rapid modification to maintain plant growth and development and evolutionary processes [31,32]. It was shown that abiotic stresses such as water deficit, cold, and salt stress induced demethylation of the repeats in the promoters of ACD6, ACO3, and GSTF14 and transcriptionally activated their expression [33]. It was also shown that the proportion of DNA methylated cytosines in Populus trichocarpa under drought stress was 10.04%, whereas it was only 7.75% in well-watered treatments [34]. Therefore, DNA methylation is considered the molecular response mechanism of plants in the face of adverse stress [8,9,31]. DNA demethylase is a bifunctional enzyme with glycosylase and depurine cleavage enzyme activities during active demethylation [35]. Currently, six DNA demethylases have been identified in plants, including transcriptional activator demeter (DME), Repressor of silencing 1 (ROS1), Demeter-like protein 2 (DML2), Demeter-like protein 3 (DML3), Demeter-like protein 4 (DML4) and Demeter-like protein 5 (DML5) [36].
Foxtail millet (Setaria italica [L.] P. Beauv), widely grown in China, Russia, India, and Africa, is one of the most important of ancient cultivated cereal crops. It was domesticated in northern China approximately 11,500 years ago [37,38]. As a C4 crop, S. italica is often used as an important model species for stress biology research [39,40,41]. DNA demethylase has already been identified and characterized in several plant species, including Arabidopsis thaliana, Oryza sativa, Daucus carota, Hordeum vulgare, Zea mays, Lycopersicon esculentum, tobacco, legumes and other plants [42,43,44,45,46,47,48,49]. However, no research has yet focused on the identification and analysis of DNA demethylase genes in S. italica, based on genome-wide analyses. Therefore, in this study, we performed a comprehensive analysis of S. italica DNA demethylase by genome-wide screening and identification, including phylogenetic analysis, gene structure analysis, conserved motifs construction, chromosomal location, and promoter sequence analysis. The expression levels of S. italica DNA demethylases under abiotic stresses such as salt, cold, and drought were also studied by qRT-PCR. These data will provide a valuable reference for further studies of DNA demethylase cascades.

2. Results

2.1. Genome-Wide Identification and Structural Analysis of DNA-deMTase Genes

Blast analysis of reported A. thaliana, O. sativa and T. aestivum DNA-deMTase proteins against the whole genome of S. italica. Subsequently, the sequences were aligned, redundant sequences were removed and combined, and Pfam and SMART structural domain searches were used to check whether the candidate genes contained ENCO3c domains, and a total of 7 candidate genes were identified in the S. italica genome (Table 1). The DNA-deMTase genes (SiDML5, SiDML3a, SiROS1b, SiROS1a, SiROS1d, SiDML3b, and SiDML4) code proteins composed of 286 to 2082 amino acids. The MW ranges from 31.52 to 230.61 kDa. The pI varies from 5.97 to 9.67, and the mean value was 7.43 which predicted that most of them were neutral proteins. The instability index of the DNA-deMTases proteins was greater than 40, and they were unstable proteins; GRAVY analysis showed that all the DNA-deMTase proteins were hydrophilic proteins. The prediction of subcellular localization showed that all DNA-deMTase proteins were localized in the nucleus, and SiDML5 was also localized in the chloroplast nucleus. Furthermore, the coding regions of DNA-deMTase genes were interrupted by 4 to 19 introns.

2.2. Phylogenetic Tree Construction of DNA-deMTases Gene Family Proteins and Analysis of Conserved Structural Domains

To reveal the evolutionary relationship of the DNA-dMTase genes in S. italica, a phylogenetic tree was constructed by comparing the protein sequences of the seven S. italica DNA-deMTase genes with those of the seven known A. thaliana DNA-deMTase genes, eight O. sativa DNA-deMTase genes, and eighteen T. aestivum DNA-deMTase genes, for a total of 40 DNA-deMTase genes using the maximum likelihood method. Based on the classification of O. sativa DNA-deMTases, the seven DNA-deMTases in cereals were classified into four subfamilies, including ROS1, DML3, DML4, and DML5. Among them, ROS1 was the most abundant, containing three DNA-deMTase genes, DML3 contained two DNA-deMTase genes, and there was one in each of the DML4 and DML5 subfamilies (Figure 1a). CDD and SMART were used to analyze the conserved domains of DNA-deMTases, and the results showed that all DNA-deMTases contained ENCO3c domain. The ENCO3c domains are DNA-deMTase-specific domains that perform the function of glycosidases. DML5 contains two domains, ENCO3c and FES; DML4 contains only one ENCO3c domain and ROS1 and DML3 subfamilies contain four common conserved domains: ENCO3c, FES, RRM_DME, and Perm-CXXC (Figure 1b).

2.3. Observation of Conserved Motif and Gene Structure Analysis of DNA-deMTases in S. italica

Genetic structure analysis is an important strategy to study genetic evolution. In this study, preliminary Pfam analysis results of the entire predicted proteins were applied to the MEME motif search tool. This analysis was used to further identify conserved motifs in the corresponding conserved domains of all seven DNA-deMTase proteins encoded by this gene family from S. italica. The search was performed for all predicted proteins and 12 conserved motifs were identified, but different subfamilies possess different motifs, consistent with their conserved domains. For example, the motifs of the ROS1 subfamily have the same composition and contain 12 conserved motifs; the DML3 subfamily contains nine motifs, including motifs 1, 2, 3, 4, 5, 6, 7, 9 and 10; the DML4 subfamily contains nine motifs, including motifs 1, 2, 3, 6 and 9; the DML5 subfamily contains nine motifs, including motifs 1, 2 and 6. (Figure 2a). Subsequently, an exon-intron structure analysis was performed to support the phylogeny reconstruction. The schematic structures revealed that each coding sequence of DNA-deMTase gene is disrupted by one or more exons. The DNA-deMTase genes within the same groups of the phylogenetic tree all showed similar exon-intron structures (Figure 2b).

2.4. Analysis of Cis-Acting Elements in DNA-deMTase Genes

Using the plant CARE database to analyze the 2000 bp cis-acting regulatory elements (CREs) upstream of seven DNA-deMTase genes, we found that the promoter of dMTase in S. italica contains types of CREs, including light response (17.5%), stress response (9.8%), hormone regulation (6%) and plant development (1.3%) (Figure 3). The DNA-deMTase promoter has various stress-responsive regulatory elements, such as the LTR involved in low-temperature stress, the MBS involved in drought induction, the TC-rich repeats involved in defense and stress responses, and the ARE regulatory element in response to anaerobic induction. There are also regulatory elements related to hormone signaling, such as the ABRE involved in the abscisic acid response, the P-box involved in the gibberellin response, and the TCA-element involved in the salicylic acid response. In addition, a number of tissue-specific cis-regulatory elements have been identified, such as the GCN4-motif involved in the GCN4-motif required for endosperm-specific expression, and AACA-motif involved in negative endosperm regulation, and AACA-motif associated with negative endosperm regulation. The motif required for endosperm-specific expression and AACA-motif associated with negative endosperm regulation, as well as the CAT-box associated with phloem-specific activation. In conclusion, the presence of multiple CREs in S. italica suggests that DNA-deMTase genes can be involved in growth and development, as in the stress response in S. italica.

2.5. Chromosomal Location and Synteny Analysis

The distribution of genes on chromosomes provides an important basis for studying the evolution and function of gene families. Combining the chromosomal information of the four S. italica genomes and the position of each DNA-deMTase gene on the chromosome, the distribution map of the S. italica demethylase gene on the chromosome was obtained (Figure 4). The seven identified genes were dispersed unevenly through nine chromosomes. Chromosome-1 contains two genes (SiROS1d and SiDML3a), whereas chromosomes 5, 7 and 8 did not contain any target genes, and the other chromosomes have one gene. Most of the DNA-deMTase genes are located in the middle of the chromosomes.

2.6. Expression Patterns of DNA-deMTase Genes in S. italica

To explore the expression pattern of the DNA-deMTase gene in different tissues in S. italica, publicly available Foxtail millet RNA-seq data (GeneAtlas v1 Tissue Sample from Phytozome 13) was used. This dataset included information from various tissues, including etiolated seedling, leaf, panicle, root, shoot, total aerial and germ shoot (Figure 5). As shown in Figure 5, SiROS1d was highly expressed in all tissues, while SiDML3b and SiDML4 were almost undetected in all tissues. SiROS1b, SiDML3a, SiDML5 and SiROS1a were highly expressed in most tissues. SiDML3a, SiDML5 and SiROS1a were all highly expressed in leaf 2–5, whereas the expression of the other genes was low, suggesting that these three genes may have similar functions in leaf development. SiDML3a was highly expressed in total aerial, indicating that probably SiDML3a has an important function in plant growth and development. Consistent with these results, the promoters of these genes have light-responsive regulatory elements such as AE-box, GT1-motif and TCT-motif. The expression of SiROS1d, SiROS1b, SiDML3a, SiDML5 and SiROS1a was higher in roots during ammonia, drought, nitrate, and urea treatments, suggesting that these genes may be involved in plant adversity stress. The analysis of cis-acting elements also showed that the promoters of these genes have multiple stress responses regulatory elements, such as MBS involved in drought induction and TC-rich repeats involved in defense and stress responses. These results suggest that the DNA-deMTase genes play regulatory roles at different stages of plant growth and development and in different tissues.

2.7. The Transcript Abundance of DNA-deMTase Genes in S. italica under Abiotic Stress

To clarify the potential roles of DNA-deMTase genes involved in abiotic stress, we used qRT-PCR to determine the abundance of DNA-deMTase genes in the transcript under cold, drought, NaCl, ABA, and MeJA stresses.
The transcript abundance of all DNA-deMTase genes changed significantly under cold treatment conditions (Figure 6). The expression patterns of SiDML5, SiDML4, SiDML3a, SiDML3b, SiROS1b, and SiROS1d were consistent with a tendency to increase and then decrease during the 24 h of cold stress. Among them, the expression of SiDML5 was significantly increased at 8 h, the expression of other genes was significantly increased at 4 h. The expression of SiROS1a was first decreased and then increased after the cold stress treatment, and then decreased again. The expression of SiROS1a was significantly increased at 8 h, and was significantly lower than 0 h at all the other times.
The transcript abundance of all DNA-deMTase genes changed significantly under drought treatment conditions (Figure 7). The expression patterns of SiDML5, SiDML4, SiDML3a, SiDML3b, SiROS1b, and SiROS1d were consistent with a tendency to increase and then decrease during the 24 h of drought stress. Among them, the expression of SiDML4 was significantly increased at 12 h, the expression of SiDML3b was significantly increased at 4 h, and the expression of other genes was significantly increased at 8 h. However, not all DNA-deMTase genes responded to drought stress, and the expression of SiROS1a was first decreased and then increased after the cold stress treatment, and then decreased again., but was always lower than 0 h.
The transcript abundance of all DNA-deMTase genes changed significantly under NaCl treatment conditions (Figure 8). The expression patterns of SiDML5, SiDML3a, and SiDML3b were consistent with a tendency to increase during the 24 h of NaCl stress. Among them, the expression of SiDML5 was significantly increased at 12 h, and the expression of other genes was significantly increased at 8 h. The expression patterns of SiDML4, SiROS1a, SiROS1b, and SiROS1d were consistent with a tendency to increase and then decrease during the 24 h of NaCl stress. The expression of SiDML4 was significantly increased at 4 h, and the expression of other genes was significantly increased at 12 h. The expression of SiDML4 was significantly increased at 4 h, the expression of SiROS1d was significantly increased at 8 h, and the expression of other genes was significantly increased at 12 h.
The transcript abundance of all DNA-deMTase genes changed significantly under ABA treatment conditions (Figure 9). The expression patterns of SiDML3a and SiROS1b were consistent with a tendency to increase and then decrease during the 24 h of ABA stress. Furthermore, the expression of SiDML3a and SiROS1b was significantly increased at 12 h; the expression of SiROS1a was first decreased and then increased after the ABA stress treatment, and then decreased again. However, not all DNA-deMTase genes responded to ABA stress. The expression of SiDML5, SiDML4, and SiROS1d was first decreased and then increased after the ABA stress treatment, and then decreased again, but was always lower than 0 h; the expression of SiDML3b was first decreased and then increased after the ABA stress treatment but was always lower than 0 h.
The transcript abundance of all DNA-deMTase genes changed significantly under MeJA treatment conditions (Figure 10). The expression patterns of SiDML3a, SiDML3b and SiROS1b were consistent with a tendency to increase and then decrease during the 24 h of MeJA stress. Furthermore, the expression of SiDML3a and SiROS1b was significantly increased at 12 h, the expression of SiDML3b was significantly increased at 4 h; the expression of SiDML5 was first increased and then decreased after the MeJA stress treatment, and then increased again; the expression of SiROS1a were consistent with a tendency to decrease and then increase during the 24 h of MeJA stress. However, not all DNA-deMTase genes responded to MeJA stress, and the expression of SiROS1d was first decreased and then increased after the MeJA stress treatment, and then decreased again, but was always lower than 0 h; the expression of SiDML4 was decreased after the MeJA stress treatment.
Temperature, drought, and salt stress are the main abiotic stresses in agriculture worldwide that threaten the growth and development of S. italica. ABA and MeJA as stress signaling factors for adversity can also improve plant stress tolerance. The above results suggest that DNA-deMTase genes from S. italica play important roles in abiotic stresses.

3. Discussion

3.1. Structural Features of DNA-deMTase in S. italica

With the release of the S. italica reference genome, the identification, classification, and prediction of gene family functions have gradually become new research hotspots. To elucidate the regulatory role of DNA demethylation in S. italica, we identified the DNA-deMTase genes at the whole genome level. The cis-acting elements, conserved motifs, phylogenetic relationships, and sequence features of these genes were also integratively analyzed. Currently, DNA-deMTase genes have been characterized in a variety of plants, such as A. thaliana, Solanum lycopersicum, Ricinus communis, Fragaria vesca, peanut [50,51,52,53,54], cotton and other plants. In this study, seven DNA-deMTases from S. italica were identified.
To identify factors in S. italica where specific epigenetic genes can regulate transcription, the promoter sequences of their putative cis-acting elements within 2000 bp upstream of the gene were analyzed. The results showed that the promoters of DNA-deMTases have four types of cis-actors: light response, phytohormone response, stress response, and growth and development. Among them, the stress-responsive elements are mainly involved in low-temperature stress, drought-induced, and stress responses. Similar CREs have been reported in the promoter regions of the DNA-deMTase gene of other plants such as Dendrobium officinale, tea plant, and kiwifruit [55,56,57].

3.2. Expression of the DNA-deMTase Gene in Abiotic Stress

Increasing evidence from recent studies suggests that DNA methylation plays an important role in regulating the stress response in plants [58,59], and change under different environmental stresses [60,61]. It was found that salt stress inhibited germination and seedling growth of S. italica seeds, and the restrain on radicel was bigger than that of on plantule, the relative radicel length was longer than relative plantule length [62]. Drought stress inhibited the rate of superoxide generation and hydrogen peroxide accumulation in S. italica seedlings, while decreasing MDA content [63]. SiARDP is a DREB-type transcription factor from S. italica. The transcript levels of SiARDP increased not only after drought, high salt, and low temperature stresses, but also after an ABA treatment of S. italica seedlings, indicating a key role in plant abiotic stress response [64]. Therefore, uncovering the mechanisms of tolerance to abiotic stress in plants and searching for related functional genes are essential tasks for successful molecular breeding. To adapt to environmental challenges, plants can make adjustments in response to various abiotic and biotic stresses [65,66]. Research has found that drought stress is associated with changes in genomic methylation of many genes, including those encoding transcription factors. In poplar, transcription factors that affect gene expression after drought treatment are influenced by methylated transposons [67]. Under salt stress, DNA methylation can control stress signals and induce plant stress responses by regulating membrane transporter genes, heavy metal transporter genes, and organic acid secretion genes. Under cold stress, genes involved in cellular metabolism, the stress response, the antioxidant system, the lysine metabolism pathway, and transcriptional regulation exhibit a correlation between methylation and their expression [68,69]. In addition, when plants are exposed to abiotic stress, the level of DNA methylation changes, and this effect can be repeatedly established [70].
In this study, we investigated the abundance of DNA-deMTase gene transcripts under cold, drought, salt, ABA, and MeJA stress. The results showed that the abundance of transcripts from all DNA-deMTase genes in S. italica under abiotic stress showed a dynamic trend at specific times from 0 to 48 h. These results suggest that gene families play different roles in abiotic stress, suggesting that demethylation of some key genes may be crucial for abiotic responses in plants.

4. Materials and Methods

4.1. Identification of DNA Demethylase Proteins in S. italica

It was shown that there are four DNA demethylase genes in the Arabidopsis thaliana dicotyledons plant: ROSl, DME, DML2 and DML3 [71], and eight DNA demethylase genes in the monocotyledon plant Oryza sativa: ROSla, ROSlb, ROSlc, ROS1d, DML3a, DML3b, DML4 and DML5 [72]. The amino acid sequences of DNA demethylases of A. thaliana, O. sativa, and Triticum aestivum [71,72,73] were used as queries in BLASTp (E value > 10−5) searches against the S. italica annotation database. To confirm the precision of these predicted genes, amino acid sequences were analyzed using CDD and Smart (http://smart.embl-heidelberg.de/, accessed on 20 August 2023), and sequences without conserved domains were excluded. S. italica genome information was downloaded from Phytozome v13 (https://phytozome-next.jgi.doe.gov/, accessed on 17 August 2023) [74]. The homologous DNA demethylase proteins in A. thaliana were obtained from TAIR (https://www.arabidopsis.org/, accessed on 20 August 2023); the homologous DNA demethylase proteins in O. sativa v7.0 were obtained from Phytozome v13 (https://phytozome-next.jgi.doe.gov/, accessed on 20 August 2023); the homologous DNA demethylase proteins in Triticum aestivum (IWGSC) were obtained from Ensembl Plants (https://plants.ensembl.org/index.html, accessed on 20 August 2023). The Plant-mPLoc 2.0 (http://www.csbio.sjtu.edu.cn/bioinf/plant-multi/, accessed on 25 September 2023) was used to predict the subcellular localization of these genes, and the ProtParam program (http://web.expasy.org/protparam/, accessed on 30 September 2023) in ExPASy was used to obtain the theoretical molecular weight (MW) and isoelectric point (pI), etc. of these proteins.

4.2. Phylogenetic Tree and Evolutionarily Conserved Protein Domain Analysis

The phylogenetic tree was established by 1000 bootstrap replications using Molecular Evolutionary Genetics Analysis (MEGA 7.0.26) software, based on the neighbor-joining (NJ) method [75,76]. DNA-deMTase were classified based on their phylogenetic relationship with the corresponding DNA-deMTase proteins from A. thaliana, O. sativa, and T. aestivum, and genes were named according to the chromosomal location of S. italica DNA-deMTase. The evolutionary tree was beautified by ITOL V6 (https://itol.embl.de, accessed on 20 August 2023) online software. The amino acid sequences and conserved domains of these proteins were analyzed using a conserved domain database (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi, accessed on 20 August 2023) [77].

4.3. Conserved Motifs and Gene Structure

The conserved motifs of all DNA-deMTase proteins were analyzed using Multiple Expectation Maximization for Motif Elicitation (MEME Suite 5.5.5, http://meme-suite.org/, accessed on 23 August 2023) by uploading the amino acid sequences following the MEME instructions. The conserved motif sites were visualized using TBtools v2.042 [78]. Gene structure Display Server 2.0 (https://gsds.gao-lab.org/, accessed on 23 August 2023) [79] was used to show the exon/intron structure of each DNA demethylase gene.

4.4. Analysis of Cis-Regulatory Elements and Chromosome Localization

The upstream sequences (2000 bp) of the start codon were retrieved from the S. italica genome database. Then the types, numbers, and functions of promoter cis-acting elements were analyzed by PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 28 August 2023) and visualized by TBtools v2.042 [80]. The locations of the genes identified on the chromosomes of S. italica were taken from the GFF (Graphics File Formats) file, which was downloaded from the genome database. The physical map of the chromosome was visualized by TBtools v2.042.

4.5. Expression of DNA Demethylases

FPKM (Fragments Per Kilobase of exon model per Million fragments mapped) expression value data (RNA-seq) of DNA demethylase genes of 7 tissues at different developmental stages of cereal were obtained from phytozome, and the FPKM values can reflect the gene expression level. The heatmap of gene expression was obtained using the TBtools v2.042.

4.6. Plant Material, Growth Condition and Abiotic Stress Treatment

The seeds of Setaria italica L. (Xiaomi) were obtained from Prof. Xingchun Wang, The Shanxi Agricultural University, China. Seeds were incubated in water for one night in the dark, evenly spread on beakers covered with gauze filled with water, and incubated in the dark for three days at 28 °C, replenishing water every six hours. Then, the uniformly grown cereal seedlings were selected and evenly placed in 96-well hydroponic boxes and cultured in Hoagland nutrient solution, which was changed every three days. The plants were placed in a plant growth chamber at 23 °C, bright fluorescent light of 150 µmol photons m−2 s−1, relative humidity of 65% and 16 h light: 8 h dark cycle.
For stress treatment, 30-day-old seedlings were subjected to six abiotic stresses, namely, cold (4 °C) [81], drought (20% PEG6000) [82], salt (150 Mm NaCl), ABA (100 µmol) [83], and MeJA (100 µmol) for a duration of 0, 4, 8, 12, and 24 h in triplicates. After 8 h of salt stress, the leaves began to wilt; at 12 h, more than half of the leaves had wilted; at 24 h, the plant wilted. No phenotypes were observed during the sampling of other treatments.

4.7. Total RNA Extraction and Transcript Abundance Analyzes of DNA-deMTase Genes

Total RNA from S. italica leaves was extracted using TransZol Up reagent (TransGen Biotech, Beijing, China). Total RNA was then reverse transcribed into the first strand cDNA using an Easyscript RT/RI Enzyme Mix (TransGen Biotech, Beijing, China). Quantitative real-time polymerase chain reaction (qRT-PCR) was performed with 2×PerfectStart Green qPCR SuperMix (TransGen Biotech, Beijing, China). The abundance of transcripts was calculated using the 2−ΔΔCt method [84], and all primers used for qRT-PCR were designed using NCBI, and both the qRT-PCR procedure and the reaction system were used in Table S1. All qRT-PCR analyzes were performed in three biological replicates. Statistical analyses were performed using IBM SPSS Statistics 21 software.

5. Conclusions

According to the reference genome of S. italica, seven DNA-deMTase genes were discovered and identified. Phylogenetic analysis demonstrated that the seven DNA-deMTase genes were divided into four categories, DML5, DML4, DML3, and ROS1. The DNA-deMTase promoter has a variety of stress-responsive, hormone signaling, and plant development-related regulatory elements, and DNA-deMTase genes play a regulatory role in different stages of plant growth and development and in different tissues. This indicates that the demethylase family plays an important role in the growth and development of S. italica. To gain further insight into the potential roles of DNA-deMTase involved in stress response, we analyzed seven DNA-deMTase genes from S. italica. Under cold and salt stress, all DNA-deMTase genes were induced with noticeably higher expression levels compared with the control (0 h). Under drought stress, SiDML5, SiDML4, SiDML3a, SiDML3b, SiROS1b, and SiROS1d were highly expressed compared with the control at 8 h and 12 h, while SiROS1a was lowly expressed. Under ABA treatments, the majority of SiDML3a, SiROS1a, and SiROS1b were upregulated. Among these genes, SiROS1b was the most preferentially expressed in response to ABA treatment. Additionally, SiDML5, SiDML3a, SiDML3b, SiROS1a, and SiROS1b were induced with higher expression levels under MeJA stress treatments. The results of this study lay the foundation for the discovery of functional genes and the further study of the stress resistance mechanisms of S. italica.

Supplementary Materials

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

Author Contributions

Conceptualization, K.L. and H.W.; Methodology, Z.S.; Software, Y.S.; Validation, X.W., Y.D. and J.L.; Formal Analysis, X.W.; Investigation, F.Z.; Resources, Z.S.; Data Curation, Y.D.; Writing—Original Draft Preparation, J.L.; Writing—Review & Editing, Y.S.; Visualization, X.W.; Supervision, Y.S.; Project Administration, Y.S.; Funding Acquisition, Z.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant numbers [32170350]), Fundamental Research Program of Shanxi Province (Grant numbers [202103021222005]), Research Project Supported by Shanxi Scholarship Council of China (Grant numbers [2022-024]), and the Open Project Program of Xinghuacun College of Shanxi University (Shanxi Institute of Brewing Technology and Industry) (Grant numbers [XCSXU-KF-202307]).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The materials and datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Goll, M.G.; Bestor, T.H. Eukaryotic cytosine methyltransferases. Ann. Rev. Biochem. 2005, 74, 481–514. [Google Scholar] [CrossRef] [PubMed]
  2. Feng, S.H.; Cokus, S.J.; Zhang, X.Y.; Chen, P.Y.; Bostick, M.; Goll, M.G.; Hetzel, J.; Jain, J.; Strauss, S.H.; Halpern, M.E.; et al. Conservation and divergence of methylation patterning in plants and animals. Proc. Natl. Acad. Sci. USA 2010, 107, 8689–8694. [Google Scholar] [CrossRef] [PubMed]
  3. Niederhuth, C.E.; Bewick, A.J.; Ji, L.X.; Alabady, M.S.; Kim, K.D.; Li, Q.; Rohr, N.A.; Rambani, A.; Burke, J.M.; Udall, J.A.; et al. Widespread natural variation of DNA methylation within angiosperms. Genome Biol. 2016, 17, 194. [Google Scholar] [CrossRef] [PubMed]
  4. Du, J.M.; Johnson, L.M.; Jacobsen, S.E.; Patel, D.J. DNA methylation pathways and their crosstalk with histone methylation. Nat. Rev. Mol. Cell Biol. 2015, 16, 519–532. [Google Scholar] [CrossRef] [PubMed]
  5. Omony, J.; Nussbaumer, T.; Gutzat, R. DNA methylation analysis in plants: Review of computational tools and future perspectives. Brief. Bioinf. 2020, 21, 906–918. [Google Scholar] [CrossRef] [PubMed]
  6. Ashapkin, V.V.; Kutueva, L.I.; Aleksandrushkina, N.I.; Vanyushin, B.F. Epigenetic Mechanisms of Plant Adaptation to Biotic and Abiotic Stresses. Int. J. Mol. Sci. 2020, 21, 7457. [Google Scholar] [CrossRef] [PubMed]
  7. Chen, W.W.; Kong, J.H.; Qin, C.; Yu, S.; Tan, J.J.; Chen, Y.R.; Wu, C.Q.; Wang, H.; Shi, Y.; Li, C.Y.; et al. Requirement of CHROMOMETHYLASE3 for somatic inheritance of the spontaneous tomato epimutation Colourless non-ripening. Sci. Rep. 2015, 5, 9192. [Google Scholar] [CrossRef] [PubMed]
  8. Manning, K.; Tör, M.; Poole, M.; Hong, Y.G.; Thompson, A.J.; King, G.J.; Giovannoni, J.J.; Seymour, G.B. A naturally occurring epigenetic mutation in a gene encoding an SBP-box transcription factor inhibits tomato fruit ripening. Nat. Genet. 2006, 38, 948–952. [Google Scholar] [CrossRef]
  9. Zhong, S.L.; Fei, Z.J.; Chen, Y.R.; Zheng, Y.; Huang, M.Y.; Vrebalov, J.; McQuinn, R.; Gapper, N.; Liu, B.; Xiang, J.; et al. Single-base resolution methylomes of tomato fruit development reveal epigenome modifications associated with ripening. Nat. Biotechnol. 2013, 31, 154–159. [Google Scholar] [CrossRef]
  10. Lee, S.; Choi, J.; Park, J.C.; Choi, J.Y.; Park, J.W.; Hong, C.P.; Choi, D.; Han, S.; Choi, K.U.; Roh, T.Y.; et al. DDM1-mediated gene body DNA methylation is associated with inducible activation of defense-related genes in Arabidopsis. Genome Biol. 2023, 24, 106. [Google Scholar] [CrossRef]
  11. Macdonald, W.A. Epigenetic mechanisms of genomic imprinting: Common themes in the regulation of imprinted regions in mammals, plants, and insects. Genet. Res. Int. 2012, 2012, 585024. [Google Scholar] [CrossRef] [PubMed]
  12. Choi, Y.; Gehring, M.; Johnson, L.; Hannon, M.; Harada, J.J.; Goldberg, R.B.; Jacobsen, S.E.; Fischer, R.L. DEMETER, a DNA glycosylase domain protein, is required for endosperm gene imprinting and seed viability in arabidopsis. Cell 2002, 110, 33–42. [Google Scholar] [CrossRef]
  13. Bartels, A.; Han, Q.; Nair, P.; Stacey, L.; Gaynier, H.; Mosley, M.; Huang, Q.Q.; Pearson, J.K.; Hsieh, T.F.; An, Y.Q.C.; et al. Dynamic DNA Methylation in Plant Growth and Development. Int. J. Mol. Sci. 2018, 19, 2144. [Google Scholar] [CrossRef]
  14. Chen, Y.Y.; Zhu, C.; Zhao, J.H.; Liu, T.; Gao, F.; Zhang, Y.C.; Duan, C.G. DNA methylation-dependent epigenetic regulation of Verticillium dahliae virulence in plants. aBIOTECH. 2023, 4, 185–201. [Google Scholar] [CrossRef]
  15. Zhang, L.; Wang, Y.; Zhang, X.; Zhang, M.; Han, D.; Qiu, C.; Han, Z. Dynamics of phytohormone and DNA methylation patterns changes during dormancy induction in strawberry (Fragaria × ananassa Duch.). Plant Cell Rep. 2012, 31, 155–165. [Google Scholar] [CrossRef]
  16. Deleris, A.; Halter, T.; Navarro, L. DNA Methylation and Demethylation in Plant Immunity. Annu. Rev. Phytopathol. 2016, 54, 579–603. [Google Scholar] [CrossRef] [PubMed]
  17. Abdulraheem, M.I.; Xiong, Y.; Moshood, A.Y.; Cadenas-Pliego, G.; Zhang, H.; Hu, J. Mechanisms of Plant Epigenetic Regulation in Response to Plant Stress: Recent Discoveries and Implications. Plants 2024, 13, 163. [Google Scholar] [CrossRef]
  18. Wang, R.H.; Xu, J.H. Genomic DNA methylation and histone methylation. Yi Chuan Hered. 2014, 36, 191–199. [Google Scholar] [CrossRef]
  19. Zhang, H.M.; Lang, Z.B.; Zhu, J.K. Dynamics and function of DNA methylation in plants. Nat. Rev. Mol. Cell Biol. 2018, 19, 489–506. [Google Scholar] [CrossRef]
  20. Cokus, S.J.; Feng, S.H.; Zhang, X.Y.; Chen, Z.G.; Merriman, B.; Haudenschild, C.D.; Pradhan, S.; Nelson, S.F.; Pellegrini, M.; Jacobsen, S.E. Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning. Nature 2008, 452, 215–219. [Google Scholar] [CrossRef]
  21. Kerr, L.; Kafetzopoulos, I.; Grima, R.; Sproul, D. Genome-wide single-molecule analysis of long-read DNA methylation reveals heterogeneous patterns at heterochromatin that reflect nucleosome organisation. PLoS Genet. 2023, 19, e1010958. [Google Scholar] [CrossRef] [PubMed]
  22. Fan, S.H.; Liu, H.F.; Liu, J.; Hua, W.; Xu, S.M.; Li, J. Systematic Analysis of the DNA Methylase and Demethylase Gene Families in Rapeseed (Brassica napus L.) and Their Expression Variations After Salt and Heat stresses. Int. J. Mol. Sci. 2020, 21, 953. [Google Scholar] [CrossRef] [PubMed]
  23. Jeltsch, A. Beyond Watson and Crick: DNA methylation and molecular enzymology of DNA methyltransferases. Chembiochem 2002, 3, 274–293. [Google Scholar] [CrossRef] [PubMed]
  24. Zhu, J.K. Active DNA demethylation mediated by DNA glycosylases. Ann. Rev. Genet. 2009, 43, 143–166. [Google Scholar] [CrossRef] [PubMed]
  25. Liu, T.K.; Li, Y.; Duan, W.K.; Huang, F.Y.; Hou, X.L. Cold acclimation alters DNA methylation patterns and confers tolerance to heat and increases growth rate in Brassica rapa. J. Exp. Bot. 2017, 68, 1213–1224. [Google Scholar] [CrossRef] [PubMed]
  26. Ganguly, D.R.; Crisp, P.A.; Eichten, S.R.; Pogson, B.J. The Arabidopsis DNA Methylome Is Stable under Transgenerational Drought Stress. Plant Physiol. 2017, 175, 1893–1912. [Google Scholar] [CrossRef] [PubMed]
  27. Ma, Y.Z.; Min, L.; Wang, M.J.; Wang, C.Z.; Zhao, Y.L.; Li, Y.Y.; Fang, Q.D.; Wu, Y.L.; Xie, S.; Ding, Y.H.; et al. Disrupted Genome Methylation in Response to High Temperature Has Distinct Affects on Microspore Abortion and Anther Indehiscence. Plant Cell. 2018, 30, 1387–1403. [Google Scholar] [CrossRef] [PubMed]
  28. Chen, X.; Schï Nberger, B.; Menz, J.; Ludewig, U. Plasticity of DNA methylation and gene expression under zinc deficiency in Arabidopsis roots. Plant Cell Physiol. 2018, 59, 1790–1802. [Google Scholar] [CrossRef] [PubMed]
  29. Chen, R.; Li, M.; Zhang, H.Y.; Duan, L.J.; Sun, X.J.; Jiang, Q.Y.; Zhang, H.; Hu, Z. Continuous salt stress-induced long non-coding RNAs and DNA methylation patterns in soybean roots. BMC Genom. 2019, 20, 730. [Google Scholar] [CrossRef]
  30. Marfil, C.; Ibañez, V.; Alonso, R.; Varela, A.; Bottini, R.; Masuelli, R.; Fontana, A.; Berli, F. Changes in grapevine DNA methylation and polyphenols content induced by solar ultraviolet-B radiation, water deficit and abscisic acid spray treatments. Plant Physiol. Biochem. PPB 2019, 135, 287–294. [Google Scholar] [CrossRef]
  31. Song, Q.X.; Zhang, T.Z.; Stelly, D.M.; Chen, Z.J. Epigenomic and functional analyses reveal roles of epialleles in the loss of photoperiod sensitivity during domestication of allotetraploid cottons. Genome Biol. 2017, 18, 99. [Google Scholar] [CrossRef]
  32. Finnegan, E.J.; Peacock, W.J.; Dennis, E.S. DNA methylation, a key regulator of plant development and other processes. Curr. Opin. Genet. Dev. 2000, 10, 217–223. [Google Scholar] [CrossRef]
  33. Yang, L.; Lang, C.; Wu, Y.; Meng, D.; Yang, T.; Li, D.; Jin, T.; Zhou, X. ROS1-mediated decrease in DNA methylation and increase in expression of defense genes and stress response genes in Arabidopsis thaliana due to abiotic stresses. BMC Plant Biol. 2022, 22, 104. [Google Scholar] [CrossRef] [PubMed]
  34. Garg, R.; Narayana Chevala, V.; Shankar, R.; Jain, M. Divergent DNA methylation patterns associated with gene expression in rice cultivars with contrasting drought and salinity stress response. Sci. Rep. 2015, 5, 14922. [Google Scholar] [CrossRef]
  35. Chen, Y.; Zhang, W.F.; Ren, H.Z.; Xiong, F.; Zhang, H.J.; Yao, L.N.; Liu, Y.; Wang, L.; Wang, X.C.; Yang, Y.J.; et al. Genome-wide Investigation and Expression Analysis of DNA Demethylase Genes in Tea Plant (Camellia Sinensis). J. Tea Sci. 2021, 41, 28–39. [Google Scholar] [CrossRef]
  36. Liu, Q.X.; Xue, Q.Z.; Xu, J.H. Phylogenetic analysis of DNA demethylase genes in angiosperm. Hereditas 2014, 36, 276–285. [Google Scholar] [CrossRef] [PubMed]
  37. Hu, H.; Mauro-Herrera, M.; Doust, A.N. Domestication and Improvement in the Model C4 Grass, Setaria. Front. Plant Sci. 2018, 9, 719. [Google Scholar] [CrossRef]
  38. Lu, H.; Zhang, J.; Liu, K.B.; Wu, N.; Li, Y.; Zhou, K.; Ye, M.; Zhang, T.; Zhang, H.; Yang, X.; et al. Earliest domestication of common millet (Panicum miliaceum) in East Asia extended to 10,000 years ago. Proc. Natl. Acad. Sci. USA 2009, 106, 7367–7372. [Google Scholar] [CrossRef] [PubMed]
  39. Li, P.H.; Brutnell, T.P. Setaria viridis and Setaria italica, model genetic systems for the Panicoid grasses. J. Exp. Bot. 2011, 62, 3031–3037. [Google Scholar] [CrossRef]
  40. Krishnamurthy, L.; Upadhyaya, H.D.; Gowda, L.L.; Kashiwagi, J.; Purushothaman, R.; Sube, S.; Vadez, V. Large variation for salinity tolerance in the core collection of foxtail millet (Setaria italica (L.) P. Beauv.) germplasm. Crop Pasture Sci. 2013, 65, 353–361. [Google Scholar] [CrossRef]
  41. Sudhakar, C.; Veeranagamallaiah, G.; Nareshkumar, A.; Sudhakarbabu, O.; Sivakumar, M.; Pandurangaiah, M.; Kiranmai, K.; Lokesh, U. Polyamine metabolism influences antioxidant defense mechanism in foxtail millet (Setaria italica L.) cultivars with different salinity tolerance. Plant Cell Rep. 2015, 34, 141–156. [Google Scholar] [CrossRef] [PubMed]
  42. Genger, R.K.; Kovac, K.A.; Dennis, E.S.; Peacock, W.J.; Finnegan, E.J. Multiple DNA methyltransferase genes in Arabidopsis thaliana. Plant Mol. Biol. 1999, 41, 269–278. [Google Scholar] [CrossRef] [PubMed]
  43. Ahmad, F.; Huang, X.; Lan, H.X.; Huma, T.; Bao, Y.M.; Huang, J.; Zhang, H.S. Comprehensive gene expression analysis of the DNA (cytosine-5) methyltransferase family in rice (Oryza sativa L.). Genet. Mol. Res. GMR 2014, 13, 5159–5172. [Google Scholar] [CrossRef] [PubMed]
  44. Palmgren, G.; Mattsson, O.; Okkels, F.T. Specific Levels of DNA Methylation in Various Tissues, Cell Lines, and Cell Types of Daucus carota. Plant Physiol. 1991, 95, 174–178. [Google Scholar] [CrossRef] [PubMed]
  45. Castiglione, M.R.; Venora, G.; Ravalli, C.; Stoilov, L.; Gecheff, K.; Cremonini, R. DNA methylation and chromosomal rearrangements in reconstructed karyotypes of Hordeum vulgare L. Protoplasma 2008, 232, 215–222. [Google Scholar] [CrossRef] [PubMed]
  46. Qian, Y.X.; Xi, Y.L.; Cheng, B.J.; Zhu, S.W. Genome-wide identification and expression profiling of DNA methyltransferase gene family in maize. Plant Cell Rep. 2014, 33, 1661–1672. [Google Scholar] [CrossRef] [PubMed]
  47. Teyssier, E.; Bernacchia, G.; Maury, S.; How Kit, A.; Stammitti-Bert, L.; Rolin, D.; Gallusci, P. Tissue dependent variations of DNA methylation and endoreduplication levels during tomato fruit development and ripening. Planta 2008, 228, 391–399. [Google Scholar] [CrossRef]
  48. Nakano, Y.; Steward, N.; Sekine, M.; Kusano, T.; Sano, H. A tobacco NtMET1 cDNA encoding a DNA methyltransferase: Molecular characterization and abnormal phenotypes of transgenic tobacco plants. Plant Cell Physiol. 2000, 41, 448–457. [Google Scholar] [CrossRef]
  49. Garg, R.; Kumari, R.; Tiwari, S.; Goyal, S. Genomic survey, gene expression analysis and structural modeling suggest diverse roles of DNA methyltransferases in legumes. PLoS ONE 2014, 9, e88947. [Google Scholar] [CrossRef]
  50. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  51. Cao, D.Y.; Ju, Z.J.; Gao, C.; Mei, X.H.; Fu, D.Q.; Zhu, H.L.; Luo, Y.B.; Zhu, B.Z. Genome-wide identification of cytosine-5 DNA methyltransferases and demethylases in Solanum lycopersicum. Gene 2014, 550, 230–237. [Google Scholar] [CrossRef] [PubMed]
  52. Victoria, D.; Aliki, K.; Venetia, K.; Georgios, M.; Zoe, H. Spatial and temporal expression of cytosine-5 DNA methyltransferase and DNA demethylase gene families of the Ricinus communis during seed development and drought stress. Plant Growth Regul. 2018, 84, 81–94. [Google Scholar] [CrossRef]
  53. Gu, T.T.; Ren, S.; Wang, Y.H.; Han, Y.H.; Li, Y. Characterization of DNA methyltransferase and demethylase genes in Fragaria vesca. Mol. Genet. Genom. MGG 2016, 291, 1333–1345. [Google Scholar] [CrossRef] [PubMed]
  54. Wang, P.F.; Gao, C.; Bian, X.T.; Zhao, S.Z.; Zhao, C.Z.; Xia, H.; Song, H.; Hou, L.; Wan, S.B.; Wang, X.J. Genome-Wide Identification and Comparative Analysis of Cytosine-5 DNA Methyltransferase and Demethylase Families in Wild and Cultivated Peanut. Front. Plant Sci. 2016, 7, 7. [Google Scholar] [CrossRef] [PubMed]
  55. Yu, Z.M.; Zhang, G.H.; Teixeira da Silva, J.A.; Li, M.Z.; Zhao, C.H.; He, C.M.; Si, C.; Zhang, M.Z.; Duan, J. Genome-wide identification and analysis of DNA methyltransferase and demethylase gene families in Dendrobium officinale reveal their potential functions in polysaccharide accumulation. BMC Plant Biol. 2021, 21, 21. [Google Scholar] [CrossRef] [PubMed]
  56. Zhu, C.; Zhang, S.T.; Zhou, C.Z.; Chen, L.; Fu, H.F.; Li, X.Z.; Lin, Y.L.; Lai, Z.X.; Guo, Y.Q. Genome-wide investigation and transcriptional analysis of cytosine-5 DNA methyltransferase and DNA demethylase gene families in tea plant (Camellia sinensis) under abiotic stress and withering processing. PeerJ. 2020, 8, e8432. [Google Scholar] [CrossRef] [PubMed]
  57. Zhang, Y.X.; He, X.Q.; Zhao, H.C.; Xu, W.C.; Deng, H.; Wang, H.; Wang, S.Y.; Su, D.; Zheng, Z.L.; Yang, B.; et al. Genome-Wide Identification of DNA Methylases and Demethylases in Kiwifruit (Actinidia chinensis). Front. Plant Sci. 2020, 11, 514993. [Google Scholar] [CrossRef] [PubMed]
  58. Sanchez, D.H.; Paszkowski, J. Heat-induced release of epigenetic silencing reveals the concealed role of an imprinted plant gene. PLoS Genet. 2014, 10, e1004806. [Google Scholar] [CrossRef]
  59. Xu, R.; Wang, Y.H.; Zheng, H.; Lu, W.; Wu, C.G.; Huang, J.G.; Yan, K.; Yang, G.D.; Zheng, C.C. Salt-induced transcription factor MYB74 is regulated by the RNA-directed DNA methylation pathway in Arabidopsis. J. Exp. Bot. 2015, 66, 5997–6008. [Google Scholar] [CrossRef]
  60. Chang, Y.N.; Zhu, C.; Jiang, J.; Zhang, H.; Zhu, J.K.; Duan, C.G. Epigenetic regulation in plant abiotic stress responses. J. Integr. Plant Biol. 2020, 62, 563–580. [Google Scholar] [CrossRef]
  61. Wang, L.; Wang, L.; Tan, M.; Wang, L.; Zhao, W.; You, J.; Wang, L.; Yan, X.; Wang, W. The pattern of alternative splicing and DNA methylation alteration and their interaction in linseed (Linum usitatissimum L.) response to repeated drought stresses. Biol. Res. 2023, 56, 12. [Google Scholar] [CrossRef]
  62. Cui, X.G.; Shi, L.R. Effect of Salt Stress on Different Seed Germination and Seedling Growth of Setaria italica (L.). Heilongjiang Agric. Sci. 2011, 14–16. [Google Scholar]
  63. Pei, S.S.; Wang, Y.G.; Yin, M.Q.; Yuan, X.Y.; Huang, M.J.; Wen, Y.Y.; Zhao, J. Effect of Compound Seed Soaking on Seed Germination and Seedling Growth in Millet(Setaria italica)under Drought Stress. J. Shanxi Agricult. Sci. 2013, 41, 676–679. [Google Scholar]
  64. Li, C.; Yue, J.; Wu, X.; Xu, C.; Yu, J. An ABA-responsive DRE-binding protein gene from Setaria italica, SiARDP, the target gene of SiAREB, plays a critical role under drought stress. J. Exp. Bot. 2014, 65, 5415–5427. [Google Scholar] [CrossRef] [PubMed]
  65. Bilichak, A.; Ilnystkyy, Y.; Hollunder, J.; Kovalchuk, I. The progeny of Arabidopsis thaliana plants exposed to salt exhibit changes in DNA methylation, histone modifications and gene expression. PLoS ONE 2012, 7, e30515. [Google Scholar] [CrossRef]
  66. Matzke, M.A.; Mosher, R.A. RNA-directed DNA methylation: An epigenetic pathway of increasing complexity. Nat. Rev. Genet. 2014, 15, 394–408. [Google Scholar] [CrossRef] [PubMed]
  67. Liang, D.; Zhang, Z.J.; Wu, H.L.; Huang, C.Y.; Shuai, P.; Ye, C.Y.; Tang, S.; Wang, Y.J.; Yang, L.; Wang, J.; et al. Single-base-resolution methylomes of Populus trichocarpa reveal the association between DNA methylation and drought stress. BMC Genet. 2014, 15 (Suppl 1), S9. [Google Scholar] [CrossRef] [PubMed]
  68. Rakei, A.; Maali-Amiri, R.; Zeinali, H.; Ranjbar, M. DNA methylation and physio-biochemical analysis of chickpea in response to cold stress. Protoplasma 2016, 253, 61–76. [Google Scholar] [CrossRef]
  69. Song, Y.; Jia, Z.; Hou, Y.; Ma, X.; Li, L.; Jin, X.; An, L. Roles of DNA Methylation in Cold Priming in Tartary Buckwheat. Front. Plant Sci. 2020, 11, 608540. [Google Scholar] [CrossRef]
  70. Sun, M.; Yang, Z.; Liu, L.; Duan, L. DNA Methylation in Plant Responses and Adaption to Abiotic Stresses. Int. J. Mol. Sci. 2022, 23, 6910. [Google Scholar] [CrossRef]
  71. Yu, A.; Lepère, G.; Jay, F.; Wang, J.; Bapaume, L.; Wang, Y.; Abraham, A.L.; Penterman, J.; Fischer, R.L.; Voinnet, O.; et al. Dynamics and biological relevance of DNA demethylation in Arabidopsis antibacterial defense. Proc. Natl. Acad. Sci. USA 2013, 110, 2389–2394. [Google Scholar] [CrossRef] [PubMed]
  72. Zemach, A.; Kim, M.Y.; Silva, P.; Rodrigues, J.A.; Dotson, B.; Brooks, M.D.; Zilberman, D. Local DNA hypomethylation activates genes in rice endosperm. Proc. Natl. Acad. Sci. USA 2010, 107, 18729–18734. [Google Scholar] [CrossRef] [PubMed]
  73. Jiang, Z.N.; Liu, D.T.; Zhang, X.; Jiang, W.; Wang, L.; Li, D.S.; Cheng, X.M.; Gao, D.R. Genome-Wide Investigation and Expression Analysis of DNA Demethylase Genes in Wheat (Triticum aestivum) during the Grain Development. J. Triticeae Crops. 2023, 43, 1361–1371. [Google Scholar] [CrossRef]
  74. Goodstein, D.M.; Shu, S.; Howson, R.; Neupane, R.; Hayes, R.D.; Fazo, J.; Mitros, T.; Dirks, W.; Hellsten, U.; Putnam, N.; et al. Phytozome: A comparative platform for green plant genomics. Nucleic Acids Res. 2012, 40, D1178–D1186. [Google Scholar] [CrossRef] [PubMed]
  75. 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]
  76. 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] [PubMed]
  77. 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] [PubMed]
  78. Chen, C.J.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.H.; Xia, R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef] [PubMed]
  79. Kayum, M.A.; Park, J.I.; Nath, U.K.; Saha, G.; Biswas, M.K.; Kim, H.T.; Nou, I.S. Genome-wide characterization and expression profiling of PDI family gene reveals function as abiotic and biotic stress tolerance in Chinese cabbage (Brassica rapa ssp. pekinensis). BMC Genom. 2017, 18, 885. [Google Scholar] [CrossRef]
  80. Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van de Peer, Y.; Rouzé, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30, 325–327. [Google Scholar] [CrossRef]
  81. Fan, Y.; Lai, D.; Yang, H.; Xue, G.; He, A.; Chen, L.; Feng, L.; Ruan, J.; Xiang, D.; Yan, J.; et al. Genome-wide identification and expression analysis of the bHLH transcription factor family and its response to abiotic stress in foxtail millet (Setaria italica L.). BMC Genom. 2021, 22, 778. [Google Scholar] [CrossRef] [PubMed]
  82. Pan, Y.; Li, J.; Jiao, L.; Li, C.; Zhu, D.; Yu, J. A Non-specific Setaria italica Lipid Transfer Protein Gene Plays a Critical Role under Abiotic Stress. Front. Plant Sci. 2016, 7, 1752. [Google Scholar] [CrossRef] [PubMed]
  83. Zhao, J.; Yu, A.; Du, Y.; Wang, G.; Li, Y.; Zhao, G.; Wang, X.; Zhang, W.; Cheng, K.; Liu, X.; et al. Foxtail millet (Setaria italica (L.) P. Beauv) CIPKs are responsive to ABA and abiotic stresses. PLoS ONE 2019, 14, e0225091. [Google Scholar] [CrossRef]
  84. Ogneva, Z.V.; Dubrovina, A.S.; Kiselev, K.V. Age-associated alterations in DNA methylation and expression of methyltransferase and demethylase genes in Arabidopsis thaliana. Biol. Plant. 2016, 60, 628–634. [Google Scholar] [CrossRef]
Ijms 25 04464 g001
Figure 1. Phylogenetic trees (a) and conserved domain (b) analysis of DNA demethylase (DNA-deMTase). The red dots in the figure are the DNA-deMTase genes from S. italica.
Figure 1. Phylogenetic trees (a) and conserved domain (b) analysis of DNA demethylase (DNA-deMTase). The red dots in the figure are the DNA-deMTase genes from S. italica.
Ijms 25 04464 g001
Ijms 25 04464 g002
Figure 2. Conserved motifs (a) and gene structure (b) of DNA-deMTase genes in Foxtail millet (Setaria italica L.).
Figure 2. Conserved motifs (a) and gene structure (b) of DNA-deMTase genes in Foxtail millet (Setaria italica L.).
Ijms 25 04464 g002
Ijms 25 04464 g003
Figure 3. Cis-regulatory elements (CREs) of DNA-deMTase genes in S. italica. (a): Heatmap showing the number of CREs belonging to the following four categories per DNA-deMTase gene. (b): The distribution of different CREs in the promoter region (2000 bp) of DNA-deMTase genes in S. italica.
Figure 3. Cis-regulatory elements (CREs) of DNA-deMTase genes in S. italica. (a): Heatmap showing the number of CREs belonging to the following four categories per DNA-deMTase gene. (b): The distribution of different CREs in the promoter region (2000 bp) of DNA-deMTase genes in S. italica.
Ijms 25 04464 g003
Ijms 25 04464 g004
Figure 4. Chromosomal locations of S. italica DNA-deMTase genes on 6 chromosomes. The gene names on the right side of each chromosome correspond to the approximate locations of each DNA-deMTase gene.
Figure 4. Chromosomal locations of S. italica DNA-deMTase genes on 6 chromosomes. The gene names on the right side of each chromosome correspond to the approximate locations of each DNA-deMTase gene.
Ijms 25 04464 g004
Ijms 25 04464 g005
Figure 5. Heat map of S. italica DNA-deMTase families’ genes expression patterns across tissues. Color scale bar at the right of heat map represents log2-transformed FPKM (Fragments Per Kilobase of exon per Million fragments mapped) values for each gene, with warmer colors denoting higher expression.
Figure 5. Heat map of S. italica DNA-deMTase families’ genes expression patterns across tissues. Color scale bar at the right of heat map represents log2-transformed FPKM (Fragments Per Kilobase of exon per Million fragments mapped) values for each gene, with warmer colors denoting higher expression.
Ijms 25 04464 g005
Ijms 25 04464 g006
Figure 6. Transcript abundance of SiDML5 (a), SiDML4 (b), SiDML3a (c), SiDML3b (d), SiROS1a (e), SiROS1b (f), SiROS1d (g) under cold stress. Data are the mean ± SE (n = 3). Different letters above bars indicate statistically significant differences between treatments (Tukey’s multiple range test, p ≤ 0.05).
Figure 6. Transcript abundance of SiDML5 (a), SiDML4 (b), SiDML3a (c), SiDML3b (d), SiROS1a (e), SiROS1b (f), SiROS1d (g) under cold stress. Data are the mean ± SE (n = 3). Different letters above bars indicate statistically significant differences between treatments (Tukey’s multiple range test, p ≤ 0.05).
Ijms 25 04464 g006
Ijms 25 04464 g007
Figure 7. Transcript abundance of SiDML5 (a), SiDML4 (b), SiDML3a (c), SiDML3b (d), SiROS1a (e), SiROS1b (f), SiROS1d (g) under drought stress. Data are the mean ± SE (n = 3). Different letters above bars indicate statistically significant differences between treatments (Tukey’s multiple range test, p ≤ 0.05).
Figure 7. Transcript abundance of SiDML5 (a), SiDML4 (b), SiDML3a (c), SiDML3b (d), SiROS1a (e), SiROS1b (f), SiROS1d (g) under drought stress. Data are the mean ± SE (n = 3). Different letters above bars indicate statistically significant differences between treatments (Tukey’s multiple range test, p ≤ 0.05).
Ijms 25 04464 g007
Ijms 25 04464 g008
Figure 8. Transcript abundance of SiDML5 (a), SiDML4 (b), SiDML3a (c), SiDML3b (d), SiROS1a (e), SiROS1b (f), SiROS1d (g) under NaCl stress. Data are the mean ± SE (n = 3). Different letters above bars indicate statistically significant differences between treatments (Tukey’s multiple range test, p ≤ 0.05).
Figure 8. Transcript abundance of SiDML5 (a), SiDML4 (b), SiDML3a (c), SiDML3b (d), SiROS1a (e), SiROS1b (f), SiROS1d (g) under NaCl stress. Data are the mean ± SE (n = 3). Different letters above bars indicate statistically significant differences between treatments (Tukey’s multiple range test, p ≤ 0.05).
Ijms 25 04464 g008
Ijms 25 04464 g009
Figure 9. Transcript abundance of SiDML5 (a), SiDML4 (b), SiDML3a (c), SiDML3b (d), SiROS1a (e), SiROS1b (f), SiROS1d (g) under ABA stress. Data are the mean ± SE (n = 3). Different letters above bars indicate statistically significant differences between treatments (Tukey’s multiple range test, p ≤ 0.05).
Figure 9. Transcript abundance of SiDML5 (a), SiDML4 (b), SiDML3a (c), SiDML3b (d), SiROS1a (e), SiROS1b (f), SiROS1d (g) under ABA stress. Data are the mean ± SE (n = 3). Different letters above bars indicate statistically significant differences between treatments (Tukey’s multiple range test, p ≤ 0.05).
Ijms 25 04464 g009
Ijms 25 04464 g010
Figure 10. Transcript abundance of SiDML5 (a), SiDML4 (b), SiDML3a (c), SiDML3b (d), SiROS1a (e), SiROS1b (f), SiROS1d (g) under MeJA stress. Data are the mean ± SE (n = 3). Different letters above bars indicate statistically significant differences between treatments (Tukey’s multiple range test, p ≤ 0.05).
Figure 10. Transcript abundance of SiDML5 (a), SiDML4 (b), SiDML3a (c), SiDML3b (d), SiROS1a (e), SiROS1b (f), SiROS1d (g) under MeJA stress. Data are the mean ± SE (n = 3). Different letters above bars indicate statistically significant differences between treatments (Tukey’s multiple range test, p ≤ 0.05).
Ijms 25 04464 g010
Table 1. Characteristic features of DNA demethylase (DNA-deMTase) gene family in Foxtail millet (Setaria italica L.).
Table 1. Characteristic features of DNA demethylase (DNA-deMTase) gene family in Foxtail millet (Setaria italica L.).
Gene NameGene IDNo. of IntronsLength (Aa)MW
(kDa)		pIInstability IndexGRAVY
Value		Subcellular Localization
SiDML5Seita.8G1043001035839.749.6741.68−0.455Chloroplas. Nucleus.
SiDML3aSeita.1G164500161060118.047.8449.08−0.549Nucleus.
SiROS1bSeita.3G223900191663185.596.2056.29−0.640Nucleus.
SiROS1aSeita.2G152900172082230.616.7946.03−0.605Nucleus.
SiROS1dSeita.1G164200161954216.946.4547.77−0.672Nucleus.
SiDML3bSeita.9G1855001157865.775.9742.44−0.665Nucleus.
SiDML4Seita.4G105300428631.529.1056.21−0.300Nucleus.
Aa: amino acid; MW: Molecular weight; pI: theoretical isoelectric point; GRAVY: Grand average of hydropathy.
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

Sun, Y.; Wang, X.; Di, Y.; Li, J.; Li, K.; Wei, H.; Zhang, F.; Su, Z. Systematic Analysis of DNA Demethylase Gene Families in Foxtail Millet (Setaria italica L.) and Their Expression Variations after Abiotic Stresses. Int. J. Mol. Sci. 2024, 25, 4464. https://doi.org/10.3390/ijms25084464

AMA Style

Sun Y, Wang X, Di Y, Li J, Li K, Wei H, Zhang F, Su Z. Systematic Analysis of DNA Demethylase Gene Families in Foxtail Millet (Setaria italica L.) and Their Expression Variations after Abiotic Stresses. International Journal of Molecular Sciences. 2024; 25(8):4464. https://doi.org/10.3390/ijms25084464

Chicago/Turabian Style

Sun, Yingying, Xin Wang, Yunfei Di, Jinxiu Li, Keyu Li, Huanhuan Wei, Fan Zhang, and Zhenxia Su. 2024. "Systematic Analysis of DNA Demethylase Gene Families in Foxtail Millet (Setaria italica L.) and Their Expression Variations after Abiotic Stresses" International Journal of Molecular Sciences 25, no. 8: 4464. https://doi.org/10.3390/ijms25084464

APA Style

Sun, Y., Wang, X., Di, Y., Li, J., Li, K., Wei, H., Zhang, F., & Su, Z. (2024). Systematic Analysis of DNA Demethylase Gene Families in Foxtail Millet (Setaria italica L.) and Their Expression Variations after Abiotic Stresses. International Journal of Molecular Sciences, 25(8), 4464. https://doi.org/10.3390/ijms25084464

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