Characterization of homeodomain-leucine zipper genes in the fern Ceratopteris richardii and the evolution of the homeodomain-leucine zipper gene family in vascular plants

Evolutionary Consequences of Functional and Regulatory Divergence of HD-Zip I Transcription Factors as a Source of Diversity in Protein Interaction Networks in Plants

The HD superfamily has been studied in detail for several decades. The plant-specific HD-Zip I subfamily attracts the most attention because of its involvement in plant development and stress responses. In this review, we provide a comprehensive insight into the evolutionary events responsible for the functional redundancy and diversification of the HD-Zip I genes in regulating various biological processes. We summarized the evolutionary history of the HD-Zip family, highlighting the important role of WGDs in its expansion and divergence of retained duplicates in the genome. To determine the relationship between the evolutionary origin and functional conservation of HD-Zip I in different species, we performed a phylogenetic analysis, compared their expression profiles in different tissues and under stress and traced the role of orthologs and paralogs in regulating developmental processes. We found that HD-Zip I from different species have similar gene structures with a highly conserved HD and Zip, bind to the same DNA sequences and are involved in similar biological processes. However, they exhibit a functional diversity, which is manifested in altered expression patterns. Some of them are involved in the regulation of species-specific leaf morphology and phenotypes. Here, we discuss the role of changes in functional domains involved in DNA binding and protein interaction of HD-Zip I and in cis-regulated regions of its target genes in promoting adaptive innovations through the formation of de novo regulatory systems. Understanding the role of the HD-Zip I subfamily in organism-environment interactions remains a challenge for evolutionary developmental biology (evo-devo).

Genome-wide analysis of the homeodomain-leucine zipper family in Lotus japonicus and the overexpression of LjHDZ7 in Arabidopsis for salt tolerance

The homeodomain-leucine zipper (HD-Zip) family participates in plant growth, development, and stress responses. Here, 40 HD-Zip transcription factors of Lotus japonicus were identified and gave an overview of the phylogeny and gene structures. The expression pattern of these candidate genes was determined in different organs and their response to abiotic stresses, including cold, heat, polyethylene glycol and salinity. The expression of the LjHDZ7 was strongly induced by abiotic stress, especially salt stress. Subsequently, LjHDZ7 gene was overexpressed in Arabidopsis. The transgenic plants grew obviously better than Col-0 plants under salt stress. Furthermore, LjHDZ7 transgenic lines accumulated higher proline contents and showed lower electrolyte leakage and MDA contents than Col-0 plants under salt stress. Antioxidant activities of the LjHDZ7 overexpression lines leaf were significantly higher than those of the Col-0 plants under salt stress. The concentration of Na⁺ ion in LjHDZ7 overexpression lines was significantly lower than that of Col-0 in leaf and root parts. The concentration of K⁺ ion in LjHDZ7 overexpression lines was significantly higher than that of Col-0 in the leaf parts. Therefore, these results showed that overexpression of LjHDZ7 increased resistance to salt stress in transgenic Arabidopsis plants, and certain genes of this family can be used as valuable tools for improving abiotic stresses.

The Homeodomain-Leucine Zipper Genes Family Regulates the Jinggangmycin Mediated Immune Response of Oryza sativa to Nilaparvata lugens, and Laodelphax striatellus

The homeodomain-leucine zipper (HDZIP) is an important transcription factor family, instrumental not only in growth but in finetuning plant responses to environmental adversaries. Despite the plethora of literature available, the role of HDZIP genes under chewing and sucking insects remains elusive. Herein, we identified 40 OsHDZIP genes from the rice genome database. The evolutionary relationship, gene structure, conserved motifs, and chemical properties highlight the key aspects of OsHDZIP genes in rice. The OsHDZIP family is divided into a further four subfamilies (i.e., HDZIP I, HDZIP II, HDZIP III, and HDZIP IV). Moreover, the protein–protein interaction and Gene Ontology (GO) analysis showed that OsHDZIP genes regulate plant growth and response to various environmental stimuli. Various microRNA (miRNA) families targeted HDZIP III subfamily genes. The microarray data analysis showed that OsHDZIP was expressed in almost all tested tissues. Additionally, the differential expression patterns of the OsHDZIP genes were found under salinity stress and hormonal treatments, whereas under brown planthopper (BPH), striped stem borer (SSB), and rice leaf folder (RLF), only OsHDZIP3, OsHDZIP4, OsHDZIP40, OsHDZIP10, and OsHDZIP20 displayed expression. The qRT-PCR analysis further validated the expression of OsHDZIP20, OsHDZIP40, and OsHDZIP10 under BPH, small brown planthopper (SBPH) infestations, and jinggangmycin (JGM) spraying applications. Our results provide detailed knowledge of the OsHDZIP gene family resistance in rice plants and will facilitate the development of stress-resilient cultivars, particularly against chewing and sucking insect pests.

Genome-Wide Association Studies for Sex Determination and Cross-Compatibility in Water Yam (Dioscorea alata L.)

Yam (Dioscorea spp.) species are predominantly dioecious, with male and female flowers borne on separate individuals. Cross-pollination is, therefore, essential for gene flow among and within yam species to achieve breeding objectives. Understanding genetic mechanisms underlying sex determination and cross-compatibility is crucial for planning a successful hybridization program. This study used the genome-wide association study (GWAS) approach for identifying genomic regions linked to sex and cross-compatibility in water yam (Dioscorea alata L.). We identified 54 markers linked to flower sex determination, among which 53 markers were on chromosome 6 and one on chromosome 11. Our result ascertained that D. alata is characterized by the male heterogametic sex determination system (XX/XY). The cross-compatibility indices, average crossability rate (ACR) and percentage high crossability (PHC), were controlled by loci on chromosomes 1, 6 and 17. Of the significant loci, SNPs located on chromosomes 1 and 17 were the most promising for ACR and PHC, respectively, and should be validated for use in D. alata hybridization activities to predict cross-compatibility success. A total of 61 putative gene/protein families with direct or indirect influence on plant reproduction were annotated in chromosomic regions controlling the target traits. This study provides valuable insights into the genetic control of D. alata sexual reproduction. It opens an avenue for developing genomic tools for predicting hybridization success in water yam breeding programs.

Genome-Wide Characterization and Expression Analysis of the HD-ZIP Gene Family in Response to Salt Stress in Pepper

