Identification, classification, and expression profiles of heat shock transcription factors in tea plant (Camellia sinensis) under temperature stress

A method for mining condition-specific co-expressed genes in Camellia sinensis based on k-means clustering

BACKGROUND

As one of the world's most important beverage crops, tea plants (Camellia sinensis) are renowned for their unique flavors and numerous beneficial secondary metabolites, attracting researchers to investigate the formation of tea quality. With the increasing availability of transcriptome data on tea plants in public databases, conducting large-scale co-expression analyses has become feasible to meet the demand for functional characterization of tea plant genes. However, as the multidimensional noise increases, larger-scale co-expression analyses are not always effective. Analyzing a subset of samples generated by effectively downsampling and reorganizing the global sample set often leads to more accurate results in co-expression analysis. Meanwhile, global-based co-expression analyses are more likely to overlook condition-specific gene interactions, which may be more important and worthy of exploration and research.

RESULTS

Here, we employed the k-means clustering method to organize and classify the global samples of tea plants, resulting in clustered samples. Metadata annotations were then performed on these clustered samples to determine the "conditions" represented by each cluster. Subsequently, we conducted gene co-expression network analysis (WGCNA) separately on the global samples and the clustered samples, resulting in global modules and cluster-specific modules. Comparative analyses of global modules and cluster-specific modules have demonstrated that cluster-specific modules exhibit higher accuracy in co-expression analysis. To measure the degree of condition specificity of genes within condition-specific clusters, we introduced the correlation difference value (CDV). By incorporating the CDV into co-expression analyses, we can assess the condition specificity of genes. This approach proved instrumental in identifying a series of high CDV transcription factor encoding genes upregulated during sustained cold treatment in Camellia sinensis leaves and buds, and pinpointing a pair of genes that participate in the antioxidant defense system of tea plants under sustained cold stress.

CONCLUSIONS

To summarize, downsampling and reorganizing the sample set improved the accuracy of co-expression analysis. Cluster-specific modules were more accurate in capturing condition-specific gene interactions. The introduction of CDV allowed for the assessment of condition specificity in gene co-expression analyses. Using this approach, we identified a series of high CDV transcription factor encoding genes related to sustained cold stress in Camellia sinensis. This study highlights the importance of considering condition specificity in co-expression analysis and provides insights into the regulation of the cold stress in Camellia sinensis.

Investigation of heat stress responses and adaptation mechanisms by integrative metabolome and transcriptome analysis in tea plants (Camellia sinensis)

Extreme high temperature has deleterious impact on the yield and quality of tea production, which has aroused the attention of growers and breeders. However, the mechanisms by which tea plant varieties respond to extreme environmental heat is not clear. In this study, we analyzed physiological indices, metabolites and transcriptome differences in three different heat-tolerant tea plant F1 hybrid progenies. Results showed that the antioxidant enzyme activity, proline, and malondialdehyde were significantly decreased in heat-sensitive 'FWS' variety, and the accumulation of reactive oxygen molecules such as H2O2 and O2- was remarkably increased during heat stress. Metabolomic analysis was used to investigate the metabolite accumulation pattern of different varieties in response to heat stress. The result showed that a total of 810 metabolites were identified and more than 300 metabolites were differentially accumulated. Transcriptional profiling of three tea varieties found that such genes encoding proteins with chaperon domains were preferentially expressed in heat-tolerant varieties under heat stress, including universal stress protein (USP32, USP-like), chaperonin-like protein 2 (CLP2), small heat shock protein (HSP18.1), and late embryogenesis abundant protein (LEA5). Combining metabolomic with transcriptomic analyses discovered that the flavonoids biosynthesis pathway was affected by heat stress and most flavonols were up-regulated in heat-tolerant varieties, which owe to the preferential expression of key FLS genes controlling flavonol biosynthesis. Take together, molecular chaperons, or chaperon-like proteins, flavonols accumulation collaboratively contributed to the heat stress adaptation in tea plant. The present study elucidated the differences in metabolite accumulation and gene expression patterns among three different heat-tolerant tea varieties under extreme ambient high temperatures, which helps to reveal the regulatory mechanisms of tea plant adaptation to heat stress, and provides a reference for the breeding of heat-tolerant tea plant varieties.

A review of changes at the phenotypic, physiological, biochemical, and molecular levels of plants due to high temperatures

MAIN CONCLUSION

This review summarizes the physiological, biochemical, and molecular regulatory network changes in plants in response to high temperature. With the continuous rise in temperature, high temperature has become an important issue limiting global plant growth and development, affecting the phenotype and physiological and biochemical processes of plants and seriously restricting crop yield and tree growth speed. As sessile organisms, plants inevitably encounter high temperatures and improve their heat tolerance by activating molecular networks related to heat stress, such as signal transduction, synthesis of metabolites, and gene expression. Heat tolerance is a polygenic trait regulated by a variety of genes, transcription factors, proteins, and metabolites. Therefore, this review summarizes the changes in physiological, biochemical and molecular regulatory networks in plants under high-temperature conditions to lay a foundation for an in-depth understanding of the mechanisms involved in plant heat tolerance responses.

