The heat shock factor gene family in Salix suchowensis: a genome-wide survey and expression profiling during development and abiotic stresses

Genome-wide characterization of COMT family and regulatory role of CsCOMT19 in melatonin synthesis in Camellia sinensis

BACKGROUND

Caffeic acid O-methyltransferase (COMT) is a key enzyme that regulates melatonin synthesis and is involved in regulating the growth, development, and response to abiotic stress in plants. Tea plant is a popular beverage consumed worldwide, has been used for centuries for its medicinal properties, including its ability to reduce inflammation, improve digestion, and boost immune function. By analyzing genetic variation within the COMT family, while helping tea plants resist adversity, it is also possible to gain a deeper understanding of how different tea varieties produce and metabolize catechins, then be used to develop new tea cultivars with desired flavor profiles and health benefits.

RESULTS

In this study, a total of 25 CsCOMT genes were identified based on the high-quality tea (Camellia sinensis) plant genome database. Phylogenetic tree analysis of CsCOMTs with COMTs from other species showed that COMTs divided into four subfamilies (Class I, II, III, IV), and CsCOMTs was distributed in Class I, Class II, Class III. CsCOMTs not only undergoes large-scale gene recombination in pairs internally in tea plant, but also shares 2 and 7 collinear genes with Arabidopsis thaliana and poplar (Populus trichocarpa), respectively. The promoter region of CsCOMTs was found to be rich in cis-acting elements associated with plant growth and stress response. By analyzing the previously transcriptome data, it was found that some members of CsCOMT family exhibited significant tissue-specific expression and differential expression under different stress treatments. Subsequently, we selected six CsCOMTs to further validated their expression levels in different tissues organ using qRT-PCR. In addition, we silenced the CsCOMT19 through virus-induced gene silencing (VIGS) method and found that CsCOMT19 positively regulates the synthesis of melatonin in tea plant.

CONCLUSION

These results will contribute to the understanding the functions of CsCOMT gene family and provide valuable information for further research on the role of CsCOMT genes in regulating tea plant growth, development, and response to abiotic stress.

An insight into Pisum sativum HSF gene family-Genome-wide identification, phylogenetic, expression, and analysis of transactivation potential of pea heat shock transcription factor

Field pea (Pisum sativum L, 2n = 14) is a popular temperate legume with high economic value. Heat shock factors (HSFs) are the core element in the regulatory mechanism of heat stress responses. HSFs in pea (P. sativum) have not been characterized and their role remains unclear in different abiotic stresses. To address this knowledge gap, the current study aimed to characterize the HSF gene family in pea. We identified 38 PsHsf members in P. sativum, which are distributed on the seven chromosomes, and based on phylogenetic analysis, we classified them into three representative classes i.e. A, B, and C. Conserved motif and gene structure analysis confirmed a high degree of similarity among the members of the same class. Additionally, identified cis-acting regulatory elements (CAREs) related to abiotic responses, development, growth, and hormone signaling provides crucial insights into the regulatory mechanisms of PsHsfs. Our research revealed instances of gene duplication in PsHsf gene family, suggesting that this mechanism could be driving the expansion of the PsHsf gene family. Moreover, Expression analysis of PsHsfs exhibited upregulation under heat stress (HS), salt stress (SS), and drought stress (DS) showing their phenomenal role in stress conditions. PsHsfs protein interaction network suggested their involvement in stress-responsive mechanisms. Further transactivation potential was checked for spliced variant of PsHsfA2a (PsHsfA2aI, PsHsfA2aII, and PsHsfA2aIII), PsHsfA3, PsHsfA6b, PsHsfA9, PsHsfB1a, and PsHsfB2a. Overall, these findings provide valuable insight into the evolutionary relationship of PsHsf gene family and their role in abiotic stress responses.

In-silico identification and characterization of O-methyltransferase gene family in peanut (Arachis hypogaea L.) reveals their putative roles in development and stress tolerance

Cultivated peanut (Arachis hypogaea) is a leading protein and oil-providing crop and food source in many countries. At the same time, it is affected by a number of biotic and abiotic stresses. O-methyltransferases (OMTs) play important roles in secondary metabolism, biotic and abiotic stress tolerance. However, the OMT genes have not been comprehensively analyzed in peanut. In this study, we performed a genome-wide investigation of A. hypogaea OMT genes (AhOMTs). Gene structure, motifs distribution, phylogenetic history, genome collinearity and duplication of AhOMTs were studied in detail. Promoter cis-elements, protein-protein interactions, and micro-RNAs targeting AhOMTs were also predicted. We also comprehensively studied their expression in different tissues and under different stresses. We identified 116 OMT genes in the genome of cultivated peanut. Phylogenetically, AhOMTs were divided into three groups. Tandem and segmental duplication events played a role in the evolution of AhOMTs, and purifying selection pressure drove the duplication process. AhOMT promoters were enriched in several key cis-elements involved in growth and development, hormones, light, and defense-related activities. Micro-RNAs from 12 different families targeted 35 AhOMTs. GO enrichment analysis indicated that AhOMTs are highly enriched in transferase and catalytic activities, cellular metabolic and biosynthesis processes. Transcriptome datasets revealed that AhOMTs possessed varying expression levels in different tissues and under hormones, water, and temperature stress. Expression profiling based on qRT-PCR results also supported the transcriptome results. This study provides the theoretical basis for further work on the biological roles of AhOMT genes for developmental and stress responses.

