Next Article in Journal
Phospholipid Membrane Protection by Sugar Molecules during Dehydration—Insights into Molecular Mechanisms Using Scattering Techniques
Next Article in Special Issue
Endoplasmic Reticulum (ER) Stress Response and Its Physiological Roles in Plants
Previous Article in Journal
Homology Models of Melatonin Receptors: Challenges and Recent Advances
Previous Article in Special Issue
Proteins Involved in Distinct Phases of Cold Hardening Process in Frost Resistant Winter Barley (Hordeum vulgare L.) cv Luxor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Role of Homeodomain Leucine Zipper (HD-Zip) IV Transcription Factors in Plant Development and Plant Protection from Deleterious Environmental Factors

Australian Centre for Plant Functional Genomics, University of Adelaide, Waite Campus, Glen Osmond, South Australia 5064, Australia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2013, 14(4), 8122-8147; https://doi.org/10.3390/ijms14048122
Submission received: 31 January 2013 / Revised: 26 March 2013 / Accepted: 3 April 2013 / Published: 12 April 2013
(This article belongs to the Special Issue Abiotic and Biotic Stress Tolerance Mechanisms in Plants)

Abstract

:
Homeobox genes comprise an important group of genes that are responsible for regulation of developmental processes. These genes determine cell differentiation and cell fate in all eukaryotic organisms, starting from the early stages of embryo development. Homeodomain leucine zipper (HD-Zip) transcription factors are unique to the plant kingdom. Members of the HD-Zip IV subfamily have a complex domain topology and can bind several cis-elements with overlapping sequences. Many of the reported HD-Zip IV genes were shown to be specifically or preferentially expressed in plant epidermal or sub-epidermal cells. HD-Zip IV TFs were found to be associated with differentiation and maintenance of outer cell layers, and regulation of lipid biosynthesis and transport. Insights about the role of these proteins in plant cuticle formation, and hence their possible involvement in plant protection from pathogens and abiotic stresses has just started to emerge. These roles make HD-Zip IV proteins an attractive tool for genetic engineering of crop plants. To this end, there is a need for in-depth studies to further clarify the function of each HD-Zip IV subfamily member in commercially important plant species.

1. Background

Transcription factors (TFs) are key components underlying mechanisms that control gene expression in all living organisms. They define phenotypic diversity and evolutionary adaptation of organisms [1,2]. TFs typically consist of at least two major domains. The first is a DNA-binding domain (BD) that recognises target DNA sequences and the second is a transcriptional activation domain (AD) that initiates transcription by interacting with general transcription factors [3].
Several TF families consist of proteins containing a DNA-binding domain known as the homeodomain (HD). The first HD-containing TF was identified in Drosophila melanogaster[4]. The homeotic gene encoding this TF was found to be responsible for the development of antennae at the position of the second leg pair of Drosophila[4]. Later, a large number of homeotic genes containing a 180 bp long homeobox (HB) were identified and their involvement in the development of Drosophila was demonstrated. Proteins with highly conserved HDs have been found in other animals, fungi and plants, indicating that these TFs play an important role in developmental pathways in species across different kingdoms [57]. Many HB-containing genes in mammals and Drosophila are homeotic genes that determine body part formation in specific positions. However, none of the plant HB-containing genes exert any such homeotic effect.
The first plant gene encoding a HD-containing protein, KNOTTED1, was isolated from maize (Zea mays L.) [8]. Since the identification of KNOTTED1, a large number of plant HB-containing TF genes have been isolated. Based on the presence of specific functional domains, plant HD-containing TFs were classified into six families: (1) homeodomain leucine zipper (HD-Zip); (2) plant homeodomain with a finger domain (PHD); (3) Bell domain (Bell); (4) zinc finger with homeodomain (ZF-HD); (5) Wuschel homeobox (WOX); and (6) Knotted homeobox (KNOX) [9]. A comprehensive classification of plant homeobox genes has been reported [10].
The HD-Zip family of TFs is unique to the plant kingdom, suggesting their involvement in developmental processes specific for plants [11,12]. The HD-Zip family consists of four subfamilies or classes. A common feature of all four subfamilies is the presence of a leucine zipper domain (Zip) adjacent to HD, which is important for homo- and hetero-dimerisation [9]. Members of the HD-Zip II subfamily can be distinguished from HD-Zip I subfamily proteins by the presence of a conserved “CPSCE” motif, which is situated downstream of the leucine zipper [6]. HD-Zip III and IV TFs contain a STeroidogenic Acute Regulatory protein-related lipid Transfer (START) domain and a conserved START-associated domain (HD-SAD), that are absent in HD-Zip I and II proteins [13,14]. It is of note that the START domains have never been associated to work in conjunction with DNA binding in mammalian systems, while the START domains have primarily been associated with the plant HD-Zip domains. Members of the HD-Zip III subfamily can be distinguished from HD-Zip IV subfamily members by an additional C-terminal MEKHLA domain [14]. Unlike other known HD proteins that are able to bind DNA sequence as monomers, HD-Zip TFs bind to DNA and consequently activate/repress downstream gene(s) only as dimers [15]. Several genes encoding members of HD-Zip proteins were shown to be involved in abiotic stress responses [1622]. Class IV of HD-Zip TFs is involved in a number of vitally important processes, including the transcriptional control of epidermal and sub-epidermal cell fate, anthocyanin accumulation, lipid transport, and cuticle biosynthesis [23]. Involvement of HD-Zip IV TFs in mediating plant responses to drought have also been reported [17].
In this review we summarise existing knowledge on structural characteristics, expression patterns, potential interacting partners and target genes of HD-Zip IV TFs, as well as describe possible roles of the HD-Zip IV genes in plant development, defence and response to abiotic stresses. Potential applications of these genes in plant biotechnology will be discussed.

2. Structure of HD-Zip IV TFs and Function of Identified Domains

HD-Zip IV TFs consist of four conserved domains. A highly conserved HD domain contains 60 or 61 amino acid residues and is responsible for binding to a specific DNA sequence by forming a structure composed of three α-helices. It was shown that HD-Zip IV proteins from Arabidopsis preferentially bind an 11 bp-long palindromic sequence 5′-GCATT(A/T)AATGC-3′, which partly overlaps with the sequence of the L1 box (5′-TAAATG(C/T)A-3′) [24]. The L1 box is responsible for specific gene expression in the epidermal L1 layer [25]. This cis-element was identified in the promoter regions of several genes, including HD-Zip IV genes [25,26]. The L1 box was found in the promoter region of the Arabidopsis PROTODERMAL FACTOR1 (PDF1) gene, which encodes a putative extracellular proline-rich protein with expression highly specific to the L1 cell layer of shoot apices and the protoderm of organ primordia. Expression of PDF1 is regulated by an Arabidopsis HD-Zip IV TF known as MERISTEM LAYER1 (AtML1), through binding of the L1 box in the promoter region of the PDF1 gene [25]. The substitution of two base pairs (AT for GG) in the L1 box located in the PDF1 promoter abolished binding of the AtML1 protein to the promoter in electrophoretic mobility shift assays (EMSA) [25]. However, HD-containing proteins often show a poor DNA-binding specificity in EMSA experiments and for strong DNA binding the HD-Zip IV proteins might require assistance from other TFs and/or modifying enzyme(s) [27]. Therefore, a final conclusion about in vivo activation of genes containing the L1 box in their promoters by AtML1 was not reached [25,27]. It is noteworthy that the plant L1 box sequence (5′-TAAATG(C/T)A-3′) resembles that of the target sequence (5′-TTAATGGCC-3′) of some homeotic proteins from Drosophila, such as Ultrabithorax [28].
A total of 53 TFs belonging to the HD-ZIP class IV are listed in Table 1. Phylogenetic relationships among individual HD-Zip IV proteins were established using either full length amino acid sequences (Figure 1A), or only HD regions of 43 TFs (Figure 1B). We have investigated phylogeny of TaGL9 from wheat (Triticum aestivum L.) in the context of 53 related TFs of the HD-Zip class IV from a range of plant species (Figure 1). Protein sequences were obtained using the BLAST tool from the NCBI database (Table 1). We constructed two types of phylogenetic trees of HD-Zip IV TFs based on either their full-length amino acid sequences (Figure 1A) or their HD domains (Figure 1B). The two types of trees using full-length and HD domain sequences were constructed to understand, which part of sequences (or domains) succumbed to selective pressures. The comparative analyses of these trees indicated that TFs from mono- and dicotyledonous species (mono- and di-cotyledonous species are indicated in bold and normal types, respectively) may have evolved independently, as distributions of TFs in individual sub-branches in both trees segregated. This observation suggested that a divergence between the monocotyledonous and dicotyledonous HD-Zip IV genes occurred early during evolution in both types of plants. It is as yet unclear, whether there are any major functional differences between the dicotyledonous and monocotyledonous HD-Zip IV proteins.
The cDNA of HD-Zip IV gene GLABRA2-like9 (TaGL9) was isolated from a wheat grain cDNA library based on strong interaction of the gene product with the bait sequence 5′-GCATTAAATGC-3′, which contains the L1 box [45]. To have a better understanding on how the HD region of HD-Zip IV proteins interacts with the L1 box, a homology molecular model of the TaGL9 HD binding the 5′-GCATTAAATGC-3′ element was constructed (Figure 2). Chain A of HD from mouse Hox-A9 along with a 20 bp-DNA duplex [50], were used as templates to generate the model. A multiple sequence alignment was performed using ProMals3D [51], and selected HDs from HD-Zip IV subfamily members (Figure 2A) including TaGL9 HD were used to assess conservation of amino acid residues in HD of TaGL9. A total of 138 publicly available amino acid sequences with HDs were found through ConSurf analysis [52] by CSI-BLAST [53] and 66 unique amino acid sequences of the HD-Zip IV proteins including TaGL9 from Triticum aestivum (Acc. No. AEI99592) were selected for analysis [52,54]. The ConSurf analysis revealed that among the three α-helices of HD, helix 3 is the region with the highest level of conservation (Figure 2B). This may be because the amino acid residues in the third helix (here designated as the recognition helix) are responsible for most interactions between the TaGL9 HD and the specific DNA element, while the other two helices are most probably responsible for the formation of the hydrophobic core, enabling the recognition helix to achieve an optimal position for formation of the protein-DNA complex [55].
The results of HD sequence analysis by ProMals3D software (Figure 2A) and the model of the TaGL9 HD, which was assigned according to the conservation values of the ConSurf Server (Figure 2B), revealed a group of five amino acid residues that are highly conserved within HDs of all known HD-Zip IV proteins. These five amino acid residues are responsible for the interaction between the HD-Zip IV protein and the 5′-GCATTAAATGC-3′ cis-element. The residues contacting DNA are indicated by black dotted lines in Figure 2C. Four of the DNA-interacting residues (Trp86, Gln88, Asn89 and Arg91, Figure 2C) were found to be located within the recognition helix and had an evolutionary conservation score of 9.9. One other amino acid residue (Arg41) was positioned at the N-terminal part of HD of TaGL9 and had an evolutionary conservation score of 8. Notably, one amino acid residue (Met93) located in the recognition helix had a relatively low evolutionary conservation score of 6 that indicates a high variability of this residue across HD-Zip IV members (Figure 2B). Nevertheless, results of the ConSurf analysis demonstrated that this Met93 can only be substituted by residues with hydrophobic side chains such as tryptophan, leucine and isoleucine. These data indicate that the specific residue in this position may be important for maintaining the structural integrity of HD or might be involved in protein-protein interaction, as residues with hydrophobic side chains are rarely exposed to a hydrophilic environment.
Comparisons of HD sequences of HD-Zip IV proteins with those of other HD-Zip classes using the homology model of TaGL9 HD revealed a low level of similarity between these sequences (28% sequence identity between either HD-Zip I or II and HD-Zip IV, and 29% sequence identity between HD-Zip III and IV). Further analysis using multiple sequence alignment generated with CLUSTALX [47] showed that the third helix of HD defines the identity of HD domains in HD-Zip TFs from different subfamilies (data not shown). These results indicated that the other two-helices, which contain highly variable amino acid residues, might be responsible for recognition of subfamily-specific cis-elements.
The second conserved domain characteristic of HD-Zip IV TFs is a leucine zipper (Zip), which is specific for the HD-Zip IV structure and in this subfamily of HD-ZIP TFs it has been described as a bipartite dimerization leucine zipper-loop-zipper (ZLZ) motif [31]. The motif is located immediately after the third helix of HD. Despite disruption of the leucine zipper of the ZLZ motif by an inserted loop sequence, dimerization via this motif can still occur [31]. Dimerisation is a prerequisite for DNA binding by HD-Zip IV proteins. Unlike other HD proteins, which can bind DNA as monomers, HD-Zip IV TFs can bind DNA efficiently only as either homo- or hetero-dimers [58].
The third domain which is present in all HD-Zip IV proteins is a START domain, which is composed of approximately 200 amino acids residues. In animals, the START domain has been found in proteins belonging to several families including the phosphatidylcholine transfer protein (PCTP), metastatic lymph node 64 (MLN64) and the StAR-related lipid transfer protein 4 (STARD4), but it has not been found to be associated with a DNA binding domain of animal TFs [13,59]. The main role of the START domain in mammals is in binding and transportation of lipid and/or sterol molecules. Approximately three quarters of Arabidopsis (26/35) and rice (22/29) proteins with a START domain contain additional protein domain(s) [13]. Among the additional domains are HD, Plekstrin Homology (PH) domain and a domain 1336 (DUF1336) with an as yet unknown function. Similarity of structural features of the START domain-containing HD-Zip IV proteins, and mammalian nuclear receptors containing DNA-binding Zinc fingers and sterol/lipid binding Ligand Binding Domain (LBD) suggest that the START domain in plants might be involved in signal transduction and direct regulation of transcription by binding and transporting steroid-type phytohormones and/or other lipid molecules [13]. The binding affinity of TFs containing a START domain to specific DNA elements may be affected either by a direct protein-lipid/sterol interaction or by an interaction of lipid/sterol with a partner protein, which is bound to the same promoter region. The function of the HD-Zip IV protein HD-GLABRA2 (HD-GL2) from Arabidopsis as a negative regulator in phospholipid signalling in roots supports the possible involvement of the START domain in such regulation [26]. It was suggested that binding of lipids by the START domain might be necessary for nuclear transport of HD-Zip IV TFs in a way similar to the model proposed for nuclear transport of the glucocorticoid receptors in mammals [60]. A nuclear localisation of the HD-GL2 protein and the maize HD-Zip IV protein ZmOCL1, have been demonstrated [61,62]. In addition, it has been shown that AD of the maize HD-Zip IV protein ZmOCL1 is located at the NH2-terminal region of the START domain [62].
One more distinguishable feature of HD-Zip IV proteins is a conserved HD-START-associated domain (HD-SAD) [10]. Little is known about functions of this domain. However, simultaneous presence of both the START and HD-SAD domains of cotton HD-Zip IV TF GbML1 is required for binding to the C-terminal domain of GbMYB25, a key regulator of cotton fibre initiation. This was demonstrated in the EMSA and the yeast 2-hybrid (Y2H) assay using truncated proteins and proteins with point mutations [63].
Using the maize Outer Cell Layer 1 (ZmOCL1) protein as bait in a Y2H screen, the SWITCH complex protein 3C1 (SWI3C1) was identified and was confirmed as an interacting partner in planta[62]. However, the domain/motif of ZmOCL1 responsible for this interaction was not mapped. SWI3C1 is a member of the SWI3 protein family and belongs to the C subfamily of such proteins in plants. It possesses all domain characteristics which are found in similar components of the SWITCH/SUCROSE NONFERMENTING (SWI/SNF) ATP-dependent chromatin re-modelling complex of yeast and mammalian origins [64]. ATP-dependent chromatin re-modelling complexes, such as SWI/SNF, temporarily disable histone interactions with DNA, allowing accessibility of DNA to TFs. However, the SWI/SNF complex binds DNA in a non-specific manner. It was shown that the recruitment of this complex in mammals and yeast requires interaction with transcriptional activators [65,66]. In maize, it is likely that ZmOCL1 specifies recruitment of the SWI/SNF complex to a particular promoter region(s), and that its association with this complex facilitates access to specific cis-element(s) [40].
The HD-Zip IV proteins are among the longest TFs in plants and have a complex domain structure. Unfortunately, attempts to understand function of these TFs through structural studies have so far produced little data.

