Arabidopsis thaliana HomeoBox 1 (AtHB1), a Homedomain-Leucine Zipper I (HD-Zip I) transcription factor, is regulated by PHYTOCHROME-INTERACTING FACTOR 1 to promote hypocotyl elongation

Cbuhdz34, a Homeodomain Leucine Zipper Transcription Factor, Positively Regulates Tension Wood Formation and Xylem Fibre Cell Elongation in Catalpa bungei

Catalpa bungei is a highly valued timber species renowned for its superior wood properties. However, the development of tension wood (TW) induced by wind and other mechanical stresses during the growing season significantly reduces its economic value. Although Homeodomain Leucine Zipper (HD-Zip), a plant-specific transcription factor family, has been reported to play various roles in plant growth, development, and stress resistance, a systematic characterisation of the HD-Zip gene family in C. bungei, particularly regarding the regulatory mechanisms involved in TW formation, is still lacking. Here, we identified a total of 48 HD-Zip genes (Cbuhdzs) in C. bungei and analysed their phylogeny, structure, and expression profiles. In particular, Cbuhdz34, a member of the HD-Zip I subfamily, was specifically upregulated during TW formation. To further explore its function, we overexpressed Cbuhdz34 (OE-Cbuhdz34) in poplar '84 K', which led to noticeable changes in plant growth and fibre cell length. Moreover, compared with wild-type plants, the OE-Cbuhdz34 plants presented increased TW formation under bending stress, as indicated by increased TW width, gelatinous layer width, and eccentric growth rate, suggesting a positive regulatory role in TW formation. Additionally, hierarchical genetic regulatory network analysis revealed the direct targets of Cbuhdz34, including CbuMYB63 and three genes involved in cell wall synthesis (CbuGATL1, CbuFLA17, and CbuLRR14). Further, yeast one-hybrid and dual-luciferase reporter assays confirmed the activation of these targets by Cbuhdz34. In conclusion, our results provide insights into the molecular mechanisms by which Cbuhdz34 regulates TW formation and lay a genetic foundation for the potential improvement of wood quality in C. bungei.

Cell Fate Determination of the Potato Shoot Apex and Stolon Tips Revealed by Single-Cell Transcriptome Analysis

Potato (Solanum tuberosum L.) is a starch-rich crop with two types of meristematic stems: the shoot and stolon. Shoots grow vertically, while stolons grow horizontally underground and produce tubers at their tips. However, transcriptional differences between shoot and stolon cells remain unclear. To address this, we performed single-cell RNA sequencing of the shoot apex and stolon tip, generating a comprehensive transcriptional landscape. We identified 23 distinct cell clusters with high cell heterogeneity, including cell-specific genes and conserved genes with cell-specific expression patterns. Hormone-related genes, particularly those involved in auxin and gibberellin pathways, exhibited distinct patterns among shoot and stolon cells. Meristematic cells were re-clustered based on the expression of StPOTH15, a homolog of SHOOT MERISTEMLESS (STM) in Arabidopsis. Co-expression networks of transcription factors identified the key transcription factors involved in stolon development. We also constructed developmental trajectories for xylem and phloem development using key vascular genes, including MP, XCP1, PP2A1 and SEOR1. Comparative analysis with Arabidopsis highlighted significant differences in cell type-specific transcript profiles. These results provide insights into the transcriptional divergence between potato shoot and stolon, and identify key transcription factors co-expressed with StPOTH15 that can be used to explore their roles in stolon development.

Two antagonistic gene regulatory networks drive Arabidopsis root hair growth at low temperature linked to a low-nutrient environment

Root hair (RH) cells can elongate to several hundred times their initial size, and are an ideal model system for investigating cell size control. Their development is influenced by both endogenous and external signals, which are combined to form an integrative response. Surprisingly, a low-temperature condition of 10°C causes increased RH growth in Arabidopsis and in several monocots, even when the development of the rest of the plant is halted. Previously, we demonstrated a strong correlation between RH growth response and a significant decrease in nutrient availability in the growth medium under low-temperature conditions. However, the molecular basis responsible for receiving and transmitting signals related to the availability of nutrients in the soil, and their relation to plant development, remain largely unknown. We have discovered two antagonic gene regulatory networks (GRNs) controlling RH early transcriptome responses to low temperature. One GNR enhances RH growth and it is commanded by the transcription factors (TFs) ROOT HAIR DEFECTIVE 6 (RHD6), HAIR DEFECTIVE 6-LIKE 2 and 4 (RSL2-RSL4) and a member of the homeodomain leucine zipper (HD-Zip I) group I 16 (AtHB16). On the other hand, a second GRN was identified as a negative regulator of RH growth at low temperature and it is composed by the trihelix TF GT2-LIKE1 (GTL1) and the associated DF1, a previously unidentified MYB-like TF (AT2G01060) and several members of HD-Zip I group (AtHB3, AtHB13, AtHB20, AtHB23). Functional analysis of both GRNs highlights a complex regulation of RH growth response to low temperature, and more importantly, these discoveries enhance our comprehension of how plants synchronize RH growth in response to variations in temperature at the cellular level.

Genome-wide identification of the bHLH gene family in Scutellaria baicalensis and their relationship with baicalin biosynthesis under drought stress

The bHLH gene family plays a critical role in regulating internal responses in plants. Although the pharmacological properties of Scutellaria baicalensis have been extensively studied, its bHLH gene family remains poorly investigated. In this study, 142 SbbHLH genes were identified using the complete genome data of S. baicalensis. Phylogenetic and conserved motif analyses were performed. Gene duplication events were analyzed, and cis-element analysis was conducted to explore regulatory factors. The expression patterns of these genes in different tissues and under drought stress were investigated using transcriptome data and qRT-PCR analysis. Phylogenetic and conserved motif analyses revealed that the gene structures within each SbbHLH clade are relatively conserved. Gene duplication analysis identified 29 duplication events in the SbbHLH gene family, most of which involved gene pairs under purifying selection. Cis-element analysis revealed that these genes are regulated by various environmental and hormonal factors. Transcriptomic data and qRT-PCR results demonstrated tissue-specific expression patterns for the 142 SbbHLH genes. Additionally, bHLH genes potentially involved in baicalin biosynthesis were identified under drought stress. The findings suggest that under drought stress, SbbHLH74, SbbHLH98, and SbbHLH142 are regulated by a network centered on SbbHLH53, which enhances baicalin biosynthesis. In conclusion, this study provides a comprehensive analysis of the bHLH gene family in S. baicalensis and identifies 4 potential SbbHLH genes involved in regulating baicalin biosynthesis under drought stress.

Passion fruit HD-ZIP genes: Characterization, expression variance, and overexpression PeHB31 enhanced drought tolerance via lignin pathway

The HD-ZIP (homeodomain-leucine zipper) genes hold significant importance in transcriptional regulation, especially in plant development and responses to abiotic stresses. However, a comprehensive study targeting HD-ZIP family members in passion fruit has been absent. In our current research, 34 HD-ZIP family members (PeHBs) were identified by bioinformatics analysis. Transcriptome analysis revealed that PeHBs exhibited distinct expression patterns when subjected to the four different abiotic stresses, and significant differential expression of PeHBs was also found among the three developmental stages of the fruit and between the purple and yellow genotype passion fruit leaves. An integrated metabolome and transcriptome analysis further revealed that the HD-ZIP III class gene PeHB31 (homologous to ATHB8), was co-upexpressed with lignans in yellow fruit P. edulis (commonly used as a resistance rootstock) when compared to purple fruit P. edulis. The transformation of Arabidopsis and yeast with the PeHB31 gene showed an enhancement in their capacity to withstand drought conditions. Notably, the transgenic Arabidopsis plants exhibited an increase in lignin content within the vascular tissues of their stems. This research lays the groundwork for future studies on the control mechanisms of lignin biosynthesis by HD-ZIP genes (especially HD-ZIP classes III and I) involved in drought tolerance.