Identification of genes associated with abiotic stress tolerance in sweetpotato using weighted gene co-expression network analysis

Sweetpotato, Ipomoea batatas (L.), a key food security crop, is negatively impacted by heat, drought, and salinity stress. The orange-fleshed sweetpotato cultivar "Beauregard" was exposed to heat, salt, and drought treatments for 24 and 48 h to identify genes responding to each stress condition in leaves. Analysis revealed both common (35 up regulated, 259 down regulated genes in the three stress conditions) and unique sets of up regulated (1337 genes by drought, 516 genes by heat, and 97 genes by salt stress) and down regulated (2445 genes by drought, 678 genes by heat, and 204 genes by salt stress) differentially expressed genes (DEGs) suggesting common, yet stress-specific transcriptional responses to these three abiotic stressors. Gene Ontology analysis of down regulated DEGs common to both heat and salt stress revealed enrichment of terms associated with "cell population proliferation" suggestive of an impact on the cell cycle by the two stress conditions. To identify shared and unique gene co-expression networks under multiple abiotic stress conditions, weighted gene co-expression network analysis was performed using gene expression profiles from heat, salt, and drought stress treated 'Beauregard' leaves yielding 18 co-expression modules. One module was enriched for "response to water deprivation," "response to abscisic acid," and "nitrate transport" indicating synergetic crosstalk between nitrogen, water, and phytohormones with genes encoding osmotin, cell expansion, and cell wall modification proteins present as key hub genes in this drought-associated module. This research lays the groundwork for exploring to a further degree, mechanisms for abiotic stress tolerance in sweetpotato.

Insights into morphological and physio-biochemical adaptive responses in mungbean (Vigna radiata L.) under heat stress

Mungbean (Vigna radiata L. Wilczek) is an important food legume crop which contributes significantly to nutritional and food security of South and Southeast Asia. The crop thrives in hot and humid weather conditions, with an optimal temperature range of 28°-35°C, and is mainly cultivated under rainfed environments. However, the rising global temperature has posed a serious threat to mungbean cultivation. Optimal temperature is a vital factor in cellular processes, and every crop species has evolved with its specific temperature tolerance ability. Moreover, variation within a crop species is inevitable, given the diverse environmental conditions under which it has evolved. For instance, various mungbean germplasm can grow and produce seeds in extreme ambient temperatures as low as 20°C or as high as 45°C. This range of variation in mungbean germplasm for heat tolerance plays a crucial role in developing heat tolerant and high yielding mungbean cultivars. However, heat tolerance is a complex mechanism which is extensively discussed in this manuscript; and at the same time individual genotypes have evolved with various ways of heat stress tolerance. Therefore, to enhance understanding towards such variability in mungbean germplasm, we studied morphological, anatomical, physiological, and biochemical traits which are responsive to heat stress in plants with more relevance to mungbean. Understanding heat stress tolerance attributing traits will help in identification of corresponding regulatory networks and associated genes, which will further help in devising suitable strategies to enhance heat tolerance in mungbean. The major pathways responsible for heat stress tolerance in plants are also discussed.

Conserved and variable heat stress responses of the Heat Shock Factor transcription factor family in maize and Setaria viridis

The Heat Shock Factor (HSF) transcription factor family is a central and required component of plant heat stress responses and acquired thermotolerance. The HSF family has dramatically expanded in plant lineages, often including a repertoire of 20 or more genes. Here we assess and compare the composition, heat responsiveness, and chromatin profiles of the HSF families in maize and Setaria viridis (Setaria), two model C4 panicoid grasses. Both species encode a similar number of HSFs, and examples of both conserved and variable expression responses to a heat stress event were observed between the two species. Chromatin accessibility and genome-wide DNA-binding profiles were generated to assess the chromatin of HSF family members with distinct responses to heat stress. We observed significant variability for both chromatin accessibility and promoter occupancy within similarly regulated sets of HSFs between Setaria and maize, as well as between syntenic pairs of maize HSFs retained following its most recent genome duplication event. Additionally, we observed the widespread presence of TF binding at HSF promoters in control conditions, even at HSFs that are only expressed in response to heat stress. TF-binding peaks were typically near putative HSF-binding sites in HSFs upregulated in response to heat stress, but not in stable or not expressed HSFs. These observations collectively support a complex scenario of expansion and subfunctionalization within this transcription factor family and suggest that within-family HSF transcriptional regulation is a conserved, defining feature of the family.

Transcriptional reprogramming of xylem cell wall biosynthesis in tension wood