Heat Shock Transcription Factor GhHSFB2a Is Crucial for Cotton Resistance to Verticillium dahliae

Heat shock transcription factors (HSFs) play a critical regulatory role in many plant disease resistance pathways. However, the molecular mechanisms of cotton HSFs involved in resistance to the soil-borne fungus Verticillium dahliae are limited. In our previous study, we identified numerous differentially expressed genes (DEGs) in the transcriptome and metabolome of V. dahliae-inoculated Arabidopsis thaliana. In this study, we identified and functionally characterized GhHSFB2a, which is a DEG belonging to HSFs and related to cotton immunity to V. dahliae. Subsequently, the phylogenetic tree of the type two of the HSFB subfamily in different species was divided into two subgroups: A. thaliana and strawberry, which have the closest evolutionary relationship to cotton. We performed promoter cis-element analysis and showed that the defense-reaction-associated cis-acting element-FC-rich motif may be involved in the plant response to V. dahliae in cotton. The expression pattern analysis of GhHSFB2a displayed that it is transcriptional in roots, stems, and leaves and significantly higher at 12 h post-inoculation (hpi). Subcellular localization of GhHSFB2a was observed, and the results showed localization to the nucleus. Virus-induced gene silencing (VIGS) analysis exhibited that GhHSFB2a silencing increased the disease index and fungal biomass and attenuated resistance against V. dahliae. Transcriptome sequencing of wild-type and GhHSFB2a-silenced plants, followed by Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, protein-protein interaction, and validation of marker genes revealed that ABA, ethylene, linoleic acid, and phenylpropanoid pathways are involved in GhHSFB2a-mediated plant disease resistance. Ectopic overexpression of the GhHSFB2a gene in Arabidopsis showed a significant increase in the disease resistance. Cumulatively, our results suggest that GhHSFB2a is required for the cotton immune response against V. dahliae-mediated ABA, ethylene, linoleic acid, and phenylpropanoid pathways, indicating its potential role in the molecular design breeding of plants.

Identification, classification, and expression profile analysis of heat shock transcription factor gene family in Salvia miltiorrhiza

In response to abiotic stresses, transcription factors are essential. Heat shock transcription factors (HSFs), which control gene expression, serve as essential regulators of plant growth, development, and stress response. As a model medicinal plant, Salvia miltiorrhiza is a crucial component in the treatment of cardiovascular illnesses. But throughout its growth cycle, S.miltiorrhiza is exposed to a series of abiotic challenges, including heat and drought. In this study, 35 HSF genes were identified based on genome sequencing of Salvia miltiorrhiza utilizing bioinformatics techniques. Additionally, 35 genes were classified into three groups by phylogeny and gene structural analysis, comprising 22 HSFA, 11 HSFB, and two HSFC. The distribution and sequence analysis of motif showed that SmHSFs were relatively conservative. In SmHSF genes, analysis of the promoter region revealed the presence of many cis-acting elements linked to stress, hormones, and growth and development, suggesting that these factors have regulatory roles. The majority of SmHSFs were expressed in response to heat and drought stress, according to combined transcriptome and real-time quantitative PCR (qRT-PCR) analyses. In conclusion, this study looked at the SmHSF gene family using genome-wide identification, evolutionary analysis, sequence characterization, and expression analysis. This research serves as a foundation for further investigations into the role of HSF genes and their molecular mechanisms in plant stress responses.

Global Integrated Genomic and Transcriptomic Analyses of MYB Transcription Factor Superfamily in C3 Model Plant Oryza sativa (L.) Unravel Potential Candidates Involved in Abiotic Stress Signaling

Plant transcription factors (TFs) are significant players in transcriptional regulations, signal transduction, and constitute an integral part of signaling networks. MYB TFs are major TF superfamilies that play pivotal roles in regulation of transcriptional reprogramming, physiological processes, and abiotic stress (AbS) responses. To explore the understanding of MYB TFs, genome and transcriptome-wide identification was performed in the C3 model plant, Oryza sativa (OsMYB). This study retrieved 114 OsMYB TFs that were computationally analyzed for their expression profiling, gene organization, cis-acting elements, and physicochemical properties. Based on the microarray datasets, six OsMYB genes which were sorted out and identified by a differential expression pattern were noted in various tissues. Systematic expression profiling of OsMYB TFs showed their meta-differential expression of different AbS treatments, spatio-temporal gene expression of various tissues and their growth in the field, and gene expression profiling in responses to phytohormones. In addition, the circular ideogram of OsMYB genes in related C4 grass plants conferred the gene synteny. Protein–protein interactions of these genes revealed the molecular crosstalk of OsMYB TFs. Transcriptional analysis (qPCR) of six OsMYB players in response to drought and salinity stress suggested the involvement in individual and combined AbS responses. To decipher how these OsMYB play functional roles in AbS dynamics, further research is a prerequisite.

Characteristics and Regulating Roles of Wheat TaHsfA2-13 in Abiotic Stresses