Functional and regulatory diversity of homeobox-leucine zipper transcription factors BnaHB6 under dehydration and salt stress in Brassica napus L

The plant-specific homeodomain-leucine zipper I subfamily is involved in the regulation of various biological processes, particularly growth, development and stress response. In the present study, we characterized four BnaHB6 homologues from Brassica napus. All BnaHB6 proteins have transcriptional activation activity. Structural and functional data indicate the complex role of BnaHB6 genes in regulating biological processes, with some functions conserved and others diverged. Transcriptional analyzes revealed that they are induced in a similar manner in different tissues but show different expression patterns in response to stress and circadian rhythm. Only the BnaA09HB6 and BnaC08HB6 genes are expressed under dehydration and salt stress, and in darkness. The partial transcriptional overlap of BnaHB6s with the evolutionarily related genes BnaHB5 and BnaHB16 was also observed. Transgenic Arabidopsis thaliana plants expressing a single proBnaHB6::GUS partially confirmed the expression results. Bioinformatic analysis allowed the identification of TF-binding sites in the BnaHB6 promoters that may control their expression under stress and circadian rhythm. ChIP-qPCR analysis revealed that BnaA09HB6 and BnaC08HB6 bind directly to the promoters of the target genes BnaABF4 and BnaDREB2A. Comparison of their expression patterns in the WT plants and the bnac08hb6 mutant showed that BnaC08HB6 positively regulates the expression of the BnaABF4 and BnaDREB2A genes under dehydration and salt stress. We conclude that four BnaHB6 homologues have distinct functions in response to stress despite high sequence similarity, possibly indicating different binding preferences with BnaABF4 and BnaDREB2A. We hypothesize that BnaC08HB6 and BnaA09HB6 function in a complex regulatory network under stress.

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.

Molecular and biochemical responses of sesame (Sesame indicum L.) to rhizobacteria inoculation under water deficit

Introduction

Water scarcity is a challenge for sesame cultivation under rainfed conditions. In this scenario, a potential strategy to alleviate the water deficit is the application of plant growth-promoting bacteria. The objective of this study was to analyze the interaction of rhizobacteria with sesame cultivation under water deficit conditions.

Methods

An experiment was conducted in pots in a greenhouse using the BRS Morena sesame cultivar. The experimental design was completely randomized in a factorial scheme: 2 (irrigation regimes - daily irrigation and water deficit by suspending irrigation until 90% stomatal closure) x 6 (treatments with nitrogen or inoculants), with 5 replications. The types of fertilization were characterized by the addition of nitrogen (ammonium sulfate; 21% N), inoculants based on Bacillus spp. (pant001, ESA 13, and ESA 402), Agrobacterium sp. (ESA 441), and without nitrogen (control). On the fifth day after the suspension of irrigation, plant material was collected for gene expression analysis (DREB1 and HDZ7), activities of antioxidant enzymes (superoxide dismutase and catalase), relative proline content, and photosynthetic pigments. At the end of the crop cycle (about 85 days), production characteristics (root dry matter, aboveground dry matter, number of capsules, and thousand seed weight), as well as leaf nitrogen (N) and phosphorus (P) content, were evaluated.

Results and Discussion

There was a positive effect on both production and biochemical characteristics (proline, superoxide dismutase, catalase, and photosynthetic pigments). Regarding gene expression, most of the inoculated treatments exhibited increased expression of the DREB1 and HDZ7 genes. These biological indicators demonstrate the potential of rhizobacteria for application in sesame cultivation, providing nutritional supply and reducing the effects of water deficit.

Multi-locus genome-wide association study and genomic prediction for flowering time in chrysanthemum

MAIN CONCLUSION

Multi-locus GWAS detected several known and candidate genes responsible for flowering time in chrysanthemum. The associations could greatly increase the predictive ability of genome selection that accelerates the possible application of GS in chrysanthemum breeding. Timely flowering is critical for successful reproduction and determines the economic value for ornamental plants. To investigate the genetic architecture of flowering time in chrysanthemum, a multi-locus genome-wide association study (GWAS) was performed using a collection of 200 accessions and 330,710 single-nucleotide polymorphisms (SNPs) via 3VmrMLM method. Five flowering time traits including budding (FBD), visible colouring (VC), early opening (EO), full-bloom (OF) and senescing (SF) stages, plus five derived conditional traits were recorded in two environments. Extensive phenotypic variations were observed for these flowering time traits with coefficients of variation ranging from 6.42 to 38.27%, and their broad-sense heritability ranged from 71.47 to 96.78%. GWAS revealed 88 stable quantitative trait nucleotides (QTNs) and 93 QTN-by-environment interactions (QEIs) associated with flowering time traits, accounting for 0.50-8.01% and 0.30-10.42% of the phenotypic variation, respectively. Amongst the genes around these stable QTNs and QEIs, 21 and 10 were homologous to known flowering genes in Arabidopsis; 20 and 11 candidate genes were mined by combining the functional annotation and transcriptomics data, respectively, such as MYB55, FRIGIDA-like, WRKY75 and ANT. Furthermore, genomic selection (GS) was assessed using three models and seven unique marker datasets. We found the prediction accuracy (PA) using significant SNPs identified by GWAS under SVM model exhibited the best performance with PA ranging from 0.90 to 0.95. Our findings provide new insights into the dynamic genetic architecture of flowering time and the identified significant SNPs and candidate genes will accelerate the future molecular improvement of chrysanthemum.

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

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

Genome-Wide Identification and Expression Analysis of the bHLH Transcription Factor Family in Wintersweet (Chimonanthus praecox)

Wintersweet (Chimonanthus praecox (L.) Link, Calycanthaceae) is an esteemed ornamental flowering shrub known for its distinct blooming period in winter, vibrant color petals, and captivating floral fragrance. Basic helix-loop-helix (bHLH) transcription factors (TFs) play pivotal roles as key regulators in secondary metabolites biosynthesis, growth, and development in plants. However, the systematic analysis of the bHLH family members and their role in the regulation of floral traits in Wintersweet remains insufficiently understood. To bridge this knowledge gap, we conducted a comprehensive genome-wide analysis of the C. praecox bHLH (CpbHLH) gene family, identifying a total of 131 CpbHLH genes across 11 chromosomes. Phylogenetic analysis classified these CpbHLH genes into 23 subfamilies, wherein most members within the same subfamily exhibited analogous intron/exon patterns and motif composition. Moreover, the expansion of the CpbHLH gene family was primarily driven by segmental duplication, with duplicated gene pairs experiencing purifying selection during evolution. Transcriptomic analysis revealed diverse expression patterns of CpbHLH genes in various tissues and distinct stages of Wintersweet flower development, thereby suggesting their involvement in a diverse array of physiological processes. Furthermore, yeast 2-hybrid assay demonstrated interaction between CpbHLH25 and CpbHLH59 (regulators of floral scent and color) as well as with CpbHLH112 and CpMYB2, suggesting potential coordinately regulation of secondary metabolites biosynthesis in Wintersweet flowers. Collectively, our comprehensive analysis provides valuable insights into the structural attributes, evolutionary dynamics, and expression profiles of the CpbHLH gene family, laying a solid foundation for further explorations of the multifaceted physiological and molecular roles of bHLH TFs in Wintersweet.

The Role of Light-Regulated Auxin Signaling in Root Development