Tension wood (TW) is a specialized xylem tissue developed under mechanical/tension stress in angiosperm trees. TW development involves transregulation of secondary cell wall genes, which leads to altered wood properties for stress adaptation. We induced TW in the stems of black cottonwood (Populus trichocarpa, Nisqually-1) and identified two significantly repressed transcription factor (TF) genes: class B3 heat-shock TF (HSFB3-1) and MYB092. Transcriptomic analysis and chromatin immunoprecipitation (ChIP) were used to identify direct TF-DNA interactions in P. trichocarpa xylem protoplasts overexpressing the TFs. This analysis established a transcriptional regulatory network in which PtrHSFB3-1 and PtrMYB092 directly activate 8 and 11 monolignol genes, respectively. The TF-DNA interactions were verified for their specificity and transactivator roles in 35 independent CRISPR-based biallelic mutants and overexpression transgenic lines of PtrHSFB3-1 and PtrMYB092 in P. trichocarpa. The gene-edited trees (mimicking the repressed PtrHSFB3-1 and PtrMYB092 under tension stress) have stem wood composition resembling that of TW during normal growth and under tension stress (i.e., low lignin and high cellulose), whereas the overexpressors showed an opposite effect (high lignin and low cellulose). Individual overexpression of the TFs impeded lignin reduction under tension stress and restored high levels of lignin biosynthesis in the TW. This study offers biological insights to further uncover how metabolism, growth, and stress adaptation are coordinately regulated in trees.

Genome-wide identification, classification, and expression analysis of the HSF gene family in pineapple (Ananas comosus)

Transcription factors (TFs), such as heat shock transcription factors (HSFs), usually play critical regulatory functions in plant development, growth, and response to environmental cues. However, no HSFs have been characterized in pineapple thus far. Here, we identified 22 AcHSF genes from the pineapple genome. Gene structure, motifs, and phylogenetic analysis showed that AcHSF families were distinctly grouped into three subfamilies (12 in Group A, seven in Group B, and four in Group C). The AcHSF promoters contained various cis-elements associated with stress, hormones, and plant development processes, for instance, STRE, WRKY, and ABRE binding sites. The majority of HSFs were expressed in diverse pineapple tissues and developmental stages. The expression of AcHSF-B4b/AcHSF-B4c and AcHSF-A7b/AcHSF-A1c were enriched in the ovules and fruits, respectively. Six genes (AcHSF-A1a , AcHSF-A2, AcHSF-A9a, AcHSF-B1a, AcHSF-B2a, and AcHSF-C1a) were transcriptionally modified by cold, heat, and ABA. Our results provide an overview and lay the foundation for future functional characterization of the pineapple HSF gene family.

Underpinning the molecular programming attributing heat stress associated thermotolerance in tea (Camellia sinensis (L.) O. Kuntze)

The most daunting issue of global climate change is the deleterious impact of extreme temperatures on tea productivity and quality, which has resulted in a quest among researchers and growers. The current study aims to unravel molecular programming underpinning thermotolerance by characterizing heat tolerance and sensitivity response in 20 tea cultivars. The significantly higher negative influence of heat stress was recorded in a sensitive cultivar with reduced water retention (47%), chlorophyll content (33.79%), oxidation potential (32.48%), and increase in membrane damage (76.4%). Transcriptional profiling of most tolerant and sensitive cultivars identified 78 differentially expressed unigenes with chaperon domains, including low and high molecular weight heat shock protein (HSP) and heat shock transcription factors (HSFs) involved in heat shock response (HSR). Further, predicted transcriptional interactome network revealed their key role in thermotolerance via well-co-ordinated transcriptional regulation of aquaporins, starch metabolism, chlorophyll biosynthesis, calcium, and ethylene mediated plant signaling system. The study identified the key role of HSPs (CsHSP90) in regulating HSR in tea, wherein, structure-based molecular docking revealed the inhibitory role of geldanamycin (GDA) on CsHSP90 by blocking ATP binding site at N-terminal domain of predicted structure. Subsequently, GDA mediated leaf disc inhibitor assay further affirmed enhanced HSR with higher expression of CsHSP17.6, CsHSP70, HSP101, and CsHSFA2 genes in tea. Through the current study, efforts were made to extrapolate a deeper understanding of chaperons mediated regulation of HSR attributing thermotolerance in tea.

Genome-wide identification of HSF family in peach and functional analysis of PpHSF5 involvement in root and aerial organ development

Background

Heat shock factors (HSFs) play important roles during normal plant growth and development and when plants respond to diverse stressors. Although most studies have focused on the involvement of HSFs in the response to abiotic stresses, especially in model plants, there is little research on their participation in plant growth and development or on the HSF (PpHSF) gene family in peach (Prunus persica).

Methods

DBD (PF00447), the HSF characteristic domain, was used to search the peach genome and identify PpHSFs. Phylogenetic, multiple alignment and motif analyses were conducted using MEGA 6.0, ClustalW and MEME, respectively. The function of PpHSF5 was confirmed by overexpression of PpHSF5 into Arabidopsis.

Results