Heat shock transcription factor (Hsf) exists widely in eukaryotes and responds to various abiotic stresses by regulating the expression of downstream transcription factors, functional enzymes, and molecular chaperones. In this study, TaHsfA2-13, a heat shock transcription factor belonging to A2 subclass, was cloned from wheat (Triticum aestivum) and its function was analyzed. TaHsfA2-13 encodes a protein containing 368 amino acids and has the basic characteristics of Hsfs. Multiple sequence alignment analysis showed that TaHsfA2-13 protein had the highest similarity with TdHsfA2c-like protein from Triticum dicoccoides, which reached 100%. The analysis of tissue expression characteristics revealed that TaHsfA2-13 was highly expressed in root, shoot, and leaf during the seedling stage of wheat. The expression of TaHsfA2-13 could be upregulated by heat stress, low temperature, H₂O₂, mannitol, salinity and multiple phytohormones. The TaHsfA2-13 protein was located in the nucleus under the normal growth conditions and showed a transcriptional activation activity in yeast. Further studies found that overexpression of TaHsfA2-13 in Arabidopsis thaliana Col-0 or athsfa2 mutant results in improved tolerance to heat stress, H₂O₂, SA and mannitol by regulating the expression of multiple heat shock protein (Hsp) genes. In summary, our study identified TaHsfA2-13 from wheat, revealed its regulatory function in varieties of abiotic stresses, and will provide a new target gene to improve stress tolerance for wheat breeding.

Ectopic Overexpression of Maize Heat Stress Transcription Factor ZmHsf05 Confers Drought Tolerance in Transgenic Rice

Drought is a key factor affecting plant growth and development. Heat shock transcription factors (Hsfs) have been reported to respond to diverse abiotic stresses, including drought stress. In the present study, functional characterization of maize heat shock transcription factor 05 (ZmHsf05) gene was conducted. Homologous analysis showed that ZmHsf05 belongs to Class A2 Hsfs. The mRNA expression level of ZmHsf05 can be affected by drought, high temperature, salt, and abscisic acid (ABA) treatment. Ectopic overexpression of ZmHsf05 in rice (Oryza sativa) could significantly enhance the drought tolerance. Faced with drought stress, transgenic rice exhibited better phenotypic performance, higher survival rate, higher proline content, and lower leaf water loss rate, compared with wild-type plant Zhonghua11. Additionally, we assessed the agronomic traits of seven transgenic rice lines overexpressing ZmHsf05 and found that ZmHsf05 altered agronomical traits in the field trials. Moreover, rice overexpressing ZmHsf05 was more sensitive to ABA and had either a lower germination rate or shorter shoot length under ABA treatment. The transcription level of key genes in the ABA synthesis and drought-related pathway were significantly improved in transgenic rice after drought stress. Collectively, our results showed that ZmHsf05 could improve drought tolerance in rice, likely in an ABA-dependent manner.

Genome-wide identification and characterization of COMT gene family during the development of blueberry fruit

BACKGROUND

Caffeic acid O-methyltransferases (COMTs) play an important role in the diversification of natural products, especially in the phenylalanine metabolic pathway of plant. The content of COMT genes in blueberry and relationship between their expression patterns and the lignin content during fruit development have not clearly investigated by now.

RESULTS

Ninety-two VcCOMTs were identified in Vaccinium corymbosum. According to phylogenetic analyses, the 92 VcCOMTs were divided into 2 groups. The gene structure and conserved motifs within groups were similar which supported the reliability of the phylogenetic structure groupings. Dispersed duplication (DSD) and whole-genome duplication (WGD) were determined to be the major forces in VcCOMTs evolution. The results showed that the results of qRT-PCR and lignin content for 22 VcCOMTs, VcCOMT40 and VcCOMT92 were related to lignin content at different stages of fruit development of blueberry.

CONCLUSION

We identified COMT gene family in blueberry, and performed comparative analyses of the phylogenetic relationships in the 15 species of land plant, and gene duplication patterns of COMT genes in 5 of the 15 species. We found 2 VcCOMTs were highly expressed and their relative contents were similar to the variation trend of lignin content during the development of blueberry fruit. These results provide a clue for further study on the roles of VcCOMTs in the development of blueberry fruit and could promisingly be foundations for breeding blueberry clutivals with higher fruit firmness and longer shelf life.

Recent Advances in the Roles of HSFs and HSPs in Heat Stress Response in Woody Plants

A continuous increase in ambient temperature caused by global warming has been considered a worldwide threat. As sessile organisms, plants have evolved sophisticated heat shock response (HSR) to respond to elevated temperatures and other abiotic stresses, thereby minimizing damage and ensuring the protection of cellular homeostasis. In particular, for perennial trees, HSR is crucial for their long life cycle and development. HSR is a cell stress response that increases the number of chaperones including heat shock proteins (HSPs) to counter the negative effects on proteins caused by heat and other stresses. There are a large number of HSPs in plants, and their expression is directly regulated by a series of heat shock transcription factors (HSFs). Therefore, understanding the detailed molecular mechanisms of woody plants in response to extreme temperature is critical for exploring how woody species will be affected by climate changes. In this review article, we summarize the latest findings of the role of HSFs and HSPs in the HSR of woody species and discuss their regulatory networks and cross talk in HSR. In addition, strategies and programs for future research studies on the functions of HSFs and HSPs in the HSR of woody species are also proposed.

Genome-Wide Analysis of Heat Shock Transcription Factors in Ziziphus jujuba Identifies Potential Candidates for Crop Improvement Under Abiotic Stress