3. Regulation of HD-Zip IV Gene Expression

Our understanding of cell-specific regulation of HD-Zip IV gene expression is very limited. In most reports expression of HD-Zip IV genes has been observed in outer cell layers, preferentially in the epidermis of developing embryo and/or other plant organs [24,40], including the endosperm, scutellum and roots [40]. The L1 box is a DNA element which directs specific gene expression in epidermal cells [25]. The promoter sequence of only one HD-Zip IV gene, AtML1, has so far been studied [67]. AtML1 expression is initially observed in the apical daughter cell of the zygote and later in epidermal cells throughout embryo and plant development [68]. A 101 bp-long sequence containing the L1 box and a putative WUSCHEL-binding site was identified using a GFP reporter gene [67]. This sequence was sufficient for embryo-specific activation of the AtML1 promoter. Interestingly, the L1 box alone was not able to activate expression of the reporter gene. However, several other promoter fragments could activate the AtML1 promoter in some but not in all epidermal cells. It was concluded that the AtML1 promoter has a complex modular structure and that several regions of the promoter are involved in spatial and temporal regulation of AtML1 expression during embryogenesis [67]. The presence of the L1 box was identified in several other HD-Zip IV promoters that might suggest feed-back regulation of HD-Zip IV genes by their own gene products or by other members of the subfamily [25,26].
Attempts were made to explain cell-specific expression of the HD-Zip genes by microRNA (miRNA) involvement. Most miRNAs found in Arabidopsis are targeted to TF genes, which are involved in plant development and particularly in development of shoot and floral meristems. One of the well characterised cases of miRNA regulation in plants is recognition of a conserved region at the beginning of the START domain in the HD-Zip III subfamily of TFs, by miR165/166 [69]. Because of tight high similarity of domain structures between HD-Zip III and HD-Zip IV proteins, this regulatory system may also be applied to HD-Zip IV genes. Regulation of HD-Zip IV genes by miRNAs might be supported by the presence of a conserved 17-nucleotide motif located at in the 3′UTR region of the HD-Zip IV genes. This conserved sequence was found in PaHB1, PaHB2 and several other HD-Zip IV genes from Arabidopsis, sunflower and cotton [35]. However, alternative mechanisms, which could involve the 17 bp long motif in regulation of gene expression, may also exist. Some metazoan HB genes expressed during embryogenesis contain regulatory elements located in their 3′-UTR regions that influence subcellular localisation, translation and stability of mRNAs [35]. Therefore, further investigation is required to test whether the conserved motif in HD-Zip IV genes has any regulatory function involving miRNAs.
Expression of the FWA/HDG6 gene was found to be controlled by epigenetic hyper-methylation of two direct repeats situated at the 5′ region of the gene around its transcription start site [70]. Ectopic FWA expression led to a late-flowering phenotype, which resulted from an in inhibition of function of the FLOWERING Locus T. (FT) gene by the FWA gene product in vegetative tissue of Arabidopsis[7072]. It was demonstrated that FWA expression is imprinted in endosperm [30,73]. FWA is silenced by DNA methylation in vegetative tissue but it is de-methylated in the central cell of the ovule and after fertilisation FWA is expressed in the endosperm from the maternal copy of the gene [74].

4. Roles of HD-Zip IV TFs

4.1. Development and Maintenance of Epidermal Cell Layers

Expression of genes encoding the HD-Zip subfamily IV proteins can typically be observed in the outer cells of embryos and in the epidermal or sub-epidermal layers of developing plant organs, where regulation of epidermal cell fate, trichome formation and anthocyanin accumulation occur [32,33,75]. The epidermis plays a critical role in plant defence against pathogens and in protection from environmental stresses [23]. Differentiation of the epidermis is one of the major events in plant development. During embryogenesis the basic plant body plan is realised by specifying position of shoot and root meristems. The meristems contain groups of stem cells, which provide initial material for the development of plant organs after seed germination. During meristem formation protoderm differentiation leads to organisation of inner and outer cell layers [76]. Epidermal-specific expression patterns of the HD-Zip IV genes have been observed in a variety of plant species including Arabidopsis[24], rice (Oryza sativa) [41], maize (Zea mays) [38], cotton (Gossypium hirsutum) [36] and Norway Spruce (Picea abies) [34].
A total of sixteen HD genes, (HDGs) namely HDG1HDG6/FWA, HDG7HDG12, GLABRA2 (GL2), ANTHOCYANINLESS2 (ANL2), AtML1, and PDF2, comprise the Arabidopsis HD-Zip IV subfamily [24]. At least some of these genes were shown to be involved in cell fate determination of the epidermal cell layer by regulating cell layer-specific gene expression (Table 1). GL2 or AtHB-10 from Arabidopsis was the first identified member of the HD-Zip IV subfamily. The expression of GL2 was detected only in the specific epidermal cells responsible for formation of trichomes and root hairless epidermal cells. Absence of trichomes was observed in the gl2 mutant suggesting that the GL2 gene may play a role in trichome formation [77]. In addition, the gl2 mutation caused root-hair formation from hairless epidermal cells, but other cellular processes that usually take place during differentiation of root-hairs and hairless cell files were not affected [78]. Arabidopsis HDG11 and HDG12 genes were reported to be involved in trichome branching. The hdg11, but not the hdg12, mutant had more branched trichomes than wild-type plants. However, the phenotype of the hdg11 mutant was enhanced in double mutants of hdg11 and hdg12[24]. Surprisingly, the characteristic phenotype of gl2 mutant was not altered in the gl2 hdg11 and gl2 hdg12 double mutants [24].
Involvement of HD-Zip IV genes in differentiation of the epidermis has been clearly demonstrated using mutants of two other HD-Zip IV genes from Arabidopsis[33]. Single gene knockout mutants of PDF2 and ATML1, which have a high level of sequence identity, retained normal shoot development, while the double mutant had severe defects in shoot epidermal cell differentiation. This suggests that PDF2 and ATML1 have the same or very similar functions in differentiation of epidermal cells, which are possibly realised by interacting with L1 boxes of downstream target-gene promoters [33].

