Role of Homeodomain leucine zipper (HD-Zip) IV transcription factors in plant development and plant protection from deleterious environmental factors

Identification of transcription factor BnHDG4-A08 as a novel candidate associated with the accumulation of oleic, linoleic, linolenic, and erucic acid in Brassica napus

KEY MESSAGE

We screened 47 significantly associated haplotype blocks for oleic, linoleic, linolenic, and erucic acid, with 17 blocks influencing multiple traits. A novel candidate of transcription factor BnHDG4 A08 influencing oleic, linoleic, linolenic, and erucic acid was identified, by a joint strategy of haplotype-based genome-wide association study, genomic resequencing, gene cloning, and co-expression network Fatty acid (FA) composition determines the quality and economic value of rapeseed oil (Brassica napus). However, the molecular network of FAs is unclear. In the current study, multi-strategies of haplotype-based genome-wide association study (GWAS), genomic resequencing, gene cloning, and co-expression network were joint to reveal novel genetic factors influencing FA accumulation in rapeseed. We identified 47 significantly associated haplotype blocks for oleic, linoleic, linolenic, and erucic acid, with 17 blocks influencing multiple traits, using a haplotype-based GWAS with phenotype data from 203 Chinese semi-winter accessions. A total of 61 rapeseed orthologs involved in acyl-lipid metabolism, carbohydrate metabolism, or photosynthesis were identified in these 17 blocks. Among these genes, BnHDG4-A08, encoding a class IV homeodomain leucine-zipper transcription factor, exhibited two single-nucleotide polymorphisms (SNPs) in the exon and intron, with significant associations with oleic, linoleic, linolenic, and erucic acid. Gene cloning further validated two SNPs in the exon of BnHDG4-A08 in a population with 75 accessions, leading to two amino acid changes (T372A and P366L) and significant variation of oleic, linoleic, linolenic, and erucic acid. A competitive allele-specific PCR (KASP) marker based on the SNPs was successfully developed and validated. Moreover, 98 genes exhibiting direct interconnections and high weight values with BnHDG4-A08 were identified through co-expression network analysis using transcriptome data from 13 accessions. Our study identified a novel FA candidate of transcription factor BnHDG4-A08 influencing oleic, linoleic, linolenic, and erucic acid. This gene provides a potential promising gene resource for the novel mechanistic understanding of transcription factors regulating FA accumulation.

ACL1-ROC4/5 complex reveals a common mechanism in rice response to brown planthopper infestation and drought

Brown planthopper (BPH) is the most destructive insect pest of rice. Drought is the most detrimental environmental stress. BPH infestation causes adaxial leaf-rolling and bulliform cells (BCs) shrinkage similar to drought. The BC-related abaxially curled leaf1 (ACL1) gene negatively regulates BPH resistance and drought tolerance, with decreased cuticular wax in the gain-of-function mutant ACL1-D. ACL1 shows an epidermis-specific expression. The TurboID system and multiple biochemical assays reveal that ACL1 interacts with the epidermal-characteristic rice outermost cell-specific (ROC) proteins. ROC4 and ROC5 positively regulate BPH resistance and drought tolerance through modulating cuticular wax and BCs, respectively. Overexpression of ROC4 and ROC5 both rescue ACL1-D mutant in various related phenotypes. ACL1 competes with ROC4/ROC5 in homo-dimer and hetero-dimer formation, and interacts with the repressive TOPLESS-related proteins. Altogether, we illustrate that ACL1-ROC4/5 complexes synergistically mediate drought tolerance and BPH resistance through regulating cuticular wax content and BC development in rice, a mechanism that might facilitate BPH-resistant breeding.

Integrated transcriptomic and pathway analyses of sorghum plants revealed the molecular mechanisms of host defense against aphids

Sugarcane aphid has emerged as a major pest of sorghum recently, and a few sorghum accessions were identified for resistance to this aphid so far. However, the molecular and genetic mechanisms underlying this resistance are still unclear. To understand these mechanisms, transcriptomics was conducted in resistant Tx2783 and susceptible BTx623 sorghum genotypes infested with sugarcane aphids. A principal component analysis revealed differences in the transcriptomic profiles of the two genotypes. The pathway analysis of the differentially expressed genes (DEGs) indicated the upregulation of a set of genes related to signal perception (nucleotide-binding, leucine-rich repeat proteins), signal transduction [mitogen-activated protein kinases signaling, salicylic acid (SA), and jasmonic acid (JA)], and plant defense (transcription factors, flavonoids, and terpenoids). The upregulation of the selected DEGs was verified by real-time quantitative PCR data analysis, performed on the resistant and susceptible genotypes. A phytohormone bioassay experiment showed a decrease in aphid population, plant mortality, and damage in the susceptible genotype when treated with JA and SA. Together, the results indicate that the set of genes, pathways, and defense compounds is involved in host plant resistance to aphids. These findings shed light on the specific role of each DEG, thus advancing our understanding of the genetic and molecular mechanisms of host plant resistance to aphids.

Genome-wide characterization and expression profiling of the HD-ZIP gene family in Acoraceae under salinity and cold stress

The Homeodomain-Leucine Zipper (HD-ZIP) transcription factors play a pivotal role in governing various aspects of plant growth, development, and responses to abiotic stress. Despite the well-established importance of HD-ZIPs in many plants, their functions in Acoraceae, the basal lineage of monocots, remain largely unexplored. Using recently published whole-genome data, we identified 137 putative HD-ZIPs in two Acoraceae species, Acorus gramineus and Acorus calamus. These HD-ZIP genes were further classified into four subfamilies (I, II, III, IV) based on phylogenetic and conserved motif analyses, showcasing notable variations in exon-intron patterns among different subfamilies. Two microRNAs, miR165/166, were found to specifically target HD-ZIP III genes with highly conserved binding sites. Most cis-acting elements identified in the promoter regions of Acoraceae HD-ZIPs are involved in modulating light and phytohormone responsiveness. Furthermore, our study revealed an independent duplication event in Ac. calamus and a one-to-multiple correspondence between HD-ZIP genes of Ac. calamus and Ac. gramineus. Expression profiles obtained from qRT-PCR demonstrated that HD-ZIP I genes are strongly induced by salinity stress, while HD-ZIP II members have contrasting stress responses in two species. HD-ZIP III and IV genes show greater sensitivity in stress-bearing roots. Taken together, these findings contribute valuable insights into the roles of HD-ZIP genes in stress adaptation and plant resilience in basal monocots, illuminating their multifaceted roles in plant growth, development, and response to abiotic stress.

Genome-Wide Identification, Characterization, and Expression Analysis of the HD-Zip Gene Family in Lagerstroemia for Regulating Plant Height

The Homeodomain leucine zipper (HD-Zip) family of transcription factors is crucial in helping plants adapt to environmental changes and promoting their growth and development. Despite research on the HD-Zip family in various plants, studies in Lagerstroemia (crape myrtle) have not been reported. This study aimed to address this gap by comprehensively analyzing the HD-Zip gene family in crape myrtle. This study identified 52 HD-Zip genes in the genome of Lagerstroemia indica, designated as LinHDZ1-LinHDZ52. These genes were distributed across 22 chromosomes and grouped into 4 clusters (HD-Zip I-IV) based on their phylogenetic relationships. Most gene structures and motifs within each cluster were conserved. Analysis of protein properties, gene structure, conserved motifs, and cis-acting regulatory elements revealed diverse roles of LinHDZs in various biological contexts. Examining the expression patterns of these 52 genes in 6 tissues (shoot apical meristem, tender shoot, and mature shoot) of non-dwarf and dwarf crape myrtles revealed that 2 LinHDZs (LinHDZ24 and LinHDZ14) and 2 LinHDZs (LinHDZ9 and LinHDZ35) were respectively upregulated in tender shoot of non-dwarf crape myrtles and tender and mature shoots of dwarf crape myrtles, which suggested the important roles of these genes in regulate the shoot development of Lagerstroemia. In addition, the expression levels of 2 LinHDZs (LinHDZ23 and LinHDZ34) were significantly upregulated in the shoot apical meristem of non-dwarf crape myrtle. These genes were identified as key candidates for regulating Lagerstroemia plant height. This study enhanced the understanding of the functions of HD-Zip family members in the growth and development processes of woody plants and provided a theoretical basis for further studies on the molecular mechanisms underlying Lagerstroemia plant height.