Plant heat shock transcription factors (Hsfs) play a significant role in adoption under abiotic stress conditions by modulating the expression of several stress-responsive genes. Analysis of the Hsf gene family will serve to understand the molecular mechanism which is involved in response to abiotic stress. The Ziziphus species grows in warm and dry regions and is inherently tolerant to abiotic stress conditions; thus, Ziziphus is a highly enriched source of genes conferring abiotic stress tolerance. Therefore, the present study provides a comprehensive genome-wide analysis of the Hsf gene family in Z. jujuba. Identified 21 non-redundant Hsf genes were grouped into three major classes (classes A, B, and C) based on the phylogenetic analysis. Promoter and gene ontology analysis suggested that ZjHsfs perform diverse functions in response to abiotic stress conditions. Two paralogous pairs resulting from tandem gene duplication events were identified. Also, physio-chemical properties of chromosomal locations, gene structure, motifs, and protein domain organization of Hsfs were analyzed. Real-time PCR expression analyses revealed that most of the Z. jujuba Hsf genes are differentially expressed in response to heat stress. The analysis suggested ZjHsf-2, ZjHsf-3, ZjHsf-5, ZjHsf-7, ZjHsf-8, ZjHsf-10, ZjHsf-12, ZjHsf-17, and ZjHsf-18 were the outstanding candidate genes for imparting heat stress tolerance and for future functional analysis. The present analysis laid the foundation for understanding the molecular mechanism of the Hsf gene family regulating Z. jujuba development and tolerance to abiotic stress conditions.

Uncovering candidate genes responsive to salt stress in Salix matsudana (Koidz) by transcriptomic analysis

Salix matsudana, a member of Salicaceae, is an important ornamental tree in China. Because of its capability to tolerate high salt conditions, S. matsudana also plays an important ecological role when grown along Chinese coastal beaches, where the salinity content is high. Here, we aimed to elucidate the mechanism of higher salt tolerance in S. matsudana variety ‘9901’ by identifying the associated genes through RNA sequencing and comparing differential gene expression between the S. matsudana salt-tolerant and salt-sensitive samples treated with 150 mM NaCl. Transcriptomic comparison of the roots of the two samples revealed 2174 and 3159 genes responsive to salt stress in salt-sensitive and salt-tolerant sample, respectively. Real-time polymerase chain reaction analysis of 9 of the responsive genes revealed a strong, positive correlation with RNA sequencing data. The genes were enriched in several pathways, including carbon metabolism pathway, plant-pathogen interaction pathway, and plant hormone signal transduction pathway. Differentially expressed genes (DEGs) encoding transcription factors associated with abiotic stress responses and salt stress response network were identified; their expression levels differed between the two samples in response to salt stress. Hub genes were also revealed by weighted gene co-expression network (WGCNA) analysis. For functional analysis of the DEG encoding sedoheptulose-1,7-bisphosphatase (SBPase), the gene was overexpressed in transgenic Arabidopsis, resulting in increased photosynthetic rates, sucrose and starch accumulation, and enhanced salt tolerance. Further functional characterization of other hub DEGs will reveal the molecular mechanism of salt tolerance in S. matsudana and allow the application of S. matsudana in coastal afforestation.

Genome-wide characterization of tea plant (Camellia sinensis) Hsf transcription factor family and role of CsHsfA2 in heat tolerance

BACKGROUND

Heat stress factors (Hsfs) play vital roles in signal transduction pathways operating in responses to environmental stresses. However, Hsf gene family has not been thoroughly explored in tea plant (Camellia sinensis L.).

RESULTS

In this study, we identified 25 CsHsf genes in C. sinensis that were separated by phylogenetic analysis into three sub-families (i.e., A, B, and C). Gene structures, conserved domains and motifs analyses indicated that the CsHsf members in each class were relatively conserved. Various cis-acting elements involved in plant growth regulation, hormone responses, stress responses, and light responses were located in the promoter regions of CsHsfs. Furthermore, degradome sequencing analysis revealed that 7 CsHsfs could be targeted by 9 miRNAs. The expression pattern of each CsHsf gene was significantly different in eight tissues. Many CsHsfs were differentially regulated by drought, salt, and heat stresses, as well as exogenous abscisic acid (ABA) and Ca²⁺. In addition, CsHsfA2 was located in the nucleus. Heterologous expression of CsHsfA2 improved thermotolerance in transgenic yeast, suggesting its potential role in the regulation of heat stress response.

CONCLUSIONS

A comprehensive genome-wide analysis of Hsf in C. sinensis present the global identification and functional prediction of CsHsfs. Most of them were implicated in a complex gene regulatory network controlling various abiotic stress responses and signal transduction pathways in tea plants. Additionally, heterologous expression of CsHsfA2 increased thermotolerance of transgenic yeast. These findings provide new insights into the functional divergence of CsHsfs and a basis for further research on CsHsfs functions.

Global transcriptome analysis of novel stress associated protein (SAP) genes expression dynamism of combined abiotic stresses in Oryza sativa (L.)

Genes encoding proteins with A20/AN1 zinc-finger domains, belonging to the stress associated protein (SAP) gene family, are present in all eukaryotes and play a decisive role in plant response to diverse physiological and molecular activities particularly on biotic and abiotic stresses (AbS). In this first and foremost study, global transcriptome analysis of members of the SAP gene family was carried out in C3 model-Oryza sativa (OsSAP) aiming at the identification of OsSAP genes activated in response to unique or Combined AbS (CAbS). Based on the available spatio-temporal and phytohormonal RNA-Seq expression profile datasets, nine OsSAP genes were filtered out and identified by a differential expression signature noted in various tissues as well as plant hormones. Comparative genome ideogram of OsSAP genes confirmed the orthologous collinearity with C4 panicoid genomes. Interactome of these genes, revealed the molecular cross-talks of OsSAP. Thus, the computational expression signature of OsSAP genes led to a better understanding of gene dynamism in diverse developmental tissues/organs. Transcriptional regulation analysis of key OsSAP genes in response to stress (drought and salinity) suggested the novel role of OsSAP1, OsSAP2, OsSAP5, OsSAP7, OsSAP8 and OsSAP11 in AbS. Altogether, the study provides deeper insights on molecular characteristics of OsSAP genes, which could be deployed further to decipher their precise functional roles in AbS responses.Communicated by Ramaswamy H. Sarma.

