Skip to main content

Genome-wide characterization of NtHD-ZIP IV: different roles in abiotic stress response and glandular Trichome induction



The plant-specific homeodomain-leucine zipper class IV (HD-ZIP IV) gene family has been involved in the regulation of epidermal development.


Fifteen genes coding for HD-ZIP IV proteins were identified (NtHD-ZIP-IV-1 to NtHD-ZIP-IV-15) based on the genome of N. tabacum. Four major domains (HD, ZIP, SAD and START) were present in these proteins. Tissue expression pattern analysis indicated that NtHD-ZIP-IV-1, − 2, − 3, − 10, and − 12 may be associated with trichome development; NtHD-ZIP-IV-8 was expressed only in cotyledons; NtHD-ZIP-IV-9 only in the leaf and stem epidermis; NtHD-ZIP-IV-11 only in leaves; and NtHD-ZIP-IV-15 only in the root and stem epidermis. We found that jasmonates may induce the generation of glandular trichomes, and that NtHD-ZIP-IV-1, − 2, − 5, and − 7 were response to MeJA treatment. Dynamic expression under abiotic stress and after application of phytohormones indicated that most NtHD-ZIP IV genes were induced by heat, cold, salt and drought. Furthermore, most of these genes were induced by gibberellic acid, 6-benzylaminopurine, and salicylic acid, but were inhibited by abscisic acid. NtHD-ZIP IV genes were sensitive to heat, but insensitive to osmotic stress.


NtHD-ZIP IV genes are implicated in a complex regulatory gene network controlling epidermal development and abiotic stress responses. The present study provides evidence to elucidate the gene functions of NtHD-ZIP IVs during epidermal development and stress response.


Plants have developed a complex regulatory network to adapt to extreme environmental stresses, in which jasmonic acid (JA), salicylic acid (SA) and abscisic acid (ABA) act as pivotal defense signal molecules [1,2,3]. Plant trichomes are involved in defense responses towards insect predation, UV damage, toxin sequestration, and excess transpiration. Trichomes are grouped into two types, glandular and non-glandular. Glandular trichomes can synthesize and secrete large numbers of specialized metabolites, including terpenes, phenylpropanoids, sucrose esters, and flavonoids [4, 5]. These natural plant compounds not only protect plants against insect pests, but also contribute to the production of industrial chemicals for use in flavors, aromas, and pharmaceuticals [6,7,8]. In Arabidopsis, it was reported that exogenous application of JA and GA induced the occurrence of non-glandular trichomes [9]. In tomato, exogenous application of JA resulted in a dramatic increase in glandular trichome density [10].