Nuclear localization of HD-Zip IV transcription factor GLABRA2 is driven by Importin α

GLABRA2 (GL2), a class IV homeodomain leucine-zipper (HD-Zip IV) transcription factor (TF) from Arabidopsis , is a developmental regulator of specialized cell types in the epidermis. GL2 contains a putative monopartite nuclear localization sequence (NLS) partially overlapping with its homeodomain (HD). We demonstrate that NLS deletion or alanine substitution of its basic residues (KRKRKK) affects nuclear localization and results in a loss-of-function phenotype. Fusion of the predicted NLS (GTNKRKRKKYHRH) to the fluorescent protein EYFP is sufficient for its nuclear localization in roots and trichomes. The functional NLS is evolutionarily conserved in a distinct subset of HD-Zip IV members including PROTODERMAL FACTOR2 (PDF2). Despite partial overlap of the NLS with the HD, genetic dissection of the NLS from PDF2 indicates that nuclear localization and DNA binding are separable functions. Affinity purification of GL2 from plant tissues followed by mass spectrometry-based proteomics identified Importin α (IMPα) isoforms as potential GL2 interactors. NLS structural prediction and molecular docking studies with IMPα-3 revealed major interacting residues. Split-ubiquitin cytosolic yeast two-hybrid assays suggest interaction between GL2 and four IMPα isoforms from Arabidopsis. Direct interactions were verified in vitro by co-immunoprecipitation with recombinant proteins. IMPα triple mutants ( impα- 1,2,3 ) exhibit defects in EYFP:GL2 nuclear localization in trichomes but not in roots, consistent with tissue-specific and redundant functions of IMPα isoforms in Arabidopsis . Taken together, our findings provide mechanistic evidence for IMPα-dependent nuclear localization of GL2 and other HD-Zip IV TFs in plants.

One sentence summary

GLABRA2, a representative HD-Zip IV transcription factor from Arabidopsis , contains an evolutionarily conserved monopartite nuclear localization sequence that is recognized by Importin α for translocation to the nucleus, a process that is necessary for cell-type differentiation of the epidermis.

HD-Zip IV transcription factors: Drivers of epidermal cell fate integrate metabolic signals

The leaf epidermis comprises the outermost layer of cells that protect plants against environmental stresses such as drought, ultraviolet radiation, and pathogen attack. Research over the past decades highlights the role of class IV homeodomain leucine-zipper (HD-Zip IV) transcription factors (TFs) in driving differentiation of various epidermal cell types, such as trichomes, guard cells, and pavement cells. Evolutionary origins of this family in the charophycean green algae and HD-Zip-specific gene expression in the maternal genome provide clues to unlocking their secrets which include ties to cell cycle regulation. A distinguishing feature of these TFs is the presence of a lipid binding pocket that integrates metabolic information with gene expression. Identities of metabolic partners are beginning to emerge, uncovering feedback loops to maintain epidermal cell specification. Discoveries of associated molecular mechanisms are revealing fascinating links to phospholipid and sphingolipid metabolism and mechanical signaling.

Genome-Wide Identification of Homeodomain Leucine Zipper (HD-ZIP) Transcription Factor, Expression Analysis, and Protein Interaction of HD-ZIP IV in Oil Palm Somatic Embryogenesis

Understanding the molecular mechanisms underlying somatic embryogenesis is essential for resolving the problems related to the long duration of the process and a low rate of somatic embryo induction in oil palm tissue culture. In this study, we conducted genome-wide identification of the oil palm homeodomain leucine zipper (EgHD-ZIP) family, which is one of the plant-specific transcription factors reported to be involved in embryogenesis. EgHD-ZIP proteins can be divided into four subfamilies, which have similarities in gene structure and protein-conserved motifs within a group. In silico expression analysis showed that the expression of EgHD-ZIP gene members in the EgHD-ZIP I and II families, as well as most members in the EgHD-ZIP IV family, were up-regulated during the zygotic and somatic embryo developmental stages. In contrast, the expression of EgHD-ZIP gene members in the EgHD-ZIP III family was down-regulated during zygotic embryo development. Moreover, the expression of EgHD-ZIP IV genes was validated in the oil palm callus and at the somatic embryo stages (globular, torpedo, and cotyledon). The results revealed that EgHD-ZIP IV genes were up-regulated at the late stages of somatic embryogenesis (torpedo and cotyledon). While BABY BOOM (BBM) gene was up-regulated at the early stage of somatic embryogenesis (globular). In addition, the Yeast-two hybrid assay revealed the direct binding between all members of the oil palm HD-ZIP IV subfamily (EgROC2, EgROC3, EgROC5, EgROC8, and EgBBM). Our findings suggested that the EgHD-ZIP IV subfamily and EgBBM work together to regulate somatic embryogenesis in oil palms. This process is important because it is widely used in plant biotechnology to produce large quantities of genetically identical plants, which can be used for oil palm tissue culture improvement.

Recent allopolyploidy alters Spartina microRNA expression in response to xenobiotic-induced stress

Environmental contamination by xenobiotics represents a major threat for natural ecosystems and public health. In response, xenobiotic detoxification is a fundamental trait of organisms for developmental plasticity and stress tolerance, but the underlying molecular mechanisms remain poorly understood in plants. To decipher this process, we explored the consequences of allopolyploidy on xenobiotic tolerance in the genus Spartina Schreb. Specifically, we focused on microRNAs (miRNAs) owing to their central function in the regulation of gene expression patterns, including responses to stress. Small RNA-Seq was conducted on the parents S. alterniflora and S. maritima, their F1 hybrid S. x townsendii and the allopolyploid S. anglica under phenanthrene-induced stress (phe), a model Polycyclic Aromatic Hydrocarbon (PAH) compound. Differentially expressed miRNAs in response to phe were specifically identified within species. In complement, the respective impacts of hybridization and genome doubling were detected, through changes in miRNA expression patterns between S. x townsendii, S. anglica and the parents. The results support the impact of allopolyploidy in miRNA-guided regulation of plant response to phe. In total, we identified 17 phe-responsive miRNAs in Spartina among up-regulated MIR156 and down-regulated MIR159. We also describe novel phe-responsive miRNAs as putative Spartina-specific gene expression regulators in response to stress. Functional validation using Arabidopsis (L.) Heynh. T-DNA lines inserted in homologous MIR genes was performed, and the divergence of phe-responsive miRNA regulatory networks between Arabidopsis and Spartina was discussed.

Genome-wide identification and systematic analysis of the HD-Zip gene family and its roles in response to pH in Panax ginseng Meyer

BACKGROUND

Ginseng, Panax ginseng Meyer, is a traditional herb that is immensely valuable both for human health and medicine and for medicinal plant research. The homeodomain leucine zipper (HD-Zip) gene family is a plant-specific transcription factor gene family indispensable in the regulation of plant growth and development and plant response to environmental stresses.

RESULTS

We identified 117 HD-Zip transcripts from the transcriptome of ginseng cv. Damaya that is widely grown in Jilin, China where approximately 60% of the world's ginseng is produced. These transcripts were positioned to 64 loci in the ginseng genome and the ginseng HD-Zip genes were designated as PgHDZ genes. Identification of 82 and 83 PgHDZ genes from the ginseng acc. IR826 and cv. ChP genomes, respectively, indicated that the PgHDZ gene family consists of approximately 80 PgHDZ genes. Phylogenetic analysis showed that the gene family originated after Angiosperm split from Gymnosperm and before Dicots split from Monocots. The gene family was classified into four subfamilies and has dramatically diverged not only in gene structure and functionality but also in expression characteristics. Nevertheless, co-expression network analysis showed that the activities of the genes in the family remain significantly correlated, suggesting their functional correlation. Five hub PgHDZ genes were identified that might have central functions in ginseng biological processes and four of them were shown to be actively involved in plant response to environmental pH stress in ginseng.

CONCLUSIONS

The PgHDZ gene family was identified from ginseng and analyzed systematically. Five potential hub genes were identified and four of them were shown to be involved in ginseng response to environmental pH stress. The results provide new insights into the characteristics, diversity, evolution, and functionality of the PgHDZ gene family in ginseng and lay a foundation for comprehensive research of the gene family in plants.

