Next Article in Journal
Seasonal Dynamics of Soil Moisture in an Integrated-Crop-Livestock-Forestry System in Central-West Brazil
Next Article in Special Issue
A New Breeding Strategy towards Introgression and Characterization of Stay-Green QTL for Drought Tolerance in Sorghum
Previous Article in Journal
Evaluating the Financial Health of Agricultural Enterprises in the Conditions of the Slovak Republic Using Bankruptcy Models
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Agriculture
Volume 11
Issue 3
10.3390/agriculture11030244
agriculture-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Comprehensive Analysis of the Histone Deacetylase Gene Family in Chinese Cabbage (Brassica rapa): From Evolution and Expression Pattern to Functional Analysis of BraHDA3

by
Seung Hee Eom
and
Tae Kyung Hyun
*
Department of Industrial Plant Science and Technology, Chungbuk National University, Cheongju 28644, Korea
*
Author to whom correspondence should be addressed.
Agriculture 2021, 11(3), 244; https://doi.org/10.3390/agriculture11030244
Submission received: 24 February 2021 / Revised: 10 March 2021 / Accepted: 10 March 2021 / Published: 12 March 2021
(This article belongs to the Special Issue Biotic and Abiotic Stresses in Crop Plants)
Browse Figures
Versions Notes

Abstract

:
Histone deacetylases (HDACs) are known as erasers that remove acetyl groups from lysine residues in histones. Although plant HDACs play essential roles in physiological processes, including various stress responses, our knowledge concerning HDAC gene families and their evolutionary relationship remains limited. In Brassica rapa genome, we identified 20 HDAC genes, which are divided into three major groups: RPD3/HDA1, HD2, and SIR2 families. In addition, seven pairs of segmental duplicated paralogs and one pair of tandem duplicated paralogs were identified in the B. rapa HDAC (BraHDAC) family, indicating that segmental duplication is predominant for the expansion of the BraHDAC genes. The expression patterns of paralogous gene pairs suggest a divergence in the function of BraHDACs under various stress conditions. Furthermore, we suggested that BraHDA3 (homologous of Arabidopsis HDA14) encodes the functional HDAC enzyme, which can be inhibited by Class I/II HDAC inhibitor SAHA. As a first step toward understanding the epigenetic responses to environmental stresses in Chinese cabbage, our results provide a solid foundation for functional analysis of the BraHDAC family.

1. Introduction

Epigenetic regulation, including DNA and histone modifications, plays an essential role in controlling chromatin structure, with the nucleosome as the basic organization unit [1]. In eukaryotes, post-transcriptional histone modifications such as acetylation, methylation, ubiquitination, and sumoylation enable the regulation of mRNA expression by recruiting histone writers, erasers, and readers [2]. Among these modifications, histone acetylation in different lysine residues of histones by histone acetyltransferases (HATs) is required for opening the chromatin structure to allow binding of RNA polymerase II and transcription factors; in contrast, histone deacetylases (HDACs) maintain the homeotic balance of histone acetylation by removing the acetyl groups from hyperacetylated histones [1,3]. The highly dynamic balance of the histone acetylation and deacetylation is mostly maintained by the physical and functional interplay between HAT and HDAC activities [3], and serves as a major epigenetic regulatory mechanism for gene transcription and in turn governs numerous physiological processes [4].
In higher plants, HDACs are generally classified into three distinct groups: the reduced potassium dependence 3/histone deacetylase 1 (RPD3/HDA1) family, silent information regulator 2 (SIR2) family, and histone deacetylase 2 (HD2) family [5]. Since 18 HDACs were identified in the Arabidopsis thaliana genome [6], several HDACs have been identified and isolated from rice [7], soybean [8], grape [9], litchi [10], maize [11], and cotton [12]. These indicate that genome-wide mining and comparative analysis are powerful approach to identify genes and understand the evolutionary relationships of genes. Genetic and physiological studies on Arabidopsis have shown that HDA6, HDA9, HDA15, HDA19, HD2C, and AtSRT1 negatively regulate stress tolerance by regulating gene expression through histone acetylation [13]. The Arabidopsis HDA6 mutants showed enhanced induction of genes involved in the acetate biosynthesis pathway, which resulted in enhanced drought tolerance [14]. Similarly, the loss-of-function of HDA9 showed increased tolerance against salt and drought stresses [15]. In addition, HDA19-deficiency enhanced tolerance to high-salinity, drought, and heat stresses in Arabidopsis [16]. In contrast, HDA19 overexpression in Arabidopsis enhanced disease resistance against Alternaria brassicicola via increased transcription of jasmonic acid and ethylene responsive genes [17]. In rice plants, overexpression of OsHDA705, a member of the RPD3/HDA1 family, resulted in decreased salt and abscisic acid (ABA) stress tolerance during seed germination, whereas its overexpression resulted in enhanced osmotic stress resistance during the seedling stage [18]. The decrease in levels of histone H4 acetylation because of the overexpression of OsHDT701 (rice HD2 family) caused reduced resistance against rice pathogens Magnaporthe oryzae and Xanthomonas oryzae pv oryzae [7], whereas the overexpression of HDT701 in transgenic rice leads to increased tolerance to NaCl and polyethylene glycol stresses [19]. Furthermore, mitogen-activated protein kinase 3 phosphorylates the histone deacetylase HD2B to control the transcription of biotic stress response genes [20]. Collectively, these indicate the importance, diversity, and specificity of HDACs in response to different stimuli.
Despite the knowledge concerning HDACs, the evolutionary relationships and functional information regarding the HDAC family in Brassica rapa, a crop species of economic importance, remain largely unknown. In this study, we aimed at genome-wide identification of the HDAC gene family in B. rapa genome. Expansion and evolutionary mechanisms of this gene family were studied to explore the evolutionary relationships among B. rapa HDAC family (BraHDAC). Furthermore, we analyzed the expression pattern of BraHDACs under different abiotic stress condition, providing valuable information for further study of BraHDAC gene functions.