The root is an important organ for obtaining nutrients and absorbing water and carbohydrates, and it depends on various endogenous and external environmental stimulations such as light, temperature, water, plant hormones, and metabolic constituents. Auxin, as an essential plant hormone, can mediate rooting under different light treatments. Therefore, this review focuses on summarizing the functions and mechanisms of light-regulated auxin signaling in root development. Some light-response components such as phytochromes (PHYs), cryptochromes (CRYs), phototropins (PHOTs), phytochrome-interacting factors (PIFs) and constitutive photo-morphorgenic 1 (COP1) regulate root development. Moreover, light mediates the primary root, lateral root, adventitious root, root hair, rhizoid, and seminal and crown root development via the auxin signaling transduction pathway. Additionally, the effect of light through the auxin signal on root negative phototropism, gravitropism, root greening and the root branching of plants is also illustrated. The review also summarizes diverse light target genes in response to auxin signaling during rooting. We conclude that the mechanism of light-mediated root development via auxin signaling is complex, and it mainly concerns in the differences in plant species, such as barley (Hordeum vulgare L.) and wheat (Triticum aestivum L.), changes of transcript levels and endogenous IAA content. Hence, the effect of light-involved auxin signaling on root growth and development is definitely a hot issue to explore in the horticultural studies now and in the future.

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.

The identification of conserved sequence features of co-translationally decayed mRNAs and upstream open reading frames in angiosperm transcriptomes

RNA turnover is essential in maintaining messenger RNA (mRNA) homeostasis during various developmental stages and stress responses. Co-translational mRNA decay (CTRD), a process in which mRNAs are degraded while still associated with translating ribosomes, has recently been discovered to function in yeast and three angiosperm transcriptomes. However, it is still unclear how prevalent CTRD across the plant lineage. Moreover, the sequence features of co-translationally decayed mRNAs have not been well-studied. Here, utilizing a collection of publicly available degradome sequencing datasets for another seven angiosperm transcriptomes, we have confirmed that CTRD is functioning in at least 10 angiosperms and likely throughout the plant lineage. Additionally, we have identified sequence features shared by the co-translationally decayed mRNAs in these species, implying a possible conserved triggering mechanism for this pathway. Given that degradome sequencing datasets can also be used to identify actively translating upstream open reading frames (uORFs), which are quite understudied in plants, we have identified numerous actively translating uORFs in the same 10 angiosperms. These findings reveal that actively translating uORFs are prevalent in plant transcriptomes, some of which are conserved across this lineage. We have also observed conserved sequence features in the regions flanking these uORFs' stop codons that might contribute to ribosome stalling at these sequences. Finally, we discovered that there were very few overlaps between the mRNAs harboring actively translating uORFs and those sorted into the co-translational decay pathway in the majority of the studied angiosperms, suggesting that these two processes might be nearly mutually exclusive in those species. In total, our findings provide the identification of CTRD and actively translating uORFs across a broad collection of plants and provide novel insights into the important sequence features associated with these collections of mRNAs and regulatory elements, respectively.

A tomato HD-zip I transcription factor, VAHOX1, acts as a negative regulator of fruit ripening

Homeodomain-leucine zipper (HD-Zip) transcription factors are only present in higher plants and are involved in plant development and stress responses. However, our understanding of their participation in the fruit ripening of economical plants, such as tomato (Solanum lycopersicum), remains largely unclear. Here, we report that VAHOX1, a member of the tomato HD-Zip I subfamily, was expressed in all tissues, was highly expressed in breaker+4 fruits, and could be induced by ethylene. RNAi repression of VAHOX1 (VAHOX1-RNAi) resulted in accelerated fruit ripening, enhanced sensitivity to ethylene, and increased total carotenoid content and ethylene production. Conversely, VAHOX1 overexpression (VAHOX1-OE) in tomato had the opposite effect. RNA-Seq results showed that altering VAHOX1 expression affected the transcript accumulation of a series of genes involved in ethylene biosynthesis and signal transduction and cell wall modification. Additionally, a dual-luciferase reporter assay, histochemical analysis of GUS activity and a yeast one-hybrid (Y1H) assay revealed that VAHOX1 could activate the expression of AP2a. Our findings may expand our knowledge about the physiological functions of HD-Zip transcription factors in tomato and highlight the diversities of transcriptional regulation during the fruit ripening process.

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.

Rapid growth of Moso bamboo (Phyllostachys edulis): Cellular roadmaps, transcriptome dynamics, and environmental factors

Moso bamboo (Phyllostachys edulis) shows remarkably rapid growth (114.5 cm/day), but the underlying biological mechanisms remain unclear. After examining more than 12,750 internodes from more than 510 culms from 17 Moso populations, we identified internode 18 as a representative internode for rapid growth. This internode includes a 2-cm cell division zone (DZ), a cell elongation zone up to 12 cm, and a secondary cell wall (SCW) thickening zone. These zones elongated 11.8 cm, produced approximately 570,000,000 cells, and deposited ∼28 mg g-1 dry weight (DW) lignin and ∼44 mg g-1 DW cellulose daily, far exceeding vegetative growth observed in other plants. We used anatomical, mathematical, physiological, and genomic data to characterize development and transcriptional networks during rapid growth in internode 18. Our results suggest that (1) gibberellin may directly trigger the rapid growth of Moso shoots, (2) decreased cytokinin and increased auxin accumulation may trigger cell DZ elongation, and (3) abscisic acid and mechanical pressure may stimulate rapid SCW thickening via MYB83L. We conclude that internode length involves a possible tradeoff mediated by mechanical pressure caused by rapid growth, possibly influenced by environmental temperature and regulated by genes related to cell division and elongation. Our results provide insight into the rapid growth of Moso bamboo.

Integrating ATAC-seq and RNA-seq Reveals the Dynamics of Chromatin Accessibility and Gene Expression in Apple Response to Drought

Drought resistance in plants is influenced by multiple signaling pathways that involve various transcription factors, many target genes, and multiple types of epigenetic modifications. Studies on epigenetic modifications of drought focus on DNA methylation and histone modifications, with fewer on chromatin remodeling. Changes in chromatin accessibility can play an important role in abiotic stress in plants by affecting RNA polymerase binding and various regulatory factors. However, the changes in chromatin accessibility during drought in apples are not well understood. In this study, the landscape of chromatin accessibility associated with the gene expression of apple (GL3) under drought conditions was analyzed by Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) and RNA-seq. Differential analysis between drought treatment and control identified 23,466 peaks of upregulated chromatin accessibility and 2447 peaks of downregulated accessibility. The drought-induced chromatin accessibility changed genes were mainly enriched in metabolism, stimulus, and binding pathways. By combining results from differential analysis of RNA-seq and ATAC-seq, we identified 240 genes with higher chromatin accessibility and increased gene expression under drought conditions that may play important functions in the drought response process. Among them, a total of nine transcription factor genes were identified, including ATHB7, HAT5, and WRKY26. These transcription factor genes are differentially expressed with different chromatin accessibility motif binding loci that may participate in apple response to drought by regulating downstream genes. Our study provides a reference for chromatin accessibility under drought stress in apples and the results will facilitate subsequent studies on chromatin remodelers and transcription factors.

Integrating ATAC-seq and RNA-seq Reveals the Dynamics of Chromatin Accessibility and Gene Expression in Apple Response to Drought