EYFP fusions to HD-Zip IV transcription factors enhance their stability and lead to phenotypic changes in Arabidopsis

Green fluorescent protein (GFP) and its derivatives are extensively used for labeling cells, monitoring gene expression and/or tracking the localization or interactions of proteins. Previous reports of detrimental effects of fluorescent protein (FP) expression include cytotoxicity and interference with fusion protein function or localization. Only a few studies have documented the fluorescent tag-specific effects in plants. Here, we show that placing an enhanced yellow FP (EYFP) tag on the amino-terminus of GLABRA2 (GL2) and PROTODERMAL FACTOR2 (PDF2), two developmentally important HD-Zip IV transcription factors from Arabidopsis, enhances their protein stability. Additionally, expression of EYFP:GL2 not only rescued the gl2 null mutant but also resulted in the abnormal development of abaxially curled leaves associated with EYFP-tag induced GL2 overexpression. Our study raises concerns on the use of FPs regarding their effects on the native properties of target proteins as well as biological consequences of fusion protein expression on morphology.

Transcriptomic analysis of bakanae disease resistant and susceptible rice genotypes in response to infection by Fusarium fujikuroi

BACKGROUND

Fusarium fujikuroi causing bakanae is one of the most significant pathogens of rice and much responsible for yield losses thereby emerging as a major risk to food security.

METHODS

In the present study transcriptomic analysis was conducted between two contrasting resistant (C101A51) and susceptible (Rasi) genotypes of rice with the combinations of C101A51 control (CC) vs. C101A51 inoculated (CI); Rasi control (RC) vs. Rasi inoculated (RI) and C101A51 inoculated (CI) vs. Rasi inoculated (RI).

RESULTS

In CC vs. CI commonly expressed genes were 12,764. Out of them 567 (4%) were significantly upregulated and 1399 (9%) genes were downregulated. For the RC vs. RI 14, 333 (79%) genes were commonly expressed. For CI vs. RI 13,662 (72%) genes were commonly expressed. Genes related to cysteine proteinase inhibitor 10, disease resistance protein TAO1-like, oleosin 16 kDa-like, pathogenesis-related protein (PR1), (PR4), BTB/POZ and MATH domain-containing protein 5-like, alpha-amylase isozyme were upregulated in resistant genotype C101A51. Whereas, genes related to GDSL esterase/lipase, serine glyoxylate aminotransferase, CASP-like protein 2C1, WAT1-related protein, Cytoplasmic linker associated proteins, xyloglucan endotransglucosylase/hydrolase protein and β-D xylosidase 7 were upregulated in susceptible genotype Rasi. Gene ontology analysis showed functions related to defence response (GO:0006952), regulation of plant hypersensitive type response (GO:0010363), Potassium ion transmembrane activity (GO:0015079), chloroplast (GO:0009507), response to wounding (GO:0009611), xylan biosynthetic process (GO:0045492) were upregulated in resistant genotype C101A51 under inoculated conditions.

CONCLUSION

Real time PCR based validation of the selected DEGs showed that the qRT-PCR was consistent with the RNA-Seq results. This is the first transcriptomic study against bakanae disease of rice in Indian genotypes. Further, functional studies on identified genes and their utilization through different methodology will be helpful for the development of bakanae disease management strategies.

The roles of HD-ZIP proteins in plant abiotic stress tolerance

Homeodomain leucine zipper (HD-ZIP) proteins are plant-specific transcription factors that contain a homeodomain (HD) and a leucine zipper (LZ) domain. The highly conserved HD binds specifically to DNA and the LZ mediates homodimer or heterodimer formation. HD-ZIP transcription factors control plant growth, development, and responses to abiotic stress by regulating downstream target genes and hormone regulatory pathways. HD-ZIP proteins are divided into four subclasses (I-IV) according to their sequence conservation and function. The genome-wide identification and expression profile analysis of HD-ZIP proteins in model plants such as Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa) have improved our understanding of the functions of the different subclasses. In this review, we mainly summarize and discuss the roles of HD-ZIP proteins in plant response to abiotic stresses such as drought, salinity, low temperature, and harmful metals. HD-ZIP proteins mainly mediate plant stress tolerance by regulating the expression of downstream stress-related genes through abscisic acid (ABA) mediated signaling pathways, and also by regulating plant growth and development. This review provides a basis for understanding the roles of HD-ZIP proteins and potential targets for breeding abiotic stress tolerance in plants.

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.

RNAi silencing of wheat gliadins alters the network of transcription factors that regulate the synthesis of seed storage proteins toward maintaining grain protein levels

Gluten proteins are responsible for the unique viscoelastic properties of wheat dough, but they also trigger the immune response in celiac disease patients. RNA interference (RNAi) wheat lines with strongly silenced gliadins were obtained to reduce the immunogenic response of wheat. The E82 line presents the highest reduction of gluten, but other grain proteins increased, maintaining a total nitrogen content comparable to that of the wild type. To better understand the regulatory mechanisms in response to gliadin silencing, we carried out a transcriptomic analysis of grain and leaf tissues of the E82 line during grain filling. A network of candidate transcription factors (TFs) that regulates the synthesis of the seed storage proteins (SSPs), α-amylase/trypsin inhibitors, lipid transfer proteins, serpins, and starch in the grain was obtained. Moreover, there were a high number of differentially expressed genes in the leaf of E82, where processes such as nutrient availability and transport were enriched. The source-sink communication between leaf and grain showed that many down-regulated genes were related to protease activity, amino acid and sugar metabolism, and their transport. In the leaf, specific proline transporters and lysine-histidine transporters were down- and up-regulated, respectively. Overall, the silencing of gliadins in the RNAi line is compensated mainly with lysine-rich globulins, which are not related to the proposed candidate network of TFs, suggesting that these proteins are regulated independently of the other SSPs. Results reported here can explain the protein compensation mechanisms and contribute to decipher the complex TF network operating during grain filling.

Genome-Wide Characterization and Expression Analysis of HD-ZIP Gene Family in Dendrobium officinale

The homeodomain-leucine zipper (HD-ZIP) gene family, as one of the plant-specific transcription factor families, plays an important role in regulating plant growth and development as well as in response to diverse stresses. Although it has been extensively characterized in many plants, the HD-ZIP family is not well-studied in Dendrobium officinale, a valuable ornamental and traditional Chinese medicinal herb. In this study, 37 HD-ZIP genes were identified in Dendrobium officinale (Dohdzs) through the in silico genome search method, and they were classified into four subfamilies based on phylogenetic analysis. Exon–intron structure and conserved protein domain analyses further supported the prediction with the same group sharing similar gene and protein structures. Furthermore, their expression patterns were investigated in nine various tissues and under cold stress based on RNA-seq datasets to obtain the tissue-specific and cold-responsive candidates. Finally, Dohdz5, Dohdz9, and Dohdz12 were selected to validate their expression through qRT-PCR analysis, and they displayed significantly differential expression under sudden chilling stress, suggesting they might be the key candidates underlying cold stress response. These findings will contribute to better understanding of the regulatory roles of the HD-ZIP family playing in cold stress and also will provide the vital targets for further functional studies of HD-ZIP genes in D. officinale.

Genome-Wide Analysis of the Homeobox Gene Family and Identification of Drought-Responsive Members in Populus trichocarpa

Homeobox (HB) genes play critical roles in the regulation of plant morphogenesis, growth and development. Here, we identified a total of 156 PtrHB genes from the Populus trichocarpa genome. According to the topologies and taxonomy of the phylogenetic tree constructed by Arabidopsis thaliana HB members, all PtrHB proteins were divided into six subgroups, namely HD-ZIP, ZF-HD, HB-PHD, TALE, WOX and HB-OTHERS. Multiple alignments of conserved homeodomains (HDs) revealed the conserved loci of each subgroup, while gene structure analysis showed similar exon–intron gene structures, and motif analysis indicated the similarity of motif number and pattern in the same subgroup. Promoter analysis indicated that the promoters of PtrHB genes contain a series of cis-acting regulatory elements involved in responding to various abiotic stresses, indicating that PtrHBs had potential functions in these processes. Collinearity analysis revealed that there are 96 pairs of 127 PtrHB genes mainly distributing on Chromosomes 1, 2, and 5. We analyzed the spatio-temporal expression patterns of PtrHB genes, and the virus-induced gene silencing (VIGS) of PtrHB3 gene resulted in the compromised tolerance of poplar seedlings to mannitol treatment. The bioinformatics on PtrHB family and preliminary exploration of drought-responsive genes can provide support for further study of the family in woody plants, especially in drought-related biological processes. It also provides a direction for developing new varieties of poplar with drought resistance. Overall, our results provided significant information for further functional analysis of PtrHB genes in poplar and demonstrated that PtrHB3 is a dominant gene regulating tolerance to water stress treatment in poplar seedlings.