Genome-wide analysis and expression profiling of the heat shock transcription factor gene family in Physic Nut (Jatropha curcas L.)

The heat shock transcription factor (Hsf) family, identified as one of the important gene families, participates in plant development process and some stress response. So far, there have been no reports on the research of the Hsf transcription factors in physic nut. In this study, seventeen putative Hsf genes identified from physic nut genome. Phylogenetic analysis manifested these genes classified into three groups: A, B and C. Chromosomal location showed that they distributed eight out of eleven linkage groups. Expression profiling indicated that fourteen JcHsf genes highly expressed in different tissues except JcHsf1, JcHsf6 and JcHsf13. In addition, induction of six and twelve JcHsf genes noted against salt stress and drought stress, respectively, which demonstrated that the JcHsf genes are involved in abiotic stress responses. Our results contribute to a better understanding of the JcHsf gene family and further study of its function.

Heat stress transcription factor OsSPL7 plays a critical role in reactive oxygen species balance and stress responses in rice

The rice spotted leaf gene, OsSPL7, induces lesion mimic (LM) spots under heat stress. Herein, we provide several lines of evidence elucidating the importance of OsSPL7 in maintaining reactive oxygen species (ROS) balance via the regulation of downstream gene expression. osspl7 knockout (spl7ko) mutants showed LM and growth retardation. Transgenic rice lines strongly overexpressing OsSPL7 (SPL7OX-S) exhibited LM accompanied by accumulated H₂O₂, whereas moderate expressers of OsSPL7 (SPL7OX-M) did not, and neither of them exhibited severe growth defects. Transient expression of OsSPL7-GFP in rice protoplasts indicated that OsSPL7 localizes predominantly in the nucleus. Transcriptional activity assay suggested its function as a transcriptional activator in rice. Disease evaluation showed that both SPL7OX and spl7ko enhanced resistance to Magnaporthe oryzae and Xanthomonas oryzae pv. oryzae, the causal agents of blast and blight diseases in rice, respectively. Additionally, SPL7OX enhanced tolerance to cold stress, whereas spl7ko showed a phenotype opposite to the overexpression lines. RNA sequencing analyses identified four major groups of differentially expressed genes associated with LM, pathogen resistance, LM-pathogen resistance, and potential direct targets of OsSPL7. Collectively, our results suggest that OsSPL7 plays a critical role in plant growth and balancing ROS during biotic and abiotic stress.

Genome-wide identification and characterization of the MADS-box gene family in Salix suchowensis

MADS-box genes encode transcription factors that participate in various plant growth and development processes, particularly floral organogenesis. To date, MADS-box genes have been reported in many species, the completion of the sequence of the willow genome provides us with the opportunity to conduct a comprehensive analysis of the willow MADS-box gene family. Here, we identified 60 willow MADS-box genes using bioinformatics-based methods and classified them into 22 M-type (11 Mα, seven Mβ and four Mγ) and 38 MIKC-type (32 MIKCc and six MIKC*) genes based on a phylogenetic analysis. Fifty-six of the 60 SsMADS genes were randomly distributed on 19 putative willow chromosomes. By combining gene structure analysis with evolutionary analysis, we found that the MIKC-type genes were more conserved and played a more important role in willow growth. Further study showed that the MIKC* type was a transition between the M-type and MIKC-type. Additionally, the number of MADS-box genes in gymnosperms was notably lower than that in angiosperms. Finally, the expression profiles of these willow MADS-box genes were analysed in five different tissues (root, stem, leave, bud and bark) and validated by RT-qPCR experiments. This study is the first genome-wide analysis of the willow MADS-box gene family, and the results establish a basis for further functional studies of willow MADS-box genes and serve as a reference for related studies of other woody plants.

Genome-Wide Identification, Classification, and Expression Analysis of the Hsf Gene Family in Carnation (Dianthus caryophyllus)

Heat shock transcription factors (Hsfs) are a class of important transcription factors (TFs) which play crucial roles in the protection of plants from damages caused by various abiotic stresses. The present study aimed to characterize the Hsf genes in carnation (Dianthus caryophyllus), which is one of the four largest cut flowers worldwide. In this study, a total of 17 non-redundant Hsf genes were identified from the D. caryophyllus genome. Specifically, the gene structure and motifs of each DcaHsf were comprehensively analyzed. Phylogenetic analysis of the DcaHsf family distinctly separated nine class A, seven class B, and one class C Hsf genes. Additionally, promoter analysis indicated that the DcaHsf promoters included various cis-acting elements that were related to stress, hormones, as well as development processes. In addition, cis-elements, such as STRE, MYB, and ABRE binding sites, were identified in the promoters of most DcaHsf genes. According to qRT-PCR data, the expression of DcaHsfs varied in eight tissues and six flowering stages and among different DcaHsfs, even in the same class. Moreover, DcaHsf-A1, A2a, A9a, B2a, B3a revealed their putative involvement in the early flowering stages. The time-course expression profile of DcaHsf during stress responses illustrated that all the DcaHsfs were heat- and drought-responsive, and almost all DcaHsfs were down-regulated by cold, salt, and abscisic acid (ABA) stress. Meanwhile, DcaHsf-A3, A7, A9a, A9b, B3a were primarily up-regulated at an early stage in response to salicylic acid (SA). This study provides an overview of the Hsf gene family in D. caryophyllus and a basis for the breeding of stress-resistant carnation.