4.2. Development of Sub-Epidermal Cell Layers and Control of Anthocyanin Pigmentation

The product of the HD-Zip IV gene ANL2 from Arabidopsis shares a high level of amino acid sequence similarity with the product of the AtML1 gene, but has different function. ANL2 controls anthocyanin pigmentation of the leaf sub-epidermal layer and cellular organisation of the primary root [32]. Plant anthocyanins are pigments that are responsible for a range of colours such as red, blue and purple. The accumulation of anthocyanins occurs during specific plant developmental stages and only in certain plant tissues. The accumulation of anthocyanins is also induced by environmental stimuli such as light, temperature, nutrients and by different stresses [32]. It is interesting that ANL2 controls pigmentation of leaves and stems, but is not involved in regulation of anthocyanin synthesis in the internal layer of the seed coat. The anl2 mutant has an abnormal radial root patterning. This mutant produces extra cells known as intervening cells located between the cortical and epidermal layers [32]. Based on these data it was proposed that ANL2 might be involved in the maintenance of sub-epidermal layer identity during plant development and that it controls the process of anthocyanin accumulation in sub-epidermal cells.

4.3. Drought Tolerance

The up-regulation of the Arabidopsis HD-Zip IV gene HDG11 was found to contribute to increased drought tolerance [17]. Such observation originates from characterisation of the enhanced drought tolerance 1 (edt1) mutant, which demonstrated better drought tolerance characteristics than control plants. The mutant had dramatically increased expression levels of the HDG11 gene. The mutant phenotype included improved root architecture, a more extensive root system, and reduced stomatal density. The mutants also had higher levels of abscisic acid, which further expanded the capacity for stress response. In addition higher levels of Superoxide Dismutase (SOD) activity and Pro in mutant plants could potentially contribute to better osmotic adjustment and Radical Oxygen Species (ROS) detoxification. The same phenotype was observed when HDG11 was constitutively over-expressed in wild-type Arabidopsis. The authors suggested that such elevated levels of expression resulted in significantly increased amounts of HDG11 protein in the nucleus, which made it a dominant TF that successfully competed for any available HD-binding site in the genome [17]. No changes in epidermal cell layers or cuticle development were reported [17]. However, higher levels of ABA in the edt1 mutant suggested that a mechanism of drought response regulation by HDG11 might be similar to that of the members of the HD-Zip I subfamily [21,22].

4.4. Function of HD-Zip IV Genes in Commercially Important Plant Species

Data obtained with Arabidopsis on the regulation of trichome development by HD-Zip IV genes, stimulated attempts to understand roles of these genes in agriculturally important plants with the aim to apply the gained knowledge for improvement of important traits. Cotton (Gossypium spp.) fibres, which are also known as seed trichomes, have a high economic value. Therefore, an understanding of the molecular mechanisms responsible for fibre initiation and control of development is important for enhancing yield and quality of fibres [37]. Cotton fibres are distinct from the Arabidopsis trichomes. They are exceptionally elongated, usually not branched and when mature, composed of almost pure cellulose [36]. The Gossypium arboretum L. HOMEOBOX genes 1 and 2 (GaHOX1 and GaHOX2) are cotton genes encoding HD-Zip IV TFs. GaHOX1 was found to be the closest homologue of GL2, while GaHOX2 was more similar to AtML1 and PDF2 from Arabidopsis. Functional complementation with the GaHOX1 gene was observed when it was expressed under the GL2 promoter in the Arabidopsis gl2 mutant. This suggested that the GaHOX1 protein has the same function as GL2 in regulation of trichome development. In contrast, GaHOX2 might be involved in seed coat development rather than in trichome growth regulation [36].
Expression levels of the HD-Zip IV gene from Gossypium hirsutum L., known as HOMEODOMAIN 1 (GhHD-1) were lower in the ovule epidermis of a fibreless mutant compared to those expressed in the ovule epidermis of wild-type plants [79]. GhHD-1 was predominantly expressed in epidermal and sub-epidermal layers of cells. The fact that this gene was up-regulated during conversion of cells into fibres might suggest its involvement in fibre differentiation [37]. A close homologue of GhHD-1 from G. barbadense (GbML1) has been shown to interact with GhMYB25 in the Y2H assay. The GhMYB25 protein in turn was demonstrated to be involved in the activation of the fibre-preferential cotton gene RD22-like1 (RDL1) by direct binding to its promoter region, and hence this gene also has a role in fibre development [63]. Because a large number of MYBs are expressed in cotton ovules, the authors speculated that GhHD-1 might interact with various MYB factors in different cell types to enable diverse functions [80]. Phenotyping of transgenic plants with silenced GhMYB25 and GhHD-1 genes indicated that these genes work in cooperation. Both types of transgenic plants had reduced trichome numbers and a delayed initiation of seed fibre differentiation. However, some differences in function of GhMYB25 and GhHD-1 were also found. For example, a decrease of fibre length was observed only in lines with silenced GhMYB25 but not in the GhHD-1-silenced lines [37,81].
A total of seventeen HD-Zip IV genes, also known as Outer Cell Layer (OCL) genes, have been identified in the maize genome. However, only a few of them have been characterised in detail (Table 1). Most are preferentially expressed in immature reproductive organs and/or the epidermal (L1) cell layer [40]. It was shown that ZmOCL1, ZmOCL3, ZmOCL4 and ZmOCL5 have epidermis-specific expression patterns, whereas expression of ZmOCL2 was restricted to the sub-epidermal layer (L2) [38,82]. It was demonstrated that ZmOCL4 is involved in anther and trichome (macrohairs in maize) development [39]. The ocl4 mutant has ectopic growth of macrohairs. Surprisingly, expression of ZmOCL4 did not rescue the gl2 mutant in Arabidopsis, and even enhanced the influence of the gl2 mutation on trichome development. The ocl4 mutant also has partial male sterility. Smaller size of anthers in mutant maize plants compared to the wild-type plants and absence of anthocyanins suggested a premature termination of anther development in the mutant [39].
Another maize gene, ZmOCL1, was found to have a role in regulation of kernel development in maize. A reduction in kernel size was observed in transgenic maize plants transformed with a construct for expression of a fusion of the repressor domain of the Drosophila Engrailed protein and the dimerisation domains of ZmOCL1, which converted ZmOCL1 from an activator to a repressor [7].
Analysis of the rice HD-Zip IV genes revealed that while some of these are involved in regulation of processes that take place in the epidermis, the function of others seems to be different. Nine GL2-type Rice outermost cell-specific (Roc) genes were found in the rice genome using a BLAST search [42]. Five of these genes were isolated as full length cDNA clones and designated Roc1 to Roc5 (Table 1). Expression of all five genes was restricted to the leaf epidermis [41,42]. Specific expression of Roc1 begins in the outermost cells, and precedes protoderm differentiation. The authors suggested two possible explanations for this phenomenon: (1) epidermal cell fate is dependent on Roc1 expression; or (2) Roc1 expression occurs before the determination of epidermal cell fate [41,42]. Expression of Roc1 was also detected in the outermost cell layer of callus before organogenesis, and was independent of whether the cells did or did not differentiate. Wounding by cutting has induced Roc1 expression in newly formed outer layers. These observations suggested that Roc1 is expressed in a position-dependent manner and that Roc1 is not involved in determination of fate of epidermal cells [41].
Another rice gene, Roc5, has been recently shown to be associated with development and formation of bulliform cells in rice leaves [83]. The bulliform cells act as specialised epidermal cells that are involved in protection against pathogen invasions and in prevention of water loss during stress. These cells are found on the adaxial surface of leaves in all monocotyledous plants except Helobiae[84,85]. Roc5 overexpression and suppression in transgenic rice plants were used to determine how Roc5 affects development of bulliform cells in rice. Overexpression of the Roc5 gene resulted in a decrease in bulliform cell number and size that in turn lead to the adaxial rolling of leaves. On the contrary, suppression of Roc5 in transgenic rice led to an increased number and size of bulliform cells and the abaxial rolling of leaves. These data indicate that Roc5 is a negative regulator of bulliform cell development and bulliform cell-dependent leaf rolling [83]. It would be also interesting to investigate the role of Roc2 in bulliform cell development since physical interaction of Roc5 and Roc2 has been demonstrated to occur in yeast [42].
Rice Transcription Factor 1 (OsTF1) is a rice HD-Zip IV gene, which shares a relatively low sequence similarity with AtML1 and ZmOCL5[86]. HD of OsTF1 is located at the N-terminal part of the protein, unlike HDs in most other HD-Zip IV members, which are situated either in the middle or at the C-terminal part of the protein [86]. However, OsTF1 contains all domains and motifs that are characteristic of the HD-Zip IV proteins [68]. Regardless of the low level of similarity between protein sequences of OsTF1 and AtML1 or ZmOCL5, the epidermis-specific expression pattern of OsTF1 in the early embryo is consistent with the expression pattern of other HD-Zip IV genes. The precise role of the OsTF1 gene remains unclear, but it may contribute to regulation of organ or tissue differentiation in the developing embryo. High levels of expression of this gene were observed shortly before organ differentiation occurred [86].
The closest wheat homologue of OsTF1, TaGL9, was found to be expressed in embryo epidermis and adjacent to the embryo portion of the syncytial endosperm during early grain development. However, later in development TaGL9 expression relocated to the main vascular bundle of the scutellum. Such patterns of expression suggest that TaGL9 could play a role in development and maintenance of inner surface layers of the vascular bundle [45]. A gene encoding a lipid transfer protein with the same pattern of expression has been identified in wheat [87]. This gene may be a potential target gene of TaGL9, or it may be that expression of both genes is regulated in a similar way by other TF(s).

4.5. Regulation of Cuticle Formation and Plant Protection