Eighteen PpHSF genes were identified within the peach genome. The PpHSF genes were nonuniformly distributed on the peach chromosomes. Seventeen of the PpHSFs (94.4%) contained one or two introns, except PpHSF18, which contained three introns. The in silico-translated PpHSFs were classified into three classes (PpHSFA, PpHSFB and PpHSFC) based on multiple alignment, motif analysis and phylogenetic comparison with HSFs from Arabidopsis thaliana and Oryza sativa. Dispersed gene duplication (DSD at 67%) mainly contributed to HSF gene family expansion in peach. Promoter analysis showed that the most common cis-elements were the MYB (abiotic stress response), ABRE (ABA-responsive) and MYC (dehydration-responsive) elements. Transcript profiling of 18 PpHSFs showed that the expression trend of PpHSF5 was consistent with shoot length changes in the cultivar ‘Zhongyoutao 14’. Further analysis of the PpHSF5 was conducted in 5-year-old peach trees, Nicotiana benthamiana and Arabidopsis thaliana, respectively. Tissue-specific expression analysis showed that PpHSF5 was expressed predominantly in young vegetative organs (leaf and apex). Subcellular localization revealed that PpHSF5 was located in the nucleus in N. benthamiana cells. Two transgenic Arabidopsis lines were obtained that overexpressed PpHSF5. The root length and the number of lateral roots in the transgenic seedlings were significantly less than in WT seedlings and after cultivation for three weeks. The transgenic rosettes were smaller than those of the WT at 2–3 weeks. The two transgenic lines exhibited a dwarf phenotype three weeks after transplanting, although there was no significant difference in the number of internodes. Moreover, the PpHSF5-OE lines exhibited enhanced thermotolerance. These results indicated that PpHSF5 might be act as a suppresser of growth and development of root and aerial organs.

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 the Heat Shock Transcription Factor Gene Family in Brassica juncea: Structure, Evolution, and Expression Profiles

Heat shock transcription factor (HSF) is ubiquitous in the whole biological world and plays an important role in regulating growth and development and responses to environment stress. In this study, a total of 60 HSF transcription factors in Brassica juncea genome were identified and analyzed. Phylogenetic analysis showed that HSF genes were divided into three groups namely: A, B, and C, of which group A was further divided into nine subgroups (A1-A9). The analysis of gene structure and conserved motifs showed that some homologous genes are highly conserved. There was strong conservative microcollinearity among Brassica rapa, B. juncea, and Brassica oleracea, which provides a basis for studying the replication of gene families. Moreover, the results revealed that the promoter regions of BjuHSF genes were rich in cis-elements related to growth and development, hormone signal, and stress response. The prediction of protein interaction results showed that HSFs could interact with multiple transcription factors and proteins in the genome, while functional annotation revealed that BjuHSF genes were involved in many biological processes. The expression patterns of BjuHSF genes were analyzed by qPCR, and the results showed that these genes were closely linked to stress response, hormones, and development process. These results are a foundation for further analysis of the regulation mechanism of HSF gene family.

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.

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.

Alternative splicing in tea plants was extensively triggered by drought, heat and their combined stresses

Drought and heat stresses can influence the expressions of genes, and thereby affect the growth and development of plants. Alternative splicing (AS) of genes plays crucial roles through increasing transcriptome diversity in plant stress responses. Tea plants, widely cultivated in the tropics and subtropics, are often simultaneously exposed to drought and heat stresses. In the present study, we performed a global transcriptome of tea leaves treated with drought, heat or their combination. In total, 19,019, 20,025 and 20,253 genes underwent AS in response to drought (DT), heat (HT) and their combined stress (HD), respectively, of which 12,178, 11,912 and 14,413 genes differentially spliced in response to DT, HT and HD, respectively. Also, 2,447 specific differentially spliced genes (DSGs) were found only in response to HD. All DSGs accounted for  48% of the annotated genes in tea tree genome. Comparison of DSGs and differentially expressive genes (DEGs) showed that the proportions of HT and HD-induced DSGs were 13.4% and 9.2%, while the proportion of DT increased to 28.1%. Moreover, the DEG-DSG overlapped genes tended to be enriched in a wide large of pathways in response to DT. The results indicated that the AS of genes in tea leaves was extensively triggered by drought, heat and their combined stresses. In addition, the AS enhanced the transcriptome adaption in response to drought and heat stresses, and the AS also provoked specific molecular functions in response to drought and heat synergy stress. The study might have practical significance for molecular genetic breeding of tea plants with stress resistance.

Development of a DNA-based real-time PCR assay for the quantification of Colletotrichum camelliae growth in tea (Camellia sinensis)

BACKGROUND

Tea, which is produced from new shoots of existing tea plants (Camellia sinensis), is one of the most popular, non-alcoholic, healthy beverages worldwide. Colletotrichum camelliae is one of the dominant fungal pathogens of tea. The interaction of C. camelliae with tea could be a useful pathosystem to elucidate various aspects of woody, medicinal plant-fungal interactions. Currently, many studies characterizing resistance or virulence and aggressiveness use lesion size at the infection sites on the leaves to quantify the growth of the pathogen. However, this method does not offer the sensitivity needed for the robust quantification of small changes in aggressiveness or the accurate quantification of pathogen growth at the early stages of infection.

RESULTS

A quantitative real-time polymerase chain reaction (qRT-PCR) assay was developed for the quantification of C. camelliae growth on tea plant. This method was based on the comparison of fungal DNA in relation to plant biomass. This assay was used to investigate the phenotypes of tea plant cultivars in response to C. camelliae infection. Two cultivars, Zhongcha 108 (ZC108) and Longjing 43 (LJ43), were tested with this method. ZC108 was previously reported as an anthracnose-resistant cultivar against C. camelliae, while LJ43 was susceptible. The traditional lesion measurement method showed that both cultivars were susceptible to a virulent strain of C. camelliae, while the qRT-PCR approach indicated that very little fungal growth occurred in the anthracnose-resistant cultivar ZC108. The observed results in this study were consistent with previously published research. In addition, the DNA-based real-time PCR method was applied for analysis of pathogenic differences in general C. camelliae isolates and among several Colletotrichum spp that infect tea.