HD-ZIP is a unique type of transcription factor in plants, which are closely linked to the regulation of plant growth and development, the response to abiotic stress, and disease resistance. However, there is little known about the HD-ZIP gene family of pepper. In this study, 40 HD-ZIP family members were analyzed in the pepper genome. The analysis indicated that the introns number of Ca-HD-ZIP varied from 1 to 17; the number of amino acids was between 119 and 841; the theoretical isoelectric point was between 4.54 and 9.85; the molecular weight was between 14.04 and 92.56; most of them were unstable proteins. The phylogenetic tree divided CaHD-ZIP into 4 subfamilies; 40 CaHD-ZIP genes were located on different chromosomes, and all of them contained the motif 1; two pairs of CaHD-ZIP parallel genes of six paralogism genes were fragment duplications which occurred in 58.28~88.24 million years ago. There were multiple pressure-related action elements upstream of the start codon of the HD-Z-IP family. Protein interaction network proved to be coexpression phenomenon between ATML1 (CaH-DZ22, CaHDZ32) and At4g048909 (CaHDZ12, CaHDZ31), and three regions of them were highly homology. The expression level of CaHD-ZIP gene was different with tissues and developmental stages, which suggested that CaHD-ZIP may be involved in biological functions during pepper progress. In addition, Pepper HD-ZIP I and II genes played a major role in salt stress. CaHDZ03, CaHDZ 10, CaHDZ17, CaHDZ25, CaHDZ34, and CaHDZ35 were significantly induced in response to salt stress. Notably, the expression of CaHDZ07, CaHDZ17, CaHDZ26, and CaHDZ30, homologs of Arabidopsis AtHB12 and AtHB7 genes, was significantly upregulated by salt stresses. CaHDZ03 possesses two closely linked ABA action elements, and its expression level increased significantly at 4 h under salt stress. qRT-P-CR and transcription analysis showed that the expression of CaHDZ03 and CaHDZ10 was upregulated under short-term salt stress, but CaHDZ10 was downregulated with long-term salt stress, which provided a theoretical basis for research the function of Ca-HDZIP in response to abiotic stress.

A moso bamboo transcription factor, Phehdz1, positively regulates the drought stress response of transgenic rice

78 HD-Zip family genes in Phyllostachys edulis were analyzed. Overexpression of Phehdz1 can improve the drought tolerance of transgenic rice and affect its secondary metabolism. Many studies suggested homeodomain-leucine zipper (HD-Zip) transcription factors are important regulators of plant growth and development, signal transduction, and responses to environmental stresses. In this study, 78 moso bamboo (Phyllostachys edulis) HD-Zip genes were investigated and classified into four subfamilies (HD-Zip I-IV). Additionally, Phehdz1 (HD-Zip I gene) was isolated and confirmed to be highly expressed in the roots. A quantitative real-time PCR analysis indicated Phehdz1 expression was significantly induced by drought, high salinity, and abscisic acid (ABA). A transient expression assay proved that Phehdz1 was localized in the nucleus of tobacco cells. Moreover, it could bind to the core region encoded by the H-box sequence (CAATAATTG) in yeast. In response to mannitol treatments, the Phehdz1-overexpressing transgenic rice had a higher germination rate and longer shoots than the wild-type controls. Moreover, Phehdz1-overexpressing rice plants had a higher survival rate as well as higher relative water and proline contents, but a lower malondialdehyde content, than the WT plants after a 30% polyethylene glycol 6000 treatment. Accordingly, the overexpression of Phehdz1 enhances the drought tolerance of transgenic rice. Many of the differentially expressed genes identified by a transcriptome analysis are involved in MAPK signal transduction and the biosynthesis of secondary metabolites. Thus, the overexpression of Phehdz1 enhances the drought stress tolerance of transgenic rice, while also potentially modulating the expression of metabolism-related genes.

HvHOX9, a novel homeobox leucine zipper transcription factor, positively regulates aluminum tolerance in Tibetan wild barley

Aluminum (Al) toxicity is the primary limiting factor of crop production on acid soils. Tibetan wild barley germplasm is a valuable source of potential genes for breeding barley with acid and Al tolerance. We performed microRNA and RNA sequencing using wild (XZ16, Al-tolerant; XZ61, Al-sensitive) and cultivated (Dayton, Al-tolerant) barley. A novel homeobox-leucine zipper transcription factor, HvHOX9, was identified as a target gene of miR166b and functionally characterized. HvHOX9 was up-regulated by Al stress in XZ16 (but unchanged in XZ61 and Dayton) and was significantly induced only in root tip. Phylogenetic analysis showed that HvHOX9 is most closely related to wheat TaHOX9 and orthologues of HvHOX9 are present in the closest algal relatives of Zygnematophyceae. Barley stripe mosaic virus-induced gene silencing of HvHOX9 in XZ16 led to significantly increased Al sensitivity but did not affect its sensitivity to other metals and low pH. Disruption of HvHOX9 did not change Al concentration in the root cell sap, but led to more Al accumulation in root cell wall after Al exposure. Silencing of HvHOX9 decreased H+ influx after Al exposure. Our findings suggest that miR166b/HvHOX9 play a critical role in Al tolerance by decreasing root cell wall Al binding and increasing apoplastic pH for Al detoxification in the root.

An ontogenetic framework for functional studies in the model fern Ceratopteris richardii

As the sister group to seed plants, ferns are a phylogenetically informative lineage. Functional studies in representatives of the fern lineage are helping bridge the knowledge gap in developmental mechanisms between angiosperms and non-vascular plants. The fern life cycle has the advantage of combining a sizable free-living haploid gametophyte, more amenable for developmental studies than the reduced seed plant gametophyte, with an indeterminate and complex diploid sporophyte. Ceratopteris richardii has long been proposed as a model fern and has recently become tractable due to stable transgenesis and increasing genomic resources, allowing researchers to test explicit questions about gene function in a fern for the first time. As with any model system, a detailed understanding of wild-type morphology and a staged ontogeny are indispensable for the characterization of mutant phenotypes resulting from genetic manipulations. Therefore, the goal of this study is to provide a unified reference ontogeny for this emerging model fern as a tool for comparative evolutionary and developmental studies. It complements earlier research by filling gaps in major stages of development of the haploid gametophyte and diploid sporophyte generations, and provides additional descriptions of the shoot apical meristem and early leaf development. This resource is meant to facilitate not only studies of candidate genes within C. richardii, but also broader ontogenetic comparisons to other model plants.