Recently, a new function of HD-Zip IV TFs has been described in the context of regulation of biosynthesis and transport of plant cuticle components. Plant cuticle has an important role in establishment of a protective hydrophobic layer on the external cell wall of epidermal cells in aerial parts of plants [88]. The cuticle layer provides protection to plants from biotic and abiotic stresses through its specific chemical, physical and optical properties [23,89,90]. Thickness, chemical composition and structure of a cuticular layer, as well as environmental parameters such as temperature at the leaf surface, are important determinants of cuticle permeability [89,91,92].
The regulation of activity of genes involved in cuticle biosynthesis is still poorly understood. Members of two families of TFs have been demonstrated to control cuticle biosynthesis and accumulation of cuticular wax during plant development and during plant responses to stresses. These TFs are the ethylene responsive factors (ERF-type transcriptional regulators) WAX INDUCER1/SHINE1 (WIN1/SHN1) and WAX PRODUCTION1 (WPX1), and several MYB TFs [9397]. Recently, a new regulatory pathway controlling cuticle development was reported. It involves a Curly Flag Leaf1 (CFL1) gene, which encodes a protein containing a WW domain [29]. As the name suggests, the WW domain contains two conserved tryptophans, which can bind proline-rich peptides [98]. Overexpression of both rice and Arabidopsis CFL1 in transgenic Arabidopsis plants showed severely impaired cuticle development. It has been demonstrated that AtCFL1 interacts directly with the HD-Zip IV protein HDG1, which negatively regulates its activity [29]. It was also recently demonstrated that HDG1 contains an N-terminal Ethylene-responsive element binding factor-associated Amphiphilic Repression (EAR) motif [99]. Transcriptional activators with attached EAR motifs usually function as dominant repressors [100]. However, over-expression of HDG1 in wild-type Arabidopsis did not produce a phenotype similar to those observed either during over-expression of AtCFL1 or repression of HDG1. It was speculated that the six amino acid residue EAR motif of HDG1 is too short to act as a dominant repressor without assistance from other cofactors [99].
Recently, it has been shown that transgenic maize with the overexpressed HD-Zip IV gene ZmOCL1 had altered wax composition in juvenile leaves [101]. This finding was supported by the fact that ZmOCL1 regulates expression of several genes involved in lipid transport and metabolism, and potentially could be involved in cuticle biosynthesis and deposition [101]. Based on microarray analysis, 14 genes were identified to be directly or indirectly regulated by over-expression of the ZmOCL1 gene. Half of these genes encode proteins that are known to be involved in biosynthesis and/or transport of lipid molecules and plant defence. Among them are the type 2 lipid transfer protein (nsLTPII), an AtCXE-18 carboxylesterase, a SEC14/phosphatidylinositol transport protein (PITP), three ABC transporters (WBC11/ABCG11-like) and a fatty acid reductase (FAR). FAR is closely related to the Arabidopsis ECERIFERUM 4 (CER4) protein, which is responsible for synthesis of long chain primary alcohols from fatty acid precursors [102].
The major function of carboxylesterase enzymes is to hydrolyse esters of short-chain fatty acids [103]. However, the majority of carboxylesterases are associated with plant defence [104]. The SEC14/PITP transport proteins, similar to the one activated by ZmOCL1, are known to be involved in trans-Golgi export pathways [105] and therefore can be potentially involved in biosynthesis of cuticular lipids. SEC14/PITP can catalyse phosphatidylinositol and phosphatidylcholine transfer in vitro. In mammals, several diseases develop due to deficiencies of PITPs [106].
As mentioned above, another group of genes activated by ZmOCL1 encode ABC transporters. The ABC transporters belonging to the White Brown Complex (WBC) family are specialised in ATP-dependent translocation of steroids and other lipid molecules in mammals [107]. The analysis of Arabidopsis mutants has shown that the ECERIFERUM 5 (CER5) protein WBC12/ABCG12 and CER5-like protein WBC11/ABCG11 are key components of cutin and wax transport from their site of biosynthesis to the cuticle layer [108,109]. It is known that CER5 is a stress-responsive gene and that it is up-regulated by the stress hormone ABA, drought and high salinity [90].
The gene encoding type 2 lipid transfer protein (nsLTPII) is among genes regulated by HD-Zip IV TF ZmOCL1. It has been proposed that LTPs play an essential role in transport of cuticular lipids through plant cell walls [110]. Recently the role of LTPs in the export of cuticular components was confirmed by analysis of Arabidopsis plants with a mutation in the LTPG1 gene. The ltpg1 mutant had a thinner wax layer; this finding suggested a role of LTPG1 in facilitation of the export of wax components through the plasma membrane [111,112]. Some NsLTPs have been shown to have a role in plant protection against fungal infection and insects [113,114]. It is also known that over-expression of HvLTP2 in barley increases tolerance to Pseudomonas syringae[115], while the ltpg1 mutant showed increased susceptibility to fungal infections [112].
Plant defencins (γ-thyonines) comprise an LTP subfamily known to have strong antifungal and antibacterial properties [116118]. A group of wheat genes encoding defensins was found to be specifically expressed in the outer cell layers of developing grain [119]. It has been recently demonstrated that a wounding-inducible HD-Zip IV TF, which is preferentially expressed in wheat grain, can activate promoters of several grain-specific defensins in a transient expression assay (our unpublished results). This findings support data [101] on activation of the LTP gene by ZmOCL1 over-expression, and suggests that HD-Zip IV TFs can be involved in plant protection from pathogens either directly, by activation of certain pathogenesis-related proteins, or indirectly, by activation of proteins that are involved in cuticle formation.

5. Conclusions and Future Perspectives

The existing data regarding HD-Zip IV TFs suggest their involvement in regulation of development and maintenance of outer cell layers of plant organs, which may play an important role in plant defence from pathogens and protection from abiotic stresses. Emerging data imply that the HD-Zip IV proteins regulate cuticle biosynthesis and wax accumulation. This function bears similarity to functions of earlier discovered TFs from the ERF and MYB families. The inducibility of cuticle-related ERF and MYB TFs by ABA, drought and high salinity has been reported and ability of these TFs to confer stress tolerance to transgenic Arabidopsis, alfalfa and rice plants has been demonstrated [9396]. It would be interesting to examine if any of the HD-Zip IV genes are involved in abiotic stress responses by modulation of size and content of cuticle.
Also fascinating is the fact that some HD-Zip IV genes were found to be induced by wounding and could regulate some pathogen responsive genes. Further studies are required to determine if ‘clever’ up-regulation of HD-Zip IV genes in transgenic plants can directly or indirectly enhance plant defence. Such ‘clever’ up-regulation would entail the use of outer cell layer-specific promoters, or modification of the promoters of HD-Zip IV genes themselves by the addition of enhancer sequences.
In contrast to Arabidopsis, where many HD-Zip IV genes are expressed in vegetative tissues, in grasses these genes are expressed mostly in reproductive parts of the plant. A role for ZmOCL1 in maize kernel development has been demonstrated [7]. Further comprehensive analysis of function needs to be done for other cereal HD-Zip IV family members, which are predominantly expressed in flowers and grain.
The mechanism of HD-Zip IV promoter regulation remains unclear and requires further attention. Nevertheless, HD-Zip IV genes are a promising source of tissue-specific promoters of moderate strength, particularly for specific expression in epidermal and/or sub-epidermal layers, as well as in specific parts of flowers and grain.

Acknowledgments

William Chew thanks to the Australian Centre for Plant Functional Genomics for financial support. This work was supported by grant DP120100900 from the Australian Research Council to MH, and by grant LP120100201 from the Australian Research Council to MH and SL. We wish to thank Dinesh Yadav, Nikolai Borisjuk and Julie Hayes for critically reading the manuscript.

Conflict of Interest

The authors declare no conflict of interest.

