Skip to main content

Genome-wide investigation of histone acetyltransferase gene family and its responses to biotic and abiotic stress in foxtail millet (Setaria italica [L.] P. Beauv)



Modification of histone acetylation is a ubiquitous and reversible process in eukaryotes and prokaryotes and plays crucial roles in the regulation of gene expression during plant development and stress responses. Histone acetylation is co-regulated by histone acetyltransferase (HAT) and histone deacetylase (HDAC). HAT plays an essential regulatory role in various growth and development processes by modifying the chromatin structure through interactions with other histone modifications and transcription factors in eukaryotic cells, affecting the transcription of genes. Comprehensive analyses of HAT genes have been performed in Arabidopsis thaliana and Oryza sativa. However, little information is available on the HAT genes in foxtail millet (Setaria italica [L.] P. Beauv).


In this study, 24 HAT genes (SiHATs) were identified and divided into four groups with conserved gene structures via motif composition analysis. Phylogenetic analysis of the genes was performed to predict functional similarities between Arabidopsis thaliana, Oryza sativa, and foxtail millet; 19 and 2 orthologous gene pairs were individually identified. Moreover, all identified HAT gene pairs likely underwent purified selection based on their non-synonymous/synonymous nucleotide substitutions. Using published transcriptome data, we found that SiHAT genes were preferentially expressed in some tissues and organs. Stress responses were also examined, and data showed that SiHAT gene transcription was influenced by drought, salt, low nitrogen, and low phosphorus stress, and that the expression of four SiHATs was altered as a result of infection by Sclerospora graminicola.


Results indicated that histone acetylation may play an important role in plant growth and development and stress adaptations. These findings suggest that SiHATs play specific roles in the response to abiotic stress and viral infection. This study lays a foundation for further analysis of the biological functions of SiHATs in foxtail millet.

Peer Review reports


Plant nuclear DNA is organized into a DNA-protein structure called chromatin. The central nucleosome consists of 147 bp DNA, usually wrapped on an octamer of histones, comprising two copies each of the core histones H2A, H2B, H3, and H4 [1]. Histone proteins can be extensively modified at their N-terminal tails, which protrude from the core structure of nucleosomes and function as preferred targets for the histone modifiers involved in a series of post-translational modifications (PTMs), including acetylation, methylation, phosphorylation, ubiquitylation, and sumoylation [2, 3]. These modifications of histones play very important roles in gene regulation, genome stability, and genome defense in eukaryotes, mainly by altering the structure of chromatin and/or recruiting regulatory factors [2, 4,5,6]. Histone acetylation is one of the most studied PTM mechanisms, which alters the physical properties of nucleosomes by weakening interactions between histones and DNA [7]. Histone methylation and other PTMs often create binding sites for other proteins that are bound by specific effector proteins. These can either be involved in the repression of transcription by compacting nucleosome arrays, or they can support transcription by recruiting chromatin remodeling complexes, modifying enzymes, or other complexes involved in elongation or splicing [8,9,10,11]. Acetylation of histone proteins at the N-terminus lysine residues plays a crucial role in regulating gene activities in eukaryotes. Acetylation of core histones is correlated with “open” chromatin configurations and is gene activated, whereas deacetylation often produces “closed” chromatin configurations and is associated with gene repression [7, 12, 13]. Similar to most PTMs, histone acetylation is reversible and is dynamically modulated by the action of histone acetyltransferases (HATs) and histone deacetylases (HDACs). Generally, HATs cause gene activation by transferring acetyl groups (CH3COO-) from acetyl-CoA onto the lysine residues of the N-terminal tails of histone proteins [14, 15]. In contrast, HDACs regulate gene repression by removing acetyl groups from these acetylated lysine residues [14, 16]. The targets of HATs and HDACs include H3K9, H3K14, H3K36, H4K5, H4K8, H4K12, and H4K16 [4, 16]. All plant HATs are subdivided into four groups: the general control non-depressible GCN5-related acetyltransferase (GNAT) family; the MOZ, YBF2, SAS3, SAS2, and TIP60 (MYST) family; the camp-responsive element-binding protein (CREB)-binding protein (CBP) family; and the TATA-binding protein-associated factor 1 (TAFII 250) family [14, 17]. Corresponding HATs, which are represented by the acronyms HAG, HAM, HAC, and HAF, respectively [14], play essential roles in regulating gene expression during plant development, exogenous hormone response, and responses to environmental stresses [18,19,20,21,22,23,24]. The GNAT group includes three subfamilies: GCN5, elongated complex protein 3 (ELP3), and HAT1-like acetyltransferases, namely HAG1, HAG2, and HAG3 [14]. The GCN5 protein is the catalytic subunit of several multi-protein HAT complexes and plays an essential role in plant development and resistance to abiotic stressors, such as heat, drought, cold, salt, and phosphate starvation [13, 17, 25, 26]. The number of HAT gene family members varies among plants; 12 HATs have been identified in Arabidopsis [14], 8 in rice [27], and 32 in tomato [28]. Recently, Kumar et al. [3] provided a comprehensive review of the regulation of histone acetylation during plant growth, development, and stress response.

Foxtail millet (Setaria italica [L.] P. Beauv), which is one of the most important and ancient cultivated cereal crops, was domesticated in northern China approximately 11,500 years ago [29, 30]. This species exhibits specific morphological features, such as root architectural modifications, small leaf area, thick cell walls, and epidermal cell arrangement, that imparts stress tolerance and high water and nutrient utilization efficiency. These features are bolstered by a small diploid genome (approximately 430 Mb), a short life cycle, and C4 photosynthesis [31, 32]. These characteristics make foxtail millet an ideal model crop for exploring basic biological features, such as plant architecture, physiology, and genome evolution [32, 33]. Additionally, stress tolerance related characteristics of foxtail millet reduce its dependence on synthetic fertilizers, pesticides, herbicides, and insecticides, highlighting it as a model crop for exploring the mechanisms of stress tolerance. With the rapid development of molecular biology, the entire genome of foxtail millet has been sequenced and published by the United States Department of Energy Joint Genomic Institute and Beijing Genomics Institute of China [31, 34], enabling research on the mechanisms of stress response and molecular regulation and providing a foundation for genome-wide analysis of the HAT family members in this species.

In this study, we examined the whole genome of foxtail millet, using bioinformatics analysis, with a focus on the physicochemical properties, chromosomal localization, systematic evolution, gene structure, and conserved domain of the HAT gene family. We performed expression profiling in foxtail millet after exposure to abiotic stresses, such as drought, high salinity, low nitrate, and low phosphate, to further identify the function of the HAT genes. Additionally, we analyzed the response of SiHATs to pathogenic Sclerospora graminicola infection. Overall, our results provide a foundation for further study of the functions of HAT genes, particularly in the responses to abiotic stresses, and will pave the way for identifying the precise role of HATs in plant growth and development.


Identification and chromosome mapping of the foxtail millet HAT gene family

