Characterization of the class IV homeodomain-Leucine Zipper gene family in Arabidopsis

GLABRA2 transcription factor integrates arsenic tolerance with epidermal cell fate determination

Arsenic poses a global threat to living organisms, compromising crop security and yield. Limited understanding of the transcriptional network integrating arsenic-tolerance mechanisms with plant developmental responses hinders the development of strategies against this toxic metalloid. Here, we conducted a high-throughput yeast one-hybrid assay using as baits the promoter region from the arsenic-inducible genes ARQ1 and ASK18 from Arabidopsis thaliana, coupled with a transcriptomic analysis, to uncover novel transcriptional regulators of the arsenic response. We identified the GLABRA2 (GL2) transcription factor as a novel regulator of arsenic tolerance, revealing a wider regulatory role beyond its established function as a repressor of root hair formation. Furthermore, we found that ANTHOCYANINLESS2 (ANL2), a GL2 subfamily member, acts redundantly with this transcription factor in the regulation of arsenic signaling. Both transcription factors act as repressors of arsenic response. gl2 and anl2 mutants exhibit enhanced tolerance and reduced arsenic accumulation. Transcriptional analysis in the gl2 mutant unveils potential regulators of arsenic tolerance. These findings highlight GL2 and ANL2 as novel integrators of the arsenic response with developmental outcomes, offering insights for developing safer crops with reduced arsenic content and increased tolerance to this hazardous metalloid.

Genome-wide analysis of HD-Zip genes in Sophora alopecuroides and their role in salt stress response

We aimed to identify HD-Zip (homologous domain leucine zipper) family genes based on the complete Sophora alopecuroides genome sequence. Eighty-six Sophora alopecuroides HD-Zip family (SaHDZ) genes were identified and categorized into four subclasses using phylogenetic analysis. Chromosome localization analysis revealed that these genes were distributed across 18 chromosomes. Gene structure and conserved motif analysis showed high similarity among members of the SaHDZ genes. Prediction analysis revealed 71 cis-acting elements in SaHDZ genes. Transcriptome and quantitative real-time polymerase chain reaction analyses showed that under salt stress, SaHDZ responded positively in S. alopecuroides, and that SaHDZ22 was significantly upregulated afterward. Functional verification experiments revealed that SaHDZ22 overexpression increased the tolerance of Arabidopsis to salt and osmotic stress. Combined with cis-acting element prediction and expression level analysis, HD-Zip family transcription factors may be involved in regulating the balance between plant growth and stress resistance under salt stress by modulating the expression of auxin and abscisic acid signaling pathway genes. The Sophora alopecuroides adenylate kinase protein (SaAKI) and S. alopecuroides tetrapeptide-like repeat protein (SaTPR; pCAMBIA1300-SaTPR-cLUC) expression levels were consistent with those of SaHDZ22, indicating that SaHDZ22 may coordinate with SaAKI and SaTPR to regulate plant salt tolerance. These results lay a foundation in understanding the salt stress response mechanisms of S. alopecuroides and provide a reference for future studies oriented toward exploring plant stress resistance.

QTL mapping for pod quality and yield traits in snap bean (Phaseolus vulgaris L.)

Pod quality and yield traits in snap bean (Phaseolus vulgaris L.) influence consumer preferences, crop adoption by farmers, and the ability of the product to be commercially competitive locally and globally. The objective of the study was to identify the quantitative trait loci (QTL) for pod quality and yield traits in a snap × dry bean recombinant inbred line (RIL) population. A total of 184 F6 RILs derived from a cross between Vanilla (snap bean) and MCM5001 (dry bean) were grown in three field sites in Kenya and one greenhouse environment in Davis, CA, USA. They were genotyped at 5,951 single nucleotide polymorphisms (SNPs), and composite interval mapping was conducted to identify QTL for 16 pod quality and yield traits, including pod wall fiber, pod string, pod size, and harvest metrics. A combined total of 44 QTL were identified in field and greenhouse trials. The QTL for pod quality were identified on chromosomes Pv01, Pv02, Pv03, Pv04, Pv06, and Pv07, and for pod yield were identified on Pv08. Co-localization of QTL was observed for pod quality and yield traits. Some identified QTL overlapped with previously mapped QTL for pod quality and yield traits, with several others identified as novel. The identified QTL can be used in future marker-assisted selection in snap bean.

Deciphering the dynamic expression network of fiber elongation and the functional role of the GhTUB5 gene for fiber length in cotton based on an introgression population of upland cotton

INTRODUCTION

Interspecific introgression between Gossypium hirsutum and G. barbadense allows breeding cotton with outstanding fiber length (FL). However, the dynamic gene regulatory network of FL-related genes has not been characterized, and the functional mechanism through which the hub gene GhTUB5 mediates fiber elongation has yet to be determined.

METHODS

Coexpression analyses of 277 developing fiber transcriptomes integrated with QTL mapping using 250 introgression lines of different FL phenotypes were conducted to identify genes related to fiber elongation. The function of GhTUB5 was determined by ectopic expression of two TUB5 alleles in Arabidopsis and knockout of GhTUB5 in upland cotton. Yeast two-hybrid, split-luciferase and pull-down assays were conducted to screen for interacting proteins, and upstream genes were identified by yeast one-hybrid, dual-LUC and electrophoretic mobility shift assays.

RESULTS

The 32,612, 30,837 and 30,277 genes expressed at 5, 10 and 15 days postanthesis (dpa) were grouped into 19 distinct coexpression modules, and 988 genes in the MEblack module were enriched in the cell wall process and exhibited significant associations with FL. A total of 20 FL-QTLs were identified, each explaining 3.34-16.04 % of the phenotypic variance in the FL. Furthermore, several FL-QTLs contained 15 genes that were differentially expressed in the MEblack module including the tubulin beta gene (TUB5). Compared with the wild type, the overexpression of GhTUB5 and GbTUB5 in Arabidopsis suppressed root cell length but promoted cellulose synthesis. Knockout of GhTUB5 resulted in longer fiber lines. Protein-based experiments revealed that GhTUB5 interacts with GhZFP6. Additionally, GhTUB5 was directly activated by GhHD-ZIP7, a homeobox-leucine zipper transcription factor, and its paralogous gene was previously reported to mediate fiber elongation.

CONCLUSION

This study opens a new avenue to dissect functional mechanism of cotton fiber elongation. Our findings provide some molecular details on how GhTUB5 mediates the FL phenotype in cotton.

The BTB-BACK-TAZ domain protein MdBT2 reduces drought resistance by weakening the positive regulatory effect of MdHDZ27 on apple drought tolerance via ubiquitination

Drought stress is one of the dominating challenges to the growth and productivity in crop plants. Elucidating the molecular mechanisms of plants responses to drought stress is fundamental to improve fruit quality. However, such molecular mechanisms are poorly understood in apple (Malus domestica Borkh.). In this study, we explored that the BTB-BACK-TAZ protein, MdBT2, negatively modulates the drought tolerance of apple plantlets. Moreover, we identified a novel Homeodomain-leucine zipper (HD-Zip) transcription factor, MdHDZ27, using a yeast two-hybrid (Y2H) screen with MdBT2 as the bait. Overexpression of MdHDZ27 in apple plantlets, calli, and tomato plantlets enhanced their drought tolerance by promoting the expression of drought tolerance-related genes [responsive to dehydration 29A (MdRD29A) and MdRD29B]. Biochemical analyses demonstrated that MdHDZ27 directly binds to and activates the promoters of MdRD29A and MdRD29B. Furthermore, in vitro and in vivo assays indicate that MdBT2 interacts with and ubiquitinates MdHDZ27, via the ubiquitin/26S proteasome pathway. This ubiquitination results in the degradation of MdHDZ27 and weakens the transcriptional activation of MdHDZ27 on MdRD29A and MdRD29B. Finally, a series of transgenic analyses in apple plantlets further clarified the role of the relationship between MdBT2 and MdHDZ27, as well as the effect of their interaction on drought resistance in apple plantlets. Collectively, our findings reveal a novel mechanism by which the MdBT2-MdHDZ27 regulatory module controls drought tolerance, which is of great significance for enhancing the drought resistance of apple and other plants.