References

  1. Wang, R.L.; Stec, A.; Hey, J.; Lukens, L.; Doebley, J. The limits of selection during maize domestication. Nature 1999, 398, 236–239. [Google Scholar]
  2. Bustamante, C.D.; Fledel-Alon, A.; Williamson, S.; Nielsen, R.; Hubisz, M.T.; Glanowski, S.; Tanenbaum, D.M.; White, T.J.; Sninsky, J.J.; Hernandez, R.D.; et al. Natural selection on protein-coding genes in the human genome. Nature 2005, 437, 1153–1157. [Google Scholar]
  3. Ma, J.; Ptashne, M. A new class of yeast transcriptional activators. Cell 1987, 51, 113–119. [Google Scholar]
  4. Garber, R.L.; Kuroiwa, A.; Gehring, W.J. Genomic and cDNA clones of the homeotic locus antennapedia in Drosophila. EMBO J 1983, 2, 2027–2036. [Google Scholar]
  5. Gehring, W.J.; Affolter, M.; Burglin, T. Homeodomain proteins. Annu. Rev. Biochem 1994, 63, 487–526. [Google Scholar]
  6. Chan, R.L.; Gago, G.M.; Palena, C.M.; Gonzalez, D.H. Homeoboxes in plant development. Biochim. Biophys. Acta 1998, 23, 1–19. [Google Scholar]
  7. Khaled, A.S.; Vernoud, V.; Ingram, G.C.; Perez, P.; Sarda, X.; Rogowsky, P.M. Engrailed-ZmOCL1 fusions cause a transient reduction of kernel size in maize. Plant Mol. Biol 2005, 58, 123–139. [Google Scholar]
  8. Vollbrecht, E.; Veit, B.; Sinha, N.; Hake, S. The developmental gene Knotted-1 is a member of a maize homeobox gene family. Nature 1991, 350, 241–243. [Google Scholar]
  9. Ariel, F.D.; Manavella, P.A.; Dezar, C.A.; Chan, R.L. The true story of the HD-Zip family. Trends Plant Sci 2007, 12, 419–426. [Google Scholar]
  10. Mukherjee, K.; Brocchieri, L.; Bürglin, T.R. A comprehensive classification and evolutionary analysis of plant homeobox genes. Mol. Biol. Evol 2009, 26, 2775–2794. [Google Scholar]
  11. 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–1791. [Google Scholar]
  12. Riechmann, J.L.; Heard, J.; Martin, G.; Reuber, L.; Jiang, C.; Keddie, J.; Adam, L.; Pineda, O.; Ratcliffe, O.J.; Samaha, R.R.; et al. Arabidopsis transcription factors: Genome-wide comparative analysis among eukaryotes. Science 2000, 290, 2105–2110. [Google Scholar]
  13. Schrick, K.; Nguyen, D.; Karlowski, W.M.; Mayer, K.F. START lipid/sterol-binding domains are amplified in plants and are predominantly associated with homeodomain transcription factors. Genome Biol. 2004, 5. [Google Scholar] [CrossRef] [Green Version]
  14. Mukherjee, K.; Bürglin, T.R. MEKHLA, a novel domain with similarity to PAS domains, is fused to plant homeodomain-leucine zipper III proteins. Plant Physiol 2006, 140, 1142–1150. [Google Scholar]
  15. Sessa, G.; Morelli, G.; Ruberti, I. The Athb-1 and -2 HD-Zip domains homodimerize forming complexes of different DNA binding specificities. EMBO J 1993, 12, 3507–3517. [Google Scholar]
  16. Harris, J.C.; Hrmova, M.; Lopato, S.; Langridge, P. Modulation of plant growth by HD-Zip class I and II transcription factors in response to environmental stimuli. New Phytol 2011, 190, 823–837. [Google Scholar]
  17. Yu, H.; Chen, X.; Hong, Y.-Y.; Wang, Y.; Xu, P.; Ke, S.-D.; Liu, H.-Y.; Zhu, J.-K.; Oliver, D.J.; Xiang, C.-B. Activated expression of an Arabidopsis HD-START protein confers drought tolerance with improved root system and reduced stomatal density. Plant Cell 2008, 20, 1134–1151. [Google Scholar]
  18. Zhang, S.; Haider, I.; Kohlen, W.; Jiang, L.; Bouwmeester, H.; Meijer, A.; Schluepmann, H.; Liu, C.-M.; Ouwerkerk, P.F. Function of the HD-Zip I gene Oshox22 in ABA-mediated drought and salt tolerances in rice. Plant Mol. Biol 2012, 80, 571–585. [Google Scholar]
  19. Cabello, J.V.; Chan, R.L. The homologous homeodomain-leucine zipper transcription factors HaHB1 and AtHB13 confer tolerance to drought and salinity stresses via the induction of proteins that stabilize membranes. Plant Biotech. J 2012, 10, 815–825. [Google Scholar]
  20. Nadakuduti, S.S.; Pollard, M.; Kosma, D.K.; Allen, C., Jr; Ohlrogge, J.B.; Barry, CS. Pleiotropic phenotypes of the sticky peel mutant provide new insight into the role of CUTIN DEFICIENT2 in epidermal cell function in tomato. Plant Physiol 2012, 159, 945–960. [Google Scholar]
  21. Valdés, A.; Övernäs, E.; Johansson, H.; Rada-Iglesias, A.; Engström, P. The homeodomain-leucine zipper (HD-Zip) class I transcription factors ATHB7 and ATHB12 modulate abscisic acid signalling by regulating protein phosphatase 2C and abscisic acid receptor gene activities. Plant Mol. Biol 2012, 80, 405–418. [Google Scholar]
  22. Lechner, E.; Leonhardt, N.; Eisler, H.; Parmentier, Y.; Alioua, M.; Jacquet, H.; Leung, J.; Genschik, P. MATH/BTB CRL3 receptors target the homeodomain-leucine zipper ATHB6 to modulate abscisic acid signaling. Dev. Cell 2011, 21, 1116–1128. [Google Scholar]
  23. Javelle, M.; Vernoud, V.; Rogowsky, P.M.; Ingram, G.C. Epidermis: The formation and functions of a fundamental plant tissue. New Phytol 2011, 189, 17–39. [Google Scholar]
  24. Nakamura, M.; Katsumata, H.; Abe, M.; Yabe, N.; Komeda, Y.; Yamamoto, K.T.; Takahashi, T. Characterization of the class IV homeodomain-leucine zipper gene family in Arabidopsis. Plant Physiol 2006, 141, 1363–7375. [Google Scholar]
  25. Abe, M.; Takahashi, T.; Komeda, Y. Identification of a cis-regulatory element for L1 layer-specific gene expression, which is targeted by an L1-specific homeodomain protein. Plant J 2001, 26, 487–494. [Google Scholar]
  26. Ohashi, Y.; Oka, A.; Rodrigues-Pousada, R.; Possenti, M.; Ruberti, I.; Morelli, G.; Aoyama, T. Modulation of phospholipid signaling by GLABRA2 in root-hair pattern formation. Science 2003, 300, 1427–1430. [Google Scholar]
  27. Graba, Y.; Aragnol, D.; Pradel, J. Drosophila Hox complex downstream targets and the function of homeotic genes. BioEssays 1997, 19, 379–388. [Google Scholar]
  28. Ekker, S.C.; von Kessler, D.P.; Beachy, P.A. Differential DNA sequence recognition is a determinant of specificity in homeotic gene action. EMBO J 1992, 11, 4059–4072. [Google Scholar]
  29. Wu, R.; Li, S.; He, S.; Wassmann, F.; Yu, C.; Qin, G.; Schreiber, L.; Qu, L.J.; Gu, H. CFL1, a WW domain protein, regulates cuticle development by modulating the function of HDG1, a class IV homeodomain transcription factor, in rice and Arabidopsis. Plant Cell 2011, 23, 3392–3411. [Google Scholar]
  30. Kinoshita, T.; Miura, A.; Choi, Y.; Kinoshita, Y.; Cao, X.; Jacobsen, S.E.; Fischer, R.L.; Kakutani, T. One-Way Control of FWA Imprinting in Arabidopsis Endosperm by DNA Methylation. Science 2004, 303, 521–523. [Google Scholar]
  31. Di Cristina, M.; Sessa, G.; Dolan, L.; Linstead, P.; Baima, S.; Ruberti, I.; Morelli, G. The Arabidopsis Athb-10 (GLABRA2) is an HD-Zip protein required for regulation of root hair development. Plant J 1996, 10, 393–402. [Google Scholar]
  32. Kubo, H.; Peeters, A.J.; Aarts, M.G.; Pereira, A.; Koornneef, M. ANTHOCYANINLESS2, a homeobox gene affecting anthocyanin distribution and root development in Arabidopsis. Plant Cell 1999, 11, 1217–1226. [Google Scholar]
  33. 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–643. [Google Scholar]
  34. Ingouff, M.; Farbos, I.; Lagercrantz, U.; von Arnold, S. PaHB1 is an evolutionary conserved HD-GL2 homeobox gene expressed in the protoderm during Norway spruce embryo development. Genesis 2001, 30, 220–230. [Google Scholar]
  35. Ingouff, M.; Farbos, I.; Wiweger, M.; von Arnold, S. The molecular characterization of PaHB2, a homeobox gene of the HD-GL2 family expressed during embryo development in Norway spruce. J. Exp. Bot 2003, 54, 1343–1350. [Google Scholar]
  36. Guan, X.Y.; Li, Q.J.; Shan, C.M.; Wang, S.; Mao, Y.B.; Wang, L.J.; Chen, X.Y. The HD-Zip IV gene GaHOX1 from cotton is a functional homologue of the Arabidopsis GLABRA2. Plant Physiol 2008, 134, 174–182. [Google Scholar]
  37. Walford, S.A.; Wu, Y.; Llewellyn, D.J.; Dennis, E.S. Epidermal cell differentiation in cotton mediated by the homeodomain leucine zipper gene, GhHD-1. Plant J 2012, 71, 464–478. [Google Scholar]
  38. Ingram, G.C.; Magnard, J.L.; Vergne, P.; Dumas, C.; Rogowsky, P.M. Expression patterns of genes encoding HD-Zip IV homeodomain proteins define specific domains in maize embryos and meristems. Plant J 2000, 22, 401–414. [Google Scholar]
  39. Vernoud, V.; Laigle, G.; Rozier, F.; Meeley, R.B.; Perez, P.; Rogowsky, P.M. The HD-ZIP IV transcription factor OCL4 is necessary for trichome patterning and anther development in maize. Plant J 2009, 59, 883–894. [Google Scholar]
  40. Javelle, M.; Klein-Cosson, C.; Vernoud, V.; Boltz, V.; Maher, C.; Timmermans, M.; Depege-Fargeix, N.; Rogowsky, P.M. Genome-wide characterization of the HD-ZIP IV transcription factor family in maize: Preferential expression in the epidermis. Plant Physiol 2011, 157, 790–803. [Google Scholar]
  41. Ito, M.; Sentoku, N.; Nishimura, A.; Hong, S.-K.; Sato, Y.; Matsuoka, M. Position dependent expression of GL2-type homeobox gene, Roc1: Significance for protoderm differentiation and radial pattern formation in early rice embryogenesis. Plant J 2002, 29, 497–507. [Google Scholar]
  42. Ito, M.; Sentoku, N.; Nishimura, A.; Hong, S.-K.; Sato, Y.; Matsuoka, M. Roles of rice GL-2 type homeobox genes in epidermis differentiation. Breed. Sci 2003, 53, 245–253. [Google Scholar]
  43. Sasaki, T.; Matsumoto, T.; Yamamoto, K.; Sakata, K.; Baba, T.; Katayose, Y.; Wu, J.; Niimura, Y.; Cheng, Z.; Nagamura, Y.; et al. The genome sequence and structure of rice chromosome 1. Nature 2002, 420, 312–316. [Google Scholar]
  44. Yu, J.; Wang, J.; Lin, W.; Li, S.; Li, H.; Zhou, J.; Ni, P.; Dong, W.; Hu, S.; Zeng, C.; et al. The genomes of Oryza sativa: A history of duplications. PLoS Biol 2005, 3, e38. [Google Scholar]
  45. Kovalchuk, N.; Wu, W.; Eini, O.; Bazanova, N.; Pallotta, M.; Shirley, N.; Singh, R.; Ismagul, A.; Eliby, S.; Johnson, A.; et al. The scutellar vascular bundle–specific promoter of the wheat HD-Zip IV transcription factor shows similar spatial and temporal activity in transgenic wheat, barley and rice. Plant Biotech. J 2012, 10, 43–53. [Google Scholar]
  46. Prigge, M.J.; Clark, S.E. Evolution of the class III HD-Zip gene family in land plants. Evol. Dev 2006, 8, 350–361. [Google Scholar]
  47. Thompson, J.D.; Gibson, T.J.; Plewniak, F.; Jeanmougin, F.; Higgins, D.G. The CLUSTAL_X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 1997, 25, 4876–4882. [Google Scholar]
  48. Saitou, N.; Nei, M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol. Biol. Evol 1987, 4, 406–425. [Google Scholar]
  49. Schultz, J.; Milpetz, F.; Bork, P.; Ponting, C.P. SMART, a simple modular architecture research tool: Identification of signaling domains. Proc. Natl. Acad. Sci. USA 1998, 95, 5857–5864. [Google Scholar]
  50. LaRonde-LeBlanc, N.A.; Wolberger, C. Structure of HoxA9 and Pbx1 bound to DNA: Hox hexapeptide and DNA recognition anterior to posterior. Genes Dev 2003, 17, 2060–2072. [Google Scholar]
  51. Pei, J.; Kim, B.H.; Grishin, N.V. PROMALS3D: A tool for multiple protein sequence and structure alignments. Nucleic Acids Res 2008, 36, 2295–2300. [Google Scholar]
  52. Landau, M.; Mayrose, I.; Rosenberg, Y.; Glaser, F.; Martz, E.; Pupko, T.; Ben-Tal, N. ConSurf 2005: The projection of evolutionary conservation scores of residues on protein structures. Nucleic Acids Res 2005, 33, W299–W302. [Google Scholar]
  53. Biegert, A.; Soding, J. Sequence context-specific profiles for homology searching. Proc. Natl. Acad. Sci. USA 2009, 106, 3770–3775. [Google Scholar]
  54. Ashkenazy, H.; Erez, E.; Martz, E.; Pupko, T.; Ben-Tal, N. ConSurf 2010: Calculating evolutionary conservation in sequence and structure of proteins and nucleic acids. Nucleic Acids Res 2010, 38, W529–W533. [Google Scholar]
  55. Tron, A.E.; Bertoncini, C.W.; Palena, C.M.; Chan, R.L.; Gonzalez, D.H. Combinatorial interactions of two amino acids with a single base pair define target site specificity in plant dimeric homeodomain proteins. Nucleic Acids Res 2001, 29, 4866–4872. [Google Scholar]
  56. Bru, C.; Courcelle, E.; Carrère, S.; Beausse, Y.; Dalmar, S.; Kahn, D. The ProDom database of protein domain families: More emphasis on 3D. Nucleic Acids Res 2005, 33, D212–D215. [Google Scholar]
  57. Cotsaftis, O.; Plett, D.; Shirley, N.; Tester, M.; Hrmova, M. A two-staged model of Na+ exclusion in rice explained by 3D modeling of HKT transporters and alternative splicing. PLoS One 2012, 7, e39865. [Google Scholar]
  58. Palena, C.M.; Tron, A.E.; Bertoncini, C.W.; Gonzalez, D.H.; Chan, R.L. Positively charged residues at the N-terminal arm of the homeodomain are required for efficient DNA binding by homeodomain-leucine zipper proteins. J. Mol. Biol 2001, 308, 39–47. [Google Scholar]
  59. Ponting, C.P.; Aravind, L. START: A lipid-binding domain in StAR, HD-ZIP and signalling proteins. Trends Biochem. Sci 1999, 24, 130–132. [Google Scholar]
  60. Kumar, R.; Thompson, E.B. Gene regulation by the glucocorticoid receptor: Structure:Function relationship. J. Steroid Biochem. Mol. Biol 2005, 94, 383–394. [Google Scholar]
  61. Szymanski, D.B.; Jilk, R.A.; Pollock, S.M.; Marks, M.D. Control of GL2 expression in Arabidopsis leaves and trichomes. Development 1998, 125, 1161–1171. [Google Scholar]
  62. Depege-Fargeix, N.; Javelle, M.; Chambrier, P.; Frangne, N.; Gerentes, D.; Perez, P.; Rogowsky, P.M.; Vernoud, V. Functional characterization of the HD-ZIP IV transcription factor OCL1 from maize. J. Exp. Bot 2011, 62, 293–305. [Google Scholar]
  63. Zhang, F.; Zuo, K.; Zhang, J.; Liu, X.; Zhang, L.; Sun, X.; Tang, K. An L1 box binding protein, GbML1, interacts with GbMYB25 to control cotton fibre development. J. Exp. Bot 2010, 61, 3599–3613. [Google Scholar]
  64. Sarnowski, T.J.; Ríos, G.; Jásik, J.; Œwieżewskia, S.; Kaczanowskia, S.; Lib, Y.; Kwiatkowskac, A.; Pawlikowskaa, K.; KoŸbiałc, M.; KoŸbiałc, P.; et al. SWI3 subunits of putative SWI/SNF chromatin-remodeling complexes play distinct roles during Arabidopsis development. Plant Cell 2005, 17, 2454–2472. [Google Scholar]
  65. Neely, K.E.; Hassan, A.H.; Brown, C.E.; Howe, L.; Workman, J.L. Transcription activator interactions with multiple SWI/SNF subunits. Mol. Cell Biol 2002, 22, 1615–1625. [Google Scholar]
  66. Simone, C. SWI/SNF: The crossroads where extracellular signaling pathways meet chromatin. J. Cell Physiol 2006, 207, 309–314. [Google Scholar]
  67. Takada, S.; Jurgens, G. Transcriptional regulation of epidermal cell fate in the Arabidopsis embryo. Development 2007, 134, 1141–1150. [Google Scholar]
  68. Lu, P.; Porat, R.; Nadeau, J.A.; O’Neill, S.D. Identification of a meristem L1 layer-specific gene in Arabidopsis that is expressed during embryonic pattern formation and defines a new class of homeobox genes. Plant Cell 1996, 8, 2155–2168. [Google Scholar]
  69. Bowman, J.L. Class III HD-Zip gene regulation, the golden fleece of ARGONAUTE activity? BioEssays 2004, 26, 938–942. [Google Scholar]
  70. Soppe, W.J.; Jacobsen, S.E.; Alonso-Blanco, C.; Jackson, J.P.; Kakutani, T.; Koornneef, M.; Peeters, A.J. The late flowering phenotype of FWA mutants is caused by gain-of-function epigenetic alleles of a homeodomain gene. Mol. Cell 2000, 6, 791–802. [Google Scholar]
  71. Kakutani, T. Genetic characterization of late-flowering traits induced by DNA hypomethylation mutation in Arabidopsis thaliana. Plant J 1997, 12, 1447–1451. [Google Scholar]
  72. Ikeda, Y.; Kobayashi, Y.; Yamaguchi, A.; Abe, M.; Araki, T. Molecular basis of late-flowering phenotype caused by dominant epi-alleles of the FWA locus in Arabidopsis. Plant Cell Physiol 2007, 48, 205–220. [Google Scholar]
  73. Fujimoto, R.; Kinoshita, Y.; Kawabe, A.; Kinoshita, T.; Takashima, K.; Nordborg, M.; Nasrallah, M.E.; Shimizu, K.K.; Kudoh, H.; Kakutani, T. Evolution and control of imprinted FWA genes in the genus Arabidopsis. PLoS Genet 2008, 4, e1000048. [Google Scholar]
  74. Fujimoto, R.; Sasaki, T.; Kudoh, H.; Taylor, J.M.; Kakutani, T.; Dennis, E.S. Epigenetic variation in the FWA gene within the genus Arabidopsis. Plant J 2011, 66, 831–843. [Google Scholar]
  75. Prigge, M.J.; Otsuga, D.; Alonso, J.M.; Ecker, J.R.; Drews, G.N.; Clark, S.E. Class III homeodomain-leucine zipper gene family members have overlapping, antagonistic, and distinct roles in Arabidopsis development. Plant Cell 2005, 17, 61–76. [Google Scholar]
  76. Bommert, P.; Werr, W. Gene expression patterns in the maize caryopsis: Clues to decisions in embryo and endosperm development. Gene 2001, 271, 131–142. [Google Scholar]
  77. Hülskamp, M.; Miséra, S.; Jürgens, G. Genetic dissection of trichome cell development in Arabidopsis. Cell 1994, 76, 555–566. [Google Scholar]
  78. Masucci, J.D.; Rerie, W.G.; Foreman, D.R.; Zhang, M.; Galway, M.E.; Marks, M.D.; Schiefelbein, J.W. The homeobox gene GLABRA2 is required for position-dependent cell differentiation in the root epidermis of Arabidopsis thaliana. Development 1996, 122, 1253–1260. [Google Scholar]
  79. Wu, Y.; Machado, A.C.; White, R.G.; Llewellyn, D.J.; Dennis, E.S. Expression profiling identifies genes expressed early during lint fibre initiation in cotton. Plant Cell Physiol 2006, 47, 107–127. [Google Scholar]
  80. Suo, J.; Liang, X.; Pu, L.; Zhang, Y.; Xue, Y. Identification of GhMYB109 encoding a R2R3 MYB transcription factor that expressed specifically in fiber initials and elongating fibers of cotton (Gossypium hirsutum L.). Biochim. Biophys. Acta 2003, 20, 25–34. [Google Scholar]
  81. Machado, A.; Wu, Y.; Yang, Y.; Llewellyn, D.J.; Dennis, E.S. The MYB transcription factor GhMYB25 regulates early fibre and trichome development. Plant J 2009, 59, 52–62. [Google Scholar]
  82. Ingram, G.C.; Magnard, J.L.; Vergne, P.; Dumas, C.; Rogowsky, P.M. ZmOCL1, an HD-GL2 family homeobox gene, is expressed in the outer cell layer throughout maize development. Plant Mol. Biol 1999, 40, 343–354. [Google Scholar]
  83. Zou, L.-P.; Sun, X.-H.; Zhang, Z.-G.; Liu, P.; Wu, J.-X.; Tian, C.-J.; Qiu, J.-L.; Lu, T.-G. Leaf rolling controlled by the homeodomain leucine zipper class IV gene Roc5 in rice. Plant Physiol 2011, 156, 1589–1602. [Google Scholar]
  84. Metcalfe, C.R. Anatomy of the Monocotyledons; Clarendon Press: Oxford, UK, 1960. [Google Scholar]
  85. Jane, W.; Chiang, S. Morphology and development of bulliform cells in arundo formosana hack. Taiwania 1991, 36, 85–97. [Google Scholar]
  86. Yang, J.Y.; Chung, M.C.; Tu, C.Y.; Leu, W.M. OSTF1: A HD-GL2 family homeobox gene is developmentally regulated during early embryogenesis in rice. Plant Cell Physiol 2002, 43, 628–638. [Google Scholar]
  87. Boutrot, F.; Meynard, D.; Guiderdoni, E.; Joudrier, P.; Gautier, M. The Triticum aestivum non-specific lipid transfer protein (TaLtp) gene family: Comparative promoter activity of six TaLtp genes in transgenic rice. Planta 2007, 225, 843–862. [Google Scholar]
  88. Jeffree, C.E. The Fine Structure of the Plant Cuticle. In Annual Plant Reviews Volume 23: Biology of the Plant Cuticle; Blackwell Publishing Ltd: Oxford, UK, 2007; pp. 11–125. [Google Scholar]
  89. Kerstiens, G. Water transport in plant cuticles: An update. J. Exp. Bot 2006, 57, 2493–2499. [Google Scholar]
  90. Kosma, D.; Jenks, M. Eco-Physiological and Molecular-Genetic Determinants of Plant Cuticle Function in Drought and Salt Stress Tolerance. In Advances in Molecular Breeding Toward Drought and Salt Tolerant Crops; Jenks, M., Hasegawa, P., Jain, S.M., Eds.; Springer: Dordrecht, The Netherlands, 2007; pp. 91–120. [Google Scholar]
  91. Burghardt, M.; Riederer, M. 9 Cuticular transpiration. In Annual Plant Reviews, Biololgy of the Plant Cuticle 23; Blackwell Publishing: Oxford, UK, 2006; Volume 292, pp. 292–312. [Google Scholar]
  92. Shepherd, T.; Wynne Griffiths, D. The effects of stress on plant cuticular waxes. New Phytol 2006, 171, 469–499. [Google Scholar]
  93. Broun, P.; Poindexter, P.; Osborne, E.; Jiang, C.-Z.; Riechmann, J.L. WIN1, a transcriptional activator of epidermal wax accumulation in Arabidopsis. Proc. Natl. Acad. Sci. USA 2004, 101, 4706–4711. [Google Scholar]
  94. Aharoni, A.; Dixit, S.; Jetter, R.; Thoenes, E.; van Arkel, G.; Pereira, A. The SHINE clade of ap2 domain transcription factors activates wax biosynthesis, alters cuticle properties, and confers drought tolerance when overexpressed in Arabidopsis. Plant Cell 2004, 16, 2463–2480. [Google Scholar]
  95. Seo, P.J.; Lee, S.B.; Suh, M.C.; Park, M.-J.; Go, Y.S.; Park, C.-M. The MYB96 transcription factor regulates cuticular wax biosynthesis under drought conditions in Arabidopsis. Plant Cell 2011, 23, 1138–1152. [Google Scholar]
  96. Zhang, J.Y.; Broeckling, C.D.; Sumner, L.W.; Wang, Z.Y. Heterologous expression of two Medicago truncatula putative ERF transcription factor genes, WXP1 and WXP2, in Arabidopsis led to increased leaf wax accumulation and improved drought tolerance, but differential response in freezing tolerance. Plant Mol. Biol 2007, 64, 265–278. [Google Scholar]
  97. Cominelli, E.; Galbiati, M.; Tonelli, C. Integration of water stress response: Cell expansion and cuticle deposition in Arabidopsis thaliana. Plant Signal Behav 2008, 3, 556–557. [Google Scholar]
  98. Chen, H.I.; Sudol, M. The WW domain of Yes-associated protein binds a proline-rich ligand that differs from the consensus established for Src homology 3-binding modules. Proc. Natl. Acad. Sci. USA 1995, 92, 7819–7823. [Google Scholar]
  99. Kagale, S.; Links, M.G.; Rozwadowski, K. Genome-wide analysis of ethylene-responsive element binding factor-associated amphiphilic repression motif-containing transcriptional regulators in Arabidopsis. Plant Physiol 2010, 152, 1109–1134. [Google Scholar]
  100. Hiratsu, K.; Matsui, K.; Koyama, T.; Ohme-Takagi, M. Dominant repression of target genes by chimeric repressors that include the EAR motif, a repression domain, in Arabidopsis. Plant J 2003, 34, 733–739. [Google Scholar]
  101. Javelle, M.; Vernoud, V.; Depege-Fargeix, N.; Arnould, C.; Oursel, D.; Domergue, F.; Sarda, X.; Rogowsky, P.M. Overexpression of the epidermis-specific homeodomain-leucine zipper IV transcription factor Outer Cell Layer1 in maize identifies target genes involved in lipid metabolism and cuticle biosynthesis. Plant Physiol 2010, 154, 273–286. [Google Scholar]
  102. Rowland, O.; Zheng, H.; Hepworth, S.R.; Lam, P.; Jetter, R.; Kunst, L. CER4 encodes an alcohol-forming fatty acyl-coenzyme a reductase involved in cuticular wax production in Arabidopsis. Plant Physiol 2006, 142, 866–877. [Google Scholar]
  103. Cummins, I.; Landrum, M.; Steel, P.; Edwards, R. Structure activity studies with xenobiotic substrates using carboxylesterases isolated from Arabidopsis thaliana. Phytochemistry 2007, 68, 811–818. [Google Scholar]
  104. Frey, M.; Stettner, C.; Paré, P.W.; Schmelz, E.A.; Tumlinson, J.H.; Gierl, A. An herbivore elicitor activates the gene for indole emission in maize. Proc. Natl. Acad. Sci. USA 2000, 97, 14801–14806. [Google Scholar]
  105. Curwin, A.J.; Fairn, G.D.; McMaster, C.R. Phospholipid transfer protein Sec14 is required for trafficking from endosomes and regulates distinct trans-Golgi export pathways. J. Biol. Chem 2009, 284, 7364–7375. [Google Scholar]
  106. Bankaitis, V.A.; Phillips, S.; Yanagisawa, L.; Li, X.; Routt, S.; Xie, Z. Phosphatidylinositol transfer protein function in the yeast Saccharomyces cerevisiae. Adv. Enzyme Regul 2005, 45, 155–170. [Google Scholar]
  107. Velamakanni, S.; Wei, S.L.; Janvilisri, T.; van Veen, H.W. ABCG transporters: Structure, substrate specificities and physiological roles: A brief overview. J. Bioenerg. Biomembr 2007, 39, 465–471. [Google Scholar]
  108. Pighin, J.A.; Zheng, H.; Balakshin, L.J.; Goodman, I.P.; Western, T.L.; Jetter, R.; Kunst, L.; Samuels, A.L. Plant cuticular lipid export requires an ABC transporter. Science 2004, 306, 702–704. [Google Scholar]
  109. Bird, D.; Beisson, F.; Brigham, A.; Shin, J.; Greer, S.; Jetter, R.; Kunst, L.; Wu, X.; Yephremov, A.; Samuels, L. Characterization of Arabidopsis ABCG11/WBC11, an ATP binding cassette (ABC) transporter that is required for cuticular lipid secretion. Plant J 2007, 52, 485–498. [Google Scholar]
  110. Kader, J. Lipid-transfer proteins in plants. Annu. Rev. Plant Physiol. Plant Mol. Biol 1996, 47, 627–654. [Google Scholar]
  111. DeBono, A.; Yeats, T.H.; Rose, J.K.C.; Bird, D.; Jetter, R.; Kunst, L.; Samuels, L. Arabidopsis LTPG is a glycosylphosphatidylinositol-anchored lipid transfer protein required for export of lipids to the plant surface. Plant Cell 2009, 21, 1230–1238. [Google Scholar]
  112. Lee, S.B.; Go, Y.S.; Bae, H.J.; Park, J.H.; Cho, S.H.; Cho, H.J.; Lee, D.S.; Park, O.K.; Hwang, I.; Suh, M.C. Disruption of glycosylphosphatidylinositol-anchored lipid transfer protein gene altered cuticular lipid composition, increased plastoglobules, and enhanced susceptibility to infection by the fungal pathogen Alternaria brassicicola. Plant Physiol 2009, 150, 42–54. [Google Scholar]
  113. Boutrot, F.; Chantret, N.; Gautier, M.-F. Genome-wide analysis of the rice and Arabidopsis non-specific lipid transfer protein (nsLtp) gene families and identification of wheat nsLtp genes by EST data mining. BMC Genomics 2008, 9, 86, :1–86:19.. [Google Scholar]
  114. Zottich, U.; da Cunha, M.; Carvalho, A.O.; Dias, G.B.; Silva, N.C.; Santos, I.S.; do Nacimento, V.V.; Miguel, E.C.; Machado, O.L.; Gomes, V.M. Purification, biochemical characterization and antifungal activity of a new lipid transfer protein (LTP) from Coffea canephora seeds with alpha-amylase inhibitor properties. Biochim. Biophys. Acta 2011, 4, 375–383. [Google Scholar]
  115. Molina, A.; Garcia-Olmedo, F. Enhanced tolerance to bacterial pathogens caused by the transgenic expression of barley lipid transfer protein LTP2. Plant J 1997, 12, 669–675. [Google Scholar]
  116. Lay, F.T.; Mills, G.D.; Poon, I.K.; Cowieson, N.P.; Kirby, N.; Baxter, A.A.; van der Weerden, N.L.; Dogovski, C.; Perugini, M.A.; Anderson, M.A.; et al. Dimerization of plant defensin NaD1 enhances its antifungal activity. J. Biol. Chem 2012, 287, 19961–19972. [Google Scholar]
  117. Rigano, M.M.; Romanelli, A.; Fulgione, A.; Nocerino, N.; D’Agostino, N.; Avitabile, C.; Frusciante, L.; Barone, A.; Capuano, F.; Capparelli, R. A novel synthetic peptide from a tomato defensin exhibits antibacterial activities against Helicobacter pylori. J. Pept. Sci 2012, 18, 755–762. [Google Scholar]
  118. Picart, P.; Pirttila, A.; Raventos, D.; Kristensen, H.-H.; Sahl, H.-G. Identification of defensin-encoding genes of Picea glauca: Characterization of PgD5, a conserved spruce defensin with strong antifungal activity. BMC Plant Biol 2012, 12, 180, :1–180:14.. [Google Scholar]
  119. Kovalchuk, N.; Li, M.; Wittek, F.; Reid, N.; Singh, R.; Shirley, N.; Ismagul, A.; Eliby, S.; Johnson, A.; Milligan, A.S.; et al. Defensin promoters as potential tools for engineering disease resistance in cereal grains. Plant Biotech. J 2010, 8, 47–64. [Google Scholar]