Genome-Wide Investigation and Expression Profiling of HD-Zip Transcription Factors in Foxtail Millet (Setaria italica L.)

HD-Zip proteins represent the major transcription factors in higher plants, playing essential roles in plant development and stress responses. Foxtail millet is a crop to investigate the systems biology of millet and biofuel grasses and the HD-Zip gene family has not been studied in foxtail millet. For further investigation of the expression profile of the HD-Zip gene family in foxtail millet, a comprehensive genome-wide expression analysis was conducted in this study. We found 47 protein-encoding genes in foxtail millet using BLAST search tools; the putative proteins were classified into four subfamilies, namely, subfamilies I, II, III, and IV. Gene structure and motif analysis indicate that the genes in one subfamily were conserved. Promotor analysis showed that HD-Zip gene was involved in abiotic stress. Duplication analysis revealed that 8 (~17%) hdz genes were tandemly duplicated and 28 (58%) were segmentally duplicated; purifying duplication plays important roles in gene expansion. Microsynteny analysis revealed the maximum relationship in foxtail millet-sorghum and foxtail millet-rice. Expression profiling upon the abiotic stresses of drought and high salinity and the biotic stress of ABA revealed that some genes regulated responses to drought and salinity stresses via an ABA-dependent process, especially sihdz29 and sihdz45. Our study provides new insight into evolutionary and functional analyses of HD-Zip genes involved in environmental stress responses in foxtail millet.

Genome-wide analysis of the HD-ZIP IV transcription factor family in Gossypium arboreum and GaHDG11 involved in osmotic tolerance in transgenic Arabidopsis

HD-ZIP IV proteins belong to the homeodomain-leucine zipper (HD-ZIP) transcription factor family and are involved in trichome development and drought stress in plants. Although some functions of the HD-ZIP IV group are well understood in Arabidopsis, little is known about their function in cotton. In this study, HD-ZIP genes were identified from three Gossypium species (G. arboreum, G. raimondii and G. hirsutum) and clustered into four families (HD-ZIP I, II, III and IV) to separate HD-ZIP IV from the other three families. Systematic analyses of phylogeny, gene structure, conserved domains, and expression profiles in different plant tissues and the expression patterns under osmotic stress in leaves were further conducted in G. arboreum. More importantly, ectopic overexpression of GaHDG11, a representative of the HD-ZIP IV family, confers enhanced osmotic tolerance in transgenic Arabidopsis plants, possibly due to elongated primary root length, lower water loss rates, high osmoprotectant proline levels, significant levels of antioxidants CAT, and/or SOD enzyme activity with reduced levels of MDA. Taken together, these observations may lay the foundation for future functional analysis of cotton HD-ZIP IV genes to unravel their biological roles in cotton.

A New Glabrous Gene (csgl3) Identified in Trichome Development in Cucumber (Cucumis sativus L.)

Spines or trichomes on the fruit of cucumbers enhance their commercial value in China. In addition, glabrous mutants exhibit resistance to aphids and therefore their use by growers can reduce pesticide residues. Previous studies have reported two glabrous mutant plants containing the genes, csgl1 and csgl2. In the present study, a new glabrous mutant, NCG157, was identified showing a gene interaction effect with csgl1 and csgl2. This mutant showed the glabrous character on stems, leaves, tendrils, receptacles and ovaries, and there were no spines or tumors on the fruit surface. Inheritance analysis showed that a single recessive gene, named csgl3, determined the glabrous trait. An F₂ population derived from the cross of two inbred lines 9930 (a fresh market type from Northern China that exhibits trichomes) and NCG157 (an American processing type with glabrous surfaces) was used for genetic mapping of the csgl3 gene. By combining bulked segregant analysis (BAS) with molecular markers, 18 markers, including two simple sequence repeats (SSR), nine insertion deletions (InDel) and seven derived cleaved amplified polymorphism sequences (dCAPs), were identified to link to the csgl3 gene. All of the linked markers were used as anchor loci to locate the csgl3 gene on cucumber chromosome 6. The csgl3 gene was mapped between the dCAPs markers dCAPs-21 and dCAPs-19, at genetic distances of 0.05 cM and 0.15 cM, respectively. The physical distance of this region was 19.6 kb. Three markers, InDel-19, dCAPs-2 and dCAPs-11, co-segregated with csgl3. There were two candidate genes in the region, Csa6M514860 and Csa6M514870. Quantitative real-time PCR showed that the expression of Csa6M514870 was higher in the tissues of 9930 than that of NCG157, and this was consistent with their phenotypic characters. Csa6M514870 is therefore postulated to be the candidate gene for the development of trichomes in cucumber. This study will facilitate marker-assisted selection (MAS) of the smooth plant trait in cucumber breeding and provide for future cloning of csgl3.

Ferns: the missing link in shoot evolution and development

Shoot development in land plants is a remarkably complex process that gives rise to an extreme diversity of forms. Our current understanding of shoot developmental mechanisms comes almost entirely from studies of angiosperms (flowering plants), the most recently diverged plant lineage. Shoot development in angiosperms is based around a layered multicellular apical meristem that produces lateral organs and/or secondary meristems from populations of founder cells at its periphery. In contrast, non-seed plant shoots develop from either single apical initials or from a small population of morphologically distinct apical cells. Although developmental and molecular information is becoming available for non-flowering plants, such as the model moss Physcomitrella patens, making valid comparisons between highly divergent lineages is extremely challenging. As sister group to the seed plants, the monilophytes (ferns and relatives) represent an excellent phylogenetic midpoint of comparison for unlocking the evolution of shoot developmental mechanisms, and recent technical advances have finally made transgenic analysis possible in the emerging model fern Ceratopteris richardii. This review compares and contrasts our current understanding of shoot development in different land plant lineages with the aim of highlighting the potential role that the fern C. richardii could play in shedding light on the evolution of underlying genetic regulatory mechanisms.

High-efficiency stable transformation of the model fern species Ceratopteris richardii via microparticle bombardment