NtHD9 modulates plant salt tolerance by regulating the formation of glandular trichome heads in Nicotiana tabacum

Salt stress is one of the main abiotic factor affecting plant growth. We have previously identified a key gene (NtHD9) in Nicotiana tabacum L. that positively regulates the formation of long glandular trichomes (LGTs). Here, we verified that both abiotic stress (aphids, drought and salt stress) could restore the phenotype lacking LGTs in NtHD9-knockout (NtHD9-KO) plants. The abiotic stress response assays indicated that NtHD9 is highly sensitive to salt stress. Compared with cultivated tobacco "K326" (CK) plants, NtHD9-overexpressing (NtHD9-OE) plants with more LGTs exhibited stronger salt tolerance, whereas NtHD9-KO with no LGTs showed weaker tolerance to salt. The densities and sizes of the glandular heads gradually increased with increasing NaCl concentrations in NtHD9-KO plants. Mineral element determination showed that leaves and trichomes of NtHD9-OE plants accumulated less Na+ but had higher K+ contents under salt stress, thus maintaining ion homeostasis in plants, which could contribute to a robust photosynthetic and antioxidant system under salt stress. Therefore, NtHD9-OE plants maintained a larger leaf area and root length under high-salt conditions than CK and NtHD9-KO plants. We verified that NtHD9 could individually interact with NtHD5, NtHD7, NtHD12, and NtJAZ10 proteins. Salt stress led to an increase in jasmonic acid (JA) levels and activated the expression of NtHDs while inhibiting the expression of NtJAZ. This study suggests that the glandular heads play an important role in plant resistance to salt stress. The activation of JA signaling leading to JAZ protein degradation may be key factors regulating the glandular heads development under salt stress.

Whole-Genome Sequencing and Analysis of Tumour-Forming Radish (Raphanus sativus L.) Line

Spontaneous tumour formation in higher plants can occur in the absence of pathogen invasion, depending on the plant genotype. Spontaneous tumour formation on the taproots is consistently observed in certain inbred lines of radish (Raphanus sativus var. radicula Pers.). In this paper, using Oxford Nanopore and Illumina technologies, we have sequenced the genomes of two closely related radish inbred lines that differ in their ability to spontaneously form tumours. We identified a large number of single nucleotide variants (amino acid substitutions, insertions or deletions, SNVs) that are likely to be associated with the spontaneous tumour formation. Among the genes involved in the trait, we have identified those that regulate the cell cycle, meristem activity, gene expression, and metabolism and signalling of phytohormones. After identifying the SNVs, we performed Sanger sequencing of amplicons corresponding to SNV-containing regions to validate our results. We then checked for the presence of SNVs in other tumour lines of the radish genetic collection and found the ERF118 gene, which had the SNVs in the majority of tumour lines. Furthermore, we performed the identification of the CLAVATA3/ESR (CLE) and WUSCHEL (WOX) genes and, as a result, identified two unique radish CLE genes which probably encode proteins with multiple CLE domains. The results obtained provide a basis for investigating the mechanisms of plant tumour formation and also for future genetic and genomic studies of radish.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genome-wide analysis and expression profiling of the HD-ZIP gene family in kiwifruit

The homeodomain-leucine zipper (HD-Zip) gene family plays a pivotal role in plant development and stress responses. Nevertheless, a comprehensive characterization of the HD-Zip gene family in kiwifruit has been lacking. In this study, we have systematically identified 70 HD-Zip genes in the Actinidia chinensis (Ac) genome and 55 in the Actinidia eriantha (Ae) genome. These genes have been categorized into four subfamilies (HD-Zip I, II, III, and IV) through rigorous phylogenetic analysis. Analysis of synteny patterns and selection pressures has provided insights into how whole-genome duplication (WGD) or segmental may have contributed to the divergence in gene numbers between these two kiwifruit species, with duplicated gene pairs undergoing purifying selection. Furthermore, our study has unveiled tissue-specific expression patterns among kiwifruit HD-Zip genes, with some genes identified as key regulators of kiwifruit responses to bacterial canker disease and postharvest processes. These findings not only offer valuable insights into the evolutionary and functional characteristics of kiwifruit HD-Zips but also shed light on their potential roles in plant growth and development.

Identification of a cis-element for long glandular trichome-specific gene expression, which is targeted by a HD-ZIP IV protein

Glandular trichomes are epidermal outgrowths that secret a variety of secondary metabolites, which not only help plants adapt to environmental stresses but also have important commercial value in fragrances, pharmaceuticals, and pesticides. In Nicotiana tabacum, it has been confirmed that a B-type cyclin, CycB2, negatively regulates the formation of long glandular trichomes (LGTs). This study aimed to identify the upstream regulatory gene involved in LGT formation by screening LGT-specific cis-elements within the NtCycB2 promoter. Using GUS as a reporter gene, the tissue-driven ability of NtCycB2 promoter showed that NtCycB2 promoter could drive GUS expression specifically in LGTs. Function analysis of a series of successive 5' truncations and synthetic segments of the NtCycB2 promoter indicated that the 87-bp region from -1221 to -1134 of the NtCycB2 promoter was required for gene expression in LGTs, and the L1-element (5'-AAAATTAATAAGAG-3') located in the 87-bp region contributed to the gene expression in the stalk of LGTs. Further Y1H and LUC assays confirmed that this L1-element exclusively binds to a HD-Zip IV protein, NtHD13. Gene function analysis revealed that NtHD13 positively controlled LGT formation, as overexpression of NtHD13 resulted in a high number of LGTs, whereas knockout of NtHD13 led to a decrease in LGTs. These findings demonstrate that NtHD13 can bind to an L1-element within the NtCycB2 promoter to regulate LGT formation.

Phosphatidic acid interacts with an HD-ZIP transcription factor GhHOX4 to influence its function in fiber elongation of cotton (Gossypium hirsutum)

Upland cotton, the mainly cultivated cotton species in the world, provides over 90% of natural raw materials (fibers) for the textile industry. The development of cotton fibers that are unicellular and highly elongated trichomes on seeds is a delicate and complex process. However, the regulatory mechanism of fiber development is still largely unclear in detail. In this study, we report that a homeodomain-leucine zipper (HD-ZIP) IV transcription factor, GhHOX4, plays an important role in fiber elongation. Overexpression of GhHOX4 in cotton resulted in longer fibers, while GhHOX4-silenced transgenic cotton displayed a "shorter fiber" phenotype compared with wild type. GhHOX4 directly activates two target genes, GhEXLB1D and GhXTH2D, for promoting fiber elongation. On the other hand, phosphatidic acid (PA), which is associated with cell signaling and metabolism, interacts with GhHOX4 to hinder fiber elongation. The basic amino acids KR-R-R in START domain of GhHOX4 protein are essential for its binding to PA that could alter the nuclear localization of GhHOX4 protein, thereby suppressing the transcriptional regulation of GhHOX4 to downstream genes in the transition from fiber elongation to secondary cell wall (SCW) thickening during fiber development. Thus, our data revealed that GhHOX4 positively regulates fiber elongation, while PA may function in the phase transition from fiber elongation to SCW formation by negatively modulating GhHOX4 in cotton.