Genome-wide analysis of poplar HD-Zip family and over-expression of PsnHDZ63 confers salt tolerance in transgenic Populus simonii × P.nigra

HD-Zip is a plant-specific HB transcription factor, which participates in plant development and stress response. In this study, we identified 63 poplar HD-Zip transcription factors, which were randomly distributed on 19 chromosomes of poplar. Based on the gene structure and phylogenetic relationship, these members are divided into four groups, which have a variety of collinear evolutionary relationships. They also have rich segmental replication events and experienced strong purification selection. Based on RNA-seq analysis, we profiled the expression pattern of the 63 HD-Zip members under salt stress. Subsequently, we carried out in-depth study on the significantly up-regulated PsnHDZ63 in the stems and leaves. The transgenic Populus simonii × P.nigra plants over-expressing PsnHDZ63 displayed better morphological and physiological indexes than WT under salt stress. In addition, PsnHDZ63 enhanced salt stress tolerance of transgenic lines by combining effective stress-resistant elements to improve reactive oxygen species scavenging ability. These studies laid a foundation for a comprehensive understanding of poplar HD-Zip family members, and revealed the important role of PsnHDZ63 in plant salt tolerance.

StAR-Related Lipid Transfer (START) Domains Across the Rice Pangenome Reveal How Ontogeny Recapitulated Selection Pressures During Rice Domestication

The StAR-related lipid transfer (START) domain containing proteins or START proteins, encoded by a plant amplified family of evolutionary conserved genes, play important roles in lipid binding, transport, signaling, and modulation of transcriptional activity in the plant kingdom, but there is limited information on their evolution, duplication, and associated sub- or neo-functionalization. Here we perform a comprehensive investigation of this family across the rice pangenome, using 10 wild and cultivated varieties. Conservation of START domains across all 10 rice genomes suggests low dispensability and critical functional roles for this family, further supported by chromosomal mapping, duplication and domain structure patterns. Analysis of synteny highlights a preponderance of segmental and dispersed duplication among STARTs, while transcriptomic investigation of the main cultivated variety Oryza sativa var. japonica reveals sub-functionalization amongst genes family members in terms of preferential expression across various developmental stages and anatomical parts, such as flowering. Ka/Ks ratios confirmed strong negative/purifying selection on START family evolution, implying that ontogeny recapitulated selection pressures during rice domestication. Our findings provide evidence for high conservation of START genes across rice varieties in numbers, as well as in their stringent regulation of Ka/Ks ratio, and showed strong functional dependency of plants on START proteins for their growth and reproductive development. We believe that our findings advance the limited knowledge about plant START domain diversity and evolution, and pave the way for more detailed assessment of individual structural classes of START proteins among plants and their domain specific substrate preferences, to complement existing studies in animals and yeast.

HD-ZIP Gene Family: Potential Roles in Improving Plant Growth and Regulating Stress-Responsive Mechanisms in Plants

Exploring the molecular foundation of the gene-regulatory systems underlying agronomic parameters or/and plant responses to both abiotic and biotic stresses is crucial for crop improvement. Thus, transcription factors, which alone or in combination directly regulated the targeted gene expression levels, are appropriate players for enlightening agronomic parameters through genetic engineering. In this regard, homeodomain leucine zipper (HD-ZIP) genes family concerned with enlightening plant growth and tolerance to environmental stresses are considered key players for crop improvement. This gene family containing HD and LZ domain belongs to the homeobox superfamily. It is further classified into four subfamilies, namely HD-ZIP I, HD-ZIP II, HD-ZIP III, and HD-ZIP IV. The first HD domain-containing gene was discovered in maize cells almost three decades ago. Since then, with advanced technologies, these genes were functionally characterized for their distinct roles in overall plant growth and development under adverse environmental conditions. This review summarized the different functions of HD-ZIP genes in plant growth and physiological-related activities from germination to fruit development. Additionally, the HD-ZIP genes also respond to various abiotic and biotic environmental stimuli by regulating defense response of plants. This review, therefore, highlighted the various significant aspects of this important gene family based on the recent findings. The practical application of HD-ZIP biomolecules in developing bioengineered plants will not only mitigate the negative effects of environmental stresses but also increase the overall production of crop plants.

Transcriptomic Analysis of Salt-Stress-Responsive Genes in Barley Roots and Leaves

Barley is characterized by a rich genetic diversity, making it an important model for studies of salinity response with great potential for crop improvement. Moreover, salt stress severely affects barley growth and development, leading to substantial yield loss. Leaf and root transcriptomes of a salt-tolerant Tunisian landrace (Boulifa) exposed to 2, 8, and 24 h salt stress were compared with pre-exposure plants to identify candidate genes and pathways underlying barley's response. Expression of 3585 genes was upregulated and 5586 downregulated in leaves, while expression of 13,200 genes was upregulated and 10,575 downregulated in roots. Regulation of gene expression was severely impacted in roots, highlighting the complexity of salt stress response mechanisms in this tissue. Functional analyses in both tissues indicated that response to salt stress is mainly achieved through sensing and signaling pathways, strong transcriptional reprograming, hormone osmolyte and ion homeostasis stabilization, increased reactive oxygen scavenging, and activation of transport and photosynthesis systems. A number of candidate genes involved in hormone and kinase signaling pathways, as well as several transcription factor families and transporters, were identified. This study provides valuable information on early salt-stress-responsive genes in roots and leaves of barley and identifies several important players in salt tolerance.

Plant Transcription Factors Involved in Drought and Associated Stresses

Transcription factors (TFs) play a significant role in signal transduction networks spanning the perception of a stress signal and the expression of corresponding stress-responsive genes. TFs are multi-functional proteins that may simultaneously control numerous pathways during stresses in plants-this makes them powerful tools for the manipulation of regulatory and stress-responsive pathways. In recent years, the structure-function relationships of numerous plant TFs involved in drought and associated stresses have been defined, which prompted devising practical strategies for engineering plants with enhanced stress tolerance. Vast data have emerged on purposely basic leucine zipper (bZIP), WRKY, homeodomain-leucine zipper (HD-Zip), myeloblastoma (MYB), drought-response elements binding proteins/C-repeat binding factor (DREB/CBF), shine (SHN), and wax production-like (WXPL) TFs that reflect the understanding of their 3D structure and how the structure relates to function. Consequently, this information is useful in the tailored design of variant TFs that enhances our understanding of their functional states, such as oligomerization, post-translational modification patterns, protein-protein interactions, and their abilities to recognize downstream target DNA sequences. Here, we report on the progress of TFs based on their interaction pathway participation in stress-responsive networks, and pinpoint strategies and applications for crops and the impact of these strategies for improving plant stress tolerance.

The underground life of homeodomain-leucine zipper transcription factors

Roots are the anchorage organs of plants, responsible for water and nutrient uptake, exhibiting high plasticity. Root architecture is driven by the interactions of biomolecules, including transcription factors and hormones that are crucial players regulating root plasticity. Multiple transcription factor families are involved in root development; some, such as ARFs and LBDs, have been well characterized, whereas others remain less well investigated. In this review, we synthesize the current knowledge about the involvement of the large family of homeodomain-leucine zipper (HD-Zip) transcription factors in root development. This family is divided into four subfamilies (I-IV), mainly according to structural features, such as additional motifs aside from HD-Zip, as well as their size, gene structure, and expression patterns. We explored and analyzed public databases and the scientific literature regarding HD-Zip transcription factors in Arabidopsis and other species. Most members of the four HD-Zip subfamilies are expressed in specific cell types and several individuals from each group have assigned functions in root development. Notably, a high proportion of the studied proteins are part of intricate regulation pathways involved in primary and lateral root growth and development.

The WRKY transcription factor AaGSW2 promotes glandular trichome initiation in Artemisia annua

