Homeodomain-Leucine Zipper II family of transcription factors to the limelight: central regulators of plant development

Genome-wide characterisation of HD-Zip transcription factors and functional analysis of PbHB24 during stone cell formation in Chinese white pear (Pyrus bretschneideri)

BACKGROUND

The homodomain-leucine zipper (HD-Zip) is a conserved transcription factor family unique to plants that regulate multiple developmental processes including lignificaion. Stone cell content is a key determinant negatively affecting pear fruit quality, which causes a grainy texture of fruit flesh, because of the lignified cell walls.

RESULTS

In this study, a comprehensive bioinformatics analysis of HD-Zip genes in Chinese white pear (Pyrus bretschneideri) (PbHBs) was performed. Genome-wide identification of the PbHB gene family revealed 67 genes encoding PbHB proteins, which could be divided into four subgroups (I, II, III, and IV). For some members, similar intron/exon structural patterns support close evolutionary relationships within the same subgroup. The functions of each subgroup of the PbHB family were predicted through comparative analysis with the HB genes in Arabidopsis and other plants. Cis-element analysis indicated that PbHB genes might be involved in plant hormone signalling and external environmental responses, such as light, stress, and temperature. Furthermore, RNA-sequencing data and quantitative real-time PCR (RT-qPCR) verification revealed the regulatory roles of PbHB genes in pear stone cell formation. Further, co-expression network analysis revealed that the eight PbHB genes could be classified into different clusters of co-expression with lignin-related genes. Besides, the biological function of PbHB24 in promoting stone cell formation has been demonstrated by overexpression in fruitlets.

CONCLUSIONS

This study provided the comprehensive analysis of PbHBs and highlighted the importance of PbHB24 during stone cell development in pear fruits.

Genome-wide identification, bioinformatics and expression analysis of HD-Zip gene family in peach

BACKGROUND

HD-Zips (Homeodomain-Leucine Zippers) are a class of plant-specific transcription factors that play multiple roles in plant growth and development. Although some functions of HD-Zip transcription factor have been reported in several plants, it has not been comprehensively studied in peach, especially during adventitious root formation of peach cuttings.

RESULTS

In this study, 23 HD-Zip genes distributed on 6 chromosomes were identified from the peach (Prunus persica) genome, and named PpHDZ01-23 according to their positions on the chromosomes. These 23 PpHDZ transcription factors all contained a homeomorphism box domain and a leucine zipper domain, were divided into 4 subfamilies(I-IV) according to the evolutionary analysis, and their promoters contained many different cis-acting elements. Spatio-temporal expression pattern showed that these genes were expressed in many tissues with different levels, and they had distinct expression pattern during adventitious root formation and development.

CONCLUSION

Our results showed the roles of PpHDZs on root formation, which is helpful to better understand the classification and function of peach HD-Zip genes.

Identification of GA20ox2 as a target of ATHB2 and TCP13 during shade response

The shade avoidance syndrome (SAS) is a collective adaptive response of plants under shade highlighted by characteristic phenotypes such as hypocotyl elongation, which is largely mediated by concerted actions of auxin and GA. We identified ATHB2, a homeodomain-leucine zipper (HD-Zip) domain transcription factor known to be rapidly induced under shade condition, as a positive regulator of GA biosynthesis necessary for the SAS by transactivating the expression of GA20ox2, a key gene in the GA biosynthesis pathway. Based on promoter deletion analysis, EMSA and ChIP assay, ATHB2 appears to regulate the GA20ox2 expression as a direct binding target. We also found that the GA20ox2 expression is under negative control by TCP13, the effect of which can be suppressed by presence of ATHB2. Considering a rapid induction kinetics of ATHB2, this relationship between ATHB2 and TCP13 may allow ATHB2 to play a shade-specific activator for GA20ox by derepressing a pre-existing activity of TCP13.

Genome-wide identification and characterization of the HD-Zip gene family and expression analysis in response to stress in Rehmannia glutinosa Libosch