Reciprocal expression of MADS-box genes and DNA methylation reconfiguration initiate bisexual cones in spruce

The naturally occurring bisexual cone of gymnosperms has long been considered a possible intermediate stage in the origin of flowers, but the mechanisms governing bisexual cone formation remain largely elusive. Here, we employed transcriptomic and DNA methylomic analyses, together with hormone measurement, to investigate the molecular mechanisms underlying bisexual cone development in the conifer Picea crassifolia. Our study reveals a "bisexual" expression profile in bisexual cones, especially in expression patterns of B-class, C-class and LEAFY genes, supporting the out of male model. GGM7 could be essential for initiating bisexual cones. DNA methylation reconfiguration in bisexual cones affects the expression of key genes in cone development, including PcDAL12, PcDAL10, PcNEEDLY, and PcHDG5. Auxin likely plays an important role in the development of female structures of bisexual cones. This study unveils the potential mechanisms responsible for bisexual cone formation in conifers and may shed light on the evolution of bisexuality.

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.

Epicuticular wax accumulation and regulation of wax pathway gene expression during bioenergy Sorghum stem development

Bioenergy sorghum is a drought-tolerant high-biomass C4 grass targeted for production on annual cropland marginal for food crops due primarily to abiotic constraints. To better understand the overall contribution of stem wax to bioenergy sorghum's resilience, the current study characterized sorghum stem cuticular wax loads, composition, morphometrics, wax pathway gene expression and regulation using vegetative phase Wray, R07020, and TX08001 genotypes. Wax loads on sorghum stems (~103-215 µg/cm2) were much higher than Arabidopsis stem and leaf wax loads. Wax on developing sorghum stem internodes was enriched in C28/30 primary alcohols (~65%) while stem wax on fully developed stems was enriched in C28/30 aldehydes (~80%). Scanning Electron Microscopy showed minimal wax on internodes prior to the onset of elongation and that wax tubules first appear associated with cork-silica cell complexes when internode cell elongation is complete. Sorghum homologs of genes involved in wax biosynthesis/transport were differentially expressed in the stem epidermis. Expression of many wax pathway genes (i.e., SbKCS6, SbCER3-1, SbWSD1, SbABCG12, SbABCG11) is low in immature apical internodes then increases at the onset of stem wax accumulation. SbCER4 is expressed relatively early in stem development consistent with accumulation of C28/30 primary alcohols on developing apical internodes. High expression of two SbCER3 homologs in fully elongated internodes is consistent with a role in production of C28/30 aldehydes. Gene regulatory network analysis aided the identification of sorghum homologs of transcription factors that regulate wax biosynthesis (i.e., SbSHN1, SbWRI1/3, SbMYB94/96/30/60, MYS1) and other transcription factors that could regulate and specify expression of the wax pathway in epidermal cells during cuticle development.

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.

Transcriptome and Metabolome Reveal Sugar and Organic Acid Accumulation in Rosa roxburghii Fruit

Sugars and organic acids significantly impact fruit sensory quality, but their accumulation patterns and regulatory mechanisms during the development of Rosa roxburghii fruit are still unclear. We utilized transcriptomics and metabolomics to investigate genes related to sugar and organic acid metabolism in Rosa roxburghii. Metabolomics data revealed that sucrose, glucose and fructose were the primary sugars, whereas citric acid and malic acid were the primary organic acids in Rosa roxburghii fruit. We constructed the metabolic pathways of major sugars and organic acids in Rosa roxburghii and identified five key genes involved in sugar and organic acid synthesis. In addition, we identified a module containing 132 transcription factors that was significantly associated with sucrose, citric acid and malic acid. Based on quantitative polymerase chain reaction (qPCR), we identified 13 transcription factors involved in sugar and organic acid metabolism, including the transcription factor RrANL2 and the sucrose synthase gene RrSUS3. Further yeast one-hybrid and dual luciferase assays showed that RrANL2 could bind to the promoter of RrSUS3 to increase its expression. These results provide new insights into the metabolism of sugars and organic acids in Rosa roxburghii fruit.

MYB-bHLH-TTG1 in a Multi-tiered Pathway Regulates Arabidopsis Seed Coat Mucilage Biosynthesis Genes Including PECTIN METHYLESTERASE INHIBITOR14 Required for Homogalacturonan Demethylesterification

MYB-bHLH-TTG1 (MBW) transcription factor (TF) complexes regulate Arabidopsis seed coat biosynthesis pathways via a multi-tiered regulatory mechanism. The MYB genes include MYB5, MYB23 and TRANSPARENT TESTA2 (TT2), which regulate GLABRA2 (GL2), HOMEODOMAIN GLABROUS2 (HDG2) and TRANSPARENT TESTA GLABRA2 (TTG2). Here, we examine the role of PECTIN METHYLESTERASE INHIBITOR14 (PMEI14) in seed coat mucilage pectin methylesterification and provide evidence in support of multi-tiered regulation of seed coat mucilage biosynthesis genes including PMEI14. The PMEI14 promoter was active in the seed coat and developing embryo. A pmei14 mutant exhibited stronger attachment of the outer layer of seed coat mucilage, increased mucilage homogalacturonan demethylesterification and reduced seed coat radial cell wall thickness, results consistent with decreased PMEI activity giving rise to increased PME activity. Reduced mucilage release from the seeds of myb5, myb23, tt2 and gl2, hdg2, ttg2 triple mutants indicated that HDG2 and MYB23 play minor roles in seed coat mucilage deposition. Chromatin immunoprecipitation analysis found that MYB5, TT8 and seven mucilage pathway structural genes are directly regulated by MYB5. Expression levels of GL2, HDG2, TTG2 and nine mucilage biosynthesis genes including PMEI14 in the combinatorial mutant seeds indicated that these genes are positively regulated by at least two of those six TFs and that TTG1 and TTG2 are major regulators of PMEI14 expression. Our results show that MYB-bHLH-TTG1 complexes regulate mucilage biosynthesis genes, including PMEI14, both directly and indirectly via a three-tiered mechanism involving GL2, HDG2 and TTG2.

Enhancer activation via TCP and HD-ZIP and repression by Dof transcription factors mediate giant cell-specific expression

Proper cell-type identity relies on highly coordinated regulation of gene expression. Regulatory elements such as enhancers can produce cell type-specific expression patterns, but the mechanisms underlying specificity are not well understood. We previously identified an enhancer region capable of driving specific expression in giant cells, which are large, highly endoreduplicated cells in the Arabidopsis thaliana sepal epidermis. In this study, we use the giant cell enhancer as a model to understand the regulatory logic that promotes cell type-specific expression. Our dissection of the enhancer revealed that giant cell specificity is mediated primarily through the combination of two activators and one repressor. HD-ZIP and TCP transcription factors are involved in the activation of expression throughout the epidermis. High expression of HD-ZIP transcription factor genes in giant cells promoted higher expression driven by the enhancer in giant cells. Dof transcription factors repressed the activity of the enhancer such that only giant cells maintained enhancer activity. Thus, our data are consistent with a conceptual model whereby cell type-specific expression emerges from the combined activities of three transcription factor families activating and repressing expression in epidermal cells.