Ferns represent the most closely related extant lineage to seed plants. The aquatic fern Ceratopteris richardii has been subject to research for a considerable period of time, but analyses of the genetic programs underpinning developmental processes have been hampered by a large genome size, a lack of available mutants, and an inability to create stable transgenic lines. In this paper, we report a protocol for efficient stable genetic transformation of C. richardii and a closely related species Ceratopteris thalictroides using microparticle bombardment. Indeterminate callus was generated and maintained from the sporophytes of both species using cytokinin treatment. In proof-of-principle experiments, a 35S::β-glucuronidase (GUS) expression cassette was introduced into callus cells via tungsten microparticles, and stable transformants were selected via a linked hygromycin B resistance marker. The presence of the transgene in regenerated plants and in subsequent generations was validated using DNA-blot analysis, reverse transcription-polymerase chain reaction, and GUS staining. GUS staining patterns in most vegetative tissues corresponded with constitutive gene expression. The protocol described in this paper yields transformation efficiencies far greater than those previously published and represents a significant step toward the establishment of a tractable fern genetic model.

Molecular evolution of the HD-ZIP I gene family in legume genomes

Homeodomain leucine zipper I (HD-ZIP I) genes were used to increase the plasticity of plants by mediating external signals and regulating growth in response to environmental conditions. The way genomic histories drove the evolution of the HD-ZIP I family in legume species was described; HD-ZIP I genes were searched in Lotus japonicus, Medicago truncatula, Cajanus cajan and Phaseolus vulgaris, and then divided into five clades through phylogenetic analysis. Microsynteny analysis was made based on genomic segments containing the HD-ZIP I genes. Some pairs turned out to conform with syntenic genome regions, while others corresponded to those that were inverted, expanded, or contracted after the divergence of legumes. Besides, we dated their duplications by Ks analysis and demonstrated that all the blocks were formed after the monocot-dicot split; we observed Ka/Ks ratios representing strong purifying selections in the four legume species which might have been followed by gene loss and rearrangement.

Comprehensive analysis of the homeodomain-leucine zipper IV transcription factor family in Cucumis sativus

The class IV homeodomain-leucine zipper (HD-Zip IV) proteins are plant-specific transcriptional factors known to play crucial roles in plant growth and development. In this study, 11 cucumber (Cucumis sativus) HD-Zip IV genes were identified in the version 2 cucumber genome and found to be distributed unevenly across the chromosomes. The CsHDZIV (Cucumis sativus homeodomain-leucine zipper IV) gene family is smaller than in other studied species (except for rice) because of the absence of gene duplication events. Phylogenetic analysis showed that HD-Zip IV genes from cucumber, Arabidopsis, tomato, cotton, maize, and rice could be classified into five subgroups. All CsHDZIV genes appear to be derived from a basic module containing 11 exons in the coding region. Two conserved motifs of 21 and 19 nucleotides were found in the 3'-untranslated regions of six CsHDZIV genes, suggesting that post-transcriptional regulation may play a role in regulation of CsHDZIV genes. In addition, 6 of 11 CsHDZIV genes were found to undergo alternative splicing events. Reverse transcription PCR analysis showed that all CsHDZIV genes (except one) were expressed and showed preferential expression in reproductive organs.

Function of the HD-Zip I gene Oshox22 in ABA-mediated drought and salt tolerances in rice

Oshox22 belongs to the homeodomain-leucine zipper (HD-Zip) family I of transcription factors, most of which have unknown functions. Here we show that the expression of Oshox22 is strongly induced by salt stress, abscisic acid (ABA), and polyethylene glycol treatment (PEG), and weakly by cold stress. Trans-activation assays in yeast and transient expression analyses in rice protoplasts demonstrated that Oshox22 is able to bind the CAAT(G/C)ATTG element and acts as a transcriptional activator that requires both the HD and Zip domains. Rice plants homozygous for a T-DNA insertion in the promoter region of Oshox22 showed reduced Oshox22 expression and ABA content, decreased sensitivity to ABA, and enhanced tolerance to drought and salt stresses at the seedling stage. In contrast, transgenic rice over-expressing Oshox22 showed increased sensitivity to ABA, increased ABA content, and decreased drought and salt tolerances. Based on these results, we conclude that Oshox22 affects ABA biosynthesis and regulates drought and salt responses through ABA-mediated signal transduction pathways.

Genome-wide identification, evolutionary expansion, and expression profile of homeodomain-leucine zipper gene family in poplar (Populus trichocarpa)

Background

Homeodomain-leucine zipper (HD-ZIP) proteins are plant-specific transcriptional factors known to play crucial roles in plant development. Although sequence phylogeny analysis of Populus HD-ZIPs was carried out in a previous study, no systematic analysis incorporating genome organization, gene structure, and expression compendium has been conducted in model tree species Populus thus far.

Principal Findings

In this study, a comprehensive analysis of Populus HD-ZIP gene family was performed. Sixty-three full-length HD-ZIP genes were found in Populus genome. These Populus HD-ZIP genes were phylogenetically clustered into four distinct subfamilies (HD-ZIP I–IV) and predominately distributed across 17 linkage groups (LG). Fifty genes from 25 Populus paralogous pairs were located in the duplicated blocks of Populus genome and then preferentially retained during the sequential evolutionary courses. Genomic organization analyses indicated that purifying selection has played a pivotal role in the retention and maintenance of Populus HD-ZIP gene family. Microarray analysis has shown that 21 Populus paralogous pairs have been differentially expressed across different tissues and under various stresses, with five paralogous pairs showing nearly identical expression patterns, 13 paralogous pairs being partially redundant and three paralogous pairs diversifying significantly. Quantitative real-time RT-PCR (qRT-PCR) analysis performed on 16 selected Populus HD-ZIP genes in different tissues and under both drought and salinity stresses confirms their tissue-specific and stress-inducible expression patterns.

Conclusions

Genomic organizations indicated that segmental duplications contributed significantly to the expansion of Populus HD-ZIP gene family. Exon/intron organization and conserved motif composition of Populus HD-ZIPs are highly conservative in the same subfamily, suggesting the members in the same subfamilies may also have conservative functionalities. Microarray and qRT-PCR analyses showed that 89% (56 out of 63) of Populus HD-ZIPs were duplicate genes that might have been retained by substantial subfunctionalization. Taken together, these observations may lay the foundation for future functional analysis of Populus HD-ZIP genes to unravel their biological roles.

