Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 22 July 2024
Sec. Plant Abiotic Stress

Integrated transcriptomics and metabolomics provides insights into the Nicotiana tabacum response to heat stress

Hao ChenHao Chen1Shaoxin Qiu,Shaoxin Qiu1,2Yuanping ChenYuanping Chen3Jiqin LiJiqin Li1Tingyu XuTingyu Xu1Pingzhan ZhongPingzhan Zhong3Xiuhong ShaoXiuhong Shao1Shihuan XuShihuan Xu3Zhuwen MaZhuwen Ma1Zhenrui HuangZhenrui Huang1Xiaoying Pan*Xiaoying Pan1*
  • 1Guangdong Key Laboratory for Crops Genetic Improvement, Crops Research Institute, Guangdong Academy of Agricultural Sciences (GAAS), Guangdong Provincial Engineering & Technology Research Center for Tobacco Breeding and Comprehensive Utilization, Guangzhou, China
  • 2College of Agronomy, South China Agricultural University, Guangzhou, China
  • 3China National Tobacco Corporation, Guangdong Company, Guangzhou, China

Heat stress is a prevalent factor that significantly damages crops, especially with the ongoing global warming and increasing frequency of extreme weather events. Tobacco is particularly sensitive to temperature fluctuations, experiencing reduced yield and quality under high temperatures. However, the underlying molecular mechanisms of heat resistance in tobacco remain poorly understood. This study comprehensively analyzed biochemical, transcriptomic, and metabolomic responses to heat stress on the root and shoot of the tobacco cultivar K326 compared to control conditions. Heat stress significantly increased the activities of antioxidant enzymes (CAT, POD, and SOD) and levels of osmotic mediators (soluble sugars, sucrose, and proline) in the shoot. Furthermore, transcriptome analysis identified 13,176 differentially expressed genes (DEGs) in the root (6,129 up-regulated and 7,047 down-regulated) and 12,283 DEGs (6,621 up-regulated and 5,662 down-regulated) in the shoot. The root had 24 enriched KEGG pathways, including phenylpropanoid metabolism, while the shoot had 32 significant pathways, such as galactose metabolism and MAPK signaling. The metabolomic data identified 647 metabolites in the root and 932 in the shoot, with carbohydrates and amino acids being the main categories. The root had 116 differentially abundant metabolites (DAMs) (107 up-regulated and 9 down-regulated), and the shoot contained 256 DAMs (251 up-regulated and 5 down-regulated). Joint transcriptome and metabolome analysis showed that galactose metabolism and starch and sucrose metabolism were co-enriched in both tissues. In contrast, amino sugar and nucleotide sugar metabolism was enriched in the root, and purine metabolism in the shoot. The purine metabolic pathway in the shoot can modulate the expression of MYB transcription factors by influencing ABA synthesis and signaling, thereby controlling the accumulation of HSPs, raffinose, sucrose, and trehalose to enhance heat tolerance. Furthermore, NtMYB78, an MYB transcription factor, enhances tolerance for heat stress in tobacco. This research offers a foundational framework for investigating and implementing heat-resistant genes and metabolic pathways in the root and shoot of tobacco seedlings.

Introduction

The escalating phenomenon of global warming has led to exceptionally high temperatures, significantly impacting crop productivity and diversity. Research suggests that every 1-degree Celsius increase in the global mean temperature could reduce the yields of wheat, rice, maize, and soybean by 6.0, 3.2, 7.4, and 3.1%, respectively, substantially threatening global food security (Shekhawat et al., 2022). High temperatures significantly reduce the yield and quality of tobacco, a plant species sensitive to temperature fluctuations (Yang et al., 2018). As a vital economic crop contributing to national tax revenue, the cultivation and utilization of heat-resistant tobacco varieties is paramount (Lu et al., 2011). Therefore, investigating the molecular mechanisms underlying heat resistance in tobacco is an important foundational framework for future breeding efforts.

The cell membrane system is the primary defense mechanism against heat stress and facilitates the transmission of heat-responsive signals. Thus, elevated temperatures can impair cell membrane integrity, heightening membrane permeability, compromising thermal stability and osmoregulation, and ultimately causing cellular demise (Mathur et al., 2014; de Pinto et al., 2015; Fortunato et al., 2023; Saeed et al., 2023). High temperatures cause an imbalance in reactive oxygen species in plant cells, leading to the accumulation of harmful substances like hydrogen peroxide and malondialdehyde, putting plants at risk of oxidative damage (Xie et al., 2019; Zhao et al., 2021; Sato et al., 2024). Plants have developed intricate heat stress response mechanisms, such as pathways for temperature sensing, signaling, gene expression regulation, protein synthesis, and metabolic regulation, in order to alleviate the detrimental effects of elevated temperatures (Ohama et al., 2017; Guihur et al., 2022; Shekhawat et al., 2022; Kan et al., 2023). When plants are subjected to high temperature stress, they detect stimuli through receptors located on the cell membrane and and subsequently relay these signals to the nucleus. Signaling molecules such as Ca2+ and MAPKs (mitogen-activated protein kinases) play key roles in heat stress signaling (Niu et al., 2020; Haider et al., 2021). The phytohormone ABA is important in response to heat stress. Overexpressing the ABA-responsive element binding protein (AREB) can improve heat tolerance in Arabidopsis (Suzuki et al., 2016). Moreover, ABA regulates ZmCDPK to improve ROS clearing and adapt to high-temperature environments (Zhao et al., 2021). Furthermore, MYB transcription factors are significant in response to stress, particularly abiotic stresses such as drought, temperature fluctuations, and salinity (Li et al., 2015; Jacob et al., 2021; Wang et al., 2021b). Heat stress transcription factors (HSFs) are essential regulators in plants that are activated in response to heat stress, binding to heat stress elements (HSEs) to stimulate the expression of heat stress proteins (HSPs) (Andrási et al., 2021; Tan et al., 2023; Wen et al., 2023). HSPs, a group of molecular chaperones, facilitate proper protein folding and prevent aggregation, playing a vital role in plant defense mechanisms against heat stress (Jacob et al., 2017; Song et al., 2021; Mondal et al., 2023).

Plant sugars serve multiple functions as energy sources, antioxidants, and essential components of cellular structures. Additionally, these sugars are crucial in regulating plant growth and development and aiding in plant adaptation to challenging environments (Chng et al., 2014; Proels and Hückelhoven, 2014). Monosaccharides such as hexose glucose, fructose, and certain glucose derivatives and disaccharides sucrose and trehalose-6-phosphate regulate various biological processes (Halford et al., 2003; Nunes et al., 2013; Tognetti et al., 2013). For instance, salinity and drought stress elevate the sucrose levels in the phloem sap of Arabidopsis, thereby maintaining the water potential (Durand et al., 2016). Trehalose interacts with various sugars, osmoprotectants, amino acids, and phytohormones to modulate metabolic reprogramming, which is essential for adaptation to heat stress (Luo et al., 2018; Zhang et al., 2022; Raza et al., 2023).

Thus, analyzing the transcriptomes and metabolomes can unravel key plant metabolic pathways and regulatory networks, improving the understanding of stress resistance in crops as sequencing costs decrease (Wang et al., 2019; Zhao et al., 2019). Many studies have used single or combined transcriptome and metabolome analyses to study cold and drought tolerance in tobacco (Jin et al., 2017; Luo et al., 2022b; Hu et al., 2023), and few on heat stress. Therefore, this study utilized an integrative analysis to investigate the transcriptome and metabolome of the root and shoot tissues of K326 tobacco. The main co-enriched pathways in both tissues were starch and sucrose metabolism and galactose metabolism. The amino acid sugar metabolism pathway was specifically and primarily enriched in the root, and purine metabolism was enriched in the shoot. Heat stress up-regulated key genes involved in sugar metabolism, including raffinose synthase (RS), sucrose synthase (SUS), trehalose 6-phosphate synthase (TPS), and trehalose 6-phosphate phosphatase (TPP), significantly increasing the levels of raffinose, sucrose, and trehalose. The MYB transcription factor, NtMYB78, discovered in this study, plays a crucial role in enhancing the heat stress response in tobacco plants. Our research can serve as a foundation for identifying key genes involved in the heat stress response in tobacco.

Materials and methods

Plant growth and heat treatment

The seeds of cultivated tobacco (Nicotiana tabacum L.) K326 (used in this study) were treated with 75% alcohol at room temperature for 5 minutes and washed five times with sterile water. Subsequently, the seeds were planted in small plugs and allowed to grow until they developed two small leaves before transplantation. At the five-leaf stage, eight similar-sized tobacco plants were heat treated at 45°C, 50% relative humidity, 16 h light/8 h dark photoperiod, and 10000 lux light intensity for 24 h. The control group was subjected to identical light and humidity conditions at 25°C. The root sample consists of all the root systems, while the leaves and stems of each plant were combined into one shoot sample. The root and shoot were collected for transcriptome sequencing.

Physiological and biochemical assays

Relative conductivity assays

Here, the leaves were washed with distilled water and dried, cut into long strips (avoiding the main veins). Next, three fresh samples (0.1 g each) were placed in a centrifuge tube containing 10 mL of deionized water, and their conductivity before boiling (R1 and R2) was determined using a DDS-307 conductivity meter. The relative conductivity was calculated as R1/R2×100%.

Enzyme activity assays

Here, 0.2 g of fresh leaves were ground into powder using liquid nitrogen and 2 mL pre-cooled PBS (50 mM pH 7.8). The mixture was homogenized in an ice bath, transferred to a centrifuge tube, washed with PBS, and centrifuged at 8000 rpm for 15 min at 4°C. The supernatant was identified as the crude enzyme solution, and the levels of malondialdehyde (MDA) and enzymatic activities (catalase (CAT), superoxide dismutase (SOD), and peroxidase (POD)) were assessed following prior protocols (Hu et al., 2016).

