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

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.

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 characterization and expression profiling of the HD-ZIP gene family in Acoraceae under salinity and cold stress

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

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

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

Identification and Characterization of the HD-Zip Gene Family and Dimerization Analysis of HB7 and HB12 in Brassica napus L

Homeodomain-leucine zipper (HD-Zip) genes encode plant-specific transcription factors, which play important roles in plant growth, development, and response to environmental stress. These genes have not been fully studied in allopolyploid Brassica napus, an important kind of oil crop. In this study, 165 HD-Zip genes were identified in B. napus and classified into four subfamilies. If proteins belong to the same subfamily, they exhibit similarities in gene structure, motifs, and domain distribution patterns. BnHD-Zip genes were unevenly distributed in the An and Cn subgenomes. Whole genome triplication (WGT) events may be major mechanisms accounting for this gene family expansion. Orthologous gene analysis showed that the process of this gene family expansion was accompanied by domain loss. We further found three genes homologous to HB7 and three genes homologous to HB12, all induced by PEG, ABA, and NaCl treatment. HB7 could not form homodimers but could form heterodimers with HB12 based on yeast two-hybrid assays. The results of this study provide valuable information for further exploration of the HD-Zip gene family in B. napus.

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.

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.

Identification and characterization of stress responsive homeodomain leucine zipper transcription factors in Medicago truncatula

BACKGROUND

Homeodomain leucine zipper (HD-ZIP) transcription factors play roles in regulating plant development and responses to abiotic stresses; however, how HD-ZIP genes in Medicago truncatula are involved in abiotic stress response remains elusive.

METHODS AND RESULTS

The HD-ZIP I genes in Medicago truncatula were identified and characterized, and their expression patterns in different tissues and under different abiotic stresses were analyzed. A total of 15 Medicago truncatula HD-ZIP I genes were identified and a phylogenetic analysis of HD-ZIP I proteins in Arabidopsis thaliana and Medicago truncatula was conducted. Fifteen HD-ZIP I genes showed diverse tissue preferences. Among them, expressions of MtHB22 and MtHB51 were specially detected in vegetative buds. In addition, they responded to various abiotic stresses, including salinity and osmotic stress and abscisic acid (ABA). For instance, MtHB7 and MtHB12 expression levels were found to be positively associated with salt, osmotic stress and ABA in both shoots and roots, while MtHB13 and MtHB23 were negatively associated with these stresses in Medicago truncatula.

CONCLUSION

The HD-ZIP I genes in Medicago truncatula are evolutionarily conserved, but also exhibit gene duplication and gene loss events. Differential expression analysis of Medicago truncatula HD-ZIP I genes indicated their crucial roles in abiotic stress responses. Our genome-wide analysis of the HD-ZIP I transcription factor family in Medicago truncatula provided a valuable reference for further research.

Genome-Wide Identification of the HD-ZIP III Subfamily in Upland Cotton Reveals the Involvement of GhHB8-5D in the Biosynthesis of Secondary Wall in Fiber and Drought Resistance

A subfamily of transcription factors known as HD-ZIP III plays distinct roles in the secondary cell wall biosynthesis, which could be attributed to the quality of cotton fiber and adaptation to drought stress. In this study, 18 HD-ZIP III genes were identified as genome wide from the upland cotton (Gossypium hirsutum). These genes are distributed on 14 different chromosomes, and all of them have undergone segmental duplications. Numerous cis-elements were identified in the promoter regions, which are related to phytohormone responses and abiotic stresses. Expression profiling of these genes by quantitative real-time (qRT)-PCR illustrated their differential spatial expression, with preferential expression in cotton fiber. Among these genes, GhHB8-5D was predicted to encode a protein that is targeted to the cell nucleus and having self-activation ability. In addition, the ectopic expression of GhHB8-5D or its synonymous mutant GhHB8-5Dm in Arabidopsis resulted in stunted plant growth, curly leaves, and twisted inflorescence stems. Microscopy examination revealed that the morphology of vascular bundles and deposition of secondary wall had substantially altered in stems, which is concomitant with the significant alteration in the transcription levels of secondary wall-related genes in these transgenic Arabidopsis. Further, ectopic expression of GhHB8-5D or GhHB8-5Dm in Arabidopsis also led to significant increase in green seedling rate and reduction in root length relative to wild type when the plants were grown under mimicked drought stress conditions. Taken together, our results may shed new light on the functional roles of GhHB8-5D that is attributable for secondary cell wall thickening in response to drought stress. Such a finding may facilitate a novel strategy for improving plant adaptations to environmental changes via regulating the biosynthesis of secondary cell wall.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

Identification and analysis of Chrysanthemum nankingense NAC transcription factors and an expression analysis of OsNAC7 subfamily members