Drought resistance in plants is influenced by multiple signaling pathways that involve various transcription factors, many target genes, and multiple types of epigenetic modifications. Studies on epigenetic modifications of drought focus on DNA methylation and histone modifications, with fewer on chromatin remodeling. Changes in chromatin accessibility can play an important role in abiotic stress in plants by affecting RNA polymerase binding and various regulatory factors. However, the changes in chromatin accessibility during drought in apples are not well understood. In this study, the landscape of chromatin accessibility associated with the gene expression of apple (GL3) under drought conditions was analyzed by Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) and RNA-seq. Differential analysis between drought treatment and control identified 23,466 peaks of upregulated chromatin accessibility and 2447 peaks of downregulated accessibility. The drought-induced chromatin accessibility changed genes were mainly enriched in metabolism, stimulus, and binding pathways. By combining results from differential analysis of RNA-seq and ATAC-seq, we identified 240 genes with higher chromatin accessibility and increased gene expression under drought conditions that may play important functions in the drought response process. Among them, a total of nine transcription factor genes were identified, including ATHB7, HAT5, and WRKY26. These transcription factor genes are differentially expressed with different chromatin accessibility motif binding loci that may participate in apple response to drought by regulating downstream genes. Our study provides a reference for chromatin accessibility under drought stress in apples and the results will facilitate subsequent studies on chromatin remodelers and transcription factors.

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.

Genome-wide characterization of homeodomain-leucine zipper genes reveals RsHDZ17 enhances the heat tolerance in radish (Raphanus sativus L.)

Homeodomain-leucine zipper (HD-Zip) transcription factors are involved in various biological processes of plant growth, development, and abiotic stress response. However, how they regulate heat stress (HS) response remains largely unclear in plants. In this study, a total of 83 RsHD-Zip genes were firstly identified from the genome of Raphanus sativus. RNA-Seq, RT-qPCR and promoter activity assays revealed that RsHDZ17 from HD-Zip Class I was highly expressed under heat, salt, and Cd stresses. RsHDZ17 is a nuclear protein with transcriptional activity at the C-terminus. Ectopic overexpression (OE) of RsHDZ17 in Arabidopsis thaliana enhanced the HS tolerance by improving the survival rate, photosynthesis capacity, and scavenging for reactive oxygen species (ROS). In addition, transient OE of RsHDZ17 in radish cotyledons impeded cell injury and augmented ROS scavenging under HS. Moreover, yeast one-hybrid, dual-luciferase assay, and electrophoretic mobility shift assay revealed that RsHDZ17 could bind to the promoter of HSFA1e. Collectively, these pieces of evidence demonstrate that RsHDZ17 could play a positive role in thermotolerance, partially through up-regulation of the expression of HSFA1e in plants. These results provide novel insights into the role of HD-Zips in radish and facilitate genetical engineering and development of heat-tolerant radish in breeding programs.

SK, Sheoran S, K. UB, K. BN, Kumar S, Singh AN, Singh HV. Seed Longevity in Legumes: Deeper Insights Into Mechanisms and Molecular Perspectives

Sustainable agricultural production largely depends upon the viability and longevity of high-quality seeds during storage. Legumes are considered as rich source of dietary protein that helps to ensure nutritional security, but associated with poor seed longevity that hinders their performance and productivity in farmer's fields. Seed longevity is the key determinant to assure proper seed plant value and crop yield. Thus, maintenance of seed longevity during storage is of prime concern and a pre-requisite for enhancing crop productivity of legumes. Seed longevity is significantly correlated with other seed quality parameters such as germination, vigor, viability and seed coat permeability that affect crop growth and development, consequently distressing crop yield. Therefore, information on genetic basis and regulatory networks associated with seed longevity, as well as molecular dissection of traits linked to longevity could help in developing crop varieties with good storability. Keeping this in view, the present review focuses towards highlighting the molecular basis of seed longevity, with special emphasis on candidate genes and proteins associated with seed longevity and their interplay with other quality parameters. Further, an attempt was made to provide information on 3D structures of various genetic loci (genes/proteins) associated to seed longevity that could facilitate in understanding the interactions taking place within the seed at molecular level. This review compiles and provides information on genetic and genomic approaches for the identification of molecular pathways and key players involved in the maintenance of seed longevity in legumes, in a holistic manner. Finally, a hypothetical fast-forward breeding pipeline has been provided, that could assist the breeders to successfully develop varieties with improved seed longevity in legumes.

Identification of the yield traits related haplotype combinations of transcription factor genes TaHDZ34 in common wheat

A predominant objective in wheat breeding is improving yield-related traits. The homeodomain-leucine zipper (HD-Zip) transcription factor plays a significant role in plant growth and development. In this study, we cloned all homeologs of TaHDZ34, which is a member of the HD-Zip class IV transcription factor family in wheat (Triticum aestivum L.). Sequence polymorphism analysis showed that TaHDZ-A34, TaHDZ-B34, and TaHDZ-D34 formed five, six, and six haplotypes, respectively, and the genes were divided into two main haplotype groups. We also developed functional molecular markers. The TaHDZ34 genes were divided into eight main haplotype combinations. Association analysis and distinct population validation preliminarily indicated that TaHDZ34 genes modulate grain number per spike, effective spikelet number per spike, thousand kernel weight, and flag leaf area per plant in wheat. Hap-ABD was the most effective haplotype combination of TaHDZ34. Subcellular localization showed that TaHDZ-A34 was localized to the nucleus. The interacting proteins of TaHDZ-A34 were involved in protein synthesis/degradation, energy production and transportation, and photosynthesis. Geographic distribution and frequencies of TaHDZ34 haplotype combinations suggested that Hap-Abd and Hap-AbD were preferentially selected in Chinese wheat breeding programs. The high-yield-related haplotype combination Hap-ABD provided beneficial genetic resources for the marker-assisted selection of new wheat cultivars.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01298-5.

The Gastrodia menghaiensis (Orchidaceae) genome provides new insights of orchid mycorrhizal interactions

BACKGROUND

To illustrate the molecular mechanism of mycoheterotrophic interactions between orchids and fungi, we assembled chromosome-level reference genome of Gastrodia menghaiensis (Orchidaceae) and analyzed the genomes of two species of Gastrodia.

RESULTS

Our analyses indicated that the genomes of Gastrodia are globally diminished in comparison to autotrophic orchids, even compared to Cuscuta (a plant parasite). Genes involved in arbuscular mycorrhizae colonization were found in genomes of Gastrodia, and many of the genes involved biological interaction between Gatrodia and symbiotic microbionts are more numerous than in photosynthetic orchids. The highly expressed genes for fatty acid and ammonium root transporters suggest that fungi receive material from orchids, although most raw materials flow from the fungi. Many nuclear genes (e.g. biosynthesis of aromatic amino acid L-tryptophan) supporting plastid functions are expanded compared to photosynthetic orchids, an indication of the importance of plastids even in totally mycoheterotrophic species.

CONCLUSION

Gastrodia menghaiensis has the smallest proteome thus far among angiosperms. Many of the genes involved biological interaction between Gatrodia and symbiotic microbionts are more numerous than in photosynthetic orchids.

The Arabidopsis transcription factors AtPHL1 and AtHB23 act together promoting carbohydrate transport from pedicel-silique nodes to seeds

Carbohydrates are produced in green tissues through photosynthesis and then transported to sink tissues. Carbon partitioning is a strategic process, fine regulated, involving specific sucrose transporters in each connecting tissue. Here we report that a screening of an Arabidopsis transcription factor (TF) library using the homeodomain-leucine zipper I member AtHB23 as bait, allowed identifying the TF AtPHL1 interacting with the former. An independent Y2H assay, and in planta by BiFC, confirmed such interaction. AtHB23 and AtPHL1 coexpressed in the pedicel-silique nodes and the funiculus. Mutant plants (phl1, and amiR23) showed a marked reduction of lipid content in seeds, although lipid composition did not change compared to the wild type. While protein and carbohydrate contents were not significantly different between mutants and control mature seeds, we observed a reduced carbohydrate content in mutant plants young siliques (7 days after pollination). Moreover, using a CFDA probe, we revealed an impaired transport to the seeds, and the gene encoding the carbohydrate transporters SWEET10 and SWEET11, usually expressed in connecting tissues, was repressed in the amiR23 and phl1 mutant plants. Altogether, the results indicated that AtHB23 and AtPHL1 act together, promoting sucrose transport, and the lack of any of them provoked a reduction in seeds lipid content.