CONCLUSIONS

This study showed that the DNA-based qRT-PCR technique is rapid, highly sensitive and easily applicable for routine experiments and could be used in screening for resistant tea plant cultivars or to identify differences in pathogen aggressiveness within and among Colletotrichum species.

Genome-Wide Characterization of the C-repeat Binding Factor (CBF) Gene Family Involved in the Response to Abiotic Stresses in Tea Plant (Camellia sinensis)

C-repeat (CRT)/dehydration responsive element (DRE)-binding factor CBFs, a small family of genes encoding transcriptional activators, play important roles in plant cold tolerance. In this study, a comprehensive genome-wide analysis was carried out to identify and characterize the functional dynamics of CsCBFs in tea plant (Camellia sinensis). A total of 6 CBF genes were obtained from the tea plant genome and named CBF1-6. All of the CsCBFs had an AP2/ERF DNA-binding domain and nuclear localization signal (NLS) sequence. CsCBF-eGFP fusion and DAPI staining analysis confirmed the nuclear localization of the CsCBFs. Transactivation assays showed that the CsCBFs, except CsCBF1, had transcriptional activity. CsCBF expression was differentially induced by cold, heat, PEG, salinity, ABA, GA, MeJA, and SA stresses. In particular, the CsCBF genes were significantly induced by cold treatments. To further characterize the functions of CsCBF genes, we overexpressed the CsCBF3 gene in Arabidopsis thaliana plants. The resulting transgenic plants showed increased cold tolerance compared with the wild-type Arabidopsis plant. The enhanced cold tolerance of the transgenic plants was potentially achieved through an ABA-independent pathway. This study will help to increase our understanding of CsCBF genes and their contributions to stress tolerance in tea plants.

Transcriptomic and metabolomic profiling of Camellia sinensis L. cv. 'Suchazao' exposed to temperature stresses reveals modification in protein synthesis and photosynthetic and anthocyanin biosynthetic pathways

To determine the mechanisms in tea plants responding to temperature stresses (heat and cold), we examined the global transcriptomic and metabolomic profiles of the tea plant cultivar 'Suchazao' under moderately low temperature stress (ML), severely low temperature stress (SL), moderately high temperature stress (MH) and severely high temperature stress (SH) using RNA-seq and high performance liquid chromatography tandem mass spectrometry/mass spectrometry (HPLC-MS/MS), respectively. The identified differentially expressed genes indicated that the synthesis of stress-resistance protein might be redirected to cope with the temperature stresses. We found that heat shock protein genes Hsp90 and Hsp70 played more critical roles in tea plants in adapting to thermal stress than cold, while late embryogenesis abundant protein genes (LEA) played a greater role under cold than heat stress, more types of zinc finger genes were induced under cold stress as well. In addition, energy metabolisms were inhibited by SH, SL and ML. Furthermore, the mechanisms of anthocyanin synthesis were different under the cold and heat stresses. Indeed, the CsUGT75C1 gene, encoding UDP-glucose:anthocyanin 5-O-glucosyl transferase, was up-regulated in the SL-treated leaves but down-regulated in SH. Metabolomics analysis also showed that anthocyanin monomer levels increased under SL. These results indicate that the tea plants share certain foundational mechanisms to adjust to both cold and heat stresses. They also developed some specific mechanisms for surviving the cold or heat stresses. Our study provides effective information about the different mechanisms tea plants employ in surviving cold and heat stresses, as well as the different mechanisms of anthocyanin synthesis, which could speed up the genetic breeding of heat- and cold-tolerant tea varieties.

Comparative transcriptome analysis to elucidate the enhanced thermotolerance of tea plants (Camellia sinensis) treated with exogenous calcium

The molecular mechanisms regulating calcium-mediated thermotolerance in Camellia sinensis were revealed by RNA-Sequencing. Heat stress is one of the most remarkable abiotic factors limiting the growth and productivity of Camellia sinensis plants. Calcium helps regulate plant responses to various adverse environmental conditions, including heat stress. In this study, the effects of exogenous calcium on the physiological characteristics of heat-stressed C. sinensis were investigated. A calcium pretreatment increased the proline, soluble sugar, Ca²⁺, and chlorophyll contents, but decreased the malondialdehyde content and relative electrical conductivity in C. sinensis leaves under heat stress. Further analysis of the ultra-structure of chloroplasts indicated that heat stress induced accumulation of starch granules and destruction of the stroma lamella in C. sinensis. However, calcium pretreatment counteracted the adverse effects of heat stress on the structure of the photosynthetic apparatus. These results imply that the calcium pretreatment increased C. sinensis thermotolerance. Moreover, RNA-sequencing was applied to characterize the calcium-mediated transcript-level responses to heat stress. A total of 923 differentially expressed genes (DEGs) including 299 up-regulated and 624 down-regulated genes were identified. Functional annotations indicated that these DEGs were primarily related to signal transduction, transcriptional regulation, and post-translational modification. In addition, a C. sinensis gene [CsCML45 (GenBank: KY652927)] encoding a calmodulin-like protein was isolated. The heterologous expression of CsCML45 enhanced the thermotolerance of transgenic Arabidopsis thaliana plants. These results may be useful for characterizing the calcium-mediated molecular mechanism responsible for C. sinensis thermotolerance.