Comprehensive analysis of OpHD-ZIP transcription factors related to the regulation of camptothecin biosynthesis in Ophiorrhiza pumila

Ophiorrhiza pumila, as a folk herb belonging to the Rubiaceae family, has become a potential source of camptothecin (CPT), which is a monoterpenoid indole alkaloid with good antitumor property. However, the camptothecin content in this herb is low, and is far from meeting the increasing clinical demand. Understanding the transcriptional regulation of camptothecin biosynthesis provides an effective strategy for improvement of camptothecin yield. Previous studies have demonstrated several transcription factors that are related to camptothecin biosynthesis, while the functions of HD-ZIP members in O. pumila have not been investigated yet. In this study, 32 OpHD-ZIP transcription factor members were genome-wide identified. Phylogenetic tree showed that these OpHD-ZIP proteins are divided into four subfamilies. Based on the transcriptome data, nine OpHD-ZIP genes were shown to be predominantly expressed in O. pumila roots, which were in line with the camptothecin biosynthetic genes. Co-expression analysis showed that OpHD-ZIP7 and OpHD-ZIP20 were potentially related to the modulation of camptothecin biosynthesis. Dual-luciferase reporter assays (Dual-LUC) showed that both OpHD-ZIP7 and OpHD-ZIP20 could activate the expression of camptothecin biosynthetic genes OpIO and OpTDC. In conclusion, this study offered the promising data for exploring the roles of OpHD-ZIP transcription factors in regulating camptothecin biosynthesis.

A conserved mechanism determines the activity of two pivotal transcription factors that control epidermal cell differentiation in Arabidopsis thaliana

The surface of plants is covered by the epidermis, which protects the plant's body from the external environment and mediates inter-cell layer signaling to regulate plant development. Therefore, the manifestation of epidermal traits at a precise location is a prerequisite for their normal growth and development. In Arabidopsis thaliana, class IV homeodomain-leucine zipper transcription factors PROTODERMAL FACTOR2 (PDF2) and ARABIDOPSIS THALIANA MERISTEM LAYER1 (ATML1) play redundant roles in epidermal cell differentiation. Nevertheless, several pieces of evidence suggest that the activity and/or function of PDF2 and ATML1 are regulated differently. The role of the steroidogenic acute regulatory protein-related lipid transfer (START) domain of ATML1 in restricting this protein's activity has been demonstrated; however, whether this lipid-dependent mechanism regulates PDF2 expression is unknown. In this study, we demonstrated that the START domains of PDF2 and ATML1, regulate protein turnover in a position-dependent manner and affect the dimeric proteins. Our results show that a conserved mechanism provides the basis for the functional redundancy of PDF2 and ATML1 in epidermal cell differentiation and that an unidentified regulatory layer specific to PDF2 or ATML1 is responsible for the difference in the activity and/or function of PDF2 and ATML1.

HD-ZIP Transcription Factors and Brassinosteroid Signaling Play a Role in Capitulum Patterning in Chrysanthemum

Chrysanthemum is a genus in the Asteraceae family containing numerous cut flower varieties with high ornamental value. It owes its beauty to the composite flower head, which resembles a compact inflorescence. This structure is also known as a capitulum, in which many ray and disc florets are densely packed. The ray florets are localized at the rim, are male sterile, and have large colorful petals. The centrally localized disc florets develop only a small petal tube but produce fertile stamens and a functional pistil. Nowadays, varieties with more ray florets are bred because of their high ornamental value, but, unfortunately, this is at the expense of their seed setting. In this study, we confirmed that the disc:ray floret ratio is highly correlated to seed set efficiency, and therefore, we further investigated the mechanisms that underlie the regulation of the disc:ray floret ratio. To this end, a comprehensive transcriptomics analysis was performed in two acquired mutants with a higher disc:ray floret ratio. Among the differentially regulated genes, various potential brassinosteroid (BR) signaling genes and HD-ZIP class IV homeodomain transcription factors stood out. Detailed follow-up functional studies confirmed that reduced BR levels and downregulation of HD-ZIP IV gene Chrysanthemum morifolium PROTODERMAL FACTOR 2 (CmPDF2) result in an increased disc:ray floret ratio, thereby providing ways to improve seed set in decorative chrysanthemum varieties in the future.

Asymmetric gene expression and cell-type-specific regulatory networks in the root of bread wheat revealed by single-cell multiomics analysis

BACKGROUND

Homoeologs are defined as homologous genes resulting from allopolyploidy. Bread wheat, Triticum aestivum, is an allohexaploid species with many homoeologs. Homoeolog expression bias, referring to the relative contribution of homoeologs to the transcriptome, is critical for determining the traits that influence wheat growth and development. Asymmetric transcription of homoeologs has been so far investigated in a tissue or organ-specific manner, which could be misleading due to a mixture of cell types.

RESULTS

Here, we perform single nuclei RNA sequencing and ATAC sequencing of wheat root to study the asymmetric gene transcription, reconstruct cell differentiation trajectories and cell-type-specific gene regulatory networks. We identify 22 cell types. We then reconstruct cell differentiation trajectories that suggest different origins between epidermis/cortex and endodermis, distinguishing bread wheat from Arabidopsis. We show that the ratio of asymmetrically transcribed triads varies greatly when analyzing at the single-cell level. Hub transcription factors determining cell type identity are also identified. In particular, we demonstrate that TaSPL14 participates in vasculature development by regulating the expression of BAM1. Combining single-cell transcription and chromatin accessibility data, we construct the pseudo-time regulatory network driving root hair differentiation. We find MYB3R4, REF6, HDG1, and GATAs as key regulators in this process.

CONCLUSIONS

Our findings reveal the transcriptional landscape of root organization and asymmetric gene transcription at single-cell resolution in polyploid wheat.

Transcription factors NtHD9 and NtHD12 control long glandular trichome formation via jasmonate signaling

Glandular trichomes are universal epidermal structures that produce abundant specialized metabolites. However, knowledge of the initiation of glandular heads in glandular trichomes is limited. Herein, we found an intrinsic link of morphogenesis between glandular trichomes and non-glandular trichomes. Two novel homeodomain leucine zipper II members in tobacco (Nicotiana tabacum), NtHD9 and NtHD12, played important roles in long glandular trichome formation: NtHD9 was responsible for glandular head formation, while NtHD12 simultaneously controlled the formation of stalks and glandular heads. DAP-seq analysis suggested that NtHD9 can bind to the KKGCATTWAWTR motif of the cytochromes P450 94C1 (NtCYP94C1) promoter, which is involved in jasmonoyl-isoleucine oxidation. RNA-seq analysis of non-transformed tobacco and nthd9 plants revealed that NtHD9 modulates the expression of jasmonate (JA) signaling- and six trichome development-related genes. Notably, MeJA treatment restored the morphogenesis of long glandular trichomes in nthd9 and nthd12 plants, and the size of glandular heads increased with increasing MeJA concentration. However, the phenotype of long glandular trichome absence in double mutants of NtHD9 and NtHD12 could not be restored by MeJA. Our data demonstrate that NtHD9 and NtHD12 have distinct major functions yet overlapping roles in long glandular trichome formation via JA signaling.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Genome-wide identification of HD-ZIP gene family and screening of genes related to prickle development in Zanthoxylum armatum