Total soluble sugar and sucrose assays

Soluble sugars were quantified using anthrone colorimetry as outlined in a previous study (Abdel-Basset et al., 2010). The sucrose content was quantified using the plant sucrose content assay kit (Shanghai Acmec Biochemical Co. Ltd, Shanghai, China), following the guidelines outlined in the instruction manual.

Proline assays

Here, 0.5 g of leaves were ground and placed in a test tube containing 5 mL of 3% sulfosalicylic acid solution. The mixture was heated in a boiling water bath for 15 minutes, filtered, and filtrate collected. Next, 2 mL of the extract was combined with equal amounts of acetic acid and ninhydrin in a new test tube and sealed. The sealed test tube was heated in boiling water for 15 minutes, then cooled, and the solution was mixed with 5 mL toluene for extraction. The toluene layer was transferred to a cuvette, and its absorbance was measured at 520 nm using a spectrophotometer (Rajametov et al., 2021).

RNA-seq library preparation and data analysis

Total RNA was isolated from tobacco root and shoot tissues using Trizol Reagent (Thermo Fisher Scientific, MA, USA). The RNA quality was evaluated using the Agilent 2100 Bioanalyzer (Agilent Technologies, CA, USA). Approximately 1 µg of RNA was utilized as input material for the transcriptome library preparation per sample. The libraries underwent sequencing on an Illumina NovaSeq platform (Tianjin, China), resulting in the generation of 150-bp paired-end reads.

The raw data was quality controlled using the FastQC software, eliminating reads containing adapter sequences, N (representing undetermined bases), and those with low-quality (> 50% of bases having a quality score < 5). Additionally, the Q20, Q30, and GC content were calculated for the clean data. Next, HISAT2 (v2.0.5), the reference genome index, was constructed with Nitab v1.0 Chr Edwards 2017 as the reference genome. The paired-end clean reads were aligned to the reference genome using HISAT2 v2.0.5, and FeatureCounts (v1.6.4) was employed to count the reads mapped to each gene (Liao et al., 2014; Kim et al., 2019). Subsequently, the FPKM (Fragments Per Kilobase of transcript per Million mapped reads) for each gene was computed based on gene length. The Pearson correlation coefficients among various samples were computed utilizing the cor function within the R package, and subsequently visualized as a heatmap through the pheatmap software. DESeq2 software (1.20.0) was utilized for differential gene expression analysis between the two comparison groups. Genes with expression levels of |log2 (fold change)| ≥ 1 and p-value ≤ 0.05 were retained for further analyses. The DEGs were further subjected to GO and KEGG enrichment analysis using clusterProfiler software (3.8.1) (Yu et al., 2012).

Metabolite extraction and data analysis

A 100 mg tissue sample was flash-frozen in liquid nitrogen, transferred to an Eppendorf tube, and mixed with 500 μL of 80% methanol solution. The samples were then subjected to a 5-minute incubation on ice, followed by centrifugation at 15,000 g and 4°C for 20 minutes. One volume of the supernatant was extracted and diluted with mass spectrometry-grade water to a final methanol content of 53%. The specimen was centrifuged again at 15,000 g and 4°C for 20 minutes, then the supernatant was collected for subsequent analysis via liquid chromatography-mass spectrometry (LC-MS/MS).

The targeted metabolomics investigation (Dunn et al., 2011) was performed on the QTRAP® 6500+ mass spectrometry platform (SCIEX, MA, USA). The detailed parameter settings of the instrument were referred to previous studies (Wang et al., 2021a).