Glandular secreting trichomes (GSTs) synthesize and secrete large quantities of secondary metabolites, some of which have well-established commercial value. An example is the anti-malarial compound artemisinin, which is synthesized in the GSTs of Artemisia annua. Accordingly, there is considerable interest in understanding the processes that regulate GST density as a strategy to increase artemisinin production. In this study, we identified a GST-specific WRKY transcription factor from A. annua, AaGSW2, which is positively regulated by the direct binding of the homeodomain proteins AaHD1 and AaHD8 to the L1-box of the AaGSW2 promoter. Overexpression of AaGSW2 in A. annua significantly increased GST density, while AaGSW2 knockdown lines showed impaired GST initiation. Ectopic expression of AaGSW2 homologs from two mint cultivars, Mentha spicata and Mentha haplocalyx, in A. annua also induced GST formation. These results reveal a molecular mechanism involving homeodomain and WRKY proteins that controls glandular trichome initiation, at least part of which is shared by A. annua and mint.

Gene Expression Responses to Sequential Nutrient Deficiency Stresses in Soybean

Throughout the growing season, crops experience a multitude of short periods of various abiotic stresses. These stress events have long-term impacts on plant performance and yield. It is imperative to improve our understanding of the genes and biological processes underlying plant stress tolerance to mitigate end of season yield loss. The majority of studies examining transcriptional changes induced by stress focus on single stress events. Few studies have been performed in model or crop species to examine transcriptional responses of plants exposed to repeated or sequential stress exposure, which better reflect field conditions. In this study, we examine the transcriptional profile of soybean plants exposed to iron deficiency stress followed by phosphate deficiency stress (-Fe-Pi). Comparing this response to previous studies, we identified a core suite of genes conserved across all repeated stress exposures (-Fe-Pi, -Fe-Fe, -Pi-Pi). Additionally, we determined transcriptional response to sequential stress exposure (-Fe-Pi) involves genes usually associated with reproduction, not stress responses. These findings highlight the plasticity of the plant transcriptome and the complexity of unraveling stress response pathways.

The HD-ZIP IV transcription factor Tril regulates fruit spine density through gene dosage effects in cucumber

Trichomes and fruit spines are important traits that directly affect the appearance quality and commercial value of cucumber (Cucumis sativus). Tril (Trichome-less), encodes a HD-Zip IV transcription factor that plays a crucial role in the initiation of trichomes and fruit spines, but little is known about the details of the regulatory mechanisms involved. In this study, analysis of tissue expression patterns indicated that Tril is expressed and functions in the early stages of organ initiation and development. Expression of Tril under the control of its own promoter (the TrilPro::Tril-3*flag fragment) could partly rescue the mutant phenotypes of tril, csgl3 (cucumber glabrous 3, an allelic mutant of tril), and fs1 (few spines 1, a fragment substitution in the Tril promoter region), providing further evidence that Tril is responsible for the initiation of trichomes and fruit spines. In lines with dense spine, fs1-type lines, and transgenic lines of different backgrounds containing the TrilPro::Tril-3*flag foreign fragment, spine density increased in conjunction with increases in Tril expression, indicating that Tril has a gene dosage effect on fruit spine density in cucumber. Numerous Spines (NS) is a negative regulatory factor of fruit spine density. Characterization of the molecular and genetic interaction between Tril and NS/ns demonstrated that Tril functions upstream of NS with respect to spine initiation. Overall, our results reveal a novel regulatory mechanism governing the effect of Tril on fruit spine development, and provide a reference for future work on breeding for physical quality in cucumber.

Genetic control of tracheid properties in Norway spruce wood

Through the use of genome-wide association studies (GWAS) mapping it is possible to establish the genetic basis of phenotypic trait variation. Our GWAS study presents the first such effort in Norway spruce (Picea abies (L). Karst.) for the traits related to wood tracheid characteristics. The study employed an exome capture genotyping approach that generated 178 101 Single Nucleotide Polymorphisms (SNPs) from 40 018 probes within a population of 517 Norway spruce mother trees. We applied a least absolute shrinkage and selection operator (LASSO) based association mapping method using a functional multi-locus mapping approach, with a stability selection probability method as the hypothesis testing approach to determine significant Quantitative Trait Loci (QTLs). The analysis has provided 30 significant associations, the majority of which show specific expression in wood-forming tissues or high ubiquitous expression, potentially controlling tracheids dimensions, their cell wall thickness and microfibril angle. Among the most promising candidates based on our results and prior information for other species are: Picea abies BIG GRAIN 2 (PabBG2) with a predicted function in auxin transport and sensitivity, and MA_373300g0010 encoding a protein similar to wall-associated receptor kinases, which were both associated with cell wall thickness. The results demonstrate feasibility of GWAS to identify novel candidate genes controlling industrially-relevant tracheid traits in Norway spruce.

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.

Genome wide identification, characterization and expression analysis of HD-ZIP gene family in Cucumis sativus L. under biotic and various abiotic stresses

Information retrieved from genomic assembly may provide important clues and various molecular aspects in plants. Our research identified 40 CsHDZ genes in the Cucumber genome database. Subsequently; we performed the conserved motif and domain analysis of CsHDZ proteins. The phylogeny of the CsHDZ proteins further divides into 4 subfamilies (HD-ZIP I, HD-ZIP II, HD-ZIP III, and HD-ZIP IV) based on the structural similarities and functional diversities. The GO (Gene ontology) analysis of CsHDZ proteins showed that they are responsive to environmental stimuli and involved in numerous growth and developmental processes. The qRT-PCR analysis of 11 CsHDZ genes showed that they are expressed in all the tested tissues of Cucumis sativus. The differential expression pattern of CsHDZ genes unfolded their possible involvement in responding to various abiotic stresses and powdery mildew stress. It has been found that the CsHDZ22 localized in the nucleus which possibly participates in the regulatory mechanisms of various biological and cellular processes. In the light of above-mentioned outcomes, it has been deducted that CsHDZ genes in the Cucumis sativus genome play an important role in mediating the resistance to various abiotic stresses and powdery mildew stress as well as provide significant clues for functional studies.

The RNA-seq transcriptomic analysis reveals genes mediating salt tolerance through rapid triggering of ion transporters in a mutant barley

Considering the complex nature of salinity tolerance mechanisms, the use of isogenic lines or mutants possessing the same genetic background albeit different tolerance to salinity is a suitable method for reduction of analytical complexity to study these mechanisms. In the present study, whole transcriptome analysis was evaluated using RNA-seq method between a salt-tolerant mutant line "M4-73-30" and its wild-type "Zarjou" cultivar at seedling stage after six hours of exposure to salt stress (300 mM NaCl). Transcriptome sequencing yielded 20 million reads for each genotype. A total number of 7116 transcripts with differential expression were identified, 1586 and 1479 of which were obtained with significantly increased expression in the mutant and the wild-type, respectively. In addition, the families of WRKY, ERF, AP2/EREBP, NAC, CTR/DRE, AP2/ERF, MAD, MIKC, HSF, and bZIP were identified as the important transcription factors with specific expression in the mutant genotype. The RNA-seq results were confirmed at several time points using qRT-PCR for some important salt-responsive genes. In general, the results revealed that the mutant accumulated higher levels of sodium ion in the root and decreased its transfer to the shoot. Also, the mutant increased the amount of potassium ion leading to the maintenance a high ratio [K+]/[Na+] in the shoot compared to its wild-type via fast stomata closure and consequently transpiration reduction under the salt stress. Moreover, a reduction in photosynthesis and respiration was observed in the mutant, resulting in utilization of the stored energy and the carbon for maintaining the plant tissues, which is considered as a mechanism of salt tolerance in plants. Up-regulation of catalase, peroxidase, and ascorbate peroxidase genes has resulted in higher accumulation of H2O2 in the wild-type compared to the mutant. Therefore, the wild-type initiated rapid ROS signals which led to less oxidative scavenging in comparison with the mutant. The mutant increased expression in the ion transporters and the channels related to the salinity to maintain the ion homeostasis. In overall, the results demonstrated that the mutant responded better to the salt stress under both osmotic and ionic stress phases and lower damage was observed in the mutant compared to its wild-type under the salt stress.

Genome-wide identification and expression profiling of HD-ZIP gene family in Medicago truncatula