Genome-wide identification and classification of the Hsf and sHsp gene families in Prunus mume, and transcriptional analysis under heat stress

The transcriptional activation of heat shock proteins (Hsps) by heat shock transcription factors (Hsfs) is presumed to have a pivotal role in plant heat stress (HS) response. Prunus mume is an ornamental woody plant with distinctive features, including rich varieties and colors. In this study, 18 Hsfs and 24 small Hsps (sHsps) were identified in P. mume. Their chromosomal locations, protein domains, conserved motifs, phylogenetic relationships, and exon–intron structures were analyzed and compared with Arabidopsis thaliana Hsfs or sHsps. A total of 18 PmHsf members were classified into three major classes, A, B, and C. A total of 24 PmsHsps were grouped into eight subfamilies (CI to CIII, P, endoplasmic reticulum, M, and CI- or P-related). Quantitative reverse transcription PCR analysis revealed that members of the A2, A7, and A9 groups became the prominent Hsfs after heat shock, suggesting their involvement in a key regulatory role of heat tolerance. Most of the PmsHsp genes were up-regulated upon exposure to HS. Overall, our data contribute to an improved understanding of the complexity of the P. mume Hsf and sHsp gene families, and provide a basis for directing future systematic studies investigating the roles of the Hsf and sHsp gene families.

Genome-wide identification of the CLAVATA3/EMBRYO SURROUNDING REGION (CLE) family in grape (Vitis vinifera L.)

Background

CLE genes play various biological roles in plant growth and development, as well as in responses to environmental stimuli.

Results

In the present study, we identified nine CLE genes in the grape genome using an effective identification method. We analyzed the expression profiles of grape CLE genes in different tissues and under environmental different stimuli. VvCLE3 was expressed in shoot apical meristem (SAM) enriched regions, and VvCLE6 was expressed in shoot tissue without SAM. When grapes were infected with bois noir, VvCLE2 was up-regulated. Under ABA treatment, VvCLE3 was down-regulated. VvCLE6 was up-regulated under high temperature stress. We found that VvCLE6 and VvCLE1 were highly expressed in root tissue. In addition, we compared the characteristics of CLEs from grape and other plant species. The CLE family in Sphagnum fallax underwent positive selection, while the CLE family in grape underwent purifying selection. The frequency of optimal codons and codon adaptation index of rice and grape CLE family members were positively correlated with GC content at the third site of synonymous codons, indicating that the dominant evolutionary pressure acting on rice and grape CLE genes was mutation pressure. We also found that closely related species had higher levels of similarity in relative synonymous codon usage in CLE genes. The rice CLE family was biased toward C and G nucleotides at third codon positions. Gene duplication and loss events were also found in grape CLE genes.

Conclusion

These results demonstrate an effective identification method for CLE motifs and increase the understanding of grape CLEs. Future research on CLE genes may have applications for grape breeding and cultivation to better understand root and nodulation development.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5944-2) contains supplementary material, which is available to authorized users.

Evolutionary Divergence of Duplicated Hsf Genes in Populus

Heat shock transcription factors (Hsfs), which function as the activator of heat shock proteins (Hsps), play multiple roles in response to environmental stress and the development of plants. The Hsf family had experienced gene expansion via whole-genome duplication from a single cell algae to higher plants. However, how the Hsf gene family went through evolutionary divergence after genome duplication is unknown. As a model wood species, Populus trichocarpa is widely distributed in North America with various ecological and climatic environments. In this study, we used P. trichocarpa as materials and identified the expression divergence of the PtHsf gene family in developmental processes, such as dormant bud formation and opening, catkins development, and in response to environments. Through the co-expression network, we further discovered the divergent co-expressed genes that related to the functional divergence of PtHsfs. Then, we studied the alternative splicing events, single nucleotide polymorphism distribution and tertiary structures of members of the PtHsf gene family. In addition to expression divergence, we uncovered the evolutionary divergence in the protein level which may be important to new function formations and for survival in changing environments. This study comprehensively analyzed the evolutionary divergence of a member of the PtHsf gene family after genome duplication, paving the way for further gene function analysis and genetic engineering.

Genome-wide analysis of the Catalpa bungei caffeic acid O-methyltransferase (COMT) gene family: identification and expression profiles in normal, tension, and opposite wood

Caffeic acid O-methyltransferase (COMT) is an important protein that participates in lignin synthesis and is associated with the ratio of G-/S-type lignin in plants. COMTs are associated with the wood properties of forest trees; however, little known about the COMT family in Catalpa bungei, a valuable timber tree species in China . We performed a comprehensive analysis of COMT genes in the C. bungei genome by describing the gene structure and phylogenetic relationships of each family member using bioinformatics-based methods. A total of 23 putative COMT genes were identified using the conserved domain sequences and amino acid sequences of COMTs from Arabidopsis thaliana and Populus trichocarpa as probes. Phylogenetic analysis showed that 23 CbuCOMTs can be divided into three groups based on their structural characteristics; five conserved domains were found in the COMT family. Promoter analysis indicated that the CbuCOMT promoters included various cis-acting elements related to growth and development. Real-time quantitative polymerase chain reaction (PCR) analysis showed differential expression among CbuCOMTs. CbuCOMT2, 7, 8, 9, 10, 12, 13, 14, 21, and 23 were mainly expressed in xylem. Only CbuCOMT23 was significantly downregulated in tension wood and upregulated in opposite wood compared to normal wood. Our study provides new information about the CbuCOMT gene family and will facilitate functional characterisation in further research.