NAC (NAM, ATAF1-2, and CUC2) transcription factors (TFs) play a vital role in plant growth and development, as well as in plant response to biotic and abiotic stressors (Duan et al., 2019; Guerin et al., 2019). Chrysanthemum is a plant with strong stress resistance and adaptability; therefore, a systematic study of NAC TFs in chrysanthemum is of great significance for plant breeding. In this study, 153 putative NAC TFs were identified based on the Chrysanthemum nankingense genome. According to the NAC family in Arabidopsis and rice, a rootless phylogenetic tree was constructed, in which the 153 CnNAC TFs were divided into two groups and 19 subfamilies. Moreover, the expression levels of 12 CnNAC TFs belonging to the OsNAC7 subfamily were analyzed in C. nankingense under osmotic and salt stresses, and different tissues were tested during different growth periods. The results showed that these 12 OsNAC7 subfamily members were involved in the regulation of root and stem growth, as well as in the regulation of drought and salt stresses. Finally, we investigated the function of the CHR00069684 gene, and the results showed that CHR00069684 could confer improved salt and low temperature resistance, enhance ABA sensitivity, and lead to early flowering in tobacco. It was proved that members of the OsNAC7 subfamily have dual functions including the regulation of resistance and the mediation of plant growth and development. This study provides comprehensive information on analyzing the function of CnNAC TFs, and also reveals the important role of OsNAC7 subfamily genes in response to abiotic stress and the regulation of plant growth. These results provide new ideas for plant breeding to control stress resistance and growth simultaneously.

Wood development regulators involved in apical growth in Pinus canariensis

The shoot apical meristem is responsible of seasonal length increase in plants. In woody plants transition from primary to secondary growth is also produced during seasonal apical growth. These processes are controlled by different families of transcription factors. Levels of transcriptomic activity during apical growth were measured by means of a cDNA microarray designed from sequences related to meristematic activity in Pinus canariensis. The identification of differentially expressed genes was performed using a time-course analysis. A total of 7170 genes were differentially expressed and grouped in six clusters according to their expression profiles. We identified master regulators, such as WUSCHEL-like HOMEOBOX (WOX), to be involved in the first stages of apical development, i.e. growth of primary tissues, while other transcription factors, such as Class III HOMEODOMAIN-LEUCINE ZIPPER (HD-ZIP III) and KNOTTED-like (KNOX) and BEL1-like (BELL) HOMEODOMAIN proteins, were found to be induced during last stages of apical seasonal development, already with secondary growth. Our results reveal the main expression patterns of these genes during apical development and the transition from primary to secondary stem growth. In particular, the regulatory factors identified play key roles in controlling stem architecture and constitute candidate genes for the study of other development processes in conifers.

Identification and Bioinformatic Analysis of the GmDOG1-Like Family in Soybean and Investigation of Their Expression in Response to Gibberellic Acid and Abscisic Acid

Seed germination is one of the most important stages during plant life cycle, and DOG1 (Delay of germination1) plays a pivotal regulatory role in seed dormancy and germination. In this study, we have identified the DOG1-Like (DOG1L) family in soybean (Glycine max), a staple oil crop worldwide, and investigated their chromosomal distribution, structure and expression patterns. The results showed that the GmDOG1L family is composed of 40 members, which can be divided into six subgroups, according to their evolutionary relationship with other known DOG1-Like genes. These GmDOG1Ls are distributed on 18 of 20 chromosomes in the soybean genome and the number of exons for all the 40 GmDOG1Ls varied greatly. Members of the different subgroups possess a similar motif structure composition. qRT-PCR assay showed that the expression patterns of different GmDOG1Ls were significantly altered in various tissues, and some GmDOG1Ls expressed primarily in soybean seeds. Gibberellic acid (GA) remarkably inhibited the expression of most of GmDOG1Ls, whereas Abscisic acid (ABA) inhibited some of the GmDOG1Ls expression while promoting others. It is speculated that some GmDOG1Ls regulate seed dormancy and germination by directly or indirectly relating to ABA and GA pathways, with complex interaction networks. This study provides an important theoretical basis for further investigation about the regulatory roles of GmDOG1L family on soybean seed germination.

Identification of Homeobox Genes Associated with Lignification and Their Expression Patterns in Bamboo Shoots

: Homeobox (HB) genes play critical roles in regulating various aspects of plant growth and development. However, little is known about HB genes in bamboo. In this study, a total of 115 HB genes (PeHB001‒PeHB115) were identified from moso bamboo (Phyllostachys edulis) and grouped into 13 distinct classes (BEL, DDT, HD-ZIP I‒IV, KNOX, NDX, PHD, PINTOX, PLINC, SAWADEE, and WOX) based on the conserved domains and phylogenetic analysis. The number of members in the different classes ranged from 2 to 24, and they usually varied in terms of exon‒intron distribution pattern and length. There were 20 conserved motifs found in 115 PeHBs, with motif 1 being the most common. Gene ontology (GO) analysis showed that PeHBs had diverse molecular functions, with 19 PeHBs being annotated as having xylem development, xylem, and phloem pattern formation functions. Co-expression network analysis showed that 10 of the 19 PeHBs had co-expression correlations, and three members of the KNOX class were hub proteins that interacted with other transcription factors (TFs) such as MYB, bHLH, and OVATE, which were associated with lignin synthesis. Yeast two-hybridization results further proved that PeHB037 (BEL class) interacted with PeHB057 (KNOX class). Transcriptome expression profiling indicated that all PeHBs except PeHB017 were expressed in at least one of the seven tissues of moso bamboo, and 90 PeHBs were expressed in all the tissues. The qRT-PCR results of the 19 PeHBs showed that most of them were upregulated in shoots as the height increased. Moreover, a KNOX binding site was found in the promoters of the key genes involved in lignin synthesis such as Pe4CL, PeC3H, PeCCR, and PeCOMT, which had positive expression correlations with five KNOX genes. Similar results were found in winter bamboo shoots with prolonged storage time, which was consistent with the degree of lignification. These results provide basic data on PeHBs in moso bamboo, which will be helpful for future functional research on PeHBs with positive regulatory roles in the process of lignification.