Zanthoxylum armatum is an important cash crop for medicinal and food purposes in Asia. However, its stems and leaves are covered with a large number of prickles, which cause many problems in the production process. The homeodomain leucine zipper (HD-ZIP) gene family is a class of transcription factors unique to plants that play an important role in biological processes such as morphogenesis, signal transduction, and secondary metabolite synthesis. However, little is known about HD-ZIP gene information that may be involved in prickle development of Z. armatum. Here, we identified 76 ZaHDZ genes from the Z. armatum genome and classified them into four subfamilies (I-IV) based on phylogenetic analysis, a classification further supported by gene structure and conserved motif analysis. Seventy-six ZaHDZ genes were unevenly distributed on chromosomes. Evolutionary analysis revealed that the expansion of ZaHDZ genes mainly were due to whole-genome duplication (WGD) or segmental duplication, and they experienced strong purifying selection pressure in the process of evolution. A total of 47 cis-elements were identified in the promoter region of ZaHDZ genes. Quantitative real-time polymerase chain reaction analysis was performed on subfamily IV ZaHDZ gene expression levels in five tissues and under four hormone treatments. Finally, ZaHDZ16 was predicted to be the candidate gene most likely to be involved in prickle development of Z. armatum. These results contribute to a better understanding of the characteristics of HD-ZIP gene family and lay a foundation for further study on the function of genes related to prickle development of Z. armatum.

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.

Genome-Wide Characterization and Expression Profiling of HD-Zip Genes in ABA-Mediated Processes in Fragaria vesca

Members of homeodomain-leucine zipper (HD-Zip) transcription factors can play their roles by modulating abscisic acid (ABA) signaling in Arabidopsis. So far, our knowledge of the functions of HD-Zips in woodland strawberries (Fragaria vesca), a model plant for studying ABA-mediated fruit ripening, is limited. Here, we identified a total of 31 HD-Zip genes (FveHDZ1-31) in F. vesca, and classified them into four subfamilies (I to IV). Promoter analyses show that the ABA-responsive element, ABRE, is prevalent in the promoters of subfamily I and II FveHDZs. RT-qPCR results demonstrate that 10 of the 14 investigated FveHDZs were consistently >1.5-fold up-regulated or down-regulated in expression in response to exogenous ABA, dehydration, and ABA-induced senescence in leaves. Five of the six consistently up-regulated genes are from subfamily I and II. Thereinto, FveHDZ4, and 20 also exhibited significantly enhanced expression along with increased ABA content during fruit ripening. In yeast one-hybrid assays, FveHDZ4 proteins could bind the promoter of an ABA signaling gene FvePP2C6. Collectively, our results strongly support that the FveHDZs, particularly those from subfamilies I and II, are involved in the ABA-mediated processes in F. vesca, providing a basis for further functional characterization of the HD-Zips in strawberry and other plants.

A novel regulatory complex mediated by Lanata (Ln) controls multicellular trichome formation in tomato

Trichomes that originate from plant aerial epidermis act as mechanical and chemical barriers against herbivores. Although several regulators have recently been identified, the regulatory pathway underlying multicellular trichome formation remains largely unknown in tomato. Here, we report a novel HD-ZIP IV transcription factor, Lanata (Ln), a missense mutation which caused the hairy phenotype. Biochemical analyses demonstrate that Ln separately interacts with two trichome regulators, Woolly (Wo) and Hair (H). Genetic and molecular evidence demonstrates that Ln directly regulates the expression of H. The interaction between Ln and Wo can increase trichome density by enhancing the expression of SlCycB2 and SlCycB3, which we previously showed are involved in tomato trichome formation. Furthermore, SlCycB2 represses the transactivation of the SlCycB3 gene by Ln and vice versa. Our findings provide new insights into the novel regulatory network controlling multicellular trichome formation in tomato.

Identification and functional validation of super-enhancers in Arabidopsis thaliana

Super-enhancers (SEs) are exceptionally large enhancers and are recognized to play prominent roles in cell identity in mammalian species. We surveyed the genomic regions containing large clusters of accessible chromatin regions (ACRs) marked by deoxyribonuclease (DNase) I hypersensitivity in Arabidopsis thaliana. We identified a set of 749 putative SEs, which have a minimum length of 1.5 kilobases and represent the top 2.5% of the largest ACR clusters. We demonstrate that the genomic regions associating with these SEs were more sensitive to DNase I than other nonpromoter ACRs. The SEs were preferentially associated with topologically associating domains. Furthermore, the SEs and their predicted cognate genes were frequently associated with organ development and tissue identity in A. thaliana. Therefore, the A. thaliana SEs and their cognate genes mirror the functional characteristics of those reported in mammalian species. We developed CRISPR/Cas-mediated deletion lines of a 3,578-bp SE associated with the thalianol biosynthetic gene cluster (BGC). Small deletions (131-157 bp) within the SE resulted in distinct phenotypic changes and transcriptional repression of all five thalianol genes. In addition, T-DNA insertions in the SE region resulted in transcriptional alteration of all five thalianol genes. Thus, this SE appears to play a central role in coordinating the operon-like expression pattern of the thalianol BGC.

The START domain mediates Arabidopsis GLABRA2 dimerization and turnover independently of homeodomain DNA binding

Class IV homeodomain leucine-zipper transcription factors (HD-Zip IV TFs) are key regulators of epidermal differentiation that are characterized by a DNA-binding HD in conjunction with a lipid-binding domain termed steroidogenic acute regulatory-related lipid transfer (START). Previous work established that the START domain of GLABRA2 (GL2), a HD-Zip IV member from Arabidopsis (Arabidopsis thaliana), is required for TF activity. Here, we addressed the functions and possible interactions of START and the HD in DNA binding, dimerization, and protein turnover. Deletion analysis of the HD and missense mutations of a conserved lysine (K146) resulted in phenotypic defects in leaf trichomes, root hairs, and seed mucilage, similar to those observed for START domain mutants, despite nuclear localization of the respective proteins. In vitro and in vivo experiments demonstrated that while HD mutations impair binding to target DNA, the START domain is dispensable for DNA binding. Vice versa, protein interaction assays revealed impaired GL2 dimerization for multiple alleles of START mutants, but not HD mutants. Using in vivo cycloheximide chase experiments, we provided evidence for the role of START, but not HD, in maintaining protein stability. This work advances our mechanistic understanding of HD-Zip TFs as multidomain regulators of epidermal development in plants.

Comparative transcriptome analysis reveals key pathways and genes involved in trichome development in tea plant (Camellia sinensis)

Trichomes, which develop from epidermal cells, are considered one of the important characteristics of the tea plant [Camellia sinensis (L.) O. Kuntze]. Many nutritional and metabolomic studies have indicated the important contributions of trichomes to tea products quality. However, understanding the regulation of trichome formation at the molecular level remains elusive in tea plants. Herein, we present a genome-wide comparative transcriptome analysis between the hairless Chuyeqi (CYQ) with fewer trichomes and the hairy Budiaomao (BDM) with more trichomes tea plant genotypes, toward the identification of biological processes and functional gene activities that occur during trichome development. In the present study, trichomes in both cultivars CYQ and BDM were unicellular, unbranched, straight, and soft-structured. The density of trichomes was the highest in the bud and tender leaf periods. Further, using the high-throughput sequencing method, we identified 48,856 unigenes, of which 31,574 were differentially expressed. In an analysis of 208 differentially expressed genes (DEGs) encoding transcription factors (TFs), five may involve in trichome development. In addition, on the basis of the Gene Ontology (GO) annotation and the weighted gene co-expression network analysis (WGCNA) results, we screened several DEGs that may contribute to trichome growth, including 66 DEGs related to plant resistance genes (PRGs), 172 DEGs related to cell wall biosynthesis pathway, 29 DEGs related to cell cycle pathway, and 45 DEGs related to cytoskeleton biosynthesis. Collectively, this study provided high-quality RNA-seq information to improve our understanding of the molecular regulatory mechanism of trichome development and lay a foundation for additional trichome studies in tea plants.

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.