Figure 1. Unrooted phylogenetic tree of selected HD-Zip IV proteins. HD-Zip IV protein sequences were retrieved from the NCBI database and aligned with CLUSTALX [47]. Unrooted phylogenetic trees were constructed based on aligned protein sequences using the Neighbor-Joining algorithm [48] with a Bootstrap value of 1000 from CLUSTALX [47]. Species of origin are indicated by two-letter prefixes. The accession numbers of the published protein sequences used in the phylogenetic trees are listed in Table 1; unpublished BnBBIP-1A has Accession ABA54874. (A) Phylogenetic tree based on full-length amino acid sequences of 43 HD-Zip IV proteins; (B) Phylogenetic tree constructed using HDs of selected HD-Zip IV proteins (60–61 amino acid residues). The HD sequences that were included in the analyses were selected by a Simple Modular Architecture Research Tool (SMART) [49]. At, Arabidopsis thaliana; Pa, Picea abies; Ga, Gossypium spp; Gh, Gossypium hirsutum; Zm, Zea mays; Os, Oryza sativa; Ta, Triticum aestivum; Td, Triticum durum; Sl, Solanum lycopersicum; Bn, Brassica napus; Pp, Physcomitrella patens; Sm, Selaginella moellendorffii. Proteins from mono- and dicotyledonous groups are indicated in bold and normal types, respectively. Picea abies, Physcomitrella patens and Selaginella moellendorffii have not been assigned to either group (underlined).
Figure 1. Unrooted phylogenetic tree of selected HD-Zip IV proteins. HD-Zip IV protein sequences were retrieved from the NCBI database and aligned with CLUSTALX [47]. Unrooted phylogenetic trees were constructed based on aligned protein sequences using the Neighbor-Joining algorithm [48] with a Bootstrap value of 1000 from CLUSTALX [47]. Species of origin are indicated by two-letter prefixes. The accession numbers of the published protein sequences used in the phylogenetic trees are listed in Table 1; unpublished BnBBIP-1A has Accession ABA54874. (A) Phylogenetic tree based on full-length amino acid sequences of 43 HD-Zip IV proteins; (B) Phylogenetic tree constructed using HDs of selected HD-Zip IV proteins (60–61 amino acid residues). The HD sequences that were included in the analyses were selected by a Simple Modular Architecture Research Tool (SMART) [49]. At, Arabidopsis thaliana; Pa, Picea abies; Ga, Gossypium spp; Gh, Gossypium hirsutum; Zm, Zea mays; Os, Oryza sativa; Ta, Triticum aestivum; Td, Triticum durum; Sl, Solanum lycopersicum; Bn, Brassica napus; Pp, Physcomitrella patens; Sm, Selaginella moellendorffii. Proteins from mono- and dicotyledonous groups are indicated in bold and normal types, respectively. Picea abies, Physcomitrella patens and Selaginella moellendorffii have not been assigned to either group (underlined).
Ijms 14 08122f1
Figure 2. Molecular modelling of the homeobox domain (HD) of TaGL9 from wheat in complex with an 11-bp long DNA fragment. TaGL9 has at least six domains as predicted by ProDom [56]. (A) A multiple sequence alignment of selected HD sequences using ProMals3D [51]. The predicted secondary structures are shown in magenta (α-helices) and black (loops). Conservation of residues on a scale of 9–5 is shown at the top of the diagram. The absolutely conserved and similar residues are shaded in green and yellow, respectively. The black box indicates the boundaries of HD domains. Vertical arrows above the alignment point to the DNA-interacting residues shown in panel C; (B) HD structure of TaGL9 showing the degree of conservation. The coloured model is based on known HD-Zip IV protein sequences found in the Consurf database [52]. 138 sequences were found by CSI-BLAST [53] but only 66 sequences were unique for the Consurf algorithm to perform calculation. Highly conserved amino acids are coloured in deep magenta, while the least conserved and average ones are coloured in cyan and white, respectively. The black arrow indicates the third helix, for which there is a particularly high level of amino acid conservation; (C) A molecular model of HD of TaGL9 in complex with 5′-GCATTAAATGC-3′/3′-GCATTTAATGC-5′; the model was constructed as described [57]. HD of Hox-A9 from mouse (Protein Data Bank accession 1puf, chain A), in complex with a 20-bp DNA duplex fragment 5′-ACTCTATGATTTACGACGCT-3′ [50] served as a template. Ribbon representation in salmon shows the disposition of secondary structure elements. Here, the α-helix 3 (perpendicular to the viewer’s plane) carries most of the residues that mediate contacts between HD and the DNA fragment. The duplex DNA is shown in cpk-green (coding strand) and cpk-yellow (complementary strand). The nucleotides interacting with HD are represented as cpk sticks. The left and right black arrows point to the NH2- and COOH-termini of HD, respectively. Separations of ≥3.5 Å between the contacting residues (1-letter codes) of HD and a DNA strand are indicated by black dotted lines. The interplay of the interacting residues within HD suggests that structural rigidity and/or flexibility could impact upon selectivity of DNA binding. It is of note that mainly the TAAA and GCAT segments of the coding and complementary strands, respectively, are interacting with the five highly conserved residues of HD. The nucleotide sequence of the DNA coding strand is shown on the left.
Figure 2. Molecular modelling of the homeobox domain (HD) of TaGL9 from wheat in complex with an 11-bp long DNA fragment. TaGL9 has at least six domains as predicted by ProDom [56]. (A) A multiple sequence alignment of selected HD sequences using ProMals3D [51]. The predicted secondary structures are shown in magenta (α-helices) and black (loops). Conservation of residues on a scale of 9–5 is shown at the top of the diagram. The absolutely conserved and similar residues are shaded in green and yellow, respectively. The black box indicates the boundaries of HD domains. Vertical arrows above the alignment point to the DNA-interacting residues shown in panel C; (B) HD structure of TaGL9 showing the degree of conservation. The coloured model is based on known HD-Zip IV protein sequences found in the Consurf database [52]. 138 sequences were found by CSI-BLAST [53] but only 66 sequences were unique for the Consurf algorithm to perform calculation. Highly conserved amino acids are coloured in deep magenta, while the least conserved and average ones are coloured in cyan and white, respectively. The black arrow indicates the third helix, for which there is a particularly high level of amino acid conservation; (C) A molecular model of HD of TaGL9 in complex with 5′-GCATTAAATGC-3′/3′-GCATTTAATGC-5′; the model was constructed as described [57]. HD of Hox-A9 from mouse (Protein Data Bank accession 1puf, chain A), in complex with a 20-bp DNA duplex fragment 5′-ACTCTATGATTTACGACGCT-3′ [50] served as a template. Ribbon representation in salmon shows the disposition of secondary structure elements. Here, the α-helix 3 (perpendicular to the viewer’s plane) carries most of the residues that mediate contacts between HD and the DNA fragment. The duplex DNA is shown in cpk-green (coding strand) and cpk-yellow (complementary strand). The nucleotides interacting with HD are represented as cpk sticks. The left and right black arrows point to the NH2- and COOH-termini of HD, respectively. Separations of ≥3.5 Å between the contacting residues (1-letter codes) of HD and a DNA strand are indicated by black dotted lines. The interplay of the interacting residues within HD suggests that structural rigidity and/or flexibility could impact upon selectivity of DNA binding. It is of note that mainly the TAAA and GCAT segments of the coding and complementary strands, respectively, are interacting with the five highly conserved residues of HD. The nucleotide sequence of the DNA coding strand is shown on the left.
Ijms 14 08122f2aIjms 14 08122f2b
Table 1. Published HD-Zip IV genes and their characteristics.
Table 1. Published HD-Zip IV genes and their characteristics.
PlantGene nameAccessionTissue localisationSuggested function*Cited
Arabidopsis thalianaHDG1NP_191674In trichomes at the base of young leaf, endodermal cell lines around emergent lateral root, stamen filamentCuticle development[24,29]
Arabidopsis thalianaHDG2NP_172015Outer cell layers of shoot apical meristems, early flower primordia, nucellus, seed coat, epidermal layers of young leaves, hairless cell files of hypocotyl epidermis, primary root tipsEmbryo development[24]
Arabidopsis thalianaHDG3NP_180796Siliques and seedlingsCotyledon development[24]
Arabidopsis thalianaHDG4NP_193506FlowersFlower organ development[24]
Arabidopsis thalianaHDG5NP_199499All plant organs except root, outer cell layers of shoot apical meristems, early flower primordia, nucellus, epidermis of the stamen filament, stomatal guard cells of the carpelNot determined[24]
Arabidopsis thalianaFWA/HDG6NP_567722EndospermEndosperm development[30]
Arabidopsis thalianaHDG7NP_200030Base of leaf primordia, apical region of the heart-stage embryo, lateral root primordia and tips, seedlingsNot determined[24]
Arabidopsis thalianaHDG8NP_186976Endosperm and embryo at early stages of development, after flower organ developmentFlower, seed development[24]
Arabidopsis thalianaHDG9NP_197234Developed flowers, chalazal, embryo sacEmbryonic development[24]
Arabidopsis thalianaHDG12NP_564041Apical meristem, stamen filament, nucellus, embryo, young epidermal tissues, lateral root tipTrichome development[24]
Arabidopsis thalianaGL2NP_565223Epidermis of leaves and roots, developing trichomes and surrounding epidermal cells, hairless cell files of hypocotyl epidermisTrichome and root hair development[31]
Arabidopsis thalianaANL2NP_567183Leaves, stems, buds, rootsAnthocyanin distribution and root development[32]
Arabidopsis thalianaATML1NP_193906Flower bud, lower parts of siliquesShoot epidermal cell differentiation[33]
Arabidopsis thalianaPDF2NP_567274Flower buds, shoot apicesShoot epidermal cell differentiation, cotyledon development[33]
Picea abiesPaHB1AAG43405EmbryoEmbryo development[34]
Picea abiesPaHB2AAL83725EmbryoEmbryo development[35]
Gossypium spp.GaHOX1ABY41242Fibre cellsCotton fibre development[36]
Gossypium spp.GaHOX2ABY67263Fibre cells at early developmental stagesSeed coat development[36]
Gossypium hirsutumGhHD-1AAFO11041Epidermal cellsEpidermal cell determination[37]
Gossypium hirsutumGhHD-1DAFO11042Epidermal cellsEpidermal cell determination[37]
Zea maysZmOCL1CAG38614Immature ears, epidermis specific expression in embryo, meristems, young organ primordia, juvenile leavesKernel development, cuticle deposition[7]
Zea maysZmOCL2CAB96422Apical meristem, sub-epidermal layer of floral meristems immature ears, young female gametophyte, kernels, rootNot determined[38]
Zea maysZmOCL3CAB96423Developing embryo and endosperm, epidermis specific expression in embryo, meristems, young organ primordial, juvenile leavesSpecification of organ identity[38]
Zea maysZmOCL4CAB96424.2Shoot epidermis, meristems and young organ primordia, immature tasselsAnther and trichome development[39]
Zea maysZmOCL5CAB96425Immature tassels, epidermis specific expression in embryo, meristems, young organ primordia,Not determined[40]
Zea maysZmOCL6DAA34955Immature tassels, immature ears, epidermal cells of juvenile leavesNot determined[40]
Zea maysZmOCL7DAA34956Immature ears, epidermal cellsNot determined[40]
Zea maysZmOCL8DAA34957Immature tassels, epidermal cellsNot determined[40]
Zea maysZmOCL9DAA34958Young developing kernels, L2 cells of shoot apical meristemNot determined[40]
Zea maysZmOCL10DAA34959Immature tassels, epidermal layer of juvenile leavesNot determined[40]
Zea maysZmOCL11DAA34960Mature tassels that include pollenNot determined[40]
Zea maysZmOCL12DAA34961Not determinedNot determined[40]
Zea maysZmOCL13DAA34962Immature tassels, epidermis of shoot apical meristemNot determined[40]
Zea maysZmOCL14DAA34963Immature tassels, epidermal cellsNot determined[40]
Zea maysZmOCL15DAA34964Immature tassels, epidermis of shoot apical meristemNot determined[40]
Zea maysZmOCL16DAA34965Young developing kernels, epidermal cellsNot determined[40]
Zea maysZmOCL17DAA34966Immature tassels, epidermis of shoot apical meristemNot determined[40]
Oryza sativaOsRoc1BAB85750Shoot apex, inflorescence apex, leaf bladeEmbryogenesis[41]
Oryza sativaOsRoc2BAC77155Shoot apex, inflorescence apex, leaf bladeNot determined[42]
Oryza sativaOsRoc3BAC77156Shoot apex, inflorescenceNot determined[42]
Oryza sativaOsRoc4BAC77157RootsNot determined[42]
Oryza sativaOsRoc5BAC77158Leaf epidermisEpidermal differentiation, bulliform cells development[42]
Oryza sativaOsROC6Q7Y0V7Not determined[42]
Oryza sativaOsRoc7BAC77160EpidermisEpidermal differentiation[42]
Oryza sativaOsRoc8BAC77161EpidermisEpidermal differentiation[42]
Oryza sativaOsRoc9BAC77162EpidermisEpidermal differentiation[42]
Oryza sativaOsTF1Q5ZAY0GrainEmbryogenesis[43]
Oryza sativaOsHOX9A2Z8L4Not determined[44]
Triticum aestivumTaGL9AEI99592Grain specific expressionNot determined[45]
Triticum durumTdGL9HAEI99593Early embryo, endosperm around embryo, later in scutellar vascular bundleMaintenance of scutellar vascular bundle[45]
Solanum lycopersicumSlCD2NP_001234657Epidermal cell and cuticleEpidermal cell and cuticle development[20]
Physcomitrella patensPpHDZ41DAA05775Not determinedNot determined[46]
Selaginella moellendorffiiSmHDZ44DAA05774Not determinedNot determined[46]
*It should be taken in consideration that the current data on HD-Zip IV function are dependent on the available data, and thus in some instances function of HD-Zip IV TFs could be speculative.