The HD-Zip family of transcription factors is unique to the plant kingdom, and play roles in modulation of plant growth and response to environmental stresses. R. glutinosa is an important Chinese medicinal material. Its yield and quality are susceptible to various stresses. The HD-Zip transcription factors is unique to the plant, and roles in modulation of plant growth and response to environmental stresses. However, there is no relevant research on the HD-ZIP of R. glutinosa. In this study, 92 HD-Zip transcription factors were identified in R. glutinosa, and denominated as RgHDZ1-RgHDZ92. Members of RgHDZ were classified into four groups (HD-ZipI-IV) based on the phylogenetic relationship of Arabidopsis HD-Zip proteins, and each group contains 38, 18, 17, and 19 members, respectively. Expression analyses of RgHDZ genes based on transcriptome data showed that the expression of these genes could be induced by the endophytic fungus of R. glutinosa. Additionally, we showed that RgHDZ genes were differentially expressed in response to drought, waterlogging, temperature, and salinity treatments. This study provides important information for different expression patterns of stress-responsive HD-Zip and may contribute to the better understanding of the different responses of plants to biotic and abiotic stresses, and provide a molecular basis for the cultivation of resistant varieties of R. glutinosa.

Coordination of biradial-to-radial symmetry and tissue polarity by HD-ZIP II proteins

Symmetry establishment is a critical process in the development of multicellular organs and requires careful coordination of polarity axes while cells actively divide within tissues. Formation of the apical style in the Arabidopsis gynoecium involves a bilateral-to-radial symmetry transition, a stepwise process underpinned by the dynamic distribution of the plant morphogen auxin. Here we show that SPATULA (SPT) and the HECATE (HEC) bHLH proteins mediate the final step in the style radialisation process and synergistically control the expression of adaxial-identity genes, HOMEOBOX ARABIDOPSIS THALIANA 3 (HAT3) and ARABIDOPSIS THALIANA HOMEOBOX 4 (ATHB4). HAT3/ATHB4 module drives radialisation of the apical style by promoting basal-to-apical auxin flow and via a negative feedback mechanism that finetune auxin distribution through repression of SPT expression and cytokinin sensitivity. Thus, this work reveals the molecular basis of axes-coordination and hormonal cross-talk during the sequential steps of symmetry transition in the Arabidopsis style.

The apical style in Arabidopsis is formed following a bilateral-to-radial symmetry transition in the gynoecium. Here the authors show that the final step in style radialization is coordinated by the adaxial regulators HAT3 and ATHB4, which are induced by the SPT and HEC transcription factors.

Genome-wide characterization and expression profiling of Eucalyptus grandis HD-Zip gene family in response to salt and temperature stress

BACKGROUND

The HD-Zip transcription factors are unique to plants and play an essential role in plant growth, development and stress responses. The HD-Zip transcription factor family consists of a highly conserved homeodomain (HD) and a leucine zipper domain (LZ) domain. Although the HD-Zip gene family has been extensively studied in many plant species, a systematic study of the Eucalyptus HD-Zip family has not been reported until today. Here, we systematically identified 40 HD-Zip genes in Eucalyptus (Eucalyptus grandis). Besides, we comprehensively analyzed the HD-Zips of Eucalyptus by studying the homology, conserved protein regions, gene structure, 3D structure of the protein, location of the genes on the chromosomes and the expression level of the genes in different tissues.

RESULTS

The HD-Zip family in Eucalyptus has four subfamilies, which is consistent with other plants such as Arabidopsis and rice. Moreover, genes that are in the same group tend to have similar exon-intron structures, motifs, and protein structures. Under salt stress and temperature stress, the Eucalyptus HD-Zip transcription factors show a differential expression pattern.

CONCLUSIONS

Our findings reveal the response of HD-Zip transcription factors under salt and temperature stresses, laying a foundation for future analysis of Eucalyptus HD-Zip transcription factors.

Genetic Contribution of Paleopolyploidy to Adaptive Evolution in Angiosperms

