Genome-wide identification and multi-dimensional functional characterization of the SIR2 family in Brassica napus L
Hanbing Hu, Ling He, Rui Chen, Ziheng Jing, Heping Wan, Zilan Xiao, Hongtao Cheng, Chunhong Wu, Xiaoyun Liu
TL;DR
This study identifies and characterizes four SIR2 family genes in rapeseed, revealing their roles in epigenetic regulation and stress responses.
Contribution
The study provides the first systematic characterization of the SIR2 family in Brassica napus, including evolutionary origins and functional roles.
Findings
Four SRT genes (BnaASRT1/2 and BnaCSRT1/2) were identified in Brassica napus with conserved SIR2 domains.
BnaSRTs show histone deacetylase activity and bind to acetylation sites like H3K9ac, influencing gene expression.
Expression of BnaSRTs is tissue-specific and responsive to abiotic stresses and hormones.
Abstract
Silent Information Regulator 2 (SIR2) family proteins are NAD+-dependent histone deacetylases that play key roles in epigenetic regulation, plant development, and stress responses; however, systematic characterization of this family (hereafter SRTs) in the important oilseed crop Brassica napus (rapeseed) remains limited. In this study, we identified 4 SRT genes in the rapeseed genome, designated BnaASRT1/2 and BnaCSRT1/2, based on their homology to Arabidopsis SRTs and chromosomal localization. These genes encode proteins with conserved SIR2 domains, sharing >80% sequence identity and >95% coverage with Arabidopsis SRTs, and exhibit similar subcellular localizations (SRT1s in the nucleus, SRT2s in mitochondria). Phylogenetic and synteny analyses confirmed the evolutionary conservation of rapeseed SRTs, which originated from whole-genome duplication events and were stably inherited from…
Genes, proteins, chemicals, diseases, species, mutations and cell lines named across the full text — each resolved to its canonical identifier and authoritative record.
Click any figure to enlarge with its caption.
Fig 1
Fig 2
Fig 3
Fig 4
Fig 5
Fig 6
Fig 7
Fig 8- —Oil Crops Research Institute of Chinese Academy of Agricultural Sciences
- —http://dx.doi.org/10.13039/501100012239Hubei Technological Innovation Special Fund
- —http://dx.doi.org/10.13039/100012554Hubei Provincial Department of Education
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsSirtuins and Resveratrol in Medicine · Plant Gene Expression Analysis · Plant Molecular Biology Research
Introduction
Plants have evolved sophisticated epigenetic regulatory mechanisms to maintain stable growth and development under diverse environmental conditions [1]. Histone modification represents one of the key epigenetic regulatory strategies employed by plants in response to environmental stresses [2,3].
Histone acetylation, among the most common histone modifications, is catalyzed by histone acetyltransferases (HATs) and histone deacetylases (HDACs), with HATs mediating acetylation and HDACs governing deacetylation [4]. The main function of HDACs is to remove acetyl groups from Nε-Acetyl-L-lysine residues on histones, thereby increasing the number of positive charges on the histones and enhancing their ability to bind DNA [5]. This process enhances the accessibility of chromatin to transcriptional regulators, thereby promoting gene expression. In contrast, the deacetylation of histones by HDACs leads to chromatin compaction, thereby inhibiting gene expression [6]. Three HDAC families have been identified in plants: Class I, designated as the Reduced Potassium Dependency 3 (RPD3)/Histone Deacetylase 1 (HDA1) family; Class II, the Silent Information Regulator 2 (Sirtuin/SIR2) family; and Class III, the HD2/HDAC tuin (HD-tuin/HDT) family [7,8]. The SIR2 gene family is a group of highly conserved NAD ⁺ -dependent histone deacetylases in eukaryotes, with vital effects on plant growth, development, and their responses to stress conditions. [9,10].
In plants, extensive research has been documented regarding the involvement of the SIR2 gene family in the regulation of energy metabolism and hormone signaling pathways, and these two processes further adjust plant growth and development along with responses to abiotic and biotic stresses. Through histone deacetylation, AtSRT1 directly suppresses the expression of FT and SOC1. Simultaneously, it directly downregulates TOR (a central hub of glucose-driven energy signal transduction), thereby promoting the deposition of the H3K27me3 marks on FLC through the TOR-FIE-PRC2 pathway, ultimately enhancing the suppression of flowering [9]. In rice, OsSRT1 inhibits glycolysis by deacetylating OsGAPDH, increases glucose content in rice plants, and regulates starch accumulation in seeds [11]. Furthermore, OsSRT1 negatively regulates leaf senescence by inhibiting the expression of biosynthetic genes involved in the MeOH-jasmonates cascade metabolic pathway [12]. AtSRT2 exhibits a primary localization in mitochondria, where it exerts a crucial effect on the fine-tuning of energy metabolism within this organelle [13]. Disruption of AtSRT2 results in decreased ATP content while causing elevated ADP levels in the leaves of Arabidopsis [14]. AtSRT1 and AtSRT2 interact with EIN2 NUCLEAR ASSOCIATED PROTEIN1 (ENAP1), and this protein complex exerts a mediating role in transcriptional repression via the regulation of H3K9 acetylation levels within the ethylene signaling transduction pathway [15]. PbSRT1 and PbSRT2 play important roles in regulating pear fruit growth and ripening by modulating hormone and sugar metabolism [16]. FaSRT2-1 plays a significant role in regulating energy metabolism, promoting the maturation, aging, and quality formation of strawberry fruits [17].
In Arabidopsis, AtSRT1 negatively controls the stress tolerance of plants as well as the glycolytic process; concurrently, it stimulates mitochondrial respiration by means of epigenetic regulation and adjusting the transcriptional activity of AtMBP-1 [18]. AtSRT2 promotes plant growth and represses osmotic stress-associated signaling [10]. Additionally, AtSRT2 may negatively regulate basal resistance against the pathogen Pseudomonas syringae pv. tomato DC3000 (Pst DC3000) and the expression of pathogenesis-related gene 1 (PR1) by inhibiting salicylic acid (SA) biosynthesis [19]. Under salt treatment, AtSRT2 modulates histone acetylation in the promoter region of Vesicle-associated membrane protein 714 (VAMP714) and represses its transcription, thereby facilitating the trafficking of H₂O₂-containing vesicles in cells and enhancing salt tolerance [20]. In Oryza sativa (rice), RNA interference (RNAi)-mediated silencing of OsSRT1 induces increased histone H3K9 acetylation and decreased H3K9 dimethylation, which in turn triggers H₂O₂ production, DNA fragmentation, cell death, and hypersensitivity reactions. By contrast, OsSRT1 overexpression enhances oxidative stress tolerance [21]. In strawberry, FaSRT1-2 positively promotes vegetative growth and fruit ripening by regulating the ABA and GA pathways, while negatively modulating resistance to heat stress and Botrytis cinerea infection via controlling the expression of relevant genes (e.g., heat shock proteins, transcription factors and stress-related genes) [22]. DgSRT2 enhances cold tolerance in chrysanthemum via the reactive oxygen species (ROS) scavenging system [23]. Recently, significant breakthroughs have been made in functional studies of TaSRT1 in wheat. Research indicates that TaSRT1 directly interacts with the key drought resistance factor TaDT-A, regulating the protein stability and transcriptional activity of TaDT-A through lysine deacetylation, thereby influencing the drought resistance of wheat [24]. Furthermore, TaSRT1 negatively regulates Fusarium crown rot (FCR) resistance by inhibiting the accumulation of TaHSP18.6 via deacetylation. TaHSP18.6 interacts with TaIAA1 to regulate resistance by altering endogenous auxin content, thus proposing a TaSRT1-TaHSP18.6 regulatory model for wheat FCR resistance [25].
Brassica napus L. (rapeseed), the second largest oilseed crop worldwide, possesses an allotetraploid genome. Extensive research indicates that the SIR2 family plays a crucial role in plant growth, development, and stress responses. However, the SIR2 family in rapeseed remains inadequately characterized. Thus, it is essential to elucidate their functional properties in this species.
In the present study, four SIR2 family members were identified in the rapeseed genome, followed by comprehensive analyses encompassing phylogeny, conserved domains, gene structures, and cis-acting elements. Furthermore, the biochemical functions of BnaSRTs were predicted via molecular docking and validated using sirtuin inhibitor. Their expression profiles were also investigated across different developmental stages of various rapeseed organs, as well as in responses to abiotic stresses and hormone treatments. In summary, this study conducted a systematic analysis of the SRTs in rapeseed, laying a solid foundation for deciphering their regulatory roles in rapeseed growth and development, stress responses, and hormonal signal transduction.
Materials and methods
Characterization and Categorization of SRTs
Following previous research, the SRTs protein sequences of Arabidopsis thaliana and Oryza sativa were downloaded from the Plant Comparative Genomics portal (Phytozome, https://phytozome-next.jgi.doe.gov/) [26] and the sequences of Homo sapiens were downloaded from the National Center for Biotechnology Information (NCBI, https://www.ncbi.nlm.nih.gov/) [27]. The Arabidopsis SRTs protein sequences were used as queries for BLAST analysis on the Brassica napus multi-omics information resource (BnIR, https://yanglab.hzau.edu.cn) [28], to preliminarily identify SRTs in the Brassica napus Zhongshuang11. Gene IDs are available in Supplementary S1 Table. Conserved protein domains of the target gene family in Brassica napus were analyzed via the NCBI Conserved Domain Database (CDD, https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) [29].
Physicochemical properties of proteins in the target gene family were analyzed via the ExPASy online software ProtParam (https://web.expasy.org/protparam/) [30], with parameters encompassing amino acid count, theoretical pI, molecular weight, instability index, aliphatic index, and grand average of hydropathicity. The CDS length, number of introns, number of exon and chromosome location of each confirmed SRTs was retrieved from the BnIR. Amino acid sequences of two AtSRTs and four BnaSRTs were used to predict the sub-cellular localization for plant proteins by ProtComp 9.0 on the Softberry website. (http://www.softberry.com/berry.phtml?topic=protcomppl&group=programs&subgroup=proloc).
Protein structures prediction
The structure of the BnaSRTs proteins were predicted and modeled using the SWISS-MODEL online software (https://swissmodel.expasy.org/interactive) [31], and then visualized using PyMOL 2.6 (https://pymol.org/) [32]. Two AtSRTs and four BnaSRTs protein sequences were simultaneously aligned with Clustal Omega (https://www.ebi.ac.uk/jdispatcher/msa/clustalo) [33], and labeled as SIR2 domains, that can be found in S1 Fig. The secondary structure was visualized in PyMOL 2.6, consisting of three structures: Helix, Sheet and Loop, and conserved domains are also displayed based on the results of the analysis on the NCBI’s CDD.
Phylogenetic analysis of SRTs proteins in four species
The identified Brassica napus SRTs protein sequences and those of three other species, including Homo sapiens, Arabidopsis thaliana and Oryza sativa, were downloaded, and the protein sequences were simultaneously aligned with MUSCLE. The phylogenetic tree was constructed using the Molecular Evolutionary Genetics Analysis 11 (MEGA11, https://www.megasoftware.net/) [34] by the neighbor-joining (NJ) method, and a bootstrap analysis was conducted using 1000 replicates at each node. The phylogenetic tree was visualized and optimized by Interactive Tree Of Life (iTOL, https://itol.embl.de/) [35]. Different colors were used to indicate the homologous genes of two AtSRTs, and different shapes were used to label the SRTs proteins of different species.
Chromosomal distribution and synteny analysis
The genomic data of Arabidopsis thaliana, Brassica napus, Brassica rapa and Brassica oleracea were obtained from the Arabidopsis Information Resource (TAIR, https://www.arabidopsis.org/) [36] and BnIR, respectively. The gene location analysis and gene density analysis were performed using the Gene Location Visualize form GTF/GFF tool and Gene Density Profile tool in TBtools-Ⅱ v2.149 (https://github.com/CJ-Chen/TBtools) [37]. The synteny analysis was performed using the One Step MCScanX tool, Dual Systeny Plot tool, File Merge For MCScanX tool and Multiple Synteny Plot tool in TBtools-Ⅱ.
Sequence analysis
The genomic data of Homo sapiens and Oryza sativa were obtained from the GenBank in NCBI (https://www.ncbi.nlm.nih.gov/genbank/) [38] and Rice Genome Annotation Project (RGAP, https://rice.uga.edu/) [39], respectively. The conserved protein domain SIR2 was searched within protein sequences of SRTs using the NCBI’s CDD. The gene structure of SRTs were assessed utilizing information from tree plant genome GFF3 file. And the motif consensus can be found in the S2 Table. The conserved features of SRTs protein sequences of four species were analyzed using the Multiple Em for Motif Elicitation v5.5.8 (MEME Suite, https://meme-suite.org/meme/tools/meme) [40]. The domains, gene structures and motifs were visualized using the Gene Structure View (Advanced) tool in TBtools-Ⅱ. The cis-elements in the promoters of four Brassica napus SRTs were identified using the PlantCare (https://bioinformatics.psb.ugent.be/webtools/plantcare/html/) [41], and visualized using the Gtf/Gff3 Sequences Extractor tool, the Fasta Stater tool, the Fasta Sequence Manipulator tool and the Simple BioSequence Viewer tool in TBtools-Ⅱ.
GO enrichment analysis and molecular docking prediction
Based on the protein sequences of BnaSRTs, protein GO enrichment analysis was performed using STRING online software (https://cn.string-db.org/) [42]. Due to the lack of relevant information on the ZS11 rapeseed gene in the STRING, information on the Darmor rapeseed gene, which is homologous to the ZS11 rapeseed gene, was used. Details of GO and KEGG pathway enrichment analysis can be found in S3 Table.
5 types of ligand-modified histone acetylation sites (H3K9ac, H3K14ac, H3K23ac, H3K27ac and H4K5ac) were downloaded from the RCSB Protein Data Bank (RCSB PDB, https://www.rcsb.org/) [43]. Receptor and ligand proteins were subjected to operations such as dehydration, hydrogenation, and energy minimization using AutoDock4 (https://autodock.scripps.edu/) [44], and molecular docking was performed using AutoDock Vina (https://vina.scripps.edu/) [45] to calculate the results. The results were visualized using PyMOL and the online software PoseEdit (https://proteins.plus/) [46].
Treatment of rapeseed with histone deacetylase inhibitors, RNA-seq and identification of histone modification levels
The ZS11 rapeseed seeds were germinated by soaking in clean water, and after 3 days, they were transferred to a modified Hogue’s nutrient solution for rapeseed cultivation in a greenhouse. Inside the greenhouse, a constant temperature of 20°C, constant humidity, 20 hours of illumination, and 4 hours of darkness. After 7 days, a portion of the rapeseed was used as a control group, while another portion was treated with Nicotinamide, one of the histone deacetylase inhibitors (with a working concentration of 50μM) [11]. On the fifth day post-treatment, significant phenotypic changes were observed. These samples are first cleaned with ultra-pure water and then preserved in liquid nitrogen.
Total RNA was extracted using TRIzol reagent (Invitrogen Life Technologies). The RNA was subsequently evaluated and quantified in two steps: 1) its purity and concentration were measured using a NanoDrop spectrophotometer (Thermo Fisher Scientific) and a Qubit 4.0 fluorometer (Thermo Fisher Scientific); 2) the integrity and quantity of RNA were assessed using the Agilent 2100/4200 system. Three micrograms of RNA were used as the starting material for subsequent sample preparation. cDNA libraries were constructed through Poly(A⁺) enrichment and using Hieff NGS® DNA Selection Beads. After library construction, the concentration was measured again with Qubit 4.0. The libraries were then sent to 563 Sci‑Service Biotechnology Co., Ltd. (Wuhan, China) for sequencing on the DNBSEQ‑T7 platform using the PE150 mode. The RNA‑seq experiment was conducted with three biological replicates. The raw FASTQ data were first subjected to quality control using fastp (v0.21.0), filtering out low-quality sequences to obtain high-quality clean reads. The clean reads were then aligned to the reference genome using HISAT2 (v2.1.0). Gene read counts were quantified using StringTie (v2.1.5), and expression levels were normalized as FPKM (fragments per kilobase of transcript per million mapped reads). Differentially expressed genes (DEGs) were identified using DESeq2 (v1.34.1) with the criteria of |log₂FoldChange| > 1 and Padj ≤ 0.05. Finally, the identified DEGs were submitted to the GO and KEGG databases for functional annotation (Padj ≤ 0.05) and their functions as well as metabolic pathways were analyzed and predicted. GO and KEGG pathway enrichment analyses for differentially expressed genes, histone‑modifying genes, and genes whose expression is regulated by these modifications were carried out using AgriGOv2. (http://systemsbiology.cau.edu.cn/agriGOv2/) [47] and KOBAS software [48] (S7 Table).
Additionally, total proteins extracted from the samples were used for Western Blotting to detect the changes in the levels of five histone modifications: H3K9ac (A7255, ABclonal), H3K14ac (A7254, ABclonal), H3K23ac (A18154, ABclonal), H3K27ac (A22264, ABclonal), and H4K5ac (A23080, ABclonal), in order to identify the biochemical functions of SRTs in rapeseed.
Expression profiling of SRTs in rapeseed
The expression pattern of SRTs in rapeseed following tissue, abiotic stress treatment and hormone stress treatment were obtained from BnIR. The study types were as follows: (1) tissues of leaf and root: bud, filament, petal,pollen, sepal, cotyledon, vegetative rosette, leaf, root, seed, silique, and so on. For detailed information, please refer to S4 Table. (2) abiotic stress treatment of leaf and root: salt, drought, freezing, cold, heat, and osmotic. For detailed information, please refer to S5 Table. (3) hormone treatment of leaf and root: Indole-3-acetic acid (IAA), 1-Aminocyclopropane-1-carboxylic Acid (ACC), Gibberellin (GA), Abscisic acid (ABA), Trans-Zeatin-Riboside (TZ), Jasmonic acid (JA), and Brassinolide (BL). For detailed information, please refer to S6 Table. RNA-seq data was downloaded from the BnIR, with the gene expression quantification method based on Kallisto, and gene expression was estimated using the transcripts per kilobase of exon model per million mapped reads (TPM) method.
Validation of transcriptome results and expression patterns via qRT-PCR
Total RNA was extracted from rapeseed materials treated with NAM (50 μM NAM, 5 days), salt stress (1.3% NaCl, 2 days), and drought stress (20% PEG6000, 1 day) using Trizol reagent. From the extracted total RNA, 1 µg was reverse‑transcribed into cDNA using the ABScript II First‑Strand cDNA Synthesis Kit (RK20400, ABclonal). The resulting first‑strand cDNA was then used as the template for the subsequent qRT‑PCR assays.
qRT‑PCR was performed on an ABI 7900 instrument. The cycling protocol began with an initial denaturation at 95°C for 15 seconds, followed by 42 amplification cycles consisting of 95°C for 5 seconds and 60°C for 40 seconds.
For melt‑curve analysis, the program included a step at 95°C for 15 seconds, then 60°C for 20 seconds, and a final step at 95°C for 15 seconds. Data acquisition and analysis were carried out using the ABI 7900 system. Relative gene expression was calculated with the 2^⁻ΔΔCT^ method, and BnActin7 (EV116054) from rapeseed was used as a reference gene. The primers employed in these analyses are presented in S8 Table.
Results
Identification of SRTs in rapeseed
To identify the SRT genes in the rapeseed (Brassica napus L.) genome, we performed a comprehensive search using Arabidopsis SRTs sequences as queries. The results revealed four SRT genes harboring the SIR2 domain in the rapeseed genome (Table 1). Based on their chromosomal locations and homology to Arabidopsis SRTs, these genes were designated as BnaASRT1/2 and BnaCSRT1/2, respectively. Physicochemical properties of the identified genes were analyzed via the ProtParam server. Our analyses showed that the lengths of rapeseed SRTs proteins ranged from 357 to 484 amino acids (aa); their molecular weights varied between 39.63 and 54.20 kDa; and their isoelectric points (pI) spanned from 6.62 to 8.96. Specifically, BnaASRT1 and BnaCSRT1 were weakly acidic proteins with pI > 6, whereas BnaASRT2 and BnaCSRT2 were basic proteins with pI > 8. The instability index ranged from 42.3 to 52.25; the aliphatic index was between 78.46 and 87.03; and the grand average of hydropathicity (GRAVY) varied from −0.246 to −0.316, indicating that BnaSRTs proteins exhibit strong hydrophilicity (Table 1). Prediction results suggested that the corresponding SRTs genes on the A and C subgenomes of rapeseed share identical subcellular localization and display similarities in isoelectric point, aliphatic index, and GRAVY. BnaASRT1, BnaCSRT1, and AtSRT1 were all localized in the nucleus, while BnaASRT2, BnaCSRT2, and AtSRT2 were all targeted to mitochondria (Table 1). Additionally, the SRTs of rapeseed and Arabidopsis showed similarities in molecular weight, instability index, aliphatic index, and GRAVY (Table 1).
Table 1: Physicochemical properties of SIR2 gene family members in Arabidopsis thaliana and Brassica napus. No.aa, Number of Amino Acid; TPI, Theoretical pI; MW, Molecular Weight; II, Instability Index; AI, Aliphatic Index; GRAVY, Grand Average of Hydropathicity; PSL, Prediction of Subcellular Localisation; N, Nuclear; M, Mitochondrial.
Protein sequence alignments were conducted between AtSRT1 and BnaASRT1/BnaCSRT1, as well as between AtSRT2 and BnaASRT2/BnaCSRT2, respectively. The alignment results demonstrated that the coverage (cov) of SRT proteins between rapeseed and Arabidopsis exceeded 95%, with a percent identity (pid) of over 80%, indicating that rapeseed and Arabidopsis SRTs share a common evolutionary origin and maintain a close sequence relationship (S1 Fig). AlphaFold was used to predict the tertiary structures of BnaSRTs proteins (Fig 1), and their secondary structural elements—including helices, sheets, and loops—were visualized using PyMOL (Fig 1A). Furthermore, the positions of the SIR2 domains, as predicted by CDD, were mapped onto the tertiary structures of BnaSRTs (Fig 1B). Structural predictions revealed a high degree of similarity in both secondary and tertiary structures between AtSRT1 and BnaASRT1/BnaCSRT1, as well as between AtSRT2 and BnaASRT2/BnaCSRT2. These findings suggest that SRT proteins are highly conserved between Brassica napus and Arabidopsis, two species belonging to the Brassicaceae family.
Protein structures prediction.(A) The secondary structure of the protein is shown within the tertiary structure. The secondary structures are distinguished by blue helices, purple sheets, and pink loops. (B) In the tertiary structure of the protein, the SIR2 conserved domain is displayed in green.
Chromosomal distribution and phylogenetic relationship analysis of BnaSRTs genes
BnaASRT1 and BnaCSRT1 are localized on chromosomes A03 and C03 of Brassica napus, respectively, while BnaASRT2 and BnaCSRT2 are situated on chromosomes A02 and C02, respectively. All these genes are positioned at the distal ends of their respective chromosomes (Fig 2A).
Chromosomal Distribution, Synteny Analysis of SRTs Genes, and Phylogenetic Analysis of SRTs.(A) Chromosomal distribution of BnaSRTs in rapeseed. The gene density heat map within the chromosome rectangles indicates gene density from low (cool colors) to high (warm colors). (B) Phylogenetic analysis of 15 SRT proteins based on Homo sapiens (7), Arabidopsis thaliana (2), Brassica napus (4), and Oryza sativa (2). A neighbor-joining phylogenetic tree was constructed, and bootstrap consensus was derived from 1000 replicates. The red branches represent proteins that are evolutionarily related to AtSRT1, while the blue branches represent proteins that are evolutionarily related to AtSRT2. Each shape represents a distinct species. The rect is Homo sapiens; the triangle is Arabidopsis thaliana; the star is Brassica napus; the circle is Oryza sativa. (C) Synteny analysis of SRTs of Arabidopsis thaliana and Brassica napus. (D) Comparative synteny of SRT genes among Brassica napus, Brassica rapa and Brassica oleracea. Background gray connections indicate whole-genome duplication events, while the collinearity of the four colors is highlighted, specifically tracking the conservative SRTs gene pairs.
To elucidate the evolutionary relationships within the SRT genes, a neighbor-joining phylogenetic tree was constructed based on multiple sequence alignments of 15 SRT proteins from Homo sapiens (7), Arabidopsis thaliana (2), Brassica napus (4), and Oryza sativa (2) (Fig 2B). The results revealed that AtSRT1, BnaASRT1, BnaCSRT1, and OsSRT701 share a close evolutionary relationship with HsSIRT6 and HsSIRT7; in contrast, AtSRT2, BnaASRT2, BnaCSRT2, and OsSRT702 are evolutionarily closely related to HsSIRT4 (S1 Table).
In order to gain a deeper understanding of the origin, evolutionary history, and functional characteristics of the SRTs genes in rapeseed, a synteny sequence analysis of the SRTs genes in Brassica napus and Arabidopsis thaliana was conducted (Fig 2C). BnaASRT1 and BnaCSRT1 are homologous to AtSRT1, while BnaASRT2 and BnaCSRT2 are homologous to AtSRT2, which is consistent with the evolutionary relationship analysis. These findings indicated that SRTs genes were not lost during homologous evolution and were stably inherited throughout the polyploidization process of Brassica napus.
To investigate the amplification mechanism of BnaSRTs genes in Brassica napus (B. napus), synteny analysis of SRTs genes was performed across the genomes of B. napus, Brassica rapa (B. rapa), and Brassica oleracea (B. oleracea). The analysis revealed that the four BnaSRTs in B. napus exhibit collinearity with two genes in B. rapa and two in B. oleracea, respectively, with all BnaSRTs mapped to chromosomes A02, A03, C02, and C03. This observation suggests that BnaSRTs are derived from the genomes of B. rapa and B. oleracea, and the SRTs gene has experienced chromosomal fragment polyploidization during its evolutionary process (Fig 2D). These results reveal a conserved intergeneric collinear relationship of BnaSRTs.
Sequence and structure analysis of SRTs
Proteins of BnaSRTs, AtSRTs, and OsSRTs all harbor the conserved SIR2 domain (Fig 3A). Gene structure analysis revealed that SRTs genes from these three plant species share similar architectural features (Fig 3B). AtSRT1, BnaASRT1, BnaCSRT1, and OsSRT701 exhibit identical numbers of exons and introns, with 14 exons and 13 introns, respectively. AtSRT2 and OsSRT702 also share consistent exon-intron counts, containing 11 exons and 10 introns. In contrast, BnaASRT2 and BnaCSRT2 possess 9 exons and 8 introns (Fig 3B).
BnaSRTs features, phylogeny and cis-element analysis.(A) SIR2 conserved domain of BnaSRTs, AtSRTs and OsSRTs. (B) exon-intron structures of BnaSRTs, AtSRTs and OsSRTs. Green rectangles correspond to the coding sequence (CDS) or exons, yellow rectangles to untranslated regions (UTRs), and grey lines to introns. (C) The phylogenetic analysis and motif composition of BnaSRTs, AtSRTs, OsSRTs and HsSIRTs, and motif logo. (D) Analysis of Cis-Elements in BnaSRTs Promoters: Distribution and Quantity of Cis-Elements per Promoter.
The protein sequence characteristics of HsSIRTs, AtSRTs, BnaSRTs, and OsSRTs were analyzed using MEME (Fig 3C). All these SRTs proteins contain the conserved motif array 2-1-5-4-3. The motif composition of SRT genes from the three plant species is consistent with the predicted results of their conserved SIR2 domains (S2 Table), indicating that these SRTs proteins may encode functional homologs with similar biological functions.
Promoter elements located in the 2000-bp region upstream of the CDS of BnaSRTs genes were predicted using PlantCARE, and the types of promoter elements were quantified (Fig 3D). Most of these promoter elements are associated with stress responses and the regulation of plant growth and development. The promoter elements of the four BnaSRTs genes display substantial divergence. Although all contain light-responsive, methyl jasmonate (MeJA)-responsive, and abscisic acid (ABA)-responsive elements, their abundances vary significantly. Some promoter elements are exclusive to a single or a specific subgroup of BnaSRT genes: for instance, low-temperature-responsive elements are restricted to BnaSRT1; circadian rhythm regulation and endosperm expression elements are specific to BnaSRT2; defense and stress-responsive elements, as well as zein metabolism regulatory elements, are unique to BnaSRT1. Hormone-responsive elements—including MeJA, gibberellin, ABA, auxin, and salicylic acid responsive ones—are more abundant in BnaSRT2 (Fig 3D).
Analysis of BnaSRTs GO enrichment and interaction with histone acetylation
GO enrichment analysis revealed that, in biological processes, BnaSRTs are involved not only in histone deacetylation but also in histone H3K9 acetylation, ethylene-activated signaling pathways, and negative regulation of biological processes (Fig 4A). Specifically, BnaSRT1s additionally participate in the negative regulation of RNA polymerase II-mediated transcription, while BnaSRT2s are involved in protein ADP-ribosylation and negative regulation of defense responses (S3 Table). These observations suggest that BnaSRTs may negatively regulate DNA transcription and contribute to the modulation of plant defense responses. In terms of molecular functions, BnaSRTs possess protein lysine deacetylase activity, NAD ⁺ -binding capacity, and metal ion-binding activity (Fig 4B). Notably, BnaSRT2s also exhibit chromatin DNA-binding activity (S3 Table), indicating their potential involvement in chromatin DNA modification. KEGG pathway enrichment analysis further demonstrated that BnaSRTs are implicated in nicotinic acid and nicotinamide metabolism (S2 Fig).
BnaSRTs GO enrichment and molecular docking to 5 histone acetylation sites.(A) The top 8 significantly enriched GO biological processes of BnaSRTs. (B) The top 7 significantly enriched GO molecular function of BnaSRTs. The horizontal axis represents the error discovery rate, while the vertical axis describes the GO terms. The size of the circles indicates the quantity of genes. (C) Each pair of binding models showed the affinity of 9 binding sites, and the binding pattern was analyzed with the binding site with the largest absolute affinity value. The different shapes and abscissa represent the 5 modification sites (H3K9ac, H3K14ac, H3K23ac, H3K27ac and H4K5ac), while the ordinate represents the affinity (kcal/mol).
Molecular docking predictions showed that BnaSRTs display strong binding activity to five common histone acetylation sites (H3K9ac, H3K14ac, H3K23ac, H3K27ac, H4K5ac) (S3 Fig). Among the predicted binding sites, the highest-affinity interactions exhibit binding energies all below −6.5 kcal/mol (Fig 4C). A binding energy < −5 kJ/mol indicates a measurable binding activity between the receptor and ligand, with lower binding energies corresponding to more stable docking conformations. These results suggest that BnaSRTs are likely to interact strongly with these histone acetylation sites in vivo and may mediate deacetylation at these positions. Among the five acetylation sites, most of the four BnaSRTs show high binding capacity to H3K9ac and H3K27ac, except for BnaCSRT2, which exhibits weak binding to H3K27ac. Additionally, the binding capacity of BnaSRT2s to H4K5ac is lower than that to other H3 sites. The active pockets of BnaSRTs exhibit high diversity with diverse residue combinations; however, certain residues are highly conserved, including Glu, Thr, Ser, Asp, and Asn (S4–S7 Figs).
Treatment of rapeseed with histone deacetylase inhibitors, differentially expressed gene analysis and biochemical functional identification
To preliminarily validate the functions of SRTs proteins in rapeseed, rapeseed plants were treated with the SIRT-specific inhibitor nicotinamide (NAM) [49] for 5 days. It was observed that the seedlings treated with NAM showed weaker growth compared to the control group (Fig 5A). Based on the assessment of histone modification levels, compared with the control (CK), NAM treatment significantly upregulated the levels of H3K9ac, H3K14ac, H3K23ac, H3K27ac, and H4K5ac modifications. This indicates a robust inhibitory effect on histone deacetylases (Fig 5B, S8 Fig), consistent with the molecular docking results.
Rapeseed NAM inhibitor treatment: phenotype, biochemistry, expression and qRT-PCR validation.(A) Phenotype of rapeseed treated with NAM for 7 days, compared with control group (CK). (B) Western blotting identified the difference in histone modification levels between the inhibitor-treated group and the control group. The modified designs were H3K9ac, H3K14ac, H3K23ac, H3K27ac and H4K5ac. Each type of modification was conducted in duplicate. (C) Volcano plots comparing the gene expression levels for the 5-day-NAM-treated rapeseed versus untreated rapeseed. The horizontal axis represents the fold change in gene expression (log2FoldChange) between the treatment and control groups, while the vertical axis indicates the significance level of the expression difference (-log10(FDR)). Red points correspond to upregulated genes, and green points correspond to downregulated genes. Perform GO functional (D) and KEGG pathway (E) enrichment analysis on the differential gene set. 10 significant terms were selected for presentation. (F) In the GO functional enrichment analysis, 10 entries related to stress response were selected for presentation. The x-axis represents the ratio of the number of differential genes annotated to each GO/ KEGG term to the total number of differential genes, while the y-axis indicates the GO/ KEGG term. The size of the points represents the number of genes annotated to the GO/ KEGG term, and the color gradient from red to blue reflects the degree of significance of the enrichment. (G) Results of RT-PCR for NAM-treated and CK rapeseed. Four genes related to stress or hormone response functions were selected from GO functional enrichment analysis for qRT-qPCR validation. The four genes are BnaA08G0189400ZS (osmotic), BnaA03G0358900ZS (salt), BnaA05G0333200ZS (heat), and BnaA06G0133300ZS (JA). The horizontal axis represents the materials, the vertical axis represents gene expression levels (2^-ΔCq average value), and the subheadings indicate the functional type of the gene tested.
RNA-seq analysis revealed that, relative to the control group, the inhibitor-treated group exhibited 12,272 upregulated genes and 9,519 downregulated genes (Fig 5C). GO functional enrichment analysis of the differentially expressed gene (DEG) set showed significant enrichment in biological processes (BP) such as photosynthesis, hydrogen peroxide catabolic process, cellular amino acid catabolic process, and chlorophyll metabolic process; in cellular components (CC) including photosystem, oxidoreductase complex, and photosystem I reaction center; and in molecular functions (MF) including antioxidant activity, oxidoreductase activity, flavin adenine dinucleotide (FAD) binding, and chlorophyll binding (Fig 5D). KEGG pathway enrichment analysis of DEGs demonstrated significant enrichment in pathways such as glycolysis/gluconeogenesis, glyoxylate and dicarboxylate metabolism, carbon fixation in photosynthetic organisms, photosynthesis, various amino acid metabolisms, and the pentose phosphate pathway (Fig 5E). Taken together, these findings suggest that BnaSRTs could participate in the regulation of photosynthesis and respiration. In addition, after NAM treatment, a large number of DEGs were enriched in stress response pathways, such as osmotic stress response, salt stress response, and high-temperature response (Fig 5F). From the results of GO functional enrichment analysis, we selected 4 genes with stress or hormone response functions (BnaA08G0189400ZS (osmotic), BnaA03G0358900ZS (salt), BnaA05G0333200ZS (heat), and BnaA06G0133300ZS (JA)) for qRT-PCR validation (Fig 5G). The validation results were consistent with those of the transcriptome. The above results indicate that the BnaSRTs gene family can regulate the stress response process of plants by influencing the level of histone deacetylation.
Expression profiles of BnaSRTs in various organs, tissues and developmental stage
This study analyzed the expression profiles of BnaSRTs across diverse tissues and developmental stages of rapeseed using the publicly available transcriptome dataset from BnIR (https://yanglab.hzau.edu.cn/BnIR/) (Fig 6, S4 Table). The results demonstrated that BnaSRTs genes are expressed in various vegetative and reproductive organs throughout different developmental phases of rapeseed. Overall, the expression levels of BnaASRT2 and BnaCSRT2 in rapeseed tissues were significantly higher than those of BnaASRT1 and BnaCSRT1. Specifically, BnaCSRT1 exhibited higher expression than BnaASRT1 in most tissues, whereas BnaCSRT2 expression was lower than that of BnaASRT2 in nearly all tissues examined. BnaASRT2 and BnaCSRT2 were predominantly expressed in buds, flowers, leaves, and siliques; in contrast, BnaASRT1 and BnaCSRT1 showed preferential expression in buds, seeds at stages before 34 days, and siliques within 12 days post-flowering. Notably, BnaASRT2 displayed higher expression levels than BnaCSRT2 in petals, sepals, and cotyledons, suggesting that BnaASRT2 may play a more critical role in regulating flowering processes. Additionally, the expression levels of all BnaSRT genes in seeds exhibited a down-regulatory trend over time.
Expression profiles of the rapeseed BnaSRTs in various organs, tissues and developmental stages.The gradient of colors from blue to yellow represents the degree of the expression value [log2(TPM + 1)] provided by the color scale.
Expression profiles of BnaSRTs under different abiotic stresses
Based on the characterized functions of SRTs in other species and promoter element predictions, BnaSRTs are hypothesized to play critical roles in stress responses. In this study, transcriptome dynamics of BnaSRTs under diverse abiotic stresses (salt, drought, freezing, cold, heat, and osmotic stress) were analyzed using publicly available datasets from BnIR (https://yanglab.hzau.edu.cn/BnIR/) (Fig 7, S5 Table, and S9 Fig).
Expression profiles of the rapeseed BnaSRTs under different abiotic stresses, and qRT-PCR validation.These expression profiles pertain respectively to the leaves (A) and the roots (B) at different time points under different abiotic stresses (salt, drought, freezing, cold, head and osmotic). The gradient of colors from blue to yellow represents the degree of the expression value [log2(TPM + 1)] provided by the color scale. qRT-PCR validation was performed using materials from high salt treatment (1.3% NaCl, 2 days) (C), drought treatment (20% PEG, 1 day) (D), and untreated rapeseed. The horizontal axis represents the treatment method and location of the material, the vertical axis represents the gene expression levels (2^-ΔCq average value), and the subtitles indicate the genes being measured.
In leaf tissues, high-salt stress induced a time-dependent downregulatory trend in the expression of BnaASRT2 and BnaCSRT2. Under drought stress, BnaASRT2 and BnaCSRT2 expression exhibited a biphasic pattern: initial upregulation followed by downregulation over time, with a critical turning point at 3 h. Cold stress triggered a similar biphasic trend (initial upregulation then downregulation) for both gene subgroups, with BnaASRT1/BnaCSRT1 and BnaASRT2/BnaCSRT2 reaching their turning points at 6 h and 12 h, respectively (Fig 7A, S5 Table). Osmotic stress consistently suppressed the expression of all BnaSRTs over time (S9 Fig, S5 Table).
In root tissues, high-salt stress led to a progressive downregulation of BnaASRT1 and BnaCSRT1. Drought stress elicited distinct responses: BnaASRT1, BnaASRT2, and BnaCSRT2 showed early upregulation (peaking at 3 h) followed by downregulation, whereas BnaCSRT1 exhibited the opposite pattern with a turning point at 6 h. Under cold stress, BnaASRT1/BnaCSRT1 expression displayed a recovery trend (initial downregulation then upregulation), in contrast to the sustained downregulation of BnaASRT2/BnaCSRT2 (Fig 7B, S5 Table). Freezing stress induced a transient surge at 6 h in BnaASRT1, BnaASRT2, and BnaCSRT2 (following initial downregulation), while BnaCSRT1 remained downregulated throughout. Heat stress triggered a triphasic response in BnaASRT1/BnaCSRT1: downregulation, upregulation,(at 6 h) then re-downregulation. Osmotic stress led to consistent downregulation of BnaASRT1 and BnaCSRT1 over time (S9 Fig, S5 Table). The detection of the expression levels of BnaSRTs under different stress conditions also validated these results (Figs 7C and 7D). These results indicated that BnaSRTs responded to different stresses with distinct patterns, suggesting that BnaSRTs may play regulatory roles under various stress conditions.
Expression profiles of BnaSRTs in response to different hormone treatment
Given the identification of multiple hormone-responsive elements in the BnaSRT promoters during our promoter element analysis, we further examined the expression dynamics of BnaSRTs in rapeseed following treatments with various hormones (IAA, ACC, GA, ABA, TZ, JA, and BL) using publicly available transcriptome datasets from BnIR (Fig 8, S6 Table and S10 Fig).
Expression profiles of the rapeseed BnaSRTs in response to different hormone treatment.These expression profiles pertain respectively to the leaves (A) and the roots (B) at different time points in response to different hormone treatment (IAA, ACC, GA, ABA, TZ, JA and BL). The gradient of colors from blue to yellow represents the degree of the expression value [log2(TPM + 1)] provided by the color scale.
In leaf tissues, IAA treatment induced a biphasic expression pattern in BnaASRT2 over time: initial upregulation followed by downregulation, with a critical turning point at 3 h. Under GA treatment, BnaASRT2 and BnaCSRT2 expression exhibited a sustained downregulatory trend over time. ABA treatment triggered a progressive upregulation of BnaASRT2 and BnaCSRT2 expression, while JA treatment led to a time-dependent upregulation specifically in BnaASRT2 (Fig 8A). ACC treatment elicited a biphasic response in BnaASRT2 expression (initial upregulation followed by downregulation) with a turning point at 1 h; similarly, TZ treatment induced a biphasic expression trend in BnaASRT2 over time (S6 Table, S10 Fig).
In root tissues, IAA treatment induced a biphasic response in both BnaASRT2 and BnaCSRT2 expression, characterized by initial upregulation followed by downregulation with a turning point at 1 h (Fig 8B, S6 Table). Overall, BnaSRT2 displayed more pronounced expression alterations in response to hormone treatments, which was consistent with the distribution of hormone-responsive elements in its promoter region.
Discussion
As key epigenetic regulators, histone deacetylases (HDACs) modulate chromatin structure and gene expression via acetyl group removal from histone tails, thereby exerting pivotal functions in plant growth, development, and stress responses [50,51]. Within the HDAC superfamily, SIR2 family is defined by its conserved NAD ⁺ -dependent deacetylase activity. Having been extensively studied in yeast, animals, and model plants, This family is involved in a variety of biological processes, such as chromatin silencing, senescence, and stress tolerance. [52,53]. However, systematic identification and functional characterization of the SIR2 family in Brassica napus L. (rapeseed)-a globally important oilseed crop with high economic value-remain limited. In this study, we performed genome-wide identification of SIR2 family members in rapeseed, analyzed their phylogenetic relationships, structural features, and expression patterns, and provided insights into their potential roles in regulating the biological processes of rapeseed.
Genome-wide screening identified 4 SIR2 family members in rapeseed, which clustered into 2 distinct subgroups based on phylogenetic analysis (Fig 2B). This number exceeds that in Arabidopsis thaliana (2 members) [54] and Oryza sativa (rice, 2 members) [55], a discrepancy presumably linked to rapeseed’s allopolyploid genome (AACC). Rapeseed originated from the hybridization and subsequent polyploidization of Brassica rapa (AA genome) and Brassica oleracea (CC genome) [56]. Synteny analysis further confirmed that most BnaSRT members have orthologous counterparts in B. rapa and B. oleracea (Fig 2D), suggesting that polyploidization likely drove the expansion of the SIR2 family in rapeseed. Such expansion may have facilitated functional divergence, enabling rapeseed to adapt to more complex environmental niches during evolution.
All identified BnaSRTs proteins contain the conserved SIR2 domain (PF02146), which harbors the catalytic core essential for NAD ⁺ -dependent deacetylase activity [57] (Fig 3A). The high conservation of this domain across subgroups indicates that the core enzymatic functions of SIR2 proteins are evolutionarily conserved in rapeseed, consistent with observations in other plant species [54,55]. Despite strong conservation in their coding sequences, Bna**SRTs genes exhibit substantial diversity in their upstream promoter regions: cis-acting elements differ markedly not only between SRT1 and SRT2 subgroups but also between homologous genes derived from different subgenomes (Fig 3D). This divergence points to functional specialization and regulatory recombination during evolution, supporting the notion that the expansion of the rapeseed SIR2 family is accompanied by both functional conservation and diversification.
Furthermore, molecular docking and inhibitor-based assays of potential modification sites revealed that rapeseed BnaSRTs target multiple deacetylation sites (Fig 4C and 5A). In contrast, most reported SRT-mediated modifications in other species are restricted to H3K9ac [11,18]. This suggests potential functional diversification of rapeseed BnaSRTs. Given the current focus on H3K9ac with limited research on other modification sites, substantial scope remains for exploring BnaSRTs' functions. Beyond histone deacetylation, emerging evidence indicates that SIR2 family proteins possess broader deacylation activities. In rice, OsSRT1 and OsSRT2 mediate the removal of lysine 2-hydroxyisobutyrylation (Khib) marks to regulate resistance against the pathogen Ustilaginoidea virens [58]. They also function in removing β-hydroxybutyrylation (Kbhb) modifications to modulate immune responses [11], while OsSRT2 is involved in crotonylation (Kcr) removal [59]. In wheat, TaSRT1 negatively regulates cold tolerance by decrotonylating TaPGK, thereby inhibiting its accumulation [60]. In cotton, GhSRT1 suppresses Verticillium wilt resistance by removing Khib and succinylation (Ksuc) marks [61]. BnSRTs may similarly possess these biochemical activities and exert specific biological functions in rapeseed.
Gene expression patterns are tightly linked to biological function [62]. Tissue-specific expression analysis revealed distinct expression profiles of BnSRTs across roots, stems, leaves, buds, flowers, siliques, and seeds (Fig 6). For example, high expression of BnaASRT2 in siliques suggests potential involvement in seed development or nutrient accumulation-a critical agronomic trait in rapeseed. Its preferential expression in roots also implies roles in root system establishment, which is vital for nutrient uptake and stress tolerance [63]. These tissue-specific patterns share both similarities and differences with those in Arabidopsis [5], reflecting conserved and species-specific regulatory roles of the SIR2 family in plant development. Notably, BnSRTs members exhibit dynamic expression responses to abiotic stresses (e.g., drought, salt, cold) and hormone treatments (e.g., ABA, JA, GA) (S7 and S8 Figs). For instance, BnaASRT1 and BnaASRT2 are significantly upregulated in roots under drought stress; BnaASRT2 and BnaCSRT1 are induced by cold treatment in leaves and roots, respectively; and ABA-a key stress hormone-strongly activates BnaASRT2, pointing to potential involvement in ABA-mediated stress signaling. Furthermore, the analysis of the promoter elements of BnaSRTs also revealed the presence of multiple cis-acting elements that respond to abiotic stresses and hormones. These findings align with previous reports that plant SIR2 regulate stress responses through epigenetic modulation of stress-related gene expression [64,65]. For rapeseed, which frequently encounters environmental stresses during cultivation, these stress-responsive BnaSRTs represent promising candidates for improving stress tolerance via epigenetic engineering.
In yeast and animals, the SIR2 family is well-characterized for roles in lifespan regulation, DNA repair, and metabolic homeostasis [66]. In plants, accumulating evidence supports their involvement in transposon silencing, flowering time control, and environmental adaptation [54,67]. This study extends these insights to rapeseed, a polyploid crop with a complex genome. The expansion of the BnaSRT family and its subgroup-specific structural features suggest adaptation to regulating polyploidy-associated processes, such as genome dosage balance or subgenome expression dominance [68], expression under stress, indicating functional redundancy, while others show subgenome-specific expression, suggesting neofunctionalization. Compared with SIR2 families in diploid relatives (B. rapa and B. oleracea), the larger number of BnaSRT members and their diverse expression patterns highlight polyploidization' s role in shaping rapeseed' s epigenetic regulatory network. This expansion likely enhanced environmental adaptation flexibility-a trait critical for rapeseed' s domestication and widespread cultivation [69].
While this study provides a comprehensive framework for understanding the BnaSRTs family, limitations remain. First, functional inferences based on expression and structural features require experimental validation. Future studies employing CRISPR/Cas9-mediated knockout or overexpression of key BnaSRT members will clarify their roles in rapeseed development and stress responses. Second, the molecular mechanisms underlying BnaSRT-mediated regulation-including target genes, interacting partners (e.g., chromatin modifiers or transcription factors), and post-translational modifications-remain unclear. Chromatin immunoprecipitation sequencing (ChIP-seq) or yeast two-hybrid assays could identify downstream targets and protein interactors, elucidating their epigenetic regulatory networks. Third, responses of BnaSRTs to biotic stresses (e.g., pathogens) or agronomic cues (e.g., nutrient deficiency) were not investigated; given their potential role in integrating environmental signals, this warrants further research.
In summary, through genome-wide identification and expression profiling, this study reveals evolutionary expansion, structural conservation, and diverse developmental/stress-responsive expression patterns of the rapeseed SIR2 family. These findings advance our understanding of epigenetic regulation in polyploid crops and provide valuable candidate genes for improving rapeseed agronomic traits via epigenetic manipulation.
Conclusions
This study is the first systematic identification and characterization of the SRTs gene family in Brassica napus, identifying four members (BnaASRT1/2, BnaCSRT1/2) and analyzing their structural, evolutionary, and functional traits. Evolutionarily, rapeseed SRTs are highly conserved with homologs in Arabidopsis, rice, and humans (retaining the NAD ⁺ -dependent deacetylase-related SIR2 domain) and were stably inherited during polyploidization. Structurally, they share conserved gene architectures and motifs, while divergent promoter cis-elements (e.g., stress/hormone-responsive) imply SRT1/SRT2 subgroup functional specialization. Functionally, BnaSRTs mediate histone deacetylation (binding H3K9/14/23ac and H4K5ac strongly) and regulate photosynthesis, metabolism, and stress pathways (supported by inhibitor and transcriptome analyses). They show tissue-specific expression (e.g., seed maturation) and respond to abiotic stresses/hormones. Overall, the study reveals SRTs' evolutionary conservation and functional diversity in rapeseed. BnaSRTs (especially BnaSRT2s) are promising for epigenetic engineering to improve stress tolerance and agronomic traits, with future studies needed to clarify mechanisms for breeding application.
Supporting information
S1 Raw imagesUncropped, original and complete Western Blot images.(PDF)
S1 FigProtein sequence alignment with labeled SIR2 domain.(TIF)
S2 FigKEGG pathway enrichment analysis results of BnaSRT proteins.(TIF)
S3 Fig3D structure diagram of histone acetylation modifications.(TIF)
S4 Fig2D molecular docking diagram of BnaASRT1 and 5 types of histone acetylation modifications.(TIF)
S5 Fig2D molecular docking diagram of BnaCSRT1 and 5 types of histone acetylation modifications.(TIF)
S6 Fig2D molecular docking diagram of BnaASRT2 and 5 types of histone acetylation modifications.(TIF)
S7 Fig2D molecular docking diagram of BnaCSRT2 and 5 types of histone acetylation modifications.(TIF)
S8 FigWB full results.The red box marks the part that is cropped and used.(TIF)
S9 FigExpression profiles of the rapeseed BnaSRTs under different abiotic stresses.These expression profiles pertain respectively to the leaves (A) and the roots (B) at different time points under different abiotic stresses (freezing, heat and osmotic).(TIF)
S10 FigExpression profiles of the rapeseed BnaSRTs in response to different hormone treatment.These expression profiles pertain respectively to the leaves (A) and the roots (B) at different time points in response to different hormone treatment (ACC, TZ and BL).(TIF)
S1 TableThe accession number for the genes used in this study.(XLSX)
S2 TableDetailed information about motifs.(XLSX)
S3 TableDetailed information about GO and KEGG pathway enrichment analysis results.(XLSX)
S4 TableGene expression information in the tissues of rapeseed.(XLSX)
S5 TableGene expression information in rapeseed during abiotic stress treatment.(XLSX)
S6 TableGene expression information in rapeseed during hormone stress treatment.(XLSX)
S7 TableGO and KEGG enrichment analysis.(XLSX)
S8 TableNucleotide sequences of primers used in this study.(XLSX)
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Tahir MS, Tian L. HD 2-type histone deacetylases: unique regulators of plant development and stress responses. Plant Cell Rep. 2021;40(9):1603–15. doi: 10.1007/s 00299-021-02688-3 34041586 · doi ↗ · pubmed ↗
- 2Chang 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(5):563–80. doi: 10.1111/jipb.12901 31872527 · doi ↗ · pubmed ↗
- 3Ma L, Xing L, Li Z, Jiang D. Epigenetic control of plant abiotic stress responses. J Genet Genomics. 2025;52(2):129–44. doi: 10.1016/j.jgg.2024.09.008 39322116 · doi ↗ · pubmed ↗
- 4Rothbart SB, Strahl BD. Interpreting the language of histone and DNA modifications. Biochim Biophys Acta. 2014;1839(8):627–43. doi: 10.1016/j.bbagrm.2014.03.001 24631868 PMC 4099259 · doi ↗ · pubmed ↗
- 5Gray SG, Ekström TJ. The human histone deacetylase family. Exp Cell Res. 2001;262(2):75–83. doi: 10.1006/excr.2000.5080 11139331 · doi ↗ · pubmed ↗
- 6Kuo MH, Allis CD. Roles of histone acetyltransferases and deacetylases in gene regulation. Bioessays. 1998;20(8):615–26. doi: 10.1002/(SICI)1521-1878(199808)20:8<615::AID-BIES 4>3.0.CO;2-H 9780836 · doi ↗ · pubmed ↗
- 7Chen X, Ding AB, Zhong X. Functions and mechanisms of plant histone deacetylases. Sci China Life Sci. 2020;63(2):206–16. doi: 10.1007/s 11427-019-1587-x 31879846 · doi ↗ · pubmed ↗
- 8Li H, Chen Y, Dai Y, Yang L, Zhang S. Genome-wide identification and expression analysis of histone deacetylase and histone acetyltransferase genes in response to drought in poplars. BMC Genomics. 2024;25(1):657. doi: 10.1186/s 12864-024-10570-1 38956453 PMC 11218084 · doi ↗ · pubmed ↗