Next Article in Journal
Applying UHPLC-HRAM MS/MS Method to Assess Host Cell Protein Clearance during the Purification Process Development of Therapeutic mAbs
Next Article in Special Issue
Identification, Cloning, and Characterization of Two Acupuncture-Injury-Inducing Promoters in Rice
Previous Article in Journal
A Transcriptomic Analysis of Laryngeal Dysplasia
Previous Article in Special Issue
Pangenome Identification and Analysis of Terpene Synthase Gene Family Members in Gossypium
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring the Roles of the Swi2/Snf2 Gene Family in Maize Abiotic Stress Responses

1
Institute of Agricultural Biotechnology/Jilin Provincial Key Laboratory of Agricultural Biotechnology, Jilin Academy of Agricultural Sciences (Northeast Innovation Center of Agricultural Science and Technology in China), Changchun 130033, China
2
College of Agriculture, Yanbian University, Yanji 133002, China
3
College of Agronomy, Jilin Agricultural University, Changchun 130118, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2024, 25(17), 9686; https://doi.org/10.3390/ijms25179686
Submission received: 17 July 2024 / Revised: 4 September 2024 / Accepted: 4 September 2024 / Published: 7 September 2024
(This article belongs to the Special Issue Physiology and Molecular Biology of Plant Stress Tolerance)

Abstract

:
The maize Snf2 gene family plays a crucial role in chromatin remodeling and response to environmental stresses. In this study, we identified and analyzed 35 members of the maize Snf2 gene family (ZmCHR1 to ZmCHR35) using the Ensembl Plants database. Each protein contained conserved SNF2-N and Helicase-C domains. Phylogenetic analysis revealed six groups among the Snf2 proteins, with an uneven distribution across subfamilies. Physicochemical analysis indicated that the Snf2 proteins are hydrophilic, with varied amino acid lengths, isoelectric points, and molecular weights, and are predominantly localized in the nucleus. Chromosomal mapping showed that these genes are distributed across all ten maize chromosomes. Gene structure analysis revealed diverse exon–intron arrangements, while motif analysis identified 20 conserved motifs. Collinearity analysis highlighted gene duplication events, suggesting purifying selection. Cis-regulatory element analysis suggested involvement in abiotic and biotic stress responses. Expression analysis indicated tissue-specific expression patterns and differential expression under various stress conditions. Specifically, qRT-PCR validation under drought stress showed that certain Snf2 genes were upregulated at 12 h and downregulated at 24 h, revealing potential roles in drought tolerance. These findings provide a foundation for further exploration of the functional roles of the maize Snf2 gene family in development and stress responses.

1. Introduction

Maize is a cornerstone of global food security and economic development, and it ranks as one of the three major crops in Asia. The variability of the global environment poses significant threats to maize production worldwide, with various stresses including biotic and abiotic challenges [1,2]. Current research confirms that abiotic stresses such as extreme temperatures, high salinity, and drought are among the primary contributors to global maize yield losses [3]. It is estimated that a temperature increase of just one degree Celsius could lead to a 7.4% reduction in global maize production [4]. Consequently, elucidating the molecular mechanisms by which maize responds to abiotic stress is imperative.
In this context, transcription factors have garnered extensive attention in the scientific community. A central question in analyzing eukaryotic transcription factors is how these DNA-binding proteins interact with the chromatin template [5]. Research has shown that Swi2/Snf2 and related proteins can disrupt nucleosome stability, utilizing ATP to alter chromatin structure and nucleosome positioning, thereby facilitating the binding of transcription factors to chromatin. This is crucial for transcription regulation, DNA replication, and DNA damage repair [6,7,8]. Genetic studies of transcription regulation in Saccharomyces cerevisiae led to the discovery of various Swi (yeast mating type switching) and Snf (sucrose non-fermenting) genes [9]. Notably, the Swi2 gene in Swi and the Snf2 gene in Snf were found to be the same; hence, the gene is referred to as Swi2/Snf2. Plant genomes contain numerous Snf2-family proteins, some of which have been demonstrated to act as key developmental regulators during different stages in model plants such as Arabidopsis and rice [10]. Despite their crucial roles, studies of these proteins in maize have been relatively sparse.
As autotrophic organisms, plants can rapidly adjust their gene expression under environmental stress to ensure adaptability, genomic integrity, and enhanced future stress tolerance [11,12]. Although DNA methylation, histone variants, and modifications play significant roles in plant stress responses, the role of Snf2 chromatin remodeling pathways remains relatively understudied. Recent research highlights the prominence of Snf2 genes under various abiotic stresses across multiple plant species [13]. In Arabidopsis, the Snf2 family member SPLAYED (SYD) regulates the expression of many genes responsive to drought, salt, and temperature stress by remodeling chromatin structure [14]; another member, BRAHMA (BRM), participates in the abscisic acid (ABA) signaling pathway to modulate plant responses to drought stress [15]. In rice, the Snf2 family member OsCHR4 enhances drought resistance by increasing the waxy content of the leaf surface cuticle [16]. Transgenic Arabidopsis overexpressing AtCHR12, an Snf2 family member, exhibit halted primary shoot growth and slowed main stem growth under drought and high-temperature stress [17]. Similarly, the tomato SlCHR1 gene, similar to Arabidopsis CHR12, regulates growth retardation under salt stress while enhancing plant salt tolerance by maintaining chromatin structure stability [13]. Furthermore, maize ZmCHB101 mutants show heightened sensitivity to salt stress due to chromatin alterations at the promoter, leading to disrupted expression of stress-responsive genes [18]. Additionally, Snf2 factors are involved in hormone pathways, including cytokinin, auxin, and gibberellin signaling [19]. Arabidopsis BRM and PKL regulate the distribution of auxin and cytokinin responses, respectively, through PIN and ARR genes [20]. Thus, Snf2 chromatin remodeling factors have emerged as critical regulators of stress responses and developmental processes in various plant species.
The Snf2 protein family is characterized by conserved SNF2_N and Helicase_C domains [21]. Based on similarities in helicase-like regions and functional variations, the Snf2 gene family is well represented in plant genomes. Arabidopsis has 41 members [22], rice has 40 [23], soybean has 66 [7], tomato has 45 [24], and barley has 38 members [25]. However, systematic identification and expression analysis under abiotic stress conditions for these genes in maize are relatively scarce. This study aims to thoroughly investigate the expression patterns of the maize Swi2/Snf2 gene family members under abiotic stress conditions, providing insights for the cloning and characterization of maize-related functional genes and molecular marker-assisted breeding.

2. Results

2.1. Identification of the Maize Snf2 Gene Family

Using the Ensembl Plants database, 35 members of the maize Snf2 gene family were identified and systematically named ZmCHR1 to ZmCHR35 (Table S1), following the nomenclature currently used for Arabidopsis and rice and based on their physical chromosomal locations. All derived Snf2 proteins were found to contain the conserved SNF2-N domain and Helicase-C domain. To explore the evolutionary relationships between maize and Arabidopsis Snf2 proteins, phylogenetic analysis was conducted using the maximum likelihood method, incorporating 35 maize Snf2 genes and 41 Arabidopsis Snf2 genes. The resulting phylogenetic tree (Figure 1) shows that Snf2 proteins can be divided into six groups: Snf2-like, Swr1-like, Rad54-like, SSO1653-like, SMARCAL1-like, and Rad5/16-like, containing 18 different subgroups (Table S2). The distribution across these subfamilies was uneven, with the DRD1 evolutionary branch being the most expanded, comprising 10 members. It was followed by Rad5/16 with nine members, Ris1 with eight members, Snf2 also with eight members, and ERCC6 with six members. In contrast, other subfamilies such as SMARCAL1, Mi-2, and Iswi each had four members, while the Swr1, Ino80, SHPRH, and Lsh subfamilies each had three members. The subfamilies ALC1, Etl1, ATRX, Rad54, and Mot1 each contained only two Snf2 members, and the Chd1 subfamily included just one Arabidopsis Snf2 member. This systematic identification and classification provide a foundational framework for further studies on the functional roles of the maize Snf2 gene family in development processes and stress responses.

2.2. Analysis of Physicochemical Properties of Maize Snf2-Family Proteins