Molecular Characterization and Expression Profile Analysis of Heat Shock Transcription Factors in Mungbean

Heat shock transcription factors (Hsfs) are essential elements in plant signal transduction pathways that mediate gene expression in response to various abiotic stresses. Mungbean (Vigna radiata) is an important crop worldwide. The emergence of a genome database now allows for functional analysis of mungbean genes. In this study, we dissect the mungbean Hsfs using genome-wide identification and expression profiles. We characterized a total of 24 VrHsf genes and classified them into three groups (A, B, and C) based on their phylogeny and conserved domain structures. All VrHsf genes exhibit highly conserved exon-intron organization, with two exons and one intron. In addition, all VrHsf proteins contain 16 distinct motifs. Chromosome location analysis revealed that VrHsf genes are located on 8 of the 11 mungbean chromosomes, and that seven duplicated gene pairs had formed among them. Moreover, transcription patterns of VrHsf genes varied in different tissues, indicating their different roles in plant growth and development. We identified multiple stress related cis-elements in VrHsf promoter regions 2 kb upstream of the translation initiation codons, and the expression of most VrHsf genes was altered under different stress conditions, suggesting their potential functions in stress resistance pathways. These molecular characterization and expression profile analyses of VrHsf genes provide essential information for further function investigation.

Heat shock transcription factor 3 regulates plant immune response through modulation of salicylic acid accumulation and signalling in cassava

As the terminal components of signal transduction, heat stress transcription factors (Hsfs) mediate the activation of multiple genes responsive to various stresses. However, the information and functional analysis are very limited in non-model plants, especially in cassava (Manihot esculenta), one of the most important crops in tropical areas. In this study, 32 MeHsfs were identified from the cassava genome; the evolutionary tree, gene structures and motifs were also analysed. Gene expression analysis found that MeHsfs were commonly regulated by Xanthomonas axonopodis pv. manihotis (Xam). Amongst these MeHsfs, MeHsf3 was specifically located in the cell nucleus and showed transcriptionally activated activity on heat stress elements (HSEs). Through transient expression in Nicotiana benthamiana leaves and virus-induced gene silencing (VIGS) in cassava, we identified the essential role of MeHsf3 in plant disease resistance, by regulating the transcripts of Enhanced Disease Susceptibility 1 (EDS1) and pathogen-related gene 4 (PR4). Notably, as regulators of defence susceptibility, MeEDS1 and MePR4 were identified as direct targets of MeHsf3. Moreover, the disease sensitivity of MeHsf3- and MeEDS1-silenced plants could be restored by exogenous salicylic acid (SA) treatment. Taken together, this study highlights the involvement of MeHsf3 in defence resistance through the transcriptional activation of MeEDS1 and MePR4.

Genomic and transcriptomic analyses of HD-Zip family transcription factors and their responses to abiotic stress in tea plant (Camellia sinensis)

Tea plant (Camellia sinensis (L.) O. Kuntze) is a perennial evergreen woody plant, and its leaves contain various beneficial ingredients and have healthy efficacy. HD-Zip (homeodomain-leucine zipper) transcription factors (TFs) are widely distributed in plants and play an important role in plant growth and environmental response. To date, knowledge on HD-Zip gene family in tea plant is still limited. In this study, 33 HD-Zip TFs were selected based on the genomic and transcriptomic databases of tea plant. The conserved domains and common motifs of these TFs were predicted and analyzed. These 33 Cshdz TFs were divided into four groups (HD-Zip I, HD-Zip II, HD-Zip III, and HD-Zip IV). The interaction network of the HD-Zip proteins of tea plant was established based on the data of Arabidopsis. In addition, the expression levels of these Cshdz genes in tea plant cv. 'Longjing43' were detected and analyzed under five abiotic stress treatments. Results showed that the different expression profiles of Cshdz genes were associated with different abiotic stress treatments. Our findings suggested a potential relationship between the resistance of tea plant and its Cshdz genes.

Transcriptome and Expression Profiling Analysis of Recalcitrant Tea (Camellia sinensis L.) Seeds Sensitive to Dehydration