Hydro-Electro Hybrid Priming Promotes Carrot (Daucus carota L.) Seed Germination by Activating Lipid Utilization and Respiratory Metabolism

Carrot (Daucus carota L.) is widely cultivated as one of the most important root crops, and developing an effective presowing treatment method can promote the development of modern mechanized precision sowing. In the present study, a novel seed priming technology, named hydro-electro hybrid priming (HEHP), was used to promote the germination of carrot seeds. Seed germination experiments showed that HEHP was able to increase the germination index (GI) and vigor index (VI) by 3.1-fold and 6.8-fold, respectively, and the effect was significantly superior to that of hydro-priming (HYD) and electrostatic field treatment (EF). The consumption and utilization rate of seed storage reserves were also greatly improved. Meanwhile, both glyoxysomes and mitochondria were found to appear ahead of time in the endosperm cells of HEHP through observations of the subcellular structure of the endosperm. Activities of isocitrate lyase (ICL), NAD-dependent malate dehydrogenase (MDH), pyruvate kinase (PK), and alcohol dehydrogenase (ADH) were significantly increased by HEHP. From transcriptome results, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to the glyoxylate cycle, glycolysis, gluconeogenesis, and the citrate cycle were significantly enriched and real-time quantitative PCR (qRT-PCR) analysis confirmed the expression pattern of 15 critical differentially expressed genes (DEGs) in these pathways. All DEGs encoding MDH, phosphoenolpyruvate carboxykinase (PEPCK), and PK were upregulated in HEHP; thus, it is reasonable to infer that the transformation of malate, oxalacetate, phosphoenolpyruvate, and pyruvate in the cytoplasm may be pivotal for the energy supply during early germination. The results suggest that the optimal effect of HEHP is achieved by initiating stored lipid utilization and respiratory metabolism pathways related to germination.

Upstream Open Reading Frame Mediated Translation of WNK8 Is Required for ABA Response in Arabidopsis

With no lysine (K) (WNK) kinases comprise a family of serine/threonine kinases belonging to an evolutionary branch of the eukaryotic kinome. These special kinases contain a unique active site and are found in a wide range of eukaryotes. The model plant Arabidopsis has been reported to have 11 WNK members, of which WNK8 functions as a negative regulator of abscisic acid (ABA) signaling. Here, we found that the expression of WNK8 is post-transcriptionally regulated through an upstream open reading frame (uORF) found in its 5′ untranslated region (5′-UTR). This uORF has been predicted to encode a conserved peptide named CPuORF58 in both monocotyledons and dicotyledons. The analysis of the published ribosome footprinting studies and the study of the frameshift CPuORF58 peptide with altered repression capability suggested that this uORF causes ribosome stalling. Plants transformed with the native WNK8 promoter driving WNK8 expression were comparable with wild-type plants, whereas the plants transformed with a similar construct with mutated CPuORF58 start codon were less sensitive to ABA. In addition, WNK8 and its downstream target RACK1 were found to synergistically coordinate ABA signaling rather than antagonistically modulating glucose response and flowering in plants. Collectively, these results suggest that the WNK8 expression must be tightly regulated to fulfill the demands of ABA response in plants.

Overexpression of the CaHB12 transcription factor in cotton (Gossypium hirsutum) improves drought tolerance

The Coffea arabica HB12 gene (CaHB12), which encodes a transcription factor belonging to the HD-Zip I subfamily, is upregulated under drought, and its constitutive overexpression (35S:CaHB12^OX) improves the Arabidopsis thaliana tolerance to drought and salinity stresses. Herein, we generated transgenic cotton events constitutively overexpressing the CaHB12 gene, characterized these events based on their increased tolerance to water deficit, and exploited the gene expression level from the CaHB12 network. The segregating events Ev8.29.1, Ev8.90.1, and Ev23.36.1 showed higher photosynthetic yield and higher water use efficiency under severe water deficit and permanent wilting point conditions compared to wild-type plants. Under well-irrigated conditions, these three promising transformed events showed an equivalent level of Abscisic acid (ABA) and decreased Indole-3-acetic acid (IAA) accumulation, and a higher putrescine/(spermidine + spermine) ratio in leaf tissues was found in the progenies of at least two transgenic cotton events compared to non-transgenic plants. In addition, genes that are considered as modulated in the A. thaliana 35S:CaHB12^OX line were also shown to be modulated in several transgenic cotton events maintained under field capacity conditions. The upregulation of GhPP2C and GhSnRK2 in transgenic cotton events maintained under permanent wilting point conditions suggested that CaHB12 might act enhancing the ABA-dependent pathway. All these data confirmed that CaHB12 overexpression improved the tolerance to water deficit, and the transcriptional modulation of genes related to the ABA signaling pathway or downstream genes might enhance the defense responses to drought. The observed decrease in IAA levels indicates that CaHB12 overexpression can prevent leaf abscission in plants under or after stress. Thus, our findings provide new insights on CaHB12 gene and identify several promising cotton events for conducting field trials on water deficit tolerance and agronomic performance.

Characterization of the development dynamics within the linear growth bamboo leaf

The leaf is the main photosynthetic organ in plants, such as bamboo. Leaves from bamboo are used as a food additive. However, according to our investigation, to date there have been no reports concerning the leaf development of bamboo. By measuring over 7500 parenchymal cells, we discovered that the linear leaf growth of Pseudosasa japonica could be divided into three growth sections. The first one is a roughly 1-cm long division zone (DZ), containing about 1580 cells, located at the bottom of the leaf; the second one is an about 3-cm long elongation zone (EZ), with ~1905 cells, located above the DZ; and the last is a mature zone (MZ) in which cell elongation is completed. The cell production rate of the DZ was 25.33-35.81 cells per hour, with an average of 29.73; and the cell division rate was ~0.45 cells per cell every 24 h. PacBio and Illumina transcriptome sequencing found 21 933 unigenes expressed in these zones. Further analysis revealed a dynamic transcriptome, with transcripts for cell division in the DZ changing to transcripts for cell elongation, photosynthetic development, secondary metabolism, stress resistance, and nutrition transport toward the leaf distal. Those transcriptome transformations correlated well with the changes of relative water content, biomass accumulation, and cellulose crystal degree and were supported by quantitative polymerase chain reaction data. These results revealed a developmental gradient of the bamboo linear growth leaf, which offers a foundation to elucidate and engineer leaf development in bamboo, an economically valuable plant.

Molecular mechanisms of mutualistic and antagonistic interactions in a plant-pollinator association

Many insects metamorphose from antagonistic larvae into mutualistic adult pollinators, with reciprocal adaptation leading to specialized insect-plant associations. It remains unknown how such interactions are established at molecular level. Here we assemble high-quality genomes of a fig species, Ficus pumila var. pumila, and its specific pollinating wasp, Wiebesia pumilae. We combine multi-omics with validation experiments to reveal molecular mechanisms underlying this specialized interaction. In the plant, we identify the specific compound attracting pollinators and validate the function of several key genes regulating its biosynthesis. In the pollinator, we find a highly reduced number of odorant-binding protein genes and an odorant-binding protein mainly binding the attractant. During antagonistic interaction, we find similar chemical profiles and turnovers throughout the development of galled ovules and seeds, and a significant contraction of detoxification-related gene families in the pollinator. Our study identifies some key genes bridging coevolved mutualists, establishing expectations for more diffuse insect-pollinator systems.