The evolution and functional significance of leaf shape in the angiosperms

Angiosperm leaves manifest a remarkable diversity of shapes that range from developmental sequences within a shoot and within crown response to microenvironment to variation among species within and between communities and among orders or families. It is generally assumed that because photosynthetic leaves are critical to plant growth and survival, variation in their shape reflects natural selection operating on function. Several non-mutually exclusive theories have been proposed to explain leaf shape diversity. These include: thermoregulation of leaves especially in arid and hot environments, hydraulic constraints, patterns of leaf expansion in deciduous species, biomechanical constraints, adaptations to avoid herbivory, adaptations to optimise light interception and even that leaf shape variation is a response to selection on flower form. However, the relative importance, or likelihood, of each of these factors is unclear. Here we review the evolutionary context of leaf shape diversification, discuss the proximal mechanisms that generate the diversity in extant systems, and consider the evidence for each the above hypotheses in the context of the functional significance of leaf shape. The synthesis of these broad ranging areas helps to identify points of conceptual convergence for ongoing discussion and integrated directions for future research.

Gene family structure, expression and functional analysis of HD-Zip III genes in angiosperm and gymnosperm forest trees

BACKGROUND

Class III Homeodomain Leucine Zipper (HD-Zip III) proteins have been implicated in the regulation of cambium identity, as well as primary and secondary vascular differentiation and patterning in herbaceous plants. They have been proposed to regulate wood formation but relatively little evidence is available to validate such a role. We characterised and compared HD-Zip III gene family in an angiosperm tree, Populus spp. (poplar), and the gymnosperm Picea glauca (white spruce), representing two highly evolutionarily divergent groups.

RESULTS

Full-length cDNA sequences were isolated from poplar and white spruce. Phylogenetic reconstruction indicated that some of the gymnosperm sequences were derived from lineages that diverged earlier than angiosperm sequences, and seem to have been lost in angiosperm lineages. Transcript accumulation profiles were assessed by RT-qPCR on tissue panels from both species and in poplar trees in response to an inhibitor of polar auxin transport. The overall transcript profiles HD-Zip III complexes in white spruce and poplar exhibited substantial differences, reflecting their evolutionary history. Furthermore, two poplar sequences homologous to HD-Zip III genes involved in xylem development in Arabidopsis and Zinnia were over-expressed in poplar plants. PtaHB1 over-expression produced noticeable effects on petiole and primary shoot fibre development, suggesting that PtaHB1 is involved in primary xylem development. We also obtained evidence indicating that expression of PtaHB1 affected the transcriptome by altering the accumulation of 48 distinct transcripts, many of which are predicted to be involved in growth and cell wall synthesis. Most of them were down-regulated, as was the case for several of the poplar HD-Zip III sequences. No visible physiological effect of over-expression was observed on PtaHB7 transgenic trees, suggesting that PtaHB1 and PtaHB7 likely have distinct roles in tree development, which is in agreement with the functions that have been assigned to close homologs in herbaceous plants.

CONCLUSIONS

This study provides an overview of HD-zip III genes related to woody plant development and identifies sequences putatively involved in secondary vascular growth in angiosperms and in gymnosperms. These gene sequences are candidate regulators of wood formation and could be a source of molecular markers for tree breeding related to wood properties.

A comprehensive classification and evolutionary analysis of plant homeobox genes

The full complement of homeobox transcription factor sequences, including genes and pseudogenes, was determined from the analysis of 10 complete genomes from flowering plants, moss, Selaginella, unicellular green algae, and red algae. Our exhaustive genome-wide searches resulted in the discovery in each class of a greater number of homeobox genes than previously reported. All homeobox genes can be unambiguously classified by sequence evolutionary analysis into 14 distinct classes also characterized by conserved intron–exon structure and by unique codomain architectures. We identified many new genes belonging to previously defined classes (HD-ZIP I to IV, BEL, KNOX, PLINC, WOX). Other newly identified genes allowed us to characterize PHD, DDT, NDX, and LD genes as members of four new evolutionary classes and to define two additional classes, which we named SAWADEE and PINTOX. Our comprehensive analysis allowed us to identify several newly characterized conserved motifs, including novel zinc finger motifs in SAWADEE and DDT. Members of the BEL and KNOX classes were found in Chlorobionta (green plants) and in Rhodophyta. We found representatives of the DDT, WOX, and PINTOX classes only in green plants, including unicellular green algae, moss, and vascular plants. All 14 homeobox gene classes were represented in flowering plants, Selaginella, and moss, suggesting that they had already differentiated in the last common ancestor of moss and vascular plants.

The evolutionary development of plant body plans

Evolutionary developmental biology, cladistic analyses, and paleontological insights make it increasingly clear that regulatory mechanisms operating during embryogenesis and early maturation tend to be highly conserved over great evolutionary time scales, which can account for the conservative nature of the body plans in the major plant and animal clades. At issue is whether morphological convergences in body plans among evolutionarily divergent lineages are the result of adaptive convergence or 'genome recall' and 'process orthology'. The body plans of multicellular photosynthetic eukaryotes ('plants') are reviewed, some of their important developmental/physiological regulatory mechanisms discussed, and the evidence that some of these mechanisms are phyletically ancient examined. We conclude that endosymbiotic lateral gene transfers, gene duplication and functional divergence, and the co-option of ancient gene networks were key to the evolutionary divergence of plant lineages.

Genome-wide identification, classification, evolutionary expansion and expression analyses of homeobox genes in rice

Homeobox genes play a critical role in regulating various aspects of plant growth and development. In the present study, we identified a total of 107 homeobox genes in the rice genome and grouped them into ten distinct subfamilies based upon their domain composition and phylogenetic analysis. A significantly large number of homeobox genes are located in the duplicated segments of the rice genome, which suggests that the expansion of homeobox gene family, in large part, might have occurred due to segmental duplications in rice. Furthermore, microarray analysis was performed to elucidate the expression profiles of these genes in different tissues and during various stages of vegetative and reproductive development. Several genes with predominant expression during various stages of panicle and seed development were identified. At least 37 homeobox genes were found to be differentially expressed significantly (more than two-fold; P < 0.05) under various abiotic stress conditions. The results of the study suggest a critical role of homeobox genes in reproductive development and abiotic stress signaling in rice, and will facilitate the selection of candidate genes of agronomic importance for functional validation.