Share and Cite

MDPI and ACS Style

Chew, W.; Hrmova, M.; Lopato, S. Role of Homeodomain Leucine Zipper (HD-Zip) IV Transcription Factors in Plant Development and Plant Protection from Deleterious Environmental Factors. Int. J. Mol. Sci. 2013, 14, 8122-8147. https://doi.org/10.3390/ijms14048122

AMA Style

Chew W, Hrmova M, Lopato S. Role of Homeodomain Leucine Zipper (HD-Zip) IV Transcription Factors in Plant Development and Plant Protection from Deleterious Environmental Factors. International Journal of Molecular Sciences. 2013; 14(4):8122-8147. https://doi.org/10.3390/ijms14048122

Chicago/Turabian Style

Chew, William, Maria Hrmova, and Sergiy Lopato. 2013. "Role of Homeodomain Leucine Zipper (HD-Zip) IV Transcription Factors in Plant Development and Plant Protection from Deleterious Environmental Factors" International Journal of Molecular Sciences 14, no. 4: 8122-8147. https://doi.org/10.3390/ijms14048122

APA Style

Chew, W., Hrmova, M., & Lopato, S. (2013). Role of Homeodomain Leucine Zipper (HD-Zip) IV Transcription Factors in Plant Development and Plant Protection from Deleterious Environmental Factors. International Journal of Molecular Sciences, 14(4), 8122-8147. https://doi.org/10.3390/ijms14048122

Article Metrics

Back to TopTop