Ancient whole-genome duplications (WGDs or polyploidy) are prevalent in plants, and some WGDs occurred during the timing of severe global environmental changes. It has been suggested that WGDs may have contributed to plant adaptation. However, this still lacks empirical evidence at the genetic level to support the hypothesis. Here, we investigated the survivors of gene duplicates from multiple ancient WGD events on the major branches of angiosperm phylogeny, and aimed to explore genetic evidence supporting the significance of polyploidy. Duplicated genes co-retained from three waves of independent WGDs (∼120 million years ago [Ma], ∼66, and <20 Ma) were investigated in 25 selected species. Gene families functioning in low temperature and darkness were commonly retained gene duplicates after the eight independently occurring WGDs in many lineages around the Cretaceous-Paleocene boundary, when the global cooling and darkness were the two main stresses. Moreover, the commonly retained duplicates could be key factors which may have contributed to the robustness of the critical stress-related pathways. In addition, genome-wide transcription factors (TFs) functioning in stresses tend to retain duplicates after waves of WGDs, and the coselected gene duplicates in many lineages may play critical roles during severe environmental stresses. Collectively, these results shed new light on the significant contribution of paleopolyploidy to plant adaptation during global environmental changes in the evolutionary history of angiosperms.

Overexpressing the HD-Zip class II transcription factor EcHB1 from Eucalyptus camaldulensis increased the leaf photosynthesis and drought tolerance of Eucalyptus

Alteration in the leaf mesophyll anatomy by genetic modification is potentially a promising tool for improving the physiological functions of trees by improving leaf photosynthesis. Homeodomain leucine zipper (HD-Zip) transcription factors are candidates for anatomical alterations of leaves through modification of cell multiplication, differentiation, and expansion. Full-length cDNA encoding a Eucalyptus camaldulensis HD-Zip class II transcription factor (EcHB1) was over-expressed in vivo in the hybrid Eucalyptus GUT5 generated from Eucalyptus grandis and Eucalyptus urophylla. Overexpression of EcHB1 induced significant modification in the mesophyll anatomy of Eucalyptus with enhancements in the number of cells and chloroplasts on a leaf-area basis. The leaf-area-based photosynthesis of Eucalyptus was improved in the EcHB1-overexpression lines, which was due to both enhanced CO₂ diffusion into chloroplasts and increased photosynthetic biochemical functions through increased number of chloroplasts per unit leaf area. Additionally, overexpression of EcHB1 suppressed defoliation and thus improved the growth of Eucalyptus trees under drought stress, which was a result of reduced water loss from trees due to the reduction in leaf area with no changes in stomatal morphology. These results gave us new insights into the role of the HD-Zip II gene.

Overexpressing the HD-Zip class II transcription factor EcHB1 from Eucalyptus camaldulensis increased the leaf photosynthesis and drought tolerance of Eucalyptus

Alteration in the leaf mesophyll anatomy by genetic modification is potentially a promising tool for improving the physiological functions of trees by improving leaf photosynthesis. Homeodomain leucine zipper (HD-Zip) transcription factors are candidates for anatomical alterations of leaves through modification of cell multiplication, differentiation, and expansion. Full-length cDNA encoding a Eucalyptus camaldulensis HD-Zip class II transcription factor (EcHB1) was over-expressed in vivo in the hybrid Eucalyptus GUT5 generated from Eucalyptus grandis and Eucalyptus urophylla. Overexpression of EcHB1 induced significant modification in the mesophyll anatomy of Eucalyptus with enhancements in the number of cells and chloroplasts on a leaf-area basis. The leaf-area-based photosynthesis of Eucalyptus was improved in the EcHB1-overexpression lines, which was due to both enhanced CO₂ diffusion into chloroplasts and increased photosynthetic biochemical functions through increased number of chloroplasts per unit leaf area. Additionally, overexpression of EcHB1 suppressed defoliation and thus improved the growth of Eucalyptus trees under drought stress, which was a result of reduced water loss from trees due to the reduction in leaf area with no changes in stomatal morphology. These results gave us new insights into the role of the HD-Zip II gene.

AtHB2, a class II HD-ZIP protein, negatively regulates the expression of CsANS, which encodes a key enzyme in Camellia sinensis catechin biosynthesis