A genome-wide survey of HD-Zip genes in rice and analysis of drought-responsive family members

The homeodomain leucine zipper (HD-Zip) genes encode transcription factors that have diverse functions in plant development and have often been implicated in stress adaptation. The HD-Zip genes are the most abundant group of homeobox (HB) genes in plants and do not occur in other eukaryotes. This paper describes the complete annotation of the HD-Zip families I, II and III from rice and compares these gene families with Arabidopsis in a phylogeny reconstruction. Orthologous pairs of rice and Arabidopsis HD-Zip genes were predicted based on neighbour joining and maximum parsimony (MP) trees with support of conserved intron-exon organization. Additionally, a number of HD-Zip genes appeared to be unique to rice. Searching of EST and cDNA databases and expression analysis using RT-PCR showed that 30 out of 31 predicted rice HD-Zip genes are expressed. Most HD-Zip genes were broadly expressed in mature plants and seedlings, but others showed more organ specific patterns. Like in Arabidopsis and other dicots, a subset of the rice HD-Zip I and II genes was found to be regulated by drought stress. We identified both drought-induced and drought-repressed HD-Zip genes and demonstrate that these genes are differentially regulated in drought-sensitive versus drought-tolerant rice cultivars. The drought-repressed HD-Zip family I gene, Oshox4, was selected for promoter-GUS analysis, showing that drought-responsiveness of Oshox4 is controlled by the promoter and that Oshox4 expression is predominantly vascular-specific. Loss-of-function analysis of Oshox4 revealed no specific phenotype, but overexpression analysis suggested a role for Oshox4 in elongation and maturation processes.

Evolution of the class III HD-Zip gene family in land plants

Arabidopsis class III homeodomain-leucine zipper (HD-Zip III) proteins play overlapping, distinct, and antagonistic roles in key aspects of development that have evolved during land plant evolution. To better understand this gene family's role in plant evolution and development as well as to address broader questions of how duplicated genes functionally diversify, we investigated the evolutionary history of this gene family. Phylogenetic analyses including homologs from diverse land plants indicate that a gene duplication event before the angiosperm--gymnosperm split gave rise to two gene lineages that diversified during angiosperm plant radiation. Heterologous expression of an HD-Zip III gene from the nonvascular plant moss within the Arabidopsis HD-zip III revoluta mutant modified but did not complement the phenotype. Comparison of the expression domains of flowering and nonflowering plant homologs indicate an ancestral role in vascular development and organ initiation but not in specifying organ polarity, a prominent role for angiosperm homologs.

A feedback regulatory module formed by LITTLE ZIPPER and HD-ZIPIII genes

The Arabidopsis thaliana REVOLUTA (REV) protein is a member of the class III homeodomain-leucine zipper (HD-ZIPIII) proteins. REV is a potent regulator of leaf polarity and vascular development. Here, we report the identification of a gene family that encodes small leucine zipper-containing proteins (LITTLE ZIPPER [ZPR] proteins) where the leucine zipper is similar to that found in REV, PHABULOSA, and PHAVOLUTA proteins. The transcript levels of the ZPR genes increase in response to activation of a steroid-inducible REV protein. We show that the ZPR proteins interact with REV in vitro and that ZPR3 prevents DNA binding by REV in vitro. Overexpression of ZPR proteins in Arabidopsis results in phenotypes similar to those seen when HD-ZIPIII function is reduced. We propose a negative feedback model in which REV promotes transcription of the ZPR genes. The ZPR proteins in turn form heterodimers with the REV protein, preventing it from binding DNA. The HD-ZIPIII/ZPR regulatory module would serve not only to dampen the effect of fluctuations in HD-ZIPIII protein levels but more importantly would provide a potential point of regulation (control over the ratio of inactive heterodimers to active homodimers) that could be influenced by other components of the pathway governing leaf polarity.

Distinct developmental mechanisms reflect the independent origins of leaves in vascular plants

Vascular plants diverged more than 400 million years ago into two lineages, the lycophytes and the euphyllophytes . Leaf-like organs evolved independently in these two groups . Microphylls in lycophytes are hypothesized to have originated as lateral outgrowths of tissue that later became vascularized (the enation theory) or through the sterilization of sporangia (the sterilization hypothesis) . Megaphylls in euphyllophytes are thought to represent modified lateral branches . The fossil record also indicates that the seed plant megaphyll evolved uniquely in the ancestor of seed plants, independent of megaphylls in ferns, because seed plants evolved from leafless progymnosperm ancestors . Surprisingly, a recent study of KNOX and ARP gene expression in a lycophyte was reported to indicate recruitment of a similar mechanism for determinacy in both types of leaves . We examined the expression of Class III HD-Zip genes in the lycophyte Selaginella kraussiana and in two gymnosperms, Ginkgo and Pseudotsuga. Our data indicate that mechanisms promoting leaf initiation, vascularization, and polarity are quite different in lycophytes and seed plants, consistent with the hypotheses that megaphylls originated as lateral branches whereas microphylls originated as tissue outgrowths.

Evolution of class III homeodomain-leucine zipper genes in streptophytes

Land plants underwent tremendous evolutionary change following the divergence of the ancestral lineage from algal relatives. Several important developmental innovations appeared as the embryophyte clade diversified, leading to the appearance of new organs and tissue types. To understand how these changes came about, we need to identify the fundamental genetic developmental programs that are responsible for growth, patterning, and differentiation and describe how these programs were modified and elaborated through time to produce novel morphologies. Class III homeodomain-leucine zipper (class III HD-Zip) genes, identified in the model plant Arabidopsis thaliana, provide good candidates for basic land plant patterning genes. We show that these genes may have evolved in a common ancestor of land plants and their algal sister group and that the gene family has diversified as land plant lineages have diversified. Phylogenetic analysis, expression data from nonflowering lineages, and evidence from Arabidopsis and other flowering plants indicate that class III HD-Zip genes acquired new functions in sporophyte apical growth, vascular patterning and differentiation, and leaf development. Modification of expression patterns that accompanied diversification of class III HD-Zip genes likely played an important role in the evolution of land plant form.

MEKHLA, a novel domain with similarity to PAS domains, is fused to plant homeodomain-leucine zipper III proteins