The homeodomain-leucine zipper (HD-ZIP) transcription factors are important regulators in various developmental processes and responses to environmental stimuli. Currently, little information is available for HD-ZIP gene family in Medicago truncatula. Here we perform a genome-wide analysis of HD-ZIP gene family in M. truncatula. Totally 52 M. truncatula HD-ZIPs (MtHDZs) were identified and classified into four distinctive subfamilies (I to IV). Members clustered in the same subfamily shared similar gene structure and protein motifs. Fifty-one MtHDZs were non-evenly distributed on eight chromosomes. Segmental duplication and purifying selection mainly contributed to the expansion and retention of M. truncatula HD-ZIP gene family. Expression profiling using the publicly available microarray data revealed that MtHDZ genes exhibited distinctive tissue-specific patterns and divergent responses to drought and salt stresses. In addition, the expression profile between each paralogous pair diverged differentially. Our results identified potential targets for the genetic improvement of abiotic stress tolerance in Medicago.

The role of HD-Zip class I transcription factors in plant response to abiotic stresses

Abiotic stresses usually affect plant growth and development, indirectly or directly causing crop production reduction and even plant death. To survive, plants utilize different mechanisms to adapt themselves to continuously changing surrounding environmental stresses. Homeodomain-leucine zipper (HD-Zip) transcription factors are unique to the plant kingdom and divided into four different subfamilies (HD-Zip I∼IV). Many HD-Zip I members have been shown to play critical roles in the regulation of plant developmental processes, signaling networks and responses to environmental stresses. This review focuses on the role of HD-Zip I transcription factors in plant responses to various abiotic stresses, including abscisic acid-mediated stress, drought and cold stress, oxidative stress, helping to identify the potential regulatory mechanisms that alleviate abiotic stress in plants.

Wheat wounding-responsive HD-Zip IV transcription factor GL7 is predominantly expressed in grain and activates genes encoding defensins

Several classes of transcription factors are involved in the activation of defensins. A new type of the transcription factor responsible for the regulation of wheat grain specific defensins was characterised in this work. HD-Zip class IV transcription factors constitute a family of multidomain proteins. A full-length cDNA of HD-Zip IV, designated TaGL7 was isolated from the developing grain of bread wheat, using a specific DNA sequence as bait in the Y1H screen. 3D models of TaGL7 HD complexed with DNA cis-elements rationalised differences that underlined accommodations of binding and non-binding DNA, while the START-like domain model predicted binding of lipidic molecules inside a concave hydrophobic cavity. The 3'-untranslated region of TaGL7 was used as a probe to isolate the genomic clone of TdGL7 from a BAC library prepared from durum wheat. The spatial and temporal activity of the TdGL7 promoter was tested in transgenic wheat, barley and rice. TdGL7 was expressed mostly in ovary at fertilisation and its promoter was active in a liquid endosperm during cellularisation and later in the endosperm transfer cells, aleurone, and starchy endosperm. The pattern of TdGL7 expression resembled that of genes that encode grain-specific lipid transfer proteins, particularly defensins. In addition, GL7 expression was upregulated by mechanical wounding, similarly to defensin genes. Co-bombardment of cultured wheat cells with TdGL7 driven by constitutive promoter and seven grain or root specific defensin promoters fused to GUS gene, revealed activation of four promoters. The data confirmed the previously proposed role of HD-Zip IV transcription factors in the regulation of genes that encode lipid transfer proteins involved in lipid transport and defence. The TdGL7 promoter could be used to engineer cereal grains with enhanced resistance to insects and fungal infections.

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.

Transcription factors involved in abiotic stress responses in Maize (Zea mays L.) and their roles in enhanced productivity in the post genomics era

BACKGROUND

Maize (Zea mays L.) is a principal cereal crop cultivated worldwide for human food, animal feed, and more recently as a source of biofuel. However, as a direct consequence of water insufficiency and climate change, frequent occurrences of both biotic and abiotic stresses have been reported in various regions around the world, and recently, this has become a constant threat in increasing global maize yields. Plants respond to abiotic stresses by utilizing the activities of transcription factors (TFs), which are families of genes coding for specific TF proteins. TF target genes form a regulon that is involved in the repression/activation of genes associated with abiotic stress responses. Therefore, it is of utmost importance to have a systematic study on each TF family, the downstream target genes they regulate, and the specific TF genes involved in multiple abiotic stress responses in maize and other staple crops.

METHOD

In this review, the main TF families, the specific TF genes and their regulons that are involved in abiotic stress regulation will be briefly discussed. Great emphasis will be given on maize abiotic stress improvement throughout this review, although other examples from different plants like rice, Arabidopsis, wheat, and barley will be used.

RESULTS

We have described in detail the main TF families in maize that take part in abiotic stress responses together with their regulons. Furthermore, we have also briefly described the utilization of high-efficiency technologies in the study and characterization of TFs involved in the abiotic stress regulatory networks in plants with an emphasis on increasing maize production. Examples of these technologies include next-generation sequencing, microarray analysis, machine learning, and RNA-Seq.

CONCLUSION

In conclusion, it is expected that all the information provided in this review will in time contribute to the use of TF genes in the research, breeding, and development of new abiotic stress tolerant maize cultivars.

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.

Genome-wide Identification, Classification, and Expression Pattern of Homeobox Gene Family in Brassica rapa under Various Stresses

Homeobox (HB) genes are crucial for plant growth and development processes. They encode transcription factors and responses to various stresses, as reported by recent emerging evidence. In this study, a total of 113 BraHB genes were identified in Brassica rapa. On the basis of domain organization and phylogenetic analysis, the BraHBs were grouped into nine subclasses, in which homeobox leucine-zipper (HB LZP-III) showed the highest number of genes (28) compared to other subclasses. The BraHBs exhibited similarities in exon-intron organization and motif composition among the members of the same subclasses. The analysis revealed that HB-Knotted was more preferentially retained than any other subclass of BraHB. Furthermore, we evaluated the impact of whole-genome triplication on the evolution of BraHBs. In order to analyze the subgenomes of B. rapa, we identified 39 paralogous pairs for which synonymous substitution values were lower than 1.00 for further purifying selection. Finally, the expression patterns of BraHBs across six tissues expressed dynamic variations combined with their responses against multiple stresses. The current study provides brief information on the homeobox gene family in B. rapa. Our findings can serve as a reference for further functional analysis of BraHBs.

Overexpression of the class I homeodomain transcription factor TaHDZipI-5 increases drought and frost tolerance in transgenic wheat

Characterization of the function of stress-related genes helps to understand the mechanisms of plant responses to environmental conditions. The findings of this work defined the role of the wheat TaHDZipI-5 gene, encoding a stress-responsive homeodomain-leucine zipper class I (HD-Zip I) transcription factor, during the development of plant tolerance to frost and drought. Strong induction of TaHDZipI-5 expression by low temperatures, and the elevated TaHDZipI-5 levels of expression in flowers and early developing grains in the absence of stress, suggests that TaHDZipI-5 is involved in the regulation of frost tolerance at flowering. The TaHDZipI-5 protein behaved as an activator in a yeast transactivation assay, and the TaHDZipI-5 activation domain was localized to its C-terminus. The TaHDZipI-5 protein homo- and hetero-dimerizes with related TaHDZipI-3, and differences between DNA interactions in both dimers were specified at 3D molecular levels. The constitutive overexpression of TaHDZipI-5 in bread wheat significantly enhanced frost and drought tolerance of transgenic wheat lines with the appearance of undesired phenotypic features, which included a reduced plant size and biomass, delayed flowering and a grain yield decrease. An attempt to improve the phenotype of transgenic wheat by the application of stress-inducible promoters with contrasting properties did not lead to the elimination of undesired phenotype, apparently due to strict spatial requirements for TaHDZipI-5 overexpression.

A novel HD-ZIP IV/MIXTA complex promotes glandular trichome initiation and cuticle development in Artemisia annua