Tea (Camellia sinensis) is an important cash crop that is beneficial to human health because of its remarkable content of catechins. The biosynthesis of catechins follows the flavonoid pathway, which is highly branched. Among the enzymes involved in catechin biosynthesis, ANTHOCYANIDIN SYNTHASE (CsANS) functions at a branch point and play a critical role. Our previous work has showed that the gene encoding CsANS is regulated by light signals; however, the molecular mechanism behind remains unclear. Here, we cloned a full-length CsANS promoter and found that it contained a cis-element recognized by Arabidopsis thaliana HOMEOBOX2 (AtHB2). AtHB2 constitutes one of the class II HOMEODOMAIN-LEUCINE ZIPPER (HD-ZIP) proteins, which accumulate in the dark and mediate the shade avoidance response in most angiosperms. To analyze the transcription of CsANS in vivo, β-glucuronidase and luciferase reporter genes driven by the obtained promoter were introduced into A. thaliana and Nicotiana attenuata, respectively. In both expression systems there were indications that the A. thaliana PRODUCTION OF ANTHOCYANIN PIGMENT1 (AtPAP1), a MYB transcription factor of flavonoid biosynthesis, increased the activity of the CsANS promoter, while AtHB2 could significantly undermine the effect of AtPAP1. Yeast two-hybrid and bimolecular fluorescence complementation assays showed that AtHB2 interacted with the A. thaliana TRANSPARENT TESTA GLABRA 1 (AtTTG1). A yeast three-hybrid assay further suggested that AtHB2 represses the expression of CsANS and regulates its response to light signals through competitive interactions with AtTTG1. These results show that HD-ZIP II proteins participate in light regulation of flavonoid biosynthesis.

Genome-wide characterization and expression profiling of the relation of the HD-Zip gene family to abiotic stress in barley (Hordeum vulgare L.)

The homeodomain-leucine zipper (HD-Zip) gene family plays an important role in plant growth and environmental responses. At present, research on the HD-Zip gene family of barley is incomplete. In this study, 32 HD-Zip genes (HvHD-Zip 1-32) were identified from the barley genome and were subsequently divided into four subfamilies according to conserved structure and motif analysis. Whole genome replication events in barley and Arabidopsis, rice, and wheat HD-Zip gene families were analyzed, yielding 3, 14 and 25 gene pairs, respectively, but no segmental or tandem duplication events were identified in the barley HD-Zip gene family. Subsequently, quantitative real-time PCR (qRT-PCR) analysis revealed that the HvHD-Zip gene is sensitive to drought stress and that members of the HD-Zip I and HD-Zip IV subfamilies are generally more sensitive to abiotic stresses. Our results suggest a relationship between barley resistance and the potential key HvHD-Zip gene, which lay the foundation for further functional studies.

Multiple Links between HD-Zip Proteins and Hormone Networks

HD-Zip proteins are unique to plants, and contain a homeodomain closely linked to a leucine zipper motif, which are involved in dimerization and DNA binding. Based on homology in the HD-Zip domain, gene structure and the presence of additional motifs, HD-Zips are divided into four families, HD-Zip I–IV. Phylogenetic analysis of HD-Zip genes using transcriptomic and genomic datasets from a wide range of plant species indicate that the HD-Zip protein class was already present in green algae. Later, HD-Zips experienced multiple duplication events that promoted neo- and sub-functionalizations. HD-Zip proteins are known to control key developmental and environmental responses, and a growing body of evidence indicates a strict link between members of the HD-Zip II and III families and the auxin machineries. Interactions of HD-Zip proteins with other hormones such as brassinolide and cytokinin have also been described. More recent data indicate that members of different HD-Zip families are directly involved in the regulation of abscisic acid (ABA) homeostasis and signaling. Considering the fundamental role of specific HD-Zip proteins in the control of key developmental pathways and in the cross-talk between auxin and cytokinin, a relevant role of these factors in adjusting plant growth and development to changing environment is emerging.

Arabidopsis HD-Zip II proteins regulate the exit from proliferation during leaf development in canopy shade

The shade avoidance response is mainly evident as increased plant elongation at the expense of leaf and root expansion. Despite the advances in understanding the mechanisms underlying shade-induced hypocotyl elongation, little is known about the responses to simulated shade in organs other than the hypocotyl. In Arabidopsis, there is evidence that shade rapidly and transiently reduces the frequency of cell division in young first and second leaf primordia through a non-cell-autonomous mechanism. However, the effects of canopy shade on leaf development are likely to be complex and need to be further investigated. Using combined methods of genetics, cell biology, and molecular biology, we uncovered an effect of prolonged canopy shade on leaf development. We show that persistent shade determines early exit from proliferation in the first and second leaves of Arabidopsis. Furthermore, we demonstrate that the early exit from proliferation in the first and second leaves under simulated shade depends at least in part on the action of the Homeodomain-leucine zipper II (HD-Zip II) transcription factors ARABIDOPSIS THALIANA HOMEOBOX2 (ATHB2) and ATHB4. Finally, we provide evidence that the ATHB2 and ATHB4 proteins work in concert. Together the data contribute new insights on the mechanisms controlling leaf development under canopy shade.