Homeodomain (HD) proteins play important roles in the development of plants, fungi, and animals. Here we identify a novel domain, MEKHLA, in the C terminus of HD-Leu zipper (HD-ZIP) III plant HD proteins that shares similarity with a group of bacterial proteins and a protein from the green alga Chlamydomonas reinhardtii. The group of bacterial MEKHLA proteins is found in cyanobacteria and other bacteria often found associated with plants. Phylogenetic analysis suggests that a MEKHLA protein transferred, possibly from a cyanobacterium or an early chloroplast, into the nuclear genome of an early plant in a first step, and attached itself to the C terminus of an HD-ZIP IV homeobox gene in a second step. Further position-specific iterated-BLAST searches with the bacterial MEKHLA proteins revealed a subregion within the MEKHLA domain that shares significant similarity with the PAS domain. The PAS domain is a sensory module found in many proteins through all kingdoms of life. It is involved in light, oxygen, and redox potential sensation. The fact that HD-ZIP III proteins are transcription factors that have this sensory domain attached to their C terminus uncovers a potential new signaling pathway in plants.

Homeodomain leucine zipper class I genes in Arabidopsis. Expression patterns and phylogenetic relationships

Members of the homeodomain leucine zipper (HDZip) family of transcription factors are present in a wide range of plants, from mosses to higher plants, but not in other eukaryotes. The HDZip genes act in developmental processes, including vascular tissue and trichome development, and several of them have been suggested to be involved in the mediation of external signals to regulate plant growth. The Arabidopsis (Arabidopsis thaliana) genome contains 47 HDZip genes, which, based on sequence criteria, have been grouped into four different classes: HDZip I to IV. In this article, we present an overview of the class I HDZip genes in Arabidopsis. We describe their expression patterns, transcriptional regulation properties, duplication history, and phylogeny. The phylogeny of HDZip class I genes is supported by data on the duplication history of the genes, as well as the intron/exon patterning of the HDZip-encoding motifs. The HDZip class I genes were found to be widely expressed and partly to have overlapping expression patterns at the organ level. Further, abscisic acid or water deficit treatments and different light conditions affected the transcript levels of a majority of the HDZip I genes. Within the gene family, our data show examples of closely related HDZip genes with similarities in the function of the gene product, but a divergence in expression pattern. In addition, six HDZip class I proteins tested were found to be activators of gene expression. In conclusion, several HDZip I genes appear to regulate similar cellular processes, although in different organs or tissues and in response to different environmental signals.

KNOX homeobox genes potentially have similar function in both diploid unicellular and multicellular meristems, but not in haploid meristems

Members of the class 1 knotted-like homeobox (KNOX) gene family are important regulators of shoot apical meristem development in angiosperms. To determine whether they function similarly in seedless plants, three KNOX genes (two class 1 genes and one class 2 gene) from the fern Ceratopteris richardii were characterized. Expression of both class 1 genes was detected in the shoot apical cell, leaf primordia, marginal part of the leaves, and vascular bundles by in situ hybridization, a pattern that closely resembles that of class 1 KNOX genes in angiosperms with compound leaves. The fern class 2 gene was expressed in all sporophyte tissues examined, which is characteristic of class 2 gene expression in angiosperms. All three CRKNOX genes were not detected in gametophyte tissues by RNA gel blot analysis. Arabidopsis plants overexpressing the fern class 1 genes resembled plants that overexpress seed plant class 1 KNOX genes in leaf morphology. Ectopic expression of the class 2 gene in Arabidopsis did not result in any unusual phenotypes. Taken together with phylogenetic analysis, our results suggest that (a) the class 1 and 2 KNOX genes diverged prior to the divergence of fern and seed plant lineages, (b) the class 1 KNOX genes function similarly in seed plant and fern sporophyte meristem development despite their differences in structure, (c) KNOX gene expression is not required for the development of the fern gametophyte, and (d) the sporophyte and gametophyte meristems of ferns are not regulated by the same developmental mechanisms at the molecular level.

Athb-12, a homeobox-leucine zipper domain protein from Arabidopsis thaliana, increases salt tolerance in yeast by regulating sodium exclusion

An Arabidopsis cDNA clone that encodes Athb-12, a homeobox-leucine zipper domain protein (HD-Zip), was isolated by functional complementation of the NaCl-sensitive phenotype of a calcineurin (CaN)-deficient yeast mutant (cnbDelta, regulatory subunit null). CaN, a Ca2+/calmodulin-dependent protein phosphatase, regulates Na+ ion homeostasis in yeast. Expression of Athb-12 increased NaCl tolerance but not osmotic stress tolerance of these cnbDelta cells. Furthermore, expression of two other HD-Zip from Arabidopsis, Athb-1 and -7, did not suppress NaCl sensitivity of cnbDelta cells. These results suggest that Athb-12 specifically functions in Na+ ion homeostasis in yeast. Consistent with these observations, expression of Athb-12 in yeast turned on transcription of the NaCl stress-inducible PMR2A, which encodes a Na+/Li+ translocating P-type ATPase, and decreased Na+ levels in yeast cells. To investigate the biological function of Athb-12 in Arabidopsis, we performed Northern blot analysis. Expression of Athb-12 was dramatically induced by NaCl and ABA treatments, but not by KCl. In vivo targeting experiments using a green fluorescent protein reporter indicated that Athb-12 was localized to the nucleus. These results suggest that Athb-12 is a putative transcription factor that may be involved in NaCl stress responses in plants.

Isolation and characterization of two homeodomain leucine zipper genes from the dioecious plant Silene latifolia

Homeodomain leucine zipper (HD-Zip) genes encode transcription factors that are characterized by both a homeodomain and a leucine zipper motif. Two HD-Zip genes were isolated from cDNA of the male flower bud of the dioecious plant Silene latifolia. The two isolated genes, SlHDL1 and SlHDL2, encode proteins with the characteristics of HD-Zip transcription factors belonging to HD-Zip classes I and II, respectively. The expression patterns of SlHDL1 and SlHDL2 throughout the floral developmental stages were studied using real-time PCR and in situ hybridization. SlHDL1 is specifically expressed in the outermost layer of the anthers and gynoeciums with a patchy pattern in the inner layers, suggesting that the product of SlHDL1 plays a role in the early developmental stage of the epidermal tissues of these floral organs. Its expression pattern in the anthers and gynoeciums suggests an involvement in differentiation of the reproductive organs. On the other hand, real-time PCR revealed accumulation of SlHDL2 transcripts in the anther and pollen grains of the male flower. These results suggest that SlHDL1 and SlHDL2 regulate specific targets in restricted regions leading to floral organ differentiation in S. latifolia.