The identified metabolites were characterized on the KEGG (https://www.genome.jp/kegg/pathway.html), HMDB (https://hmdb.ca/metabolites), and LIPIDMaps (http://www.lipidmaps.org/) databases. Principal components analysis (PCA) was conducted using the metaX software. T-test was used to determine the statistical significance (P-value) and fold change (FC) of metabolites between two groups. Metabolites meeting the criteria of a < 0.05 p-value, > 1 VIP (Variable Importance in Projection), and |log2fold change|≥1 were considered as differentially expressed. A volcano plot displaying the differentially abundant metabolites (DAMs) was generated using the ggplot2 package in R. Subsequently, Pearson correlation analysis was conducted on the DAMs using the cor () function in R. The statistical significance was determined using the cor.mtest () function in R, with a p-value of < 0.05 indicating statistical significance. The corrplot package in R was employed to visualize the correlation matrix, and a bubble plot was generated using the ggplot2 package in R. The functions and metabolic pathways of the metabolites were explored using the KEGG database. Pathway enrichment was determined by comparing the x/n > y/n ratio, and pathways with p-value < 0.05 were considered as significantly enriched.

Mapping the network of correlations between DEGs and DAMs

Pearson’s correlation coefficients were calculated using Cor in R (Version 4.2.3) to assess the relationship between genes and metabolites within each common pathway. Next, a network diagram was created using Cytoscape (Version 3.10.1) to show the metabolite-gene using Pearson’s correlation coefficients above 0.9 in all groups.

Quantitative real-time polymerase chain reaction

High-quality RNA was extracted and reverse transcribed with the HiScript II 1st Strand cDNA Synthesis Kit (Vazyme, Nanjing, China). The gene-specific primers were designed using the Primer5.0 software, and qRT-PCR was performed using the CFX96 real-time PCR detection system (Bio-Rad, CA, USA). The 10 μL reaction included 1 μL cDNA template, 0.25 μL of each primer, 5 μL SYBR Green Mix, and 3.5 μL ddH2O. The thermal cycling program was performed as follows: 95°C for 30 s, 40 cycles of 95°C for 5 s, and 60°C for 30 s. Relative gene expression levels were calculated using the cycling threshold (Ct) 2–ΔΔCt method with NtACTIN as the internal reference gene. Each group had three biological and three technical replicates. The information of genes for qRT-PCR analysis is listed in Supplementary Table 1.

Virus-induced gene silencing of the target genes

Based on the online website VIGS tool (https://vigs.solgenomics.net/), specific primers were designed to amplify a 300 bp silencing fragment using the tobacco seedling root cDNA as a template, and the amplified fragment was ligated into the pCE2 TA/Blunt-Zero Vector (Vazyme Biotech, Nanjing, China) to obtain a positive clone by sequencing. The designed homology arm was cloned into the linearized pTRV2 vector using the homologous recombination method using EcoRI and KpnI restriction endonucleases, and the positive plasmid was obtained by sequencing. The pTRV1, pTRV2-empty vector, pTRV2-PDS (positive control), and pTRV2-NtMYB78 plasmids were transformed into Agrobacterium tumefaciens strain GV3010. Next, the positive single clones were picked out and cultured in 500 μL LB liquid medium containing 25 μg/mL rifampicin and 50 μg/mL kanamycin, with constant shaking at 200 rpm, 28°C overnight. The constructs were cultivated further for 16 to 24 hours at 28°C in 50 mL of the same LB medium with constant shaking until their OD600 reached 0.6-0.8. The mixture was centrifuged at 3000×g for 15 min, the supernatant was removed, and the OD was adjusted to 0.8-1 using the infection buffer (10 mM MES; 10 mM MgCl2; 200 μM AS). An Agrobacterium solution with equal concentration of pTRV1 and pTRV2-empty vector, pTRV2-PDS, and pTRV2-NtMYB78 were separately mixed and left at room temperature for 2 hours. Syringes were used to inject 5-8 uniformly-growing 4-week-old tobacco (Nicotiana benthamiana) plants in each group. The seedlings were maintained at 18°C and darkness for 2 d. Then, they were cultured at 22°C/18°C, 16h day/8h night for approximately 10 d, when the pTRV2-PDS plants showed an albino phenotype. Samples were collected for RNA extraction and analyzed for NtMYB78 expression. The information of primers was listed in Supplementary Table 1.

Results

Heat stress significantly altered the phenotype and physiological and biochemical indexes of K326 seedlings

After 24h of heat stress treatment at 45°C, the five-leaf stage K326 tobacco plants showed severe leaf wilting (Figure 1A). The notable rise in leaf relative conductivity following exposure to heat stress suggests a substantial impairment of leaf cell membrane integrity (Figure 1B). Moreover, heat treatment significantly elevated the content of MDA, the end product of lipid peroxidation in plant membranes and an important response parameter to the antioxidant capacity of the plant (Wang et al., 2016). Heat stress significantly increased the levels of antioxidant enzymes CAT, POD, and SOD, as did the concentrations of soluble sugar, sucrose, and proline (Figure 1B).

Figure 1
www.frontiersin.org

Figure 1 The phenotypic, physiological, and biochemical indicators of heat stress in tobacco seedlings. (A) The control group (25°C for 24h (CK left)), heat stress group (45°C for 24h. (HS, right)), Bar=10 cm. (B) Leaf relative conductivity, MDA content, CAT, POD, and SOD enzyme activity, soluble sugar content, sucrose content, and proline contents. Significance was assessed using more than four biological replicates, with asterisks denoting statistically significant differences between treatments as determined by the paired Student’s t-test (***P < 0.001).

Transcriptome analysis of tobacco seedlings in response to heat stress

Three biological replicates per tissue (collected at 0 and 24 hours) generated 112.13G of raw transcriptome sequencing data. Following quality control, a minimum of 7.89G clean data were obtained for each sample. The Q20 base percentage surpassed 98%, Q30 exceeded 96%, and the GC content was 42.56-44.94% (Supplementary Table 2). The three biological replicates of different tissues and time points were significantly correlated, with correlation coefficients exceeding 0.9 (Supplementary Figure 1A). Furthermore, the PCA results showed that the 12 samples were grouped into four categories, with similar biological replicates clustering based on tissue and treatment (Supplementary Figure 1B). The sequencing data were of high quality and consistent among replicates, making it suitable for further analysis.

Heat stress resulted in the up-regulation of 6129 genes and the down-regulation of 7047 genes in the root when compared to R24 and R0. In contrast, 6621 genes were up-regulated, and 5662 genes were down-regulated compared to S24 and S0 (Supplementary Figures 2A, B; Supplementary Table 3). Venn analysis revealed 2326 up-regulated and 2382 down-regulated genes shared between the root and shoot (Supplementary Figures 2C, D). More genes belonged to specific differentially expressed genes. It is implied that the root and shoot may have a specialized response mechanism to heat stress.

The Gene Ontology (GO) analysis identified 103 significantly enriched GO terms in the root. The enriched Biological Processes (BP) predominantly encompassed stress response-related GO terms, including response to oxidative stress (GO:0006979), DNA integration (GO:0015074), response to stress (GO:0006950), and DNA metabolic process (GO:0006259). The significantly enriched Molecular Functions (MF) includes peroxidase activity (GO:0004601), oxidoreductase activity acting on peroxide as acceptor (GO:0016684), antioxidant activity (GO:0016209), and heme binding (GO:0020037). Finally, the significantly enriched Cellular Components (CC) comprised of the nucleosome (GO:0000786), chromatin (GO:0000785), protein-DNA complex (GO:0032993), and DNA packaging complex (GO:0044815) (Figure 2A; Supplementary Table 4).

Figure 2
www.frontiersin.org

Figure 2 The GO and KEGG enrichment of the DEGs (before and after heat stress) in both the root and shoot. The top ten GO classifications under various stress conditions in each category include biological process (BP), molecular function (MF), and cellular component (CC) in the root tissue (A) and shoot (B). (C) shows the enriched KEGG pathways in the root and shoot. The size of the circle indicates the number of genes, and the color indicates the enrichment significance (p-value).

A total of 110 GO terms were significantly enriched in the shoot tissue. The enriched BP category included protein folding (GO:0006457), cellular response to stimulus (GO:0051716), small molecule metabolic process (GO:0044281), and DNA metabolic process (GO:0006259). The enriched MF category included catalytic activity acting on DNA (GO:0140097), DNA helicase activity (GO:0003678), unfolded protein binding (GO:0051082), and signal transducer activity (GO:0004871). In the CC category, prominent terms comprised mitochondrial matrix (GO:0005759), photosystem II oxygen-evolving complex (GO:0009654), mitochondrial part (GO:0044429), and thylakoid membrane (GO:0042651). The GO terms for DNA metabolic processes were enriched in both tissues (under BP) and mitochondria (under CC). The GO terms enriched in the root predominantly pertain to antioxidant processes, including response to oxidative stress and stress (under BP), as well as peroxidase and antioxidant activities (under MF). In contrast, the GO terms enriched in shoot tissue are primarily associated with signaling responses and transduction, such as cellular response to stimuli and signal transduction (under BP), as well as signal transducer and signaling receptor activities (under MF) (Figure 2B; Supplementary Table 4).

A mapping analysis was conducted using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://genome.jp/kegg/) to examine the regulatory pathways of DEGs induced by heat stress. Heat stress enriched 126 pathways in the root, of which 24 were significantly enriched. The shoot had 126 enriched pathways, including 32 significantly enriched ones (Figure 2C; Supplementary Table 5). Half of the pathways significantly enriched in the root following heat stress were the same as those in the shoot, suggesting a common alteration in the metabolic pathways. These shared pathways include plant hormone signal transduction, glutathione metabolism, starch and sucrose metabolism, MAPK signaling pathway, and glucose metabolism pathway, which are crucial in plant stress response. Metabolic pathways that were notably enriched in the root consist of phenylpropanoid biosynthesis, amino sugar and nucleotide sugar metabolism, and nitrogen metabolism. Conversely, the pathways enriched in the shoot were specific to carbon metabolism, peroxisome, circadian rhythm-plant, and photosynthesis-antenna protein, all associated with the response to heat stress. These findings indicate that distinct adaptive mechanisms mitigate the effects of heat stress on various tissues.

Transcription factors expressed in heat-stressed tobacco seedlings

Transcription factors mediate plant responses to heat stress by regulating target gene expression. Thus, a statistical analysis of the diversity and abundance of TFs within the DEGs revealed 79 classes of TFs in the root, the top five classes being MYB (55), bHLH (49), AP2/ERF-ERF (47), NAC (42), and C2H2 (41). The shoot had 80 classes of transcription factors, and the top five most abundant ones included AP2/ERF-ERF (42), MYB (40), bHLH (38), C2H2 (36), and MYB-related (31). These TFs are involved in plant stress tolerance processes (Mizoi et al., 2012; Li et al., 2015; Guo et al., 2021) (Figures 3A, B; Supplementary Table 6).

Figure 3
www.frontiersin.org

Figure 3 Distribution of transcriptome factors (TFs) in DEGs. (A) Statistical analysis of transcription factors in the root before and after heat treatment. (B) Statistical analysis of TFs in the shoot before and after heat treatment. (C) Expression of NMYB78 in the RNA-seq data. The relative expression of NtMYB78 in the root (D) and shoot (E) following heat treatments for 0, 3, 6, 12, and 24 hours. Significance analysis was conducted using a Student’s t-test with three biological replicates, each with a control treatment (0 hours), *P < 0.05; **P < 0.01; ***P < 0.001. (F) Phenotypic analysis of pTRV : NtMYB78 plants after treatment at 45°C and 24h. pTRV-PDS was the positive control, and the vector was the negative control.

Moreover, MYB and MYB-related transcription factors were the most prominent in the root and shoot, consistent with previous findings that MYB is crucial for enhancing plant stress tolerance (Li et al., 2019; Wang et al., 2021b). Therefore, we analyzed one of the MYB transcription factors, Nitab4.5_0000895g0070, which is homologous to Arabidopsis thaliana AtMYB78, hereafter referred to as NtMYB78. NtMYB78 significantly decreased expression in both tissues after heat stress, with higher levels in the root than shoot (Figure 3C). Heat stress initially decreased NtMYB78, followed by an increase and peaking at 12 h before gradually declining in the root (Figure 3D). In the shoot, NtMYB78 initially increased, peaking at 3 hours, followed by a decrease (Figure 3E). These findings suggest that the shoot responded more rapidly to heat stress than the root. To further verify whether NtMYB78 was involved in heat stress response, we conducted VIGS-mediated knock-down of gene expression. The results revealed that NtMYB78 expression level was significantly down-regulated in the pTRV : NtMYB78 plants (Supplementary Figure 3). The pTRV : NtMYB78 exhibited increased sensitivity following exposure to heat stress, suggesting a role for NtMYB78 in the positive regulation of heat stress response (Figure 3F).

Heat-responsive HSPs in heat-stressed tobacco seedlings

The heat shock proteins (HSPs), significant in the cellular response to heat stress (Tian et al., 2021; Kang et al., 2022; Mondal et al., 2023), were differentially expressed in tobacco roots and shoots under heat stress. In this study, heat stress differentially regulated 106 HSPs in the root, including HSP20 (39), HSP40 (12), HSP70 (26), and HSP90 (29). The shoot had 115 significant differentially expressed HSPs, including HSP20 (37), HSP40 (12), HSP70 (31), and HSP90 (36) (Supplementary Figure 4A; Supplementary Table 6). Using two genes from each HSP family (HSP20, HSP40, HSP70, and HSP90), qRT-PCR analysis revealed consistent expression with the transcriptome data (Supplementary Figure 4B).

Differential metabolomes in heat-stressed tobacco seedlings

A comprehensive LC-MS/MS identified 647 metabolites in the root samples. The metabolites are predominantly grouped into categories of amino acids and their derivatives (21.51%), carbohydrates and their derivatives (12.02%), and lipids (11.57%) (Supplementary Figure 5A; Supplementary Table 7). The shoot contained 932 metabolites, including amino acids and their derivatives (19.31%), carbohydrates and their derivatives (10.41%), and organic acids, lipids, and their derivatives (collectively accounting for 10.19%) (Supplementary Figure 5B; Supplementary Table 7). The PCA of the root data revealed that PC1 accounted for 50.87% of the total variation, while PC2 explained 17.7% of the variation (Supplementary Figure 5C). The PCA of the shoot showed that PC1 and PC2 accounted for 51.91% and 15.42% of the total variation, respectively (Supplementary Figure 5D). Furthermore, the biological replicates within each group clustered in the PCA results, providing additional evidence of the group distinctions and indicating the reliability and suitability of the metabolite data for subsequent analysis.

The root contained 116 DAMs (107 up-regulated and 9 down-regulated). The top 20 up-regulated compounds identified in the study included alpha-trehalose, isomaltulose, melezitose, coniferin, turanose, beta-D-Lactose, raffinose, sucrose, 1,4-dihydro-1-Methyl-4-oxo-3-pyridinecarboxamid, selgin O-hexosyl-O-hexoside, petunidin 3-O-glucoside, anthranilate O-hexosyl-O-hexoside, 2,6-Dihydroxypurine, iP7G, lactose, selgin 5-O-hexoside, maltose, methylQuercetin O-hexoside, 3-O-p-Coumaroylquinic acid, and D-(+)-cellobiose. Notably, sugar compounds accounted for most (11) up-regulated compounds. The nine most significantly down-regulated DAMs included LysoPC 10:0, isomangiferolic acid, L-threo-3-methylaspartate, 2’-deoxyguanosine, ligustrazine, 4-hydroxymandelonitrile, N-acetyl-L-methionine, eucommiol, and aminophylline. Terpenoids and amino acids were the most abundant classes, represented by two compounds (Figures 4A, C; Supplementary Table 7). The S24 vs S0 comparison group had 256 DAMs, including 251 up-regulated and 5 down-regulated compounds (Figure 4B). Further, the top twenty DAMs were raffinose, sn-glycero-3-phosphocholine, melezitose, 1,4-dihydro-1-methyl-4-oxo-3-pyridinecarboxamide, and 2,6-Dihydroxypurine, aspartic acid di-O-glucoside, selgin O-hexosyl-O-hexoside, alpha-trehalose, syringetin 5-O-hexoside, galactinol, beta-D-Lactose, sucrose, coniferin, melibiose, isomaltulose, lactose, isomaltose, maltose, 1-caffeoylquinic acid, and lysoPC 20:4. Notably, sugar compounds were the most abundant (12) among the top twenty DAMs. The top five down-regulated DAMs included 2-aminoadipic acid, L-2-aminoadipic acid, D-lactic acid, 2-methyladenosine, and lactic acid (Figure 4D). There were 51 co-up-regulated DAMs, carbohydrates and their derivatives being abundant (17), including lactose, raffinose, trehalose, sucrose, and sorbose. Amino acids and their derivatives were the second most abundant (10), including asparagine, histidine, nepsilon-acetyl-L-lysine, and proline (Supplementary Figures 6A, C). There were no shared and down-regulated DAMs (Supplementary Figure 6B). The shoot had significantly more DAMs than the root, suggesting substantial alterations in shoot metabolite following heat treatment. This disparity may be attributed to the heightened susceptibility of the shoot to heat stress relative to the root.

Figure 4
www.frontiersin.org

Figure 4 Differential metabolite analysis. (A, B) show a volcano map of DAMs in R24_VS_R0 and S24_VS_S0. Red dots indicate the up-regulation of DAMs, green dots indicate the down-regulation of DAMs, gray dots indicate no DAMs, and the circle size indicates the VIP value. (C, D) show the top 20 DAMs that increased and decreased in R24_VS_R0 and S24_VS_S0.

Furthermore, the root contained 33 significantly enriched KEGG pathways, including starch and sucrose metabolism, glycerolipid metabolism, and caffeine metabolism (Figure 5; Supplementary Table 8). The shoot DAMs enriched 41 pathways, the most significantly enriched being galactose metabolism, gluconeogenesis, starch and sucrose metabolism, and zeatin biosynthesis (Figure 5; Supplementary Table 8). The two tissues were co-enriched with various metabolic processes such as starch and sucrose metabolism, glycerolipid metabolism, galactose metabolism, ABC transporter, and purine metabolism. The pathways enriched in the root included caffeine metabolism, nicotinate and nicotinamide metabolism, lysine degradation, and the pentose phosphate pathway. In contrast, the shoot was primarily enriched with benzoxazinoid biosynthesis, phenylalanine, tyrosine, and tryptophan biosynthesis, pentose and glucuronate interconversions, and photosynthesis (Figure 5; Supplementary Table 8). In conclusion, shoots may demonstrate heightened complexity in their response mechanisms to heat stress.

Figure 5
www.frontiersin.org

Figure 5 The differential metabolite KEGG enrichment analysis was conducted, with bubble plots illustrating the top 22 significantly enriched pathways in both the root and shoot tissues. The size of each bubble corresponds to the number of metabolites, while the color of the bubble represents the significance level as indicated by the p-value.

Integrated transcriptomics and metabolomics identified key heat-regulated metabolic pathways

A combined analysis of the transcriptome and metabolites was performed on both the root and shoot tissues to further explore critical pathways of heat stress response. Heat stress regulated (DEGs and DAMs) and enriched several functional pathways in both tissues, including sucrose starch metabolism, glycerolipid metabolism, and galactose metabolism. Pathways for amino sugar and nucleotide sugar metabolism, glyoxylate and dicarboxylate metabolism, and alpha-linolenic acid metabolism were enriched exclusively in the root. In contrast, carbon metabolism, plant hormone signal transduction, pyruvate metabolism, and purine metabolism were specifically enriched in the shoot (Figure 6).

Figure 6
www.frontiersin.org

Figure 6 The (A, B) column diagrams show the KEGG co-enrichment pathway of DEGs (red column) and DAMs (blue column) in the root and shoot.

Gene-metabolite correlation networks revealed a strong correlation (cor > 0.9 and p-value < 0.01) between the differential genes and metabolites in critical metabolic pathways. The root had major metabolic pathways, amino sugar and nucleotide sugar metabolism, enriched by 17 and 3 significantly correlated DEGs and DAMs, respectively. The corresponding enriched metabolic pathways included UDP-galactose, which is involved in the galactose metabolism pathway for response to heat stress (Figure 7A; Supplementary Figure 7A). Heat stress significantly enriched the purine metabolic pathway in the shoot, with 28 DEGs that were strongly correlated (predominantly displaying a positive correlation) to 11 DAMs. Heat treatment also generated DEGs and DAMs that enriched the purine metabolic pathway in the shoot (Figure 7B; Supplementary Figure 7B). Furthermore, heat stress up-regulated adenine phosphoribosyltransferase (APT), adenylate kinase (AK), and nucleoside-diphosphate kinase (NDPK) (Supplementary Table 3). Heat stress notably elevated the concentrations of guanine, xanthine, hypoxanthine, and adenine, pivotal intermediates in purine metabolism (Supplementary Table 7).

Figure 7
www.frontiersin.org

Figure 7 Major enriched metabolic pathways in the root and shoot after heat treatment. (A) DAMs and DEGs that significantly enriched the amino sugar and nucleotide sugar metabolism pathways in the root. (B) DAMs and DEGs that significantly enriched the purine metabolism pathway in the shoot. (C) DAMs and DEGs that co-enriched the galactose metabolism and starch and sucrose metabolism pathways in the root and shoot after heat treatment. Red represents DAMs, and blue represents DEGs.

Heat stress significantly altered the abundance of sugars in both tissues. Specifically, the enriched galactose metabolic pathway in the root contained 12 DEGs, significantly associated with 7 DAMs. In the shoot, 24 DEGs were correlated with 15 DAMs (Supplementary Figure 8A). Similarly, the enriched sucrose starch metabolism pathway in the root contained 30 DEGs, which was significantly correlated with 7 DAMs, and the shoot had 38 DEGs, which was significantly associated with 7 DAMs (Supplementary Figure 8B). Sugar metabolism-related DEGs and DAMs were more abundant in the shoot than in the root. Subsequently, we integrated our RNA-seq and metabolomics datasets to contrust a gene-metabolite network (Figure 7C). The qRT-PCR results indicated a significant increase in the expression of raffinose synthase (RS), sucrose synthase (SUS), trehalose 6-phosphate synthase (TPS), and trehalose 6-phosphate phosphatase (TPP), consistent with the transcriptome data (Supplementary Figure 8C; Supplementary Table 3). Metabolomic analysis revealed significant increases in raffinose, sucrose, and trehalose contents in both root and shoot, with 52 - 419 fold changes (Supplementary Table 7). These results suggested that these genes may directly or indirectly regulate the accumulation of sugar metabolites in response to high-temperature stress.

Discussion

The shoot of tobacco seedlings are more sensitive to heat stress than the root

Extensive research has been conducted on the effects of heat stress on plant shoots; however, there remains a notable gap in the literature regarding the response of plant roots to heat stress (Wang et al., 2019; Liu et al., 2020; Xie et al., 2023). Thus, this study simultaneously exposed tobacco seedling roots and shoots to elevated temperatures to explore the comparative responses of roots and shoots to heat stress. Subsequent transcriptome analyses on the both tissues identified 13176 and 12283 differentially expressed genes (DEGs) in the root and shoot, respectively. The enriched GO terms in the root were associated with stress response and antioxidant pathways, whereas those in the shoot were predominantly related to plant signaling pathways. These findings are consistent with the GO-enriched terms observed in roots and leaves during a previous study on drought stress in switchgrass (Tiedge et al., 2022).

The root and shoot exhibited significant enrichment of KEGG pathways, with 24 pathways identified in the root and 32 in the shoot. More than half of the metabolic pathways were common to both tissues, including plant hormone signal transduction, MAPK signaling pathway, and starch and sucrose metabolism. The significant metabolic pathways enriched in response to heat stress align with previous reports on Arabidopsis and pepper (Wang et al., 2019; Olas et al., 2021). Heat stress mainly enriched phenylpropane metabolism, amino acid sugar and nucleic acid sugar metabolism in the roots. Phenylpropane metabolism is a key secondary pathway in plants, producing important metabolites like lignin, anthocyanins, and organic acids that regulate plant growth and stress response (Dong and Lin, 2021). Heat stress significantly enriched pathways related to carbon metabolism, purine metabolism, and peroxisomes in the shoot tissue. Furthermore, heat stress significantly elevated antioxidant enzyme activities in the shoot (Figure 1B), consistent with prior research (Silva et al., 2018).

The root and shoot had 647 and 932 metabolites, including 115 DAMs and 255 DAMs, respectively, suggesting that heat stress may affect the root less. There were 33 pathways identified for DAMs enrichment in the root and 41 pathways for DAMs enrichment in the shoot, indicating a greater complexity of metabolic pathways in the shoot response to heat stress.This result aligns with previous studies showing wheat shoots experience more metabolite changes during drought tolerance (Kang et al., 2019). The shoot exhibited a greater sensitivity to heat stress, as evidenced by more metabolite alterations and enriched KEGG pathways.

Purine metabolism plays a key role in responding to heat stress

The integrated transcriptomic and metabolomic data indicated that the purine metabolism pathway was predominantly enriched in the shoot, and gene-metabolite correlations identified 28 and 11 significantly linked DEGs and DAMs. Therefore, purine metabolism may play a crucial role in the response to heat stress. Purines are necessary metabolites in all organisms since they regulate various cellular functions, such as cell signaling, redox metabolism, and energy metabolism (Kokina et al., 2019). Elevated temperatures notably increased purine metabolic compounds in the high-altitude fish Triplophysa siluroides (Chen et al., 2023). In quinoa, heat stress enriched the purine metabolic pathway and up-regulated substances related to purine metabolism in both heat-sensitive and heat-resistant varieties (Xie et al., 2023). The purine metabolic compounds adenine, xanthine, hypoxanthine, and guanine increased notably within the shoot (Figure 7B). A recent study found that cAMP, the first identified second messenger, influences heat stress response by regulating various cellular processes such as protein processing, ion balance, and the ubiquitin-proteasome system (Liang et al., 2022). Interestingly, the metabolome results demonstrated a significant elevation of ATP, which is a precursor substance of cAMP, and a 1.2-fold elevation of cAMP (Supplementary Table 7). It has been reported that the accumulation of allantoin, an intermediate product of purine metabolism, promotes ABA synthesis and participates in the response to adversity stress (Watanabe et al., 2014; Kaur et al., 2023). However, the content of allantoin in the shoot did not change significantly in this study, probably because of the early heat stress response. Consequently, this study analyzed genes associated with ABA synthesis and signaling, revealing a substantial up-regulation of NCED1 (Nitab4.5_0001924g0060) (Supplementary Figures 9A, C). This pivotal gene implicated in ABA synthesis increased by 87-fold in the shoot expression, and ABA levels elevated by 8-fold (Supplementary Tables 3, 7). Furthermore, heat stress increased the expression of the ABA-responsive element-binding protein (ABF), indicating their role in modulating the expression of ABA-responsive genes under heat stress (Supplementary Figure 9B; Supplementary Table 3).

MYB transcription factors play important roles in ABA signaling, and they can be regulated by ABA to influence plant responses to abiotic stresses (Ma et al., 2023; Du et al., 2024). These MYB transcription factors can also interact with MYB cis-elements in the promoters of target genes in response to abiotic stress (Park et al., 2011; Zheng et al., 2012; Wang et al., 2021b). Heat stress significantly altered 71 MYB and MYB-related transcription factors in the shoot (20 up-regulated and 51 down-regulated (Supplementary Table 6). AtMYB78 in Arabidopsis thaliana is crucial for resistance to drought, salt, and oxidative stresses (Mengiste et al., 2003). The homologs of AtMYB78 are equally significant in mediating responses to both biotic and abiotic stresses in wheat, tomato, and oilseed rape (Abuqamar et al., 2009; Liu et al., 2011; Chen et al., 2016). We therefore analyzed the expression profile of NtMYB78, exhibited a pattern of initial increase followed by decrease in both root and shoot tissues under varying durations of heat stress treatment. Notably, the shoot displayed a more rapid response, reaching peak expression levels after 3 hours of heat stress treatment (Figures 3D, E). This observation aligns with findings from a previous study on heat-stressed wheat, where MYB-associated transcription factors rapidly accumulated within one hour of heat stress (Zhao et al., 2017). Conversely, certain MYB transcription factors, such as AtMYBS1, have been shown to exert a negative regulatory effect on heat stress by inhibiting the expression of the MAX1 gene (Li et al., 2023). This was why more MYB transcription factors were down-regulated after 24 hours of heat stress.

MYB transcription factors can interact with HSFs to regulate HSP protein expression (Jacob et al., 2017; Priya et al., 2019). These HSPs are molecular chaperones for the assembly, stabilization, and maturation of proteins and protein complexes. They are also important in plant development and responses to abiotic and biotic stresses (ul Haq et al., 2019). The up-regulation of BcHSP70 and NtHSP70-8b notably enhances heat tolerance in tobacco (Wang et al., 2016; Zhang et al., 2023). In this study, heat stress differentially expressed 106 (80 up-regulated, 26 down-regulated) and 115 (83 up-regulated, 32 down-regulated) HSPs in the root and shoot, respectively (Supplementary Figure 2). Thus, increasing the expression of HSPs could improve heat tolerance tobacco.

Sugar metabolism pathways have a significant impact on heat tolerance

Sugars, including raffinose, trehalose, and sucrose, are osmoregulatory compounds that aid in plant adaptation to molecular, cellular, and physiological alterations (Sami et al., 2016; Raza et al., 2023). Specifically, UDP-galactose and myo-inositol use HSFA1 to modulate GALACTITOL SYNTHASE (GOLS1) expression, synthesizing galactitol. Galactitol produces raffinose using RAFFINOSE SYNTHASE (RS) in a temperature-dependent condition (Panikulangara et al., 2004; Cortijo et al., 2017). Therefore, heat-stressed plants improve heat tolerance by increasing the content of raffinose (Serrano et al., 2019; Yan et al., 2022). This study identified notable alterations in GOLS and RS expression levels in the root and shoot of heat-stressed tobacco. Specifically, heat stress significantly increased RS (Nitab4.5_0002682g0060) expression in both tissues by 6- and 9-fold, respectively. (Figure 7C; Supplementary Table 3). Metabolomic analyses further revealed a substantial 71- and 356-fold rise in the raffinose content in the root and shoot, respectively (Figure 7C; Supplementary Table 7). A recent study showed that trehalose-6-phosphate synthase 1 can regulate heat tolerance in A. thaliana by inducing raffinose accumulation (Reichelt et al., 2023). Trehalose biosynthesis in plants involves converting UDP-glucose into trehalose-6P using TPS and then converting TPS to trehalose using TPP (Avonce et al., 2006). Trehalose is also an osmoprotectant that prevents cell dehydration, stabilizes macromolecules, scavenges ROS, and improves cellular antioxidant properties (Raza et al., 2023). Applying a small quantity of trehalose can mitigate the detrimental effects of extreme temperatures on plants, enhance growth and development, and positively enhance plant resistance (Luo et al., 2018, 2022a). The heightened expression of sucrose and trehalose-related genes observed in this study notably augmented sucrose and trehalose levels (Figure 7C; Supplementary Figure 8). Prior research has demonstrated that trehalose up-regulates stress response genes, prompting the accumulation of osmolytes like proline, betaine, and soluble sugars, thereby influencing stress tolerance (Hassan et al., 2023). Moreover, heat stress significantly elevated proline levels in the root and shoot (Figure 1B; Supplementary Table 7). Soluble sugars can also reduce membrane osmotic potential and maintain cell expansion, possibly explaining why higher levels of sugars like glucose, sucrose, and trehalose can improve plant tolerance to abiotic stress.

Collectively, heat stress-related transcription factors (MYB, NAC, bHLH, and HSF) can regulate response to heat stress in the root and shoot through the plant hormone signal transduction pathway, which regulates HSP expression. Heat shock proteins can interact with osmotic regulators (trehalose and proline) and antioxidant enzymes (CAT, POD, and SOD) to regulate heat stress response. The root can engage in the pathways for sugar metabolism via the amino sugar and nucleotide sugar metabolism pathway to enhance raffinose, sucrose, and trehalose contents. The shoot can partake in synthesizing and signaling the phytohormone ABA through purine metabolism. This process can regulate sugar accumulation and ultimately enhance the tolerance of tobacco seedlings to high heat stress (Figure 8).

Figure 8
www.frontiersin.org

Figure 8 Model of heat stress response pathways in the root and shoot of tobacco seedlings.

Conclusion

This study integrated physiological and biochemical analyses with transcriptomic and metabolomic approaches to examine the molecular pathways involved in the response of tobacco seedling roots and shoots to heat stress. The primary findings of the study are outlined below. Heat stress significantly elevated the activities of antioxidant enzymes, including CAT, POD, and SOD in the shoot. Additionally, heat stress significantly increased the levels of osmotic mediators, such as soluble sugars and proline. Amino and nucleotide sugar metabolism were the primary metabolic pathways enriched in the root of heat-stressed tobacco, while purine metabolism was the predominant pathway in the shoot. These pathways can collectively enhance sugar metabolism, increasing the levels of sucrose, trehalose, raffinose, and other compounds that improve tolerance to heat stress in tobacco. While we have outlined the impact of the shoot on ABA levels via purine metabolism, leading to the regulation of HSP proteins and osmoregulatory substances through MYB and HSF transcription factors, as well as identified the positive role of the MYB transcription factor NtMYB78 in heat stress regulation, the specific regulatory mechanisms remain unclear. Further validation of these molecular functions is required.

Data availability statement

The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/, the project number is PRJNA1093408.

Author contributions

HC: Writing – review & editing, Writing – original draft. SQ: Writing – review & editing, Investigation, Methodology. YC: Resources, Validation, Writing – review & editing. JL: Writing – review & editing. TX: Resources, Writing – review & editing. PZ: Resources, Writing – review & editing. XS: Writing – review & editing. SX: Resources, Writing – review & editing. ZM: Resources, Writing – review & editing. ZH: Supervision, Writing – review & editing. XP: Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was financially supported by the NSFC grant (32000214), the talent training program of Guangdong Academy of Agricultural Sciences (R2022PY-QY003), the foundation of director of Crops Research Institute, Guangdong Academy of Agricultural Sciences (202303), Guangdong S&T programme (No. 2023B1212060038), and the science and technology project of Meizhou Tobacco Corporation (202301).

Conflict of interest

YC, PZ and SX were employed by China National Tobacco Corporation, Guangdong Company.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1425944/full#supplementary-material

Supplementary Figure 1 | Principle component analysis (PCA) of the transcriptome data. (A) The correlation heatmap of the 12 samples. R0 and S0 denote the root and shoot samples subjected to control conditions (25°C, 24 h), while R24 and S24 denote the samples subjected to heat stress (45°C, 24 h). The numbers 1, 2, and 3 correspond to the three biological replicates. (B) Principle components of differentially expressed genes (DEGs) in the root and shoot before and after heat treatment.

Supplementary Figure 2 | Identification of DEGs. Root (A) and shoot (B) volcano plots of DEGs after heat treatments (before and after 24h). Red dots denote up-regulated genes, green dots denote down-regulated genes, and blue dots denote genes that were not differentially expressed. (C) The Venn diagram illustrates the up-regulated genes between the root and shoot tissues (D). The Venn diagram illustrates the down-regulated genes in roots versus shoots.

Supplementary Figure 3 | Relative expression levels of VIGS-mediated knock-down of gene NtMYB78.

Supplementary Figure 4 | Expression of heat shock proteins (HSPs). Heat map of HSPs expression in the root (A) and shoot (B) before and after heat treatment with three biological replicates per treatment. (C) qRT-PCR validated the expression of HSP20 (2), HSP40 (2), HSP70 (2) and HSP90 (2). The expression level was normalized to that of NtACTIN. Data are given as means ± SD of three biological replicates. A significance analysis was performed using the Student’s t-test. *P < 0.05, ***P < 0.001.

Supplementary Figure 5 | Distribution of identified metabolite species and PCA analysis. (A, B) show the classification and percentage of metabolites identified in the roots and shoots. PCA plots of total detected metabolites of R24 and R0 (C), S24 and S0 (D).

Supplementary Figure 6 | Differential metabolite veen analysis. (A, B) indicates the veen graph that up-regulated DAMs and down-regulated DAMs in the root and shoot, respectively. (C) Heatmap of co-up-regulated DAMs in the root and shoot.

Supplementary Figure 7 | Network diagram of key metabolic pathway genes and metabolite correlation. (A) The network diagram illustrates the correlations between genes and metabolites within the root-specific amino sugar and nucleotide sugar metabolism pathway that has been significantly enriched. (B) The network diagram illustrates the gene-metabolite correlations within the shoot-specific significantly enriched purine metabolic pathway. Green circles represent metabolites, while pink circles represent genes. Red lines denote positive correlations, while blue lines denote negative correlations.

Supplementary Figure 8 | A network diagram of key metabolic genes and metabolite correlation. (A) A network of gene-metabolite correlations in the root and shoot featuring the significantly enriched galactose metabolism pathway. (B) A network of gene-metabolite correlations in the root and shoot featuring the significantly enriched starch and sugar metabolism pathway. Green circles indicate metabolites. Pink circles indicate genes. The red line indicates a positive correlation, and the blue line indicates a negative correlation. (C) Validation of RS, SUS, TPS, and TPP expression. Data was given as means ± SD of three biological replicates. A significance analysis was performed using the Student’s t-test. **P < 0.01, ***P < 0.001.

Supplementary Figure 9 | Heat map displaying the expression of ABA synthesis genes (A) and signaling genes (B), with log2FC values representing alterations in gene expression. (C) The relative expression of NCED1 in tobacco root and shoot. (D) The relative ABA content in the shoot.

References

Abdel-Basset, R., Ozuka, S., Demiral, T., Furuichi, T., Sawatani, I., Baskin, T. I., et al. (2010). Aluminium reduces sugar uptake in tobacco cell cultures: a potential cause of inhibited elongation but not of toxicity. J. Exp. Bot. 61, 1597–1610. doi: 10.1093/jxb/erq027

PubMed Abstract | CrossRef Full Text | Google Scholar

Abuqamar, S., Luo, H., Laluk, K., Mickelbart, M. V., Mengiste, T. (2009). Crosstalk between biotic and abiotic stress responses in tomato is mediated by the AIM1 transcription factor. Plant J. 58, 347–360. doi: 10.1111/j.1365-313X.2008.03783.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Andrási, N., Pettkó-Szandtner, A., Szabados, L. (2021). Diversity of plant heat shock factors: regulation, interactions, and functions. J. Exp. Bot. 72, 1558–1575. doi: 10.1093/jxb/eraa576

PubMed Abstract | CrossRef Full Text | Google Scholar

Avonce, N., Mendoza-Vargas, A., Morett, E., Iturriaga, G. (2006). Insights on the evolution of trehalose biosynthesis. BMC Evol. Biol. 6, 109. doi: 10.1186/1471-2148-6-109

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, B., Niu, F., Liu, W. Z., Yang, B., Zhang, J., Ma, J., et al. (2016). Identification, cloning and characterization of R2R3-MYB gene family in canola (Brassica napus L.) identify a novel member modulating ROS accumulation and hypersensitive-like cell death. DNA Res. 23, 101–114. doi: 10.1093/dnares/dsv040

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Wu, X., Li, P., Liu, Y., Song, M., Li, F., et al. (2023). Integrated metabolomic and transcriptomic responses to heat stress in a high-altitude fish, Triplophysa siluroides. Fish Shellfish Immun. 142, 109118. doi: 10.1016/j.fsi.2023.109118

CrossRef Full Text | Google Scholar

Chng, W. B. A., Sleiman, M. S. B., Schüpfer, F., Lemaitre, B. (2014). Transforming growth factor β/activin signaling functions as a sugar-sensing feedback loop to regulate digestive enzyme expression. Cell Rep. 9, 336–348. doi: 10.1016/j.celrep.2014.08.064

PubMed Abstract | CrossRef Full Text | Google Scholar

Cortijo, S., Charoensawan, V., Brestovitsky, A., Buning, R., Ravarani, C., Rhodes, D., et al. (2017). Transcriptional regulation of the ambient temperature response by H2A.Z nucleosomes and HSF1 transcription factors in Arabidopsis. Mol. Plant 10, 1258–1273. doi: 10.1016/j.molp.2017.08.014

PubMed Abstract | CrossRef Full Text | Google Scholar

de Pinto, M. C., Locato, V., Paradiso, A., De, Gara., L. (2015). Role of redox homeostasis in thermo-tolerance under a climate change scenario. Ann. Bot. 116, 487–496. doi: 10.1093/aob/mcv071

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, N. Q., Lin, H. X. (2021). Contribution of phenylpropanoid metabolism to plant development and plant-environment interactions. J. Integr. Plant Biol. 63, 180–209. doi: 10.1111/jipb.13054

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, S. X., Wang, L. L., Yu, W. P., Xu, S. X., Chen, L., Huang, W. (2024). Appropriate induction of TOC1 ensures optimal MYB44 expression in ABA signaling and stress response in Arabidopsis. Plant Cell Environ. 47, 3046–3062. doi: 10.1111/pce.14922

PubMed Abstract | CrossRef Full Text | Google Scholar

Dunn, W. B., Broadhurst, D., Begley, P., Zelena, E., Francis-McIntyre, S., Anderson, N., et al. (2011). Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry. Nat. Protoc. 6, 1060–1083. doi: 10.1038/nprot.2011.335

PubMed Abstract | CrossRef Full Text | Google Scholar

Durand, M., Porcheron, B., Hennion, N., Maurousset, L., Lemoine, R., Pourtau, N. (2016). Water deficit enhances c export to the roots in Arabidopsis thaliana plants with contribution of sucrose transporters in both shoot and roots. Plant Physiol. 170, 1460–1479. doi: 10.1104/pp.15.01926

PubMed Abstract | CrossRef Full Text | Google Scholar

Fortunato, S., Lasorella, C., Dipierro, N., Vita, F., de Pinto, M. C. (2023). Redox signaling in plant heat stress response. Antioxidants. 12, 605. doi: 10.3390/antiox12030605

CrossRef Full Text | Google Scholar

Guihur, A., Rebeaud, M. E., Goloubinoff, P. (2022). How do plants feel the heat and survive? Trends Biochem. Sci. 47, 824–838. doi: 10.1016/j.tibs.2022.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, J., Sun, B., He, H., Zhang, Y., Tian, H., Wang, B. (2021). Current understanding of bHLH transcription factors in plant abiotic stress tolerance. Int. J. Mol. Sci. 22, 4921. doi: 10.3390/ijms22094921

PubMed Abstract | CrossRef Full Text | Google Scholar

Haider, S., Iqbal, J., Naseer, S., Yaseen, T., Shaukat, M., Bibi, H., et al. (2021). Molecular mechanisms of plant tolerance to heat stress: current landscape and future perspectives. Plant Cell Rep. 40, 2247–2271. doi: 10.1007/s00299-021-02696-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Halford, N. G., Hey, S., Jhurreea, D., Laurie, S., McKibbin, R. S., Paul, M., et al. (2003). Metabolic signalling and carbon partitioning: role of Snf1-related (SnRK1) protein kinase. J. Exp. Bot. 54, 467–475. doi: 10.1093/jxb/erg038

PubMed Abstract | CrossRef Full Text | Google Scholar

Hassan, M. U., Nawaz, M., Shah, A. N., Raza, A., Barbanti, L., Skalicky, M., et al. (2023). Trehalose: a key player in plant growth regulation and tolerance to abiotic stresses. J. Plant Growth Regul. 42, 4935–4957. doi: 10.1007/s00344-022-10851-7

CrossRef Full Text | Google Scholar

Hu, Z., Fan, J., Xie, Y., Amombo, E., Liu, A., Gitau, M. M., et al. (2016). Comparative photosynthetic and metabolic analyses reveal mechanism of improved cold stress tolerance in Bermudagrass by exogenous melatonin. Plant Physiol. Biochem. 100, 94–104. doi: 10.1016/j.plaphy.2016.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, Z., He, Z., Li, Y., Wang, Q., Yi, P., Yang, J., et al. (2023). Transcriptomic and metabolic regulatory network characterization of drought responses in tobacco. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.1067076

CrossRef Full Text | Google Scholar

Jacob, P., Brisou, G., Dalmais, M., Thévenin, J., van der Wal, F., Latrasse, D., et al. (2021). The seed development factors TT2 and MYB5 regulate heat stress response in Arabidopsis. Genes. 12, 746. doi: 10.3390/genes12050746

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacob, P., Hirt, H., Bendahmane, A. (2017). The heat-shock protein/chaperone network and multiple stress resistance. Plant Biotechnol. J. 15, 405–414. doi: 10.1111/pbi.12659

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, J., Zhang, H., Zhang, J., Liu, P., Chen, X., Li, Z., et al. (2017). Integrated transcriptomics and metabolomics analysis to characterize cold stress responses in Nicotiana tabacum. BMC Genom. 18, 496. doi: 10.1186/s12864-017-3871-7

CrossRef Full Text | Google Scholar

Kan, Y., Mu, X. R., Gao, J., Lin, H. X., Lin, Y. (2023). The molecular basis of heat stress responses in plants. Mol. Plant 16, 1612–1634. doi: 10.1016/j.molp.2023.09.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, Z., Babar, M. A., Khan, N., Guo, J., Khan, J., Islam, S., et al. (2019). Comparative metabolomic profiling in the roots and leaves in contrasting genotypes reveals complex mechanisms involved in post-anthesis drought tolerance in wheat. PloS One 14, e0213502. doi: 10.1371/journal.pone.0213502

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, Y., Lee, K., Hoshikawa, K., Kang, M., Jang, S. (2022). Molecular bases of heat stress responses in vegetable crops with focusing on heat shock factors and heat shock proteins. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.837152

CrossRef Full Text | Google Scholar

Kaur, R., Chandra, J., Varghese, B., Keshavkant, S. (2023). Allantoin: a potential compound for the mitigation of adverse effects of abiotic stresses in plants. Plants. 12, 3059. doi: 10.3390/plants12173059

CrossRef Full Text | Google Scholar

Kim, D., Paggi, J. M., Park, C., Bennett, C., Salzberg, S. L. (2019). Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915. doi: 10.1038/s41587-019-0201-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Kokina, A., Ozolina, Z., Liepins, J. (2019). Purine auxotrophy: Possible applications beyond genetic marker. Yeast. 36, 649–656. doi: 10.1002/yea.3434

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Han, G., Sun, C., Sui, N. (2019). Research advances of MYB transcription factors in plant stress resistance and breeding. Plant Signal Behav. 14, 1613131. doi: 10.1080/15592324.2019.1613131

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Lu, J., Zhu, X., Dong, Y., Liu, Y., Chu, S., et al. (2023). AtMYBS1 negatively regulates heat tolerance by directly repressing the expression of MAX1 required for strigolactone biosynthesis in Arabidopsis. Plant Commun. 4, 100675. doi: 10.1016/j.xplc.2023.100675

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Ng, C. K. Y., Fan, L. M. (2015). MYB transcription factors, active players in abiotic stress signaling. Environ. Exp. Bot. 114, 80–91. doi: 10.1016/j.envexpbot.2014.06.014

CrossRef Full Text | Google Scholar

Liang, S., Sun, J., Luo, Y., Lv, S., Chen, J., Liu, Y., et al. (2022). cAMP is a promising regulatory molecule for plant adaptation to heat stress. Life. 12, 885. doi: 10.3390/life12060885

CrossRef Full Text | Google Scholar

Liao, Y., Smyth, G. K., Shi, W. (2014). FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 30, 923–930. doi: 10.1093/bioinformatics/btt656

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, M., Ju, Y., Min, Z., Fang, Y., Meng, J. (2020). Transcriptome analysis of grape leaves reveals insights into response to heat acclimation. Sci. Hortic. 272, 109554. doi: 10.1016/j.scienta.2020.109554

CrossRef Full Text | Google Scholar

Liu, H., Zhou, X., Dong, N., Liu, X., Zhang, H., Zhang, Z. (2011). Expression of a wheat MYB gene in transgenic tobacco enhances resistance to Ralstonia solanacearum, and to drought and salt stresses. Funct. Integr. Genomics 11, 431–443. doi: 10.1007/s10142-011-0228-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, J., Xing, X. J., Zhu, L. Q., Wang, Y., Yin, H., Yuan, J. J. (2011). Effects of exogenous glycine betaine and CaCl2 on physiological responses of tobacco plants under stresses of heat and drought. J. Plant Nutr. 17, 1437–1443. doi: 10.11674/zwyf.2011.1054

CrossRef Full Text | Google Scholar

Luo, Y., Wang, W., Fan, Y. Z., Gao, Y. M., Wang, D. (2018). Exogenously-supplied trehalose provides better protection for d1 protein in winter wheat under heat stress. Russ J. Plant Physiol. 65, 115–122. doi: 10.1134/S1021443718010168

CrossRef Full Text | Google Scholar

Luo, Y., Wang, Y., Xie, Y., Gao, Y., Li, W., Lang, S. (2022a). Transcriptomic and metabolomic analyses of the effects of exogenous trehalose on heat tolerance in wheat. Inter. J. Mol. Sci. 23, 5194. doi: 10.3390/ijms23095194

CrossRef Full Text | Google Scholar

Luo, Z., Zhou, Z., Li, Y., Tao, S., Hu, Z. R., Yang, J. S., et al. (2022b). Transcriptome-based gene regulatory network analyses of differential cold tolerance of two tobacco cultivars. BMC Plant Biol. 22, 369. doi: 10.1186/s12870-022-03767-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, R., Liu, B., Geng, X., Ding, X., Yan, N., Sun, X., et al. (2023). Biological function and stress response mechanism of MYB transcription factor family genes. J. Plant Growth Regul. 42, 83–95. doi: 10.1007/s00344-021-10557-2

CrossRef Full Text | Google Scholar

Mathur, S., Agrawal, D., Jajoo, A. (2014). Photosynthesis: response to high temperature stress. J. Photochem. Photobiol. B. 137, 116–126. doi: 10.1016/j.jphotobiol.2014.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Mengiste, T., Chen, X., Salmeron, J., Dietrich, R. (2003). The BOTRYTIS SUSCEPTIBLE1 gene encodes an R2R3MYB transcription factor protein that is required for biotic and abiotic stress responses in Arabidopsis. Plant Cell. 15, 2551–2565. doi: 10.1105/tpc.014167

PubMed Abstract | CrossRef Full Text | Google Scholar

Mizoi, J., Shinozaki, K., Yamaguchi-Shinozaki, K. (2012). AP2/ERF family transcription factors in plant abiotic stress responses. Biochim. Biophys. Acta 1819, 86–96. doi: 10.1016/j.bbagrm.2011.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Mondal, S., Karmakar, S., Panda, D., Pramanik, K., Bose, B., Singhal, R. ,. K. (2023). Crucial plant processes under heat stress and tolerance through heat shock proteins. Plant Stress. 10, 100227. doi: 10.1016/j.stress.2023.100227

CrossRef Full Text | Google Scholar

Niu, W. T., Han, X. W., Wei, S. S., Shang, Z. L., Wang, J., Yang, D. W., et al. (2020). Arabidopsis cyclic nucleotide-gated channel 6 is negatively modulated by multiple calmodulin isoforms during heat shock. J. Exp. Bot. 71, 90–104. doi: 10.1093/jxb/erz445

PubMed Abstract | CrossRef Full Text | Google Scholar

Nunes, C., Primavesi, L. F., Patel, M. K., Martinez-Barajas, E., Powers, S. J., Sagar, R., et al. (2013). Inhibition of SnRK1 by metabolites: tissue-dependent effects and cooperative inhibition by glucose 1-phosphate in combination with trehalose 6-phosphate. Plant Physiol. Biochem. 63, 89–98. doi: 10.1016/j.plaphy.2012.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Ohama, N., Sato, H., Shinozaki, K., Yamaguchi-Shinozaki, K. (2017). Transcriptional regulatory network of plant heat stress response. Trends Plant Sci. 22, 53–65. doi: 10.1016/j.tplants.2016.08.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Olas, J. J., Apelt, F., Annunziata, M. G., John, S., Richard, S. I., Gupta, S., et al. (2021). Primary carbohydrate metabolism genes participate in heat-stress memory at the shoot apical meristem of Arabidopsis thaliana. Mol. Plant 14, 1508–1524. doi: 10.1016/j.molp.2021.05.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Panikulangara, T. J., Eggers-Schumacher, G., Wunderlich, M., Stransky, H., Schöffl, F. (2004). Galactinol synthase1. A novel heat shock factor target gene responsible for heat-induced synthesis of raffinose family oligosaccharides in Arabidopsis. Plant Physiol. 136, 3148–3158. doi: 10.1104/pp.104.042606

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, M. Y., Kang, J. Y., Kim, S. Y. (2011). Overexpression of AtMYB52 confers ABA hypersensitivity and drought tolerance. Mol. Cells 31, 447–454. doi: 10.1007/s10059-011-0300-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Priya, M., Dhanker, O. P., Siddique, K. H. M., HanumanthaRao, B., Nair, R. M., Pandey, S., et al. (2019). Drought and heat stress-related proteins: an update about their functional relevance in imparting stress tolerance in agricultural crops. Theor. Appl. Genet. 132, 1607–1638. doi: 10.1007/s00122-019-03331-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Proels, R. K., Hückelhoven, R. (2014). Cell-wall invertases, key enzymes in the modulation of plant metabolism during defence responses. Mol. Plant Pathol. 15, 858–864. doi: 10.1111/mpp.12139

PubMed Abstract | CrossRef Full Text | Google Scholar

Rajametov, S. N., Yang, E. Y., Cho, M. C., Chae, S. Y., Jeong, H. B., Chae, W. B. (2021). Heat-tolerant hot pepper exhibits constant photosynthesis via increased transpiration rate, high proline content and fast recovery in heat stress condition. Sci. Rep. 11, 14328. doi: 10.1038/s41598-021-93697-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Raza, A., Bhardwaj, S., Atikur Rahman, M., García-Caparrós, P., Habib, M., Saeed, F., et al. (2023). Trehalose: A sugar molecule involved in temperature stress management in plants. Crop J. 12, 1–16. doi: 10.1016/j.cj.2023.09.010

CrossRef Full Text | Google Scholar

Reichelt, N., Korte, A., Krischke, M., Mueller, M. J., Maag, D. (2023). Natural variation of warm temperature-induced raffinose accumulation identifies TREHALOSE-6-PHOSPHATE SYNTHASE 1 as a modulator of thermotolerance. Plant Cell Environ. 46, 3392–3404. doi: 10.1111/pce.14664

PubMed Abstract | CrossRef Full Text | Google Scholar

Saeed, F., Chaudhry, U. K., Raza, A., Charagh, S., Bakhsh, A., Bohra, A., et al. (2023). Developing future heat-resilient vegetable crops. Funct. Integr. Genomics 23, 47. doi: 10.1007/s10142-023-00967-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Sami, F., Yusuf, M., Faizan, M., Faraz, A., Hayat, S. (2016). Role of sugars under abiotic stress. Plant Physiol. Biochem. 109, 54–61. doi: 10.1016/j.plaphy.2016.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Sato, H., Mizoi, J., Shinozaki, K., Yamaguchi-Shinozaki, K. (2024). Complex plant responses to drought and heat stress under climate change. Plant J. 117, 1873–1892. doi: 10.1111/tpj.16612

PubMed Abstract | CrossRef Full Text | Google Scholar

Serrano, N., Ling, Y., Bahieldin, A., Mahfouz, M. M. (2019). Thermopriming reprograms metabolic homeostasis to confer heat tolerance. Sci. Rep. 9, 181. doi: 10.1038/s41598-018-36484-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Shekhawat, K., Almeida-Trapp, M., García-Ramírez, G. X., Hirt, H. (2022). Beat the heat: plant- and microbe-mediated strategies for crop thermotolerance. Trends Plant Sci. 27, 802–813. doi: 10.1016/j.tplants.2022.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Silva, E. N., Silveira, J. A. G., Ribeiro, R. V., Oliveira, J.Ê., Cardoso, R. A. (2018). Photosynthetic and antioxidant responses of jatropha curcas plants to heat stress: on the relative sensitivity of shoots and roots. J. Plant Growth Regul. 37, 255–265. doi: 10.1007/s00344-017-9723-5

CrossRef Full Text | Google Scholar

Song, P., Jia, Q., Xiao, X., Tang, Y., Liu, C., Li, W., et al. (2021). HSP70-3 interacts with phospholipase Dδ and participates in heat stress defense. Plant Physiol. 185, 1148–1165. doi: 10.1093/plphys/kiaa083

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, N., Bassil, E., Hamilton, J. S., Inupakutika, M. A., Zandalinas, S. I., Tripathy, D., et al. (2016). ABA is required for plant acclimation to a combination of salt and heat stress. PloS One 11, e0147625. doi: 10.1371/journal.pone.0147625

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, W., Chen, J., Yue, X., Chai, S., Liu, W., Li, C., et al. (2023). The heat response regulators HSFA1s promote Arabidopsis thermomorphogenesis via stabilizing PIF4 during the day. Sci. Adv. 9, eadh1738. doi: 10.1126/sciadv.adh1738

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, F., Hu, X. L., Yao, T., Yang, X., Chen, J. G., Lu, M. Z., et al. (2021). Recent advances in the roles of HSFs and HSPs in heat stress response in woody plants. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.704905

PubMed Abstract | CrossRef Full Text | Google Scholar

Tiedge, K., Li, X., Merrill, A. T., Davisson, D., Chen, Y., Yu, P., et al. (2022). Comparative transcriptomics and metabolomics reveal specialized metabolite drought stress responses in switchgrass (Panicum virgatum). New Phytol. 236, 1393–1408. doi: 10.1111/nph.18443

PubMed Abstract | CrossRef Full Text | Google Scholar

Tognetti, J. A., Pontis, H. G., Martínez Noël, G. M. A. (2013). Sucrose signaling in plants: a world yet to be explored. Plant Signal Behav. 8, e23316. doi: 10.4161/psb.23316

PubMed Abstract | CrossRef Full Text | Google Scholar

ul Haq, S., Khan, A., Ali, M., Khattak, A. M., Gai, W. X., Zhang, H. X., et al. (2019). Heat shock proteins: dynamic biomolecules to counter plant biotic and abiotic stresses. Int. J. Mol. Sci. 20, 5321. doi: 10.3390/ijms20215321

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Liu, Y., Han, Z., Chen, Y., Huai, D., Kang, Y., et al. (2021a). Integrated transcriptomics and metabolomics analysis reveal key metabolism pathways contributing to cold tolerance in peanut. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.752474

CrossRef Full Text | Google Scholar

Wang, J., Lv, J., Liu, Z., Liu, Y., Song, J., Ma, Y., et al. (2019). Integration of transcriptomics and metabolomics for Pepper (Capsicum annuum L.) in response to heat stress. Int. J. Mol. Sci. 20, 5042. doi: 10.3390/ijms20205042

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Niu, Y., Zheng, Y. (2021b). Multiple functions of MYB transcription factors in abiotic stress responses. Int. J. Mol. Sci. 22, 6125. doi: 10.3390/ijms22116125

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Yan, B., Shi, M., Zhou, W., Zekria, D., Wang, H., et al. (2016). Overexpression of a Brassica campestris HSP70 in tobacco confers enhanced tolerance to heat stress. Protoplasma. 253, 637–645. doi: 10.1007/s00709-015-0867-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Watanabe, S., Matsumoto, M., Hakomori, Y., Takagi, H., Shimada, H., Sakamoto, A. (2014). The purine metabolite allantoin enhances abiotic stress tolerance through synergistic activation of abscisic acid metabolism. Plant Cell Environ. 37, 1022–1036. doi: 10.1111/pce.12218

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen, J., Qin, Z., Sun, L., Zhang, Y., Wang, D., Peng, H., et al. (2023). Alternative splicing of TaHSFA6e modulates heat shock protein-mediated translational regulation in response to heat stress in wheat. New Phytol. 239, 2235–2247. doi: 10.1111/nph.19100

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, X., He, Z., Chen, N., Tang, Z., Wang, Q., Cai, Y. (2019). The roles of environmental factors in regulation of oxidative stress in plant. BioMed. Res. Int. 2019, 9732325. doi: 10.1155/2019/9732325

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, H., Zhang, P., Jiang, C., Wang, Q., Guo, Y., Zhang, X., et al. (2023). Combined transcriptomic and metabolomic analyses of high temperature stress response of quinoa seedlings. BMC Plant Biol. 23, 292. doi: 10.1186/s12870-023-04310-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, S., Liu, Q., Li, W., Yan, J., Fernie, A. R. (2022). Raffinose family oligosaccharides: crucial regulators of plant development and stress responses. Crit. Rev. Plant Sci. 41, 286–303. doi: 10.1080/07352689.2022.2111756

CrossRef Full Text | Google Scholar

Yang, L. Y., Yang, S. L., Li, J. Y., Ma, J. H., Pang, T., Zou, C. M., et al. (2018). Effects of different growth temperatures on growth, development, and plastid pigments metabolism of tobacco (Nicotiana tabacum L.) plants. Bot. Stud. 59, 5. doi: 10.1186/s40529-018-0221-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, G., Wang, L. G., Han, Y., He, Q. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 16, 284–287. doi: 10.1089/omi.2011.0118

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Li, J., Li, M., Zhang, S., Song, S., Wang, W., et al. (2023). NtHSP70-8b positively regulates heat tolerance and seed size in Nicotiana tabacum. Plant Physiol. Biochem. 201, 107901. doi: 10.1016/j.plaphy.2023.107901

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Sun, M., Gao, Y., Luo, Y. (2022). Exogenous trehalose differently improves photosynthetic carbon assimilation capacities in maize and wheat under heat stress. J. Plant Interact. 17, 361–370. doi: 10.1080/17429145.2022.2041119

CrossRef Full Text | Google Scholar

Zhao, Y., Du, H., Wang, Y., Wang, H., Yang, S., Li, C., et al. (2021). The calcium-dependent protein kinase ZmCDPK7 functions in heat-stress tolerance in maize. J. Integr. Plant Biol. 63, 510–527. doi: 10.1111/jipb.13056

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Tian, X., Wang, F., Zhang, L., Xin, M., Hu, Z., et al. (2017). Characterization of wheat MYB genes responsive to high temperatures. BMC Plant Biol. 17, 208. doi: 10.1186/s12870-017-1158-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Zhou, M., Xu, K., Li, J., Li, S., Zhang, S., et al. (2019). Integrated transcriptomics and metabolomics analyses provide insights into cold stress response in wheat. Crop J. 7, 857–866. doi: 10.1016/j.cj.2019.09.002

CrossRef Full Text | Google Scholar

Zheng, Y., Schumaker, K. S., Guo, Y. (2012). Sumoylation of transcription factor MYB30 by the small ubiquitin-like modifier E3 ligase SIZ1 mediates abscisic acid response in Arabidopsis thaliana. Proc. Natl. Acad. Sci. U.S.A. 109, 12822–12827. doi: 10.1073/pnas.1202630109

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: tobacco, heat stress, transcriptome, metabolome, sugar metabolism, purine metabolism

Citation: Chen H, Qiu S, Chen Y, Li J, Xu T, Zhong P, Shao X, Xu S, Ma Z, Huang Z and Pan X (2024) Integrated transcriptomics and metabolomics provides insights into the Nicotiana tabacum response to heat stress. Front. Plant Sci. 15:1425944. doi: 10.3389/fpls.2024.1425944

Received: 30 April 2024; Accepted: 03 July 2024;
Published: 22 July 2024.

Edited by:

Junya Mizoi, The University of Tokyo, Japan

Reviewed by:

Daisuke Todaka, RIKEN Yokohama, Japan
Anthony Guihur, Université de Lausanne, Switzerland

Copyright © 2024 Chen, Qiu, Chen, Li, Xu, Zhong, Shao, Xu, Ma, Huang and Pan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xiaoying Pan, cGFueGlhb3lpbmcwNzA3QDEyNi5jb20=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.