The physicochemical properties of ZmCHR-family proteins were comprehensively analyzed using the ProtParam tool and the CELLO online platform. Although there are differences among members, these variations in physicochemical properties are relatively minor (Table 1). Analysis of protein structural characteristics revealed that the number of amino acids in these proteins ranged from 411 to 3828. The isoelectric point (pI) spanned a range from 5.07 (ZmCHR20) to 9.66 (ZmCHR24). The instability index varied between 40.88 and 64.6. The aliphatic index, which reflects the relative volume occupied by aliphatic side chains (alanine, valine, isoleucine, and leucine), ranged from 66.3 to 97.49, serving as an indicator of protein thermostability. Molecular weights were observed to be between 46.6 kDa and 415 kDa. Negative GRAVY values (from −0.853 to −0.153) indicated that all Snf2 proteins are hydrophilic. Subcellular localization predictions showed that most of the family members are predominantly localized in the nucleus, with only a few found in the cytoplasm and mitochondria. This analysis provides insights into the key structural features of maize Snf2 genes and proteins. Furthermore, the goal of describing protein domain families is to better understand the functions of the proteins.

2.3. Chromosomal Distribution of Maize Snf2 Genes

The chromosomal distribution of the maize Snf2 gene family was mapped across the maize genome, revealing that these genes are distributed on all ten pairs of chromosomes, albeit unevenly (Figure 2). Chromosomes 4 and 6 each host six ZmCHR genes, while chromosomes 2 and 5 each contain four genes. Chromosomes 1, 9, and 10 each have three genes, with the remaining chromosomes containing two genes each. This widespread chromosomal distribution underscores the evolutionary significance and potential functional diversity of the ZmCHR gene family in maize.

2.4. Gene Structure and Protein Domain Analysis of the Maize Snf2 Gene Family Members

Variations in exon–intron arrangement are a significant source of gene family diversity and contribute to plant diversity by causing changes in gene expression and function. Analysis of the exon/intron structures of the 35 Snf2 family genes revealed considerable variation in exon counts, ranging from forty (ZmCHR20) to four (ZmCHR14) (Figure 3A). Further analysis indicated that members of the DRD1 subfamily generally possess fewer exons than members of other subfamilies. Snf2-family genes typically feature long sequences, with ten genes exceeding 20 kb and five genes exceeding 30 kb.
To study the homology among members of the maize Snf2 gene family, 20 conserved motifs were identified (Figure 3B, Table S3). The number of motifs within each ZmCHR varied from six (ZmCHR13) to eighteen (ZmCHR9 and ZmCHR28). Thirty members contained motifs 1 and 7, while thirty-four members, excluding ZmCHR31, contained motif 2, indicating that motif 2 is a relatively conserved pattern, likely involved in multiple cellular processes. Additionally, the composition of motifs is quite conservative within the same subfamily.
The goal of describing protein domain families is to better understand protein functions. To elucidate the roles of Snf2 family proteins, sequences were analyzed using Pfams and the Conserved Domain Database (CDD), identifying 20 conserved domains (Figure 3B). Analyses showed that nearly all members of the DRD1, SMARCAL1, SHPRH, Rad5/16, and Ris1 subfamilies, except ZmCHR12, possess the Helicase_C and SNF2-rel_dom superfamily domains. The Helicase_C and SNF2-rel_dom domains were found to be present in all subfamilies except Mi-2, DRD1, SMARCAL1, SHPRH, Rad5/16, and Ris1. The Chromo domain was exclusively found in the Mi-2 subfamily of maize Snf2 genes.

2.5. Collinearity Analysis of the Maize Snf2 Gene Family

To elucidate the collinearity among Snf2 family genes, intraspecific analyses within maize and interspecific analyses among maize, Arabidopsis, rice, and Sorghum were conducted using MCScanX. This analysis, facilitated by TBtools, assessed the thirty-five identified ZmCHR genes and identified three pairs of duplicated genes: ZmCHR1-ZmCHR31, ZmCHR35-ZmCHR10, and ZmCHR15-ZmCHR18 (Figure 4). Each pair of duplicated genes belongs to the same subgroup, suggesting that these chromosomal segments may have undergone duplication events without complete diversification, potentially leading to functional redundancy. The synonymous and nonsynonymous substitution rates for these duplicated gene pairs were calculated using the Simple Ka/Ks Calculator in TBtools. Results indicated that all six ZmCHR tandem duplicate sequences had a Ka/Ks ratio less than one (Table S4). It is widely accepted that a Ka/Ks ratio less than one indicates negative or purifying selection, a ratio of one indicates neutral selection, and a ratio greater than one suggests positive selection. These findings suggest that the duplicated genes in maize have undergone purifying selection aimed at eliminating deleterious mutations. This purifying selection underscores the evolutionary pressure to maintain the fundamental functions of the Snf2 gene family in maize and reduce redundancy, further revealing the genetic stability and functional integrity of these genes.
Additionally, the orthologous relationships of Snf2 genes between Zea mays and two other Poaceae species (Oryza sativa and Sorghum bicolor) were investigated to explore the origins and evolutionary relationships of the Snf2 genes (Figure 5). There are 31 orthologous gene pairs between Z. mays and S. bicolor, and 26 between Z. mays and O. sativa. No orthologous gene pairs were found with Arabidopsis thaliana, and the extensive homology between Z. mays, O. sativa, and S. bicolor supports their close evolutionary relationships within the Poaceae family.

2.6. Analysis of Cis-Regulatory Elements in the Promoters of Maize Snf2 Genes

To explore the potential regulatory functions of Snf2 genes under abiotic and biotic stresses, the 2 Kb promoter sequences of Snf2 genes were analyzed to detect cis-regulatory elements (CREs). Eight plant-hormone-responsive cis-regulatory elements were identified in the promoter regions of Snf2 genes, including ABRE elements, TCA element, CGTCA motif, TGACG motif, TGA element, AuxRR-core, GARE-motif, and P-box (Figure 6). ABRE (abscisic acid-responsive) elements were identified in 32 Snf2 genes, while CGTCA (MeJA-responsive) and TGACG (MeJA-responsive) motifs were found in 30 Snf2 genes each. Additionally, 33, 31, 22, 16, 19, and 33 ZmCHR genes contained G-Box (abscisic acid- and light-responsive), ARE (antioxidant response element), DRE (dehydration-responsive element), LTR (low-temperature responsive), MBS (drought-inducible element), and CAT-box, respectively. The visualization of cis-regulatory elements in plant promoters was facilitated using TBtools, revealing that the cis-regulatory elements in 35 ZmCHRs are related to light response, abscisic acid response, methyl jasmonate response, anaerobic induction, drought induction via MYB-binding sites, and low-temperature response (Figure 6). These results suggest that the maize Snf2 gene family may be involved in regulating responses to abiotic and biotic stresses, as well as hormone signaling pathways.

2.7. Expression Analysis of Maize Snf2 Genes in Various Tissues

In examining the expression patterns of Snf2 genes across different tissues and developmental stages, we explored the dynamic changes within the maize Snf2 gene family. Snf2 genes exhibited tissue-specific expression profiles. Members of the maize Snf2 family were predominantly expressed during the reproductive growth stages, specifically in immature tassels and cobs. Minimal to no expression was observed during early developmental stages (coleoptile 6 days after sowing GH and primary root 6 days after sowing GH), vegetative growth stages (stem apical meristem, young leaves, and early whole root system), grain development stages (endosperm and whole seed at 24 days), and other critical tissues (anthers, rachises before pollination, and leaves and internodes before and after pollination). The immature rachis represents a crucial stage in the formation of maize reproductive organs, characterized by rapid cell proliferation and differentiation. The Snf2 gene family typically encodes core subunits of chromatin remodeling complexes that play key roles in regulating gene expression, DNA repair, and chromatin conformation changes. Therefore, the elevated expression of Snf2 genes during these periods of high proliferation and differentiation is logical. Notably, ZmCHR27 was highly expressed in early primary roots, ZmCHR5 and ZmCHR8 showed higher expression levels in mid-vegetative-growth leaves, ZmCHR21 was highly expressed in anthers, and ZmCHR35 exhibited increased expression in internodes post-pollination. These findings suggest their potential significance in plant growth and development. On the other hand, the high expression of ZmCHR33 in embryos at 24 days suggests their involvement in seed germination (Figure 7, Table S5).