Evolutionary relationship analysis of STARD gene family and VvSTARD5 improves tolerance of salt stress in transgenic tomatoes

The steroidogenic acute regulatory protein-related lipid transfer domain (STARD) forms a protein that can bind membrane-derived phospholipid second messengers and plasma membranes. Although it has been reported in many plants, the evolutionary relationship of the STARD gene family has not been systematically analyzed, and functions of the HD-START and HD-START-MEKHLA domain subgroup genes under hormone and abiotic stress are also unclear in grapes. This study identified and analyzed 23 VvSTARD genes, which were distinctly divided into five subgroups according to five conserved domain types. The analyses of codon preference, selective pressure, and synteny relationship revealed that grape had higher homology with Arabidopsis compared with rice. Interestingly, the expression levels of VvSTARD genes in subgroups 1, 2, and 3 exhibited significant upregulation under NaCl treatment at 24 h, but VvSTARD genes in subgroups 4 and 5 were upregulated under methyl jasmonate (MeJA) treatment at 24 h. The subcellular localization showed that VvSTARD5 was localized in the nucleus. Additionally, under NaCl treatment at 24 h, there were an obvious decrease in the relative electrical leakages and the content of malondialdehyde (MDA), while the relative expression level of VvSTARD5 and content of proline were obviously enhanced in three transgenic lines. Therefore, the overexpression of VvSTARD5 greatly increased the salt tolerance of transgenic tomatoes. Collectively, this study preliminarily explores the comprehensive function of the STARD gene family in grapes and verifies the function of VvSTARD5 in response to salt.

CRISPR/Cas9 mediated gene-editing of GmHdz4 transcription factor enhances drought tolerance in soybean (Glycine max [L.] Merr.)

The HD-Zip transcription factors play a crucial role in plant development, secondary metabolism, and abiotic stress responses, but little is known about HD-Zip I genes in soybean. Here, a homeodomain-leucine zipper gene designated GmHdz4 was isolated. Chimeric soybean plants, GmHdz4 overexpressing (GmHdz4-oe), and gene-editing via CRISPR/Cas9 (gmhdz4) in hairy roots, were generated to examine the GmHdz4 gene response to polyethylene glycol (PEG)-simulated drought stress. Bioinformatic analysis showed GmHdz4 belonged to clade δ, and was closely related to other drought tolerance-related HD-Zip I family genes such as AtHB12, Oshox12, and Gshdz4. The GmHdz4 was located in the plant nucleus and showed transcriptional activation activity by yeast hybrid assay. Quantitative real-time PCR analysis revealed that GmHdz4 expression varied in tissues and was induced by PEG-simulated drought stress. The gmhdz4 showed promoted growth of aboveground parts, and its root system architecture, including the total root length, the root superficial area, and the number of root tips were significantly higher than those of GmHdz4-oe even the non-transgenic line (NT) on root tips number. The better maintenance of turgor pressure by osmolyte accumulation, and the higher activity of antioxidant enzymes to scavenge reactive oxygen species, ultimately suppressed the accumulation of hydrogen peroxide (H₂O₂), superoxide anion (O^2-), and malondialdehyde (MDA), conferring higher drought tolerance in gmhdz4 compared with both GmHdz4-oe and NT. Together, our results provide new insights for future research on the mechanisms by which GmHdz4 gene-editing via CRISPR/Cas9 system could promote drought stress and provide a potential target for molecular breeding in soybean.

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.

Identification and analysis of HD-Zip genes involved in the leaf development of Liriodendron chinense using multidimensional analysis

Homeodomain-leucine zipper (HD-Zip) proteins are plant-specific transcription factors that play important roles in different biological processes, especially leaf development. However, no studies to date have identified the HD-Zip genes in Liriodendron chinense nor characterized their functions. We identified the HD-Zip genes in L. chinense by analysing the phylogeny, chromosome location, structure, conserved motif, cis-regulatory elements, synteny, post-transcriptional regulation and expression patterns of these genes during leaf development. A total of 36 LcHD-Zip genes were identified and divided into four subfamilies (HD-Zip I to IV). Synteny analysis revealed that segmental duplication was the main force driving the expansion of LcHD-Zip genes. These 36 LcHD-Zip genes exhibited 11 different expression patterns. Pattern 1, 2, 3, 4, 6, 7, 8 and 9 genes may play important roles in leaf development, such as leaf initiation, leaf polarity establishment, leaf shape development, phytohormone-mediated leaf growth and leaf epidermal structure formation. Four HD-Zip III genes were targeted by microRNAs (miRNAs), and the miR165/166a-HD-Zip regulatory module formed regulated leaf initiation and leaf polarity establishment. Overall, LcHD-Zip genes play key roles in leaf development of L. chinense. This work provides a foundation for the functional verification of HD-Zip genes identified in this study.

A Chromosome Level Genome Assembly of a Winter Turnip Rape (Brassica rapa L.) to Explore the Genetic Basis of Cold Tolerance

Winter rapeseed (Brassica rapa L.) is an important overwintering oilseed crop that is widely planted in northwest China and suffers chronic low temperatures in winter. So the cold stress becomes one of the major constraints that limit its production. The currently existing genomes limit the understanding of the cold-tolerant genetic basis of rapeseed. Here we assembled a high-quality long-read genome of B. rapa "Longyou-7" cultivar, which has a cold-tolerant phenotype, and constructed a graph-based pan-genome to detect the structural variations within homologs of currently reported cold-tolerant related genes in the "Longyou-7" genome, which provides an additional elucidation of the cold-tolerant genetic basis of "Longyou-7" cultivar and promotes the development of cold-tolerant breeding in B. rapa.

CROWN ROOTLESS1 binds DNA with a relaxed specificity and activates OsROP and OsbHLH044 genes involved in crown root formation in rice

In cereals, the root system is mainly composed of post-embryonic shoot-borne roots, named crown roots. The CROWN ROOTLESS1 (CRL1) transcription factor, belonging to the ASYMMETRIC LEAVES2-LIKE/LATERAL ORGAN BOUNDARIES DOMAIN (ASL/LBD) family, is a key regulator of crown root initiation in rice (Oryza sativa). Here, we show that CRL1 can bind, both in vitro and in vivo, not only the LBD-box, a DNA sequence recognized by several ASL/LBD transcription factors, but also another not previously identified DNA motif that was named CRL1-box. Using rice protoplast transient transactivation assays and a set of previously identified CRL1-regulated genes, we confirm that CRL1 transactivates these genes if they possess at least a CRL1-box or an LBD-box in their promoters. In planta, ChIP-qPCR experiments targeting two of these genes that include both a CRL1- and an LBD-box in their promoter show that CRL1 binds preferentially to the LBD-box in these promoter contexts. CRISPR/Cas9-targeted mutation of these two CRL1-regulated genes, which encode a plant Rho GTPase (OsROP) and a basic helix-loop-helix transcription factor (OsbHLH044), show that both promote crown root development. Finally, we show that OsbHLH044 represses a regulatory module, uncovering how CRL1 regulates specific processes during crown root formation.