OSTF1: a HD-GL2 family homeobox gene is developmentally regulated during early embryogenesis in rice

In many eukaryotic organisms, homeobox genes are important regulators that specify the cell fate and body plan in early embryogenesis. In this study, a gene designated OSTF1 (Oryza sativa transcription factor 1) encoding a homeodomain protein in rice was isolated and characterized. The encoded OSTF1, although sharing only approximately 51% sequence identity with other HD-GL2 members, contains four characteristic motifs (an N-terminal acidic region, a homeodomain, a truncated leucine zipper, and a START domain). OSTF1 was detected as a single copy gene in rice. The transcripts were absent in young panicle or mature spikelet before anthesis, but appeared very early in the pollinated grain with a transient profile. In vegetative tissues examined, expression was only detectable in root. In situ hybridization analysis on developing grains revealed that OSTF1 was strongly and uniformly expressed in the embryo at the globular stage and preferentially localized to the protoderm at 3-6 d after pollination. Expression was also detectable in the integument and throughout the endosperm. Although OSTF1 is not closely related to the remaining HD-GL2 members in sequences, this gene exhibits an analogous epidermis-preferential expression pattern.

Primary phloem-specific expression of a Zinnia elegans homeobox gene

Some plant homeobox genes are expressed specifically in vascular cells and are assumed to function in the differentiation of specific types of vascular cells. However, homeobox genes exhibiting primary phloem-specific expression have not been reported. To elucidate the molecular mechanisms of vascular development, we undertook to isolate from Zinnia elegans primary phloem-specific homeobox genes that may function in phloem development. An HD-Zip type homeobox gene, ZeHB3, was isolated. This gene encodes a class I HD-Zip protein, and constitutes a gene subfamily with the Daucus carota gene CHB6, and Arabidopsis thaliana genes Athb-5, Athb-6, and Athb-16. In situ hybridization of 1-, 14- and 50-day-old plants demonstrated that ZeHB3 mRNA accumulation is restricted to a few cells destined to differentiate into phloem cells and to the immature phloem cells surrounding the sieve elements and companion cells. ZeHB3 protein was also localized to immature phloem cells. These findings clearly indicate that ZeHB3 is a novel homeobox gene that marks, and may function in, the early stages of phloem differentiation.

PaHB1 is an evolutionary conserved HD-GL2 homeobox gene expressed in the protoderm during Norway spruce embryo development

In angiosperms, the protoderm or outer cell layer is the first tissue to differentiate in the embryo proper. In gymnosperms, it is not known whether a protoderm is defined and similarly differentiated. Here, we report a cDNA designated PaHB1 (for Picea abies Homeobox1), which is expressed during somatic embryogenesis in Norway spruce. PaHB1 exon/intron organization and its corresponding protein are highly similar to those of the HD-GL2 angiosperm counterparts. A phylogenetic analysis reveals that PaHB1 is strongly associated with one subclass consisting of protoderm/epiderm-specific genes. Moreover, PaHB1 expression switches from a ubiquitous expression in proembryogenic masses to an outer cell layer-specific localization during somatic embryo development. Ectopic expression of PaHB1 in somatic embryos leads to an early developmental block. The transformed embryos lack a smooth surface. These findings show that the PaHB1 expression pattern is highly analogous to angiosperm HD-GL2 homologues, suggesting similarities in the definition of the outer cell layer in seed plants.

Isolation of homeodomain-leucine zipper genes from the moss Physcomitrella patens and the evolution of homeodomain-leucine zipper genes in land plants

Homeobox genes encode transcription factors involved in many aspects of developmental processes. The homeodomain-leucine zipper (HD-Zip) genes, which are characterized by the presence of both a homeodomain and a leucine zipper motif, form a clade within the homeobox superfamily and were previously reported only from vascular plants. Here we report the isolation of 10 HD-Zip genes (named PPHB:1-PPHB:10) from the moss Physcomitrella patens. Based on a phylogenetic analysis of the 10 PPHB: genes and previously reported vascular plant HD-Zip genes, all of the PPHB: genes except Pphb3 belong to three of the four HD-Zip subfamilies (HD-Zip I, II, and III), indicating that these subfamilies originated before the divergence of the vascular plant and moss lineages. Pphb3 is sister to the HD-Zip II subfamily and has some distinctive characteristics, including the difference of the a(1) and d(1) sites of its leucine zipper motif, which are well conserved in each HD-Zip subfamily. Comparison of the genetic divergence of representative HD-Zip I and II genes showed that the evolutionary rate of HD-Zip I genes was faster than that of HD-Zip II genes.

Negative autoregulation of the Arabidopsis homeobox gene ATHB-2

The Arabidopsis homeobox gene ATHB-2 is tightly regulated by light signals, and is thought to direct morphological changes during shade avoidance responses. To understand how ATHB-2 mediates light signals in plant morphogenesis, we investigated its transcriptional network. We constructed a gene encoding a chimeric transcription factor (HD-Zip-2-V-G) that is expected to activate target genes of ATHB-2 in a glucocorticoid-dependent manner. In transgenic Arabidopsis plants expressing HD-Zip-2-V-G, glucocorticoid treatment activates the ATHB-2 gene itself, independent of de novo protein synthesis. An in vitro DNase I-footprinting experiment showed that recombinant ATHB-2 protein specifically bound to an ATHB-2 promoter region. These complementary results indicate that ATHB-2 recognizes its own promoter. Consistent with the fact that ATHB-2 itself has been shown to act as a repressor, expression of the endogenous ATHB-2 gene was repressed in transgenic plants overexpressing an ATHB-2 transgene. Moreover, target-gene analysis using the HD-Zip-2-V-G suggested that ATHB-2 recognizes other HD-Zip II subfamily genes. We conclude that ATHB-2 has a negative autoregulatory loop and may be involved in a complicated transcriptional network involving paralogous genes, as is the case with animal homeobox genes.