2. Materials and Methods

2.1. Identification and Characterization of HDAC Family in B. rapa

Phytozome v12.1 (https://phytozome.jgi.doe.gov; Brassica rapa FPsc v1.3) was used to identify HDACs in B. rapa genome using BLASTp with Arabidopsis HDACs as queries. The conserved domains, molecular weight, phylogenetic, and isoelectric point (pI) analyses and subcellular localization of putative BraHDACs has been conducted by Eom and Hyun [21]. A nomenclature system for BraHDACs (generic name: BraHDA1 to 16 for the RPD3/HDA1 family, BraHDT1 and 2 for the HD2 family, and BraSRT1 and 2 for the SIR2 family) was used as described by Hollender and Liu [6].

2.2. Gene Duplication and Calculating Duplication Times of BraHDAC Genes

The Plant Genome Duplication Database (http://chibba.agtec.uga.edu/duplication/) was used to identify the homologous pairs among BraHDAC genes. Using PAL2NAL (http://www.bork.embl.de/pal2nal/), the synonymous rate (Ks), nonsynonymous rate (Ka), and evolutionary constraint (Ka/Ks) of duplicated BraHDAC genes were analyzed. The divergence time (T) was calculated using the formula T = Ks/(2 × 1.5 × 10−8) × 10−6 million years ago (MYA; [22]).

2.3. Plant Growth and Treatment Conditions

Chinese cabbage (cultivar Chunkwang) seeds were grown in a growth chamber under long-day conditions at 22 °C and 60% relative humidity. After 4 weeks, Chinese cabbage plants were treated at 40 °C, and samples were harvested at the times indicated in Figure 2. Heat treatment was carried out three times, resulting in three biological replicates.

2.4. The Analysis of BraHDAC Expression

Total RNA was extracted, quantified, and qualified using nanodrop. cDNA was synthesized using ReverTra Ace qPCR RT Master Mix with gDNA Remover. Quantitative RT-PCR (qRT-PCR) was performed using the SYBR Green Real-time PCR Master Mix. The expression levels of each gene were normalized to the internal reference gene actin, and then was expressed relative to its value at 0 h. The specific primer pairs were designed using GenScript Real-time PCR Primer Design tool, and are listed in Table S1.

2.5. Production of Recombinant BraHDA3 in E. coli, and the Analysis of HDAC Activity

The PCR product of BraHDA3 was ligated into the pCR-Blunt, and subcloned into the Escherichia coli expression vector pPROEX-Htb vector that produces a histidine-fused translation product. To produce his tagged-BraHDA3 fusion proteins, E. coli strain BL21 harboring pPROEX-BraHDA3 plasmid was incubated at 37 °C. Once the optical density at 600 nm reached 0.5, fusion proteins were produced by 0.1 mM isopropyl β- d-1-thiogalactopyranoside (IPTG) treatment for 16 h at 30 °C. The fusion protein was purified using Ni-NTA affinity purification, and the HDAC activity was determined using HDAC Activity Colorimetric Assay Kit (BioVision, Milpitas, CA, USA). HDAC activity was calculated based on a standard curve generated with Ac-Lys-pNA (deacetylated standard).

2.6. Statistical Analyses

The experimental results are expressed as means ± SE. SPSS Statistics version 25 (IBM Microsoft, New York, NY, USA) was used to carry out statistical analyses (Duncan multiple range tests). P values < 0.05 were regarded as statistically significant difference.

3. Results and Discussion

3.1. Identification of Histone Deacetylases in Chinese Cabbage

Since HDACs have been introduced as erasers of epigenetic acetylation marks on lysine residues in histones, a number of plant HDAC genes have been identified from a variety of plant species [3]. From the B. rapa genome, we identified 20 HDAC members (Table 1). Among them, BraHDA6 was the largest protein, with 651 amino acids (AA), whereas BraHDA1 was identified as the smallest protein with 254 AA. The molecular weight of BraHDACs varied according to protein size, ranging from 28.8 KDa to 72.1 KDa, and their pI varied from 4.52 (BraHDT2) to 8.96 (BraSRT1). Subcellular locations as predicted by WoLF PSORT indicated that BraHDACs were potentially localized in the cytosol, nucleus, chloroplast, mitochondria, and peroxisome (Table 1). Similar to the Arabidopsis and soybean HD2 proteins [23,24], the B. rapa HD2 family were predicted as nuclear proteins (Table 1). In addition, BraSRT1 and BraSRT2 appeared to be localized in the mitochondria and nucleus, respectively, similar to the rice SRT (silent information regulator) family [25]. Taken together, the various localization patterns of BraHDACs indicated that they may have distinct roles.
Plant RPD3/HDA1, HD2, and SIR2 families have been characterized by the presence of conserved domains such as the Hist_deacetyl domain (PF00850), SIR2 domain (PF02146), and ZnF_C2H2 domain (SM000355) [11,12]. As shown in Figure 1A, BraHDACs were divided into three major classes: RPD3/HDA1, HD2, and SIR2 families. The B. rapa genome encodes 16 proteins (Figure 1 and Figure S1) clustered in the RPD3/HDA1 family contained the conserved Hist_deacetyl domain. In addition, two SIR2 family HDACs with highly conserved SIR2 domain and two HD2 family HDACs with conserved zinc finger C2H2 type domain (ZnF_C2H2) were identified. Another type of zinc finger domain, ZnF_RTZ (SM000547), has been found in BraHDA07, indicating that BraHDT1, BraHDT1, and BraHDA07 have high binding affinity to DNA or modulate protein-protein interactions mediated by zinc finger domains [26]. The presence of domains similar to those in other plant HDACs suggested that all putative BraHDACs belong to the histone deacetylase family.

3.2. Gene Duplication Analysis of BraHDACs

Gene duplication has long been regarded as a key evolutionary process for acquiring new genes and genetic novelty [27]. It could occur on any of these scales: whole genome duplication, segmental duplication, and single gene duplication including tandem duplication [28]. It has been shown that some HDACs expanded in the Gossypium and soybean genomes partly because of segmental duplication events [12,24]. In the B. rapa genome, 20 BraHDACs were asymmetrically distributed on 7 chromosomes (Table 1). Chromosome A06 contains the highest density of HDACs with six members (BraHDA8 to 13), followed by chromosome A01 with four genes (BraHDA1 to 3, BraSRT1), whereas only one BraHDAC was present on chromosome A05 (Table 1). A tandem duplication is characterized as two or more adjacent homologous genes located physically “close” to one another; tandem duplicated pairs are paralogous genes on the same chromosome within a 50 kb physical distance [29], indicating that only one set (BraHDA8 and 9) was tandemly duplicated on chromosome A06 separated by 14,097 bp (Table 1) with 85.7% and 75.4% of sequence identity at the nucleotide and protein levels. To further determine whether segmental duplication events also contributed to the expansion of BraHDACs, we analyzed the PGDD database combined with phylogenetic analysis. As shown in Figure 1B, seven duplicated BraHDAC gene pairs were identified, indicating that segmental duplication is a major contributory event leading to the expansion of BraHDAC genes in B. rapa genome. The highest frequency of segmental duplication events occurred between chromosomes A06 and 09, which contained two segmentally duplicated gene pairs.
To examine the evolutionary selection process of eight duplicated gene pairs in B. rapa genome, we analyzed the evolutionary constraint (Ka/Ks). Ka/Ks ratios ranged from 0.0895 to 0.8713, and the average Ka/Ks value of the duplicated gene pairs was 0.404 (Table 2). This indicates that these duplicated gene pairs were subjected to negative selection (purifying selection), similar to duplicated HDAC gene pairs in Gossypium and soybean genomes [12,24]. The segmental duplications of HDAC genes in B. rapa genome occurred between 12.73 MYA and 42.04 MYA, with the average value being 29.19 MYA. The Ks of tandem duplications of BraHDAC genes was 0.1640, dating the duplication event at 5.47 MYA (Table 2). Brassicaceae, including Arabidopsis and Brassica, have been shown to have undergone three rounds of whole genome duplication, including γ triplication (135 MYA) and β (90–100 MYA) and α (24–40 MYA) duplications [30,31,32]. In addition, B. rapa exhibits this complex history, with the addition of a whole-genome triplication occurring between 13 and 17 MYA [33], indicating that the segmental duplications of BraHDAC genes, except BraHDA 5/7 and BraHDA 12/15 gene pairs, occurred before whole-genome triplication of the Brassica genome.

3.3. Effects of Heat Stress on Acetylation Status of Histone H3 in Chinese Cabbage

Heat stress induces excessive generation of reactive oxygen species, resulting in alterations in plant growth, development, photosynthesis, physiological processes, and production yield [34]. To evaluate the efficacy of heat treatment (40 °C), the physiological response of Chinese cabbage including lipid peroxidation and proline content was analyzed. As shown in Figure 2A, heat stress induced the accumulation of malondialdehyde (MDA), which is an indicator of lipid peroxidation [35]. In higher plants, the accumulation of proline, an amino acid and a compatible solute, is a common phenomenon in response to various environmental stresses [36]. We observed an increase in proline levels in response to heat stress (Figure 2B). The accumulated proline purportedly acts as a compatible osmolyte, molecular chaperone, cytosolic pH buffer, cell redox balancer, free radical scavenger, and stabilizer for subcellular structures during various stresses, especially osmotic and salt stresses [36,37], and improves resistance to various stresses [38]. However, proline accumulation in response to heat stress leads to increased production of reactive oxygen species in the mitochondria and inhibits ABA and ethylene biosynthesis, resulting in increased sensitivity to heat stress [36]. Taken together, these physiological changes suggested that the heat stress treatment was effective, and the leaves of heat-treated plants were harvested for further analyses.
To gain insight into the possible involvement of BraHDACs during their response to heat stress, the expression profiles of BraHDACs were analyzed. As shown in Figure 2C, exposure of Chinese cabbage plants to heat stress induced or repressed the expression of BraHDACs. After 48-h exposure to heat stress, expression levels of BraHDA2, BraHDA4, BraHDA5, BraHDA6, BraHDA7, BraHDA11, BraHDA12, BraHDA15, BraHDA16, BraHDT2, and BraSRT2 had increased, whereas those of other BraHDACs decreased, remained unchanged, or were undetected. In addition, the transcript level of BraHDA13, BraHDA15, BraHDA16, BraHDT2, and BraSRT2 was up-regulated in response to drought, salt, and wounding, whereas the decreasing level of BraHDA2, BraHDA3, BraHDA9, and BraHDA14 transcripts was observed (Figure S2). Under heat-stress condition, BraHDA4 and BraHDA14 exhibited a different expression pattern, although BraHDA4 and BraHDA14 were clustered into same sub-group in RPD3/HDA1 family, indicating that the duplication of those genes resulted in divergence of their expression pattern and function. Similar with BraHDACs, several Arabidopsis HDACs and maize HDACs in the class II of RPD3/HDA1 family were induced by heat treatment [11,39]. In addition, Arabidopsis, maize and cotton HDACs were differentially expressed in response to hormones and stresses, suggesting the functional divergence and specific regulation of HDACs in response to a particular stress [11,12,39]. According to the abiotic stress induced expression pattern of BraHDACs, we found that BraHDA8 was not expressed in leaves and on abiotic stresses. Based on transcriptome analysis, transcript level of BraHDA8 has been detected only in the silique [40], indicating that BraHDA8 should be expressed in specific tissue. In addition, the expression pattern of paralogous gene pair (BraHDA4/14) diverged, whereas three paralogous gene pairs (BraHDA4/6, BraHDA5/7, and BraHDA12/15) exhibited similar expression patterns (Figure 2C) under heat stress condition. However, these expression patterns were altered in response to different stresses including drought, salt, and wounding (Figure S2), indicating that BraHDACs had diversified substantially, which might be due to the different divergent fates after duplications.

3.4. Histone Deacetylase Activity of BarHDA3

Although many fundamental questions remain unanswered regarding the physiological function of plant HDACs in response to environmental stresses, it has been shown that RPD3/HDA1 family including Arabidopsis HDA5/14/15/18/19 positively or negatively regulated to the environmental-stress tolerances [41]. Among them, AtHDA14 controlled the activation state of RuBisCO under low-light condition, and is involved in biosynthesis of melatonin, which plays an important role in plant response to environmental stresses [42,43]. Under heat stress condition, exogenous melatonin improved the antioxidant defense systems, resulting in the suppression of heat stress-induced damage [44]. In Chinese cabbage, BraHDA3, homologous of AtHDA14 (Figure S1), dramatically down-regulated in response to heat stress (Figure 2C). This indicates that down-regulation of BraHDA3 might be required for the expression of heat response genes, although the physiological function of BraHDA3 needs to be characterized. To analyze the function of BraHDA3, recombinant BraHDA3 (BraHDA3-His) was prepared using E. coli expression system. As shown in Figure 3, purified recombinant BraHDA3 protein gave an activity of 18.1 µM/100 µg of protein. When HDAC inhibitor SAHA (5 µM) was mixed with purified recombinant BraHDA3 protein, the HDAC activity decreased to 3.8 µM/100 µg of protein. This indicates that BraHDA3 encodes the functional HDAC enzyme, which can be inhibited by Class I/II HDAC inhibitor SAHA.

4. Conclusions

Despite the knowledge concerning HDACs, evolutionary and functional information regarding HDACs in B. rapa remains relatively unknown. In this study, we identified 20 putative BraHDAC genes divided into three major classes: RPD3/HDA1, HD2, and SIR2 families. Most HDAC gene duplications in B. rapa appeared to have been caused by segmental duplication, which occurred between 12.73 MYA and 42.04 MYA. Several BraHDACs were up- or down-regulated by abiotic stress treatments, indicating that the variation of BraHDAC transcription levels may contribute to alteration in histone H3 acetylation levels in response to various stresses. Our results provide a solid foundation for further understanding of the underlying evolutionary mechanisms in the HDAC family and will provide the basis for future research on functional characterization of the BraHDAC family in response to environmental stresses.

Supplementary Materials

The following are available online at https://www.mdpi.com/2077-0472/11/3/244/s1, Table S1: Primer sequences for PCR analysis, Figure S1. Phylogenetic relationships of Arabidopsis, maize, and Brassica rapa HDACs, Figure S2: Expression pattern of BraHDAC genes in response to wounding, salt, and drought stresses.

Author Contributions

Conceptualization, S.H.E. and T.K.H.; investigation, S.H.E.; writing—original draft preparation, S.H.E. and T.K.H.; writing—review and editing, S.H.E. and T.K.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was carried out with the support of “Cooperative Research Program for Agriculture Science and Technology Development (Project No. PJ01501905)” Rural Development Administration, Republic of Korea.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available within the article.

Conflicts of Interest

The authors confirm that there is no conflict of interest to declare.

References

  1. Lämke, J.; Bäurle, I. Epigenetic and chromatin-based mechanisms in environmental stress adaptation and stress memory in plants. Genome Biol. 2017, 18, 1–11. [Google Scholar] [CrossRef] [PubMed]
  2. Berr, A.; Shafiq, S.; Shen, W.H. Histone modifications in transcriptional activation during plant development. Biochim. Biophys. Acta Gene Regul. Mech. 2011, 1809, 567–576. [Google Scholar] [CrossRef]
  3. Chen, X.; Ding, A.B.; Zhong, X. Functions and mechanisms of plant histone deacetylases. Sci. China Life Sci. 2020, 63, 206–216. [Google Scholar] [CrossRef] [Green Version]
  4. Smith, L.M.; Pontes, O.; Searle, I.; Yelina, N.; Yousafzai, F.K.; Herr, A.J.; Pikaard, C.S.; Baulcombe, D.C. An SNF2 protein associated with nuclear RNA silencing and the spread of a silencing signal between cells in Arabidopsis. Plant Cell 2007, 19, 1507–1521. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Pandey, R.; Müller, A.; Napoli, C.A.; Selinger, D.A.; Pikaard, C.S.; Richards, E.J.; Bender, J.; Mount, D.W.; Jorgensen, R.A. Analysis of histone acetyltransferase and histone deacetylase families of Arabidopsis thaliana suggests functional diversification of chromatin modification among multicellular eukaryotes. Nucleic Acids Res. 2002, 30, 5036–5055. [Google Scholar] [CrossRef] [Green Version]
  6. Hollender, C.; Liu, Z. Histone deacetylase genes in arabidopsis development. J. Integr. Plant Biol. 2008, 57, 875–885. [Google Scholar] [CrossRef] [PubMed]
  7. Hu, Y.; Qin, F.; Huang, L.; Sun, Q.; Li, C.; Zhao, Y.; Zhou, D.X. Rice histone deacetylase genes display specific expression patterns and developmental functions. Biochem. Biophys. Res. Commun. 2009, 388, 266–271. [Google Scholar] [CrossRef] [PubMed]
  8. Liew, L.C.; Singh, M.B.; Bhalla, P.L. An RNA-Seq Transcriptome Analysis of Histone Modifiers and RNA Silencing Genes in Soybean during Floral Initiation Process. PLoS ONE 2013, 8, e77502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Aquea, F.; Timmermann, T.; Arce-Johnson, P. Analysis of histone acetyltransferase and deacetylase families of Vitis vinifera. Plant Physiol. Biochem. 2010, 48, 194–199. [Google Scholar] [CrossRef] [PubMed]
  10. 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. [Google Scholar] [CrossRef]
  11. Zhang, K.; Yu, L.; Pang, X.; Cao, H.; Si, H.; Zang, J.; Xing, J.; Dong, J. In silico analysis of maize HDACs with an emphasis on their response to biotic and abiotic stresses. PeerJ 2020, 8, e8539. [Google Scholar] [CrossRef] [Green Version]
  12. Imran, M.; Shafiq, S.; Naeem, M.K.; Widemann, E.; Munir, M.Z.; Jensen, K.B.; Wang, R.R.C. Histone deacetylase (HDAC) gene family in allotetraploid cotton and its diploid progenitors: In silico identification, molecular characterization, and gene expression analysis under multiple abiotic stresses, DNA damage and phytohormone treatments. Int. J. Mol. Sci. 2020, 21, 321. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Hu, Y.; Lu, Y.; Zhao, Y.; Zhou, D.X. Histone Acetylation Dynamics Integrates Metabolic Activity to Regulate Plant Response to Stress. Front. Plant Sci. 2019, 10, 1236. [Google Scholar] [CrossRef] [PubMed]
  14. Kim, J.M.; To, T.K.; Matsui, A.; Tanoi, K.; Kobayashi, N.I.; Matsuda, F.; Habu, Y.; Ogawa, D.; Sakamoto, T.; Matsunaga, S.; et al. Acetate-mediated novel survival strategy against drought in plants. Nat. Plants 2017, 3. [Google Scholar] [CrossRef]
  15. Zheng, Y.; Ding, Y.; Sun, X.; Xie, S.; Wang, D.; Liu, X.; Su, L.; Wei, W.; Pan, L.; Zhou, D.X. Histone deacetylase HDA9 negatively regulates salt and drought stress responsiveness in Arabidopsis. J. Exp. Bot. 2016, 67, 1703–1713. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Ueda, M.; Matsui, A.; Nakamura, T.; Abe, T.; Sunaoshi, Y.; Shimada, H.; Seki, M. Versatility of HDA19-deficiency in increasing the tolerance of Arabidopsis to different environmental stresses. Plant Signal. Behav. 2018, 13, e1475808. [Google Scholar] [CrossRef] [PubMed]
  17. Zhou, C.; Zhang, L.; Duan, J.; Miki, B.; Wu, K. Histone Deacetylase19 is involved in jasmonic acid and ethylene signaling of pathogen response in Arabidopsis. Plant Cell 2005, 17, 1196–1204. [Google Scholar] [CrossRef] [Green Version]
  18. Zhao, J.; Li, M.; Gu, D.; Liu, X.; Zhang, J.; Wu, K.; Zhang, X.; Teixeira Da Silva, J.A.; Duan, J. Involvement of rice histone deacetylase HDA705 in seed germination and in response to ABA and abiotic stresses. Biochem. Biophys. Res. Commun. 2016, 470, 439–444. [Google Scholar] [CrossRef]
  19. Zhao, J.; Zhang, J.; Zhang, W.; Wu, K.; Zheng, F.; Tian, L.; Liu, X.; Duan, J. Expression and functional analysis of the plant-specific histone deacetylase HDT701 in rice. Front. Plant Sci. 2015, 5. [Google Scholar] [CrossRef] [Green Version]
  20. Latrasse, D.; Jégu, T.; Li, H.; de Zelicourt, A.; Raynaud, C.; Legras, S.; Gust, A.; Samajova, O.; Veluchamy, A.; Rayapuram, N.; et al. MAPK-triggered chromatin reprogramming by histone deacetylase in plant innate immunity. Genome Biol. 2017, 18, 131. [Google Scholar] [CrossRef]
  21. Eom, S.H.; Hyun, T.K. Histone acetyltransferases (HATs) in chinese cabbage: Insights from histone H3 acetylation and expression profiling of HATs in response to abiotic stresses. J. Am. Soc. Hortic. Sci. 2018, 143, 296–303. [Google Scholar] [CrossRef] [Green Version]
  22. Koch, M.A.; Haubold, B.; Mitchell-Olds, T. Comparative evolutionary analysis of chalcone synthase and alcohol dehydrogenase loci in Arabidopsis, Arabis, and related genera (Brassicaceae). Mol. Biol. Evol. 2000, 17, 1483–1498. [Google Scholar] [CrossRef] [Green Version]
  23. Zhou, C.; Labbe, H.; Sridha, S.; Wang, L.; Tian, L.; Latoszek-Green, M.; Yang, Z.; Brown, D.; Miki, B.; Wu, K. Expression and function of HD2-type histone deacetylases in Arabidopsis development. Plant J. 2004, 38, 715–724. [Google Scholar] [CrossRef]
  24. Yang, C.; Shen, W.; Chen, H.; Chu, L.; Xu, Y.; Zhou, X.; Liu, C.; Chen, C.; Zeng, J.; Liu, J.; et al. Characterization and subcellular localization of histone deacetylases and their roles in response to abiotic stresses in soybean. BMC Plant Biol. 2018, 18, 226. [Google Scholar] [CrossRef] [PubMed]
  25. Huang, L.; Sun, Q.; Qin, F.; Li, C.; Zhao, Y.; Zhou, D.X. Down-regulation of a Silent Information Regulator2-related histone deacetylase gene, OsSRT1, induces DNA fragmentation and cell death in rice. Plant Physiol. 2007, 144, 1508–1519. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Giesecke, A.V.; Fang, R.; Joung, J.K. Synthetic protein-protein interaction domains created by shuffling Cys 2His2 zinc-fingers. Mol. Syst. Biol. 2006, 2, 0011. [Google Scholar] [CrossRef] [Green Version]
  27. Panchy, N.; Lehti-Shiu, M.; Shiu, S.H. Evolution of gene duplication in plants. Plant Physiol. 2016, 171, 2294–2316. [Google Scholar] [CrossRef] [Green Version]
  28. Freeling, M. Bias in plant gene content following different sorts of duplication: Tandem, whole-genome, segmental, or by transposition. Annu. Rev. Plant Biol. 2009, 60, 433–453. [Google Scholar] [CrossRef]
  29. Cannon, S.B.; Mitra, A.; Baumgarten, A.; Young, N.D.; May, G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 2004, 4, 10. [Google Scholar] [CrossRef] [Green Version]
  30. Blanc, G.; Wolfe, K.H. Widespread paleopolyploidy in model plant species inferred from age distributions of duplicate genes. Plant Cell 2004, 16, 1667–1678. [Google Scholar] [CrossRef] [Green Version]
  31. Lysak, M.A.; Koch, M.A.; Pecinka, A.; Schubert, I. Chromosome triplication found across the tribe Brassiceae. Genome Res. 2005, 15, 516–525. [Google Scholar] [CrossRef] [Green Version]
  32. Du, J.; Hu, S.; Yu, Q.; Wang, C.; Yang, Y.; Sun, H.; Yang, Y.; Sun, X. Genome-wide identification and characterization of BrrTCP transcription factors in Brassica rapa ssp. rapa. Front. Plant Sci. 2017, 8, 1588. [Google Scholar] [CrossRef] [Green Version]
  33. Cheng, F.; Wu, J.; Fang, L.; Sun, S.; Liu, B.; Lin, K.; Bonnema, G.; Wang, X. Biased gene fractionation and dominant gene expression among the subgenomes of Brassica rapa. PLoS ONE 2012, 7, e36442. [Google Scholar] [CrossRef] [Green Version]
  34. Hasanuzzaman, M.; Nahar, K.; Alam, M.M.; Roychowdhury, R.; Fujita, M. Physiological, biochemical, and molecular mechanisms of heat stress tolerance in plants. Int. J. Mol. Sci. 2013, 14, 9643–9684. [Google Scholar] [CrossRef]
  35. Morales, M.; Munné-Bosch, S. Malondialdehyde: Facts and artifacts. Plant Physiol. 2019, 180, 1246–1250. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Lv, W.T.; Lin, B.; Zhang, M.; Hua, X.J. Proline accumulation is inhibitory to arabidopsis seedlings during heat stress. Plant Physiol. 2011, 156, 1921–1933. [Google Scholar] [CrossRef] [Green Version]
  37. Szabados, L.; Savouré, A. Proline: A multifunctional amino acid. Trends Plant Sci. 2010, 15, 89–97. [Google Scholar] [CrossRef]
  38. Hayat, S.; Hayat, Q.; Alyemeni, M.N.; Wani, A.S.; Pichtel, J.; Ahmad, A. Role of proline under changing environments: A review. Plant Signal. Behav. 2012, 7, 1456–1466. [Google Scholar] [CrossRef] [Green Version]
  39. Alinsug, M.V.; Yu, C.W.; Wu, K. Phylogenetic analysis, subcellular localization, and expression patterns of RPD3/HDA1 family histone deacetylases in plants. BMC Plant Biol. 2009, 9, 37. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Tong, C.; Wang, X.; Yu, J.; Wu, J.; Li, W.; Huang, J.; Dong, C.; Hua, W.; Liu, S. Comprehensive analysis of RNA-seq data reveals the complexity of the transcriptome in Brassica rapa. BMC Genom. 2013, 14, 689. [Google Scholar] [CrossRef] [Green Version]
  41. Ueda, M.; Matsui, A.; Watanabe, S.; Kobayashi, M.; Saito, K.; Tanaka, M.; Ishida, J.; Kusano, M.; Seo, M.; Seki, M. Transcriptome analysis of the hierarchical response of histone deacetylase proteins that respond in an antagonistic manner to salinity stress. Front. Plant. Sci. 2019, 10, 1323. [Google Scholar] [CrossRef]
  42. Hartl, M.; Füßl, M.; Boersema, P.J.; Jost, J.O.; Kramer, K.; Bakirbas, A.; Sindlinger, J.; Plöchinger, M.; Leister, D.; Uhrig, G.; et al. Lysine acetylome profiling uncovers novel histone deacetylase substrate proteins in Arabidopsis. Mol. Syst. Biol. 2017, 13, 949. [Google Scholar] [CrossRef] [PubMed]
  43. Lee, K.; Lee, H.Y.; Back, K. Rice histone deacetylase 10 and Arabidopsis histone deacetylase 14 genes encode N-acetylserotonin deacetylase, which catalyzes conversion of N-acetylserotonin into serotonin, a reverse reaction for melatonin biosynthesis in plants. J. Pineal Res. 2018, 64, e12460. [Google Scholar] [CrossRef] [PubMed]
  44. Buttar, Z.A.; Wu, S.N.; Arnao, M.B.; Wang, C.; Ullah, I.; Wang, C. Melatonin suppressed the heat stress-induced damage in wheat seedlings by modulating the antioxidant machinery. Plants 2020, 9, 809. [Google Scholar] [CrossRef] [PubMed]
Agriculture 11 00244 g001 550
Figure 1. Phylogenetic tree, domain architectures (A) and synteny analysis (B) of the Brassica rapa HDAC family. Using the neighbor-joining method, phylogenetic tree was generated. The conserved domains in the putative BraHDACs were analyzed with the SMART program. Paralogous gene pairs were analyzed using plant genome duplication database.
Figure 1. Phylogenetic tree, domain architectures (A) and synteny analysis (B) of the Brassica rapa HDAC family. Using the neighbor-joining method, phylogenetic tree was generated. The conserved domains in the putative BraHDACs were analyzed with the SMART program. Paralogous gene pairs were analyzed using plant genome duplication database.
Agriculture 11 00244 g001
Agriculture 11 00244 g002 550
Figure 2. Physiological response to heat stress in Chinese cabbage. Changes in malondialdehyde (A) and proline (B) levels after exposure to heat stress (40 °C) were determined. (C) Expression patterns of Brassica rapa histone deacetylases (BraHDACs) in response to heat stress. Expression is indicated as a log2 ratio of experimental treatments relative to the values at 0 h. The bars represent the mean ± standard error across three independent experiments. Different letters indicate statistically significant differences (p < 0.05).
Figure 2. Physiological response to heat stress in Chinese cabbage. Changes in malondialdehyde (A) and proline (B) levels after exposure to heat stress (40 °C) were determined. (C) Expression patterns of Brassica rapa histone deacetylases (BraHDACs) in response to heat stress. Expression is indicated as a log2 ratio of experimental treatments relative to the values at 0 h. The bars represent the mean ± standard error across three independent experiments. Different letters indicate statistically significant differences (p < 0.05).
Agriculture 11 00244 g002
Agriculture 11 00244 g003 550
Figure 3. HDAC activity of recombinant BraHDA3. SDS-PAGE exhibited the purified recombinant BraHDA3 proteins. BraHDA3 enzyme activity was analyzed using HDAC Activity Colorimetric Assay Kit. Different letters indicate statistically significant differences (p < 0.05).
Figure 3. HDAC activity of recombinant BraHDA3. SDS-PAGE exhibited the purified recombinant BraHDA3 proteins. BraHDA3 enzyme activity was analyzed using HDAC Activity Colorimetric Assay Kit. Different letters indicate statistically significant differences (p < 0.05).
Agriculture 11 00244 g003
Table
Table 1. Histone deacetylase gene family in Brassica rapa.
Table 1. Histone deacetylase gene family in Brassica rapa.
NameGene IDLocationCDS (bp)AAIntron Nr.pIkDaSubcellular Localization
BraHDA1Brara.A02927.1A01:24968034..2497035576525474.9828.8Cytosol
BraHDA2Brara.A03767.1A01:29972774..29976555150950255.1756.0Nucleus
BraHDA3Brara.A00453.1A01:2203688..2206445127842586.0645.9Chloroplast
BraHDA4Brara.B03745.1A02:30976246..30978653166855585.9362.0Cytosol
BraHDA5Brara.C03690.1A03:19011956..190149351797598136.0065.8Nucleus
BraHDA6Brara.C04226.1A03:22262637..222660321956651125.2372.1Cytosol
BraHDA7Brara.E02343.1A05:21170590..211746131638545165.6059.8Nucleus
BraHDA8Brara.F01687.1A06:10370998..10372699123040944.9745.4Cytosol
BraHDA9Brara.F01690.1A06:10386796..10389169120039945.8244.4Cytosol
BraHDA10Brara.F01971.1A06:14950658..14953055127242348.7047.6Mitochondrial
BraHDA11Brara.F02002.1A06:16050873..160538591281426135.3248.7Cytosol
BraHDA12Brara.F02210.1A06:18915988..1891823792130644.9634.3Nucleus
BraHDA13Brara.F02773.1A06:22749037..22751166111937296.1742.0Cytosol
BraHDA14Brara.I00647.1A09:3621688..36254701959652125.3271.5Peroxisome
BraHDA15Brara.I00727.1A09:4109663..4112303141347054.9452.7Nucleus
BraHDA16Brara.I05307.1A09:43269577..43271821116138625.6842.2Cytosol
BraHDT1Brara.A02332.1A01:17662954..1767106178626176.5228.0Nucleus
BraHDT2Brara.J02801.1A10:19273043..1927496785828564.5231.2Nucleus
BraSRT1Brara.B00303.1A02:1417334..1419795104735788.9639.7Mitochondrial
BraSRT2Brara.C01390.1A03:6592244..65955151407468136.4952.4Nucleus
Table
Table 2. Divergence between paralogous histone acetyltransferase gene pairs in Brassica rapa.
Table 2. Divergence between paralogous histone acetyltransferase gene pairs in Brassica rapa.
Gene 1Gene 2SNKsKaKa/KsMya
BraHDA8BraHDA9298.2898.80.16400.14290.87135.47
BraHDA1BraHDA5176.0574.01.26120.62440.495142.04
BraHDA1BraHDA7165.2545.81.20660.41680.345440.22
BraHDA4BraHDA6371.31194.70.80110.31650.395126.70
BraHDA4BraHDA14399.21229.80.88700.38670.436029.57
BraHDA5BraHDA7319.5925.50.40430.13440.332413.48
BraHDA15BraHDA8302.9906.11.18710.31720.267239.57
BraHDA15BraHDA12215.5693.50.38180.03420.089512.73
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Eom, S.H.; Hyun, T.K. Comprehensive Analysis of the Histone Deacetylase Gene Family in Chinese Cabbage (Brassica rapa): From Evolution and Expression Pattern to Functional Analysis of BraHDA3. Agriculture 2021, 11, 244. https://doi.org/10.3390/agriculture11030244

AMA Style

Eom SH, Hyun TK. Comprehensive Analysis of the Histone Deacetylase Gene Family in Chinese Cabbage (Brassica rapa): From Evolution and Expression Pattern to Functional Analysis of BraHDA3. Agriculture. 2021; 11(3):244. https://doi.org/10.3390/agriculture11030244

Chicago/Turabian Style

Eom, Seung Hee, and Tae Kyung Hyun. 2021. "Comprehensive Analysis of the Histone Deacetylase Gene Family in Chinese Cabbage (Brassica rapa): From Evolution and Expression Pattern to Functional Analysis of BraHDA3" Agriculture 11, no. 3: 244. https://doi.org/10.3390/agriculture11030244

APA Style

Eom, S. H., & Hyun, T. K. (2021). Comprehensive Analysis of the Histone Deacetylase Gene Family in Chinese Cabbage (Brassica rapa): From Evolution and Expression Pattern to Functional Analysis of BraHDA3. Agriculture, 11(3), 244. https://doi.org/10.3390/agriculture11030244

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop