- 1Key Laboratory of Aquatic Biodiversity Research in HeiLongjiang Province, Harbin Normal University, Harbin, China
- 2Heilongjiang Provincial Key Laboratory of Plant Biology in Ordinary Colleges and Universities, Harbin Normal University, Harbin, China
- 3Institute of Botany, Chinese Academy of Sciences, Beijing, China
- 4Botanical Garden, Institute of Botany, Chinese Academy of Sciences, Beijing, China
Lavender (Lavandula angustifolia Mill.) is a valuable aromatic plant with significant commercial importance. However, cold stress–one of the primary abiotic factors impacting sugar metabolism–adversely affects its agricultural productivity in Northeast China. To investigate the mechanisms underlying cold tolerance in L. angustifolia and support economic development, we measured the sugar content and performed transcriptome analysis at temperatures of 30°C (control), 20°C, 10°C, and 0°C. The results revealed that when the temperature dropped from 30°C to 0°C, the amylase activities and the content of maltose and glucose increased, while the starch content decreased. During the process, the up-regulation of LaAMY and LaBAM1/3 suggests an adaptive response in L. angustifolia to cold stress by promoting the breakdown of starch. Meanwhile, the up-regulation of sugar metabolism genes LaRHM1, LaMUR4, LaUGD4, alongside the downregulation of photosynthesis-related genes LaPSAD1, LaPSAN, LaPSBQ2, LaLHCB4.2, and LaPSB27-1 illustrate a strong connection to soluble sugar metabolism. These key genes exhibit significant correlations with starch content and amylase activities, specifically in the decomposition of starch into soluble sugars. The results indicate the decomposition of starch into soluble sugars plays a crucial role in osmotic regulation, facilitating subsequent sugar metabolism in L. angustifolia under cold stress. The correlation between gene expression and physiological indicators suggests that genes can potentially mitigate light-induced damage while promoting cellular homeostasis. Molecular docking analyses between the proteins PSAN and RHM1, MUR4 and UGD4, as well as between LHCB4.2 and RHM1, MUR4, and UGD4 predict that these protein interactions involved in sugar metabolism and photosynthesis contribute to enhancing cold resistance in L. angustifolia.
1 Introduction
Cold stress is caused by sudden temperature drops in autumn and unexpected cold snaps in late spring. For most plants, cold stress can be categorized into two types: chilling stress (below 20°C) and freezing stress (below 0°C) (Dou et al., 2016). The damage induced by cold stress significantly affects growth, health, and distribution limits of plants and causes more worldwide economic losses to agriculture than any other climate-related hazard (Lamichhane, 2021). Plants have undergone intricate evolutionary processes, equipping them adaptive strategies to survive in cold environments (Meng et al., 2025). It is well-documented that starch metabolism can facilitate adaptive changes in source-sink carbon allocation, thereby helping plants withstand cold stress (Dong and Beckles, 2019). Thus, investigating the regulatory mechanisms of starch metabolism underlying plants responses to cold stress has significant implications for agricultural productivity and biodiversity conservation.
The predominant soluble sugars, including glucose, maltose, and sucrose, are crucial for transporting and storing organic materials in both “source” and “sink” tissues of plants. Importantly, accumulation of sugars such as sucrose and glucose under cold stress is thought to protect cellular structures from freezing damage and regulate osmotic pressure. For instance, sucrose acts as a cryoprotectant by reducing the damage of cell membranes at low temperatures and scavenging reactive oxygen species (ROS) (Bai et al., 2023). Meanwhile, the starch and sucrose metabolism involves various enzymes, including β-amylase (BAM), α-amylase (AMY), sucrose synthase (SUS), and starch synthase (SS). Modifying the activities of these enzymes can accelerate the starch-sucrose pathway and enhance soluble sugars formation at low temperatures (Zhang et al., 2017a). Additionally, a study has revealed a significant increase in the levels of glucose and rhamnose in plants when subjected to cold stress (Ghabel and Karamian, 2020). Thus, an investigation of changes into the soluble sugar metabolism under cold stress conditions is beneficial for enhancing plants cold tolerance.
Sugar signals are multifunctional molecules in plants that connect photosynthetic metabolism with cold stress regulation. The regulation of photosynthesis is intricately linked to the metabolism of sugars in plants (Rozpądek et al., 2019). Photo inhibition and light-induced damage can adversely affect energy synthesis and metabolism in plants, especially under cold stress (Guo et al., 2024). This provides a substantial explanation for the effects of low temperatures on sugar metabolism (Hassan et al., 2021). Thus, the stability of the photosynthetic system (PS) at low temperatures is of significant importance. During the process of photosynthesis, energy can be harvested by the light-harvesting chlorophyll a/b binding protein (LHCB) in photosystem II (PSII) (Zhang et al., 2022). The proteins photosystem II subunit Q2 (PSBQ2) and LHCB4.2 in PSII are also crucial for capturing light energy and maintaining the stability of PSII (Wu et al., 2021; Ojosnegros et al., 2022). The photoremediation function of photosystem II subunit Psb27 (PSB27) in PSII is associated with a plant’s ability to adapt to low temperatures. Meanwhile, proteins such as photosystem I subunit D (PSAD) and photosystem I subunit N (PSAN) in photosystem I (PSI) also enhance the efficiency of photosynthesis, with PSAN helping plants acclimate to low temperatures by reducing the formation of ROS (Blanc et al., 2012; Renard et al., 2020).
The starch-sucrose metabolism makes key functions of providing vital substrates for soluble sugar metabolism (Yu et al., 2024). Key genes exert crucial regulatory roles in these metabolic pathway. The synthesis of UDP-rhamnose (UDP-Rha) is catalyzed by the rhamnose synthase (RHM), utilizing UDP-glucose (UDP-Glc), as observed in the substrate of peaches and petunias (Gu et al., 2020; Zhao et al., 2020; Sang et al., 2022). UDP-glucose dehydrogenase (UGD) has been found to convert UDP-Glc to UDP-glucuronic acid (UDP-GlcA) within the UDP-glucose pathway, leading to increased sugar accumulation levels in bamboo and sugarcane (Yang et al., 2020; Mason et al., 2023). UDP-GlcA also serves as a precursor for the synthesis of arabinose (Ara), which can be synthesized from UDP-xylose (UDP-Xyl) through the regulation of the UDP-D-xylose 4-epimerase1 (MUR4) (Dugard et al., 2016; Saez-Aguayo et al., 2017). Additionally, under cold stress conditions, LHCBs positively regulate guard cell signaling in response to abscisic acid (ABA) by modulating ROS homeostasis in Arabidopsis (Xu et al., 2012). Low temperatures increase ABA levels, activating the downstream ABA response promoter ABA-responsive element binding factor 2 ABF2 gene and initiating the ABA-dependent pathway (Chen et al., 2020). Moreover, sugar signals integrate multiple hormone signaling pathways to regulate the overall growth of plants (Chen et al., 2021). The MAP kinase kinase (MKK3) cascade is also involved in mediating the negative feedback regulation of ABA signaling (Mao et al., 2024). Low temperatures trigger the transport of Ca2+ and activates the MAPK cascade, facilitating soluble sugars accumulation and enhancing plant resilience under cold stress conditions (Zhang et al., 2017b).
Lavandula angustifolia, a shrub belonging to the Lamiaceae family, is native to the Mediterranean region (Dobros et al., 2023). The flowers and leaves of L. angustifolia are rich in polysaccharides and other antioxidants usually used in the production of high-value cosmetics, essential oils, and medicines (Kazeminia et al., 2020). Given the diverse economic significance of this species, it has been widely introduced from Yili, Xinjiang, to northeastern China. However, L. angustifolia may be affected by cold stress in these regions. Large temperature fluctuations due to late spring frosts and their effects on L. angustifolia have become more frequent in the expanding cultivation areas in recent years (Supplementary Table S1).
While numerous studies have explored the effects of cold stress on other crops, the specific mechanisms of soluble sugar metabolism and cold tolerance in L. angustifolia remain poorly understood. Previous study indicates that the activation of self-protective mechanisms serves as the key process through which L. angustifolia acclimates to cold conditions (Kazeminia et al., 2020). This involves the activation of specific transcriptional programs, allowing the plant to establish a new balance between development and defense against stressful environments (Zhang et al., 2020). Among these, the BAHD acyltransferase family (BAHD) acyltransferase superfamily and the transcription factor MYC7 (LaMYC7) promoter in L. angustifolia exhibit responsiveness to low-temperature conditions (Dong et al., 2024; Zhang et al., 2024). Our recent research on the transcriptome also suggests that regulating photosynthetic genes may increase soluble sugar content in L. angustifolia, enhancing its cold tolerance (Li et al., 2023a). However, significant gaps still exist in our understanding of how genes related to photosynthesis and sugar metabolism function to enhance resistance in L. angustifolia under low temperature stress.
Little is currently understood regarding the molecular mechanisms of sugar metabolism in lavender, particularly under cold stress conditions. Therefore, the objective of this study was to investigate the key genes and their regulatory networks associated with sugar metabolism in L. angustifolia when subjected to cold stress. Initially, WGCNA alongside gene function annotation methods were utilized to screen key differentially expressed genes (DEGs). A co-expression network of DEGs was established using weighted correlation network analysis (WGCNA) to link physiological changes in soluble sugar metabolism to the transcriptome under cold stress. Through the functional annotation of PFAM genes, we also focused on the Glycosyl hydrolase family and the Glycosyl transferases group. Gene modules related to cold resistance were then identified and analyzed for their primary metabolic pathways and potential functions. Next, pivotal genes were identified based on their connectivity within the corresponding networks, followed by 3D structure prediction and protein docking analysis. Finally, a cold resistance network was constructed for hub genes associated with sugar metabolism in L. angustifolia, providing a theoretical foundation for exploring the molecular mechanisms underlying cold resistance in L. angustifolia.
2 Materials and methods
2.1 Material cultivation
Healthy L. angustifolia plants in the lavender garden at Harbin Normal University (126°32’49″-126°33’E, 45°51’50″-45°52’N) were selected for this study. When these plants are cultivated in the introduced field, temperatures range from 30°C in July (vigorous growth stage) to 0°C in late October (late growth stage). During the seasonal transition from summer to autumn, northwest winds frequently induce rapid temperature drops, with dramatic declines of up to10°C occurring within a single day. As shown in Supplementary Table S1, in areas focusing on lavender cultivation, temperatures have notably declined due to an unexpected late spring, while frost occurrences have risen. To clarify the mechanisms by which lavender responds to decreasing temperatures, we focused our investigation on the physiological and molecular changes that take place as the temperature decreases from 30°C to 0°C. For the above reasons, these plants were cultivated under gradient chilling treatments at 30°C (control), 20°C, 10°C, and 0°C. The plants were cultured in an artificial climate chamber (RGX250E, Tianjin, China) at 30°C for 24 hours (light/dark:12h/12h), then gradually exposed to 20°C, followed by 10°C, and finally to 0°C at a rate of 2°C/hour. They were cultured for 24 hours at each temperature, with three replicates and six plants per replicate for each treatment.
Plants were grouped to ensure that transcriptome sequencing samples were collected from the same plants. Mature leaves were weighed using an electronic analytical balance, quickly placed in liquid nitrogen, and then stored in a freezer (-80°C) for the physiological index measurement and transcriptome sequencing.
2.2 Measurement of soluble sugar metabolism index of L. angustifolia
The soluble sugar metabolism index was assessed using test kits for starch (Sta), glucose (Glu), sucrose (Suc), α-amylase activity (α-AMY), and β-amylase activity (β-AMY) from Solarbio, China, as well as maltose (Mal) from OKA, China. This assessment was conducted following the provided instructions.
2.3 Transcriptome analysis of L. angustifolia
Sample sequencing was performed using the IlluminaHiSeq2500 system by Beijing Ruibo Xingke Biotechnology. The original sequencing data can be accessed in the NCBI database under the accession number PRJNA765132. The dataset was comprised of 12 samples that were exposed to 4 treatment temperatures, with each treatment being biologically replicated three times.
To ensure the quality of data analysis, contaminated and low-quality sequences in the original RawReads were eliminated. The Agilent 2100 RNA Nano 6000 Assay Kit (Agilent Technologies, CA, USA) was used to evaluate RNA sample integrity and concentration. During quality control, reads with longer than 5 bp bases contaminated by linkers, those with over 5% N, and low-quality reads were removed. The GC content of Unigene is 41.02%, while the sequencing error rate for each base was calculated using the Phred score formula (Qphred). Notably, Q30 bases made up more than 90% of the samples. The CleanData obtained after quality control underwent bioinformatic analysis. Trinity was employed for assembling and splicing the remaining sequences to quantitatively analyze the abundance of Unigenes. Based on the genomic data of L. angustifolia (Li et al., 2023b), gene functional annotation analyses were conducted on the transcriptome. The genome and transcriptome originated from the same variety of lavender, ensuring that the comparative results were reliable. The Unigene sequence and predicted ORF sequence, assembled based on Trinity, were compared using Blast. Additional comparisons were made with UniProt, PFAM, and other databases, while Trinotate was used to integrate the comprehensive function annotation results. Through the functional annotation of PFAM genes, we specifically focused on the Glycosyl hydrolase family and the Glycosyl transferases group. Then, five key genes closely associated with sugar metabolism were identified in L. angustifolia.
2.4 Co-expression network analysis of L. angustifolia
The co-expression network of differentially expressed genes (DEGs) in L. angustifolia was analyzed using the WGCNA package on the Annoyun online website (https://c.solargenomics.com). The transcription data was categorized to construct a gene expression network for L. angustifolia under cold stress. Gene modules were identified based on clustering relationships, and co-expression modules specifically responsive to low temperatures were found by analyzing their correlation with sugar metabolism indices. Utilizing the OmicStudio tool (https://www.omicstudio.cn) and the Pearson correlation calculation method, a cluster correlation heat map with markers was generated to identify the core hub genes of L. angustifolia in response to low temperatures. The functional annotation analysis of DEGs with a Q-value ≤ 0.05 was conducted using the Gene Ontology (GO) database (http://geneontology.org/) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (https://www.genome.jp/kegg/).
2.5 Protein interaction network analysis and protein structure prediction of L. angustifolia
The protein sequences were compared with those of Arabidopsis thaliana using the STRING Protein Interaction Database (http://string-db.org/). After prediction by STRING, key protein functions were validated against transcriptome Blast results, confirming consistency with STRING predictions. Key proteins were identified using Cytoscape (version 3.7.2) and the CytoHubba plug-in MCC (Maximal Clique Centrality).
ProtParam was utilized for predicting basic information about the proteins. Prabi was employed for the prediction of protein secondary structure. SMART was used to annotate the protein domains. SWISS-MODEL (https://swissmodel.expasy.org/) was applied for predicting the tertiary structure of the proteins.
We used the Ramachandran Plot from the SAVESv2.0 server to validate the three-dimensional protein model. We also employed ZDOCK3.0.2 software for protein-protein molecular docking and PDBePISA for analyzing protein-protein interaction sites. Visualization of protein-protein docking was carried out using Pymol, and Lig-plot+ software was used to analyze protein interaction forces (Laskowski and Swindells, 2011).
2.6 Validation of RNA-seq data by qRT-PCR
A total of fifteen DEGs, including those closely associated with photosynthesis and carbohydrate metabolism, were selected for qRT-PCR using β-ACTIN as a reference gene. The Primer Premier 6.0 software was utilized for primer design, and the primers used for qRT-PCR can be found in Supplementary Table S2.
Following the manufacturer’s instructions, total RNA was extracted from the L. angustifolia leaves using the MagZol™Reagent (Trizol Reagent) (Magen Biotechnology, China) and reverse transcribed using the Reverse Transciptase KIT(M-MLV) (Beijing Zomen Biotechnology, China). Subsequently, the qRT-PCR was carried out on the Mx3000P system (Agilent Stratagene) using 2×T5 Fast qPCR Mix (SYBRGreenI). The relative expression of each gene was calculated using the 2-ΔΔCT method, and the experiment was conducted in triplicate.
2.7 Statistical analysis
The physiological data of soluble sugar metabolism were statistically processed using Excel (2016) and one-way analysis of variance (one-way ANOVA), and Student’s t-test were performed using SPSS software (IBM SPSS Statistics 25.0) to determine the differences in soluble sugar contents and their statistical significance. The statistical results of the data conform to the normal distribution. The indices of sugar metabolism were determined using a standard curve, ensuring the data obtained is accurate and reliable. The mean ± standard deviation (S.D.) was calculated from three separate experiments for statistical analysis (P<0.05).
The OmicStudio tool and Pearson correlation analysis (https://www.omicstudio.cn) were used to generate a cluster heat map with markers. Then, cluster heat map analyses of key genes involved in sugar metabolism and physiological indexes for soluble sugar metabolism were conducted. Correlation analyses between physiological indexes were performed using CorrPlot within the OmicStudio tools.
3 Results and analysis
3.1 Soluble sugar metabolism in responses to cold stress in L. angustifolia
As the temperature dropped from 30°C to 0°C, starch began to decompose and its contents decreased in L. angustifolia. There was an increase in the contents of the intermediate products glucose and maltose and a decrease in the contents of end product sucrose. These were related to an observed increase in amylase activity, leading to starch decomposition (Figure 1A). The regression analyses between physiological indices and temperatures revealed that the sugar metabolism of L. angustifolia appears to be closely linked to low-temperature stress. The changes in glucose levels under stressful temperatures follow a binomial equation trend, while the changes in other sugar metabolic indexes exhibit a linear decrease with decreased temperatures (Figure 1A).

Figure 1. Analysis of soluble sugar metabolism and DEGs based on WGCNA in L. angustifolia. (A) The changes of soluble sugar contents and amylase activities in responses to cold stress in L. angustifolia. (B) Sample Cluster Tree; (C) Determination of soft threshold of gene co-expression network; (D) Gene clustering and module cutting in co-expression network. (E) Heatmap of association between co-expression network modules of genes and sugar metabolism indices. Sta, starch; Glu, glucose; Mal, maltose; Suc, sucrose; α-AMY, α-amylase activity; β-AMY, β-amylase activity; AMY, amylase. Bars represented means ± S.D. (n = 3).
3.2 Transcriptome analysis of L. angustifolia
3.2.1 Co-expression network analysis of unigenes in cold stress
After performing data analysis and normalization processing, a total of 30,746 high-expression genes were identified. Sample clustering was applied to analyze 12 groups of gene samples, and the results indicated strong gene clustering under low temperatures, validating the reliability of the data (Figure 1B).
Further screening of weight values led us to select β=15 in the construction of the network (Figure 1C). The dynamic cut tree method was then applied to merge the modules with similar expressions. This resulted in a total of 21 co-expression modules (Figure 1D), each represented by a different color. The gene expression level was analyzed to obtain sample clustering and trait correlation through calculation. As shown in Supplementary Figure S1, the results demonstrated a strong correlation between the gene clustering tree and each tissue sample.
The number of genes per module was determined based on expression levels, and the genes with higher clustering degrees were assigned to specific modules. The MEbisque4 and MEdarkgreen contained the largest number of genes, with 1965 and 1280 genes, respectively. The MEgrey consisted of genes that were not assigned to any other modules.
Twenty modules and seven differential glycometabolism physiological indices (starch, glucose, maltose, sucrose, α-amylase activity, β-amylase activity, and amylase) were used in conjunction with the Pearson correlation analysis to examine any existing relationships. The findings revealed that two modules (MEbisque4 and MEdarkgreen) were correlated with the glycometabolism content in L. angustifolia (Figure 1E).
Notably, the MEbisque4 exhibited the strongest correlation with Mal (P<0.01) and demonstrated a significant correlation with starch, sucrose, α-amylase activity, β-amylase activity, and amylase (P<0.01) (Figure 1E). The MEdarkgreen exhibited the strongest correlation with starch, maltose, sucrose, α-amylase activity, β-amylase activity, and amylase (P<0.01), as well as a significant correlation with glucose (P<0.05) (Figure 1E).
3.2.2 The function annotation of DEGs in L. angustifolia
The expression patterns of the MEbisque4 presented notable down-regulated expression under cold stress, whereas those within the MEdarkgreen exhibited significant up-regulated expression (Figures 2A, B). To determine the function of the acquired unigenes, we analyzed the enrichment of GO and KEGG in MEbisque4 and MEdarkgreen, providing annotation information about its functions. Following GO enrichment analysis (Q-value ≤ 0.05), we identified the most enriched GO terms. In the biological process category, the DEGs in the MEdarkgreen were significantly enriched in the regulation of transcription, the oxidation-reduction process (GO:0055114), and the metabolic process (GO:0008152) (Figure 2E). In the biological process category, the MEbisque4 DEGs were enriched in the metabolic process (GO:0008152) and in photosynthesis (GO:0015979) (Figure 2D). Notably, protein binding (GO:0005515) was the predominant molecular function category in these two modules.

Figure 2. MEbisque4 and MEdarkgreen genes function analysis based on WGCNA analysis in L. angustifolia. (A, B) represented heatmap and histogram of key gene expression within the MEbisque4 and MEdarkgreen, respectively. The horizontal axis represented 3 replicates corresponding to temperatures of 30, 20, 10, and 10°C. (C) represented KEGG pathway within the MEbisque4 and MEdarkgreen. The KEGG enrichment factor is depicted by the size of the words in the image; a larger font size corresponds to a higher level of KEGG enrichment, thereby intuitively reflecting the significance of each KEGG pathway. (D, E) represented GO function enrichment within the MEbisque4 and MEdarkgreen.
In KEGG enrichment analysis, it was found that the two modules were enriched in the plant-pathogen interaction pathway and the plant hormone signal transduction pathway. Specifically, the photosynthetic pathway was enriched in the MEbisque4, while the starch, sucrose metabolism, and glycolysis/gluconeogenesis pathways were enriched in the MEdarkgreen (Figure 2C). The KEGG pathway results were similar to GO analysis.
Based on GO analysis, the KEGG pathway, and weight values in the co-expression network, pairs of network relationships between the two modules were selected to generate a co-expression network. The DEGs in the MEbisque4 and MEdarkgreen groups correspond to specific proteins, respectively. We have selected protein names that align with the model organism Arabidopsis thaliana from the STRING database for further investigation of protein interactions (Figures 3A, C). Hub genes were identified using the CytoHubba plugin, and the top 10 DEGs were selected for analysis based on their MCC score (Figures 3B, D). Meanwhile, we screened 5 key DEGs (LaBAM1/3, LaAMY3, LaSUS3, and LaSS4) that are integral to the metabolic pathways of starch and sucrose based on transcriptomic data using gene function annotation methods. We subsequently analyzed their protein interactions with hub genes in the MEbisque4, MEdarkgreen networks and starch metabolic pathway. In total, proteins were included in this analysis, which existed 92 interactions within the protein network (Figure 3E).

Figure 3. Protein network interactions and hub gene screening in L. angustifolia (according to Arabidopsis thaliana database). (A, B) were protein network interactions and hub gene screening in MEbisque4, while (C, D) were in MEdarkgreen. (E) Hub proteins network interactions from the MEbisque4, the MEdarkgreen and starch metabolic pathway in L. angustifolia. AMY3, Alpha-amylase 3; BAM3, Beta-amylase 3; SUS3, Sucrose synthase 3; SS4, Starch synthase 4; BAM1, Beta-amylase 1; LHCB4.2, Light-harvesting complex II chlorophyll a/b binding protein 4.2; MPK3, Mitogen-activated protein kinase 3; TUBB1, Tubulin beta-1 chain; PSBR, Photosystem II reaction center protein R; TIC, Translocon at the inner envelope membrane of chloroplasts; psaD1, Photosystem I reaction center subunit Il D1; UXS4, UDP-glucuronate decarboxylase 4; NSF, N-ethylmaleimide-sensitive factor; HSP70-4, Heat shock protein 70-4; RHM1, Rhamnose biosynthesis 1; PSBW, Photosystem Il reaction center W protein; STR4, Sterol 3-glucosyltransferase 4; UGD4-UDP-glucose dehydrogenase 4; PSB27-1 - Photosystem Il subunit Psb27-1; Uxs2, UDP-glucuronate decarboxylase 2; PSAF, Photosystem I subunit F.
3.2.3 Predictive analysis of key proteins functionality in L. angustifolia
The proteins LaPSAD1, LaPSAN, LaPSBQ2, LaLHCB4.2, LaPSB27-1, LaRHM1, LaMUR4, LaUGD4, LaMPK3, and LaABF2 were determined to be acidic and hydrophilic. Among them, the proteins LaPSB27–1 and LaMPK3 exhibited strong stability, while the stability of LaUGD4 was comparatively lower (Supplementary Table S3).
Protein domains were annotated using SMART (Supplementary Figure S2), revealing that LaPSAD1, LaPSAN, LaPSBQ2, LaLHCB4.2, and LABSB27–1 contain corresponding photosynthetic domains involved in photosynthesis. Furthermore, LaRHM1 and LaMUR4 were found to possess an Epimerase domain, indicating their involvement in important metabolic pathways. LaUGD4 contained three central domains of UDP-glucose (UDPG_MGDP_dh). LaMPK3 possessed a serine/threonine protein kinase (STK) domain, while LaABF2 had a BRLZ (Basic Region and Leucine Zipper, Bzip) domain. We determined that LaABF2 contained two LaPSAD1 and three low-complexity domains (LCDs), which were involved in various cellular processes, ranging from gene expression to signal transduction.
3.3 Analyses of key protein docking
The prediction and construction of the tertiary structures of the above proteins were performed using SWISS-MODEL. The tertiary structure model with the highest confidence (based on GMQE size) was selected. The reliability of the three-level structure model was established using SAVES v6.0 Raman spectrum analysis. This showed that more than 80% of amino acid residues fell within the most acceptable region, indicating that the ten protein models were reliable (Figure 4A).

Figure 4. Three-dimensional structural model and protein docking diagram in L. angustifolia. (A) Three-dimensional structural model and Raman spectrum in the MEbisque4 and MEdarkgreen. Every protein corresponded to Raman spectrum below it. (B–G) Schematic diagram of the 3D protein docking diagram and action structure of amino acid residues.
The ten key proteins were categorized into photosynthetic proteins, glycometabolism proteins, and hormone proteins, with protein molecular docking predictions performed for each protein. The free energy calculations for the docking of PSAN: RHM1, PSAN: MUR4, and PSAN: UGD4, as well as LHCB4.2: RHM1, LHCB4.2: MUR4, and LHCN4.2: UGD4 showed that these interactions had the lower free energy values (Supplementary Table S4). This observation implies that they exhibit a higher likelihood of successful binding. The interaction sites of the protein docking complexes were analyzed using the PDBePISA online program. This program has the capability to predict the formation of hydrogen bonds and salt bridges between amino acid residues in protein-protein interactions.
The interaction sites between PSAN and RHM1, MUR4 and UGD4, as well as between LHCB4.2 and RHM1, MUR4, and UGD4 were predicted. Protein-protein docking was performed using the ZDOCK3.0.2 tool, resulting in the respective structures of the protein docking complexes. Our analysis showed that the best binding sites for these proteins were in the hydrophobic cavities of glycometabolic proteins. We also found that these glycometabolic protein molecules could provide enough space for PSAN and LHCB4.2 (Figures 4B–G).
Upon further examination, it was noted that a unique interaction involving distinct amino acid residues exists between PSAN and RHM1, MUR4 and UGD4 as well as LHCB4.2 and RHM1, and MUR4 and UGD4, (Figures 4B–G). Because of this interaction, the water presence was minimized, facilitating the formation of stable complexes. The molecular docking was bound not only by van der Waals forces and hydrogen bonds, but also by hydrophobic forces (Supplementary Figures S3, S4). The docking of PSAN, glycobotany proteins, and the amino acid residues Ile43, Asn18, and Arg96 were all predicted. We speculated that these sites were the main sites when glycobotany proteins docked with PSAN (Figures 4B–D). Similarly, we found that Arg68 was also an important locus in LHCB 4.2: RHM1 and LHCB 4.2: MUR4 (Figures 4E–G). The theoretical analysis of binding sites and interaction types provided structural information for our experiment.
3.4 Physiological and gene expression verification in L. angustifolia
The hub genes were classified into two distinct categories. LaPSAD1, LaPSAN, LaPSBQ2, LaLHCB4.2, and LaPSB27–1 exhibited significant positive correlations with starch and sucrose but had negative correlation with glucose, maltose, amylase, α-amylase activity, and β-amylase activity. However, LaAMY3, LaBAM1, LaRHM1, LaMUR4, LaUGD4, LaMPK3, and LaABF2 exhibited an opposite correlation with sugar metabolites (Figure 5A). Concurrently, the starch metabolism pathway and the MEbisque4 genes displayed a similar relationship with glucose metabolic processes, suggesting that the MEbisque4 genes are closely associated with glucose metabolism. These findings are consistent with the WGCNA module and effectively demonstrate the reliability of the gene clustering analysis.

Figure 5. Analysis of key genes and sugar metabolism physiological indexes in L. angustifolia. (A) Heat map analysis of key genes in the MEbisque4, the MEdarkgreen and the starch metabolicbased pathway clustered with physiological indexes. Sta, starch; Glu, glucose; Mal, maltose; Suc, sucrose; α-AMY, α-amylase activity; β-AMY, β-amylase activity; AMY, amylase; (B) Heat map analysis between physiological indexes. *P<0.05, **P<0.01. (C) RNA sequencing data accuracy was verified by qRT-PCR in L. angustifolia. Both FPKM and qRT-PCR repeated three times, and the average value is taken for drawing.
Through the correlations among physiological indexes, we observed that starch content was strongly correlated with sucrose and maltose contents, as well as with amylase, α-amylase activity, and β-amylase activity under low-temperature stress (P<0.01, P<0.05) (Figure 5B). Specifically, there was a positive correlation between starch and sucrose, and a negative correlation between starch and maltose, as well as between starch and amylase activities (P<0.01, P<0.05). Sucrose showed a significant correlation with maltose and amylase activities, and maltose was positively correlated with amylase activities, suggesting that the decomposition of starch and the synthesis of sucrose were closely related to AMY activity under low-temperature stress in L. angustifolia. In contrast, glucose, as an intermediate, had little effect under low-temperature stress.
To validate the accuracy of the RNA-seq results, fifteen DEGs were selected for qRT-PCR verification—five genes from the starch metabolic pathway and five each from MEbisque4 and MEdarkgreen. The expression levels of LaAMY3 and LaBAM1/3 exhibited an increasing trend, which facilitated the decomposition and formation of starch. The overall expression of LaSUS3 and LaSS4 were increased. It is speculated that LaSUS3 plays a role in the reversible conversion process of sucrose to glucose and fructose in L.angustifolia, thus promoting an increase in glucose levels in response to cold stress. The quantitative detection findings for these genes were generally consistent with the sequencing data, indicating a high level of reliability in our sequencing results (Figure 5C).
4 Discussion
Plants possess the capacity to accumulate compatible solutes at low temperatures as a result of adaptive evolutionary processes (Zuther et al., 2012). When plants are exposed to cold stress, the total sugar contents can increase by changes of the enzyme activities and acceleration of the starch-sucrose pathway (Li et al., 2023a). These processes cause starch to decompose into glucose and maltose, resulting in sugar accumulation (Zhang et al., 2017a). Meanwhile, the sugars serve as a protective agent for plant cells at low temperature by reducing cell osmotic potential and scavenging active oxygens to increase the cold resistance (Bai et al., 2023). Our study found that as the temperature dropped from 30 °C to 0 °C, the amylase activities significantly increased, while starch contents reduced (Figure 1A). This drop in temperature corresponded with elevated levels of maltose and glucose levels, which were closely associated with heightened amylase activity. This accelerated starch breakdown, providing plants with the energy to adapt to cold conditions and maintain cellular homeostasis and physiological processes (Zhang et al., 2017a).
Plant adaptation to low-temperature stress is governed by intricate physiological and molecular regulatory mechanisms (Liu et al., 2021). Using WGCNA analysis, we divided the transcriptome data into distinct gene co-expression modules, with genes in the same module exhibiting similar expression patterns and carrying out comparable biological functions. We obtained specific MEbisque4 and MEdarkgreen in L. angustifolia at low temperatures, and these modules showed close associations with sugar metabolism and demonstrated high levels of gene clustering (Figure 1E). Five key genes involved in the starch metabolism pathway were also identified, along ten key genes present in MEbisque4 and MEdarkgreen modules. Our findings demonstrated a significant correlation between these key genes and glycometabolism content in L. angustifolia (Figure 5A). This established an important foundation for investigating the key genes associated with soluble sugar metabolism and constructing the regulatory network of soluble sugar metabolism. Various molecular pathways could aid plants in acclimating to cold environments by regulating the expression of a series of cold-response genes (Bian et al., 2022). We specifically identified LaPSB27, LaPSBQ, LaPSAD1, LaPSAN, LaLHCB4.2, LaRHM1, LaMUR4, and LaUGD4 as potential core genes involved in sugar metabolism for adapting to cold stress in L. angustifolia. These genes are involved in the transcription regulation, metabolic processes, photosystem I and II, protein binding, DNA-dependent processes, oxidation-reduction processes, and membrane functions.
The down-regulation of PSB27 has been observed in Arabidopsis under low-temperature stress, resulting in reduced efficiency of PSII while simultaneously providing protection to this complex (Hou et al., 2015). In this study, the expression of LaPSB27 decreased under cold stress, suggesting its protective role in L. angustifolia. Like PSB27, PSBQ has a four-helix bundle tertiary structure and shares a similar binding site, further signifying their close relationship (Gisriel and Brudvig, 2022). The removal of PSBQ alone in Arabidopsis reduced light-harvesting complexes and altered state transitions (Allahverdiyeva et al., 2013). This suggests that L. angustifolia could modulate the energy flow within the PSII system through regulating LaPSBQ expression, thereby contributing to photoprotection (Ojosnegros et al., 2022). Under cold stress, PSAD could down-regulate and block the electron transfer chain, causing the inhibition of PSI (Chitnis et al., 1997). This led to a reduction in photosynthetic activity and triggered a series of protective mechanisms to maintain biological activities (Dong et al., 2022). Similarly, our study observed a down-regulation of LaPSAD1 associated with PSI, which suggests its potential protective role in L. angustifolia under cold stress. Meanwhile, the expressions of LaPSB27, LaPSBQ, and LaPSAD1 were negatively correlated with the activities of glucose, maltose, amylase, α-amylase activity, and β-amylase activity in our study (Figure 5A). The down-regulation in gene expression due to lower temperatures enhanced starch synthesis and its conversion into soluble sugars, boosting L. angustifolia’s resistance to cold.
LHCB is a crucial component of the photosynthetic antenna system, playing a pivotal role in both photosynthesis and stress response regulation (Liu et al., 2021). Under low-temperature stress (4°C), down-regulation of PbrLhcbs resulted in a reduced chlorophyll content and light energy utilization, allowing plants to partially adapt to cold environments (Li et al., 2023c). We found that the expression of LaLHCB4.2 decreased under cold stress, suggesting its protective role in L. angustifolia’s adaptation to low temperatures. LaPSAN were also found to be down-regulated under cold stress, consistent with the results in pepper. PSAN is responsible for mediating LHCII energy transfer to the PSI core, and it is speculated that the LaPSAN protein in PSI plays a significant role in the response of L. angustifolia to cold stress (Vanselow et al., 2009; Liu et al., 2021).
A recent study has revealed that RHMs can directly regulate the expression of key enzyme genes in the glycolytic pathway, impacting glucose decomposition and energy production (Crowe et al., 2024). Our findings indicate that the up-regulating LaRHM1 may enhance energy production at low temperatures in L. angustifolia. UGD is an essential enzyme involved in the biosynthesis of UDP-GlcA (Siddique et al., 2014), an important precursor for plant cell wall polysaccharides. In this study, we proposed that the upregulation of LaUGD4 contributed to the protection of L. angustifolia’s cell wall under cold stress. MUR4 could be involved in converting UDP-D-xylose to L-Ara, thus regulating the content of Ara in the cell wall (Burget et al., 2003). In this study, the expression of the LaMUR4 was significantly up-regulated as temperatures decreased, promoting the formation of protective cell wall substances at low temperatures (Figure 5C).
As previously demonstrated, MUR4 is involved in ABA signal transduction in Arabidopsis (Yan et al., 2020). ABF2 contributes to stress tolerance of mature plants in the ABA signaling pathway (Yoshida et al., 2010). Our results suggest that LaMUR4 and LaABF2 function synergistically in L. angustifolia response to low temperatures. We also identified the involvement of MAPK pathways, which are indirectly associated with sugar metabolism during adaptation to cold. The MAPK cascade mediates plant signal transduction under environmental stress. Overexpression of SIMPK3 significantly enhances the soluble sugar content and cold tolerance of transgenic tomatoes (Yu et al., 2015). Additionally, we observed an increase of LaMPK3 expression after exposure to low temperatures, providing further protection for L. angustifolia.
Protein-protein interactions are crucial for energy, substance metabolism, and gene expression regulation. Our interaction analysis revealed a strong genetic correlation between the MEbisque4 and MEdarkgreen modules from WGCNA (Figure 3E), supporting finding from previous research (Chen et al., 2021). We then utilized computational techniques to generate atomic structure models of protein-protein complexes and efficiently deduced molecular functions in cells through the docking models (Sanyal et al., 2021).
In the protein docking analysis, we identified specific interaction sites involving multiple amino acid residues between LaPSAN and LaRHM1, LaMUR4 and LaUGD4 (Figures 4B–D), as well as between LaLHCB4.2 and LaRHM1, LaMUR4, and LaUGD4 (Figures 4E–G). These interactions contributed to protein stability and may also impacted biological function and activity. Among the six protein-docking results, the residues Ile43, Asn18, Arg96, Arg68, and Tyr81 were frequently found at the binding site. The ion pairing between these residues may decrease the rate constant of intermolecular charge transfer, thereby reducing energy conversion in plants (Qian et al., 2019). This reduction is beneficial for mitigating light damage at low temperatures (Yuan et al., 2021). The proteins interacting with LaPSAN and LaLHCB4.2 play a pivotal role in regulating glucose metabolism and cellular protection, highlighting the significance of photosynthetic regulation in enabling plants to withstand low temperatures. These docking results also facilitate the prediction of protein complex structures and the development of detailed metabolic pathways.
Photosynthesis is crucial for plants to produce soluble sugar, but cold stress can hinder this process and impact sugar metabolism (Ding et al., 2023). Our findings align with previous research that highlight the role of transcriptional regulation of glycometabolism genes in plants adaptation to cold stress (Liu et al., 2021). As shown in Figure 6, when L. angustifolia is exposed to cold stress, starch decomposes into glucose and maltose through amylase LaAMY and LaBAM1/3. Glucose then circulates through G6P and TCA pathways to yield G1P and sucrose (Yuan et al., 2024). Simultaneously, G1P generates UDP-Rha and UDP-GlcA via the intermediate UDP-Glc, facilitated by RHM and UDG enzymes that influence cell wall composition (Zhao et al., 2020; Mason et al., 2023).Furthermore, Ara synthesis occurs from UDP-Xyl under MUR catalysis (Dugard et al., 2016; Saez-Aguayo et al., 2017).

Figure 6. Soluble sugar metabolic pathway in L. angustifolia’s adaptation to low temperatures. Sta, starch; Glu, glucose; Mal, maltose; Suc, sucrose; α-AMY, α-amylase activity; β-AMY, β-amylase activity; AMY, amylase. The blue fonts represented the starch and sucrose metabolic pathway. The red fonts represented key genes that regulated sugar metabolism. The black solid and dashed lines represented direct and indirect actions, respectively. The dashed gray lines depicted molecular docking interactions between proteins. “T” solid line represented inhibitory action.
Plants optimize sugar metabolism for energy production, cellular homeostasis, and enhanced cold tolerance by modulating gene expression. Our analysis of protein-protein interactions revealed that hub genes form complex regulatory networks with other proteins, indicating a coordinated process involving multiple genes and proteins in sugar metabolism under cold stress. Our experiment also verified the correlation between gene expression and physiological indicators, providing deeper insights into L. angustifolia’s adaptability to cold stress.
5 Conclusions
In summary, this study shows that L. angustifolia adapts to cold stress by regulating its sugar metabolism and suppressing photosynthesis. As temperatures decreased from 30°C to 0°C, the degradation of starch into soluble sugars (maltose and glucose) significantly increased (P<0.05). This is attributed to elevated activities of α-amylase and β-amylase, closely associated with the upregulation of LaAMY and LaBAM1/3 genes, which aid in osmotic regulation and energy supply. Transcriptomic analyses revealed an up-regulation of sugar metabolism genes (LaRHM1, LaMUR4, and LaUGD4) alongside a down-regulation of photosynthesis-related genes (LaPSAD1, LaPSAN, LaPSBQ2, LaLHCB4.2, LaPSB27-1), indicating a metabolic shift toward stress mitigation. Molecular docking predictions suggest interactions between sugar metabolism proteins (RHM1, MUR4, UGD4) and photosynthesis components (PSAN, LHCB4.2), implying a coordinated mechanism for balancing energy allocation while minimizing light-induced damage under cold conditions. These findings highlight the importance of starch-to-sugar conversion and transcriptional regulation in enhancing cold tolerance in L. angustifolia. This mechanistic insight is vital for breeding cold-resistant lavender varieties and optimizing cultivation strategies in low-temperature regions such as Northeast China, thus supporting agricultural productivity and economic development.
Data availability statement
The original sequencing data can be accessed in the NCBI database under the accession number PRJNA765132.
Author contributions
YL: Formal Analysis, Investigation, Methodology, Visualization, Writing – original draft. YL: Project administration, Resources, Writing – review & editing. YW: Formal Analysis, Investigation, Writing – review & editing. RY: Writing – review & editing. ZY: Conceptualization, Funding acquisition, Supervision, Writing – review & editing. HL: Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Project of expenditures for basic research in Heilongjiang provincial universities and Digital economy, bio-economy, carbon peaking and carbon neutrality strategy, and creative design special projects of Harbin Normal University (No.XGT202307).
Conflict of interest
The 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.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
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.2025.1537516/full#supplementary-material
References
Allahverdiyeva, Y., Suorsa, M., Rossi, F., Pavesi, A., Kater, M. M., Antonacci, A., et al. (2013). Arabidopsis plants lacking PsbQ and PsbR subunits of the oxygen-evolving complex show altered PSII super-complex organization and short-term adaptive mechanisms. Plant J. 75, 671–684. doi: 10.1111/tpj.12230, PMID: 23647309
Bai, L., Zhu, H., Shi, Y., Li, Y., Miao, Y., Yu, X., et al. (2023). Antisense overexpression of Gγ Subunit csGG3. 1–2 reduces soluble sugar content and chilling tolerance in cucumber. Horticulturae. 9, 240. doi: 10.3390/horticulturae9020240
Bian, X., Jiang, H., Meng, Y., Li, Y. P., Fang, J., and Lu, Z. (2022). Regulation of gene expression by glycolytic and gluconeogenic enzymes. Trends Cell Biol. 32, 786–799. doi: 10.1016/j.tcb.2022.02.003, PMID: 35300892
Blanc, G., Agarkova, I., Grimwood, J., Kuo, A., Brueggeman, A., Dunigan, D. D., et al. (2012). The genome of the polar eukaryotic microalga Coccomyxa subellipsoidea reveals traits of cold adaptation. Genome Biol. 13, 1–12. doi: 10.1186/gb-2012-13-5-r39, PMID: 22630137
Burget, E. G., Verma, R., Mølhøj, M., and Reiter, W. D. (2003). The biosynthesis of L-arabinose in plants: molecular cloning and characterization of a Golgi-localized UDP-D-xylose 4-epimerase encoded by the MUR4 gene of Arabidopsis. Plant Cell 15, 523–531. doi: 10.1105/tpc.008425, PMID: 12566589
Chen, K., Li, G. J., Bressan, R., Song, C. P., Zhu, J. K., and Zhao, Y. (2020). Abscisic acid dynamics, signaling, and functions in plants. J. Integr. Plant Biol. 62, 25–54. doi: 10.1111/jipb.12899, PMID: 31850654
Chen, Q. S., Zhang, J., and Li, G. (2021). Dynamic epigenetic modifications in plant sugar signal transduction. Trends Plant Sci. 27, 379–390. doi: 10.1016/j.tplants.2021.10.009, PMID: 34865981
Chitnis, V. P., Ke, A., and Chitnis, P. R. (1997). The PsaD subunit of photosystem I. (mutations in the basic domain reduce the level of PsaD in the membranes). Plant Physiol. 115, 1699–1705. doi: 10.1104/pp.115.4.1699, PMID: 9414569
Crowe, A. S., Liu, Y., Zhao, X., Scheller, H. V., and Keasling, J. D. (2024). Advances in engineering nucleotide sugar metabolism for natural product glycosylation in Saccharomyces cerevisiae. ACS Synthetic Biol. 13, 1589–1599. doi: 10.1021/acssynbio.3c00737, PMID: 38820348
Ding, Y., Wang, X. T., Wang, F., Shao, Y. L., Zhang, A. M., and Chang, W. (2023). The effects of chilling stress on antioxidant enzymes activities and proline, malondialdehyde, soluble sugar contents in three Paphiopedilum species. Russian J. Plant Physiol. 70, 61. doi: 10.1134/S1021443722603184
Dobros, N., Zawada, K. D., and Paradowska, K. (2023). Phytochemical profiling, antioxidant and anti-inflammatory activity of plants belonging to the Lavandula genus. Molecules 28, 256. doi: 10.3390/molecules28010256, PMID: 36615453
Dong, S. and Beckles, D. M. (2019). Dynamic changes in the starch-sugar interconversion within plant source and sink tissues promote a better abiotic stress response. J. Plant Physiol. 234–235), 80–93. doi: 10.1016/j.jplph.2019.01.007, PMID: 30685652
Dong, S., Jiang, Y., Bu, Y., Wang, S., Zhang, H., and Wang, R. (2022). Growth, photosynthetic pigment proteins responses and transcriptome analysis provide insights into survival strategies against short-term cold stress in the blue-green algae, Arthrospira. Aquaculture Rep. 27, 101403. doi: 10.1016/j.aqrep.2022.101403
Dong, Y., Wei, Z., Zhang, W., Li, J., Han, M., Bai, H., et al. (2024). LaMYC7, a positive regulator of linalool and caryophyllene biosynthesis, confers plant resistance to Pseudomonas syringae. Horticulture Res. 11, uhae044. doi: 10.1093/hr/uhae044, PMID: 38623075
Dou, M., Cheng, S., Zhao, B., Xuan, Y., and Shao, M. (2016). The indeterminate domain protein ROC1 regulates chilling tolerance via activation of DREB1B/CBF1 in rice. Int. J. Mol. Sci. 17, 233. doi: 10.3390/ijms17030233, PMID: 26927068
Dugard, C. K., Mertz, R. A., Rayon, C., Mercadante, D., Hart, C., Benatti, M. R., et al. (2016). The cell wall arabinose-deficient Arabidopsis thaliana mutant murus5 encodes a defective allele of REVERSIBLY GLYCOSYLATED POLYPEPTIDE2. Plant Physiol. 171, 1905–1920. doi: 10.1104/pp.15.01922, PMID: 27217494
Ghabel, V. K. and Karamian, R. (2020). Effects of TiO2 nanoparticles and spermine on antioxidant responses of Glycyrrhiza glabra L. @ to cold stress. Acta Botanica Croatica 79, 137–147. doi: 10.37427/botcro-2020-025
Gisriel, C. J. and Brudvig, G. W. (2022). Comparison of PsbQ and Psb27 in photosystem II provides insight into their roles. Photosynthesis Res. 152, 177–191. doi: 10.1007/s11120-021-00888-2, PMID: 35001227
Gu, N., Qiu, C., Zhao, L., Zhang, L., and Pei, J. (2020). Enhancing UDP-rhamnose supply for rhamnosylation of flavonoids in Escherichia coli by regulating the modular pathway and improving NADPH availability. J. Agric. Food Chem. 68, 9513–9523. doi: 10.1021/acs.jafc.0c03689, PMID: 32693583
Guo, M., Li, Z., Wang, L., Xu, T., Huang, H., Kanwar, M. K., et al. (2024). BAG8 positively regulates cold stress tolerance by modulating photosystem, antioxidant system and protein protection in Solanum lycopersicum. Plant Physiol. Biochem. 206, 108267. doi: 10.1016/j.plaphy.2023.108267, PMID: 38091937
Hassan, M. A., Xiang, C., Farooq, M., Muhammad, N., Yan, Z., Hui, X., et al. (2021). Cold stress in wheat: plant acclimation responses and management strategies. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.676884, PMID: 34305976
Hou, X., Fu, A., Garcia, V. J., Buchanan, B. B., and Luan, S. (2015). PSB27: a thylakoid protein enabling Arabidopsis to adapt to changing light intensity. Proc. Natl. Acad. Sci. 112, 1613–1618. doi: 10.1073/pnas.1424040112, PMID: 25605904
Kazeminia, M., Abdi, A., Vaisi-Raygani, A., Jalali, R., Shohaimi, S., Daneshkhah, A., et al. (2020). The effect of lavender (Lavandula stoechas L.) on reducing labor pain: A systematic review and meta-analysis. Evidence-Based Complementary Altern. Med. 2020, 4384350. doi: 10.1155/2020/4384350, PMID: 33224252
Lamichhane, J. R. (2021). Rising risks of late-spring frosts in a changing climate. Nat. Climate Change. 11, 554–555. doi: 10.1038/s41558-021-01090-x
Laskowski, R. A. and Swindells, M. B. (2011). LigPlot+: multiple ligand–protein interaction diagrams for drug discovery. J. Chem. Inf. Modeling 51, 2778–2786. doi: 10.1021/ci200227u, PMID: 21919503
Li, J. R., Li, H., Wang, Y. M., Zhang, W. Y., Wang, D., Dong, Y. M., et al. (2023b). Decoupling subgenomes within hybrid lavandin provide new insights into speciation and monoterpenoid diversification of Lavandula. Plant Biotechnol. J. 21, 2084–2099. doi: 10.1111/pbi.14115, PMID: 37399213
Li, L., Liang, Y. C., Liu, Y. N., Sun, Z. Y., Liu, Y. N., Yuan, Z. N., et al. (2023a). Transcriptome analyses reveal photosynthesis-related genes involved in photosynthetic regulation under low temperature stress in Lavandula angustifolia Mill. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1268666, PMID: 38107014
Li, X. L., Jiang, Z. Y., Zhang, C. F., Cai, K. F., Wang, H., Pan, W. Y., et al. (2023c). Comparative genomics analysis provide insights into evolution and stress responses of Lhcb genes in Rosaceae fruit crops. BMC Plant Biol. 23, 484–484. doi: 10.1186/s12870-023-04438-x, PMID: 37817059
Liu, Z. B., Song, J. S., Miao, W., Yang, B. Z., Zhang, Z. Q., Chen, W. C., et al. (2021). Comprehensive proteome and lysine acetylome analysis reveals the widespread involvement of acetylation in cold resistance of pepper (Capsicum annuum L.). Front. Plant Sci. 12, 730489. doi: 10.3389/fpls.2021.730489, PMID: 34512705
Mao, X., Zheng, X., Sun, B., Jiang, L., Zhang, J., Lyu, S., et al. (2024). MKK3 cascade regulates seed dormancy through a negative feedback loop modulating ABA signal in rice. Rice 17, 2. doi: 10.1186/s12284-023-00679-4, PMID: 38170405
Mason, P. J., Hoang, N. V., Botha, F. C., Furtado, A., Marquardt, A., and Henry, R. J. (2023). Organ-specific expression of genes associated with the UDP-glucose metabolism in sugarcane (Saccharum spp. hybrids). BMC Genomics 24, 18. doi: 10.1186/s12864-023-09124-8, PMID: 36639618
Meng, L., Zhou, H. M., Tan, L., Li, Q. Y., Hou, Y. J., Li, W. J., et al. (2025). VaWRKY65 contributes to cold tolerance through dual regulation of soluble sugar accumulation and reactive oxygen species scavenging in Vitis amurensis. Horticulture Res. 12, uhae367. doi: 10.1093/hr/uhae367, PMID: 40078721
Ojosnegros, S., Alvarez, J. M., Grossmann, J., Gagliardini, V., Quintanilla, L. G., Grossniklaus, U., et al. (2022). The shared proteome of the apomictic fern Dryopteris affinis ssp. affinis and its sexual relative Dryopteris oreades. Int. J. Mol. Sci. 23, 14027. doi: 10.3390/ijms232214027, PMID: 36430514
Qian, C., Ivan, A., Sergey, O., and Baker, D. (2019). Protein interaction networks revealed by proteome coevolution. Sci. (New York). 365, 185–189. doi: 10.1126/science.aaw6718, PMID: 31296772
Renard, J., Niñoles, R., Martínez-Almonacid, I., Gayubas, B., Mateos-Fernández, R., Bissoli, G., et al. (2020). Identification of novel seed longevity genes related to oxidative stress and seed coat by genome-wide association studies and reverse genetics. Plant Cell Environ. 43, 2523–2539. doi: 10.1111/pce.13822, PMID: 32519347
Rozpądek, P., Nosek, M., Domka, A., Ważny, R., Jędrzejczyk, R., Tokarz, K., et al. (2019). Acclimation of the photosynthetic apparatus and alterations in sugar metabolism in response to inoculation with endophytic fungi. Plant Cell Environ. 42, 1408–1423. doi: 10.1111/pce.13485, PMID: 30516827
Saez-Aguayo, S., Rautengarten, C., Temple, H., Sanhueza, D., Ejsmentewicz, T., Sandoval-Ibañez, O., et al. (2017). UUAT1 is a Golgi-localized UDP-uronic acid transporter that modulates the polysaccharide composition of Arabidopsis seed mucilage. Plant Cell 29, 129–143. doi: 10.1105/tpc.16.00465, PMID: 28062750
Sang, L., Chen, G., Cao, J., Liu, J., and Yu, Y. (2022). PhRHMs play important roles in leaf and flower development and anthocyanin synthesis in petunia. Physiologia Plantarum 174, e13773. doi: 10.1111/ppl.13773, PMID: 36066309
Sanyal, A., Zbornik, E. A., Watson, B. G., Christoffer, C., Ma, J., Kihara, D., et al. (2021). Kinetic and structural parameters governing Fic-mediated adenylylation/AMPylation of the Hsp70 chaperone, BiP/GRP78. Cell Stress Chaperones 26, 639–656. doi: 10.1007/s12192-021-01208-2, PMID: 33942205
Siddique, S., Endres, S., Sobczak, M., Radakovic, Z. S., Fragner, L., Grundler, F. M., et al. (2014). Myo-inositol oxygenase is important for the removal of excess myo-inositol from syncytia induced by Heterodera schachtii in Arabidopsis roots. New Phytol. 201, 476–485. doi: 10.1111/nph.12535, PMID: 24117492
Vanselow, C., Weber, A. P., Krause, K., and Fromme, P. (2009). Genetic analysis of the Photosystem I subunits from the red alga, Galdieria sulphuraria. Biochim. Biophys. Acta (BBA)-Bioenergetics 1787, 46–59. doi: 10.1016/j.bbabio.2008.10.004, PMID: 19007746
Wu, X., Khan, R., Gao, H., Liu, H., Zhang, J., and Ma, X. (2021). Low light alters the photosynthesis process in cigar tobacco via modulation of the chlorophyll content, chlorophyll fluorescence, and gene expression. Agriculture 11, 755. doi: 10.3390/agriculture11080755
Xu, Y. H., Liu, R., Yan, L., Liu, Z. Q., Jiang, S. C., Shen, Y. Y., et al. (2012). Light-harvesting chlorophyll a/b-binding proteins are required for stomatal response to abscisic acid in Arabidopsis. J. Exp. Bot. 63, 1095–1106. doi: 10.1093/jxb/err315, PMID: 22143917
Yan, J., Fang, L., Yang, L., He, H., Huang, Y., Liu, Y., et al. (2020). Abscisic acid positively regulates l-arabinose metabolism to inhibit seed germination through ABSCISIC ACID INSENSITIVE4-mediated transcriptional promotions of MUR4 in Arabidopsis thaliana. New Phytol. 225, 823–834. doi: 10.1111/nph.16149, PMID: 31461539
Yang, Y., Kang, L., Wu, R., Chen, Y., and Lu, C. (2020). Genome-wide identification and characterization of UDP-glucose dehydrogenase family genes in moso bamboo and functional analysis of PeUGDH4 in hemicellulose synthesis. Sci. Rep. 10, 10124. doi: 10.1038/s41598-020-67227-8, PMID: 32576917
Yoshida, T., Fujita, Y., Sayama, H., Kidokoro, S., Maruyama, K., Mizoi, J., et al. (2010). AREB1, AREB2, and ABF3 are master transcription factors that cooperatively regulate ABRE-dependent ABA signaling involved in drought stress tolerance and require ABA for full activation. Plant J. 61, 672–685. doi: 10.1111/j.1365-313x.2009.04092.x, PMID: 19947981
Yu, J. Q., Li, Z. T., Chen, S., Gao, H. S., and Sheng, L. X. (2024). Analysis of ethylene signal regulating sucrose metabolism in strawberry fruits based on RNA-seq. Plants 13, 1121. doi: 10.3390/plants13081121, PMID: 38674530
Yu, L., Yan, J., Yang, Y., and Zhu, W. (2015). Overexpression of tomato mitogen-activated protein kinase SlMPK3 in tobacco increases tolerance to low temperature stress. Plant Cell Tissue Organ Culture 121, 21–34. doi: 10.1007/s11240-014-0675-1
Yuan, X., Stephanie, C. S., and Thomas, D. S. (2024). The oxidative pentose phosphate pathway in photosynthesis: a tale of two shunts. New Phytol. 6, 2453–2463. doi: 10.1111/nph.19730, PMID: 38567702
Yuan, L., Zhang, L., Wu, Y., Zheng, Y., Nie, L., Zhang, S., et al. (2021). Comparative transcriptome analysis reveals that chlorophyll metabolism contributes to leaf color changes in wucai (Brassica campestris L.) in response to cold. BMC Plant Biol. 21, 438–438. doi: 10.1186/s12870-021-03218-9, PMID: 34583634
Zhang, H., Hou, J., Liu, J., Zhang, J., Song, B., and Xie, C. (2017a). The roles of starch metabolic pathways in the cold-induced sweetening process in potatoes. Starch - Stärke 69, 1600194. doi: 10.1002/star.201600194
Zhang, W., Li, J., Dong, Y., Huang, Y., Qi, Y., Bai, H., et al. (2024). Genome-wide identification and expression of BAHD acyltransferase gene family shed novel insights into the regulation of linalyl acetate and lavandulyl acetate in lavender. J. Plant Physiol. 292, 154143. doi: 10.1016/j.jplph.2023.154143, PMID: 38064887
Zhang, Z., Li, J., Li, F., Liu, H., Yang, W., Chong, K., et al. (2017b). OsMAPK3 phosphorylates OsbHLH002/OsICE1 and inhibits its ubiquitination to activate OsTPP1 and enhances rice chilling tolerance. Dev. Cell 43, 731–743. doi: 10.1016/j.devcel.2017.11.016, PMID: 29257952
Zhang, Y., Raza, A., Huang., H., Su, W., Luo., D., Zeng., L., et al. (2022). Analysis of Lhcb gene family in rapeseed (Brassica napus L.) identifies a novel member “BnLhcb3.4” modulating cold tolerance. Environ. Exp. Bot. 198, 104848. doi: 10.1016/j.envexpbot.2022.104848
Zhang, H., Zhao, Y., and Zhu, J. K. (2020). Thriving under stress: how plants balance growth and the stress response. Dev. Cell 55, 529–543. doi: 10.1016/j.devcel.2020.10.012, PMID: 33290694
Zhao, Z., Ren, C., Xie, L., Xing, M., Zhu, C., Ji, R., et al. (2020). Functional analysis of PpRHM1 and PpRHM2 involved in UDP-l-rhamnose biosynthesis in Prunus persica. Plant Physiol. Biochem. 155, 658–666. doi: 10.1016/j.plaphy.2020.08.011, PMID: 32861032
Keywords: Lavandula angustifolia Mill., cold stress, sugar metabolism, key genes, molecular docking
Citation: Liang Y, Liu Y, Wang Y, Yang R, Yuan Z and Li H (2025) Investigation of the physiological and molecular regulatory mechanism of soluble sugar metabolism in Lavandula angustifolia Mill. under cold stress. Front. Plant Sci. 16:1537516. doi: 10.3389/fpls.2025.1537516
Received: 01 December 2024; Accepted: 16 June 2025;
Published: 08 July 2025.
Edited by:
Klára Kosová, Crop Research Institute (CRI), CzechiaReviewed by:
Muhammad Ahmad Hassan, Anhui Agricultural University, ChinaZhen Li, International Center for Bamboo and Rattan, China
Copyright © 2025 Liang, Liu, Wang, Yang, Yuan and Li. 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: Zening Yuan, eGlhb3l1YW4xNjhva0AxNjMuY29t; Hui Li, bGlodWliamZ1QDEyNi5jb20=