- Research article
- Open Access
Genome-wide characterization of NtHD-ZIP IV: different roles in abiotic stress response and glandular Trichome induction
BMC Plant Biology volume 19, Article number: 444 (2019)
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 . In tomato, exogenous application of JA resulted in a dramatic increase in glandular trichome density .
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 . Two closely-related and functionally-redundant AtHD-ZIP IVs, HDG11 and HDG12, regulate branching of the trichome . 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 . 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.
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).
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.
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.
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.
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.
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 . 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 . In plants, segmental duplication is the most frequent mechanism due to the property of diploidized polyploid . 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 . 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 . 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 .
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 (https://solgenomics.net/) and the Rice Genome Database (http://rice.plantbiology.msu.edu/), respectively. The A. thaliana HD-ZIP IV proteins were obtained using the Arabidopsis Information Resource (http://www.arabidopsis.org/). The Arabidopsis HD-ZIP IV proteins were used as query seeds to identify the N. tabacum HD-ZIP IV proteins (https://solgenomics.net/), 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 (http://circos.ca/).
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 .
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.
- HD-ZIP IV:
homeodomain-leucine zipper class IV family
Open reading frame
Quantitative real-time PCR
Reverse transcription-polymerase chain reaction
Lorenzo O, Solano R. Molecular players regulating the jasmonate signaling network. Curr Opin Plant Biol. 2005;8:532–40.
Mauch-Mani B, Mauch F. The role of abscisic acid in plant-pathogen interactions. Curr Opin Plant Biol. 2005;8:409–14.
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.
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.
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.
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.
Schilmiller AL, Last RL, Pichersky E. Harnessing plant trichome biochemistry for the production of useful compounds. Plant J. 2008;54:702–11.
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.
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.
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.
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.
Roodbarkelari F, Groot EP. Regulatory function of homeodomain-leucine zipper (HD-ZIP) family proteins during embryogenesis. New Phytol. 2016;213(1):95–104.
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.
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.
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.
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.
Ariel FD, Manavella PA, Dezar CA, Chan RL. The true story of the HD-zip family. Trends Plant Sci. 2007;12(9):419–26.
Javelle M, Vernoud V, Rogowsky PM, Ingram GC. Epidermis: the formation and functions of a fundamental plant tissue. New Phytol. 2011;189:17–39.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Wagner GJ. Secreting glandular trichomes: more than just hairs. Plant Physiol. 1991;96(3):675–9.
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.
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.
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.
El Sayed KA, Sylvester PW. Biocatalytic and semisynthetic studies of the anticancer tobacco cembranoids. Expert Opin Investig Drugs. 2007;16:877–87.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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)
About this article
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). https://doi.org/10.1186/s12870-019-2023-4
- Nicotiana tabacum
- Expression pattern