Genome-Wide Identification and Expression Analysis of Homeodomain Leucine Zipper Subfamily IV (HD-ZIP IV) Gene Family in Cannabis sativa L

The plant-specific homeodomain zipper family (HD-ZIP) of transcription factors plays central roles in regulating plant development and environmental resistance. HD-ZIP transcription factors IV (HDZ IV) have been involved primarily in the regulation of epidermal structure development, such as stomata and trichomes. In our study, we identified nine HDZ IV-encoding genes in Cannabis sativa L. by conducting a computational analysis of cannabis genome resources. Our analysis suggests that these genes putatively encode proteins that have all the conserved domains of HDZ IV transcription factors. The phylogenetic analysis of HDZ IV gene family members of cannabis, rice (Oryza sativa), and Arabidopsis further implies that they might have followed distinct evolutionary paths after divergence from a common ancestor. All the identified cannabis HDZ IV gene promoter sequences have multiple regulation motifs, such as light- and hormone-responsive elements. Furthermore, experimental evidence shows that different HDZ IV genes have different expression patterns in root, stem, leaf, and flower tissues. Four genes were primarily expressed in flowers, and the expression of CsHDG5 (XP_030501222.1) was also correlated with flower maturity. Fifty-nine genes were predicted as targets of HDZ IV transcription factors. Some of these genes play central roles in pathogen response, flower development, and brassinosteroid signaling. A subcellular localization assay indicated that one gene of this family is localized in the Arabidopsis protoplast nucleus. Taken together, our work lays fundamental groundwork to illuminate the function of cannabis HDZ IV genes and their possible future uses in increasing cannabis trichome morphogenesis and secondary metabolite production.

Multi-Dimensional Molecular Regulation of Trichome Development in Arabidopsis and Cotton

Plant trichomes are specialized epidermal cells that are widely distributed on plant aerial tissues. The initiation and progression of trichomes are controlled in a coordinated sequence of multiple molecular events. During the past decade, major breakthroughs in the molecular understanding of trichome development were achieved through the characterization of various trichomes defective mutants and trichome-associated genes, which revealed a highly complex molecular regulatory network underlying plant trichome development. This review focuses on the recent millstone in plant trichomes research obtained using genetic and molecular studies, as well as 'omics' analyses in model plant Arabidopsis and fiber crop cotton. In particular, we discuss the latest understanding and insights into the underlying molecular mechanisms of trichomes formation at multiple dimensions, including at the chromatin, transcriptional, post-transcriptional, and post-translational levels. We summarize that the integration of multi-dimensional trichome-associated genes will enable us to systematically understand the molecular regulation network that landscapes the development of the plant trichomes. These advances will enable us to address the unresolved questions regarding the molecular crosstalk that coordinate concurrent and ordered the changes in cotton fiber initiation and progression, together with their possible implications for genetic improvement of cotton fiber.

Two R2R3-MYB genes cooperatively control trichome development and cuticular wax biosynthesis in Prunus persica

The fruit surface has an enormous impact on the external appearance and postharvest shelf-life of fruit. Here, we report two functionally redundant genes, PpMYB25 and PpMYB26, involved in regulation of fruit skin texture in peach. PpMYB25 can activate transcription of PpMYB26 and they both induce trichome development and cuticular wax accumulation, resulting in peach fruit with a fuzzy and dull appearance. By contrast, nonfunctional mutation of PpMYB25 caused by an insertional retrotransposon in the last exon in nectarine fails to activate transcription of PpMYB26, resulting in nectarine fruit with a smooth and shiny appearance due to loss of trichome initiation and decreased cuticular wax accumulation. Secondary cell wall biosynthesis in peach fruit pubescence is controlled by a transcriptional regulatory network, including the master regulator PpNAC43 and its downstream MYB transcription factors such as PpMYB42, PpMYB46 and PpMYB83. Our results show that PpMYB25 and PpMYB26 coordinately regulate fruit pubescence and cuticular wax accumulation and their simultaneous perturbation results in the origin of nectarine, which is botanically classified as a subspecies of peach.

Function and Regulation of microRNA171 in Plant Stem Cell Homeostasis and Developmental Programing

MicroRNA171 (miR171), a group of 21-nucleotide single-strand small RNAs, is one ancient and conserved microRNA family in land plants. This review focuses on the recent progress in understanding the role of miR171 in plant stem cell homeostasis and developmental patterning, and the regulation of miR171 by developmental cues and environmental signals. Specifically, miR171 regulates shoot meristem activity and phase transition through repressing the HAIRYMERISTEM (HAM) family genes. In the model species Arabidopsis, miR171 serves as a short-range mobile signal, which initiates in the epidermal layer of shoot meristems and moves downwards within a limited distance, to pattern the apical-basal polarity of gene expression and drive stem cell dynamics. miR171 levels are regulated by light and various abiotic stresses, suggesting miR171 may serve as a linkage between environmental factors and cell fate decisions. Furthermore, miR171 family members also demonstrate both conserved and lineage-specific functions in land plants, which are summarized and discussed here.

A Survey of Enhanced Cold Tolerance and Low-Temperature-Induced Anthocyanin Accumulation in a Novel Zoysia japonica Biotype

Zoysia japonica is a warm-season turfgrass that is extensively used in landscaping, sports fields, and golf courses worldwide. Uncovering the low-temperature response mechanism of Z. japonica can help to accelerate the development of new cold-tolerant cultivars, which could be used to prolong the ornamental and usage duration of turf. A novel Z. japonica biotype, YueNong-9 (YN-9), was collected from northeastern China for this study. Phenotypic measurements, cold-tolerance investigation, and whole-transcriptome surveys were performed on YN-9 and LanYin-3 (LY-3), the most popular Z. japonica cultivar in Southern China. The results indicated the following: YN-9 has longer second and third leaves than LY-3; when exposed to the natural low temperature during winter in Guangzhou, YN-9 accumulated 4.74 times more anthocyanin than LY-3; after cold acclimation and freezing treatment, 83.25 ± 9.55% of YN-9 survived while all LY-3 leaves died, and the dark green color index (DGCI) value of YN-9 was 1.78 times that of LY-3; in YN-9, there was a unique up-regulation of Phenylalanine ammonia-lyase (PAL), Homeobox-leucine Zipper IV (HD-ZIP), and ATP-Binding Cassette transporter B8 (ABCB8) expressions, as well as a unique down-regulation of zinc-regulated transporters and iron-regulated transporter-like proteins (ZIPs) expression, which may promote anthocyanin biosynthesis, transport, and accumulation. In conclusion, YN-9 exhibited enhanced cold tolerance and is thus an excellent candidate for breeding cold-tolerant Z. japonica variety, and its unique low-temperature-induced anthocyanin accumulation and gene responses provide ideas and candidate genes for the study of low-temperature tolerance mechanisms and genetic engineering breeding.

Developmental regulation of lupulin gland-associated genes in aromatic and bitter hops (Humulus lupulus L.)

BACKGROUND

Hop (Humulus lupulus L.) bitter acids are valuable metabolites for the brewing industry. They are biosynthesized and accumulate in glandular trichomes of the female inflorescence (hop cone). The content of alpha bitter acids, such as humulones, in hop cones can differentiate aromatic from bitter hop cultivars. These contents are subject to genetic and environmental control but significantly correlate with the number and size of glandular trichomes (lupulin glands).

RESULTS