Multiple Pathways in the Control of the Shade Avoidance Response

To detect the presence of neighboring vegetation, shade-avoiding plants have evolved the ability to perceive and integrate multiple signals. Among them, changes in light quality and quantity are central to elicit and regulate the shade avoidance response. Here, we describe recent progresses in the comprehension of the signaling mechanisms underlying the shade avoidance response, focusing on Arabidopsis, because most of our knowledge derives from studies conducted on this model plant. Shade avoidance is an adaptive response that results in phenotypes with a high relative fitness in individual plants growing within dense vegetation. However, it affects the growth, development, and yield of crops, and the design of new strategies aimed at attenuating shade avoidance at defined developmental stages and/or in specific organs in high-density crop plantings is a major challenge for the future. For this reason, in this review, we also report on recent advances in the molecular description of the shade avoidance response in crops, such as maize and tomato, and discuss their similarities and differences with Arabidopsis.

Variation in cucumber (Cucumis sativus L.) fruit size and shape results from multiple components acting pre-anthesis and post-pollination

Morphological, QTL, and gene expression analyses indicate variation in cucumber fruit size and shape results from orientation, timing, and extent of cell division and expansion, and suggest candidate gene factors. Variation in cucumber (Cucumis sativus L.) fruit size and shape is highly quantitative, implicating interplay of multiple components. Recent studies have identified numerous fruit size and shape quantitative trait loci (QTL); however, underlying factors remain to be determined. We examined ovary and fruit development of two sequenced cucumber genotypes with extreme differences in fruit size and shape, Chinese Long '9930' (CL9930), and pickling type 'Gy14'. Differences were observed in several independent factors that can influence size and shape: ovule number, rate and period of cell division in longitudinal and cross section in ovaries and fruit, timing and rate of fruit expansion in length and diameter, and cell shape. Level and timing of expression of select fruit growth stage marker genes and candidate fruit size gene homologs associated with cucumber fruit size and shape QTL were examined from 5-day pre-anthesis to 20-day post-pollination. Our results indicate that variation in fruit size and shape results from differences in cell number and shape in longitudinal and cross section, driven in turn by differences in orientation, timing, and duration of cell division and expansion, both pre- and post-anthesis, and suggest candidate genes contributing to determination of cucumber fruit size and shape.

Role of Transcription Factor HAT1 in Modulating Arabidopsis thaliana Response to Cucumber mosaic virus

Arabidopsis thaliana homeodomain-leucine zipper protein 1 (HAT1) belongs to the homeodomain-leucine zipper (HD-Zip) family class II that plays important roles in plant growth and development as a transcription factor. To elucidate further the role of HD-Zip II transcription factors in plant defense, the A. thaliana hat1, hat1hat3 and hat1hat2hat3 mutants and HAT1 overexpression plants (HAT1OX) were challenged with Cucumber mosaic virus (CMV). HAT1OX displayed more susceptibility, while loss-of-function mutants of HAT1 exhibited less susceptibility to CMV infection. HAT1 and its close homologs HAT2 and HAT3 function redundantly, as the triple mutant hat1hat2hat3 displayed increased virus resistance compared with the hat1 and hat1hat3 mutants. Furthermore, the induction of the antioxidant system (the activities and expression of enzymatic antioxidants) and the expression of defense-associated genes were down-regulated in HAT1OX but up-regulated in hat1hat2hat3 when compared with Col-0 after CMV infection. Further evidence showed that the involvement of HAT1 in the anti-CMV defense response might be dependent on salicylic acid (SA) but not jasmonic acid (JA). The SA level or expression of SA synthesis-related genes was decreased in HAT1OX but increased in hat1hat2hat3 compared with Col-0 after CMV infection, but there were little difference in JA level or JA synthesis-related gene expression among HAT1OX or defective plants. In addition, HAT1 expression is dependent on SA accumulation. Taken together, our study indicated that HAT1 negatively regulates plant defense responses to CMV.