Genome-Wide Characterization of the sHsp Gene Family in Salix suchowensis Reveals Its Functions under Different Abiotic Stresses

Small heat shock proteins (sHsps) function mainly as molecular chaperones that play vital roles in response to diverse stresses, especially high temperature. However, little is known about the molecular characteristics and evolutionary history of the sHsp family in Salix suchowensis, an important bioenergy woody plant. In this study, 35 non-redundant sHsp genes were identified in S. suchowensis, and they were divided into four subfamilies (C, CP, PX, and MT) based on their phylogenetic relationships and predicted subcellular localization. Though the gene structure and conserved motif were relatively conserved, the sequences of the Hsp20 domain were diversified. Eight paralogous pairs were identified in the Ssu-sHsp family, in which five pairs were generated by tandem duplication events. Ka/Ks analysis indicated that Ssu-sHsps had undergone purifying selection. The expression profiles analysis showed Ssu-Hsps tissue-specific expression patterns, and they were induced by at least one abiotic stress. The expression correlation between two paralogous pairs (Ssu-sHsp22.2-CV/23.0-CV and 23.8-MT/25.6-MT) were less than 0.6, indicating that they were divergent during the evolution. Various cis-acting elements related to stress responses, hormone or development, were detected in the promoter of Ssu-sHsps. Furthermore, the co-expression network revealed the potential mechanism of Ssu-sHsps under stress tolerance and development. These results provide a foundation for further functional research on the Ssu-sHsp gene family in S. suchowensis.

Characterization and Functional Analysis of FaHsfC1b from Festuca arundinacea Conferring Heat Tolerance in Arabidopsis

Heat transcription factors (Hsfs) belong to a large gene family classified into A, B, and C groups, with classes A and B Hsfs being well-characterized and known for their roles in plant tolerance to abiotic stresses. The functions and roles of Class C Hsfs are not well-documented. The objectives of this study were to characterize a class C Hsf gene (FaHsfC1b) cloned from tall fescue (Festuca arundinacea), a perennial grass species, and to determine the physiological functions of FaHsfC1b in regulating heat tolerance by overexpressing FaHsfC1b in Arabidopsis thaliana. Full length cDNA of FaHsfC1b was cloned and the sequence alignment showed that it had high similarity to OsHsfC1b with typical DNA binding domain, hydrophobic oligomerization domain, and a nucleus localization signal. Transient expression with FaHsfC1b-eGFP in protoplasts of Arabidopsis leaves indicated its nucleus localization. qRT-PCR analysis showed that FaHsfC1b responded to heat, osmotic, salt, and cold stress in leaves and roots during 48-h treatment. Physiological analysis showed that FaHsfC1b overexpression enhanced plant survival rate, chlorophyll content, and photochemical efficiency, while it resulted in decreases in electrolyte leakage, H₂O₂ and O²⁻ content under heat stress. qRT-PCR showed that endogenous HsfC1 was induced in transgenic plants and the expression levels of heat protection protein genes, including several HSPs, AtGalSyn1, AtRof1, and AtHSA32, as well as ABA-synthesizing gene (NCED3) were significantly upregulated in transgenic plants overexpressing FaHsfC1b under heat stress. Our results first demonstrate that HsfC1b plays positive roles in plant tolerance to heat stress in association with the induction and upregulation of heat-protective genes. HsfC1b may be used as a candidate gene for genetic modification of cool-season plant species for improving heat tolerance.

Selection of suitable reference genes for quantitative real-time PCR gene expression analysis in Salix matsudana under different abiotic stresses

Salix matsudana is a deciduous, rapidly growing willow species commonly cultivated in China, which can tolerate drought, salt, and heavy metal stress conditions. Selection of suitable reference genes for quantitative real-time PCR is important for normalizing the expression of the key genes associated with various stresses. To validate suitable reference genes, we selected 11 candidate reference genes (five traditional housekeeping genes and six novel genes) and analyzed their expression stability in various samples, including different tissues and under different abiotic stress treatments. The expression of these genes was determined using five programs-geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder. The results showed that α-TUB2 (alpha-tubulin 2) and DnaJ (chaperone protein DnaJ 49) were the most stable reference genes across all the tested samples. We measured the expression profiles of the defense response gene SmCAT (catalase) using the two most stable and one least stable reference genes in all samples of S. matsudana. The relative quantification of SmCAT varied greatly according to the different reference genes. We propose that α-TUB2 and DnaJ should be the preferred reference genes for normalization and quantification of transcript levels in future gene expression studies in willow species under various abiotic stress conditions.

Divergent Expression Patterns in Two Vernicia Species Revealed the Potential Role of the Hub Gene VmAP2/ERF036 in Resistance to Fusarium oxysporum in Vernicia montana