The tea plant (Camellia sinensis (L.) O. Kuntze) is an economically important woody perennial nonalcoholic health beverage crop. Tea seeds are categorized as recalcitrant and are sensitive to dehydration treatment. However, the molecular basis of this phenomenon has not been investigated. Thus, we analyzed the genome-wide expression profiles of three dehydration stages using RNA-Seq and digital gene expression (DGE) technologies. We performed de novo assembly and obtained a total of 91,925 nonredundant unigenes, of which 58,472 were extensively annotated. By a hierarchical clustering of differentially expressed genes (DEGs), we found that 8929 DEGs were downregulated and 5875 DEGs were upregulated during dehydration treatment. A series of genes related to ABA biosynthesis and signal transduction, transcription factor, antioxidant enzyme, LEA protein, and proline metabolism that have been reported to function in dehydration process were found to be downregulated. Additionally, the expression profiles of 12 selected genes related to tea seed dehydration treatment were confirmed by qRT-PCR analysis. To our knowledge, this is the first genome-wide study elucidating the possible molecular mechanisms of sensitivity of recalcitrant tea seeds to dehydration. The results obtained in this study contribute to the preservation of tea seeds as genetic resources and can also be used to explore the mechanism of dehydration sensitivity of other recalcitrant seeds.

Transcriptome-wide identification and expression profile analysis of the bHLH family genes in Camellia sinensis

The tea plant is an important commercial horticulture crop cultivated worldwide. Yield and quality of this plant are influenced by abiotic stress. The bHLH family transcription factors play a pivotal role in the growth and development, including abiotic stress response, of plants. A growing number of bHLH proteins have been functionally characterized in plants. However, few studies have focused on the bHLH proteins in tea plants. In this study, 120 CsbHLH TFs were identified from tea plants using computational prediction method. Structural analysis detected 23 conservative residues, with over 50% identities in the bHLH domain. Moreover, 103 CsbHLH proteins were assumed to bind DNA and encompassed 98 E-Box binders and 85 G-Box binders. The CsbHLH proteins were grouped into 20 subfamilies based on phylogenetic analysis and a previous classification system. A survey of transcriptome profiling screened 22 and 39 CsbHLH genes that were upregulated under heat and drought stress. Nine CsbHLH genes were validated using qRT-PCR. Results were approximately in accordance with transcriptome data. These genes could be induced by one or more abiotic stresses.

Reliable reference genes for normalization of gene expression data in tea plants (Camellia sinensis) exposed to metal stresses

Tea plants [Camellia sinensis (L.) O. Kuntze] are an important leaf-type crop that are widely used for the production of non-alcoholic beverages in the world. Exposure to excessive amounts of heavy metals adversely affects the quality and yield of tea leaves. To analyze the molecular responses of tea plants to heavy metals, a reliable quantification of gene expression is important and of major importance herein is the normalization of the measured expression levels for the target genes. Ideally, stably expressed reference genes should be evaluated in all experimental systems. In this study, 12 candidate reference genes (i.e., 18S rRNA, Actin, CYP, EF-1α, eIF-4α, GAPDH, MON1, PP2AA3, TBP, TIP41, TUA, and UBC) were cloned from tea plants, and the stability of their expression was examined systematically in 60 samples exposed to diverse heavy metals (i.e., manganese, aluminum, copper, iron, and zinc). Three Excel-based algorithms (geNorm, NormFinder, and BestKeeper) were used to evaluate the expression stability of these genes. PP2AA3 and 18S rRNA were the most stably expressed genes, even though their expression profiles exhibited some variability. Moreover, commonly used reference genes (i.e., GAPDH and TBP) were the least appropriate reference genes for most samples. To further validate the suitability of the analyzed reference genes, the expression level of a phytochelatin synthase gene (i.e., CsPCS1) was determined using the putative reference genes for data normalizations. Our results may be beneficial for future studies involving the quantification of relative gene expression levels in tea plants.

The CsHSP17.2 molecular chaperone is essential for thermotolerance in Camellia sinensis

Small heat shock proteins (sHSPs) play important roles in responses to heat stress. However, the functions of sHSPs in tea plants (Camellia sinensis) remain uncharacterized. A novel sHSP gene, designated CsHSP17.2, was isolated from tea plants. Subcellular localization analyses indicated that the CsHSP17.2 protein was present in the cytosol and the nucleus. CsHSP17.2 expression was significantly up-regulated by heat stress but was unaffected by low temperature. The CsHSP17.2 transcript levels increased following salt and polyethylene glycol 6000 treatments but decreased in the presence of abscisic acid. The molecular chaperone activity of CsHSP17.2 was demonstrated in vitro. Transgenic Escherichia coli and Pichia pastoris expressing CsHSP17.2 exhibited enhanced thermotolerance. The transgenic Arabidopsis thaliana exhibited higher maximum photochemical efficiencies, greater soluble protein proline contents, higher germination rates and higher hypocotyl elongation length than the wild-type controls. The expression levels of several HS-responsive genes increased in transgenic A. thaliana plants. Additionally, the CsHSP17.2 promoter is highly responsive to high-temperature stress in A. thaliana. Our results suggest that CsHSP17.2 may act as a molecular chaperone to mediate heat tolerance by maintaining maximum photochemical efficiency and protein synthesis, enhancing the scavenging of reactive oxygen species and inducing the expression of HS-responsive genes.

Heat shock transcription factors in banana: genome-wide characterization and expression profile analysis during development and stress response