A SNP of HD-ZIP I transcription factor leads to distortion of trichome morphology in cucumber (Cucumis sativus L.)

BACKGROUND

Trichomes are excellent model systems for the analysis of cell differentiation and play essential roles in plant protection. From cucumber inbred line 'WD1', we identified an EMS-induced trichome abnormally developing mutant, nps, which exhibited smaller, denser and no pyramid-shaped head trichomes.

RESULTS

Using F₂ and BC₁ populations constructed from a cross between nps and '9930', the genetic analysis showed that the nps trait is controlled by a single recessive nuclear gene. We identified CsNps by map-based cloning with 576 individuals of the F₂ population generated from the cross of nps and inbred line '9930'. The CsNps was located at a 13.4-kb genomic region on chromosome 3, which region contains three predicted genes. Sequence analysis showed that only one single nucleotide mutation (C → T) between 9930 and nps was found in the second exon of Csa3G748220, a plant-specific class I HD-Zip gene. The result of allelism test also indicated that nps is a novel allelic mutant of Mict (Micro-trichome). Thus, nps was renamed mict-L130F. By comparing the transcriptome of mict-L130F vs WD1 and 06-2 (mict) vs 06-1 (wildtype, near-isogenic line of 06-2), several potential target genes that may be related to trichome development were identified.

CONCLUSIONS

Our results demonstrate that Mict-L130F is involved in the morphogenesis of trichomes. Map-based cloning of the Mict-L130F gene could promote the study of trichome development in cucumber.

Promoter Architecture and Transcriptional Regulation of Genes Upregulated in Germination and Coleoptile Elongation of Diverse Rice Genotypes Tolerant to Submergence

Rice has the natural morphological adaptation to germinate and elongate its coleoptile under submerged flooding conditions. The phenotypic deviation associated with the tolerance to submergence at the germination stage could be due to natural variation. However, the molecular basis of this variation is still largely unknown. A comprehensive understanding of gene regulation of different genotypes that have diverse rates of coleoptile elongation can provide significant insights into improved rice varieties. To do so, publicly available transcriptome data of five rice genotypes, which have different lengths of coleoptile elongation under submergence tolerance, were analyzed. The aim was to identify the correlation between promoter architecture, associated with transcriptional and hormonal regulation, in diverse genotype groups of rice that have different rates of coleoptile elongation. This was achieved by identifying the putative cis-elements present in the promoter sequences of genes upregulated in each group of genotypes (tolerant, highly tolerant, and extremely tolerant genotypes). Promoter analysis identified transcription factors (TFs) that are common and unique to each group of genotypes. The candidate TFs that are common in all genotypes are MYB, bZIP, AP2/ERF, ARF, WRKY, ZnF, MADS-box, NAC, AS2, DOF, E2F, ARR-B, and HSF. However, the highly tolerant genotypes interestingly possess binding sites associated with HY5 (bZIP), GBF3, GBF4 and GBF5 (bZIP), DPBF-3 (bZIP), ABF2, ABI5, bHLH, and BES/BZR, in addition to the common TFs. Besides, the extremely tolerant genotypes possess binding sites associated with bHLH TFs such as BEE2, BIM1, BIM3, BM8 and BAM8, and ABF1, in addition to the TFs identified in the tolerant and highly tolerant genotypes. The transcriptional regulation of these TFs could be linked to phenotypic variation in coleoptile elongation in response to submergence tolerance. Moreover, the results indicate a cross-talk between the key TFs and phytohormones such as gibberellic acid, abscisic acid, ethylene, auxin, jasmonic acid, and brassinosteroids, for an altered transcriptional regulation leading to differences in germination and coleoptile elongation under submergence. The information derived from the current in silico analysis can potentially assist in developing new rice breeding targets for direct seeding.

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

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

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.

Understanding salt tolerance mechanism using transcriptome profiling and de novo assembly of wild tomato Solanum chilense

Soil salinity affects the plant growth and productivity detrimentally, but Solanum chilense, a wild relative of cultivated tomato (Solanum lycopersicum L.), is known to have exceptional salt tolerance. It has precise adaptations against direct exposure to salt stress conditions. Hence, a better understanding of the mechanism to salinity stress tolerance by S. chilense can be accomplished by comprehensive gene expression studies. In this study 1-month-old seedlings of S. chilense and S. lycopersicum were subjected to salinity stress through application of sodium chloride (NaCl) solution. Through RNA-sequencing here we have studied the differences in the gene expression patterns. A total of 386 million clean reads were obtained through RNAseq analysis using the Illumina HiSeq 2000 platform. Clean reads were further assembled de novo into a transcriptome dataset comprising of 514,747 unigenes with N50 length of 578 bp and were further aligned to the public databases. Genebank non-redundant (Nr), Viridiplantae, Gene Ontology (GO), KOG, and KEGG databases classification suggested enrichment of these unigenes in 30 GO categories, 26 KOG, and 127 pathways, respectively. Out of 265,158 genes that were differentially expressed in response to salt treatment, 134,566 and 130,592 genes were significantly up and down-regulated, respectively. Upon placing all the differentially expressed genes (DEG) in known signaling pathways, it was evident that most of the DEGs involved in cytokinin, ethylene, auxin, abscisic acid, gibberellin, and Ca²⁺ mediated signaling pathways were up-regulated. Furthermore, GO enrichment analysis was performed using REVIGO and up-regulation of multiple genes involved in various biological processes in chilense under salinity were identified. Through pathway analysis of DEGs, “Wnt signaling pathway” was identified as a novel pathway for the response to the salinity stress. Moreover, key genes for salinity tolerance, such as genes encoding proline and arginine metabolism, ROS scavenging system, transporters, osmotic regulation, defense and stress response, homeostasis and transcription factors were not only salt-induced but also showed higher expression in S. chilense as compared to S. lycopersicum. Thus indicating that these genes may have an important role in salinity tolerance in S. chilense. Overall, the results of this study improve our understanding on possible molecular mechanisms underlying salt tolerance in plants in general and tomato in particular.

uORFs: Important Cis-Regulatory Elements in Plants

Gene expression is regulated at many levels, including mRNA transcription, translation, and post-translational modification. Compared with transcriptional regulation, mRNA translational control is a more critical step in gene expression and allows for more rapid changes of encoded protein concentrations in cells. Translation is highly regulated by complex interactions between cis-acting elements and trans-acting factors. Initiation is not only the first phase of translation, but also the core of translational regulation, because it limits the rate of protein synthesis. As potent cis-regulatory elements in eukaryotic mRNAs, upstream open reading frames (uORFs) generally inhibit the translation initiation of downstream major ORFs (mORFs) through ribosome stalling. During the past few years, with the development of RNA-seq and ribosome profiling, functional uORFs have been identified and characterized in many organisms. Here, we review uORF identification, uORF classification, and uORF-mediated translation initiation. More importantly, we summarize the translational regulation of uORFs in plant metabolic pathways, morphogenesis, disease resistance, and nutrient absorption, which open up an avenue for precisely modulating the plant growth and development, as well as environmental adaption. Additionally, we also discuss prospective applications of uORFs in plant breeding.

Beyond Trees: Regulons and Regulatory Motif Characterization