2.8. Expression Analysis of Maize Snf2 Genes under Various Stress Conditions

To investigate the expression dynamics of maize Snf2 genes under abiotic stress conditions, RNA-seq data for the inbred line B73 were retrieved from the GDB database (https://maizegdb.org/ (accessed on 6 June 2024)). We analyzed the expression patterns of maize Snf2 family members in aerial tissues of seedlings under cold, heat, high-salinity, and UV stress; in leaves under drought, salinity, and combined stress; and in primary roots of seedlings under varying degrees of drought stress, as documented by Makarevitch [26] et al. (2015), Forestan [27] et al. (2016), and Opitz [28] et al. (2014). Our findings indicated significant downregulation of ZmCHR19, ZmCHR5, ZmCHR14, ZmCHR3, ZmCHR7, and ZmCHR16 under cold stress, while ZmCHR8, ZmCHR28, ZmCHR2, and ZmCHR24 were notably upregulated. All studied genes were downregulated under high-temperature stress. Under salt stress, ZmCHR19, ZmCHR5, ZmCHR14, ZmCHR3, and ZmCHR7 showed significant downregulation, whereas ZmCHR33, ZmCHR27, ZmCHR4, ZmCHR16, ZmCHR28, ZmCHR6, and ZmCHR31 were markedly upregulated. Under UV stress, ZmCHR5 and ZmCHR14 were significantly downregulated, while other genes (except ZmCHR19, ZmCHR33, ZmCHR8, ZmCHR2, ZmCHR6) showed significant upregulation. T0 represents ten days after stress treatment, and T7 represents seven days of normal conditions following the ten-day treatment, marking recovery from stress. Ten days after combined stress, ZmCHR24 to ZmCHR28 showed marked upregulation, while ZmCHR19 to ZmCHR22 were notably downregulated. Under mild drought conditions, ZmCHR23 to ZmCHR32 were significantly upregulated, whereas ZmCHR35, ZmCHR7, and ZmCHR12 were downregulated. Under severe drought, ZmCHR25, ZmCHR24, ZmCHR27, and ZmCHR33 were markedly upregulated, and ZmCHR35, ZmCHR7, and ZmCHR12 were significantly downregulated (Figure 8, Tables S6–S8).

2.9. Analysis of Maize Snf2 Gene Expression under Drought Stress and qPCR Validation

Based on the gene expression profiles of maize Snf2 genes under drought stress, four genes were randomly selected among the genes upregulated and downregulated at 24 h, respectively, and analyzed by qRT-PCR at 0, 6, 12, and 24 h of drought stress. The results showed that ZmCHR1 and ZmCHR23 exhibited an upward trend at 12 h, followed by a decrease at 24 h, although they remained upregulated compared to the control group. ZmCHR25 displayed a rapid increase in expression levels after 24 h of drought stress. ZmCHR31 rose sharply at 12 h after stress and remained constant at 24 h. In contrast, ZmCHR35 and ZmCHR12 exhibited a downward trend at 12 h. ZmCHR24 and ZmCHR33 showed significant decreases starting at 24 h (Figure 9). The expression patterns of the selected genes were consistent with transcriptome data. Overall, certain Snf2 genes in maize respond to drought stress, and the observed divergent trends may indicate varying degrees of drought tolerance in maize.

3. Discussion

3.1. Functional Domains and Evolutionary Dynamics of the Maize Snf2 Gene Family

Snf2 proteins, acting as core ‘motors’, have been characterized at the whole-genome level in various crops including Arabidopsis, wheat, rice, and sorghum. In this study, we systematically identified and characterized 35 Snf2 family genes on 10 chromosomes in maize. In previous studies, 38 genes encoding Snf2 proteins were identified in barley and sorghum [25,29], 41 in Arabidopsis [22], 45 in tomato [24], and 40 in rice [23], and the number of genes encoding Snf2 proteins in maize was less than the numbers in the above crops, suggesting that the differences in species and genome sizes are not directly related to the number of Snf2 family genes, which lays the foundation for revealing multiple chromatin remodeling mechanisms mediated by different Snf2 protein subfamilies. As catalytic ATPase subunits of chromatin remodeling complexes, Snf2 proteins contain conserved helicase-C and SNF2-N domains [30], enabling their involvement in multiple DNA processes, including replication [31], transcription [32], homologous recombination [33], and repair [34]. These conserved domains mediate the varied chromatin remodeling activities of Snf2 proteins in these fundamental nuclear events. Such functional diversification allows Snf2 proteins to participate in regulating various plant-specific developmental processes and stress response pathways [35].
Phylogenetic analysis of the maize Snf2 protein family classified the thirty-five maize Snf2 members into six major groups and seventeen subfamilies, with each major group containing between two and ten maize members. Interestingly, most subfamilies (Snf2, Iswi, ALC1, EtI1, ATRX, Rad54, ERCC6, Mot1, and SMARCAL1) exhibit a 1:1 orthologous conservation pattern between maize and Arabidopsis, suggesting that these subfamilies are evolutionarily more conserved than others. This was similarly found in wheat, Arabidopsis, and rice [36].

3.2. Function and Diversity: The Role of the Maize Snf2 Gene Family in Chromatin Remodeling

The Snf2 gene family in maize plays a crucial role as a component of the chromatin remodeling complex, which is essential for plant growth and environmental adaptability [37]. Due to the sedentary growth nature of plants, they must adapt in situ to external changes [38]. This unique mode of growth means that plants exhibit distinct epigenetic regulatory mechanisms compared to animals. Historically, there has been a relative lack of systematic research into the classification, subunit composition, and role of chromatin remodeling complexes in regulating transcription in plants. However, with the integration of technologies such as mass spectrometry, protein interaction network analysis, high-throughput sequencing, and protein structure resolution, significant progress has been made in this field [39]. The SWI/SNF complex, first identified through genetic screening in yeast, includes the Swi2/Snf2 protein, which contains evolutionarily conserved Helicase_C and SNF2_N domains. These domains alter chromatin structure through their catalytic ATPase activities, impacting DNA accessibility and associated biological processes. Additionally, maize Snf2 proteins include supplementary domains that demonstrate subfamily-dependent functional diversity. In contrast to all others, members of the DRD1, SMARCAL1, Ris1, Rad5/16 (except ZmCHR12), and SHPRH subfamilies contain a SNF2-related domain superfamily member upstream of the Helicase_C domain, whose ATPase activity regulates the chromatin remodeling process, thereby influencing the accessibility of DNA and its related biological processes.

3.3. Differential Expression Patterns of Maize Snf2 Gene in Various Organs

Spatio-temporal expression patterns of genes during reproductive development may provide clues to elucidate their potential functions. Most of the Snf2 genes in maize are expressed during the reproductive growth stage and are particularly highly expressed in immature tassels and cobs. This suggests that Snf2 family genes may activate cell division and differentiation. Yu et al. found that the core subunit of the SWI/SNF chromatin remodeling complex, ZmCHB101, is involved in vesicle division and development and that transgenic lines of the ZmCHB101 RNA interference (RNAi) construct had impaired spike and rachis development [40]. The Snf2 gene ZmDDM1 has been reported in maize, and its homologs are abundantly expressed in the embryo; once the gene is deleted, abnormal cell proliferation occurs during the early stages of kernel development [41]. There is also evidence for this in other crops. Chen et al. found that the Snf2 gene family has multiple regulatory roles in spike development in barley [25]. Several Snf2 genes in rice have been reported to be specifically expressed in the early development of the pistil [42]. OsCHR721 interacts with OsRPA1a, which is involved in DNA repair and plays a key role in male and female reproductive development [43]. Several previous studies in Arabidopsis have also shown that the Snf2 gene is essential for regulating flower developmental processes such as floral transition [44], floral organ characterization [45], and inflorescence structure [46]. These findings suggest that the Snf2 gene has potential research implications in maize growth and development.

3.4. Role of the Maize Snf2 Gene Family in Response to Abiotic Stress

When plants are subjected to abiotic stresses, they undergo a variety of physiological and biochemical responses that alter the gene expression program to rapidly adapt to adverse conditions. To investigate the expression pattern of the Snf2 gene under abiotic stress, we performed a computer-simulated cis-regulatory element analysis of the putative promoter region of the maize Snf2 gene, and we identified several cis-regulatory elements associated with abiotic stresses such as gibberellic acid and drought. Subsequent analysis of the Snf2 gene in different tissues of maize under different stresses using publicly available gene expression profiling databases revealed a significant upregulation/downregulation trend of the maize Snf2 gene under high temperature, low temperature, salt, UV, and drought stresses. This gives a preliminary validation of the Snf2 gene family in response to abiotic stresses. Finally, we conducted qRT-PCR experiments using eight genes selected in the expression profiles for further validation. The results showed that all eight selected genes responded to drought stress.
The response of the Snf2 gene family to abiotic stresses has been reported in several plant species. In Arabidopsis, the SWI/SNF subunits are involved in many biotic and abiotic stress responses; the PKL mutant pkl is highly sensitive to cold stress, resulting in reduced cotyledon greening and reduced primary root elongation at high salinity [47]; and BRM negatively regulates drought tolerance in Arabidopsis [48]. Both AtLFR in Arabidopsis and its homologous gene GmLFR1 in soybean are negative regulators of drought tolerance [49]. A previous report showed that overexpression of the SWI/SNF complex subunit OsSWI3C in rice resulted in reduced drought tolerance [50], and OsCHR710 was upregulated under drought stress [42]. The involvement of the SWI2/SNF2 genes in UV-B tolerance and acclimatization was mentioned in an analysis showing that chromatin remodeling proteins are responsive to UV-B [51]. These findings emphasize the close link between at least some of the Snf2 genes and plant stress responses, as well as their possible function in drought tolerance, highlighting their key role in addressing related environmental challenges.

4. Materials and Methods

4.1. Identification and Physicochemical Property Prediction of the Maize Snf2 Gene

The maize genome and GFF annotation files were downloaded from the Ensembl Plants database (https://plants.ensembl.org/index.html, accessed on 1 June 2024). The amino acid sequence of the Arabidopsis Snf2 gene was retrieved from The Arabidopsis Information Resource (TAIR) (https://www.arabidopsis.org/, accessed on 1 June 2024). Local BLAST searches were conducted using the ‘Blast Compare Two Seqs’ module within TBtools (v2.097), setting the E-value threshold to 10−5. Furthermore, characteristic structural domains of the maize Snf2 gene family (PF00176 and PF00271) were identified using the Pfam database (http://pfam.xfam.org/, accessed on 4 June 2024) through a Simple HMM Search. The intersection of sequence IDs from two screening steps was taken to obtain the longest amino acid sequence variant within different transcripts using the ‘Fasta Get Representative’ module in TBtools (v2.097). Finally, sequences with incomplete conserved domains were excluded using the Conserved Domain Database (CDD) search tool provided by the National Center for Biotechnology Information (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi, accessed on 4 June 2024).
The physicochemical properties of the predicted maize Snf2 protein were calculated using the ProtParam tool (https://web.expasy.org/protparam/, accessed on 4 June 2024). These properties included amino acid sequence length, theoretical isoelectric point (pI), molecular weight (Mw), instability index, aliphatic index, and grand average of hydropathicity (GRAVY) index. Subcellular localization of the maize Snf2 gene was predicted using the CELLO (v.2.5) online tool (http://cello.life.nctu.edu.tw/, accessed on 4 June 2024).

4.2. Chromosomal Localization of the Maize Snf2 Gene

Chromosomal position information for the maize Snf2 gene was extracted from the maize genome annotation file. Visualization of this chromosomal localization was performed using the ‘Gene Location Visualize from GTF/GFF’ module in TBtools (v2.097).

4.3. Phylogenetic, Gene Structure, and Conserved Motif Analysis of the Maize Snf2 Gene

Amino acid sequences of the Snf2 gene from maize and Arabidopsis were selected for phylogenetic analysis. Multiple sequence alignments were conducted using MUSCLE, followed by trimming with trimAL. A phylogenetic tree was constructed using the maximum likelihood (ML) method in the IQ-TREE software (1.6.12), with a bootstrap value set to 5000. The resulting phylogenetic tree was visualized and embellished using iTOL (https://itol.embl.de/, accessed on 5 June 2024).
The motif architecture of the maize Snf2 gene family members was analyzed using the online MEME suite (http://www.OMIcsclass.com/article/67, accessed on 5 June 2024). The maximum number of motifs was set to 20, with minimum and maximum motif widths set at 6 and 50, respectively. The gene structure of maize Snf2, including coding sequences (CDSs) and untranslated regions (UTRs), was extracted from the maize genome annotation file. Visualization of these features was performed using TBtools (v2.097).

4.4. Collinearity Analysis and Gene Duplication of the Maize Snf2 Gene

Homologous Snf2 genes within maize were identified using the BLASTP program, with an E-value threshold set to <10−5. Collinearity among maize Snf2 genes was analyzed using the default parameters in the MCScanX program. Homologous blocks and duplicated Snf2 gene pairs were visualized using TBtools. Genome sequences and annotation files for Arabidopsis and rice were downloaded from the Phytozome v13 website (https://phytozome-next.jgi.doe.gov/, accessed on 5 June 2024), and those for sorghum were obtained from the NCBI website (https://www.ncbi.nlm.nih.gov/, accessed on 5 June 2024). Collinearity between maize and Arabidopsis, rice, and sorghum was analyzed using the MCScanX program.

4.5. Analysis of Cis-Acting Elements in the Promoter Regions of the Maize Snf2 Gene

The 2000 bp region upstream of the start codon (ATG) of each maize Snf2 gene was considered as the promoter sequence. These promoter sequences were extracted using TBtools (v2.097) and analyzed using the PlantCARE database (https://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 6 June 2024) to predict cis-regulatory elements. Visualization and summary of cis-elements related to stress response, plant growth and development, plant hormone responsiveness, and light responsiveness were performed using the Basic Biosequence View and HeatMap modules in TBtools (v2.097).

4.6. Expression Analysis of the Snf2 Gene Family in Maize

Transcriptomic data for the maize Snf2 gene in various tissues were analyzed using the qTeller tool (https://qTeller.MaizeGDB.org, accessed on 6 June 2024) MaizeGDB, which included tissues such as the immature tassel and rachis, embryo, endosperm, unpollinated internodes, unpollinated leaves, pollinated internodes, pollinated leaves, stem apical meristem, complete root system, and whole seed (https://qTeller.MaizeGDB.org, accessed on 6 June 2024). These analyses were conducted using RNA-Seq data from Stelpflug et al. [52]. Additionally, the transcriptional expression levels of the maize Snf2 gene under various abiotic stress conditions were examined using RNA-Seq data from studies by Makarevitch [26], Forestan [27], and Opitz [28], among others.

4.7. Plant Material and Stress Treatment Procedures

Maize inbred line B73 was used as the plant material. Seeds were sterilized with 75% alcohol and 2% sodium hypochlorite before being planted in autoclaved nutrient soil. The plants were subsequently grown in a controlled greenhouse environment, maintained at 25 °C with a photoperiod of 8 h of light followed by 16 h of darkness. At the three-leaf stage, seedlings with uniform growth were selected for drought stress experiments. To simulate drought conditions, these seedlings were treated with 15% PEG6000. The second fully expanded leaves from three seedlings were collected at 0, 6, 12, and 24 h post-treatment and immediately frozen in liquid nitrogen at −80 °C for subsequent RNA extraction.

4.8. RNA Extraction and Quantitative Real-Time PCR Analysis

A 150–300 mg mass of leaves was collected at each stress time point, then added to liquid nitrogen and quickly ground into powder. The resulting powder was divided into triplicates, and then RNA was extracted using an RNA extraction kit (CW Bio, Taizhou, China). The quality of the extracted RNA was assessed through gel electrophoresis and quantified using a NanoDrop ND-2000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA). RNA samples that met the quality standards were pooled and reverse-transcribed into cDNA, with concentrations subsequently measured. The cDNA from each sample was then uniformly diluted to 200 ng/μL for use as templates in real-time quantitative PCR (qRT-PCR) experiments.
To explore the expression patterns of maize Snf2 genes under abiotic stress, eight ZmCHR genes were randomly selected for qRT-PCR analysis. Specific primers for these genes were designed using NCBI’s Primer-BLAST tool (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 6 June 2024), and qRT-PCR was performed using SYBR™ Green PCR Master Mix (Thermo Fisher Scientific, Waltham, CA, USA). Each sample was analyzed in triplicate, with the 0 h stress time point serving as the control. The results were assessed using the relative quantification (RQ) method. The reaction system and conditions for RT-qPCR were based on the method outlined by Mo et al. [53]. All primer sequences utilized in this study are provided in Supplementary Table S9.

5. Conclusions

This study identified and comprehensively analyzed 35 ZmCHR genes in the maize genome. The structure and chromosomal locations of these genes were detailed, a phylogenetic tree was constructed, and the properties of their encoded proteins were analyzed. Cis-regulatory element and transcriptional expression analyses suggest that the SWI/SNF complex may play a significant role in maize’s response to environmental stress. This research provides a critical foundation for understanding the composition of the SWI/SNF complex in maize and paves the way for further investigation into its role in abiotic stress responses.

Supplementary Materials

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

Author Contributions

Y.Y. and C.Z. designed the experiments and methods. J.H., Q.W. and B.Q. performed the experiments and data analysis. Q.L., Z.W., Y.L., Z.C. and W.W. assisted with the experiments. Q.W., W.W. and Z.C. participated in the discussion for the experimental direction. Q.W. and Y.Y. supervised and guided the experiments. J.H. and Q.W. drafted the manuscript. Y.Y., Q.W. and C.Z. proofread the manuscript. C.Z. and Y.Y. reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Jilin Science and Technology Development Program Project (20240404056ZP), the China Postdoctoral Science Foundation (2024M751102), the Jilin Agricultural Science and Technology Innovation Project (CXGC2024ZD001), and the Jilin Scientific and Technological Development Program (20220101330JC).

Informed Consent Statement

All co-authors involved in the paper consent to publish this article in the International Journal of Molecular Sciences.

Data Availability Statement

All relevant data are provided as figures or tables within the paper.

Acknowledgments

We thank our contributors for their dedication and compliance through the many stages of this research, as well as the editors and anonymous reviewers whose comments helped to greatly improve this paper.

Conflicts of Interest

The authors declare that they have no conflicts of interest or competing interests.

References

  1. Ntiamoah, E.B.; Li, D.; Appiah-Otoo, I.; Twumasi, M.A.; Yeboah, E.N. Towards a sustainable food production: Modelling the impacts of climate change on maize and soybean production in Ghana. Environ. Sci. Pollut. Res. 2022, 29, 72777–72796. [Google Scholar] [CrossRef] [PubMed]
  2. Shiferaw, B.; Prasanna, B.M.; Hellin, J.; Bänziger, M. Crops that feed the world 6. Past successes and future challenges to the role played by maize in global food security. Food Secur. 2011, 3, 307–327. [Google Scholar] [CrossRef]
  3. Zhang, H.; Zhu, J.; Gong, Z.; Zhu, J.K. Abiotic stress responses in plants. Nat. Rev. Genet. 2022, 23, 104–119. [Google Scholar] [CrossRef] [PubMed]
  4. Zhao, C.; Liu, B.; Piao, S.; Wang, X.; Lobell, D.B.; Huang, Y.; Huang, M.T.; Yao, Y.T.; Bassu, S.; Ciais, P.; et al. Temperature increase reduces global yields of major crops in four independent estimates. Proc. Natl. Acad. Sci. USA 2017, 114, 9326–9331. [Google Scholar] [CrossRef] [PubMed]
  5. van Bakel, H.; Hughes, T. (Eds.) Interactions of Transcription Factors with Chromatin; Springer: Dordrecht, The Netherlands, 2011; Volume 52. [Google Scholar]
  6. Singh, A.; Modak, S.B.; Chaturvedi, M.M.; Purohit, J.S. SWI/SNF chromatin remodelers: Structural, functional and mechanistic implications. Cell Biochem. Biophys. 2023, 81, 167–187. [Google Scholar] [CrossRef]
  7. Wang, J.; Sun, Z.; Liu, H.; Yue, L.; Wang, F.; Liu, S.; Su, B.; Liu, B.; Kong, F.; Fang, C. Genome-wide identification and characterization of the soybean Snf2 gene family and expression response to rhizobia. Int. J. Mol. Sci. 2023, 24, 7250. [Google Scholar] [CrossRef]
  8. Ryan, D.P.; Owen-Hughes, T. Snf2-family proteins: Chromatin remodellers for any occasion. Curr. Opin. Chem. Biol. 2011, 15, 649–656. [Google Scholar] [CrossRef]
  9. Cairns, B.R.; Kim, Y.J.; Sayre, M.H.; Laurent, B.C.; Kornberg, R.D. A multisubunit complex containing the SWI1/ADR6, SWI2/SNF2, SWI3, SNF5, and SNF6 gene products isolated from yeast. Proc. Natl. Acad. Sci. USA 1994, 91, 1950–1954. [Google Scholar] [CrossRef]
  10. Ojolo, S.P.; Cao, S.; Priyadarshani, S.; Li, W.; Yan, M.; Aslam, M.; Zhao, H.; Qin, Y. Regulation of plant growth and development: A review from a chromatin remodeling perspective. Front. Plant Sci. 2018, 9, 1232. [Google Scholar] [CrossRef]
  11. Arnholdt-Schmitt, B. Stress-induced cell reprogramming. A role for global genome regulation. Plant Physiol. 2004, 136, 2579–2586. [Google Scholar] [CrossRef]
  12. 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]
  13. Song, Z.T.; Liu, J.X.; Han, J.J. Chromatin remodeling factors regulate environmental stress responses in plants. J. Integr. Plant Biol. 2021, 63, 438–450. [Google Scholar] [CrossRef]
  14. Sarnowska, E.; Gratkowska, D.M.; Sacharowski, S.P.; Cwiek, P.; Tohge, T.; Fernie, A.R.; Siedlecki, J.A.; Koncz, C.; Sarnowski, T.J. The role of SWI/SNF chromatin remodeling complexes in hormone crosstalk. Trends Plant Sci. 2016, 21, 594–608. [Google Scholar] [CrossRef]
  15. Walley, J.W.; Rowe, H.C.; Xiao, Y.; Chehab, E.W.; Kliebenstein, D.J.; Wagner, D.; Dehesh, K. The chromatin remodeler SPLAYED regulates specific stress signaling pathways. PLoS Pathog. 2008, 4, e1000237. [Google Scholar] [CrossRef]
  16. Guo, T.; Wang, D.; Fang, J.; Zhao, J.; Yuan, S.; Xiao, L.; Li, X. Mutations in the rice OsCHR4 gene, encoding a CHD3 family chromatin remodeler, induce narrow and rolled leaves with increased cuticular wax. Int. J. Mol. Sci. 2019, 20, 2567. [Google Scholar] [CrossRef]
  17. Mlynárová, L.; Nap, J.P.; Bisseling, T. The SWI/SNF chromatin-remodeling gene AtCHR12 mediates temporary growth arrest in Arabidopsis thaliana upon perceiving environmental stress. Plant J. 2007, 51, 874–885. [Google Scholar] [CrossRef]
  18. Yu, X.; Meng, X.; Liu, Y.; Li, N.; Zhang, A.; Wang, T.-J.; Jiang, L.; Pang, J.; Zhao, X.; Qi, X.; et al. The chromatin remodeler ZmCHB101 impacts expression of osmotic stress-responsive genes in maize. Plant Mol. Biol. 2018, 97, 451–465. [Google Scholar] [CrossRef]
  19. Wang, Y.; Wang, D.; Gan, T.; Liu, L.; Long, W.; Wang, Y.; Niu, M.; Li, X.; Zheng, M.; Jiang, L.; et al. CRL6, a member of the CHD protein family, is required for crown root development in rice. Plant Physiol. Biochem. 2016, 105, 185–194. [Google Scholar] [CrossRef]
  20. Campos-Rivero, G.; Osorio-Montalvo, P.; Sánchez-Borges, R.; Us-Camas, R.; Duarte-Aké, F.; De-La-Peña, C. Plant hormone signaling in flowering: An epigenetic point of view. J. Plant Physiol. 2017, 214, 16–27. [Google Scholar] [CrossRef]
  21. Flaus, A.; Martin, D.M.A.; Barton, G.J.; Owen-Hughes, T. Identification of multiple distinct Snf2 subfamilies with conserved structural motifs. Nucleic Acids Res. 2006, 34, 2887–2905. [Google Scholar] [CrossRef]
  22. Knizewski, L.; Ginalski, K.; Jerzmanowski, A. Snf2 proteins in plants: Gene silencing and beyond. Trends Plant Sci. 2008, 13, 557–565. [Google Scholar] [CrossRef] [PubMed]
  23. Hu, Y.; Zhu, N.; Wang, X.; Yi, Q.; Zhu, D.; Lai, Y.; Zhao, Y. Analysis of rice Snf2 family proteins and their potential roles in epigenetic regulation. Plant Physiol. Bioch. 2013, 70, 33–42. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, D.; Gao, S.; Yang, P.; Yang, J.; Yang, S.; Wu, K. Identification and expression analysis of Snf2 family proteins in tomato (Solanum lycopersicum). Int. J. Genom. 2019, 2019, 5080935. [Google Scholar] [CrossRef] [PubMed]
  25. Chen, G.; Mishina, K.; Zhu, H.; Kikuchi, S.; Sassa, H.; Oono, Y.; Komatsuda, T. Genome-wide analysis of Snf2 gene family reveals potential role in regulation of spike development in barley. Int. J. Mol. Sci. 2022, 24, 457. [Google Scholar] [CrossRef] [PubMed]
  26. Makarevitch, I.; Waters, A.J.; West, P.T.; Stitzer, M.; Hirsch, C.N.; Ross-Ibarra, J.; Springer, N.M. Transposable elements contribute to activation of maize genes in response to abiotic stress. PLoS Genet. 2015, 11, e1004915. [Google Scholar]
  27. Forestan, C.; Aiese Cigliano, R.; Farinati, S.; Lunardon, A.; Sanseverino, W.; Varotto, S. Stress-induced and epigenetic-mediated maize transcriptome regulation study by means of transcriptome reannotation and differential expression analysis. Sci. Rep. 2016, 6, 30446. [Google Scholar] [CrossRef]
  28. Opitz, N.; Paschold, A.; Marcon, C.; Malik, W.A.; Lanz, C.; Piepho, H.-P.; Hochholdinger, F. Transcriptomic complexity in young maize primary roots in response to low water potentials. BMC Genom. 2014, 15, 741. [Google Scholar] [CrossRef]
  29. Hu, Y.; Chen, X.; Zhou, C.; He, Z.; Shen, X. Genome-wide identification of chromatin regulators in Sorghum bicolor. 3 Biotech 2022, 12, 117. [Google Scholar] [CrossRef]
  30. Kwon, C.S.; Wagner, D. Unwinding chromatin for development and growth: A few genes at a time. Trends Genet. 2007, 23, 403–412. [Google Scholar] [CrossRef]
  31. Vincent, J.A.; Kwong, T.J.; Tsukiyama, T. ATP-dependent chromatin remodeling shapes the DNA replication landscape. Nat. Struct. Mol. Biol. 2008, 15, 477–484. [Google Scholar] [CrossRef]
  32. Raab, J.R.; Resnick, S.; Magnuson, T. Genome-wide transcriptional regulation mediated by biochemically distinct SWI/SNF complexes. PLoS Genet. 2015, 11, e1005748. [Google Scholar] [CrossRef] [PubMed]
  33. Ceballos, S.J.; Heyer, W.-D. Functions of the Snf2/Swi2 family Rad54 motor protein in homologous recombination. BBA-Gene Regul. Mech. 2011, 1809, 509–523. [Google Scholar] [CrossRef]
  34. Lans, H.; Marteijn, J.A.; Vermeulen, W. ATP-dependent chromatin remodeling in the DNA-damage response. Epigenet. Chromatin 2012, 5, 4. [Google Scholar] [CrossRef] [PubMed]
  35. Peng, R.; Wang, J.; Yang, S. Functions of chromatin remodeling factors during plant development. J. Trop. Subtrop. Bot. 2019, 27, 565–579. [Google Scholar]
  36. Han, S.-K.; Wagner, D. Role of chromatin in water stress responses in plants. J. Exp. Bot. 2014, 65, 2785–2799. [Google Scholar] [CrossRef]
  37. Mohammadi, M.A.; Han, X.; Zhang, Z.; Xi, Y.; Boorboori, M.; Wang-Pruski, G. Phosphite application alleviates pythophthora infestans by modulation of photosynthetic and physio-biochemical metabolites in potato leaves. Pathogens 2020, 9, 170. [Google Scholar] [CrossRef]
  38. Fahad, M.; Liu, C.; Shen, Y.; Sajid, M.; Wu, L. enomic exploration and functional insights into the diversification of the snf2 gene family in subgenomes of Triticum aestivum. Crop Des. 2024, 3, 100050. [Google Scholar]
  39. Dai, X.; Shen, L. Advances and trends in omics technology development. Front. Med. 2022, 9, 911861. [Google Scholar] [CrossRef]
  40. Yu, X.; Jiang, L.; Wu, R.; Meng, X.; Zhang, A.; Li, N.; Xia, Q.; Qi, X.; Pang, J.; Xu, Z.-Y.; et al. The core subunit of a chromatin-remodeling complex, ZmCHB101, plays essential roles in maize growth and development. Sci. Rep. 2016, 6, 38504. [Google Scholar] [CrossRef]
  41. Long, J.C.; Xia, A.A.; Liu, J.H.; Jing, J.L.; Wang, Y.Z.; Qi, C.Y.; He, Y. Decrease in DNA methylation 1 (DDM1) is required for the formation of mCHH islands in maize. J. Integr. Plant Biol. 2019, 61, 749–764. [Google Scholar] [CrossRef]
  42. Guo, M.; Zhao, H.; He, Z.; Zhang, W.; She, Z.; Mohammadi, M.A.; Shi, C.; Yan, M.; Tian, D.; Qin, Y. Comparative expression profiling of Snf2 family genes during reproductive development and stress responses in rice. Front. Plant Sci. 2022, 13, 910663. [Google Scholar] [CrossRef] [PubMed]
  43. Zhang, Y.; Chen, Q.; Zhu, G.; Zhang, F.; Fang, X.; Ren, H.; Jiang, J.; Yang, F.; Zhang, D.; Chen, F. CHR721, interacting with OsRPA1a, is essential for both male and female reproductive development in rice. Plant Mol. Biol. 2020, 103, 473–487. [Google Scholar] [CrossRef] [PubMed]
  44. Jing, Y.; Guo, Q.; Lin, R. The chromatin-remodeling factor PICKLE antagonizes polycomb repression of FT to promote flowering. Plant Physiol. 2019, 181, 656–668. [Google Scholar] [CrossRef] [PubMed]
  45. Li, C.; Chen, C.; Gao, L.; Yang, S.; Nguyen, V.; Shi, X.; Siminovitch, K.; Kohalmi, S.E.; Huang, S.; Wu, K.; et al. The Arabidopsis SWI2/SNF2 chromatin remodeler BRAHMA regulates polycomb function during vegetative development and directly activates the flowering repressor gene SVP. PLoS Genet. 2015, 11, e1004944. [Google Scholar] [CrossRef]
  46. Zhao, M.; Yang, S.; Chen, C.Y.; Li, C.; Shan, W.; Lu, W.; Cui, Y.; Liu, X.; Wu, K. Arabidopsis BREVIPEDICELLUS interacts with the SWI2/SNF2 chromatin remodeling ATPase BRAHMA to regulate KNAT2 and KNAT6 expression in control of inflorescence architecture. PLoS Genet. 2015, 11, e1005125. [Google Scholar] [CrossRef]
  47. Yang, R.; Hong, Y.; Ren, Z.; Tang, K.; Zhang, H.; Zhu, J.-K.; Zhao, C. A role for PICKLE in the regulation of cold and salt stress tolerance in Arabidopsis. Front. Plant Sci. 2019, 10, 900. [Google Scholar] [CrossRef]
  48. Han, S.K.; Sang, Y.; Rodrigues, A.; Wu, M.F.; Rodriguez, P.L.; Wagner, D. The SWI2/SNF2 chromatin remodeling ATPase BRAHMA represses abscisic acid responses in the absence of the stress stimulus in Arabidopsis. Plant Cell 2012, 24, 4892–4906. [Google Scholar] [CrossRef]
  49. Chen, Q.; Shi, X.; Ai, L.; Tian, X.; Zhang, H.; Tian, J.; Wang, Q.; Zhang, M.; Cui, S.; Yang, C.; et al. Genome-wide identification of genes encoding SWI/SNF components in soybean and the functional characterization of GmLFR1 in drought-stressed plants. Front. Plant Sci. 2023, 14, 1176376. [Google Scholar] [CrossRef]
  50. Yang, J.; Chang, Y.; Qin, Y.; Chen, D.; Zhu, T.; Peng, K.; Wang, H.; Tang, N.; Li, X.; Wang, Y.; et al. A lamin-like protein OsNMCP1 regulates drought resistance and root growth through chromatin accessibility modulation by interacting with a chromatin remodeller OsSWI3C in rice. New Phytol. 2020, 227, 65–83. [Google Scholar] [CrossRef]
  51. Casati, P.; Stapleton, A.E.; Blum, J.E.; Walbot, V. Genome-wide analysis of high-altitude maize and gene knockdown stocks implicates chromatin remodeling proteins in response to UV-B. Plant J. 2006, 46, 613–627. [Google Scholar] [CrossRef]
  52. Stelpflug, S.C.; Sekhon, R.S.; Vaillancourt, B.; Hirsch, C.N.; Buell, C.R.; de Leon, N.; Kaeppler, S.M. An Expanded Maize Gene Expression Atlas based on RNA Sequencing and its Use to Explore Root Development. Plant Genome 2016, 9. [Google Scholar] [CrossRef] [PubMed]
  53. Mo, F.; Zhang, N.; Qiu, Y.; Meng, L.; Cheng, M.; Liu, J.; Yao, L.; Lv, R.; Liu, Y.; Zhang, Y.; et al. Molecular Characterization, Gene Evolution and Expression Analysis of the F-Box Gene Family in Tomato (Solanum lycopersicum). Genes 2021, 12, 417. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Phylogenetic tree of maize and Arabidopsis, in which ZmCHR genes are categorized into four main classes: Snf2-like, Swr1-like, Rad54-like, SSO1653-like, SMARCAL1-like, and Rad5/16-like. Within this phylogenetic tree, the AtCHR gene of Arabidopsis is specifically marked with a square and the ZmCHR gene of maize is specifically marked with a circle for easy identification.
Figure 1. Phylogenetic tree of maize and Arabidopsis, in which ZmCHR genes are categorized into four main classes: Snf2-like, Swr1-like, Rad54-like, SSO1653-like, SMARCAL1-like, and Rad5/16-like. Within this phylogenetic tree, the AtCHR gene of Arabidopsis is specifically marked with a square and the ZmCHR gene of maize is specifically marked with a circle for easy identification.
Ijms 25 09686 g001
Figure 2. Chromosomal localization of the maize ZmCHR genes. Chromosome numbers are positioned on the left, with the scale on the left measured in megabases (Mb). An estimated genetic interval of 350 kb was set to accurately depict gene distribution. The color gradient from red to blue across the chromosomes represents gene density from high to low, where dark signifies areas of higher gene density and light indicates areas of lower gene density. Unmarked regions on the chromosomes represent genetic areas lacking gene distribution data.
Figure 2. Chromosomal localization of the maize ZmCHR genes. Chromosome numbers are positioned on the left, with the scale on the left measured in megabases (Mb). An estimated genetic interval of 350 kb was set to accurately depict gene distribution. The color gradient from red to blue across the chromosomes represents gene density from high to low, where dark signifies areas of higher gene density and light indicates areas of lower gene density. Unmarked regions on the chromosomes represent genetic areas lacking gene distribution data.
Ijms 25 09686 g002
Figure 3. The structure, conserved domains, and motif patterns of maize ZmCHR genes. In part (A), the structure of ZmCHR genes is shown, where green boxes represent the untranslated regions (UTRs) at the 5′ and 3′ ends; yellow boxes indicate coding sequences (CDSs); and black lines denote introns, with numbers (0, 1, 2) indicating the phase of introns. Different conserved domains are marked with boxes of various colors. Part (B) depicts the motif patterns of ZmCHR members; Part C describes twenty conserved structural domains, with different colors corresponding to different structural domains. Additionally, the lengths of related gene structures and motif components can be estimated using the corresponding scale at the bottom of the panel.
Figure 3. The structure, conserved domains, and motif patterns of maize ZmCHR genes. In part (A), the structure of ZmCHR genes is shown, where green boxes represent the untranslated regions (UTRs) at the 5′ and 3′ ends; yellow boxes indicate coding sequences (CDSs); and black lines denote introns, with numbers (0, 1, 2) indicating the phase of introns. Different conserved domains are marked with boxes of various colors. Part (B) depicts the motif patterns of ZmCHR members; Part C describes twenty conserved structural domains, with different colors corresponding to different structural domains. Additionally, the lengths of related gene structures and motif components can be estimated using the corresponding scale at the bottom of the panel.
Ijms 25 09686 g003
Figure 4. Collinearity among homologous maize ZmCHR genes. In the diagram, collinear blocks across the whole genome are depicted with a gray background, and red curves are used to connect pairs of duplicated ZmCHR genes.
Figure 4. Collinearity among homologous maize ZmCHR genes. In the diagram, collinear blocks across the whole genome are depicted with a gray background, and red curves are used to connect pairs of duplicated ZmCHR genes.
Ijms 25 09686 g004
Figure 5. The collinearity diagram of Snf2 genes between Z. mays, A. thaliana, O. sativa, and S. bicolor. The gray line in the background represents the collinear gene pairs between the two species, and the Snf2 collinear genes are connected by red lines.
Figure 5. The collinearity diagram of Snf2 genes between Z. mays, A. thaliana, O. sativa, and S. bicolor. The gray line in the background represents the collinear gene pairs between the two species, and the Snf2 collinear genes are connected by red lines.
Ijms 25 09686 g005
Figure 6. The ZmCHR gene family promoter sequences were used to predict cis-acting elements; the predicted cis-acting elements were organized and visualized.
Figure 6. The ZmCHR gene family promoter sequences were used to predict cis-acting elements; the predicted cis-acting elements were organized and visualized.
Ijms 25 09686 g006
Figure 7. The transcriptomic expression levels of the ZmCHR gene family across various tissues, including the embryos, endosperm, immature tassel and cob, internodes and leaves before and after pollination, anthers, silks, complete root system, and whole seeds. Note: DAP, days after pollination; GH, greenhouse; R1, reproductive 1; SAM, shoot apical meristem; 7 d, 7 days after sowing.
Figure 7. The transcriptomic expression levels of the ZmCHR gene family across various tissues, including the embryos, endosperm, immature tassel and cob, internodes and leaves before and after pollination, anthers, silks, complete root system, and whole seeds. Note: DAP, days after pollination; GH, greenhouse; R1, reproductive 1; SAM, shoot apical meristem; 7 d, 7 days after sowing.
Ijms 25 09686 g007
Figure 8. The expression levels of maize ZmCHR gene family members under abiotic stresses are depicted in (A). (A) depicts the expression levels of ZmCHR in the whole above-ground tissues of maize seedlings under cold, heat, salt, and ultraviolet stresses; (B) depicts the expression levels of ZmCHR in maize leaves under different drought and salt stresses; and (C) depicts the expression levels of ZmCHR in the primary roots of maize seedlings under different water deficit stresses. Note: for cold stress, seedlings were incubated at 5 °C for 16 h; for heat stress, seedlings were incubated at 50 °C for 4 h; for high-salt stress, plants were watered with 300 mM NaCl 20 h before tissue collection; for UV stress, the UV stressor was applied for 2 h using a UV-B lamp under growth chamber conditions before tissue collection; T0, after ten days of continuous treatment; T7, after ten days of treatment with normal watering for seven days to recover from stress; Mild_water_deficit indicates moderate water deficit treatment (−0.2 Mp water potential); Severe_water_deficit indicates severe water deficit treatment (−0.8 Mp water potential).
Figure 8. The expression levels of maize ZmCHR gene family members under abiotic stresses are depicted in (A). (A) depicts the expression levels of ZmCHR in the whole above-ground tissues of maize seedlings under cold, heat, salt, and ultraviolet stresses; (B) depicts the expression levels of ZmCHR in maize leaves under different drought and salt stresses; and (C) depicts the expression levels of ZmCHR in the primary roots of maize seedlings under different water deficit stresses. Note: for cold stress, seedlings were incubated at 5 °C for 16 h; for heat stress, seedlings were incubated at 50 °C for 4 h; for high-salt stress, plants were watered with 300 mM NaCl 20 h before tissue collection; for UV stress, the UV stressor was applied for 2 h using a UV-B lamp under growth chamber conditions before tissue collection; T0, after ten days of continuous treatment; T7, after ten days of treatment with normal watering for seven days to recover from stress; Mild_water_deficit indicates moderate water deficit treatment (−0.2 Mp water potential); Severe_water_deficit indicates severe water deficit treatment (−0.8 Mp water potential).
Ijms 25 09686 g008
Figure 9. Quantitative RT-PCR analysis of selected representative genes in the maize ZmCHR family. Vertical bars indicated the standard deviations. The values represent the mean ± standard deviation (SD) of three independent replicates. An asterisk indicates that the corresponding gene is significantly upregulated or downregulated compared with the 0 h status (* p < 0.05, ** p < 0.01, t-test).
Figure 9. Quantitative RT-PCR analysis of selected representative genes in the maize ZmCHR family. Vertical bars indicated the standard deviations. The values represent the mean ± standard deviation (SD) of three independent replicates. An asterisk indicates that the corresponding gene is significantly upregulated or downregulated compared with the 0 h status (* p < 0.05, ** p < 0.01, t-test).
Ijms 25 09686 g009
Table 1. Characteristics of gene structures and protein properties of maize ZmCHR family members.
Table 1. Characteristics of gene structures and protein properties of maize ZmCHR family members.
Gene NameSequence IDLength
(aa)
MW
(kDa)
pIInstability IndexAliphatic IndexGRAVYPredicted Location(s)
ZmCHR1Zm00001eb017000_T0011398159.898.1148.8883.53−0.529Nuclear
ZmCHR2Zm00001eb046650_T0011481168.185.4651.976.91−0.674Nuclear
ZmCHR3Zm00001eb055640_T00283794.165.4645.0792.19−0.417Nuclear
ZmCHR4Zm00001eb073870_T001973109.046.1252.0581.01−0.525Nuclear
ZmCHR5Zm00001eb076670_T0011048116.955.8542.1183.11−0.533Nuclear
ZmCHR6Zm00001eb108690_T0011129128.886.7852.4980.88−0.629Nuclear
ZmCHR7Zm00001eb117870_T00177988.406.2144.790.8−0.461Cytoplasmic
ZmCHR8Zm00001eb127010_T0011198133.407.6448.279.47−0.609Nuclear
ZmCHR9Zm00001eb131790_T0021114128.435.2249.0372.56−0.851Nuclear
ZmCHR10Zm00001eb171780_T0031033114.167.7547.1980.92−0.455Nuclear
ZmCHR11Zm00001eb173050_T002952106.756.6140.8886.13−0.418Nuclear
ZmCHR12Zm00001eb176980_T0011193132.338.5248.0187.04−0.248Nuclear
ZmCHR13Zm00001eb181350_T00141146.278.5664.689.93−0.185Nuclear
ZmCHR14Zm00001eb186260_T0011360151.366.4351.8574.21−0.675Nuclear
ZmCHR15Zm00001eb187690_T0032018231.235.3850.3277.33−0.654Nuclear
ZmCHR16Zm00001eb220010_T00287999.648.6548.9181.95−0.501Nuclear
ZmCHR17Zm00001eb229760_T0012208246.568.8556.4166.3−0.815Nuclear
ZmCHR18Zm00001eb251100_T0031975226.475.2952.1776.95−0.653Nuclear
ZmCHR19Zm00001eb255780_T002968107.386.4249.1781.27−0.46Nuclear
ZmCHR20Zm00001eb280740_T0013828414.735.0752.1971.07−0.647Nuclear
ZmCHR21Zm00001eb297280_T0012036225.905.8744.0197.49−0.188Nuclear
ZmCHR22Zm00001eb322510_T00185095.878.6351.5176.42−0.48Nuclear
ZmCHR23Zm00001eb326710_T0011202135.368.3950.388.78−0.38Nuclear
ZmCHR24Zm00001eb328680_T00171778.849.6650.4894.3−0.153Mitochondrial
ZmCHR25Zm00001eb329060_T00385495.168.6841.6990.91−0.265Nuclear
ZmCHR26Zm00001eb330720_T0031709193.326.6848.4784.96−0.375Nuclear
ZmCHR27Zm00001eb331510_T0011214137.938.0543.4982.06−0.404Nuclear
ZmCHR28Zm00001eb341650_T0021102127.345.2548.0972.99−0.853Nuclear
ZmCHR29Zm00001eb369880_T0031239138.325.1552.2872.32−0.513Nuclear
ZmCHR30Zm00001eb379660_T0022254249.356.0453.384.63−0.444Nuclear
ZmCHR31Zm00001eb394610_T0021380157.948.7946.6482.78−0.553Cytoplasmic
ZmCHR32Zm00001eb404510_T002904101.656.5546.1394.21−0.214Nuclear
ZmCHR33Zm00001eb416780_T001950106.546.6143.0684.57−0.442Nuclear
ZmCHR34Zm00001eb417900_T0021110126.986.0452.8681.82−0.624Nuclear
ZmCHR35Zm00001eb431470_T001979109.846.2951.0978.77−0.464Nuclear
ZmCHR1Zm00001eb017000_T0011398159.898.1148.8883.53−0.529Nuclear
ZmCHR2Zm00001eb046650_T0011481168.185.4651.976.91−0.674Nuclear
ZmCHR3Zm00001eb055640_T00283794.165.4645.0792.19−0.417Nuclear
ZmCHR4Zm00001eb073870_T001973109.046.1252.0581.01−0.525Nuclear
ZmCHR5Zm00001eb076670_T0011048116.955.8542.1183.11−0.533Nuclear
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

Han, J.; Wang, Q.; Qian, B.; Liu, Q.; Wang, Z.; Liu, Y.; Chen, Z.; Wu, W.; Zhang, C.; Yin, Y. Exploring the Roles of the Swi2/Snf2 Gene Family in Maize Abiotic Stress Responses. Int. J. Mol. Sci. 2024, 25, 9686. https://doi.org/10.3390/ijms25179686

AMA Style

Han J, Wang Q, Qian B, Liu Q, Wang Z, Liu Y, Chen Z, Wu W, Zhang C, Yin Y. Exploring the Roles of the Swi2/Snf2 Gene Family in Maize Abiotic Stress Responses. International Journal of Molecular Sciences. 2024; 25(17):9686. https://doi.org/10.3390/ijms25179686

Chicago/Turabian Style

Han, Jiarui, Qi Wang, Buxuan Qian, Qing Liu, Ziyu Wang, Yang Liu, Ziqi Chen, Weilin Wu, Chuang Zhang, and Yuejia Yin. 2024. "Exploring the Roles of the Swi2/Snf2 Gene Family in Maize Abiotic Stress Responses" International Journal of Molecular Sciences 25, no. 17: 9686. https://doi.org/10.3390/ijms25179686

APA Style

Han, J., Wang, Q., Qian, B., Liu, Q., Wang, Z., Liu, Y., Chen, Z., Wu, W., Zhang, C., & Yin, Y. (2024). Exploring the Roles of the Swi2/Snf2 Gene Family in Maize Abiotic Stress Responses. International Journal of Molecular Sciences, 25(17), 9686. https://doi.org/10.3390/ijms25179686

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