Genome Wide Association Mapping in Arabidopsis thaliana Identifies Novel Genes Involved in Linking Allyl Glucosinolate to Altered Biomass and Defense

A key limitation in modern biology is the ability to rapidly identify genes underlying newly identified complex phenotypes. Genome wide association studies (GWAS) have become an increasingly important approach for dissecting natural variation by associating phenotypes with genotypes at a genome wide level. Recent work is showing that the Arabidopsis thaliana defense metabolite, allyl glucosinolate (GSL), may provide direct feedback regulation, linking defense metabolism outputs to the growth, and defense responses of the plant. However, there is still a need to identify genes that underlie this process. To start developing a deeper understanding of the mechanism(s) that modulate the ability of exogenous allyl GSL to alter growth and defense, we measured changes in plant biomass and defense metabolites in a collection of natural 96 A. thaliana accessions fed with 50 μM of allyl GSL. Exogenous allyl GSL was introduced exclusively to the roots and the compound transported to the leaf leading to a wide range of heritable effects upon plant biomass and endogenous GSL accumulation. Using natural variation we conducted GWAS to identify a number of new genes which potentially control allyl responses in various plant processes. This is one of the first instances in which this approach has been successfully utilized to begin dissecting a novel phenotype to the underlying molecular/polygenic basis.

Interplay of HD-Zip II and III transcription factors in auxin-regulated plant development

The homeodomain-leucine zipper (HD-Zip) class of transcription factors is unique to plants. HD-Zip proteins bind to DNA exclusively as dimers recognizing dyad symmetric sequences and act as positive or negative regulators of gene expression. On the basis of sequence homology in the HD-Zip DNA-binding domain, HD-Zip proteins have been grouped into four families (HD-Zip I-IV). Each HD-Zip family can be further divided into subfamilies containing paralogous genes that have arisen through genome duplication. Remarkably, all the members of the HD-Zip IIγ and -δ clades are regulated by light quality changes that induce in the majority of the angiosperms the shade-avoidance response, a process regulated at multiple levels by auxin. Intriguingly, it has recently emerged that, apart from their function in shade avoidance, the HD-Zip IIγ and -δ transcription factors control several auxin-regulated developmental processes, including apical embryo patterning, lateral organ polarity, and gynoecium development, in a white-light environment. This review presents recent advances in our understanding of HD-Zip II protein function in plant development, with particular emphasis on the impact of loss-of-function HD-Zip II mutations on auxin distribution and response. The review also describes evidence demonstrating that HD-Zip IIγ and -δ genes are directly and positively regulated by HD-Zip III transcription factors, primary determinants of apical shoot development, known to control the expression of several auxin biosynthesis, transport, and response genes. Finally, the interplay between HD-Zip II and III transcription factors in embryo apical patterning and organ polarity is discussed.

Transcription factors involved in brassinosteroid repressed gene expression and their regulation by BIN2 kinase

Brassinosteroids (BRs) play important roles in plant growth, development and responses to environmental cues. BRs signal through plasma membrane-localized receptors to control BES1/BZR1 family transcription factors, which mediate the expression of thousands of genes. BRs activate and repress approximately equal numbers of genes. BES1/BZR1 interact with other transcription factors, histone-modifying enzymes, and transcription elongation factors to activate BR-induced genes. However, the mechanisms by which BES1/BZR1 mediate the BR-repressed gene expression are not well understood. Recent studies revealed that 2 BR-repressed transcription factors, MYBL2 and HAT1, cooperate with BES1 to downregulate BR-repressed genes expression. Moreover, BIN2 kinase, a well-established negative regulator in the BR signaling pathway, phosphorylates MYBL2 and HAT1. While BIN2 phosphorylates and destabilizes BES1/BZR1, BIN2 phosphorylated MYBL2 and HAT1 appear to be stabilized. These results not only extended our understanding of BR-repressed gene expression, but also revealed multiple inputs of BR signaling into BR transcriptional networks.