Extensive searches of public and proprietary transcripts and genomic databases, with all previously reported HAT proteins (containing GNAT, MYST, P300/CBP, and TAF1) of rice and Arabidopsis, were conducted. A total of 24 HATs were identified in foxtail millet from the Yugu1 genome after excluding redundant genes (Additional file 1). In addition, the position and direction of transcription of each gene were determined on foxtail millet chromosome pseudomolecules available on Phytozome (v12.1) (, as shown in Fig. 1. The 24 foxtail millet HAT genes were found to be distributed on nine chromosomes: eight on chromosome 2; three each on chromosomes 1, 4, and 5; two each on chromosomes 6 and 9; and one each on chromosomes 3, 7, and 8 (Additional file 1, Fig. 1).

Fig. 1
figure 1

Chromosome locations of histone acetylation genes (HATs) in Setaria italica. Chromosomal location was performed on 24 histone acetylation gene family members in S. italica

Additionally, we analyzed the physical and chemical properties of all HAT family genes and encoded proteins, including the number of amino acids, molecular weight (Mw), isoelectric point (pI), and subcellular location in foxtail millet. The sizes of the 24 predicted SiHAT proteins ranged from 425 aa (SiHAT9) to 5068 aa (SiHAT5), with molecular weights ranging from 38.08 (SiHAT16) to 563.55 kDa (SiHAT5). The pI values ranged from 4.96 (SiHAT20) to 9.86 (SiHAT24). SiHAT21, SiHAT17, and SiHAT5 were determined to be neutral proteins (–0.5 < index of GRAVY values < 0.5), whereas the GRAVY values of the remaining proteins were < 0, indicating hydrophilic properties. Subcellular localization showed that most of the HAT genes were located in the nucleus, while two were localized to the mitochondria (SiHAT6 and SiHAT24), and two were cytoplasmic (SiHAT11 and SiHAT21). Only one protein (SiHAT5) was localized to the endoplasmic reticulum. More detailed information, including sequence, aliphatic index, instability index, and subcellular localization, are listed in Additional file 1. Prediction of the secondary structure of SiHAT proteins indicated that every member contained α-helix, extended chain, β-folding, and irregular curl structures. The irregular curl and α-spiral structures were the main secondary components, accounting for 30–50% of the secondary structure, while β-folding accounted for only about 5% (Additional file 1).

A chromosome region containing more than two genes within 200 kb is defined as tandem duplication [35]. Homology analysis of SiHATs showed that there were two tandem duplication events in the foxtail millet chromosome sequences, each containing SiHAT9 and SiHAT10 on chromosome 2 and SiHAT17 and SiHAT18 on chromosome 5 (Fig. 1).

Phylogenetic analysis, motif composition, and structure analysis of SiHATs

Neighbor-joining phylogenetic analysis was performed, and the 24 SiHAT proteins were clustered into Groups I, II, III, and IV, with 12, 3, 5, and 4 members, respectively (Fig. 2a). Notably, most SiHAT proteins fell into sister pairs (SiHAT3 and SiHAT20, SiHAT9 and SiHAT22, SiHAT12 and SiHAT14), triplets (SiHAT8, SiHAT10, and SiHAT24) or quadruplets (SiHAT2, SiHAT15, SiHAT4, and SiHAT18) in the joint phylogenetic tree (Fig. 2a).

Fig. 2
figure 2

Maximum likelihood phylogenetic trees and structure and conserved domains of histone acetylation gene (HAT) in Setaria italica (Si). a Phylogenetic tree and subfamily of SiHATs, which are further divided into four groups. b The exon-intron organization of SiHATs. c Conserved domain of the HAT protein in foxtail millet. Bromo domains are conserved domain of GNAT subfamily; PHD, ZnF and ZZ are conserved domains of CBP subfamily, and TBP-binding domains are conserved domain of TAF subfamily

To obtain more insights into gene evolution, the exon-intron organization of SiHAT genes was investigated by aligning predicted coding sequences (CDS) against corresponding genomic sequences using the online service Gene Structure Display Server (GSDS). The number of introns in the SiHAT family ranged from 2 to 23. Overall, highly similar gene structures and domains were observed for the four HAT subfamilies. In contrast, SiHAT18 and SiHAT14 did not contain both upstream and downstream untranslated regulatory regions (UTR), and SiHAT24 did not contain upstream regulatory regions. The other 21 genes exhibited upstream and downstream regulatory regions (Fig. 2b, Additional file 2). Noticeably, the closest members from the same subgroups had highly similar intron/exon structure (intron number and exon length; Fig. 2b).

To further study the characteristic regions of SiHAT proteins, the motifs of 24 SiHAT proteins were analyzed using Multiple Expectation maximizations for Motif Elicitation (MEME). The results showed that 12 SiHAT genes from group I belonged to the GNAT family and possessed the bromodomain. (Fig. 2c). The genes in Group II contained the motifs of the CBP (ZnF_TAZ), MYST(PLN00104), and TAF1(Bromo_AAA) family, individually. While the quadruplets genes in Group III belonged to the CBP family and contained the typical motif of the HAT_KAT11 domain, PHD_SF domain, zf-TAZ, and ZZ domains. The five genes in Group IV belonged to the GNAT, CBP, and TAF families. Several HAT proteins possessed unique conserved domains, such as ELP3 in SiHAT21, HAT1 chromodomain, and Znf-C2H2 in the GNAT/MYST family. PHD (Plant Homeodomain), Znf-ZZ, and Znf-TAZ domains were observed in the CBP family. These four subfamily proteins were also compared with other typical GNAT/CBP/TAF/MYST conserved domains in Arabidopsis and O. sativa (Additional files 3, 4, 5 and 6). These conserved domains allowed SiHATs to interact with RNA Pol II during transcript elongation, bind with the transactivation domain of transcription factors and acetylated histone lysine residues, and interact with co-factors (Additional file 2). Interestingly, the sister pair genes also had the same structure and conserved domain, indicating they may have the same function in foxtail millet.

Phylogenetic relationship and collinearity analysis of HATs in Setaria italica, Oryza sativa, and Arabidopsis thaliana

To better understand the phylogeny of the foxtail millet HAT gene family, the SiHATs were subjected to synteny analysis with HAT genes of the typical model plants: the dicot Arabidopsis thaliana and monocot Oryza sativa. A total of 19 SiHAT genes were synchronized with those in O. sativa, thus, a phylogenetic tree was constructed using the protein sequences of all 24 SiHATs, 19 OsHATs, and 12 AtHATs. These 55 HAT proteins were divided into four clades (Fig. 3a). The results showed the phylogenetic relationship of HAT proteins between dicots and monocots. Apart from Clade II, which is the unique group of foxtail millet, containing 11 HAT gene family members, the other clades included HAT proteins from the three species, suggesting that these genes existed before the divergence of monocots and dicots. Clade I was further subdivided into 4 classes, namely a, b, c, and d, with 4, 3, 5, and 4 members respectively. Clade III was further divided into classes e, f, and g (Fig. 3a). Further, we analyzed the 11 SiHATs in Clade II and found they all belonged to the GNAT family, possessed a Bromo domain, and had structures obviously different from those of the other family members in other clades (Additional files 3 and 7).

Fig. 3
figure 3

Phylogenetic analysis and Collinearity analysis of histone acetylation proteins (HATs) in Arabidopsis thaliana (At), Oryza sativa (Os), and Setaria italica (Si). a Phylogenetic tree of HAT genes. b Collinearity of foxtail millet HAT and related species. The green rectangular color block represents the foxtail millet chromosome, yellow rectangular color block represents the rice chromosome and red rectangular block represents Arabidopsis chromosome. The number represents the chromosome number

The ratio of non-synonymous to synonymous (Ka/Ks) nucleotide substitutions was calculated to investigate the selection pressure on SiHATs [36]. The SiHAT genes of foxtail millet were subjected to one-to-one orthologous analysis with its homologous genes in Arabidopsis thaliana and Oryza sativa. (Fig. 3b, Additional file 8). In total, 19 SiHAT genes displayed a syntenic relationship with those in Oryza sativa, while there were two homologous pairings in Arabidopsis thaliana. The results indicated that the foxtail millet and rice HAT genes were genetically similar. We also found that some SiHAT genes were evolved from rice and Arabidopsis, respectively (Fig. 3b, Additional file 8).

Cis-elements analysis of SiHAT promoters

To further investigate the putative functions of SiHAT genes, a plant promoter database (PlantCARE) search was conducted in the promoter regions at 2000-bp upstream of the transcription initiation site of SiHAT genes. As shown in Fig. 4 and Additional file 9, three main categories of cis elements were found in the promoter sequences of SiHAT genes. The first category was involved in phytohormones, such as abscisic acid (ABA), methyl jasmonate (MeJA), auxin, and salicylic acid (SA). The second category was associated with stresses, such as anaerobic induction, drought inducibility, low-temperature responsiveness, pathogen infection, wound responsiveness, and salt inducibility. The last category was related to plant growth and development, such as zein metabolism regulation. Meristem, root, endosperm inducibility (GC-motif), and abscisic acid responsive elements were found in almost all gene promoters. Importantly, all the 24 SiHAT genes contained the light responsive element, while the MeJA-responsive element (TGACG-motif and CGTCA-motif), anoxic specific inducibility element (GC-motif), and the abscisic acid responsive element (ABRE) were found in almost all gene promoters. Interestingly, distinct differences in cis elements between the sister pair genes, including the GA and MeJA response elements, were found in the promoter of SiHAT3, whereas ABA and defense and stress response elements were found in SiHAT20. These results showed that SiHATs may have affected hormone signal responsiveness, stress adaptation, and development. No cytokinin-responsive elements were identified in these promoter regions.

Fig. 4
figure 4

Maximum likelihood phylogenetic trees and prediction of cis-acting elements in promoter of histone acetylation genes ( HAT s) in Setaria italica (Si). a Phylogenetic trees of foxtail millet HAT genes. b Promoter analysis of foxtail millet HAT genes. The 2 Kb promoter sequences of corresponding HAT genes were used to analyze hormone-related cis-elements, plant growth and development cis-elements and stress-related elements. Different cis-elements were indicated by different colored symbols and placed in their relative position on the promoter of SiHATs

Spatial and temporal expression of SiHAT genes

To obtain insight into the expression patterns of SiHAT genes in various tissues, a heat map was generated using the gene expression data in the foxtail millet Exp database. The results showed complex specific and overlapping SiHAT expression in various tissues and organs. The expression level of the same gene varied among tissues and organs; for example, SiHAT17 was highly expressed in top leaves 2–3 days after the heading stage. In contrast, low or no expression signals were detected in panicles. On the other hand, the expressions of different genes were also notably different in the same tissues and organs. For example, in the root at the filling stage, the expression of SiHAT3, SiHAT9, SiHAT13, and SiHAT22 from the GNAT gene family was significantly higher than that of other genes. Some genes were exclusively expressed in single tissues or organs; for example, SiHAT17 was expressed in leaf top after 2–3 days, and SiHAT16 expression was observed in immature ears. In addition, SiHAT3, SiHAT13, and SiHAT22 were highly expressed in all tissues at different developmental stages, whereas two SiHAT15 and SiHAT12 exhibited almost no expression in any of the tested tissues (Fig. 5). These results demonstrated that the expression patterns of SiHATs differed among tissues and were associated with plant growth and development.

Fig. 5
figure 5

Relative expression patterns of histone acetylation genes (HATs) in different tissues of Setaria italica. Heat maps reflect the fragments per kilobase of transcript per million mapped fragments (FPKM) of HATs. Color from red to blue indicates high to low expression

Expression analysis of SiHATs under stress

To confirm whether the expression of SiHAT genes could be regulated by abiotic and biotic stress, we tested the effects of abiotic stresses such as nitrate deficiency and phosphate deficiency, salt-alkali, and drought.

Under low nitrate conditions, the expressions of most SiHAT genes were either slightly upregulated or downregulated. Most genes were downregulated after 2 h, which was the inverse of what was seen at 24 h. The expression of several genes (SiHAT17, SiHAT8, and SiHAT5) were clearly different between shoots and roots. SiHAT3 expression was continuously upregulated under low nitrate conditions in the shoot, but was only upregulated at 2 h in the root (Fig. 6). This suggests that SiHAT3 likely performs different control functions during nitrate absorption and transport.

Fig. 6
figure 6

Expression analysis of histone acetylation genes (HATs) in Setaria italica under low nitrogen stress. Low nitrate stress time represented by 10 min, 30 min, 2, 8, 24, and 72 h, NN and LN represent normal and low nitrogen treatment, respectively. Leaf and Root describe tissues sampled

Most SiHAT genes were not strongly upregulated under low phosphate conditions, while five SiHAT genes (SiHAT1, SiHAT6, SiHAT7, SiHAT19, and SiHAT21) were strongly upregulated in roots and weakly upregulated in shoots. In contrast, SiHAT17 was highly expressed in the shoot, and its expression pattern showed a sharp decline initially, then a gradual increase until returning to its original level at 24 h (Fig. 7).

Fig. 7
figure 7

Expression analysis of histone acetylation genes (HATs) in Setaria italica under low phosphorus stress. Low phosphorus stress time represented by 0.5, 2, 6, 12, 24, and 72 h. NP and LP represent normal and low phosphorus treatments, respectively. Leaf and Root describe tissues sampled

Previous studies have reported that the expression of HAT genes was induced by drought [37]. Therefore, we analyzed the expression of the 24 SiHATs in our RNA-seq data. The results showed that except for SiHAT15, SiHAT3, SiHAT21, and SiHAT18, the other SiHAT genes responded to drought stress with different expression patterns. Nine genes (SiHAT1, SiHAT5, SiHAT7, SiHAT8, SiHAT9, SiHAT10, SiHAT13, SiHAT19, and SiHAT22) were upregulated, while SiHAT17 and SiHAT15 were downregulated under drought conditions. SiHAT9 was highly expressed in drought-sensitive foxtail millet varieties. In response to circadian and drought treatments, SiHAT3, SiHAT9, SiHAT13, and SiHAT20 showed higher expression levels under dark conditions; meanwhile, the remaining five genes (SiHAT2, SiHAT8, SiHAT11, SiHAT23, and SiHAT24) showed lower expression (Fig. 8).

Fig. 8
figure 8

Expression analysis of HATs genes in foxtail millet under drought stress. Control and Drought represented Control and Drought treatment groups, respectively. R (AN04) and S (Yugu1) represent drought-resistant and drought-sensitive varieties, respectively. Light and Dark represent different sampling times and illumination: (light) illumination, (dark) no illumination

Under salt and alkali stress, SiHAT6 was upregulated at low levels, while the other genes were significantly downregulated. The transcript levels of most genes were lower at the germination stage than at the two-leaf one-heart stage. Only SiHAT6 and SiHAT19 showed high expression in the germination stage. Additionally, sensitive and resistant varieties showed differences in gene expression. At the T2 stage, the expression of 9 genes (SiHAT2, SiHAT3, SiHAT4, SiHAT5, SiHAT8, SiHAT10, SiHAT14, SiHAT18, and SiHAT20) was higher in the salt-resistant variety than in the sensitive variety, while the expression of three genes (SiHAT13, SiHAT17, and SiHAT19) showed the opposite trend (Fig. 9). SiHAT12 and SiHAT15 were the only two genes not or minimally expressed under all conditions (Figs. 6, 7, 8 and 9).

Fig. 9
figure 9

Expression analysis of HATs genes in foxtail millet under salt-alkali stress. CK and SAS represent control group and salt-alkali treatment group, respectively. R (B103) and S (B355) represent salt-resistant and salt-sensitive varieties, respectively. T1 and T2 represent different sampling tissues: T1 is Seedlings germinating for 3 days and T2 is one-tip-two-leaf Seedlings

SiHATs involved in Sclerospora graminicola infection

The oomycete S. graminicola (Sacc.) causes 20–30% of downey mildew cases in foxtail millet cultivated in China and results in the deterioration of yield and quality [37]. It is also prevalent in India, Japan, and Russia [37, 38]. We investigated the expression of SiHATs in our transcriptome sequencing data obtained after S. graminicola infection. Four SiHAT genes were detected in response to the infection, and their expression patterns were different. During the three-leaf-one-heart stage, SiHAT16 and SiHAT24 expressions were downregulated in the pathogen-resistant variety but upregulated in the sensitive variety, while SiHAT6 and SiHAT17 expressions were upregulated in both the resistant and sensitive varieties. At the five-leaf-one-heart stage, SiHAT6 and SiHAT24 expressions were upregulated after infection, and there was no difference between the sensitive and resistant varieties. At the seven-leaf-one-heart stage, the expression of all four genes were downregulated in the pathogen-resistant variety; however, SiHAT16 and SiHAT17 expressions were upregulated in the sensitive variety. The other genes were not or minimally expressed (Fig. 10).

Fig. 10
figure 10

Expression analysis histone acetylation genes (HATs) in Setaria italica under Sclerospora graminicola infection stress. CK and T represent the control and infection treatment groups, respectively. R and S represent resistant and sensitive varieties, respectively. Numbers 3, 5, and 7 represent leaf stages at different sampling times


Characterization of the expanded HAT gene family in foxtail millet

Multiple members of a specific gene family in a particular organism are the natural products generated of a long evolutionary history of that organism [39]. Key chromosome reshuffling events have occurred between foxtail millet, rice, and sorghum [31]. In this study, we identified and characterized 24 foxtail millet HAT genes using genome-wide analysis. We found that the foxtail millet HAT family is larger than that of Arabidopsis [14] and rice with 12 and 19 members, respectively [14, 27], but is smaller than that of wheat with 31 members [40]. The phylogenetic analysis shows that three sister pairs, one triplet, and one quadruplet were identified within the SiHAT family. However, none of these pairs were genetically linked to each other, as compared to their corresponding chromosomal locations. Conversely, all closely linked SiHAT loci, such as SiHAT9 and SiHAT10 on chromosome 2, and SiHAT17 and SiHAT18 on chromosome 5, were not paired together into the same sister groups. Moreover, there were no sister pairs mapped on the same duplicated chromosomal blocks (Additional file 1), as described previously. Domains and motifs have been shown to be involved in various activities, including protein interaction, transcriptional activity, and DNA binding; in this study, we found that homologous genes, such as sister pairs and quadruplets, contained conserved motifs. We also found that most SiHAT genes shared a similar exon/intron structure within the same phylogenetic group, although some differences were also observed (Fig. 3b). Gene structure analysis can provide important information about gene function and evolution. It has been revealed that intron gain or loss in plants is the result of selection pressures during evolution, and genes tend to evolve into diverse exon-intron structures and perform distinct functions [41]. Therefore, our results suggest that gene differentiation might have occurred in the SiHAT family to accomplish different biological functions under selection pressure during foxtail millet genome formation and evolution. Interestingly, clade II is composed of eleven foxtail millet HATs, all belonging to the GNAT subfamily, indicating expanding HAT gene events, particularly in the GNAT subgroup of foxtail millet. Genes in this group showed tissue-specific expression, indicating that it may be involved in important processes of growth and development in foxtail millet.

Expression divergence between duplicated SiHAT genes

The presence of duplicated SiHATs raises questions about their functional redundancy. According to evolutionary models, duplicated genes may undergo different selection processes: nonfunctionalization, where one copy loses function; hypofunctionalization, where one copy decreases in expression or function; neofunctionalization, where one copy gains a novel function; or subfunctionalization, where the two copies partition or specialize into distinct functions [42,43,44]. These evolutionary fates may be indicated by the divergence in expression patterns or protein structure. All SiHATs in this study appeared to be functional because they contained a credible and complete open reading frame, and their corresponding cDNA/ESTs are available in the NCBI database. Evidence for the divergence between the duplicate genes can be inferred from the expression pattern for the HAT quadruplet gene set. SiHAT2 was highly expressed in leaves and stems; however, the transcript levels of SiHAT4, SiHAT15, and SiHAT18 were very low. In addition, possible subfunctionalization trends were clear in the expression pattern shifts of gene pairs. For example, the mRNA abundance of SiHAT3 peaked in the root at the filling stage; however, SiHAT20 was highly expressed in the neck panicle internode. Under nitrogen deficiency, phosphorus deficiency, and salt and alkali stress, the expression patterns of SiHAT3 and SiHAT20 also differed possibly due to the large differences between the cis elements in their promoter.

Regulation of SiHAT gene expression

Accumulating evidence has shown that HAT family genes play pivotal roles in regulating plant growth and development. Arabidopsis AtHAG1/GCN5 is one of the most widely studied and functionally characterized acetyltransferases and regulates cell differentiation, leaf and floral meristem patterning, and plant defense pathways [45]. In foxtail millet, two SiHATs, SiHAT17 and SiHAT23, were grouped into the same cluster as AtHAG1/GCN5. Of these, SiHAT17 had a similar expression pattern to AtHAG1, which is involved in low phosphorus stress, drought, and salt and alkali stress responses. Many identified HATs are involved in seed development; in this study, we found that SiHAT23 was highly expressed in germinated seeds and contained the endosperm expressive element (GCN4_motif) in its promoter. This indicates the potential role of HAT in seed development and dormancy. A review by Nguyen et al. [20] demonstrated that HAT can collaborate with ARFs to trigger auxin response and regulate root growth and development. Additionally, some SiHATs (SiHAT3, SiHAT9, SiHAT13, and SiHAT22) are highly expressed in roots and may play important roles in root development.

Histone acetylation and abiotic stress responses

The relative expression levels of HAT genes change significantly under various biotic and abiotic stresses [16]. GCN5, which is involved in salt stress response, was first characterized in maize roots, where its expression increased under salt stress conditions [46]. Similar patterns of expression were detected in Arabidopsis under salt stress conditions [47], and the expression of GCN5 homologous genes SiHAT17 and SiHAT23 were also upregulated under salt stress in our experiment. OsHAC701 has been reported to respond to salt treatment in rice [27], and its homologs, SiHAT4 and SiHAT18, were also differentially expressed under salt and alkali stress in this study, especially at the seeding of one-heart-two-leaf. This suggests that it may play an important role in salt stress response in foxtail millet. In rice, four histone acetyltransferases, OsHAC703, OsHAG703, OsHAF701, and OsHAM701, were involved in drought responses by hyperacetylating lysine residues evident by their upregulation [48]. In our study, their closest homologs also responded to the drought stress, suggesting that the two likely perform the same functions against drought stress in foxtail millet.

Histone acetylation in nitrogen and phosphorous stress responses

Nitrogen and phosphorus are two major mineral nutrients required for plant growth and development. Histone modifications have been reported in responses to nitrogen levels. AtHNI9, which encodes a component of the RNA polymerase II complex (IWS1), plays an important role in modulating the low nitrogen response in Arabidopsis [49]. ZmCHB101 encodes a core subunit of the SWI/SNF-type ATP-dependent chromatin remodeling complex, which can bind to the promoter of the nitrogen transporter ZmNRT2.1/2.2 through nitrate-responsive cis-elements and regulate nitrogen transport. Additionally, its RNAi lines display accelerated root growth and increased biomass under low nitrogen conditions compared to that in wild-type plants [50]. This suggests that histone modification plays an important role in regulating plant growth and development under low nitrogen conditions. In foxtail millet, we found that some HAT genes were responsive to low nitrogen stress, and the expression patterns differed between the different groups, indicating that the mechanism of histone regulation is complex, and further investigation is needed to clarify the function of SiHATs under low nitrate responses.

The histone acetyltransferase AtGCN5 is required for the activation of several genes under low Pi conditions, including At4 and AtWRKY6. Mutation of AtGCN5 reduces Pi concentration in Arabidopsis plant and impairs Pi accumulation between roots and shoots. This indicates that GCN5-mediated histone acetylation regulates the phosphate starvation response through the At4-miR399-PHO2 pathway [51]. We found that the GCN5 homolog gene SiHAT17 responded to low phosphorus stress and was highly expressed in roots under low phosphorus conditions. However, the expression pattern was reversed in the shoot, suggesting a difference in histone modification between shoot and root under low phosphate conditions in foxtail millet. These findings prove that the mechanisms regulating the responses of shoots and roots to low phosphate stress are different, and further investigation should be conducted on the role of SiHAT17 in the regulation of low phosphorus stress.

Histone acetylation and biotic stress responses

Histone acetylation dynamics by HATs and HDACs is a key regulatory epigenetic mechanism that ultimately regulates plant development, hormone homeostasis, stress response, and defense responses in plants [3, 52,53,54,55]. In this study, we identified 24 HATs in foxtail millet and found more abiotic and biotic responsive elements in their promoters; and most SiHAT genes responded to nitrogen deficiency, phosphorus deficiency, drought, salt alkali stress, and S. graminicola infection. All these favorable conditions indicate that these proteins may be involved in plant responses to abiotic and biotic stress; however, the mechanisms behind this stress response, and how HATs interact with other histone modifiers, remain unclear. Future studies should include transgenic research and protein interaction studies on key genes related to abiotic and biotic stress, and the use of new CRISPR-Cas9 technology.


Recent studies show that histone acetylation and deacetylation play essential roles in the regulation of plant growth and development, as well as in responses to stress. Analyses of various mutations of HAT and HDAC genes in Arabidopsis have revealed the function of histone acetylation/deacetylation in plant development. Here, we screened 24 HAT genes that may be involved in abiotic and biotic stress responses in foxtail millet. We identified their chromosomal locations, protein structures, gene duplications, promoters, and conserved motifs. Phylogenetic and synteny comparisons between AtHATs, OsHATs, and SiHATs were performed, and the potential roles of SiHATs in the development and growth of foxtail millet were investigated based on our previously published RNA-seq data. Candidate SiHATs involved in responses to low nitrate, low phosphate, drought, salt response, and S. graminicola infection were examined and assessed through the analysis of promoter elements. Our results lay a foundation for research regarding histone modification mediated stress response regulation in foxtail millet and similar crops. Further studies, including biological experiments, will be required to confirm the functions of the candidate genes, identify the interacting HATs and HDACs, and clarify the molecular mechanisms by which histone acetylation/deacetylation affects various biological processes. This study serves as a reference to improve the stress resistance of plants.


Sequence retrieval and identification of HAT genes

We retrieved data containing sequence IDs, protein sequences, genomic sequences, and conserved domain database (CDD) of foxtail millet and rice from Phytozome ( V12.1). The AtHAT gene sequence from Arabidopsis thaliana was retrieved from Uniprot (, and the Hidden Markov Model was used to identify foxtail millet genes in a protein database with the BLASTP program (p value = 0.001). The obtained proteins were run through the Pfam database and SMART ( to eliminate the sequences not containing complete HAT domains. Then, the conserved domains were checked by the conserved domain database (CDD) program ( in order to confirm the presence of the complete HAT domain. Then, all candidate sequences were verified using BLASTn ( and HMMER ( Finally, the cis-acting elements were predicted using PlantCARE (

Phylogenetic analysis of HATs

The exon/intron structures of the SiHATs were deduced from alignments of cDNA and BAC genomic sequences using the gene structure displayer ( Multiple-sequence alignments of SiHAT proteins were carried out using the Clustal W (version 2.0) program [56]. The protein sequences of histone acetyltransferase proteins in Arabidopsis and rice were obtained from the TIGR database. Then, phylogenetic analysis was performed with MEGA7.0 [57] using the neighbor-joining method (1000 bootstrap replications), and the results were visualized in Itol [58].

Swiss-Prot, physicochemical properties, subcellular localization, pI, and Mw of the putative SiHATs were calculated using the ExPASy online tool ( Secondary protein structures were predicted using SOPMA software (

Chromosome information was obtained from Phytozome (, and a chromosomal location map of the genes was generated by TBtools ( A distinctive name was given to each SiHAT according to its initial position on the chromosomes. Conserved motif analysis of foxtail millet HAT protein sequences was conducted using the MEME suite 4.11.1 software ( [59], with the motif width set to 6–300 and the maximum number of motifs set to 20. The results were visualized using TBtools.

Analysis of duplication and synteny of HAT family genes in Arabidopsis, foxtail millet, and rice

To confirm the gene duplication events of HATs, we investigated the ancient duplication events between A. thaliana, O. sativa, and S. italica. The non-synonymous substitution rate (Ka), synonymous substitution rate (Ks), and Ka/Ks ratio for each pair of duplicated genes among A. thaliana, O. sativa, and S. italica were computed between pairs of homologous genes using the Ka/Ks calculator in TBtools with default settings. Multiple collinear scanning was used to simultaneously detect the homologous genetic relationships between different species.

Promoter analysis of SiHAT genes

The cis-acting elements of promoters are essential for determining tissue-specific expression and are involved in the regulation of gene expression under abiotic stress. A region approximately 2000 bp upstream of the start codon (ATG) was investigated from the reference genome sequence of foxtail millet. The online software PlantCARE ( was used to search for cis-acting regulatory elements in the promoters of SiHATs.

Analysis of gene expression profiles

Data on the expression of SiHAT genes in 23 tissue types of the foxtail millet cultivar JG21 were obtained from a published foxtail millet database ( [32]. These tissues included leaves at the seeding and different filling stages, seeds at various maturation stages, roots at the filling stage, panicles, and neck panicles, as detailed in Additional file 10. Detailed information can be found in Additional file 9. All RNA sequencing library preparations were conducted as described in previous study [32]. The time and spatial expression data for SiHATs were collected from our Electronic Fluorescent Pictograph (xEFP) program ( and the aforementioned database [32]. Heatmap Illustrator HemiI v.1.0, was used to plot the heat map of gene expression [60].

Tissue preparation

Seeds from the foxtail millet cultivar Jingu21 were germinated on paper rolls for three days, then incubated at 25 ℃ (day) or 22 ℃ (night) on a 16: 8 h light: dark schedule for three weeks. Seedlings at the five-leaf stage were transferred to Hoagland solution with 1 mM KH2PO4 (Control, CK) or 5 µM KH2PO4 (low phosphorus treatment). The shoot and root tissues were harvested at 0.5, 2, 6, 12, 24, and 72 h after treatment.

For the nitrate experiment, the plants were prepared as before. The nitrate concentration was 2.0 mM in the control treatment, (normal nitrate; NN), and 0.2 mM in low nitrate (LN). After 10 min, 30 min, and 2, 8, 24, and 72 h, the shoot and root tissue were harvested. Special care was taken to characterize the materials used, and three replicates were harvested for RNA extraction.

The seedlings were grown in distilled water for 14 d, then subjected to remaining stress treatments. For the drought treatment, two different varieties of foxtail millet were used: AN04 (drought resistant) and Yugu 1 (drought sensitive). Conditions included different circadian stages as morning (8:00–9:00), noon (12:00–13:00), evening (18:00–19:00). Leaf tissues from plants under the control and drought conditions were collected in the morning, noon, and evening after three days of treatment.

For the salt and saline treatments, salt-sensitive (B355) and salt-tolerant (B103) seeds were germinated in Hoagland nutrient solution supplemented with 200 mM NaCl and a 30 mM mixture of Na2CO3 and NaHCO3 to represent salt and saline-alkaline stress; they were germinated in distilled water as the control treatment. Plants were harvested after 3 d (at the germination stage) and two weeks (two-leaf one-heart stage).

For the S. graminicola infection experiments, Jingu21 (pathogen-sensitive) and Jingu42 (pathogen-resistant) cultivars were germinated as described for the phosphate experiments. They were then infected at three growth periods: the 3-leaf, 5-leaf, and 7-leaf stages. The uninfected plant was set as the control, and three replicates were harvested for each genotype between 8:00 and 9:00 am.

Availability of data and materials

All datasets supporting the results of this study are included within the article and its supplementary information.


  1. Luger K, Mäder AW, Richmond RK, Sargent DF, Richmond TJ. Crystal structure of the nucleosome core particle at 2.8 Å resolution. Nature. 1997;389:251–60.

    CAS  PubMed  Article  Google Scholar 

  2. Kouzarides T. Chromatin modifications and their function. Cell. 2007;128:693–705.

    CAS  PubMed  Article  Google Scholar 

  3. Kumar V, Thakur JK, Prasad M. Histone acetylation dynamics regulating plant development and stress responses. Cell Mol Life Sci. 2021;78:4467–86.

    CAS  PubMed  Article  Google Scholar 

  4. Zhang K, Sridhar VV, Zhu J, Kapoor A, Zhu JK. Distinctive core histone post-translational modification patterns in Arabidopsis thaliana. PLoS One. 2007;11:e1210.

    Article  CAS  Google Scholar 

  5. Chinnusamy V, Zhu JK. Epigenetic regulation of stress responses in plants. Curr Opin Plant Biol. 2009;12:133–39.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  6. Portela A, Esteller M. Epigenetic modifications and human disease. Nat Biotechnol. 2010;28:1057–68.

    CAS  PubMed  Article  Google Scholar 

  7. Grunstein M. Histone acetylation in chromatin structure and transcription. Nature. 1997;389:349–52.

    CAS  PubMed  Article  Google Scholar 

  8. Lin G, Zhou Y, Li M, Fang Y. Histone 3 lysine 36 to methionine mutations stably interact with and sequester SDG8 in Arabidopsis thaliana. Sci China Life Sci. 2018;61:225–34.

    CAS  PubMed  Article  Google Scholar 

  9. Liu R, Li X, Chen W, Du J. Structure and mechanism of plant histone mark readers. Sci China Life Sci. 2018;61:170–7.

    CAS  PubMed  Article  Google Scholar 

  10. Qian S, Lv X, Scheid RN, Lu L, Yang Z, Chen W, et al. Dual recognition of H3K4me3 and H3K27me3 by a plant histone reader SHL. Nat Commun. 2018;9:2425.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  11. Yang Z, Qian S, Scheid RN, Lu L, Chen X, Liu R, et al. EBS is a bivalent histone reader that regulates floral phase transition in Arabidopsis. Nat Genet. 2018;50:1247–53.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  12. Shahbazian MD, Grunstein M. Functions of site-specific histone acetylation and deacetylation. Annu Rev Biochem. 2007;76:75–100.

    CAS  PubMed  Article  Google Scholar 

  13. Yuan LY, Liu XC, Luo M, Yang SG, Wu KQ. Involvement of histone modification in plant abiotic stress responses. J Integr Plant Biol. 2013;55:892–901.

    CAS  PubMed  Article  Google Scholar 

  14. Pandey R, Müller A, Napoli CA, Selinger DA, Pikaard CS, Richards EJ, et al. 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–55.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  15. Aquea F, Timmermann T, Arce-Johnson P. Analysis of histone acetyltransferase and deacetylase families of Vitis vinifera. Plant Physiol Biochem. 2010;48:194–9.

    CAS  PubMed  Article  Google Scholar 

  16. Bjerling P, Silverstein RA, Thon G, Caudy A, Grewal S, Ekwall K. Functional divergence between histone deacetylases in fission yeast by distinct cellular localization and in vivo specificity. Mol Cell Biol. 2002;22:2170–81.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  17. Perrella G, Consiglio MF, Aiese-Cigliano R, Cremona G, Sanchez-Moran E, Barra L, et al. Histone hyperacetylation affects meiotic recombination and chromosome segregation in Arabidopsis. Plant J. 2010;62:796–806.

    CAS  PubMed  Article  Google Scholar 

  18. Liu X, Yang S, Yu CW, Chen CY, Wu K. Histone acetylation and plant development. Enzymes. 2016;40:173–99.

    CAS  PubMed  Article  Google Scholar 

  19. Papaefthimiou D, Likotrafiti E, Kapazoglou A, Bladenopoulos K, Tsaftaris A. Epigenetic chromatin modifiers in barley: III. Isolation and characterization of the barley GNAT-MYST family of histone acetyltransferases and responses to exogenous ABA. Plant Physiol Biochem. 2010;48:98–107.

    CAS  PubMed  Article  Google Scholar 

  20. Nguyen CT, Tran GB, Nguyen NH. Homeostasis of histone acetylation is critical for auxin signaling and root morphogenesis. Plant Mol Biol. 2020;103:1–7.

    CAS  PubMed  Article  Google Scholar 

  21. Wang Z, Cao H, Sun Y, Li X, Chen F, Carles A, et al. Arabidopsis paired amphipathic helix proteins SNL1 and SNL2 redundantly regulate primary seed dormancy via abscisic acid-ethylene antagonism mediated by histone deacetylation. Plant Cell. 2013;25:149–66.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  22. Kaldis A, Tsementzi D, Tanriverdi O, Vlachonasios KE. Arabidopsis thaliana transcriptional co-activators ADA2b and SGF29a are implicated in salt stress responses. Planta. 2011;233:749–62.

    CAS  PubMed  Article  Google Scholar 

  23. Zhou Y, Tan B, Luo M, Li Y, Liu C, Chen C, et al. Histone deacetylase19 interacts with HSL1 and participates in the repression of seed maturation genes in Arabidopsis seedlings. Plant Cell. 2013;25:134–48.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  24. Li C, Liu D, Lin Z, Guan B, Liu D, Yang L, et al. Histone acetylation modification affects cell wall degradation and aerenchyma formation in wheat seminal roots under waterlogging. Plant Growth Regul. 2019;87:149–63.

    Article  CAS  Google Scholar 

  25. Hu Z, Song N, Zheng M, Liu X, Liu Z, Xing J, et al. Histone acetyltransferase GCN5 is essential for heat stress-responsive gene activation and thermotolerance in Arabidopsis. Plant J. 2015;84:1178–91.

    CAS  PubMed  Article  Google Scholar 

  26. Mao YP, Pavangadkar KA, Thomashow MF, Triezenberg SJ. Physical and functional interactions of Arabidopsis ADA2 transcriptional coactivator proteins with the acetyltransferase GCN5 and with the cold-induced transcription factor CBF1. Biochim Biophys Acta. 2006:1759: 69–79.

    CAS  PubMed  Article  Google Scholar 

  27. Liu X, Luo M, Zhang W, Zhao JH, Zhang JX, Wu KQ, et al. Histone acetyltransferases in rice (Oryza sativa L.): phylogenetic analysis, subcellular localization and expression. BMC Plant Biol. 2012;12:145.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. Cigliano RA, 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 Genom. 2013;14:57.

    CAS  Article  Google Scholar 

  29. Hu H, Mauro-Herrera M, Doust AN. Domestication and improvement in the model C4 grass Setaria. Front Plant Sci. 2018;9:719.

    PubMed  PubMed Central  Article  Google Scholar 

  30. Lu H, Zhang J, Liu KB, Wu N, Li Y, Zhou K, et al. Earliest domestication of common millet (Panicum miliaceum) in East Asia extended to 10,000 years ago. Proc Natl Acad Sci U S A. 2009;106: 7367–72.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  31. Zhang GY, Liu X, Quan ZW, Cheng SF, Xu X, Pan SK, et al. Genome sequence of foxtail millet (Setaria italica) provides insights into grass evolution and biofuel potential. Nat Biotechnol. 2012;30:549–54.

    CAS  PubMed  Article  Google Scholar 

  32. Yang Z, Zhang H, Li X, Shen H, Gao J, Hou S, et al. A mini foxtail millet with an Arabidopsis-like life cycle as a C4 model system. Nat Plants. 2020;6:1167–78.

    CAS  PubMed  Article  Google Scholar 

  33. Doust AN, Kellogg EA, Devos KM, Bennetzen JL. Foxtail millet: a sequence-driven grass model system. Plant Physiol. 2009;149:137–41.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. Bennetzen JL, Schmutz J, Wang H, Percifield R, Hawkins J, Pontaroli AC, et al. Reference genome sequence of the model plant Setaria. Nat Biotechnol. 2012;30:555–61.

    CAS  PubMed  Article  Google Scholar 

  35. Holub EB. The arms race is ancient history in Arabidopsis, the wildflower. Nature Rev Genet. 2001;12:516–27.

    Article  CAS  Google Scholar 

  36. Nekrutenko A, Makova KD, Li WH. The Ka/Ks ratio test for assessing the protein-coding potential of genomic regions: An empirical and simulation study. Genome Res. 2002;12:198–202.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  37. Kobayashi M, Hiraka Y, Abe A, Yaegashi H, Natsume S, Kikuchi H, et al. Genome analysis of the foxtail millet pathogen Sclerospora graminicola reveals the complex effector repertoire of graminicolous downy mildews. BMC Genom. 2017;18:897.

    Article  CAS  Google Scholar 

  38. Li RJ, Han YQ, Zhang Q, Chang GR, Han YH, Li XK, et al. Transcriptome profiling analysis reveals co-regulation of hormone pathways in foxtail millet during Sclerospora graminicola infection. Int J Mol Sci. 2020;21:1226.

    CAS  PubMed Central  Article  Google Scholar 

  39. Xing HY, Pudake RN, Guo GG, Xing GF, Hu ZR, Zhang YR, et al. Genome-wide identification and expression profiling of Auxin Response Factor (ARF) gene family in maize. BMC Genom, 2011,12:178

    CAS  Article  Google Scholar 

  40. Gao SQ, Li LZ, Han XL, Liu TT, Jin P, Cai LN, et al. Genome-wide identification of the histone acetyltransferase gene family in Triticum aestivum. BMC Genom. 2021;22:49.

    CAS  Article  Google Scholar 

  41. Wang M, Yue H, Feng KW, Deng PC, Song WN, Nie XJ. Genome-wide identification, phylogeny and expressional profiles of mitrogen activated protein kinase kinase kinase (MAPKKK) gene family in bread wheat (Triticum aestivum L.). BMC Genom. 2016;17:668.

    Article  CAS  Google Scholar 

  42. Lynch M, Conery JS. The evolutionary fate and consequences of duplicate genes. Science. 2000;290:1151–5

    CAS  PubMed  Article  Google Scholar 

  43. Otto SP, Yong P. The evolution of gene duplicates. Adv Genet. 2002;46:451–83.

    CAS  PubMed  Article  Google Scholar 

  44. Duarte JM, Cui L, Wall PK, Zhang Q, Zhang X, Leebens-Mack J, et al. Expression pattern shifts following duplication indicative of subfunctionalization and neofunctionalization in regulatory genes of Arabidopsis. Mol Biol Evol. 2006;23:469–78.

    CAS  PubMed  Article  Google Scholar 

  45. Servet C, e Silva NC, Zhou DX. Histone acetyltransferase AtGCN5/HAG1 is a versatile regulator of developmental and inducible gene expression in Arabidopsis. Mol Plant. 2010;3:670–7.

    CAS  PubMed  Article  Google Scholar 

  46. Li H, Yan S, Zhao L, Tan J, Zhang Q, Gao F, et al. Histone acetylation associated up-regulation of the cell wall related genes is involved in salt stress induced maize root swelling. BMC Plant Biol. 2014;14:1–14.

    Article  Google Scholar 

  47. Zheng M, Liu XB, Lin JC, Liu XY, Wang Z, Xin M, et al. Histone acetyltransferase GCN5 contributes to cell wall integrity and salt stress tolerance by altering the expression of cellulose synthesis genes. Plant J. 2019;97:587–602.

    CAS  PubMed  Google Scholar 

  48. Fang H, Liu X, Thorn G, Duan J, Tian L. Expression analysis of histone acetyltransferases in rice under drought stress. Biochem Biophys Res Commun. 2014;443:400–5.

    CAS  PubMed  Article  Google Scholar 

  49. Widiez T, El Kafafi S, Girin T, Berr A, Ruffel S, Krouk G, et al. High nitrogen insensitive 9 (HNI9)-mediated systemic repression of root NO3 uptake is associated with changes in histone methylation. Proc Natl Acad Sci U S A. 2011;108:13329–34.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  50. Meng X, Yu X, Wu Y, Kim DH, Nan N, Cong W, et al. Chromatin remodeling protein ZmCHB101 regulates nitrate-responsive gene expression in maize. Front Plant Sci. 2020;11:52.

    PubMed  PubMed Central  Article  Google Scholar 

  51. Wang T, Xing J, Liu Z, Zheng M, Yao Y, Hu Z, et al. Histone acetyltransferase GCN5-mediated regulation of long non-coding RNA At4 contributes to phosphate starvation response in Arabidopsis. J Exp Bot. 2019;70:6337–48.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  52. Kim S, Piquerez SJM, Ramirez-Prado JS, Mastorakis E, Veluchamy A, Latrasse D, et al. GCN5 modulates salicylic acid homeostasis by regulating H3K14ac levels at the 5′ and 3′ ends of its target genes. Nucleic Acids Res. 2020;48:5953–66.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  53. Wang T, Xing J, Liu X, Yao Y, Hu Z, Peng H, et al. GCN5 contributes to stem cuticular wax biosynthesis by histone acetylation of CER3 in Arabidopsis. J Exp Bot. 2018;69:291122.

    Google Scholar 

  54. Kong L, Zhi P, Liu J, Li H, Zhang X, Xu J, et al. Epigenetic activation of Enoyl CoA reductase by an acetyltransferase complex triggers wheat wax biosynthesis. Plant Physiol. 2020;183:1250–67.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  55. Zhi P, Kong L, Liu J, Zhang X, Wang X, Li H, et al. Histone deacetylase TaHDT701 functions in TaHDA6-TaHOS15 complex to regulate wheat defense responses to Blumeria graminis f.sp. tritici. Int J Mol Sci. 2020;21:2640.

    CAS  PubMed Central  Article  Google Scholar 

  56. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–48.

    CAS  PubMed  Article  Google Scholar 

  57. Kumar S, Stecher G, Tamura K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Molec Biol Evol. 2016;33:1870–74.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  58. Letunic I, Bork P. Interactive tree of life (iTOL) v4: recent updated and new developments. Nucleic Acids Res. 2019;47: 256–259.

    Article  CAS  Google Scholar 

  59. Bailey TL, Williams N, Misleh C, Li WW. MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006;34:369–73.

    Article  CAS  Google Scholar 

  60. Deng WK, Wang YB, Liu ZX, Cheng H, Xue Y. HemI: a toolkit for illustrating heatmaps. Plos One. 2014;9:e111988.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

Download references


We would like to thank Doctor Siyu Hou and Li Zhang, University of Shanxi Agricultural University, China, for discussions and help with the manuscript.


This work was supported by Shanxi Scholarship Council of China (Grant No. 2021-071); the Natural Science Foundation of Shanxi Province (Grant No. 20210302123364); National Science Foundation of China (Grant No.31901598, 32070366).

Author information

Authors and Affiliations



MJ and RQ performed the bioinformatics analyses; GX drafted the manuscript; XZ, YH, and FM analyzed the transcriptomic data; JZ, YH, and MJ planted the crops and collected materials; XZ, XL, and XW supervised the research and modified the manuscript. GX and MJ contributed equally. All authors have read and approved the final manuscript.

Corresponding authors

Correspondence to Fangfang Ma or Xiongwei Zhao.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

 Table S1. Basic physicochemical properties, Secondary structure prediction and subcellular localization of HATs in Setaria italica.

Additional file 2: 

Table S2. Information of histone acetylation gene (HAT) structure and conserved domain in foxtail millet.

Additional file 3: Fig. S1. 

Phylogenetic trees and domain composition of GNAT subfamily. Phylogenetic tree and domain composition of GNAT subfamily predicted proteins from Arabidopsis thaliana (At),Oryza sativa (Os) and Setaria italica (Si). Conservative domains include Bromo_plant1/Bromodomain superfamily/Bromodomain/Bromo_gcn5_like, BET/BET superfamily, ZipA superfamily, PHA03247 superfamily, COG5076 superfamily, Hat1_N, NAT_SF and ELP3 superfamily.

Additional file 4: Fig. S2.

Phylogenetic trees and domain composition of CBP subfamily. Phylogenetic tree and domain composition of CBP subfamily predicted proteins from Arabidopsis thaliana (At),Oryza sativa (Os) and Setaria italica (Si). Conservative domains include HAT_KAT11 superfamily, PHD_HAC_like/PHD_SF superfamily, zf-TAZ/ZnF_TAZ/ZnF_UBR1, ZZ superfamily/ZZ_CBP/ZZ_dah/ZZ, Med15 superfamily, BTB_POZ superfamily, BACK and E3_UbLigase_R4.

Additional file 5: Fig. S3.

Phylogenetic trees and domain composition of MYST subfamily. Phylogenetic tree and domain composition of MYST subfamily predicted proteins from Arabidopsis thaliana (At), Oryza sativa (Os) and Setaria italica (Si). All members of MYST subfamily have a conserved domain PLN00104.

Additional file 6: Fig. S4. 

Phylogenetic trees and domain composition of TAF subfamily. Phylogenetic tree and domain composition of TAF subfamily predicted proteins from Arabidopsis thaliana (At), Oryza sativa (Os) and Setaria italica (Si). Conservative domains include DUF3591/DUF3591 superfamily, Bromodomain superfamily/Bromodomain/Bromo_AAA, Ubiquitin_like_fold superfamily, TBP-binding superfamily, zf-CCHC_6 superfamily, P-loop_NTPase superfamily and SpoVK.

Additional file 7: Fig. S5. 

Phylogenetic trees and three dimensional structures of GNAT proteins in Setaria italica. Bootstrap values higher than 50% are shown. In the same subfamily, the three dimensional structure is similar. The higher the bootstrap values, the closer the kinship and the more similar the three dimensional structure.

Additional file 8: Table S3.

One-to-one orthologous relationships between Setaria italica (Si) and Oryza sativa(Os).

Additional file 9: Table S4.

Prediction of cis-acting elements in promoters ofhistone acetylation genes (HATs) infoxtail millet.

Additional file 10: Table S5.

Details of 23tissues sampled in the spatial and temporal expression experiment.

Additional file 11.  

HATs Protein sequence of  Arabidopsis thaliana, Oryza sativa, and Setaria italica.

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Xing, G., Jin, M., Qu, R. et al. 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 22, 292 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • SiHAT genes
  • Expression analysis
  • Nitrogen deficiency
  • Phosphorus deficiency
  • Abiotic stress
  • Sclerospora graminicola infection