Glandular trichomes and cuticles are both specialized structures that cover the epidermis of aerial plant organs. The former are commonly regarded as 'biofactories' for producing valuable natural products. The latter are generally considered as natural barriers for defending plants against abiotic and biotic stresses. However, the regulatory network for their formation and relationship remains largely elusive. Here we identify a homeodomain-leucine zipper (HD-ZIP) IV transcription factor, AaHD8, directly promoting the expression of AaHD1 for glandular trichome initiation in Artemisia annua. We found that AaHD8 positively regulated leaf cuticle development in A. annua via controlling the expression of cuticle-related enzyme genes. Furthermore, AaHD8 interacted with a MIXTA-like protein AaMIXTA1, a positive regulator of trichome initiation and cuticle development, forming a regulatory complex and leading to enhanced transcriptional activity in regulating the expression of AaHD1 and cuticle development genes. Our results reveal a molecular mechanism by which a novel HD-ZIP IV/MIXTA complex plays a significant role in regulating epidermal development, including glandular trichome initiation and cuticle formation.

Arabidopsis ACYL-COA-BINDING PROTEIN1 interacts with STEROL C4-METHYL OXIDASE1-2 to modulate gene expression of homeodomain-leucine zipper IV transcription factors

Fatty acids (FAs) and sterols constitute building blocks of eukaryotic membranes and lipid signals. Co-regulation of FA and sterol synthesis is mediated by sterol regulatory element-binding proteins in animals but remains elusive in plants. We reported recently that Arabidopsis ACYL-COA-BINDING PROTEIN1 (ACBP1) modulates sterol synthesis via protein-protein interaction with STEROL C4-METHYL OXIDASE1-1 (SMO1-1). Herein, ACBP1 was demonstrated to co-express and interact with SMO1-2 by yeast two-hybrid, co-localization, pull-down, co-immunoprecipitation and β-glucuronidase assays. SMO1-2 silenced in acbp1 was used in phenotyping, GC-MS and expression profiling. ACBP1 co-expressed with SMO1-2 in embryo sacs, pollen and trichomes, corroborating with cooperative tissue-specific functions unseen with SMO1-1. SMO1-2 silencing in acbp1 impaired seed development, male and female gamete transmission, and pollen function. Genes encoding homeodomain-leucine zipper IV transcription factors (HDG5, HDG10, HDG11 and GLABRA2), which potentially bind phospholipids/sterols, were transcribed aberrantly. GLABRA2 targets (MYB23, MUM4 and PLDα1) were misregulated, causing glabra2-resembling trichome, seed coat mucilage and oil-accumulating phenotypes. Together with altered sterol and FA compositions upon ACBP1 mutation and/or SMO1-2 silencing, ACBP1-SMO1 interaction appears to mediate homeostatic co-regulation of FAs and sterols, which serve as lipid modulators for gene expression of homeodomain-leucine zipper IV transcription factors.

Survey of High Throughput RNA-Seq Data Reveals Potential Roles for lncRNAs during Development and Stress Response in Bread Wheat

Long non-coding RNAs (lncRNAs) are a family of regulatory RNAs that play essential role in the various developmental processes and stress responses. Recent advances in sequencing technology and computational methods enabled identification and characterization of lncRNAs in certain plant species, but they are less known in Triticum aestivum (bread wheat). Herein, we analyzed 52 RNA seq data (>30 billion reads) and identified 44,698 lncRNAs in T. aestivum genome, which were characterized in comparison to the coding sequences (mRNAs). Similar to the mRNAs, lncRNAs were also derived from each sub-genome and chromosome, and showed tissue developmental stage specific and differential expression, as well. The modulated expression of lncRNAs during abiotic stresses like heat, drought, and salt indicated their putative role in stress response. The co-expression of lncRNAs with vital mRNAs including various transcription factors and enzymes involved in Abscisic acid (ABA) biosynthesis, and gene ontology mapping inferred their regulatory roles in numerous biological processes. A few lncRNAs were predicted as precursor (19 lncRNAs), while some as target mimics (1,047 lncRNAs) of known miRNAs involved in various regulatory functions. The results suggested numerous functions of lncRNAs in T. aestivum, and unfolded the opportunities for functional characterization of individual lncRNA in future studies.

The Vascular Pathogen Verticillium longisporum Does Not Affect Water Relations and Plant Responses to Drought Stress of Its Host, Brassica napus

Verticillium longisporum is a host-specific vascular pathogen of oilseed rape (Brassica napus L.) that causes economic crop losses by impairing plant growth and inducing premature senescence. This study investigates whether plant damage through Verticillium stem striping is due to impaired plant water relations, whether V. longisporum affects responses of a susceptible B. napus variety to drought stress, and whether drought stress, in turn, affects plant responses to V. longisporum. Two-factorial experiments on a susceptible cultivar of B. napus infected or noninfected with V. longisporum and exposed to three watering levels (30, 60, and 100% field capacity) revealed that drought stress and V. longisporum impaired plant growth by entirely different mechanisms. Although both stresses similarly affected plant growth parameters (plant height, hypocotyl diameter, and shoot and root dry matter), infection of B. napus with V. longisporum did not affect any drought-related physiological or molecular genetic plant parameters, including transpiration rate, stomatal conductance, photosynthesis rate, water use efficiency, relative leaf water content, leaf proline content, or the expression of drought-responsive genes. Thus, this study provides comprehensive physiological and molecular genetic evidence explaining the lack of wilt symptoms in B. napus infected with V. longisporum. Likewise, drought tolerance of B. napus was unaffected by V. longisporum, as was the level of disease by drought conditions, thus excluding a concerted action of both stresses in the field. Although it is evident that drought and vascular infection with V. longisporum impair plant growth by different mechanisms, it remains to be determined by which other factors V. longisporum causes crop loss.

Genome-wide characterization and expression profiling of HD-Zip gene family related to abiotic stress in cassava

Homeodomain-leucine zipper (HD-Zip) gene family plays important roles in various abiotic stresses and hormone signaling in plants. However, no information is currently available regarding this family in cassava (Manihot esculenta), an important drought-tolerant crop in tropical and sub-tropical areas. Here, 57 HD-Zip genes (MeHDZ01-57) were identified in the cassava genome, and they were classified into four subfamilies based on phylogenetic analysis, which was further supported by their gene structure and conserved motif characteristics. Of which five gene pairs were involved in segmental duplication but none for tandem duplication, suggesting that segmental duplication was the main cause for the expansion of MeHDZ gene family in cassava. Global expression profiles revealed that MeHDZ genes were constitutively expressed, or not expressed, or tissue-specific expressed in examined tissues in both cultivated and wild subspecies. Transcriptomic analysis of three genotypes showed that most of MeHDZ genes responded differently to drought and polyethylene glycol treatments. Subsequently, quantitative RT-PCR analysis revealed comprehensive responses of twelve selected MeHDZ genes to various stimuli including cold, salt, and ABA treatments. These findings will increase our understanding of HD-Zip gene family involved in abiotic stresses and signaling transduction, and will provide a solid base for further functional characterization of MeHDZ genes in cassava.

HOMEODOMAIN PROTEIN 1 is required for jasmonate-mediated glandular trichome initiation in Artemisia annua

Glandular trichomes are generally considered biofactories that produce valuable chemicals. Increasing glandular trichome density is a very suitable way to improve the productivity of these valuable metabolites, but little is known about the regulation of glandular trichome formation. Phytohormone jasmonate (JA) promotes glandular trichome initiation in various plants, but its mechanism is also unknown. By searching transcription factors regulated by JA in Artemisia annua, we identified a novel homeodomain-leucine zipper transcription factor, HOMEODOMAIN PROTEIN 1 (AaHD1), which positively controls both glandular and nonglandular trichome initiations. Overexpression of AaHD1 in A. annua significantly increased glandular trichome density without harming plant growth. Consequently, the artemisinin content was improved. AaHD1 interacts with A. annua jasmonate ZIM-domain 8 (AaJAZ8), which is a repressor of JA, thereby resulting in decreased transcriptional activity. AaHD1 knockdown lines show decreased sensitivity to JA on glandular trichome initiation, which indicates that AaHD1 plays an important role in JA-mediated glandular trichome initiation. We identified a new transcription factor that promotes A. annua glandular trichome initiation and revealed a novel molecular mechanism by which a homeodomain protein transduces JA signal to promote glandular trichome initiation. Our results also suggested a connection between glandular and nonglandular trichome formations.

Over-Expression of OsHOX24 Confers Enhanced Susceptibility to Abiotic Stresses in Transgenic Rice via Modulating Stress-Responsive Gene Expression