Tung oil tree (Vernicia fordii) is a promising industrial oil crop; however, this tree is highly susceptible to Fusarium wilt disease. Conversely, Vernicia montana is resistant to the pathogen. The APETALA2/ethylene-responsive element binding factor (AP2/ERF) transcription factor superfamily has been reported to play a significant role in resistance to Fusarium oxysporum. In this study, comprehensive analysis identified 75 and 81 putative Vf/VmAP2/ERF transcription factor-encoding genes in V. fordii and V. montana, respectively, which were divided into AP2, ERF, related to ABI3 and VP1 (RAV) and Soloist families. After F. oxysporum infection, a majority of AP2/ERF superfamily genes showed strong patterns of repression in both V. fordii and V. montana. We then identified 53 pairs of one-to-one orthologs in V. fordii and V. montana, with most pairs of orthologous genes exhibiting similar expression in response to the pathogen. Further investigation of Vf/VmAP2/ERF gene expression in plant tissues indicated that the pairs of genes with different expression patterns in response to F. oxysporum tended to exhibit different tissue profiles in the two species. In addition, VmAP2/ERF036, showing the strongest interactions with 666 genes, was identified as a core hub gene mediating resistance. Moreover, qRT-PCR results indicated VmAP2/ERF036 showed repressed expression while its orthologous gene VfAP2/ERF036 had the opposite expression pattern during pathogen infection. Overall, comparative analysis of the Vf/VmAP2/ERF superfamily and indication of a potential hub resistance gene in resistant and susceptible Vernicia species provides valuable information for understanding the molecular basis and selection of essential functional genes for V. fordii genetic engineering to control Fusarium wilt disease.

Molecular evolution and expression divergence of the Populus euphratica Hsf genes provide insight into the stress acclimation of desert poplar

Heat shock transcription factor (Hsf) family is one of the most important regulators in the plant kingdom. Hsf has been demonstrated to be involved in various processes associated with plant growth, development as well as in response to hormone and abiotic stresses. In this study, we carried out a comprehensive analysis of Hsf family in desert poplar, Populus euphratica. Total of 32 genes encoding Hsf were identified and they were classified into three main classes (A, B, and C). Gene structure and conserved motif analyses indicated that the members in each class were relatively conserved. Total of 10 paralogous pairs were identified in PeuHsf family, in which nine pairs were generated by whole genome duplication events. Ka/Ks analysis showed that PeuHsfs underwent purifying selection pressure. In addition, various cis-acting elements involved in hormone and stress responses located in the promoter regions of PeuHsfs. Gene expression analysis indicated that several PeuHsfs were tissue-specific expression. Compared to Arabidopsis, more PeuHsf genes were significantly induced by heat, drought, and salt stresses (21, 19, and 22 PeuHsfs, respectively). Our findings are helpful in understanding the distinguished adaptability of P. euphratica to extreme environment and providing a basis for functional analysis of PeuHsfs in the future.

The Plant Heat Stress Transcription Factors (HSFs): Structure, Regulation, and Function in Response to Abiotic Stresses

Abiotic stresses such as high temperature, salinity, and drought adversely affect the survival, growth, and reproduction of plants. Plants respond to such unfavorable changes through developmental, physiological, and biochemical ways, and these responses require expression of stress-responsive genes, which are regulated by a network of transcription factors (TFs), including heat stress transcription factors (HSFs). HSFs play a crucial role in plants response to several abiotic stresses by regulating the expression of stress-responsive genes, such as heat shock proteins (Hsps). In this review, we describe the conserved structure of plant HSFs, the identification of HSF gene families from various plant species, their expression profiling under abiotic stress conditions, regulation at different levels and function in abiotic stresses. Despite plant HSFs share highly conserved structure, their remarkable diversification across plants reflects their numerous functions as well as their integration into the complex stress signaling and response networks, which can be employed in crop improvement strategies via biotechnological intervention.

Selection of Reliable Reference Genes for Gene Expression Analysis under Abiotic Stresses in the Desert Biomass Willow, Salix psammophila

Salix psammophila is a desert shrub willow that has extraordinary adaptation to abiotic stresses and plays an important role in maintaining local ecosystems. Moreover, S. psammophila is regarded as a promising biomass feedstock because of its high biomass yields and short rotation coppice cycle. However, few suitable reference genes (RGs) for quantitative real-time polymerase chain reaction (qRT-PCR) constrain the study on normalization of gene expression in S. psammophila until now. Here, we investigated the expression stabilities of 14 candidate RGs across tissue types and under four abiotic stress treatments, including heat, cold, salt, and drought treatments. After calculation of PCR efficiencies, three different software, NormFinder, geNorm, and BestKeeper were employed to analyze systematically the qRT-PCR data, and the outputs were merged by RankAggreg software. The optimal RGs selected for gene expression analysis were EF1α (Elongation factor-1 alpha) and OTU (OTU-like cysteine protease family protein) for different tissue types, UBC (Ubiquitin-conjugating enzyme E2) and LTA4H (Leukotriene A-4 hydrolase homolog) for heat treatment, HIS (Histone superfamily protein H3) and ARF2 (ADP-ribosylation factor 2) for cold treatment, OTU and ACT7 (Actin 7) for salt treatment, UBC and LTA4H for drought treatment. The expression of UBC, ARF2, and VHAC (V-type proton ATPase subunit C) varied the least across tissue types and under abiotic stresses. Furthermore, the relative genes expression profiles of one tissue-specific gene WOX1a (WUSCHEL-related homeobox 1a), and four stress-inducible genes, including Hsf-A2 (Heat shock transcription factors A2), CBF3 (C-repeat binding factor 3), HKT1 (High-Affinity K⁺ Transporter 1), and GST (Glutathione S-transferase), were conducted to confirm the validity of the RGs in this study. These results provided an important RGs application guideline for gene expression characterization in S. psammophila.