Trees and their seeds regulate their germination, growth, and reproduction in response to environmental stimuli. These stimuli, through signal transduction, trigger transcription factors that alter the expression of various genes leading to the unfolding of the genetic program. A regulon is conceptually defined as a set of target genes regulated by a transcription factor by physically binding to regulatory motifs to accomplish a specific biological function, such as the CO-FT regulon for flowering timing and fall growth cessation in trees. Only with a clear characterization of regulatory motifs, can candidate target genes be experimentally validated, but motif characterization represents the weakest feature of regulon research, especially in tree genetics. I review here relevant experimental and bioinformatics approaches in characterizing transcription factors and their binding sites, outline problems in tree regulon research, and demonstrate how transcription factor databases can be effectively used to aid the characterization of tree regulons.

OsbHLH073 Negatively Regulates Internode Elongation and Plant Height by Modulating GA Homeostasis in Rice

Internode elongation is one of the key agronomic traits determining a plant’s height and biomass. However, our understanding of the molecular mechanisms controlling internode elongation is still limited in crop plant species. Here, we report the functional identification of an atypical basic helix-loop-helix transcription factor (OsbHLH073) through gain-of-function studies using overexpression (OsbHLH073-OX) and activation tagging (osbhlh073-D) lines of rice. The expression of OsbHLH073 was significantly increased in the osbhlh073-D line. The phenotype of osbhlh073-D showed semi-dwarfism due to deficient elongation of the first internode and poor panicle exsertion. Transgenic lines overexpressing OsbHLH073 confirmed the phenotype of the osbhlh073-D line. Exogenous gibberellic acid (GA₃) treatment recovered the semi-dwarf phenotype of osbhlh073-D plants at the seedling stage. In addition, quantitative expression analysis of genes involving in GA biosynthetic and signaling pathway revealed that the transcripts of rice ent-kaurene oxidases 1 and 2 (OsKO1 and OsKO2) encoding the GA biosynthetic enzyme were significantly downregulated in osbhlh073-D and OsbHLH073-OX lines. Yeast two-hybrid and localization assays showed that the OsbHLH073 protein is a nuclear localized-transcriptional activator. We report that OsbHLH073 participates in regulating plant height, internode elongation, and panicle exsertion by regulating GA biosynthesis associated with the OsKO1 and OsKO2 genes.

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 AtHB1 Transcription Factor Controls the miR164-CUC2 Regulatory Node to Modulate Leaf Development

The presence of small tooth-like indentations, or serrations, characterizes leaf margins of Arabidopsis thaliana plants. The NAC family member CUP-SHAPED COTYLEDON 2 (CUC2), which undergoes post-transcriptional gene silencing by three micro-RNA genes (MIR164A, B and C), controls the extension of leaf serration. Here, we analyzed the role of AtHB1, a transcription factor (TF) belonging to the homeodomain-leucine zipper subfamily I, in shaping leaf margins. Using mutants with an impaired silencing pathway as background, we obtained transgenic plants expressing AtHB1 over 100 times compared to controls. These plants presented an atypical developmental phenotype characterized by leaves with deep serration. Transcript measurements revealed that CUC2 expression was induced in plants overexpressing AtHB1 and repressed in athb1 mutants, indicating a positive regulation exerted by this TF. Moreover, molecular analyses of AtHB1 overexpressing and mutant plants revealed that AtHB1 represses MIR164 transcription. We found that overexpression of MIR164B was able to reverse the serration phenotype of plants overexpressing AtHB1. Finally, chromatin immunoprecipitation assays revealed that AtHB1 was able to bind in vivo the promoter regions of all three MIR164 encoding loci. Altogether, our results indicate that AtHB1 directly represses MIR164 expression to enhance leaf serration by increasing CUC2 levels.

Characterization of wheat homeodomain-leucine zipper family genes and functional analysis of TaHDZ5-6A in drought tolerance in transgenic Arabidopsis

BACKGROUND

Many studies in Arabidopsis and rice have demonstrated that HD-Zip transcription factors play important roles in plant development and responses to abiotic stresses. Although common wheat (Triticum aestivum L.) is one of the most widely cultivated and consumed food crops in the world, the function of the HD-Zip proteins in wheat is still largely unknown.

RESULTS

To explore the potential biological functions of HD-Zip genes in wheat, we performed a bioinformatics and gene expression analysis of the HD-Zip family. We identified 113 HD-Zip members from wheat and classified them into four subfamilies (I-IV) based on phylogenic analysis against proteins from Arabidopsis, rice, and maize. Most HD-Zip genes are represented by two to three homeoalleles in wheat, which are named as TaHDZX_ZA, TaHDZX_ZB, or TaHDZX_ZD, where X denotes the gene number and Z the wheat chromosome on which it is located. TaHDZs in the same subfamily have similar protein motifs and intron/exon structures. The expression profiles of TaHDZ genes were analysed in different tissues, at different stages of vegetative growth, during seed development, and under drought stress. We found that most TaHDZ genes, especially those in subfamilies I and II, were induced by drought stress, suggesting the potential importance of subfamily I and II TaHDZ members in the responses to abiotic stress. Compared with wild-type (WT) plants, transgenic Arabidopsis plants overexpressing TaHDZ5-6A displayed enhanced drought tolerance, lower water loss rates, higher survival rates, and higher proline content under drought conditions. Additionally, the transcriptome analysis identified a number of differentially expressed genes between 35S::TaHDZ5-6A transgenic and wild-type plants, many of which are involved in stress response.

CONCLUSIONS

Our results will facilitate further functional analysis of wheat HD-Zip genes, and also indicate that TaHDZ5-6A may participate in regulating the plant response to drought stress. Our experiments show that TaHDZ5-6A holds great potential for genetic improvement of abiotic stress tolerance in crops.

Genome-Wide DNA Methylation Analysis of Soybean Curled-Cotyledons Mutant and Functional Evaluation of a Homeodomain-Leucine Zipper (HD-Zip) I Gene GmHDZ20

DNA methylation is a major, conserved epigenetic modification that influences many biological processes. Cotyledons are specialized tissues that provide nutrition for seedlings at the early developmental stage. To investigate the patterns of genomic DNA methylation of germinated cotyledons in soybean (Glycine max) and its effect on cotyledon development, we performed a genome-wide comparative analysis of DNA methylation between the soybean curled-cotyledons (cco) mutant, which has abnormal cotyledons, and its corresponding wild type (WT) by whole-genome bisulfite sequencing. The cco mutant was methylated at more sites but at a slightly lower level overall than the WT on the whole-genome level. A total of 46 CG-, 92 CHG-, and 9723 CHH- (H = A, C, or T) differentially methylated genes (DMGs) were identified in cotyledons. Notably, hypomethylated CHH-DMGs were enriched in the gene ontology term "sequence-specific DNA binding transcription factor activity." We selected a DMG encoding a homeodomain-leucine zipper (HD-Zip) I subgroup transcription factor (GmHDZ20) for further functional characterization. GmHDZ20 localized to the nucleus and was highly expressed in leaf and cotyledon tissues. Constitutive expression of GmHDZ20 in Arabidopsis thaliana led to serrated rosette leaves, shorter siliques, and reduced seed number per silique. A yeast two-hybrid assay revealed that GmHDZ20 physically interacted with three proteins associated with multiple aspects of plant growth. Collectively, our results provide a comprehensive study of soybean DNA methylation in normal and aberrant cotyledons, which will be useful for the identification of specific DMGs that participate in cotyledon development, and also provide a foundation for future in-depth functional study of GmHDZ20 in soybean.

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.

Genome-wide identification, characterization and expression analysis of the HD-Zip gene family in the stem development of the woody plant Prunus mume