Homeobox transcription factors play critical roles in plant development and abiotic stress responses. In the present study, we raised rice transgenics over-expressing stress-responsive OsHOX24 gene (rice homeodomain-leucine zipper I sub-family member) and analyzed their response to various abiotic stresses at different stages of development. At the seed germination stage, rice transgenics over-expressing OsHOX24 exhibited enhanced sensitivity to abiotic stress conditions and abscisic acid as compared to wild-type (WT). OsHOX24 over-expression rice seedlings showed reduced root and shoot growth under salinity and desiccation stress (DS) conditions. Various physiological and phenotypic assays confirmed higher susceptibility of rice transgenics toward abiotic stresses as compared to WT at mature and reproductive stages of rice development too. Global gene expression profiling revealed differential regulation of several genes in the transgenic plants under control and DS conditions. Many of these differentially expressed genes were found to be involved in transcriptional regulatory activities, besides carbohydrate, nucleic acid and lipid metabolic processes and response to abiotic stress and hormones. Taken together, our findings highlighted the role of OsHOX24 in regulation of abiotic stress responses via modulating the expression of stress-responsive genes in rice.

GL2-type homeobox gene Roc4 in rice promotes flowering time preferentially under long days by repressing Ghd7

Under long day (LD) lengths, flowering can be delayed in rice by modulating several regulatory genes. We found activation tagging lines that showed an early flowering phenotype preferentially under LD conditions. Expression of Rice outermost cell-specific gene 4 (Roc4), encoding a homeodomain Leu-zipper class IV family protein, was significantly increased. Transcript levels of Grain number, plant height, and heading date7 (Ghd7) were significantly reduced while those of Ghd7 downstream genes were increased. However, other flowering regulators were unaffected. Whereas constitutive overexpression of Roc4 in 'Dongjin' japonica rice, which carries active Ghd7, also caused LD-preferential early flowering, its overexpression in 'Longjing27' rice, which is defective in functional Ghd7, did not produce the same result. This confirmed that Roc4 regulates flowering time mainly through Ghd7. Phytochromes and O. sativa GIGANTEA (OsGI) function upstream of Roc4. Transgenic plants showed ubiquitous expression of the β-glucuronidase reporter gene under the Roc4 promoter. Furthermore, Roc4 had transcriptional activation activity in the N-terminal region of the StAR-related lipid-transfer domain. All of these findings are evidence that Roc4 is an LD-preferential flowering enhancer that functions downstream of phytochromes and OsGI, but upstream of Ghd7.

Are We on the Right Track: Can Our Understanding of Abscission in Model Systems Promote or Derail Making Improvements in Less Studied Crops?

As the world population grows and resources and climate conditions change, crop improvement continues to be one of the most important challenges for agriculturalists. The yield and quality of many crops is affected by abscission or shattering, and environmental stresses often hasten or alter the abscission process. Understanding this process can not only lead to genetic improvement, but also changes in cultural practices and management that will contribute to higher yields, improved quality and greater sustainability. As plant scientists, we have learned significant amounts about this process through the study of model plants such as Arabidopsis, tomato, rice, and maize. While these model systems have provided significant valuable information, we are sometimes challenged to use this knowledge effectively as variables including the economic value of the crop, the uniformity of the crop, ploidy levels, flowering and crossing mechanisms, ethylene responses, cultural requirements, responses to changes in environment, and cellular and tissue specific morphological differences can significantly influence outcomes. The value of genomic resources for lesser-studied crops such as cranberries and grapes and the orphan crop fonio will also be considered.

Transcriptional Responses in Root and Leaf of Prunus persica under Drought Stress Using RNA Sequencing

Prunus persica L. Batsch, or peach, is one of the most important crops and it is widely established in irrigated arid and semi-arid regions. However, due to variations in the climate and the increased aridity, drought has become a major constraint, causing crop losses worldwide. The use of drought-tolerant rootstocks in modern fruit production appears to be a useful method of alleviating water deficit problems. However, the transcriptomic variation and the major molecular mechanisms that underlie the adaptation of drought-tolerant rootstocks to water shortage remain unclear. Hence, in this study, high-throughput sequencing (RNA-seq) was performed to assess the transcriptomic changes and the key genes involved in the response to drought in root tissues (GF677 rootstock) and leaf tissues (graft, var. Catherina) subjected to 16 days of drought stress. In total, 12 RNA libraries were constructed and sequenced. This generated a total of 315 M raw reads from both tissues, which allowed the assembly of 22,079 and 17,854 genes associated with the root and leaf tissues, respectively. Subsets of 500 differentially expressed genes (DEGs) in roots and 236 in leaves were identified and functionally annotated with 56 gene ontology (GO) terms and 99 metabolic pathways, which were mostly associated with aminobenzoate degradation and phenylpropanoid biosynthesis. The GO analysis highlighted the biological functions that were exclusive to the root tissue, such as "locomotion," "hormone metabolic process," and "detection of stimulus," indicating the stress-buffering role of the GF677 rootstock. Furthermore, the complex regulatory network involved in the drought response was revealed, involving proteins that are associated with signaling transduction, transcription and hormone regulation, redox homeostasis, and frontline barriers. We identified two poorly characterized genes in P. persica: growth-regulating factor 5 (GRF5), which may be involved in cellular expansion, and AtHB12, which may be involved in root elongation. The reliability of the RNA-seq experiment was validated by analyzing the expression patterns of 34 DEGs potentially involved in drought tolerance using quantitative reverse transcription polymerase chain reaction. The transcriptomic resources generated in this study provide a broad characterization of the acclimation of P. persica to drought, shedding light on the major molecular responses to the most important environmental stressor.

Genome-Wide Identification and Expression Analysis of Homeodomain Leucine Zipper Subfamily IV (HDZ IV) Gene Family from Musa accuminata

The homeodomain zipper family (HD-ZIP) of transcription factors is present only in plants and plays important role in the regulation of plant-specific processes. The subfamily IV of HDZ transcription factors (HD-ZIP IV) has primarily been implicated in the regulation of epidermal structure development. Though this gene family is present in all lineages of land plants, members of this gene family have not been identified in banana, which is one of the major staple fruit crops. In the present work, we identified 21 HDZIV encoding genes in banana by the computational analysis of banana genome resource. Our analysis suggested that these genes putatively encode proteins having all the characteristic domains of HDZIV transcription factors. The phylogenetic analysis of the banana HDZIV family genes further confirmed that after separation from a common ancestor, the banana, and poales lineages might have followed distinct evolutionary paths. Further, we conclude that segmental duplication played a major role in the evolution of banana HDZIV encoding genes. All the identified banana HDZIV genes expresses in different banana tissue, however at varying levels. The transcript levels of some of the banana HDZIV genes were also detected in banana fruit pulp, suggesting their putative role in fruit attributes. A large number of genes of this family showed modulated expression under drought and salinity stress. Taken together, the present work lays a foundation for elucidation of functional aspects of the banana HDZIV encoding genes and for their possible use in the banana improvement programs.

Comprehensive analysis and expression profile of the homeodomain leucine zipper IV transcription factor family in tomato

Homeodomain leucine zipper IV (HD-ZIP IV) proteins are plant-specific transcription factors that play important roles in development of epidermal cell layers and cuticle formation. The functions of two HD-ZIP IV family genes, CD2 and Wo, have been well characterized in tomato (Solanum lycopersicum). CD2 and Wo are involved in cuticle biosynthesis and trichome formation, respectively. In this study, we identified 13 novel tomato HD-ZIP IV (SlHDZIV) genes. We analyzed the structures, chromosome locations, phylogeny, protein motifs, and expression profiles of these SlHDZIV genes. Gene structure analysis revealed that a module of 11 exons and 10 introns existed in the SlHDZIV genes. These genes were asymmetrically distributed on chromosomes, except on chromosome 4 and 5. Segmental duplication possibly contributed to the expansion of tomato HD-ZIP IV genes. The expression profiles of these genes revealed their broad expression pattern and high expression in young leaves and flowers. Each gene responded to more than one of different phytohormones [abscisic acid, ethephon, 4-(indolyl)-butyric acid, jasmonic acid, salicylic acid, gibberellic acid, and 6-benzylaminopurine] and four abiotic stress treatments (cold, heat, salt, and drought). This study provided significant insights into the diverse roles of SlHDZIV genes in tomato growth and development.