The plant-specific homeodomain-leucine zipper (HD-Zip) gene family plays a crucial role in abiotic stress response and plant development [11,12,13]. These proteins can be further grouped into 4 subfamilies according to their structural features, conserved domains, and physiological functions [[14,15,16,17]. The Class IV HD-Zip (hereafter “HD-Zip IV”) gene family is associated with lipid transport, epidermal development, cuticle biosynthesis, and anthocyanin deposition [18,19,20]. HD-ZIP IVs are also implicated in mediating plant defense to osmotic stress [21, 22]. In Arabidopsis, the HD-ZIP IV family comprises 16 genes; the first identified HD-ZIP IV gene (GL2) was implicated in root hair differentiation and trichome development [23, 24]. Two AtHD-ZIP IVs, ML1 and PDF2, have been involved in regulating epidermis and embryo development and determining floral organ identity [25, 26]. One AtHD-ZIP IV gene, AtANL2, controls epidermal cell proliferation, root development, and anthocyanin accumulation [27]. Two closely-related and functionally-redundant AtHD-ZIP IVs, HDG11 and HDG12, regulate branching of the trichome [19]. HD-ZIP IVs has been characterized in various groups other than Arabidopsis, namely maize, rice, soybean, and cucumber [18, 19, 28,29,30]. It was found that HD-ZIP IVs are primarily expressed in the epidermal tissue. Moreover, Arabidopsis, maize, rice, soybean, and cucumber possess only non-glandular trichomes. The recently published expression profile of HD ZIP IVs in tomatoes suggests that each member may fulfill distinct functions in plant development [31]. Up to now, the specific roles of HD-ZIP IVs in the induction of glandular trichomes has remained enigmatic.

The common tobacco (Nicotiana tabacum), a broadleaf crop with large yields and planting areas, has glandular trichomes on the surface of its leaves. These trichomes produce various terpenoids, alkaloids and defensive proteins, together representing up to 30% dry weight of the leaf [32,33,34]. Diterpenoids, including labdanoids and cembranoids, are more abundant in Pinus and Nicotiana than in other genera [35, 36]. In addition, cembranoids have neuroprotective, anti-microbial, and anti-tumor properties, and can help in the treatment of human immunodeficiency virus [37,38,39,40]. However, knowledge concerning the occurrence of glandular trichomes is fragmentary.

N. tabacum is an excellent model to clarify the gene functions of HD-ZIP IVs in dicotyledons. To elucidate the potential functions of NtHD-ZIP IVs in abiotic stress response and plant development, N. tabacum HD-ZIP IV genes were identified by the computational analysis of N. tabacum genome resource. We analyze gene structure, synteny, phylogeny, tissue expression pattern, and the expression profile under various exogenous hormones and abiotic stresses. In particular, we compare the transcript level of HD-ZIP IVs in the sub-epidermal and epidermal layers. Our study lays the foundation for characterization of HD-ZIP IVs in epidermis-related functions.


Identification and analysis of HD-ZIP IV genes in N. tabacum

Based on the latest genome data of tobacco, 32 HD-ZIP IV genes were identified in N. tabacum genomes. These HD-ZIP IV proteins had conserved domains namely HD, LZ, SAD and START. The positions of the HD-ZIP IV genes showed a scattered distribution pattern in the tobacco chromosome (Table 1). Chromosome 4 had three HD-ZIP IV gene copies, chromosomes 1, 11, 13, and 23 contained two copies, and chromosomes 2, 6, 8, 10, 12, 14, 17, and 22 individually had one copy. Moreover, nine pairs of HD-ZIP IVs were duplicated in tobacco genome (Fig. 1). The molecular weight of HD-ZIP IV proteins ranged from 49.66 to 91.77 kDa, the predicted full-length amino acid sequences ranged from 359 to 828, and the number of exons ranged from 4 to 11.

Table 1 HD ZIP IV gene family in N. tabacum
Fig. 1

Chromosome distributions and synteny relationships of HD-ZIP IVs in N. tabacum. Gray lines show segmental duplications

Evolutionary analysis showed that 74 HD-ZIP IV proteins (32 from tobacco, 13 from tomato, 16 from Arabidopsis, and 13 from rice) were clustered into 5 groups (Fig. 2). Each group contained HD-ZIP IVs from the four species. The result of the phylogenetic analysis was consistent with the taxonomic classification: the HD-ZIP IV genes from the solanaceous plants (tobacco and tomato) had highly homologous sequences; and the HD-ZIP IVs of eudicots (tobacco, tomato and Arabidopsis) were more closely clustered than were those of the monocot (rice).

Fig. 2

Phylogenetic tree of HD-ZIP IV proteins in different plant species. The protein sequences of 32 N. tabacum, 13 tomato, 15 rice, and 16 Arabidopsis HD-ZIP IVs were used for the phylogenetic analysis. ▲, Nicotiana tabacum; , Arabidopsis thaliana; ■, Oryza sativa; , Solanum lycopersicum

Gene structure analysis can give insights into the origin and evolution of the HD-ZIP IV gene family in tobacco. A phylogenetic tree was constructed to verify the consistency of the exon-intron pattern and the phylogenetic classification. The tobacco HD-ZIP IV genes were divided into 15 categories, which we designate with the prefix “Nt”: NtHD-ZIP-IV-1 to NtHD-ZIP-IV-15 (Fig. 3a). The closely related NtHD-ZIP IV genes had a similar gene structure. Similarly to the situation found in other plants, the features of the NtHD-ZIP IV gene family varied substantially, and the exon number varied from 4 to 11 (Fig. 3b). It is noteworthy that non intron sequence was inserted in the conserved domain. Analysis of conserved motifs found that 20 motifs were present in the 15 NtHD-ZIP IV proteins (Fig. 3c, Additional file 1: Figure S1). There were usually similar motif patterns in closely-related proteins in the phylogenetic tree, thus indicating evolutionary and functional conservation within a clade.

Fig. 3

Multiple sequence alignment, gene structure, and conserved motif analysis of NtHD-ZIP IVs. a Multiple sequence alignment of NtHD-ZIP Vs in N. tabacum. b Exon-intron structure analysis of NtHD-ZIP IV genes. Introns and exons are indicated by black lines and rectangles, respectively. c Analysis of the conserved motifs. Conserved motifs are labeled with different colored frames

Spatial gene expression of HD-ZIP IVs

The expression pattern of 15 NtHD-ZIP IVs in five tobacco tissues was investigated to assess their role in the epidermal development. There were no trichomes on the cotyledons, but many glandular and non-glandular trichomes occurred on the outer surface of leaves and stems (Fig. 4a). As shown in Fig. 4b, NtHD-ZIP-IV-1 and NtHD-ZIP-IV-2 were specifically expressed in the leaf, root, and stem epidermis. No expression of these genes was detected in cotyledons and in stems without epidermis. This indicates that NtHD-ZIP-IV-1 and NtHD-ZIP-IV-2 are trichome-specific genes. The expression of NtHD-ZIP-IV-3, NtHD-ZIP-IV-10, and NtHD-ZIP-IV-12 was weak in stems without epidermis; this suggests that these three genes may relate to trichome development. Five NtHD-ZIP IV genes (NtHD-ZIP-IV-4, − 5, − 6, − 13, and − 14), had similar expression patterns: not expressed in cotyledons, but expressed in leaves, roots, stem epidermis, and stems without epidermis. NtHD-ZIP-IV-7 showed a consistent transcript level in five tissues. These genes may have complex roles in tobacco development. Notably, NtHD-ZIP-IV-8 was expressed only in the cotyledons, NtHD-ZIP-IV-9 only in the leaf and stem epidermis, NtHD-ZIP-IV-11 only in the leaf, and NtHD-ZIP-IV-15 only in the root and stem epidermis. These results indicated that each NtHD-ZIP IV gene may be associated with the development of different plant organs.

Fig. 4

Spatial expressional analysis of NtHD-ZIP IVs. a Morphological features of the epidermis on various tissues. Scale bar = 100 μm. b Gene transcript levels in various tobacco tissues. The lowest transcription for each gene was regarded as a standard, and L25 gene was taken as endogenous control. Gels: upper, NtHD-ZIP IV gene segments amplified by semi-quantitive RT-PCR; lower, L25 gene segments amplified by semi-quantitive RT-PCR. Data was analyzed using one-way ANOVA followed by least significant difference (LSD) to determine the significance of differences between means using SPSS version 11.0. Each bar represents the average of three biological replicates. Different letters in the same gene indicate significant differences (P < 0.05)

MeJA application induced the initiation of long-stalk glandular trichomes

Only non-glandular and short-stalked glandular trichomes are present on the surface of tobacco T.I.1112 plants. After MeJA application, long-stalk glandular trichomes were observed, and the density was significantly increased with the increase of the MeJA concentration; this was not the case for the non-glandular and short-stalked glandular trichomes (Fig. 5a, b). These results indicated that the morphogenesis of different trichome types was regulated by different networks.

Fig. 5

Effects of MeJA on long-stalk glandular trichomes. a Exogenous MeJA application induced the initiation of long-stalk glandular trichomes. Scale bar = 100 μm. b Trichome density affected by MeJA application. Different letters show significant differences (P ≤ 0.05). c Transcript levels of NtHD-ZIP IVs upon 5.0 mM MeJA application. The lowest transcript level for each gene was regarded as a standard. L25 was selected as a control gene. The results were calculated by the 2-ΔΔCT method

Detecting the transcript level of 15 NtHD-ZIP-IVs in the epidermis found that most NtHD-ZIP IV genes were not response to MeJA treatment except NtHD-ZIP-IV-1, − 2, − 5, and − 7. NtHD-ZIP-IV-1, − 2, and − 7 were inhibited under MeJA application, whereas the transcription level of NtHD-ZIP-IV-5 increased after MeJA treatment (Fig. 5c).

Expression pattern of NtHD-ZIP IV genes under abiotic stress and hormone treatments

Plant hormones are key regulators in plant growth and development, as are various environmental stimuli. The NtHD-ZIP IV genes had diverse responses to the various hormone treatments (Fig. 6). Following ABA treatment, NtHD-ZIP-IV-1, 2, − 3, − 5, − 7, − 9, − 10, and − 13 were inhibited, NtHD-ZIP-IV-6, − 11, and − 12 were slightly induced, whereas NtHD-ZIP-IV-4, − 8, − 13, and − 14 showed no response. GA treatment induced expression of NtHD-ZIP-IV-4, − 5, − 6, − 9, − 10, − 12, and − 13, whereas the remaining NtHD-ZIP IV genes showed no response. Similarly, most NtHD-ZIP IV genes were activated by 6-BA treatment (but not NtHD-ZIP-IV-8, − 10, and − 15). Following SA treatment, NtHD-ZIP-IV-1, − 2, and − 3 were inhibited, NtHD-ZIP-IV-6, − 8, and − 15 did not respond, and the transcript level of the remaining NtHD-ZIP IV genes increased. Compared with other NtHD-ZIP IV genes, NtHD-ZIP-IV-9 and -14 could be upregulated at a constant rate by exogenous SA. The findings suggested that NtHD-ZIP IVs might be implicated in complex networks, with each member having distinct funtions.

Fig. 6

Heat maps of NtHD-ZIP IVs gene expression under different hormones and abiotic stresses. qRT-PCR was performed to analyze the transcript level of NtHD-ZIP IVs, and the results were calculated by the 2-ΔΔCT method. L25 was selected as a control gene. The lowest transcript level at each treatment for each gene was set as 1

We found that that most NtHD-ZIP IVs could be activated by abiotic stress, to varying degrees. NtHD-ZIP IV genes were more sensitive to heat stress than to salt, drought, and cold stress. Following high salinity treatment, the expression of NtHD-ZIP-IV-1, − 2, − 3, − 9, − 11, − 12, and − 13 was upregulated, whereas the remaining genes showed no clear changes. Under drought stress, most of NtHD-ZIP IV genes were up-regulated, except for NtHD-ZIP-IV-8 and -15. Conversely, most NtHD-ZIP IV genes were not obviously activated by cold, except for NtHD-ZIP-IV-6, − 7, − 10, and − 11. Among the four genes, NtHD-ZIP-IV-7 and -11 showed the strongest response to cold stress. Under heat stress, most NtHD-ZIP IV genes were significantly activated, except for NtHD-ZIP-IV-5, − 8, − 9, − 13, and − 15; in those that were activated, the expression levels were high. Moreover, the transcription levels of NtHD-ZIP-IV-1, − 2, − 4, − 6, − 7, and − 11 were high at each sampling occasion after application of heat stress.


HD-ZIP IV genes are conserved during evolution

Plant-specific HD-ZIP IVs have been involved in the regulation of epidermal development, including root hairs, trichomes, cuticles and stomates [17]. Sequence analysis indicates that HD-ZIP IVs are highly conserved during the evolution or separation processes of various plant species. However, there were different epidermal characteristics in different plant species. Here, we characterized the 15 NtHD-ZIP IV genes in the tobacco genome.

There are mainly three principal evolutionary mechanisms of gene duplications: tandem duplication, segmental duplication, and transposition events [41]. In plants, segmental duplication is the most frequent mechanism due to the property of diploidized polyploid [42]. In the present study, some NtHD-ZIP IV genes were distributed in duplicated blocks, indicating that segmental duplications have contributed to the gene duplication of NtHD-ZIP IVs. The phylogenetic analysis of the HD-ZIP IVs from different species suggested that the HD-ZIP IV duplications within species were first clustered into the same clade, and then grouped together with other species. This finding indicates that HD-ZIP IVs diversified and expanded after the radiation of species.

Some HD-ZIP IVs may play crucial roles in the trichome formation

The regulatory network that controls unicellular trichome formation has been well studied in Arabidopsis [43, 44]. Similarly to the situation found in the tomato and potato, trichomes in tobacco are typically multicellular structures. Prior to now, there has been a fragmented understanding of the molecular mechanisms underlying multicellular trichome formation. Two paralogous HD-ZIP IVs (HDG11 and HDG12) were involved in the trichome branching. Specifically, hdg11 mutants had more branched trichomes in the leaves, and hdg12 mutants had more normal trichomes than occurred in the wild type. The excessive-branching morphology of the trichome in hdg11 mutants was enhanced by hdg12, revealing a synergistic effect on the trichome development. In our study, the HDG11 and HDG12 homologous genes were HD-ZIP-IV-11 and HD-ZIP-IV-12. NtHD-ZIP-IV-11 was expressed only in the leaf. NtHD-ZIP-IV-12 was strongly expressed in the cotyledons and stem epidermis, whereas weak expression was detected in stems without epidermis. From this, we deduce that NtHD-ZIP-IV-12 may be related to epidermal development. In tomato, an HD-ZIP IV gene (Wo) is involved in the initiation of multicellular trichomes [45, 46]. Suppression of Wo expression by RNA interference decreases the density of type I trichomes. The homologous gene of Wo in Arabidopsis, PDF2, may regulate shoot epidermal cell differentiation [47]. These results indicate that the formation of multicellular trichomes might be regulated by a distinct network unlike the unicellular trichomes. Further, these HD-ZIP IVs may act different roles in the initiation of the unicellular and multicellular trichomes. Here, the predicted proteins coded by the Wo gene showed 73, 75, 78, and 79% amino acid sequence identity to the four Wo homologs in tobacco, which were further clustered as NtHD-ZIP-IV-1 and NtHD-ZIP-IV-2 (Additional file 1: Figure S2). Tissue-preferential expression pattern is an indication of the specific gene function. We found that NtHD-ZIP-IV-1 and NtHD-ZIP-IV-2 were trichome-specific genes. Moreover, NtHD-ZIP-IV-1 and NtHD-ZIP-IV-2 were strongly upregulated under MeJA application, followed by the initiation of secreting trichomes. Our results indicate that NtHD-ZIP-IV-1 and -2 may act crucial roles in the induction of the secreting trichome, similar to the role of Wo in tomato.

Diverse HD-ZIP IVs were implicated in hormone and abiotic stress response

In the present study, JA could induce the generation of glandular trichomes. Recently, HDG11 in Arabidopsis, a homologous gene of NtHD-ZIP-IV-11, has been reported to control the JA biosynthesis [48]. However, NtHD-ZIP-IV-11 was not responsive to JA in our present study. The transcripts of most NtHD-ZIP-IVs were not respond to MeJA treatment, except for NtHD-ZIP-IV-1, − 2, − 5, and − 7, which may play important roles in the induction of secreting trichomes. Surprisingly, most NtHD-ZIP IVs responded to ABA, GA, 6-BA, and SA. These hormones are key signaling regulators in plant responses to abiotic stresses [49].

This study primarily focused on determining the dynamic transcriptional changes in NtHD-ZIP IV genes under various abiotic stresses. The results indicate that most NtHD-ZIP IVs were sensitive to heat, but insensitive to cold and osmotic stress; each NtHD-ZIP IV gene had distinct functions; and NtHD-ZIP IVs were implicated in a complex network of responses to abiotic stress. The NtHD-ZIP IV genes might be good target genes for improving abiotic-stress tolerance in crop plants.


Fifteen HD-ZIP IV genes were identified from N. tabacum genome. These NtHD-ZIP-IVs showed differential tissue-specific expression patterns. Jasmonates could induce the generation of glandular trichome, and four NtHD-ZIP-IVs were implicated in glandular trichome induction. Each NtHD-ZIP IV gene had a distinct role in abiotic stress and phytohormone response. The present study provides evidence to elucidate the gene functions of NtHD-ZIP IVs in epidermal development and stress responses.


Analysis of the HD-ZIP IV gene family in N. tabacum

The sequence of the Solanum lycopersicum and Oryza sativa HD-ZIP IV gene family was obtained at the Solanaceae Genome Network ( and the Rice Genome Database (, respectively. The A. thaliana HD-ZIP IV proteins were obtained using the Arabidopsis Information Resource ( The Arabidopsis HD-ZIP IV proteins were used as query seeds to identify the N. tabacum HD-ZIP IV proteins (, via a BlastP search (e < 1− 10). These predicted HD-ZIP IV proteins were further confirmed and analyzed using the Pfam tool and SMART web server. The biophysical properties of the HD-ZIP IVs were estimated with the ExPASy ProtParam tool.

To estimate the phylogeny of the HD-ZIP IV genes, phylogenetic analysis was carried out using MEGA 7.0 with 1000 replicates, using the HD-ZIP IVs in tobacco, tomato, rice and Arabidopsis. Sequences were aligned with ClustalW program. Gene structure was visualized with the Gene Structure Display Server 2.0. The Multiple Expectation Maximization for Motif Elicitation tool was performed to identify the conserved motif. To determine synteny, the synteny blocks containing HD-ZIP IVs in the N. tabacum genome were scaned using the MCScanX project. The position of each gene in the corresponding chromosome and its synteny relationship were generated using Circos (

Tissue-specific expression analysis

N. tabacum ‘K326’ seedlings were raised in a growth chamber at 22 °C with a 12/12 h light-dark photoperiod. For the tissue-specific expression analysis, cotyledons were sampled from one-week-old seedlings, and the leaf, root, stem epidermis, and the stem with its epidermis removed were sampled from three-week-old seedlings.

Total RNA was extracted and removed the residual DNA using DNase I. Quantitative real-time PCR (qRT-PCR) and semi-quantitative RT-PCR were employed to determine the relative mRNA transcriptions of HD-ZIP IVs in five tobacco tissues using the gene-specific primers (Additional file 1: Table S1). L25 gene was selected as an internal control. q-PCR reaction was performed on an ABI PRISM 7000 system (Applied Biosystems, USA) with the SYBR Green RT-PCR Kit (Takara, China). Each reaction was run in triplicate, and analysis was performed using the 2-ΔΔCT method [50].

Induction of long-stalk glandular trichomes by MeJA

N. tabacum T.I.1112 without long-stalked glandular trichomes was developed by the Oxford Tobacco Research Station. Seedlings at the four-leaf stage were sprayed with 5.0 mM methyl jasmonate (MeJA). Plants were sprayed until all plants were saturated. Three applications were repeated every one week. Three weeks later, three plants from each treatment were selected, and the youngest terminal leaflet at least 5 cm in length on each plant was sampled for the trichome morphology observation. The area of glandular head, and trichome density on the upper leaf surface were counted using an Axioplan 2 microscope (Zeiss, Oberkochen, Germany). The morphological data were analyzed using one-way ANOVA. Moreover, the leaf epidermis of plants exposed to the 5.0 mM MeJA treatment, and of the control, was removed to analyze the expression level of HD-ZIP IVs.

Abiotic stress and hormone treatments

To test the effects of abiotic stress, K326 tobacco seedlings at the four-leaf stage were stressed by placing the plants under one of four treatments: application of 300 mM NaCl or PEG-6000 (− 0.5 MPa) solutions; and exposure to low (4 °C) or high (42 °C) temperatures. In preliminary studies, we found that these treatments caused significant stress to the plants. Control plants were cultured normally without treatment.

To test the effects of exogenous hormone treatment, seedlings at the four-leaf stage were sprayed separately with 100 μM abscisic acid (ABA), 100 μM 6-benzylaminopurine (6-BA), 2.0 mM salicylic acid (SA), and 150 μM gibberellic acid (GA). Control seedlings were sprayed with distilled water. True leaves were collected at 0, 1, 3, 6, 12, 24, 48, and 72 h post treatment for q-PCR analysis.

Availability of data and materials

All data generated in this study is available as Additional files.





Abscisic acid




homeodomain-leucine zipper class IV family


Methyl jasmonate


Open reading frame


Quantitative real-time PCR


Reverse transcription-polymerase chain reaction


Salicylic acid


  1. 1.

    Lorenzo O, Solano R. Molecular players regulating the jasmonate signaling network. Curr Opin Plant Biol. 2005;8:532–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Mauch-Mani B, Mauch F. The role of abscisic acid in plant-pathogen interactions. Curr Opin Plant Biol. 2005;8:409–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Qu Y, Wang YY, Yin QS, Huang LL, Jiang YG, Li GZ, Hao L. Multiple biological processes involved in the regulation of salicylic acid in Arabidopsis response to NO2 exposure. Environ Exp Bot. 2018;149:9–16.

    Article  CAS  Google Scholar 

  4. 4.

    Gang DR, Wang J, Dudareva N, Nam KH, Simon JE, Lewinsohn E, Pichersky E. An investigation of the storage and biosynthesis of phenylpropenes in sweet basil. Plant Physiol. 2001;125:539–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Kroumova AB, Wagner GJ. Different elongation pathways in the biosynthesis of acyl groups of trichome exudates sugar esters from various solanaceous plants. Planta. 2003;216:1013–21.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Zhang H, Zhang S, Yang Y, Jia H, Cui H. Metabolic flux engineering of cembratrien-ol production in both the glandular trichome and leaf mesophyll in Nicotiana tabacum. Plant Cell Physiol. 2018;59(3):566–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Schilmiller AL, Last RL, Pichersky E. Harnessing plant trichome biochemistry for the production of useful compounds. Plant J. 2008;54:702–11.

    Article  CAS  Google Scholar 

  8. 8.

    Rios-Estepa R, Turner GW, Lee JM, Croteau RB, Lange BM. A systems biology approach identifies the biochemical mechanisms regulating monoterpenoid essential oil composition in peppermint. Proc Natl Acad Sci U S A. 2008;105:2818–23.

    Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Traw MB, Bergelson J. Interactive effects of jasmonic acid, salicylic acid, and gibberellin on induction of trichomes in Arabidopsis. Plant Physiol. 2003;133:1367–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Boughton AJ, Hoover K, Felton GW. Methyl jasmonate application induces increased densities of glandular trichomes on tomato, Lycopersicon esculentum. J Chem Ecol. 2005;31(9):2211–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Perotti MF, Ribone PA, Chan RL. Plant transcription factors from the homeodomain-leucine zipper family I. role in development and stress responses. IUBMB Life. 2017;69(5):280–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Roodbarkelari F, Groot EP. Regulatory function of homeodomain-leucine zipper (HD-ZIP) family proteins during embryogenesis. New Phytol. 2016;213(1):95–104.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Yan T, Li L, Xie L, Chen M, Shen Q, Pan Q, Fu X, Shi P, Tang Y, Huang H. A novel HD-ZIP IV/MIXTA complex promotes glandular trichome initiation and cuticle development in Artemisia annua. New Phytol. 2018;218(2):567–78.

    Article  CAS  Google Scholar 

  14. 14.

    Ruberti I, Sessa G, Lucchetti S, Morelli G. A novel class of plant proteins containing a homeodomain with a closely linked leucine zipper motif. EMBO J. 1991;10:1787–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Johannesson H, Wang Y, Engström P. DNA-binding and dimerization preferences of Arabidopsis homeodomain-leucine zipper transcription factors in vitro. Plant Mol Biol. 2001;45:63–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Sessa G, Carabelli M, Ruberti I, Lucchetti S, Baima S, Morelli G. Identification of distinct families of HD-ZIP proteins in Arabidopsis thaliana. Plant Mol Biol. 1994;81:411–26.

    Article  CAS  Google Scholar 

  17. 17.

    Ariel FD, Manavella PA, Dezar CA, Chan RL. The true story of the HD-zip family. Trends Plant Sci. 2007;12(9):419–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Javelle M, Vernoud V, Rogowsky PM, Ingram GC. Epidermis: the formation and functions of a fundamental plant tissue. New Phytol. 2011;189:17–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Nakamura M, Katsumata H, Abe M, Yabe N, Komeda Y, Yamamoto KT, Takahashi T. Characterization of the class IV homeodomain-leucine zipper gene family in Arabidopsis. Plant Physiol. 2006;141(4):1363–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Zalewski CS, Floyd SK, Furumizu C, Sakakibara K, Stevenson DW, Bowman JL. Evolution of the class IV HD-zip gene family in Streptophytes. Mol Biol Evol. 2013;30(10):2347–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Yu H, Chen X, Hong YY, Wang Y, Xu P, Ke SD, Liu HY, Zhu JK, Oliver DJ, Xiang CB. Activated expression of an Arabidopsis HD-START protein confers drought tolerance with improved root system and reduced stomatal density. Plant Cell. 2008;20(4):1134–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Chen E, Zhang X, Yang Z, Wang X, Yang Z, Zhang C, Wu Z, Kong D, Liu Z, Zhao G. Genome-wide analysis of the HD-ZIP IV transcription factor family in Gossypium arboreum and GaHDG11 involved in osmotic tolerance in transgenic Arabidopsis. Mol Gen Genomics. 2017;292(3):1–17.

    Article  CAS  Google Scholar 

  23. 23.

    Wang S, Kwak SH, Zeng Q, Ellis BE, Chen XY, Schiefelbein J, Chen JG. TRICHOMELESS1 regulates trichome patterning by suppressing GLABRA1 in Arabidopsis. Development. 2007;134(21):3873.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Masucci JD, Rerie WG, Foreman DR, Zhang M, Galway ME, Marks MD, Schiefelbein JW. The homeobox gene GLABRA2 is required for position-dependent cell differentiation in the root epidermis of Arabidopsis thaliana. Development. 1996;122:1253–60.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Kamata N, Okada H, Komeda Y, Takahashi T. Mutations in epidermis-specific HD-ZIP IV genes affect floral organ identity in Arabidopsis thaliana. Plant J. 2013;75:430–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Ogawa E, Yamada Y, Sezaki N, Kosaka S, Kondo H, Kamata N, Abe M, Komeda Y, Takahashi T. ATML1 and PDF2 play a redundant and essential role in Arabidopsis embryo development. Plant Cell Physiol. 2015;56(6):1183.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Kubo H, Peeters AJ, Aarts MG, Pereira A, Koornneef M. ANTHOCYANINLESS2, a homeobox gene affecting anthocyanin distribution and root development in Arabidopsis. Plant Cell. 1999;11:1217–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Fu R, Liu W, Li Q, Li J, Wang L, Ren Z. Comprehensive analysis of the homeodomain-leucine zipper IV transcription factor family in Cucumis sativus. Genome. 2013;56:395–405.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Belamkar V, Weeks NT, Bharti AK, Farmer AD, Graham MA, Cannon SB. Comprehensive characterization and RNA-Seq profiling of the HD-zip transcription factor family in soybean (Glycine max) during dehydration and salt stress. BMC Genomics. 2014;15:950.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Javelle M, KleinCosson C, Vernoud V, Boltz V, Maher C, Timmermans M, DepègeFargeix N, Rogowsky PM. Genome-wide characterization of the HD-ZIP IV transcription factor family in maize: preferential expression in the epidermis. Plant Physiol. 2011;157(2):790–803.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Gao Y, Gao S, Xiong C, Yu G, Chang J, Ye Z, Yang C. Comprehensive analysis and expression profile of the homeodomain leucine zipper IV transcription factor family in tomato. Plant Physiol Bioch. 2015;96:141–53.

    Article  CAS  Google Scholar 

  32. 32.

    Wagner GJ. Secreting glandular trichomes: more than just hairs. Plant Physiol. 1991;96(3):675–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Kennedy BS, Nielsen MT, Severson RF. Biorationals from nicotiana protect cucumbers against Colletotrichum-Lagenarium (pass) ell & halst disease development. J Chem Ecol. 1995;21:221–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Lin Y, Wagner GJ. Surface disposition and stability of pest-interactive, trichome-exuded diterpenes and sucrose esters of tobacco. J Chem Ecol. 1994;20(8):1907–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Liu X, Zhang J, Liu Q, Tang G, Wang H, Fan C, Yin S. Bioactive cembranoids from the South China Sea soft coral Sarcophyton elegans. Molecules. 2015:13324–35.

  36. 36.

    El Sayed KA, Sylvester PW. Biocatalytic and semisynthetic studies of the anticancer tobacco cembranoids. Expert Opin Investig Drugs. 2007;16:877–87.

    Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Duan S, Du Y, Hou X, Li D, Ren X, Dong W, Zhao W, Zhang Z. Inhibitory effects of tobacco extracts on eleven phytopathogenic fungi. Nat Prod Res Dev. 2015;27:470–80.

    Google Scholar 

  38. 38.

    Zubair MS, Anam S, Al-Footy KO, Abdel-Lateef A, Alarif WM. Cembranoid diterpenes as antitumour: molecular docking study to several protein receptor targets. Proc Int Conf Comput Sci Technol. 2014;2:121–5.

    Google Scholar 

  39. 39.

    Martins AH, Hu J, Xu Z, Mu C, Alvarez P, Ford BD, El Sayed K, Eterovic VA, Ferchmin PA, Hao J. Neuroprotective activity of (1S,2E,4R,6R,-7E,11E)-2,7,11-cembratriene-4,6-diol (4R) in vitro and in vivo in rodent models of brain ischemia. Neuroscience. 2015;291:250–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Ning Y, Du Y, Liu X, Zhang H, Liu Y, Shi J, Xue SJ, Zhang Z. Analyses of effects of α-cembratrien-diol on cell morphology and transcriptome of Valsa Mali var. Mali. Food Chem. 2017;214:110–8.

    Article  CAS  Google Scholar 

  41. 41.

    Gu Z, Steinmetz LM, Gu X, Scharfe C, Davis RW, Li W. Role of duplicate genes in genetic robustness against null mutations. Nature. 2003;421:63–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Cannon SB, Mitra A, Baumgarten A, Young ND, 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.

    Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Szymanski DB, Lloyd AM, Marks MD. Rogress in the molecular genetic analysis of trichome initiation and morphogenesis in Arabidopsis. Trends Plant Sci. 2000;5(5):214–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Larkin JC, Brown ML, Schiefelbein J. How do cells know what they want to be when they grow up? Lessons from epidermal patterning in Arabidopsis. Annu Rev Plant Biol. 2003;54:403–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Yang C, Li H, Zhang J, Luo Z, Gong P, Zhang C, Li J, Wang T, Zhang Y, Lu YE. A regulatory gene induces trichome formation and embryo lethality in tomato. Proc Natl Acad Sci U S A. 2011;108:11836–41.

    Article  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Gao S, Gao Y, Xiong C, Yu G, Chang J, Yang Q, Yang C, Ye Z. The tomato B-type cyclin gene, SlCycB2, plays key roles in reproductive organ development, trichome initiation, terpenoids biosynthesis and Prodenia litura defense. Plant Sci. 2017:103–14.

  47. 47.

    Abe M, Katsumata H, Komeda Y, Takahashi T. Regulation of shoot epidermal cell differentiation by a pair of homeodomain proteins in Arabidopsis. Development. 2003;130:635–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Cai XT, Xu P, Wang Y, Xiang CB. Activated expression of AtEDT1/HDG11 promotes lateral root formation in Arabidopsis mutant edt1 by upregulating jasmonate biosynthesis. J Integr Plant Biol. 2015;57:1017–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. 49.

    Gómez-Cadenas A, Ollas CD, Manzi M, Arbona V. Phytohormonal crosstalk under abiotic stress. In: Tran LS, Pal S, editors. Phytohormones: a window to metabolism, signaling and biotechnological applications. New York: Academic; 2014. p. 289–321.

    Google Scholar 

  50. 50.

    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8.

    Article  CAS  Google Scholar 

Download references


We thank the Oxford Tobacco Research Station in Oxford, North Carolina, USA for the T.I.1112 and K326 seeds.


The research was financially supported by the National Science Foundation of Henan province [Grant No. 182300410094], Technology Center, China Tobacco Henan Industrial Co., Ltd. [Grant No. ZW2014004], State Tobacco Monopoly Administration of China [Grant No. 110201401003 (JY-03)], and the Key Research Project of Institution of Higher Education of Henan Province [Grant No. 19A210003]. These funding bodies had no role in the design of the study; in the collection, analyses, and interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Author information




HYZ and XDM conceived and write the manuscript. WJL and DXN performed the expression pattern experiment. ZJW and XXY performed the genome-wide characterization. XLY and YFY participated to the data analysis. HC projected design and supervision. All authors carefully checked and approved this version of the manuscript.

Corresponding author

Correspondence to Hong Cui.

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.

Specific primers of HD-ZIP IV in qRT-PCR. Figure S1. Motif analysis of the NtHD-ZIP IV proteins. The 20 motifs were analyzed using the MEME online tool. Different letters represent the abbreviation of various amino acids. The higher the letter height, the stronger the conservatism of the amino acid at that position. Figure S2. Sequence alignment of NtHD-ZIP IV proteins and Wo from S. lycopersicum. Alignments were performed using Megalign program of DNAStar. Identical amino acid residues are shared in black background. Dashed lines represent gaps that were introduced to maximize alignment. (DOCX 568 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhang, H., Ma, X., Li, W. et al. Genome-wide characterization of NtHD-ZIP IV: different roles in abiotic stress response and glandular Trichome induction. BMC Plant Biol 19, 444 (2019).

Download citation


  • Nicotiana tabacum
  • Expression pattern
  • MeJA
  • ABA
  • GA
  • SA