The homeodomain-leucine zipper (HD-Zip) gene family, a group of plant-specific transcriptional factors (TFs), participates in regulating growth, development, and environmental responses. However, the characteristics and biological functions of HD-Zip genes in Prunus mume, which blooms in late winter or early spring, have not been reported. In this study, 32 HD-Zip genes, named PmHB1-PmHB32 based on their chromosomal positions, were identified in the genome of P. mume. These genes are distributed among seven chromosomes and are phylogenetically clustered into four major groups. Gene structure and motif composition were mostly conserved in each group. The Ka/Ks ratios showed that purifying selection has played a leading role in the long-term evolution of the genes, which maintained the function of this family. MicroRNA target site prediction indicated that the genes of the HD-Zip III subfamily may be regulated by miR165/166. Expression pattern analysis showed that the 32 genes were differentially expressed across five different tissues (leaf, flower bud, stem, fruit, and root) and at different stages of stem and leaf-bud development, suggesting that 10 of the genes may play important roles in stem development. Protein-protein interaction predictions showed that the subfamily III genes may regulate vascular development and shoot apical meristem (SAM) maintenance. Promoter analysis showed that the HD-Zip III genes might be involved in responses to light, hormones, and abiotic stressors and stem development. Taken together, our results provide an overview of the HD-Zip family in P. mume and lay the foundation for the molecular breeding of woody ornamental plants.

Genome-Wide Identification and Expression Analysis of HD-ZIP I Gene Subfamily in Nicotiana tabacum

The homeodomain-leucine zipper (HD-Zip) gene family, whose members play vital roles in plant growth and development, and participate in responding to various stresses, is an important class of transcription factors currently only found in plants. Although the HD-Zip gene family, especially the HD-Zip I subfamily, has been extensively studied in many plant species, the systematic report on HD-Zip I subfamily in cultivated tobacco (Nicotiana tabacum) is lacking. In this study, 39 HD-Zip I genes were systematically identified in N. tabacum (Nt). Interestingly, that 64.5% of the 31 genes with definite chromosome location information were found to originate from N. tomentosoformis, one of the two ancestral species of allotetraploid N. tabacum. Phylogenetic analysis divided the NtHD-Zip I subfamily into eight clades. Analysis of gene structures showed that NtHD-Zip I proteins contained conserved homeodomain and leucine-zipper domains. Three-dimensional structure analysis revealed that most NtHD-Zip I proteins in each clade, except for those in clade η, share a similar structure to their counterparts in Arabidopsis. Prediction of cis-regulatory elements showed that a number of elements responding to abscisic acid and different abiotic stresses, including low temperature, drought, and salinity, existed in the promoter region of NtHD-Zip I genes. The prediction of Arabidopsis ortholog-based protein-protein interaction network implied that NtHD-Zip I proteins have complex connections. The expression profile of these genes showed that different NtHD-Zip I genes were highly expressed in different tissues and could respond to abscisic acid and low-temperature treatments. Our study provides insights into the evolution and expression patterns of NtHD-Zip I genes in N. tabacum and will be useful for further functional characterization of NtHD-Zip I genes in the future.

Shade delays flowering in Medicago sativa

Shade-intolerant plants respond to the decrease in the red (R) to far-red (FR) light ratio (R:FR) occurring under shade by elongating stems and petioles and by re-positioning leaves, in a race to outcompete neighbors for the sunlight resource. In some annual species, the shade avoidance syndrome (SAS) is accompanied by the early induction of flowering. Anticipated flowering is viewed as a strategy to set seeds before the resources become severely limiting. Little is known about the molecular mechanisms of SAS in perennial forage crops like alfalfa (Medicago sativa). To study SAS in alfalfa, we exposed alfalfa plants to simulated shade by supplementing with FR light. Low R:FR light produced a classical SAS, with increased internode and petiole lengths, but unexpectedly also with delayed flowering. To understand the molecular mechanisms involved in uncoupling SAS from early flowering, we used a transcriptomic approach. The SAS is likely to be mediated by increased expression of msPIF3 and msHB2 in low R:FR light. Constitutive expression of these genes in Arabidopsis led to SAS, including early flowering, strongly suggesting that their roles are conserved. Delayed flowering was likely to be mediated by the downregulation of msSPL3, which promotes flowering in both Arabidopsis and alfalfa. Shade-delayed flowering in alfalfa may be important to extend the vegetative phase under suboptimal light conditions, and thus assure the accumulation of reserves necessary to resume growth after the next season.

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.

Transcriptome Analysis of Intrusively Growing Flax Fibers Isolated by Laser Microdissection

The intrusive growth, a type of plant cell elongation occurring in the depths of plant tissues, is characterized by the invasion of a growing cell between its neighbours due to a higher rate of elongation. In order to reveal the largely unknown molecular mechanisms of intrusive growth, we isolated primary flax phloem fibers specifically at the stage of intrusive growth by laser microdissection. The comparison of the RNA-Seq data from several flax stem parts enabled the characterization of those processes occurring specifically during the fiber intrusive elongation. The revealed molecular players are summarized as those involved in the supply of assimilates and support of turgor pressure, cell wall enlargement and modification, regulation by transcription factors and hormones, and responses to abiotic stress factors. The data obtained in this study provide a solid basis for developing approaches to manipulate fiber intrusive elongation, which is of importance both for plant biology and the yield of fiber crops.

De novo transcriptome assembly and identification of salt-responsive genes in sugar beet M14

Sugar beet (Beta vulgaris) is an important crop of sugar production in the world. Previous studies reported that sugar beet monosomic addition line M14 obtained from the intercross between Beta vulgaris L. (cultivated species) and B. corolliflora Zoss (wild species) exhibited tolerance to salt (up to 0.5 M NaCl) stress. To estimate a broad spectrum of genes involved in the M14 salt tolerance will help elucidate the molecular mechanisms underlying salt stress. Comparative transcriptomics was performed to monitor genes differentially expressed in the leaf and root samples of the sugar beet M14 seedlings treated with 0, 200 and 400 mM NaCl, respectively. Digital gene expression revealed that 3856 unigenes in leaves and 7157 unigenes in roots were differentially expressed under salt stress. Enrichment analysis of the differentially expressed genes based on GO and KEGG databases showed that in both leaves and roots genes related to regulation of redox balance, signal transduction, and protein phosphorylation were differentially expressed. Comparison of gene expression in the leaf and root samples treated with 200 and 400 mM NaCl revealed different mechanisms for coping with salt stress. In addition, the expression levels of nine unigenes in the reactive oxygen species (ROS) scavenging system exhibited significant differences in the leaves and roots. Our transcriptomics results have provided new insights into the salt-stress responses in the leaves and roots of sugar beet.

Arabidopsis thaliana homeodomain-leucine zipper type I transcription factors contribute to control leaf venation patterning

Venation patterning is a taxonomic attribute for classification of plants and it also plays a role in the interaction of plants with the environment. Despite its importance, the molecular physiology controlling this aspect of plant development is still poorly understood. Auxin plays a central role modulating the final vein network and patterning. This addendum discusses recent findings on the role of homeodomain-leucine zipper (HD-Zip) transcription factors on the regulation of leaf venation patterning. Moreno-Piovano et al. reported that ectopic expression of a sunflower HD-Zip I gene, HaHB4, increased the asymmetry of leaf venation. Even more, this work showed that auxin transport in the leaf through LAX carriers controls venation patterning. Here, we provide evidence indicating that some Arabidopsis thaliana HD-Zip I genes play a role in the determination of the final leaf venation patterning. We propose that these genes contribute to regulate vein patterning, likely controlling auxin homeostasis.