- Open Access
Genome-wide identification of Brassicaceae histone modification genes and their responses to abiotic stresses in allotetraploid rapeseed
BMC Plant Biology volume 23, Article number: 248 (2023)
Histone modification is an important epigenetic regulatory mechanism and essential for stress adaptation in plants. However, systematic analysis of histone modification genes (HMs) in Brassicaceae species is lacking, and their roles in response to abiotic stress have not yet been identified.
In this study, we identified 102 AtHMs, 280 BnaHMs, 251 BcHMs, 251 BjHMs, 144 BnHMs, 155 BoHMs, 137 BrHMs, 122 CrHMs, and 356 CsHMs in nine Brassicaceae species, respectively. Their chromosomal locations, protein/gene structures, phylogenetic trees, and syntenies were determined. Specific domains were identified in several Brassicaceae HMs, indicating an association with diverse functions. Syntenic analysis showed that the expansion of Brassicaceae HMs may be due to segmental and whole-genome duplications. Nine key BnaHMs in allotetraploid rapeseed may be responsible for ammonium, salt, boron, cadmium, nitrate, and potassium stress based on co-expression network analysis. According to weighted gene co-expression network analysis (WGCNA), 12 BnaHMs were associated with stress adaptation. Among the above genes, BnaPRMT11 simultaneously responded to four different stresses based on differential expression analysis, while BnaSDG46, BnaHDT10, and BnaHDA1 participated in five stresses. BnaSDG46 was also involved in four different stresses based on WGCNA, while BnaSDG10 and BnaJMJ58 were differentially expressed in response to six different stresses. In summary, six candidate genes for stress resistance (BnaPRMT11, BnaSDG46, BnaSDG10, BnaJMJ58, BnaHDT10, and BnaHDA1) were identified.
Taken together, these findings help clarify the biological roles of Brassicaceae HMs. The identified candidate genes provide an important reference for the potential development of stress-tolerant oilseed plants.
Histone modification (HM) is an epigenetic regulatory mechanism that plays crucial roles in various aspects of plant growth and stress response by activating or silencing gene expression [1,2,3,4]. HM genes (HMs) include histone methyltransferases (HMTs), histone demethylases (HDMs), histone acetylases (HATs), and histone deacetylases (HDACs) [5,6,7,8].
HMTs are encoded by the SET-domain group (SDG) and protein arginine methyltransferase (PRMT) genes and catalyze HM . Methylation and environmental factors are related to stress, which affects gene expression by changing methylation levels and stress resistance . Several processes, such as fungal pathogen resistance, shoot and root branching, circadian cycle, hormone regulation, abscisic acid (ABA) morphogenesis, and salt stress, are affected by HMTs [11, 12]. For example, AtSDG8 is involved in the regulation of shoot meristem activity, while AtSDG26 and AtPRMT10 are involved in the regulation of flowering in Arabidopsis [13,14,15]. Histone modification can be erased by HDMs, including lysine-specific demethylase 1 (LSD1) and Jumonji C (JmjC) domain containing proteins [16,17,18]. HDMs function in brassinosteroid (BR) signaling, pollen development, chromatin regulation, floral induction, and floral organ formation [19, 20]. In Arabidopsis, JMJ30 expression changes in response to environmental stimuli, e.g., enhancement by salt and heat stress [21, 22], and flower repressor JMJ13 can be affected by temperature and photoperiod . HATs and HDACs catalyze the transfer of acetyl groups from acetyl-CoA to lysine residues [24, 25]. HATs and HDACs participate in the regulation of developmental transition, environmental signal responses, reproductive development, and gene silencing [26,27,28]. For example, HAC1 inactivation affects both vegetative and reproductive development in Arabidopsis , AtSRT2 regulates salt tolerance during seed germination , and AtHDT4 participates in abiotic stress responses . Previous studies have suggested that modifications affect functions, including transcriptional regulation of other genes in yeast .
Abiotic stress, such salinity, inappropriate nutrition, and metal toxicity, can adversely affect crop growth and yield [33, 34]. Nutrient imbalances, membrane damage, and dysfunctional antioxidant system can occur under soil salinization . Various nutrients are essential for optimal plant growth and yield. Nitrogen (N) is a macronutrient “life element” that strongly affects plant growth and development , while excess ammonium (NH4+), an inorganic N nutrient, is toxic to plants . Phosphorus (P) shares essential roles in regulating plant energy metabolism, and its deficiency can reduce cell division and elongation in grass leaves . Potassium (K) is a vital macronutrient for plant growth and organ development, and participates in many physiological processes, such as osmoregulation. Moreover, K + transport participates in abiotic stress responses [39, 40]. Boron (B) is a micronutrient essential for the transport of carbohydrates, although both excess and deficiency can adversely impact crop growth and yield [41, 42]. Plants can also be affected by non-essential heavy metals, such as cadmium (Cd), which is highly biotoxic and easily absorbed by plants through sewage effluent, industrial waste, and agricultural run-off .
Brassicaceae plants are important and economically valuable crops, noted for their oil production [44, 45]. Given their immobility, plants are unable to avoid abiotic and biotic stresses, which can impair growth, development, and production. However, plants can adapt to stress by activating a series of physiological and molecular mechanisms, such as HM [46, 47]. Therefore, improving stress resistance and yield in Brassicaceae plants is a key goal of breeding . To date, however, few studies have explored the regulation of gene expression related to stress resistance or conducted systematic study of HMs in Brassicaceae species. Here, we conducted a comprehensive study of HMs in nine Brassicaceae species, including Arabidopsis thaliana, Brassica napus, Brassica carinata, Brassica juncea, Brassica nigra, Brassica oleracea, Brassica rapa, Capsella rubella, and Camelina sativa. We further determined their chromosomal locations, conserved domains, gene structures, phylogenetic relationships, and syntenies. The responses of HMs to NH4+ toxicity, B deficiency and excess, Cd exposure, K shortage, N limitation, P starvation, and salt stress were explored in allotetraploid rapeseed. Potential candidate BnaHMs that responded to the above stresses were also identified. This study provides important clues for understanding the Brassicaceae HM gene family.
Genome-wide identification of Brassicaceae HMs and their phylogenetic analysis
In the present study, we identified 1 798 HMs, including 102, 280, 251, 251, 144, 155, 137, 122, and 356 in (A) thaliana, (B) napus, B. carinata, B. juncea, B. nigra, B. oleracea, B. rapa, (C) rubella, and C. sativa (Figure S1 and Table S1). The number of HMTs, HDMs, HATs, and HDACs varied among species, with 2.7-, 2.5-, 2.5-, and 3.5-fold as many BnaHMs, BcaHMs, BjuHMs, and CsHMs as AtHMs, respectively (Figure S1a). There were 47–159 SDGs, 7–27 PRMTs, 2–8 HDMAs, 20–77 JMJs, 3–10 HAGs, 1–7 HAMs, 4–10 HACs, 1–4 HAFs, 12–40 HDAs, 2–8 SRTs, and 4–16 HDTs in the above Brassicaceae species, respectively (Figure S1b), named according on their chromosomal position in each species (Figure S2).
To elucidate the evolutionary relationships among HMs, unrooted phylogenetic trees were constructed. Generally, each type of HAT, HDAC, HDM, and HMT shared relatively close relationships in distinct groups, with some exceptions (Figure S3). For example, in terms of HATs, all HACs were in group a, most HAGs were in group b, and HAFs and HAMs were in groups c and d, respectively (Figure S3-1).
Conserved domain, structure, and synteny analysis of HMs
Diverse conserved domains were identified in the different HMs (Figure S4) and the number of conserved motifs was determined in the Arabidopsis HMs (Figure S4-1). Most conserved domains in the Arabidopsis HMs were also present in the non-model plants (B. napus, B. carinata, B. juncea, B. nigra, B. oleracea, B. rapa, C. rubella, and C. sativa). However, several distinct domains were identified in the non-model Brassicaceae HMs (Figure S4), including the SHOCT domain in BjHDA10 and BjHDT11, which may bind to itself to perform important functions as an oligomerization domain or bind to other protein domains/motifs and nucleic acids . BjHDA24 shares a domain with the CYCLIN superfamily, which functions in the cell cycle and transcriptional control (Figure S4-8). In general, each class of HM shared a similar gene structure. Of note, several HMs, including BoJMJ23, BoSDG24, and BcJMJ54, contained long introns (Figure S4).
To determine the expansion patterns of HMs, duplication events within gene pairs were investigated in duplicated blocks of each Brassicaceae genome. In total, 1 176 gene pairs were identified, including 11, 256, 194, 215, 49, 55, 42, 15, and 339 pairs in (A) thaliana, (B) napus, B. carinata, B. juncea, B. nigra, B. oleracea, B. rapa, (C) rubella, and C. sativa, respectively (Figure S5 and Table S2). To understand the potential roles of unknown Brassicaceae HM genes, collinearity analysis was performed between Arabidopsis and non-model Brassicaceae species. In total, 151, 157, 178, 101, 89, 216, 83, and 109 gene pairs were identified in A. thaliana-B. carinata, A. thaliana-B. juncea, A. thaliana-B. napus, A. thaliana-B. oleracea, A. thaliana-B. rapa, A. thaliana-C. sativa, A. thaliana-C. rubella, and A. thaliana-B. nigra, respectively (Figure S6 and Table S3).
Effects of NH4+, salt, B, and Cd on expression patterns of BnaHMs
Although NH4+ is the main N source for plants, excess can cause toxicity to crops and reduce grain yields [50, 51]. Here, the expression profiles of BnaHM genes were investigated to predict their potential involvement in NH4+ toxicity resistance. In roots, 12 BnaHMs were differentially expressed after excess NH4+ treatment, half of which were up-regulated (Fig. 1a). In shoots, 37 BnaHMs were differentially expressed, six of which showed low levels in the NH4+-treated group (Fig. 1b). Among these differentially expressed genes (DEGs), based on gene co-expression network analysis (GCNA), BnaPRMT11 and HDT10 may be critical genes in response to NH4+ toxicity (Fig. 1c). In roots, 15 and 38 BnaHMs were suppressed and induced by salt treatment, respectively (Fig. 1d), with BnaHDA11 and BnaPRMT8 potentially playing roles in salt adaptation (Fig. 1e). In shoots, 48 BnaHMs, especially BnaSDG58, were markedly regulated by salt exposure (Fig. 1f). According to GCNA, BnaHDT10 was identified as a hub gene in response to salt stress (Fig. 1g).
Both B deficiency and toxicity can have adverse effects on plant growth and development . However, whether BnaHMs are involved in B-mediated plant growth is unclear. Our results identified several BnaHMs that were differentially expressed after B treatment (Fig. 2). In roots, BnaHDA3 and BnaSDG46 were inhibited by B deficiency, while five BnaHMs were induced (Fig. 2a). B toxicity also altered the expression patterns of BnaHMs (Fig. 2b, e). In the B deficiency group, BnaJMJ18, BnaSDG82, and BnaJMJ9 were up-regulated in shoots, while 69 BnaHMs were down-regulated (Fig. 2c). BnaSDG4 was identified as a key gene (Fig. 2d). In shoots, only BnaHDA12 increased in response to excess B, while the remaining BnaHMs were reduced (Fig. 2e). Among them, BnaSDG94 was identified as a potential hub gene (Fig. 2f).
Cd is a non-essential heavy metal toxic for plant growth . In roots, 15 and six BnaHMs exhibited higher and lower expression, respectively, in the Cd-treated group compared with the control group (Fig. 3a). In shoots, BnaSDG30 and BnaSDG75 were significantly inhibited by Cd, while BnaHDT2 was induced (Fig. 3d). BnaSDG75 was also identified as a key gene in the co-expression network (Fig. 3e).
Effects of N, K, and P on expression patterns of BnaHMs
As an essential macronutrient, N is required for rapeseed growth and development . To investigate the response of BnaHMs to N limitation, we identified their expression profiles. BnaSDG4, BnaJMJ9, and BnaJMJ43 were up-regulated in the N-treated roots, while nine other genes were down-regulated (Fig. 3b). In shoots, BnaPRMT10, BnaHAF1, BnaHDA27, BnaHDA11, and BnaSDG23 were substantially induced by N deficiency, while BnaSDG43, BnaJMJ13, BnaSDG102, and BnaJMJ61 were repressed (Fig. 3c).
Previous studies have shown that K can also cause stress to plants [55, 56]. Our results showed that limited K induced 11 BnaHMs and suppressed seven BnaHMs in the roots, especially BnaSDG81 (Fig. 4a). In shoots, 10 BnaHMs (e.g., BnaHDA15, BnaSDG46, and BnaSDG1) were decreased after K treatment, while 52 BnaHMs, especially BnaJMJ47, BnaSDG86, and BnaSDG88, were increased (Fig. 4d). BnaHDA15 was identified as a key gene according to GCNA (Fig. 4e). Given its close involvement in photosynthesis, P is an essential nutrient for plant growth and development . Here, in response to P stress, the expression levels of several BnaHMs, especially BnaJMJ6, increased in roots, whereas five BnaHMs were markedly suppressed (Fig. 4b). In shoots, 14 BnaHMs showed higher expression levels after P treatment, while 29 were inhibited by P stress (Fig. 4c).
Identification of weighted gene co-expression network analysis (WGCNA) modules and hub genes associated with target traits
All genes in the RNA sequencing (RNA-seq) data, not just DEGs, were analyzed for significant associations with phenotypes using WGCNA based on previous methods . WGCNA was established to analyze hub genes in response to A, salt, Cd, N, and K stress.
The “lightyellow” (r = -0.64, p < 0.01) and “turquoise” (r = -0.92, p < 0.01) modules were negatively correlated with chlorophyll content (SPAD) after NH4+ toxicity treatment (Fig. 5a). Two co-expression networks were constructed to identify core genes. In the “lightyellow” and “turquoise” modules, BnaPRMT15, and BnaSDG64, BnaSDG53, and BnaSDG36 were respectively identified in response to NH4+ exposure (Fig. 5b, c). In total, 37 genes in the “lightyellow” module and 40 genes in the “turquoise” module are involved in various stresses, such as oxidative stress, and interact with core genes (Table S5-1).
WGCNA was also performed to evaluate the relationship between modules and salinity (Fig. 6). The “salmon” (r = -0.83, p < 0.05) and “blue” (r = -0.91, p < 0.05) modules were negatively correlated with biomass and leaf area, respectively (Fig. 6a). Four BnaHMs (BnaSDG53, BnaSDG36, BnaSDG46, and BnaSDG64) were identified as important genes in the “blue” module (Fig. 6b). BnaHAG3, BnaHDA12, BnaHDA8, and BnaHAG7 were identified as hub genes in the “salmon” module (Fig. 6c). In addition, seven and nine genes in the “salmon” and “blue” modules, respectively, were salt-responsive and associated with core genes (Table S5-2).
Using WGCNA, core genes associated with Cd stress were identified. As shown in Fig. 7, both “green” and “purple” module were negatively correlated with SPAD and positively correlated with biomass (Fig. 7a). “Yellow” module was too, while “dark turquoise” was negatively correlated with SPAD and “purple” was positively correlated with biomass (Fig. 7a). Gene interaction networks were established for these two modules, and two key genes were identified (BnaSDG46 and BnaPRMT4, respectively) (Fig. 7b, c). In both modules, several Cd-resistance genes were identified and were associated with core genes (Table S5-3).
The relationship between WGCNA modules and N shortage was also explored. All genes were clustered into seven modules, and genes in the “green” module (r = -0.83, p < 0.05) were significantly correlated with SPAD (Fig. 8a). Three hub genes (BnaSDG53, BnaHDA1, and BnaSDG46) were screened from co-expression gene network mapping (Fig. 8b) and may play additional roles in adaptation to various stresses. Furthermore, several genes in the “green” module play roles in stress adaptation and interact with the three hub genes (Table S5-4).
In response to K stress, eight WGCNA modules were obtained. The “turquoise” module (r = -0.92, p < 0.05) showed a negative correlation with SPAD (Fig. 9a). BnaSDG60 and BnaSDG46 were identified as critical genes in this module (Fig. 9b). In addition, stress-related and K-transport genes in the ‘turquoise” module were associated with BnaSDG60 and BnaSDG46 (Table S5-5).
Diverse responses of BnaHMs to nutrient stresses
To investigate whether BnaHMs responded to diverse stresses simultaneously, we constructed a Venn diagram. Results showed that most BnaHMs were affected by more than one stress (Fig. 10 and Table S4). For example, 27 BnaHMs were simultaneously under the control of two stresses; 31 BnaHMs simultaneously responded to three stress signals; 32 BnaHMs simultaneously responded to four stresses; 11 BnaHMs were controlled by five stresses; and two genes responded to six stresses.
HMs play essential roles in plant growth and stress responses and have been successfully identified in many plants, such as Arabidopsis, wheat, and maize . However, information on Brassicaceae HMs remains limited. In this study, we systematically characterized HMs in nine Brassicaceae species and identified 1 798 HMs, including 102 AtHMs, 280 BnaHMs, 251 BcHMs, 251 BjHMs, 144 BnHMs, 155 BoHMs, 137 BrHMs, 122 CrHMs, and 356 CsHMs. We further analyzed their phylogeny, conserved domains, gene structure, and synteny, as well as their expression profiles in response to NH4+, B, salt, Cd, N, P, and K stress. These results will contribute to a comprehensive understanding of Brassicaceae HM genes.
Comparison of HMs among nine Brassicaceae species
We identified 280, 251, 251, 144, 155, 137, 122, and 356 HMs in B. napus, B. carinata, B. juncea, B. nigra, B. oleracea, B. rapa, Capsella rubella, and Camelina sativa, respectively (Figure S1 and Table S1). We also found significantly more BnaHMs, BcHMs, BjHMs, BnHMs, BoHMs, BrHMs, CrHMs, and CsHMs than AtHMs (2.7-, 2.4-, 2.4-, 1.4-, 1.5-,1.3-, 1.1-, and 3.4-fold higher, respectively) (Figure S1 and Table S1). Orthologous HMs were found based on synteny analysis. We identified 11 AtHM, 256 BnaHM, 194 BcHM, 215 BjHM, 49 BnHM, 55 BoHM, 42 BrHM,15 CrHM, and 339 CsHM pairs (Figure S5 and Table S2). Results showed that more segmental duplications of HMs were found in non-model Brassicaceae species than in Arabidopsis, which may induce the expression of non-model Brassicaceae HMs. Whole-genome replication is known to occur in B. napus, B. carinata, B. juncea, and Camelina sativa [60,61,62,63,64]. Therefore, segmental and whole-genome duplications may have contributed to the expansion and evolution of HMs in the above species.
Synteny analysis between duplicated blocks of Arabidopsis-B. carinata, Arabidopsis-B. juncea, Arabidopsis-B. napus, Arabidopsis-B. oleracea, Arabidopsis-B. rapa, Arabidopsis-C.sativa, Arabidopsis-C. rubella, and Arabidopsis-B. nigra was also performed, yielding 151, 157, 178, 101, 89, 216, 83 and 109 gene pairs, respectively (Figure S6 and Table S3). These gene pairs are considered to have originated from common ancestors with AtHMs [60,61,62,63,64], suggesting that they may have similar functions to the corresponding Arabidopsis genes. Thus, the functions of non-model Brassicaceae HMs were predicted based on homologous Arabidopsis HMs. Several AtHMs are involved in stress responses [10, 65, 66]. Although many unknown non-model Brassicaceae HMs could be inferred from orthologous Arabidopsis genes, these comparisons require further experiments.
Conserved domains are associated with gene function . We identified typical domains in the HMs (Figure S4). Most Brassicaceae HMs with conserved domains shared similar functions, but several distinct domains were identified in several non-model Brassicaceae HMs, such as the FYVE_like_SF superfamily domain in BnaJMJ65, which plays an important role in vesicular traffic and signal transduction (Figure S4-2). Novel functions may be predicted from unique domains, and thus greater attention should be paid to genes with special elements in the future.
Putative functions of BnaHMs in stress response
HMs are important in plant defense. Here, the expression patterns of BnaHMs were determined to explore their function under various stresses. In roots and shoots, 79 and 81 BnaHMs were up-regulated or down-regulated by B deficiency and toxicity (Fig. 2) and BnaHM expression patterns were changed by NH4+ and N deficiency (Figs. 1a-c and 3b-c). More than 50 BnaHMs showed differential expression in response to P shortage, and many BnaHMs were influenced by K deficiency stress (Fig. 4). These findings indicate that BnaHMs play essential roles in the stress response. Various abiotic stresses, including drought, salinity, and cold, adversely affect plant growth and development. HMs share important roles in regulating stress adaptation. For example, AtHDA6 and AtHDA19 are involved in ABA responses and are required for salt tolerance [59, 68]. Here, many BnaHMs responded to Cd and salt stress, with altered expression in the roots or shoots (Figs. 1d-g and 3d-e). These findings suggest that BnaHMs and methylation play essential roles in rapeseed resistance to diverse stresses.
Candidate BnaHMs were determined through DEG co-expression analysis and WGCNA. BnaPRMT15, BnaSDG36, BnaSDG53, BnaSDG64, BnaHDT10, and BnaPRMT11 were identified in response to NH4+ toxicity (Figs. 1c and 5). BnaHDA11, BnaHDT10, BnaPRMT8, BnaHAG3, BnaHAG7, BnaSDG36, BnaSDG46, BnaSDG53, BnaHDA12, BnaHDA8, and BnaSDG64 were associated with plant survival under salt stress (Fig. 1e and g, and Fig. 6). BnaPRMT4, BnaSDG46, and BnaSDG75 were identified as Cd-related genes (Figs. 3e and 7). BnaSDG4 and BnaSDG94 were identified as B stress candidate genes (Fig. 2d and f). BnaSDG46, BnaSDG53, and BnaHDA1 were identified as N-deficiency candidate genes (Fig. 8). BnaHDA15, BnaSDG46, and BnaSDG60 were identified as K limitation-related genes (Fig. 9). Based on orthologous gene analysis, the ortholog of AtHDA6, which responds to drought stress , was identified as BnaHDA8, and the ortholog of AtHDA14, which functions in regulating stress responses [70,71,72,73,74], was identified as BnaHDA1. In addition, according to WGCNA, several downstream genes identified in modules that may be involved in various stresses, such as low temperature and salt, interacted with the core genes, indicating that these core genes may participate in stress tolerance by interacting with downstream stress-related genes (Table S5). These results suggest that HMs play an important role in stress response. As such, future studies should pay attention to the above candidate genes.
This study also found that many differentially expressed BnaHMs responded to different stresses at the same time (Table S4). For example, two BnaHMs (BnaSDG10 and BnaJMJ58) were simultaneously regulated by six stresses, and 11 BnaHMs (e.g., BnaHDT10, BnaSDG46, BnaPRMT10) were simultaneously regulated by five stresses. However, certain genes were only impacted by a single stress signal, implying that many BnaHMs may participant in different stresses, while others only play a core role under a specific stress.
Previous studies have shown that several HMs in rice may participate in stress adaptations. For example, OsHDT701 and OsHDT702 in rice are repressed by drought and salt simultaneously [75, 76]. Here, several key genes identified by co-expression analysis or WGCNA also responded to more than three different types of stress. For instance, BnaPRMT11and BnaHDA1 were differentially expressed under four and five types of stress, respectively (Fig. 10 and Table S4) and BnaSDG10 and BnaJMJ58 simultaneously responded to six different stresses. The salt stress-correlated core gene BnaHDT10 also responded to four other stresses (i.e., A, B, Cd, and P). In addition, BnaSDG46 was identified as a salt-, B-, Cd- and K-related key gene by WGCNA. These results suggest that the above hub BnaHMs may play critical roles in resistance to multiple stressors, and that they may show different functions under different stress. Therefore, future studies should focus on the potential functions of these genes.
HM gene identification, phylogenetic relationship, chromosomal location, conserved domains, gene structure, and synteny
Known AtHM protein sequences were used as queries and the B. napus, B. carinata, B. juncea, B. nigra, B. oleracea, B. rapa, C. rubella, and C. sativa protein databases were searched using “Blast Several Sequences to a Big Database” in TBtools  with an e-value of e-5. After aligning the full-length protein sequences by ClustalW with default parameters, MEGA X was used to construct the phylogenetic tree with the maximum-likelihood method .
Using chromosome length and gene position files, the chromosomal distributions of HMs were acquired and visualized using “Gene Location Visualize (Advanced)” in TBtools. The conserved domains in HMs were confirmed using the Batch Web CD-Search Tool (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi) . The conserved domains were visualized using “Visualize NCBI CDD Domain Pattern” in TBtools . The Visualize Gene Structure (Basic) tool was used to draw the gene structure map based on generic feature format v3 (gff3) files of the HMs.
We used “One Step MCScanX” in TBtools to analyze HM duplication events with genome sequences and gff3 files. “Table Row Extract or Filter”, “File Transformat for Microsynteny Viewer and Advanced Circos”, “Fasta stats”, and “File Merge for MCScanX” in TBtools were used to visualize the syntenic relationships of HM genes based on previous studies .
Transcriptome analysis, GCNA, and WGCNA of BnaHMs
The transcriptome data can be found in previously published papers [79,80,81,82,83]. All data required to reproduce these findings can be obtained by contacting the correlation authors. Fastp software (v0.20.1) was used to evaluate the overall sequencing quality of the raw reads and low-quality reads were removed. Alignment of high-quality reads with B. napus reference genome sequences (http://cbi.hzau.edu.cn/cgi-bin/rape/download_ext, accessed on 15 May 2022) was performed using Hisat2 (v2.1.0) and SAMtools (v1.6) software. Stringtie (v1.3.3b) was used to calculate the expression levels of high-confidence genes in each sample. The R package “edgeR”, with p < 0.05, false-discovery rate (FDR) < 0.05, and |log2(fold-change)| ≥ 1, was used to define DEGs. GCNA was performed using the cor.test function in R (v4.1), and network visualized using Cytoscape (v3.8.2, https://cytoscape.org/download.html, accessed on 13 April 2022) . The R WGCNA package (v1.51) was used to complete WGCNA with high-quality genes. Significant module-trait relationships with target traits were determined by calculating modular trait gene values. Gene co-expression network maps were generated using Cytoscape (v3.8.2, https://cytoscape.org/download.html, accessed on 13 April 2022). The gene with high | log2 (a fold - change) | and degree are selected as hub gene, and was placed in the middle of the network.
Plant materials and treatments
Uniform 7-day-old B. napus (Zhongshuang 11) seedings were transplanted into black plastic containers containing Hoagland nutrient solution (5.0 mM KNO3, 1.0 mM KH2PO4, 2.0 mM MgSO4·7H2O, 5.0 mM Ca (NO3)2·4H2O, 0.10 µM Na2MoO4·2H2O, 0.050 mM EDTA-Fe, 0.80 µM ZnSO4·7H2O, 9.0 µM MnCl2·4H2O, 0.30 µM CuSO4·5H2O, and 46 µM H3BO3). Before treatments, the B. napus seedlings were cultivated for 10 days (d) in a chamber under 25 °C daytime/22°C night-time temperature, 300–320 µmol m− 2 s− 1 light intensity, 16-h light/8-h dark photoperiod, and 70% relative humidity. B deficiency and toxicity treatments: We cultivated 17-day-old seedlings for 10 d in B-deficient (0.25 µM H3BO3) and B-excess (1 500 µM H3BO3) treatment groups; N, P, and K depletion treatments: We cultivated 17-day-old B. napus seedlings in Hoagland nutrient solution (consisting of 0.30 mM N, 5 mΜ P, and 0.30 mM K) for 3 d; NH4+toxicity treatment: We cultivated 17-day-old uniform Zhongshuang 11 seedlings in Hoagland nutrient solution (consisting of normal nitrate) for 10 d, followed by transfer to a N-free solution for 3 d, and final exposure to 9.0 mM NH4+ (excess NH4+) for 6 h; Cd toxicity and salt treatments: For Cd- and salt-treatment, we cultivated 17-day-old Zhongshuang 11 seedlings in 10 µM CdCl2 and 200 mM NaCl for 12 h and 1 d, respectively. The seedlings in the control groups were cultivated in a normal solution for the appropriate times based on the aforementioned treatments. Transcriptome sequencing was performed using roots and shoots from control and stress-treated plants as described above [84,85,86].
In this study, 1 798 HM genes were systematically identified in nine Brassicaceae species. Their chromosomal locations, protein/gene structure, and phylogenetic and syntenic relationships were characterized. The BnaHMs responding to A, salt, Cd, N, and K stress were investigated through differential expression analysis (GCNA and WGCNA). Taken together, BnaPRMT11, BnaJMJ58, BnaSDG46, BnaHDA1, BnaSDG10, and BnaHDT10, were identified as potential hub genes, especially BnaSDG46 and BnaHDT10. Our results suggest that BnaHMs may be crucial for regulating stress adaptation in rapeseed. The candidate genes identified here should be validated in future studies.
The raw transcriptome sequencing data were submitted to the National Centre for Biotechnology Information (NCBI) (http://www.ncbi.nlm.nih.gov/) under BioProject PRJNA340053, PRJNA718104, and PRJCA001323. The datasets used and/or analyzed in the current study are available from the corresponding author upon reasonable request.
- HMs :
Histone modification genes
- NH4 + :
Differentially expressed genes
Weighted gene co-expression network analysis
Gene co-expression network analysis.
Klose RJ, Zhang Y. Regulation of histone methylation by demethylimination and demethylation. NAT REV MOL CELL BIO. 2007;8(4):307–18.
Fan S, Liu H, Liu J, Hua W, Xu S, Li J. Systematic analysis of the DNA methylase and demethylase gene families in rapeseed (Brassica napus L.) and their expression variations after salt and heat stresses. Int J Mol Sci. 2020;21(3):953.
Wang Z, Cao H, Chen F, Liu Y. The roles of histone acetylation in seed performance and plant development. PLANT PHYSIOL BIOCH. 2014;84:125–33.
An W. Histone acetylation and methylation: combinatorial players for transcriptional regulation. Subcell Biochem. 2007;41:351–69.
Xu J, Xu H, Liu Y, Wang X, Xu Q, Deng X. Genome-wide identification of sweet orange (Citrus sinensis) histone modification gene families and their expression analysis during the fruit development and fruit-blue mold infection process. FRONT PLANT SCI. 2015;6:607.
Li S, He X, Gao Y, Zhou C, Chiang VL, Li W. Histone Acetylation Changes in Plant Response to Drought Stress. GENES-BASEL 2021, 12(9).
Peng M, Ying P, Liu X, Li C, Xia R, Li J, Zhao M. Genome-wide identification of histone modifiers and their expression patterns during Fruit Abscission in Litchi. FRONT PLANT SCI. 2017;8:639.
Aiese CR, Sanseverino W, Cremona G, Ercolano MR, Conicella C, Consiglio FM. Genome-wide analysis of histone modifiers in tomato: gaining an insight into their developmental roles. BMC Genomics. 2013;14:57.
Chen DH, Qiu HL, Huang Y, Zhang L, Si JP. Genome-wide identification and expression profiling of SET DOMAIN GROUP family in Dendrobium catenatum. BMC PLANT BIOL. 2020;20(1):40.
Chinnusamy V, Zhu JK. Epigenetic regulation of stress responses in plants. CURR OPIN PLANT BIOL. 2009;12(2):133–9.
Ahmad A, Cao X. Plant PRMTs broaden the scope of arginine methylation. J GENET GENOMICS. 2012;39(5):195–208.
Dong G, Ma DP, Li J. The histone methyltransferase SDG8 regulates shoot branching in Arabidopsis. BIOCHEM BIOPH RES CO. 2008;373(4):659–64.
Cazzonelli CI, Cuttriss AJ, Cossetto SB, Pye W, Crisp P, Whelan J, Finnegan EJ, Turnbull C, Pogson BJ. Regulation of carotenoid composition and shoot branching in Arabidopsis by a chromatin modifying histone methyltransferase, SDG8. PLANT CELL. 2009;21(1):39–53.
Qi PL, Zhou HR, Zhao QQ, Feng C, Ning YQ, Su YN, Cai XW, Yuan DY, Zhang ZC, Su XM, et al. Characterization of an autonomous pathway complex that promotes flowering in Arabidopsis. NUCLEIC ACIDS RES. 2022;50(13):7380–95.
Niu L, Lu F, Zhao T, Liu C, Cao X. The enzymatic activity of Arabidopsis protein arginine methyltransferase 10 is essential for flowering time regulation. PROTEIN CELL. 2012;3(6):450–9.
Liu C, Lu F, Cui X, Cao X. Histone methylation in higher plants. ANNU REV PLANT BIOL. 2010;61:395–420.
He K, Cao X, Deng X. Histone methylation in epigenetic regulation and temperature responses. CURR OPIN PLANT BIOL. 2021;61:102001.
Jiang D, Yang W, He Y, Amasino RM. Arabidopsis relatives of the human lysine-specific Demethylase1 repress the expression of FWA and FLOWERING LOCUS C and thus promote the floral transition. PLANT CELL. 2007;19(10):2975–87.
Chen X, Hu Y, Zhou DX. Epigenetic gene regulation by plant Jumonji group of histone demethylase. Biochim Biophys Acta. 2011;1809(8):421–6.
Xing L, Qi S, Zhou H, Zhang W, Zhang C, Ma W, Zhang Q, Shah K, Han M, Zhao J. Epigenomic Regulatory mechanism in vegetative phase transition of Malus hupehensis. J AGR FOOD CHEM. 2020;68(17):4812–29.
Yamaguchi N, Matsubara S, Yoshimizu K, Seki M, Hamada K, Kamitani M, Kurita Y, Nomura Y, Nagashima K, Inagaki S, et al. H3K27me3 demethylases alter HSP22 and HSP17.6 C expression in response to recurring heat in Arabidopsis. NAT COMMUN. 2021;12(1):3480.
Yamaguchi N. Removal of H3K27me3 by JMJ Proteins Controls Plant Development and environmental responses in Arabidopsis. FRONT PLANT SCI. 2021;12:687416.
He KX, Mei HL, Zhu JP, Qiu Q, Cao XF, Deng X. The histone H3K27 demethylase REF6/JMJ12 promotes thermomorphogenesis in Arabidopsis. NATL SCI REV 2022, 9(5).
Boycheva I, Vassileva V, Iantcheva A. Histone acetyltransferases in plant development and plasticity. CURR GENOMICS. 2014;15(1):28–37.
Chen X, Ding AB, Zhong X. Functions and mechanisms of plant histone deacetylases. SCI CHINA LIFE SCI. 2020;63(2):206–16.
Zheng Y, Ding Y, Sun X, Xie S, Wang D, Liu X, Su L, Wei W, Pan L, Zhou DX. Histone deacetylase HDA9 negatively regulates salt and drought stress responsiveness in Arabidopsis. J EXP BOT. 2016;67(6):1703–13.
Hou Y, Lu Q, Su J, Jin X, Jia C, An L, Tian Y, Song Y. Genome-Wide Analysis of the HDAC Gene Family and Its Functional Characterization at Low Temperatures in Tartary Buckwheat (Fagopyrum tataricum). INT J MOL SCI 2022, 23(14).
Xing G, Jin M, Qu R, Zhang J, Han Y, Han Y, Wang X, Li X, Ma F, Zhao X. Genome-wide investigation of histone acetyltransferase gene family and its responses to biotic and abiotic stress in foxtail millet (Setaria italica [L.] P. Beauv). BMC PLANT BIOL. 2022;22(1):292.
Longo C, Lepri A, Paciolla A, Messore A, De Vita D, di Patti M, Amadei M, Madia VN, Ialongo D, Di Santo R et al. New Inhibitors of the Human p300/CBP Acetyltransferase Are Selectively Active against the Arabidopsis HAC Proteins. INTERNATIONAL JOURNAL OF MOLECULAR SCIENCES 2022, 23(18).
Tang WS, Zhong L, Ding QQ, Dou YN, Li WW, Xu ZS, Zhou YB, Chen J, Chen M, Ma YZ. Histone deacetylase AtSRT2 regulates salt tolerance during seed germination via repression of vesicle-associated membrane protein 714 (VAMP714) in Arabidopsis. NEW PHYTOL. 2022;234(4):1278–93.
Han ZF, Yu HM, Zhao Z, Hunter D, Luo XJ, Duan J, Tian LN. AtHD2D Gene Plays a Role in Plant Growth, Development, and Response to Abiotic Stresses in Arabidopsis thaliana. FRONT PLANT SCI 2016, 7.
Zhang XY, Clarenz O, Cokus S, Bernatavichute YV, Pellegrini M, Goodrich J, Jacobsen SE. Whole-genome analysis of histone H3 lysine 27 trimethylation in Arabidopsis. PLOS BIOL. 2007;5(5):1026–35.
Kamal KY, Khodaeiaminjan M, Yahya G, El-Tantawy AA, Abdel ED, El-Esawi MA, Abd-Elaziz M, Nassrallah AA. Modulation of cell cycle progression and chromatin dynamic as tolerance mechanisms to salinity and drought stress in maize. PHYSIOL Plant. 2021;172(2):684–95.
Gong Z, Xiong L, Shi H, Yang S, Herrera-Estrella LR, Xu G, Chao DY, Li J, Wang PY, Qin F, et al. Plant abiotic stress response and nutrient use efficiency. SCI CHINA LIFE SCI. 2020;63(5):635–74.
Munns R, Gilliham M. Salinity tolerance of crops - what is the cost? NEW PHYTOL. 2015;208(3):668–73.
Sun Y, Wang M, Mur L, Shen Q, Guo S. Unravelling the Roles of Nitrogen Nutrition in Plant Disease Defences. INT J MOL SCI 2020, 21(2).
Xun Z, Guo X, Li Y, Wen X, Wang C, Wang Y. Quantitative proteomics analysis of tomato growth inhibition by ammonium nitrogen. PLANT PHYSIOL BIOCH. 2020;154:129–41.
Chutia R, Scharfenberg S, Neumann S, Abel S, Ziegler J. Modulation of Phosphate Deficiency-Induced Metabolic Changes by Iron Availability in Arabidopsis thaliana. INT J MOL SCI 2021, 22(14).
Cui J, Tcherkez G. Potassium dependency of enzymes in plant primary metabolism. PLANT PHYSIOL BIOCH. 2021;166:522–30.
Wang Y, Chen YF, Wu WH. Potassium and phosphorus transport and signaling in plants. J INTEGR PLANT BIOL. 2021;63(1):34–52.
Song G, Li X, Munir R, Khan AR, Azhar W, Khan S, Gan Y. BnaA02.NIP6;1a encodes a boron transporter required for plant development under boron deficiency in Brassica napus. PLANT PHYSIOL BIOCH. 2021;161:36–45.
Pereira GL, Siqueira JA, Batista-Silva W, Cardoso FB, Nunes-Nesi A, Araujo WL. Boron: more than an essential element for land plants? FRONT PLANT SCI. 2020;11:610307.
Xue D, Jiang H, Deng X, Zhang X, Wang H, Xu X, Hu J, Zeng D, Guo L, Qian Q. Comparative proteomic analysis provides new insights into cadmium accumulation in rice grain under cadmium stress. J HAZARD MATER. 2014;280:269–78.
Essoh AP, Monteiro F, Pena AR, Pais MS, Moura M, Romeiras MM. Exploring glucosinolates diversity in Brassicaceae: a genomic and chemical assessment for deciphering abiotic stress tolerance. PLANT PHYSIOL BIOCH. 2020;150:151–61.
Ramirez D, Abellan-Victorio A, Beretta V, Camargo A, Moreno DA. Functional Ingredients From Brassicaceae Species: Overview and Perspectives. INT J MOL SCI 2020, 21(6).
Kumar V, Thakur JK, Prasad M. Histone acetylation dynamics regulating plant development and stress responses. CELL MOL LIFE SCI. 2021;78(10):4467–86.
Joseph JT, Shah JM. Biotic stress-induced epigenetic changes and transgenerational memory in plants. BIOLOGIA. 2022;77(8):2007–21.
Wang JX, Wang XM, Geng S, Singh SK, Wang YH, Pattanaik S, Yuan L. Genome-wide identification of hexokinase gene family in Brassica napus: structure, phylogenetic analysis, expression, and functional characterization (vol 248, pg 171, 2018). Planta. 2018;248(1):183.
Eberhardt RY, Bartholdson SJ, Punta M, Bateman A. The SHOCT domain: a widespread domain under-represented in model organisms. PLoS ONE. 2013;8(2):e57848.
Liu Y, von Wiren N. Ammonium as a signal for physiological and morphological responses in plants. J EXP BOT. 2017;68(10):2581–92.
Yang S, Hao D, Jin M, Li Y, Liu Z, Huang Y, Chen T, Su Y. Internal ammonium excess induces ROS-mediated reactions and causes carbon scarcity in rice. BMC PLANT BIOL. 2020;20(1):143.
Wimmer MA, Eichert T. Review: mechanisms for boron deficiency-mediated changes in plant water relations. PLANT SCI. 2013;203–204:25–32.
Zhang ZH, Zhou T, Tang TJ, Song HX, Guan CY, Huang JY, Hua YP. A multiomics approach reveals the pivotal role of subcellular reallocation in determining rapeseed resistance to cadmium toxicity. J EXP BOT. 2019;70(19):5437–55.
Zhang GB, Meng S, Gong JM. The Expected and Unexpected Roles of Nitrate Transporters in Plant Abiotic Stress Resistance and Their Regulation. INT J MOL SCI 2018, 19(11).
Mostofa MG, Rahman MM, Ghosh TK, Kabir AH, Abdelrahman M, Rahman KM, Mochida K, Tran LP. Potassium in plant physiological adaptation to abiotic stresses. PLANT PHYSIOL BIOCH. 2022;186:279–89.
Wang M, Zheng Q, Shen Q, Guo S. The critical role of potassium in plant stress response. INT J MOL SCI. 2013;14(4):7370–90.
Vance CP, Uhde-Stone C, Allan DL. Phosphorus acquisition and use: critical adaptations by plants for securing a nonrenewable resource. NEW PHYTOL. 2003;157(3):423–47.
Ma S, Zheng L, Liu X, Zhang K, Hu L, Hua Y, Huang J. Genome-Wide Identification of Brassicaceae Hormone-Related Transcription Factors and Their Roles in Stress Adaptation and Plant Height Regulation in Allotetraploid Rapeseed. INT J MOL SCI 2022, 23(15).
Zheng L, Ma S, Shen D, Fu H, Wang Y, Liu Y, Shah K, Yue C, Huang J. Genome-wide identification of Gramineae histone modification genes and their potential roles in regulating wheat and maize growth and stress responses. BMC PLANT BIOL. 2021;21(1):543.
Chalhoub B, Denoeud F, Liu S, Parkin IA, Tang H, Wang X, Chiquet J, Belcram H, Tong C, Samans B, et al. Plant genetics. Early allopolyploid evolution in the post-neolithic Brassica napus oilseed genome. Science. 2014;345(6199):950–3.
Song X, Wei Y, Xiao D, Gong K, Sun P, Ren Y, Yuan J, Wu T, Yang Q, Li X, et al. Brassica carinata genome characterization clarifies U’s triangle model of evolution and polyploidy inBrassica. PLANT PHYSIOL. 2021;186(1):388–406.
Liu S, Liu Y, Yang X, Tong C, Edwards D, Parkin IAP, Zhao M, Ma J, Yu J, Huang S et al. The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. NAT COMMUN 2014, 5(1).
Parkin IA, Koh C, Tang H, Robinson SJ, Kagale S, Clarke WE, Town CD, Nixon J, Krishnakumar V, Bidwell SL, et al. Transcriptome and methylome profiling reveals relics of genome dominance in the mesopolyploid Brassica oleracea. GENOME BIOL. 2014;15(6):R77.
Halldorsson BV, Hardarson MT, Kehr B, Styrkarsdottir U, Gylfason A, Thorleifsson G, Zink F, Jonasdottir A, Jonasdottir A, Sulem P, et al. Author correction: the rate of meiotic gene conversion varies by sex and age. NAT GENET. 2018;50(11):1616.
Wang L, Ahmad B, Liang C, Shi X, Sun R, Zhang S, Du G. Bioinformatics and expression analysis of histone modification genes in grapevine predict their involvement in seed development, powdery mildew resistance, and hormonal signaling. BMC PLANT BIOL. 2020;20(1):412.
Uga Y, Sugimoto K, Ogawa S, Rane J, Ishitani M, Hara N, Kitomi Y, Inukai Y, Ono K, Kanno N, et al. Control of root system architecture by DEEPER ROOTING 1 increases rice yield under drought conditions. NAT GENET. 2013;45(9):1097–102.
Xu GX, Guo CC, Shan HY, Kong HZ. Divergence of duplicate genes in exon-intron structure. P NATL ACAD SCI USA. 2012;109(4):1187–92.
Huang J, Ma S, Zhang K, Liu X, Hu L, Wang W, Zheng L. Genome-Wide Identification of Gramineae Brassinosteroid-Related Genes and Their Roles in Plant Architecture and Salt Stress Adaptation. INT J MOL SCI 2022, 23(10).
Kurita K, Sakamoto Y, Naruse S, Matsunaga TM, Arata H, Higashiyama T, Habu Y, Utsumi Y, Utsumi C, Tanaka M, et al. Intracellular localization of histone deacetylase HDA6 in plants. J PLANT RES. 2019;132(5):629–40.
Wu K, Zhang L, Zhou C, Yu CW, Chaikam V. HDA6 is required for jasmonate response, senescence and flowering in Arabidopsis. J EXP BOT. 2008;59(2):225–34.
Chen LT, Luo M, Wang YY, Wu K. Involvement of Arabidopsis histone deacetylase HDA6 in ABA and salt stress response. J EXP BOT. 2010;61(12):3345–53.
Luo M, Wang YY, Liu X, Yang S, Lu Q, Cui Y, Wu K. HD2C interacts with HDA6 and is involved in ABA and salt stress response in Arabidopsis. J EXP BOT. 2012;63(8):3297–306.
Hollender C, Liu Z. Histone deacetylase genes in Arabidopsis development. J INTEGR PLANT BIOL. 2008;50(7):875–85.
Hartl M, Fussl M, Boersema PJ, Jost JO, Kramer K, Bakirbas A, Sindlinger J, Plochinger M, Leister D, Uhrig G, et al. Lysine acetylome profiling uncovers novel histone deacetylase substrate proteins in Arabidopsis. MOL SYST BIOL. 2017;13(10):949.
Hu Y, Qin F, Huang L, Sun Q, Li C, Zhao Y, Zhou DX. Rice histone deacetylase genes display specific expression patterns and developmental functions. BIOCHEM BIOPH RES CO. 2009;388(2):266–71.
Ma X, Lv S, Zhang C, Yang C. Histone deacetylases and their functions in plants. PLANT CELL REP. 2013;32(4):465–78.
Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R. TBtools: an integrative Toolkit developed for interactive analyses of big Biological Data. MOL PLANT. 2020;13(8):1194–202.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. MOL BIOL EVOL. 2018;35(6):1547–9.
Zhu K, Xu S, Li K, Chen S, Zafar S, Cao W, Wang Z, Ding L, Yang Y, Li Y et al. Transcriptome analysis of the irregular shape of shoot apical meristem in dt (dou tou) mutant of Brassica napus L. MOL Breed 2019, 39(3).
Zhou T, Yue CP, Huang JY, Cui JQ, Liu Y, Wang WM, Tian C, Hua YP. Genome-wide identification of the amino acid permease genes and molecular characterization of their transcriptional responses to various nutrient stresses in allotetraploid rapeseed. BMC PLANT BIOL. 2020;20(1):151.
Zheng M, Hu M, Yang H, Tang M, Zhang L, Liu H, Li X, Liu J, Sun X, Fan S, et al. Three BnaIAA7 homologs are involved in auxin/brassinosteroid-mediated plant morphogenesis in rapeseed (Brassica napus L). PLANT CELL REP. 2019;38(8):883–97.
Wang X, Zheng M, Liu H, Zhang L, Chen F, Zhang W, Fan S, Peng M, Hu M, Wang H, et al. Fine-mapping and transcriptome analysis of a candidate gene controlling plant height in Brassica napus L. BIOTECHNOL BIOFUELS. 2020;13:42.
Li H, Cheng X, Zhang L, Hu J, Zhang F, Chen B, Xu K, Gao G, Li H, Li L, et al. An integration of genome-wide Association study and gene co-expression network analysis identifies candidate genes of stem lodging-related traits in Brassica napus. FRONT PLANT SCI. 2018;9:796.
Cui JQ, Hua YP, Zhou T, Liu Y, Huang JY, Yue CP. Global Landscapes of the Na+/H + Antiporter (NHX) Family Members Uncover their Potential Roles in Regulating the Rapeseed Resistance to Salt Stress. Int J Mol Sci 2020, 21(10).
Zhou T, Hua YP, Zhang BC, Zhang XQ, Zhou YH, Shi L, Xu FS. Low-boron tolerance strategies involving pectin-mediated cell Wall Mechanical Properties in Brassica napus. PLANT CELL PHYSIOL. 2017;58(11):1991–2005.
Hua YP, Feng YN, Zhou T, Xu FS. Genome-scale mRNA transcriptomic insights into the responses of oilseed rape (Brassica napus L.) to varying boron availabilities. PLANT SOIL. 2017;416(1–2):205–25.
This work was supported by the Chinese Postdoctoral Science Foundation (2021M692944), Research Start-Up Project (32212399 and 32213006), Famous Teachers in Central Plains (22610002), and Application of Molecular Design Breeding of Oil Crops and Intelligent Auxiliary Information System in Supercomputing Ecology, Henan Key Project of Science and Technology (202102110006).
Ethics approval and consent to participate
This article does not contain any studies with human participants or animals performed by any of the authors. No specific permits were required. All methods were in compliance with relevant institutional, national, and international guidelines and legislation.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Hu, LL., Zheng, LW., Zhu, XL. et al. Genome-wide identification of Brassicaceae histone modification genes and their responses to abiotic stresses in allotetraploid rapeseed. BMC Plant Biol 23, 248 (2023). https://doi.org/10.1186/s12870-023-04256-1
- Allotetraploid rapeseed
- Histone modification
- Abiotic stress