Banana (Musa acuminata) is one of the most popular fresh fruits. However, the rapid spread of fungal pathogen Fusarium oxysporum f. sp. cubense (Foc) in tropical areas severely affected banana growth and production. Thus, it is very important to identify candidate genes involved in banana response to abiotic stress and pathogen infection, as well as the molecular mechanism and possible utilization for genetic breeding. Heat stress transcription factors (Hsfs) are widely known for their common involvement in various abiotic stresses and plant-pathogen interaction. However, no MaHsf has been identified in banana, as well as its possible role. In this study, genome-wide identification and further analyses of evolution, gene structure and conserved motifs showed closer relationship of them in every subgroup. The comprehensive expression profiles of MaHsfs revealed the tissue- and developmental stage-specific or dependent, as well as abiotic and biotic stress-responsive expressions of them. The common regulation of several MaHsfs by abiotic and biotic stress indicated the possible roles of them in plant stress responses. Taken together, this study extended our understanding of MaHsf gene family and identified some candidate MaHsfs with specific expression profiles, which may be used as potential candidates for genetic breeding in banana.

Transcriptome-Wide Identification and Expression Analysis of the NAC Gene Family in Tea Plant [Camellia sinensis (L.) O. Kuntze]

In plants, the NAC (NAM-ATAF1/2-CUC) family of proteins constitutes several transcription factors and plays vital roles in diverse biological processes, such as growth, development, and adaption to adverse factors. Tea, as a non-alcoholic drink, is known for its bioactive ingredients and health efficacy. Currently, knowledge about NAC gene family in tea plant remains very limited. In this study, a total of 45 CsNAC genes encoding NAC proteins including three membrane-bound members were identified in tea plant through transcriptome analysis. CsNAC factors and Arabidopsis counterparts were clustered into 17 subgroups after phylogenetic analysis. Conserved motif analysis revealed that CsNAC proteins with a close evolutionary relationship possessed uniform or similar motif compositions. The distribution of NAC family MTFs (membrane-associated transcription factors) among higher plants of whose genome-wide has been completed revealed that the existence of doubled TMs (transmembrane motifs) may be specific to fabids. Transcriptome analysis exhibited the expression profiles of CsNAC genes in different tea plant cultivars under non-stress conditions. Nine CsNAC genes, including the predicted stress-related and membrane-bound genes, were examined through qRT-PCR (quantitative real time polymerase chain reaction) in two tea plant cultivars, namely, 'Huangjinya' and 'Yingshuang'. The expression patterns of these genes were investigated in different tissues (root, stem, mature leaf, young leaf and bud) and under diverse environmental stresses (drought, salt, heat, cold and abscisic acid). Several CsNAC genes, including CsNAC17 and CsNAC30 that are highly orthologous to known stress-responsive ANAC072/RD26 were identified as highly responsive to abiotic stress. This study provides a global survey of tea plant NAC proteins, and would be helpful for the improvement of stress resistance in tea plant via genetic engineering.

Transcriptome-Based Analysis of Dof Family Transcription Factors and Their Responses to Abiotic Stress in Tea Plant (Camellia sinensis)

Tea plant (Camellia sinensis (L.) O. Kuntze) is affected by abiotic stress during its growth and development. DNA-binding with one finger (Dof) transcription factors (TFs) play important roles in abiotic stress tolerance of plants. In this study, a total of 29 putative Dof TFs were identified based on transcriptome of tea plant, and the conserved domains and common motifs of these CsDof TFs were predicted and analyzed. The 29 CsDof proteins were divided into 7 groups (A, B1, B2, C1, C2.1, C2.2, and D2), and the interaction networks of Dof proteins in C. sinensis were established according to the data in Arabidopsis. Gene expression was analyzed in “Yingshuang” and “Huangjinya” under four experimental stresses by qRT-PCR. CsDof genes were expressed differentially and related to different abiotic stress conditions. In total, our results might suggest that there is a potential relationship between CsDof factors and tea plant stress resistance.

A Comparative Analysis of Machine Learning with WorldView-2 Pan-Sharpened Imagery for Tea Crop Mapping

Tea is an important but vulnerable economic crop in East Asia, highly impacted by climate change. This study attempts to interpret tea land use/land cover (LULC) using very high resolution WorldView-2 imagery of central Taiwan with both pixel and object-based approaches. A total of 80 variables derived from each WorldView-2 band with pan-sharpening, standardization, principal components and gray level co-occurrence matrix (GLCM) texture indices transformation, were set as the input variables. For pixel-based image analysis (PBIA), 34 variables were selected, including seven principal components, 21 GLCM texture indices and six original WorldView-2 bands. Results showed that support vector machine (SVM) had the highest tea crop classification accuracy (OA = 84.70% and KIA = 0.690), followed by random forest (RF), maximum likelihood algorithm (ML), and logistic regression analysis (LR). However, the ML classifier achieved the highest classification accuracy (OA = 96.04% and KIA = 0.887) in object-based image analysis (OBIA) using only six variables. The contribution of this study is to create a new framework for accurately identifying tea crops in a subtropical region with real-time high-resolution WorldView-2 imagery without field survey, which could further aid agriculture land management and a sustainable agricultural product supply.

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.