We evaluated the expression levels of 37 genes involved in bitter acid biosynthesis and morphological and developmental differentiation of glandular trichomes to identify key regulatory factors involved in bitter acid content differences. For bitter acid biosynthesis genes, upregulation of humulone synthase genes, which are important for the biosynthesis of alpha bitter acids in lupulin glands, could explain the higher accumulation of alpha bitter acids in bitter hops. Several transcription factors, including HlETC1, HlMYB61 and HlMYB5 from the MYB family, as well as HlGLABRA2, HlCYCB2-4, HlZFP8 and HlYABBY1, were also more highly expressed in the bitter hop cultivars; therefore, these factors may be important for the higher density of lupulin glands also seen in the bitter hop cultivars.

CONCLUSIONS

Gene expression analyses enabled us to investigate the differences between aromatic and bitter hops. This study confirmed that the bitter acid content in glandular trichomes (lupulin glands) is dependent on the last step of alpha bitter acid biosynthesis and glandular trichome density.

Inhibition of Arabidopsis stomatal development by plastoquinone oxidation

Stomata are the pores in the epidermal surface of plant leaves that regulate the exchange of water and CO₂ with the environment thus controlling leaf gas exchange.¹ In the model dicot plant Arabidopsis thaliana, the transcription factors SPEECHLESS (SPCH) and MUTE sequentially control formative divisions in the stomatal lineage by forming heterodimers with ICE1.² SPCH regulates entry into the stomatal lineage and its stability or activity is regulated by a mitogen-activated protein kinase (MAPK) signaling cascade, mediated by its interaction with ICE1.^3-6 This MAPK pathway is regulated by extracellular epidermal patterning factor (EPFs) peptides, which bind a transmembrane receptor complex to inhibit (EPF1 and EPF2) or promote (STOMAGEN/EPFL9) stomatal development.^7-9 MUTE controls the transition to guard mother cell identity and is regulated by the HD-ZIP transcription factor HDG2, which is expressed exclusively in stomatal lineage cells.^10,11 Light signals acting through phytochrome and cryptochrome photoreceptors positively regulate stomatal development in response to increased irradiance.^12,13 Here we report that stomatal development is also regulated by the redox state of the photosynthetic electron transport chain (PETC). Oxidation of the plastoquinone (PQ) pool inhibits stomatal development by negatively regulating SPCH and MUTE expression. This mechanism is dependent on MPK6 and forms part of the response to lowering irradiance, which is distinct to the photoreceptor dependent response to increasing irradiance. Our results show that environmental signals can act through the PETC, demonstrating that photosynthetic signals regulate the development of the pores through which CO2 enters the leaf.

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.

The lipid-binding START domain regulates the dimerization of ATML1 via modulating the ZIP motif activity in Arabidopsis thaliana

In Arabidopsis thaliana, the epidermis is the outermost cell layer composed of many specialized types of epidermal cells, such as pavement cells, trichomes, and guard cells. The homeodomain-leucine zipper (HD-ZIP) class Ⅳ transcription factors (TFs), which are unique to the plant kingdom, have been recognized as key regulators of epidermis development. Unlike animal HD proteins, which can bind to DNA as monomers, plant HD-ZIP class Ⅳ TFs bind to DNA as dimers, although little is known about the regulation of their dimerization process. Here, we show that the homodimerization of ARABIDOPSIS THALIANA MERISTEM LAYER 1 (ATML1) - HD-ZIP class Ⅳ TF that is required for protoderm development - is regulated by the lipid-binding steroidogenic acute regulatory protein-related lipid transfer (START) domain. We found that ATML1 forms homodimer through interaction via its ZIP motif in yeast and plant cells, although the interaction is abolished by generating a mutation into the lipid-binding START domain to disrupt the lipid-binding ability. These results suggest that lipidic ligands function as key regulators of protoderm development via modulating the dimerization of ATML1.

Cytosine methylation of the FWA promoter promotes direct in vitro shoot regeneration in Arabidopsis thaliana

Epigenetic modifications within promoter sequences can act as regulators of gene expression. Shoot regeneration is influenced by both DNA methylation and histone methylation, but the mechanistic basis of this regulation is obscure. Here, we identified 218 genes related to the regeneration capacity of callus that were differentially transcribed between regenerable calli (RC) and non-regenerable calli (NRC) in Arabidopsis thaliana. An analysis of the promoters of five of the differentially expressed genes (FWA, ACC1, TFL1, MAX3, and GRP3) pointed to an inverse relationship between cytosine methylation and transcription. The FWA promoter was demethylated and highly expressed in NRC, whereas it was methylated and expressed at low levels in RC. Explants of the hypomethylation mutants fwa-1 and fwa-2 showed strong levels of FWA expression and regenerated less readily than the wild type, suggesting that FWA inhibits direct in vitro shoot regeneration. WUSCHEL-RELATED HOMEOBOX 9 (WOX9), which is required for shoot apical meristem formation, was directly repressed by FWA. Overexpressing WOX9 partly rescued the shoot regeneration defect of fwa-2 plants. These findings suggest that cytosine methylation of the FWA promoter forms part of the regulatory system governing callus regenerability and direct in vitro shoot regeneration.

Genomic editing of intronic enhancers unveils their role in fine-tuning tissue-specific gene expression in Arabidopsis thaliana

Enhancers located in introns are abundant and play a major role in the regulation of gene expression in mammalian species. By contrast, the functions of intronic enhancers in plants have largely been unexplored and only a handful of plant intronic enhancers have been reported. We performed a genome-wide prediction of intronic enhancers in Arabidopsis thaliana using open chromatin signatures based on DNase I sequencing. We identified 941 candidate intronic enhancers associated with 806 genes in seedling tissue and 1,271 intronic enhancers associated with 1,069 genes in floral tissue. We validated the function of 15 of 21 (71%) of the predicted intronic enhancers in transgenic assays using a reporter gene. We also created deletion lines of three intronic enhancers associated with two different genes using CRISPR/Cas. Deletion of these enhancers, which span key transcription factor binding sites, did not abolish gene expression but caused varying levels of transcriptional repression of their cognate genes. Remarkably, the transcriptional repression of the deletion lines occurred at specific developmental stages and resulted in distinct phenotypic effects on plant morphology and development. Clearly, these three intronic enhancers are important in fine-tuning tissue- and development-specific expression of their cognate genes.

Proteomics of Homeobox7 Enhanced Salt Tolerance in Mesembryanthemum crystallinum

Mesembryanthemum crystallinum (common ice plant) is a halophyte species that has adapted to extreme conditions. In this study, we cloned a McHB7 transcription factor gene from the ice plant. The expression of McHB7 was significantly induced by 500 mM NaCl and it reached the peak under salt treatment for 7 days. The McHB7 protein was targeted to the nucleus. McHB7-overexpressing in ice plant leaves through Agrobacterium-mediated transformation led to 25 times more McHB7 transcripts than the non-transformed wild type (WT). After 500 mM NaCl treatment for 7 days, the activities of superoxide dismutase (SOD) and peroxidase (POD) and water content of the transgenic plants were higher than the WT, while malondialdehyde (MDA) was decreased in the transgenic plants. A total of 1082 and 1072 proteins were profiled by proteomics under control and salt treatment, respectively, with 22 and 11 proteins uniquely identified under control and salt stress, respectively. Among the 11 proteins, 7 were increased and 4 were decreased after salt treatment. Most of the proteins whose expression increased in the McHB7 overexpression (OE) ice plants under high salinity were involved in transport regulation, catalytic activities, biosynthesis of secondary metabolites, and response to stimulus. The results demonstrate that the McHB7 transcription factor plays a positive role in improving plant salt tolerance.